Abstract
Postharvest processing critically determines the quality of Codonopsis pilosula, yet the molecular mechanisms underlying the traditional “rubbing–sweating” technique remain unclear. We hypothesized that rubbing–sweating imposes stronger abiotic stress than shade drying, thereby activating stress-responsive pathways and enhancing the accumulation of bioactive constituents. To test this, freshly harvested roots were processed by shade drying (SD) and rubbing–sweating drying (RD), and compared with fresh controls (FC) in terms of chemical composition, antioxidant enzyme activity, and transcriptomic profiles. After 6 days, RD significantly increased lobetyolin content by 15.3% relative to FC and 9.7% relative to SD (p<0.01), while polysaccharides reached 19.5% in RD versus 10.6% in FC (p<0.05). Antioxidant enzymes also exhibited marked increases under RD, with catalase activity elevated by 235% compared to FC. Transcriptome sequencing revealed 17,338 DEGs in RD vs. SD and 11,007 in RD vs. FC, enriched in MAPK signaling, hormone transduction, and aromatic amino acid biosynthesis. These findings support our hypothesis that rubbing–sweating enhances the medicinal quality of C. pilosula through stress-induced activation of metabolic pathways. This work provides the first transcriptomic evidence for the molecular basis of this traditional technique, offering new insights for optimizing and modernizing postharvest processing.
1 Introduction
In recent years, rising living standards and advancements in the pharmaceutical industry have led to a growing annual demand for high-quality traditional Chinese medicinal (TCM) materials with health-promoting properties. The quality of TCM is determined by a complex interplay of factors, including the growing environment, cultivation practices, and duration of growth, as well as appropriate post-harvest processing methods at the production site, which are critical to the final quality of the medicinal materials (; ). Common traditional post-harvest processing methods for Chinese medicinal materials include steaming, boiling, blanching, smoking and “Rubbing-sweating” (; ). However, some of these methods are labor-intensive, time-consuming, and technically complex, rendering them increasingly incompatible with the modernization needs of the traditional Chinese medicine (TCM) industry. Therefore, it is essential to conduct in-depth investigations into traditional processing techniques to elucidate their underlying scientific mechanisms. Such research can offer valuable insights for developing processing methods that not only improve the quality of medicinal materials but also align with the demands of modern industrial production.
Codonopsis Radix (Dangshen) refers to the dried roots of perennial species in the Campanulaceae family, including Codonopsis pilosula (Franch.) Nannf., C. pilosula Nannf. var. modesta (Nannf.) L.T. Shen, and C. tangshen Oliv (). Due to its functions in strengthening the spleen, moistening the lungs, nourishing the blood, promoting the generation of bodily fluids, modulating immune function, and exhibiting antitumor properties (; ), C. pilosula has been used for centuries in both food and traditional folk medicine across Asian countries, including China, Japan, and Korea (). In the traditional production areas of C. pilosula, a distinctive post-harvest processing technique known as “Rubbing-sweating” has gradually evolved through long-term practice. Locally, it is widely believed that C. pilosula processed using this method possesses a superior commercial appearance and enhanced quality. The traditional rubbing and sweating process consists of the following steps: (1) Freshly harvested roots are cleaned to remove impurities and shade-dried for 2–3 days until their texture changes from hard to pliable; (2) The roots are bundled and vigorously rubbed by hand or with mechanical assistance to make the tissues denser and firmer; (3) The rubbed roots are heaped and covered with straw mats to allow surface moisture to condense. This sweating step lasts 1–2 days, with careful monitoring to prevent mold or spoilage; (4) The roots are then sun-dried for 1–2 days. Steps (2) and (3) are typically repeated 2–3 times until the roots are fully dried ().
Freshly harvested plant tissues—particularly roots and rhizomes—retain substantial physiological activity for a period following harvest. During this time, they can respond to external stimuli by initiating a series of physiological and biochemical reactions, thereby influencing metabolic processes and the accumulation of secondary metabolites within the plant (). Studies have demonstrated that appropriate postharvest processing methods can significantly enhance the accumulation of volatile C6-compounds in oolong tea (). Similarly, postharvest high-temperature treatment of jujube fruit (soaking in 50°C water for 4 minutes) has been shown to effectively delay senescence and decay (). In addition, the traditional “sweating” process has been employed in the postharvest handling of several Chinese medicinal materials, including Magnolia officinalis, Salvia miltiorrhiza, and Gentiana macrophylla. Some researchers have suggested that sweating alters the composition of microbial communities in M. officinalis, thereby increasing the content of active compounds such as magnolol and honokiol in the final product (). Others have found that sweating significantly promotes the production of intermediate metabolites involved in the biosynthesis of tanshinones and salvianolic acids in S. miltiorrhiza, contributing to the accumulation of these bioactive components and improving the quality of the dried material (). Furthermore, sweating has been reported to enhance the biosynthesis of phenylalanine, tyrosine, and tryptophan, as well as polyphenol and terpene biosynthesis in G. macrophylla, thus promoting the accumulation of key active ingredients and improving the overall quality of the herb ().
At present, it remains unclear how the traditional “rubbing–sweating” process affects the active components of C. pilosula and its molecular regulatory mechanism. We hypothesized that the “rubbing–sweating” treatment might enhance the quality of C. pilosula by inducing the stress response signaling pathway and promoting the accumulation of bioactive compounds. To verify this hypothesis, this study compared the differences in chemical components, antioxidant enzyme activities, and transcriptome profiles among fresh C. pilosula, shade-dried C. pilosula, and C. pilosula treated with the “rubbing–sweating” process. The aim of this study was to clarify the molecular mechanism by which the “rubbing–sweating” treatment improves the quality of C. pilosula, providing a theoretical basis for the optimization of traditional processing techniques and the development of modern drying technologies.
2 Materials and methods
2.1 Plant material and experimental design
Root samples of C. pilosula were collected from Wen County, Gansu Province, China (32.944° N, 104.683° E), and were of uniform age (3 years). Botanical identification was performed by Professor Chen Yuan from the College of Agronomy, Gansu Agricultural University, who confirmed the samples as members of the Campanulaceae family, specifically Codonopsis pilosula (Franch.) Nannf.
The collected C. pilosula samples, selected for uniform size, were divided into three groups, each consisting of 12 roots. One group was designated as the control (fresh C. pilosula, FC), while the other two groups were subjected to shade-drying (SD) and “Rubbing-sweating” drying (RD) treatments, respectively. For subsequent analyses, three biological replicates were established per treatment, with each replicate comprising pooled tissues from four randomly selected roots.
2.1.1 Fresh C. pilosula
The freshly harvested C. pilosula samples were immediately stored at −80°C to preserve them for subsequent analysis.
2.1.2 “Rubbing-sweating” dried C. pilosula
The samples were kept indoors and manually rubbed every three days. After each rubbing session, the roots were covered with a breathable tarp for 12 hours to facilitate sweating, followed by air-drying for two days to promote moisture loss and prevent mold growth. This cycle was repeated three times. The ambient temperature was maintained at 2–10°C throughout the process.
2.1.3 Shade-dried C. pilosula
The samples were placed in a cool, shaded environment under the same conditions as the RD group, except without rubbing or sweating. They were dried continuously for 12 days, with air-drying every three days to facilitate moisture loss and prevent mold formation.
To investigate transcriptomic changes during the drying process, samples were collected after four days of treatment. Each sample weighed 9 g and was immediately stored in an ultra-low-temperature freezer (−80°C) for further analysis. Voucher specimens of C. pilosula used in this study were deposited in the Herbarium of Gansu Agricultural University (College of Agronomy) to facilitate future research and verification.
2.2 Detection of chemical composition
Freeze-dried C. pilosula tissue (2 g) was ground into a fine powder and extracted with 20 mL of 50% ethanol using ultrasonic assistance at 60°C for 60 minutes. The extract was then centrifuged at 1000 rpm for 5 minutes, and the resulting supernatant was filtered through a 0.22 μm organic membrane filter. The filtrate was subsequently analyzed using a high-performance liquid chromatography with photodiode array detection (HPLC-PDA) system. Quantitative determination of five major chemical constituents—adenosine, protocatechuic acid, tryptophan, syringin, and lobetyolin—was performed using an Agilent 1260 liquid chromatograph equipped with a Kromasil C18 reversed-phase column (4.6 mm × 250 mm). The mobile phase consisted of 0.5% formic acid aqueous solution (solvent A) and acetonitrile (solvent B), delivered at a flow rate of 1 mL·min-¹. The gradient elution program was as follows: 0–5 min, 5% B; 5–10 min, 5–8% B; 10–15 min, 8% B; 15–25 min, 8–14% B; 25–30 min, 14–20% B; 30–40 min, 20–30% B; and 40–50 min, 30–50% B. The injection volume was 10 μL, and detection was carried out at 260 nm.
The contents of polysaccharides, proteins, and ash were determined in accordance with AOAC standard methods (; ).
2.3 Superoxide dismutase, peroxidase, catalase, and ascorbate peroxidase activity assay
A total of 1.00 g of fresh C. pilosula tissue powder was mixed with 9 mL of normal saline for enzymatic extraction. The mixture was centrifuged at 2,500 rpm for 10 minutes at 4°C, and the resulting supernatant was collected for the determination of superoxide dismutase (SOD), peroxidase (POD), catalase (CAT), and ascorbate peroxidase (APX) activities. Enzyme activity assay kits were purchased from Nanjing Jiancheng Bioengineering Institute.
2.4 Transcriptome data determination
2.4.1 RNA abstraction and quality evaluation
Samples from FC, SD (6 days), and RD (6 days) groups (0.5 g fresh weight each) were selected for total RNA extraction using TRIzol reagent (Invitrogen, Carlsbad, CA, USA), following the manufacturer’s instructions. The integrity and concentration of the extracted RNA were assessed using a 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) and quantified using a NanoDrop ND-2000 spectrophotometer (Thermo Scientific, Wilmington, DE, USA). High-quality RNA was subsequently used for library construction.
2.4.2 Complementary DNA library transcription group fabrication and sequencing
RNA purification, reverse transcription, library construction, and sequencing were performed by EMajorbio Bio-Pharm Biotechnology Co., Ltd. (Shanghai, China). RNA-seq libraries were prepared using the TruSeq™ RNA Sample Preparation Kit (Illumina, USA). Poly(A) mRNA was isolated from total RNA using oligo(dT)-attached magnetic beads, followed by fragmentation using a fragmentation buffer. The resulting short RNA fragments were used as templates to synthesize double-stranded complementary DNA (ds-cDNA) using the SuperScript Double-Stranded cDNA Synthesis Kit (Invitrogen) and RHP reagents (Illumina). Subsequently, the cDNA underwent end repair, phosphorylation, and addition of a single ‘A’ base according to Illumina’s library preparation protocol. Libraries were size-selected for cDNA fragments of 200–300 bp using 2% low-range ultra agarose (LRUA) gel and amplified via PCR using Phusion High-Fidelity DNA Polymerase for 15 cycles. After quantification using the TBS-380 fluorometer, the libraries were sequenced on an Illumina HiSeq X Ten or NovaSeq 6000 platform to generate 2 × 150 bp paired-end reads.
2.4.3 De novo assembling and annotation
The original paired end reads were treated with trimming and QC by SeqPrep and Sickle. Afterwards, the cleaning data from the samples were adopted to perform de novo assembling with Trinity (). The entire assembling transcription products were retrieved against the National Center for Biotechnology Information (NCBI) protein NR (https://www.ncbi.nlm.nih.gov/public/, July 2023), GO (http://www.geneontology.org/, July 2023), and KEGG (http://www.genome.jp/kegg/, September 2023) data centers via BLASTX (2.9.0) to determine the proteins with the greatest sequential resemblance compared to the aforementioned transcription products to search the functional notations and a representative cut-off E-values<1.0 × 10−5 was set. BLAST2GO software (Agu 2018) was adopted to acquire gene ontology (GO) notations of distinctive assembling transcription products for the description of biology activities, molecule roles and cell constituents. Metabolism pathway assay was completed via the KEGG method ().
2.4.4 Differential expression analysis and functional enrichment
To identify differentially expressed genes (DEGs) between groups, transcript expression levels were calculated using the transcripts per million (TPM) method. Gene abundances were quantified using RSEM (). Differential expression analysis was conducted using either DESeq2 () or DEGseq (). Genes with |log2 fold change| ≥ 1 and false discovery rate (FDR)<0.05 (DESeq2) or FDR<0.001 (DEGseq) were considered significantly differentially expressed. Furthermore, functional enrichment analyses, including Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, were performed to identify significantly enriched GO terms and metabolic pathways. Enrichment was evaluated using a Bonferroni-corrected p-value<0.05 compared with the whole transcriptome background. GO and KEGG analyses were conducted using GOATOOLS (https://pypi.org/project/goatools/, Version 1.4.4) and the Python SciPy library (https://docs.scipy.org/doc/, Version 1.10.0), respectively.
2.5 Quantitative real-time PCR validation
Overall RNA was abstracted via Plant RNA Kit II (OMEGA, Norcross, USA). The first normal cDNA was prepared from overall RNA via a PrimeScript™ RT Reagent Kit with gDNA Eraser (Takara, Dalian, China). Actin was chosen to be an inner control. The heat cycle procedurefor quantitative reverse transcription polymerase chain reaction (qRT-PCR): incipient denaturating at 95°Cfor 2min, denaturating at 94°C for 30 s, annealing at 58°C for 30 s and elongation at 72°Cfor 38s, for an overall 40 cycles. The primers adopted for qRT-PCR herein were presented by Supporting Information (Supplementary Table S2). Each assay was completed three times. The comparative expression level of every unigene was calculated via the 2−ΔΔCt approach.
2.6 Statistical analysis
All experiments were performed with three independent biological replicates per treatment, each consisting of pooled tissues from four roots, and with three technical replicates for each assay. Data are presented as mean±standard error (SE). Statistical analysis was performed using SPSS version 26.0 (SPSS Inc., Chicago, IL, USA). Duncan’s multiple range test was used to determine significant differences among groups, with different uppercase letters indicating significance at p<0.01 and different lowercase letters indicating significance at p<0.05. Graphs were generated using OriginPro 2021 (OriginLab Corporation, Northampton, MA, USA).
3 Results
3.1 Effect of dry method on the chemical and nutritional compositions of C. pilosula
As shown in Figure 1, the levels of active ingredients in C. pilosula changed significantly after a 6-day drying period (p<0.05). The total content of five major active components was significantly higher in the RD and SD samples compared to the FC samples, with increases of 2.01% and 3.06%, respectively (p<0.05; Figure 1A). The adenosine content in FC (38.30 μg/g) was significantly higher than in RD (34.21 μg/g; p<0.05), but not significantly different from SD (36.07 μg/g; p > 0.05; Figure 1B). Similarly, the protocatechuic acid content in FC (92.90 μg/g) was significantly higher than in RD (91.42 μg/g) and SD (92.19 μg/g) (p<0.05; Figure 1C). In contrast, tryptophan content in SD (238.22 μg/g) was significantly higher than in both FC (215.08 μg/g) and RD (216.75 μg/g) (p<0.05; Figure 1D). For syringin, the content in FC (165.45 μg/g) and RD (163.37 μg/g) was significantly higher than in SD (158.46 μg/g) (p<0.05; Figure 1E). Notably, lobetyolin content was highest in RD (141.18 μg/g), which was significantly higher than in FC (122.44 μg/g) and SD (128.67 μg/g) (p<0.01; Figure 1F), with increases of 15.31% and 9.72%, respectively. These findings suggest that different drying treatments exert distinct effects on the accumulation of active components in C. pilosula. Among them, the “Rubbing-sweating” method appears to be particularly effective in promoting lobetyolin accumulation.
Figure 1
The contents of polysaccharides (POL), protein, ash, total polyphenols (TP), total flavonoids (TF), ethanol-soluble extract (ASE), and water-soluble extract (WSE) in the samples are presented in Figure 2 (Supplementary Table S1). The results demonstrated that different drying treatments significantly reduced the contents of protein, total phenolic compounds, ethanol-soluble extracts, and water-soluble extracts in C. pilosula, while significantly increasing the polysaccharide content.
Figure 2
Compared with the FC samples, the protein, ash, TP, ASE, and WSE contents were significantly higher in FC (9.86%, 3.08%, 0.41%, 75.16%, and 76.40%, respectively) than in RD (7.61%, 2.09%, 0.38%, 62.51%, and 65.68%, respectively) and SD (8.94%, 2.91%, 0.33%, 59.59%, and 62.48%, respectively) (p<0.05).
Polysaccharides, which are considered one of the primary bioactive constituents of C. pilosula, were notably increased by both drying treatments. The POL content in RD (19.48%) and SD (15.30%) was significantly higher than in FC (10.56%) (p<0.05). These results indicate that both RD and SD can enhance the accumulation of polysaccharides in C. pilosula, with the RD method being particularly effective.
3.2 Effect of dry method on SOD, POD, CAT, and APX activity
In “Rubbing-sweating” dried C. pilosula tissues, the activities of SOD and CAT reached 126.05 U/g FW and 227.42 U/g/min, respectively, which were significantly higher than those in fresh samples (49.09 U/g FW and 67.35 U/g/min) and shade-dried samples (100.38 U/g FW and 161.44 U/g/min) (p<0.01). The APX activities in “Rubbing-sweating” dried and shade-dried samples were 586.12 U/g FW and 578.46 U/g FW, respectively, both of which were significantly higher than that in fresh samples (333.12 U/g FW) (p<0.01). In addition, POD activity in “Rubbing-sweating” dried tissues reached 1093.83 U/g/min, which was significantly higher than that in fresh (365.07 U/g/min) and shade-dried (856.53 U/g/min) samples (p<0.01). SOD, APX, CAT, and POD are critical antioxidant enzymes involved in reactive oxygen species (ROS) scavenging. The significant enhancement in their activities following the RD treatment suggests that this traditional method induces more intense external stress, likely triggering increased ROS accumulation and, consequently, a stronger antioxidative response in C. pilosula (Figure 3).
Figure 3
3.3 Sequencing and sequence assembly
During the drying process, plant tissues undergo gradual senescence, accompanied by the degradation of proteins and RNA. Therefore, analyzing C. pilosula tissues at the early stages of drying is essential for elucidating the mechanisms underlying the accumulation of active compounds. After 6 days of drying, both RD and SD samples exhibited significant increases in active ingredient content and antioxidant enzyme activities. Consequently, transcriptomic analysis was performed on C. pilosula samples subjected to 6 days of RD and SD treatment to further investigate the molecular mechanisms involved in the accumulation of active compounds during the drying process.
In this study, 44,687,908, 49,800,810, and 47,962,490 clean reads were obtained from the FC, RD, and SD samples, respectively. After removing low-quality sequences, adapters, and ambiguous reads, a total of 399.11 million high-quality clean reads were acquired. These reads were assembled into 194,009 transcripts with an average length of 1,440.03 bp and an N50 length of 2,636 bp. Subsequently, the transcripts were further assembled into 102,614 unigenes, with an average length of 1,033.70 bp and an N50 length of 1,904 bp. For all samples, the mapping rate of clean reads to the reference database ranged from 85.58% to 87.56% (Table 1).
Table 1
| Sample | Raw read | Clean read | Mapped reads | Mapped ratio (%) | Error (%) | Q20 (%) | Q30 (%) |
|---|---|---|---|---|---|---|---|
| FC1 | 43651440 | 43292084 | 18897607 | 87.56 | 0.01 | 98.86 | 96.35 |
| FC2 | 43624858 | 43262026 | 20421843 | 87.23 | 0.01 | 98.84 | 96.29 |
| FC3 | 47892620 | 47509614 | 19204790 | 87.56 | 0.01 | 98.83 | 96.24 |
| RD1 | 45431518 | 45053218 | 18876840 | 87.21 | 0.01 | 98.85 | 96.34 |
| RD2 | 43277238 | 42904240 | 18801078 | 86.92 | 0.01 | 98.85 | 96.32 |
| RD3 | 43599356 | 43231112 | 20551352 | 86.51 | 0.01 | 98.86 | 96.37 |
| SD1 | 43539938 | 43167336 | 19326767 | 85.80 | 0.01 | 98.84 | 96.31 |
| SD2 | 47243318 | 46825102 | 18412532 | 85.83 | 0.01 | 98.79 | 96.15 |
| SD3 | 44250164 | 43867122 | 18499041 | 85.58 | 0.01 | 98.84 | 96.31 |
| Total | 402510450 | 399111854 | |||||
| Type | Unigene | Transcript | |||||
| Total number | 102614 | 194009 | |||||
| Total base | 106071857 | 279379336 | |||||
| Largest length (bp) | 18547 | 18547 | |||||
| Smallest length (bp) | 201 | 201 | |||||
| Average length (bp) | 1033.7 | 1440.03 | |||||
| N50 length (bp) | 1904 | 2636 | |||||
| E90N50 length (bp) | 3487 | 2975 | |||||
| Fragment mapped percent(%) | 61.65 | 87.689 | |||||
| GC percent (%) | 38.84 | 39.53 | |||||
| TransRate score | 0.3118 | 0.44594 | |||||
Summary of transcriptome sequencing data and transcriptome assembly.
Raw reads: total number of sequencing reads before quality control; Clean reads: total number of sequencing reads after quality control; Mapped reads: number of clean reads successfully aligned to the assembled transcripts; Mapped ratio: percentage of paired clean reads mapped to the assembled transcripts; Error (%): average base error rate of quality-controlled data; Q20 (%): percentage of bases with a sequencing quality score higher than 99%; Q30 (%): percentage of bases with a sequencing quality score higher than 99.9%.
The contigs from the nine transcriptome sequencing datasets were integrated and assembled into a total of 102,614 unigenes. Functional annotation of these unigenes was performed using BLAST against seven public databases: GO, KEGG, eggNOG, NR, Swiss-Prot, and Pfam. A total of 30,521, 13,041, 31,191, 38,642, 21,750, and 19,871 unigenes were successfully aligned to these databases, respectively. In total, 39,037 unigenes (38.04% of all unigenes) were annotated in at least one of the functional databases (Table 2). This high annotation rate highlights the overall quality of the sequencing data and the effectiveness of the assembly and annotation process ().
Table 2
| Database | Unigene | Transcript | ||
|---|---|---|---|---|
| Number | Percent (%) | Number | Percent (%) | |
| GO | 30521 | 29.74% | 88372 | 45.55% |
| KEGG | 13041 | 12.71% | 45639 | 23.52% |
| eggNOG | 31191 | 30.40% | 94486 | 48.70% |
| NR | 38642 | 37.66% | 107869 | 55.60% |
| Swiss-Prot | 21750 | 21.20% | 73421 | 37.84% |
| Pfam | 19871 | 19.36% | 66588 | 34.32% |
| Total number of annotated Unigenes/Transcript | 39037 | 38.04% | 108540 | 55.95% |
| Total number of Unigenes/Transcript | 102614 | 100.00% | 194009 | 100.00% |
Numbers of unigenes/transcript annotated using different databases.
3.4 Identification of DEGs and cluster analysis
Principal component analysis (PCA) revealed that the transcriptomes of FC, RD, and SD samples were clearly distinct from one another (Figure 4A). The correlation coefficients (R²) among the transcriptomes ranged from 0.352 to 1.000 (Figure 4B), indicating relatively consistent expression patterns among the nine samples and indirectly confirming the reliability of the sequencing and sampling procedures.
Figure 4
To identify DEGs associated with the different drying treatments, three pairwise comparisons were performed: RD vs. SD, RD vs. FC, and SD vs. FC. DEGs were defined as those with an adjusted p-value ≤ 0.05 and |log2 fold change| ≥ 1. A total of 17,338 DEGs (9,390 upregulated and 7,948 downregulated) were identified in RD vs. SD; 11,007 DEGs (4,500 upregulated and 6,507 downregulated) in RD vs. FC; and 28,138 DEGs (14,963 upregulated and 13,175 downregulated) in SD vs. FC (Figure 4C).
Among the 35,211 unique DEGs identified across all comparisons, 3,779 DEGs were shared among all three comparisons, 16,983 DEGs were shared between two comparisons, and 14,449 DEGs were specific to only one comparison (Figure 4D).
3.5 Functional analysis of DEGs with different dry process
To further elucidate the primary functions of the identified DEGs, Gene Ontology (GO) term annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted. According to the GO enrichment analysis results, DEGs were categorized into three main GO domains: biological processes (BP), cellular components (CC), and molecular functions (MF) (Figure 5). In the RD vs. SD comparison, DEGs were significantly enriched in 7 biological processes, 11 cellular components, and 2 molecular functions (Figure 5A). For the RD vs. FC comparison, DEGs were primarily enriched in 7 biological processes, 5 cellular components, and 8 molecular functions (Figure 5B). In the SD vs. FC comparison, DEGs were enriched in 8 biological processes, 5 cellular components, and 7 molecular functions (Figure 5C). GO enrichment analysis indicated that processes such as the regulation of RNA biosynthetic process, hormone-mediated signaling pathway, and defense response; components including the extracellular region and intracellular membrane-bounded organelles; and functions such as DNA-binding transcription factor activity, oxidoreductase activity, transcription regulator activity, and molecular function inhibitor activity were significantly associated with the DEGs. These findings suggest that these biological pathways and molecular functions may play crucial roles in metabolite transformation and the underlying molecular mechanisms of C. pilosula during different postharvest drying treatments.
Figure 5
The top 20 KEGG pathways identified from the three pairwise comparisons are illustrated in Figure 6. In the RD vs. SD comparison, a total of 3,554 DEGs were mapped to 134 KEGG pathways, with two pathways—DNA replication and cysteine and methionine metabolism—showing significant enrichment (Figure 6A). In the RD vs. FC comparison, 2,418 DEGs were enriched across 137 pathways, among which sixteen pathways exhibited significant enrichment (Figure 6B). These included biosynthesis of various plant secondary metabolites, starch and sucrose metabolism, plant hormone signal transduction, pentose and glucuronate interconversions, monoterpenoid biosynthesis, MAPK signaling pathway, phenylpropanoid biosynthesis, fructose and mannose metabolism, cysteine and methionine metabolism, among others.
Figure 6
For the SD vs. FC comparison, 5,010 DEGs were mapped to 139 pathways, with ten pathways significantly enriched (Figure 6C). Key enriched pathways included plant hormone signal transduction, starch and sucrose metabolism, MAPK signaling pathway, and biosynthesis of various plant secondary metabolites. These findings suggest that metabolic and signaling pathways related to hormone regulation, carbohydrate metabolism, and secondary metabolite biosynthesis play pivotal roles in mediating the effects of traditional postharvest drying methods on the quality and bioactivity of C. pilosula. These pathways likely contribute to the enhanced accumulation of active compounds observed following rubbing-sweating and shade drying treatments.
3.6 qRT-PCR-based verification
To validate the RNA-sequencing (RNA-Seq) results, six DEGs RP-1 (pathogenesis-related protein 1), ETR (Ethylene receptor), GPX (glutathione peroxidase), PDHB (putative pyruvate dehydrogenase), FAB2 (acyl-[acyl-carrier protein] desaturase), SUS (sucrose synthase) related to C. pilosula substance stress response, signal transduction and metabolism were selected for qRT-PCR analysis. The results showed that The mRNA-Seq and RT-PCR data were very closely correlated, and there was high consistency in the up- and down-regulated expression of DEGs. This result supports the reliability of the RNA-Seq analysis, indicating that transcription data were accurate and effective and can be used for the gene expression profile analysis in different drying processes of C. pilosula (Figure 7).
Figure 7
4 Discussion
4.1 Differences in chemical composition and antioxidant enzyme activity of C. pilosula under different drying treatments
Significant differences in the chemical and nutritional compositions of C. pilosula were observed under different drying methods (Figures 1 and 2). Compared with FC, both SD and RD significantly increased the contents of bioactive constituents (p<0.01), with the increase being more pronounced in RD. In particular, the polysaccharide and lobetyolin contents in RD reached 19.50% and 141.18 μg/g, respectively, which were significantly higher than those in FC (10.56% and 122.44 μg/g) and SD (15.30% and 128.67 μg/g) (p<0.01). These results indicate that RD facilitates the accumulation of major active constituents, which may be associated with stress-induced activation of secondary metabolism. This finding is consistent with previous reports that postharvest processing methods markedly improve the quality of other medicinal materials, such as rhubarb and Radix Gentianae (, ; ), further supporting the notion that abiotic stress can enhance the biosynthesis and accumulation of key active compounds in medicinal plants.
Analysis of antioxidant enzyme activities showed that after 6 days of RD and SD treatment, the activities of SOD, CAT, APX, and POD were all significantly elevated (p<0.01), indicating that both drying methods triggered oxidative stress responses in C. pilosula tissues. Notably, the activities of SOD, CAT, and POD in RD were significantly higher than those in SD (p<0.01), suggesting that RD induced a stronger stress stimulus, leading to more pronounced oxidative stress and activation of antioxidant defense systems. This enhanced stress response may be closely related to the more substantial accumulation of secondary metabolites observed under RD.
In summary, rubbing–sweating drying significantly increased the contents of active constituents such as polysaccharides and lobetyolin and markedly enhanced antioxidant enzyme activities. These findings indicate that this drying method improves the medicinal quality of C. pilosula by inducing stress responses and activating relevant metabolic pathways. To systematically elucidate the molecular mechanisms by which different drying methods affect the quality of medicinal materials, transcriptomic analysis was performed on day 6, when changes in active constituents and antioxidant enzyme activities were most pronounced.
4.2 Differential gene expression related to stress response
Reactive oxygen species (ROS) are oxygen-containing molecules with higher chemical reactivity than molecular oxygen and act as double-edged regulators in plants (). Excessive ROS accumulation damages proteins, lipids, and nucleic acids, leading to cell death, whereas controlled ROS production functions as an essential signal to coordinate responses to biotic and abiotic stress (). The major ROS in plant cells—singlet oxygen, superoxide anions, hydrogen peroxide, and hydroxyl radicals—are primarily generated in chloroplasts, mitochondria, and peroxisomes (; ). Antioxidant defenses, including enzymatic (SOD, CAT, POD, APX) and non-enzymatic pathways such as the ascorbate-glutathione (AsA-GSH) cycle, maintain redox homeostasis (; ; ).
Peroxisomes play a particularly important role in regulating oxidative metabolism. Peroxins, which are essential for peroxisome biogenesis and function, maintain peroxisomal integrity and mediate stress responses (; ). In this study, the transcriptional levels of PEX12 and PEX14 were significantly upregulated in RD vs. SD, suggesting enhanced peroxisomal activity under rubbing-sweating treatment. Meanwhile, MPV17, a gene encoding a peroxisomal membrane protein involved in ROS generation (), was upregulated in RD vs. SD but downregulated in SD vs. FC, indicating that RD imposes stronger oxidative stress than SD. Furthermore, genes encoding peroxisomal enzymes such as glyoxylate aminotransferase (AGXT) and hydroxyacid oxidase (HAO), both of which contribute to adaptation under drought and salt stress (; ; ), were significantly upregulated in RD vs. SD, but downregulated in SD vs. FC. In addition, PIPOX, a key enzyme regulating pipecolic acid metabolism and plant immunity (; ), was markedly upregulated in both RD vs. SD and RD vs. FC, while no changes were observed in SD vs. FC. Together, the coordinated upregulation of PEX12, PEX14, MPV17, AGXT, HAO, and PIPOX in RD vs. SD supports the conclusion that rubbing-sweating imposes greater abiotic stress on C. pilosula tissues, triggering peroxisome-related pathways and promoting ROS generation that likely activates secondary metabolism (Figure 8).
Figure 8
Plants employ a range of antioxidant systems to mitigate ROS toxicity. SOD catalyzes the dismutation of superoxide radicals to H2O2, which is further decomposed by CAT and POD (). In this study, genes encoding SOD were upregulated in RD vs. SD but downregulated in SD vs. FC, indicating a stronger oxidative response under RD. POD and CAT showed similar trends, supporting the notion of enhanced ROS detoxification capacity under RD treatment. The AsA-GSH cycle represents a critical pathway for maintaining redox balance in chloroplasts, mitochondria, and cytosol (). Ascorbate peroxidase (APX) plays a complementary role by catalyzing the reduction of H2O2 using ascorbate (AsA) as an electron donor, producing monodehydroascorbate (MDHA) and H2O (; ). Within this cycle, glutathione peroxidase (GPX) reduces H2O2 using GSH, which is regenerated via NADPH-dependent reactions. In this study, GPX-related genes were significantly upregulated in both RD vs. FC and SD vs. FC, suggesting activation under both drying conditions. However, APX genes were particularly upregulated in RD vs. FC, reinforcing the specific contribution of the AsA-GSH cycle to ROS scavenging under rubbing-sweating treatment. Importantly, genes encoding NADPH-generating enzymes such as IDH1, PGD, and G6PD were significantly upregulated in RD vs. SD. These enzymes are not only critical for regenerating reduced glutathione and ascorbate (; ), but also central to the TCA cycle and pentose phosphate pathway, thereby contributing to enhanced energy production and precursor supply. Their upregulation indicates that RD treatment boosts metabolic capacity to support both redox homeostasis and secondary metabolite biosynthesis (Figure 8).
Taken together, these results demonstrate that RD imposes stronger oxidative stress on C. pilosula than SD, as evidenced by elevated ROS generation and coordinated activation of peroxisome-associated metabolic pathways, antioxidant enzymes, and the AsA-GSH cycle. This enhanced stress response not only maintains cellular redox balance but also drives increased metabolic activity and the accumulation of bioactive compounds, thereby improving the medicinal quality of C. pilosula.
4.3 Differential gene expression related to signal transduction
Signal transduction is essential for plants to perceive environmental stimuli and activate appropriate physiological and molecular responses, thereby enhancing stress tolerance (). Among the various pathways, the mitogen-activated protein kinase (MAPK) cascade plays a central role in mediating responses to both biotic and abiotic stresses. In this study, DEGs within the MAPK pathway were mainly associated with defense activation, H2O2 production, and the regulation of ROS homeostasis.
In plant immunity, the receptor FLS2 specifically recognizes pathogen-associated peptides such as flg22, with BAK1 serving as a co-receptor to initiate downstream immune signaling (; ). PR-1, a key defense protein, subsequently accumulates to inhibit pathogen proliferation (). Here, BAK1 and PR-1 were significantly upregulated in RD vs. SD, while FLS2 and PR-1 were enhanced in RD vs. FC, suggesting that mechanical stress during rubbing-sweating (RD) may disrupt cell walls, enhance PAMP (pathogen-associated molecular patterns) perception, and activate immune signaling. NDPK2 and ANP1 further regulate the MAPK cascade. NDPK2 promotes H2O2 accumulation and programmed cell death (), while ANP1, as a MAPKKK, responds to oxidative stress to activate downstream MAPK (). Both were upregulated in RD vs. SD but downregulated in SD vs. FC, indicating that RD strongly induced oxidative signaling and stress-related cell death in C. pilosula tissues. Calcium-mediated signaling also plays a crucial role. Calmodulin (CALM) functions as a Ca²+ sensor, regulating RbohD activity, which transfers electrons from NADPH to oxygen to generate ROS (; ; ). In this study, CALM and RbohD were significantly upregulated in RD vs. SD, implying that RD enhanced Ca²+ influx and ROS production. Meanwhile, MYC2-related genes were downregulated in SD vs. FC, indicating differential activation of jasmonate signaling under shade drying. Together, these results suggest that RD treatment imposed stronger abiotic stress, triggering Ca²+–ROS signaling and potentially regulating secondary metabolism (Figure 8).
Plant hormone signaling further contributes to stress responses. In the abscisic acid (ABA) pathway, SnRK2 mediates responses to drought, salinity, and other stresses (). In this study, SnRK2, MKK3, and MPK1 were significantly upregulated in RD vs. SD, while MKK3 and MPK1 were downregulated in SD vs. FC. As MKK3 activates downstream MAPKs such as MPK1 and MPK7 to regulate disease resistance (; ), these results indicate that RD imposed stronger stress, thereby enhancing ABA- and MAPK-mediated signaling. Ethylene signaling was also activated under RD. RTE1 regulates the receptor ETR1, which controls CTR1 activity to initiate downstream signaling (). Both RTE1 and ETR1 were significantly upregulated in RD vs. SD, suggesting enhanced ethylene signaling. Ethylene mediates responses to diverse stresses, including drought, temperature, salinity, and mechanical damage (). Moreover, ChiB, a chitinase gene related to ethylene-mediated defense, was upregulated in both RD vs. SD and RD vs. FC, further supporting the view that the RD process significantly enhances the defense capacity of C. pilosula tissues (). Jasmonate (JA) signaling also responded strongly to RD. In this pathway, mechanical stress increases JA levels, promoting COI1–JA receptor complex formation, degradation of JAZ repressors, and release of MYC2 to activate defense-related and senescence-associated genes (). In our study, COI1 and MYC2 were upregulated in RD vs. SD, while MYC2 was suppressed in SD vs. FC. These findings suggest that RD activated JA signaling more effectively than SD, enhancing stress resistance but potentially accelerating senescence (Figure 8).
Taken together, transcriptome data demonstrate that rubbing-sweating more strongly activates signal transduction pathways than shade drying. RD treatment upregulated multiple genes in MAPK signaling and hormone pathways, including BAK1, PR-1, NDPK2, ANP1, CALM, RbohD, SnRK2, MKK3, MPK1, RTE1, ETR1, ChiB, COI1, and MYC2. This indicates that RD imposes stronger abiotic stresses—mechanical damage, hypoxia, and elevated temperature—thereby stimulating immune responses, ROS production, and hormone signaling. Such activation likely enhances stress adaptation and secondary metabolite accumulation, contributing to improved medicinal quality of C. pilosula. Nevertheless, the pronounced stress response may also accelerate senescence and programmed cell death, reflecting a trade-off between enhanced quality and tissue longevity under traditional rubbing-sweating processing.
4.4 Differential gene expression related to phenylalanine, tyrosine and tryptophan biosynthesis
C. pilosula, a traditional medicinal and edible herb, is rich in amino acids that contribute significantly to its nutritional and pharmacological value (). Among them, the aromatic amino acids (AAAs)—L-tryptophan (Trp), L-phenylalanine (Phe), and L-tyrosine (Tyr)—are not only indispensable for protein synthesis but also serve as precursors for diverse natural products, influencing plant growth, defense, and stress responses (). In this study, multiple genes involved in AAA biosynthesis were differentially expressed under different drying treatments. Their biosynthesis occurs in two main stages: the conversion of D-erythrose 4-phosphate into chorismate via the shikimate pathway, followed by chorismate-mediated branching into distinct downstream pathways. The shikimate pathway also generates antioxidants such as flavonoids, phenolics, and lignin, which mitigate oxidative stress by scavenging ROS (; ). Within this pathway, aroF, aroB, and aroDE encode key enzymes that catalyze the formation of intermediates including 3-dehydroshikimate and chorismate. Notably, 3-dehydroquinate synthase (aroB) regulates protocatechuic acid biosynthesis, a phenolic compound that reduces ROS levels and enhances stress resilience (; ). In wheat, aroDE is upregulated under drought and salinity stress (). Here, aroF, aroB, and aroDE were significantly upregulated in RD vs. SD, suggesting that rubbing-drying imposes stronger mechanical, hypoxic, and dehydration stresses than SD. This activation likely promotes chorismate synthesis, driving downstream AAA biosynthesis and contributing to secondary metabolite accumulation.
Phe and Tyr are not only structural amino acids but also precursors of diverse polyphenols with antioxidant properties (; ). Both are implicated in activating ROS-scavenging systems and alleviating oxidative damage (). In this study, genes including ADT, hisC, GOT1, and TAT—encoding key enzymes for Phe and Tyr biosynthesis—were significantly upregulated in RD vs. SD, unchanged in RD vs. FC, and downregulated in SD vs. FC. These patterns indicate that RD more effectively enhances Phe/Tyr biosynthesis, potentially improving ROS detoxification and stress tolerance in C. pilosula (Figure 9).
Figure 9
For Trp biosynthesis, anthranilate synthase (trpG and trpE) catalyzes the conversion of chorismate to anthranilate, initiating the pathway. These genes are known to regulate root stress responses and development (). Shade treatment upregulates trpG in walnut, increasing protein accumulation (). Other genes such as trpF, trpC, and trpB also participate in Trp synthesis (). In our study, trpG, trpE, trpF, trpC, and trpB were significantly upregulated in RD vs. SD but downregulated in RD vs. FC, suggesting enhanced transcriptional activation of the Trp pathway under rubbing stress. However, despite this upregulation, Trp content was higher in SD than RD. This discrepancy may reflect accelerated Trp turnover under stress: as a precursor for indole-derived metabolites such as auxin (IAA), serotonin, and alkaloids, Trp is rapidly consumed during stress responses (; ). Thus, RD may promote higher flux through the Trp pathway but reduce net accumulation due to enhanced conversion into secondary metabolites (Figure 9).
Chorismate and its derivatives Phe, Tyr, and Trp are central intermediates in numerous physiological processes, serving as precursors of lignin, flavonoids, and auxins while sustaining metabolic homeostasis and stress adaptation (; ). Our findings show that RD strongly upregulated genes across this pathway, indicating that mechanical and hypoxic stresses activate AAA biosynthesis and promote secondary metabolite production. This enhanced metabolic activity may increase the nutritional and medicinal value of C. pilosula. Overall, the results highlight the dual role of RD-induced stress: while it activates key biosynthetic pathways, leading to enhanced accumulation of secondary metabolites and potential improvement in medicinal quality, it also accelerates metabolic turnover of certain compounds such as Trp. Therefore, the RD process should be carefully optimized—by adjusting kneading cycles and sweating duration—to balance stress-induced quality enhancement with the risk of excessive nutrient loss.
4.5 Differential gene expression related to lobetyolin biosynthesis
Lobetyolin, a characteristic polyacetylene compound in C. pilosula, exhibits diverse pharmacological activities, including antioxidant, anti-inflammatory, and immunomodulatory effects, and is therefore regarded as an important quality marker (; ). However, its biosynthetic pathway remains poorly characterized. Previous studies have suggested that lobetyolin may derive from pyranose-form glucose metabolism () or originate from oleic acid through fatty acid metabolism, with intermediates such as citric acid and linoleic acid serving as potential precursors (). Comparative analyses among C. pilosula varieties further indicate that lobetyolin likely consists of a fatty acid chain conjugated with a glycosyl moiety, reinforcing the role of fatty acid metabolism in its biosynthesis (; ). In the putative pathway, pyruvate is converted into acetyl-CoA by PDHB and pdhC, which then enters fatty acid metabolism. In this study, the genes encoding PDHB and pdhC were significantly upregulated in the RD vs. FC comparison, suggesting enhanced acetyl-CoA supply under rubbing-drying. Acetyl-CoA is subsequently elongated and modified through enzymes such as CEM1 and FATB, leading to palmitic acid formation, which can be further desaturated by FAB2 to oleic acid. Oleic acid is converted by FAD2 to linoleic acid, a key intermediate for polyacetylene biosynthesis (). Here, CEM1 was significantly upregulated in RD vs. FC, while both CEM1 and FATB showed strong induction in SD vs. FC, indicating that fatty acid metabolism is activated by both drying methods, though via slightly different patterns. FAB2 expression was significantly upregulated in both RD vs. FC and SD vs. FC, while FAD2 was higher in RD vs. SD, suggesting that RD more strongly promotes the conversion of oleic acid to linoleic acid. Taken together, these results suggest that lobetyolin biosynthesis in C. pilosula is closely linked to fatty acid metabolism, particularly the conversion of acetyl-CoA through sequential desaturation steps leading to linoleic acid. Both SD and RD enhanced the expression of key genes in this pathway, but RD imposed stronger stress stimuli, thereby exerting a more pronounced promotive effect on lobetyolin accumulation (Figure 9).
In addition to unsaturated fatty acid chains, lobetyolin biosynthesis requires glycosylation mediated by glycosidic compounds (). Within the starch and sucrose metabolism pathway, sucrose is converted into UDP-glucose through the actions of INV, scrK, pgm, UGP2, and SUS, providing the essential glycosyl donor. UDP-glucose is subsequently conjugated to lobetyol by glycosyltransferases (GTs), yielding lobetyolin. In this study, pgm and SUS were significantly upregulated in both RD vs. FC and SD vs. FC, while INV, scrK, and UGP2 showed higher expression in RD vs. SD. These results indicate that RD more strongly stimulates UDP-glucose biosynthesis and related glycosylation processes than SD, thereby facilitating lobetyolin formation (Figure 9).
Overall, lobetyolin biosynthesis in C. pilosula is closely associated with genes from the TCA cycle, fatty acid metabolism, and starch and sucrose metabolism. Transcriptomic analysis demonstrated that drying methods differentially regulate these pathways, with rubbing-drying exerting a stronger promotive effect. These transcriptomic results were consistent with quantitative content determination (Figure 1), confirming that rubbing-drying enhances lobetyolin accumulation in C. pilosula.
4.6 Differential gene expression related to C. pilosula polysaccharide biosynthesis
The biosynthesis of polysaccharides in plants mainly involves starch/sucrose and amino sugar/nucleotide sugar metabolism. Sucrose is first converted into UDP-glucose by sucrose synthase (SUS), which provides precursors for multiple nucleotide sugars. For example, GalE forms UDP-D-galacturonate, UGDH produces UDP-glucuronic acid, AXS generates UDP-D-xylose, RGP forms UDP-L-arabinofuranose, and RHM synthesizes UDP-L-rhamnose. These activated sugars are incorporated into polysaccharide chains by GTs, contributing to cell wall architecture, intracellular signaling, and stress responses (; ; ). Alternatively, sucrose can also be metabolized by SUS and HK into D-mannose-6-phosphate, then converted into GDP-mannose by GMPP, which also functions as a sugar donor in polysaccharide biosynthesis.
In this study, genes encoding RGP and RHM were significantly upregulated in RD vs. FC, while AXS, HK, and GMPP showed an upward trend. Furthermore, galE, UGDH, AXS, RGP, RHM, and GMPP were significantly upregulated in RD vs. SD, suggesting that RD promotes polysaccharide biosynthesis more effectively than SD or FC (Figure 9). Previous studies support this view: UGDH is induced by drought in barley (), AXS enhances oxidative stress resistance in rice (), RGP contributes to stress defense (), and GMPP improves salt tolerance in rice (). These results indicate that RD imposes stronger abiotic stress, including mechanical injury, pathogen exposure, hypoxia, and elevated temperature, thereby stimulating polysaccharide accumulation in C. pilosula.
5 Conclusion
In this study, the differences in key active compounds in C. pilosula under shade drying and rubbing–sweating drying were systematically analyzed. For the first time, the molecular mechanisms underlying the accumulation of these compounds during postharvest processing were explored using transcriptomic approaches. The results demonstrated that both drying methods enhanced the quality of C. pilosula to varying degrees, with rubbing–sweating being more conducive to the accumulation of major active constituents such as lobetyolin and polysaccharides.
Transcriptomic analysis further revealed that rubbing–sweating more strongly activated the expression of genes involved in MAPK signaling and hormone transduction, thereby inducing stress responses, defense mechanisms, and programmed cell death. In addition, peroxisome-related pathways and the antioxidant enzyme system were markedly upregulated, promoting redox homeostasis and the biosynthesis of secondary metabolites. Compared with SD, RD also enhanced the expression of genes related to phenylalanine, tyrosine, tryptophan, lobetyolin, and polysaccharide biosynthesis, suggesting that the enhanced quality results from stress-induced activation of multiple metabolic pathways.
Although this study, for the first time, provides transcriptomic evidence linking the rubbing and sweating treatment with the improvement of C. pilosula quality, further in-depth investigation is still required. Future research will integrate multi-omics data to conduct joint analyses of DEGs and metabolites during the drying process of C. pilosula, as well as perform functional validation of key genes involved in critical pathways. On this basis, systematic optimization of key processing parameters in producing areas (such as rubbing intensity, rubbing frequency, sweating duration, and environmental conditions) will be undertaken to establish technical standards for C. pilosula processing, thereby standardizing the processing workflow. This will provide both scientific evidence and practical guidance to promote the modernization and mechanization of C. pilosula processing.
Statements
Data availability statement
The data presented in the study are deposited in the CNCB-NGDC Genome Sequence Archive (GSA), BioProject accession number PRJCA039682 (https://ngdc.cncb.ac.cn/gsa/browse/CRA025395).
Author contributions
WL: Writing – original draft, Formal Analysis, Funding acquisition, Project administration, Writing – review & editing, Investigation. GB: Software, Investigation, Writing – review & editing, Formal Analysis, Methodology. JS: Writing – original draft, Methodology, Writing – review & editing. WT: Methodology, Writing – review & editing, Formal Analysis. QL: Writing – review & editing. PD: Writing – review & editing, Data curation, Formal Analysis. HW: Formal Analysis, Writing – original draft, Data curation. JC: Formal Analysis, Data curation, Writing – original draft. FG: Formal Analysis, Writing – review & editing, Methodology, Investigation, Conceptualization. YC: Formal Analysis, Writing – review & editing, Methodology, Investigation, Conceptualization.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. This work was supported by Gansu Agricultural University Public Recruitment Doctoral Research Initiation Fund Project (No. GAU-KYQD-2022-07), Gansu Provincial Department of Education: Young Doctoral Support Project (No. 2023-QB-125), Gansu Provincial Youth Science and Technology Fund (No. 23JRRA1444), Lanzhou youth science and technology talent innovation project (No. 2023-QN-162), State Key Laboratory of Aridland Crop Science, Gansu Agricultural University (No. GSCS-2023-13).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2025.1650787/full#supplementary-material
References
1
AghkandA. R.DilmaghaniK. A.GhezelbashZ. D.AsghariB. (2019). Synthesis and biological evaluation of some novel S-β-D-glucosides of 4-amino-5-alkyl-1,2,4-triazole-3-thiones derivatives. Acta Chimica Slovenica66, 344–350. doi: 10.17344/acsi.2018.4841
2
BaiR. B.LiW. Y.LiY. D.MaM.WangY. P.ZhangJ.et al. (2018). Cytotoxicity of two water-soluble polysaccharides from Codonopsis pilosula Nannf. var. modesta (Nannf.) LTShen against human hepatocellular carcinoma HepG2 cells and its mechanism. Int. J. Biol. Macromolecules120, 1544–1550. doi: 10.1016/j.ijbiomac.2018.09.123
3
BaillyC. (2021). Anticancer properties of lobetyolin, an essential component of Radix Codonopsis (Dangshen). Natural products bioprospecting11, 143–153. doi: 10.1007/s13659-020-00283-9
4
BarqawiA. A.AbulfarajA. A. (2023). Salt stress-related mechanisms in leaves of the wild barley hordeum spontaneum generated from RNA-seq datasets. Life-Basel13, 1–17. doi: 10.3390/life13071454
5
CaccamoA.LazzarottoF.Margis-PinheiroM.MessensJ.RemacleC. (2024). The ascorbate peroxidase-related protein: insights into its functioning in Chlamydomonas and Arabidopsis. Front. Plant Sci.15. doi: 10.3389/fpls.2024.1487328
6
CaoM. Y.LiuY. Y.JiangW. M.MengX. X.ZhangW.ChenW. D.et al. (2020). UPLC/MS-based untargeted metabolomics reveals the changes of metabolites profile of Salvia miltiorrhiza bunge during Sweating processing. Sci. Rep.10, 1–10. doi: 10.1038/s41598-020-76650-w
7
CastroB.CittericoM.KimuraS.StevensD. M.WrzaczekM.CoakerG. (2021). Stress-induced reactive oxygen species compartmentalization, perception and signalling. Nat. Plants7, 403–412. doi: 10.1038/s41477-021-00887-0
8
ChenH.HanC.CuiL. M.LiuZ. M.YuF. Y. (2024). Transcriptome analysis of antioxidant system response in Styrax tonkinensis seedlings under flood-drought abrupt alternation. BMC Plant Biol.24, 1–18. doi: 10.1186/s12870-024-05130-4
9
ChenG.HanH. M.YangX. L.DuR. Y.WangX. (2022a). Salt tolerance of rice is enhanced by the SS3 gene, which regulates ascorbic acid synthesis and ROS scavenging. Int. J. Mol. Sci.23, 1–15. doi: 10.3390/ijms231810338
10
ChenY. C.HolmesE. C.RajniakJ.KimJ. G.TangS.FischerC. R.et al. (2018). N-hydroxy-pipecolic acid is a mobile metabolite that induces systemic disease resistance in Arabidopsis. Proc. Natl. Acad. Sci. United States America115, E4920–E4929. doi: 10.1073/pnas.1805291115
11
ChenJ.LiL. F.LinZ. Z.ChengX. L.WeiF.MaS. C. (2023). A quality-comprehensive-evaluation-index-based model for evaluating traditional Chinese medicine quality. Chin. Med.18, 1–14. doi: 10.1186/s13020-023-00782-0
12
ChenJ. C.SuiX. Y.MaB. R.LiY. T.LiN.QiaoL. F.et al. (2022b). Arabidopsis CPR5 plays a role in regulating nucleocytoplasmic transport of mRNAs in ethylene signaling pathway. Plant Cell Rep.41, 1075–1085. doi: 10.1007/s00299-022-02838-1
13
ChoudhuryF. K.RiveroR. M.BlumwaldE.MittlerR. (2017). Reactive oxygen species, abiotic stress and stress combination. Plant J.90, 856–867. doi: 10.1111/tpj.13299
14
CollinA.Daszkowska-GolecA. (2025). Navigating plant stress response: from nucleus to peroxisomes. Trends Plant Sci.30, 815–817. doi: 10.1016/j.tplants.2025.03.016
15
ConesaA.GotzS.Garcia-GomezJ. M.TerolJ.TalonM.RoblesM. (2005). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinf. (Oxford England)21, 3674–3676. doi: 10.1093/bioinformatics/bti610
16
DubiellaU.SeyboldH.DurianG.KomanderE.LassigR.WitteC. P.et al. (2013). Calcium-dependent protein kinase/NADPH oxidase activation circuit is required for rapid defense signal propagation. Proc. Natl. Acad. Sci. U.S.A.110, 8744–8749. doi: 10.1073/pnas.1221294110
17
DugasaM. T.ChalaI. G.WuF. B. (2020). Genotypic difference in secondary metabolism-related enzyme activities and their relative gene expression patterns, osmolyte and plant hormones in wheat. Physiologia Plantarum168, 921–933. doi: 10.1111/ppl.13032
18
El-AzazJ.MaedaH. A. (2025). The multilayered regulation of aromatic amino acid biosynthesis in plants. Trends Biochem. Sci.7, 1–21. doi: 10.1016/j.tibs.2025.07.008
19
GechevT. S.BeninaM.ObataT.TohgeT.SujeethN.MinkovI.et al. (2013). Molecular mechanisms of desiccation tolerance in the resurrection glacial relic Haberlea rhodopensis. Cell. Mol. Life Sci.70, 689–709. doi: 10.1007/s00018-012-1155-6
20
GrabherrM. G.HaasB. J.YassourM.LevinJ. Z.ThompsonD. A.AmitI.et al. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol.29, 644–U130. doi: 10.1038/nbt.1883
21
HaoR. L.LiM. Q.LiF.Sun-WaterhouseD. X.LiD. P. (2022). Protective effects of the phenolic compounds from mung bean hull against H2O2-induced skin aging through alleviating oxidative injury and autophagy in HaCaT cells and HSF cells. Sci. Total Environ.841, 1–11. doi: 10.1016/j.scitotenv.2022.156669
22
HasanuzzamanM.BhuyanM.ZulfiqarF.RazaA.MohsinS. M.Al MahmudJ.et al. (2020). Reactive oxygen species and antioxidant defense in plants under abiotic stress: revisiting the crucial role of a universal defense regulator. Antioxidants9, 1–52. doi: 10.3390/antiox9080681
23
HeW.TaoW. W.ZhangF.JieQ.HeY.ZhuW.et al. (2020). Lobetyolin induces apoptosis of colon cancer cells by inhibiting glutamine metabolism. J. Cell. Mol. Med.24, 3359–3369. doi: 10.1111/jcmm.15009
24
HuJ. P.BakerA.BartelB.LinkaN.MullenR. T.ReumannS.et al. (2012). Plant peroxisomes: biogenesis and function. Plant Cell24, 2279–2303. doi: 10.1105/tpc.112.096586
25
HuangH. L.UllahF.ZhouD. X.YiM.ZhaoY. (2019). Mechanisms of ROS regulation of plant development and stress responses. Front. Plant Sci.10. doi: 10.3389/fpls.2019.00800
26
IydaJ. H.FernandesA.FerreiraF. D.AlvesM. J.PiresT.BarrosL.et al. (2019). Chemical composition and bioactive properties of the wild edible plant Raphanus raphanistrum L. Food Res. Int.121, 714–722. doi: 10.1016/j.foodres.2018.12.046
27
LiS. C. (2023). Novel insight into functions of ascorbate peroxidase in higher plants: More than a simple antioxidant enzyme. Redox Biol.64, 1–6. doi: 10.1016/j.redox.2023.102789
28
LiB.DeweyC. N. (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinf.12, 1–16. doi: 10.1186/1471-2105-12-323
29
LiY. Q.KongD. X.FuY.SussmanM. R.WuH. (2020). The effect of developmental and environmental factors on secondary metabolites in medicinal plants. Plant Physiol. Biochem.148, 80–89. doi: 10.1016/j.plaphy.2020.01.006
30
LiL.LiM.YuL. P.ZhouZ. Y.LiangX. X.LiuZ. X.et al. (2014). The FLS2-associated kinase BIK1 directly phosphorylates the NADPH oxidase rbohD to control plant immunity. Cell Host Microbe15, 329–338. doi: 10.1016/j.chom.2014.02.009
31
LiT. T.ZhangH. Z.XuL. W.ChenX. K.FengJ. S.WuW. J.et al. (2022). StMPK7 phosphorylates and stabilizes a potato RNA-binding protein StUBA2a/b to enhance plant defence responses. Horticulture Res.9, 1–13. doi: 10.1093/hr/uhac177
32
LiangW.ChenY.LiX.GuoF. X.SunJ. C.ZhangX. M.et al. (2021). Label-free proteomic analysis of smoke-drying and shade-drying processes of postharvest rhubarb: A comparative study. Front. Plant Sci.12. doi: 10.3389/fpls.2021.663180
33
LiangM. M.DongQ. L.ZhangX. M.LiuY.LiH.GuoS. P.et al. (2023). Metabolomics and transcriptomics analyses reveals the molecular regulatory mechanisms of walnut (Juglans regia L.) embryos in response to shade treatment. Int. J. Mol. Sci.24, 1–25. doi: 10.3390/ijms241310871
34
LiangW.SunJ. C.BaiG.QiuD. Y.LiQ.DongP. B.et al. (2024). Codonopsis radix: a review of resource utilisation, postharvest processing, quality assessment, and its polysaccharide composition. Front. Pharmacol.15. doi: 10.3389/fphar.2024.1366556
35
LiangW.ZhangW. M.ChenY.GuoF. X.SunJ. C.ZhangX. M.et al. (2022). Accumulation of functional metabolites and transcriptomics in postharvest fume-drying and air-drying process in rhubarb. J. Sci. Food Agric.102, 5628–5641. doi: 10.1002/jsfa.11910
36
LiuG.DuY.ChenM.LuY.LiuL.YangW.et al. (2025). A Colletotrichum-unique effector with the Cx(11)NC motif enhances plant NDPK2 kinase activity to suppress plant immunity. Sci. Adv.11, eadt7970. doi: 10.1126/sciadv.adt7970
37
LjungK. (2013). Auxin metabolism and homeostasis during plant development. Development140, 943–950. doi: 10.1242/dev.086363
38
LoveM. I.HuberW.AndersS. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol.15, 1–21. doi: 10.1186/s13059-014-0550-8
39
LuanF.JiY. F.PengL. X.LiuQ.CaoH. J.YangY.et al. (2021). Extraction, purification, structural characteristics and biological properties of the polysaccharides from Codonopsis pilosula: A review. Carbohydr. Polymers261, 1–24. doi: 10.1016/j.carbpol.2021.117863
40
MaY. T.MaL.XuL.WeiR. A.ChenG. P.DangJ. H.et al. (2024). Research on the mechanism of growth of codonopsis pilosula (Franch.) nannf. Root responding to phenolic stress induced by benzoic acid. Int. J. Mol. Sci.25, 1–21. doi: 10.3390/ijms252011007
41
MaedaH.DudarevaN. (2012). The shikimate pathway and aromatic amino acid biosynthesis in plants. Annu. Rev. Plant Biol.63, 73–105. doi: 10.1146/annurev-arplant-042811-105439
42
MannaM.RengasamyB.SinhaA. K. (2023). Revisiting the role of MAPK signalling pathway in plants and its manipulation for crop improvement. Plant Cell Environ.46, 2277–2295. doi: 10.1111/pce.14606
43
MartiL.SavatinD. V.Gigli-BiscegliaN.De TurrisV.CervoneF.De LorenzoG. (2021). The intracellular ROS accumulation in elicitor-induced immunity requires the multiple organelle-targeted Arabidopsis NPK1-related protein kinases. Plant Cell Environ.44, 931–947. doi: 10.1111/pce.13978
44
MittlerR. (2017). ROS are good. Trends Plant Sci.22, 11–19. doi: 10.1016/j.tplants.2016.08.002
45
MittlerR.ZandalinasS. I.FichmanY.Van BreusegemF. (2022). Reactive oxygen species signalling in plant stress responses. Nat. Rev. Mol. Cell Biol.23, 663–679. doi: 10.1038/s41580-022-00499-2
46
NiL.WangQ. W.ChenC.WangS.ShenT.JiangJ. J.et al. (2022). OsDMI3-mediated OsUXS3 phosphorylation improves oxidative stress tolerance by modulating OsCATB protein abundance in rice. J. Integr. Plant Biol.64, 1087–1101. doi: 10.1111/jipb.13255
47
NiuJ. F.ZhaoG. M.MiZ. Y.ChenL. J.LiuS.WangS. Q.et al. (2020). De novo sequencing of Bletilla striata (Orchidaceae) transcriptome and identification of genes involved in polysaccharide biosynthesis. Genet. Mol. Biol.43, 1–12. doi: 10.1590/1678-4685-gmb-2019-0417
48
Pérez-LlorcaM.PollmannS.MüllerM. (2023). Ethylene and jasmonates signaling network mediating secondary metabolites under abiotic stress. Int. J. Mol. Sci.24, 1–29. doi: 10.3390/ijms24065990
49
QiD. D.ShiY. L.LuM.MaC. Y.DongC. W. (2024). Effect of withering/spreading on the physical and chemical properties of tea: A review. Compr. Rev. Food Sci. Food Saf.23, 1–26. doi: 10.1111/1541-4337.70010
50
RamzanT.ShahbazM.MaqsoodM. F.ZulfiqarU.SamanR. U.LiliN.et al. (2023). Phenylalanine supply alleviates the drought stress in mustard (Brassica campestris) by modulating plant growth, photosynthesis, and antioxidant defense system. Plant Physiol. Biochem.201, 1–14. doi: 10.1016/j.plaphy.2023.107828
51
RenC. Y.ChengT.JiaJ. R.CaoL.ZhangW. J.ZhangS. Z.et al. (2025). Exogenous tryptophan enhances cold resistance of soybean seedlings by promoting melatonin biosynthesis. Physiologia Plantarum177, 1–15. doi: 10.1111/ppl.70189
52
SaqibA.SchellerH. V.FredslundF.WelnerD. H. (2019). Molecular characteristics of plant UDP-arabinopyranose mutases. Glycobiology29, 839–846. doi: 10.1093/glycob/cwz067
53
SeyboldH.TrempelF.RanfS.ScheelD.RomeisT.LeeJ. (2014). Ca2+ signalling in plant immune response: from pattern recognition receptors to Ca2+ decoding mechanisms. New Phytol.204, 782–790. doi: 10.1111/nph.13031
54
ShendeV. V.BaumanK. D.MooreB. S. (2024). The shikimate pathway: gateway to metabolic diversity. Natural Product Rep.41, 604–648. doi: 10.1039/d3np00037k
55
SunJ. C.DuL.LiangW.QuZ.ZhaoH.GuanW. Q. (2023). Effect of postharvest processing on quality traits of Radix Gentianae Macrophyllae: A integrative analysis of metabolomics and proteomics. Plant Physiol. Biochem.204, 1–12. doi: 10.1016/j.plaphy.2023.108099
56
TaiF. J.WangS.LiangB. S.LiY.WuJ. K.FanC. J.et al. (2022). Quaternary ammonium iminofullerenes improve root growth of oxidative-stress maize through ASA-GSH cycle modulating redox homeostasis of roots and ROS-mediated root-hair elongation. J. Nanobiotechnol.20, 1–15. doi: 10.1186/s12951-021-01222-7
57
TohgeT.WatanabeM.HoefgenR.FernieA. R. (2013). Shikimate and phenylalanine biosynthesis in the green lineage. Front. Plant Sci.4. doi: 10.3389/fpls.2013.00062
58
TraverM. S.BradfordS. E.OlmosJ. L.WrightZ. J.MillerM. D.XuW. J.et al. (2022). The structure of the arabidopsis PEX4-PEX22 peroxin complex-insights into ubiquitination at the peroxisomal membrane. Front. Cell Dev. Biol.10. doi: 10.3389/fcell.2022.838923
59
TuT. L.ZhengS. S.RenP. R.MengX. W.ZhaoJ. H.ChenQ.et al. (2021). Coordinated cytokinin signaling and auxin biosynthesis mediates arsenate-induced root growth inhibition. Plant Physiol.185, 1166–1181. doi: 10.1093/plphys/kiaa072
60
VaghelaB.VashiR.RajputK.JoshiR. (2022). Plant chitinases and their role in plant defense: A comprehensive review. Enzyme Microbial. Technol.159, 1–11. doi: 10.1016/j.enzmictec.2022.110055
61
VitámvásP.UrbanM. O.SkodácekZ.KosováK.PitelkováI.VítámvásJ.et al. (2015). Quantitative analysis of proteome extracted from barley crowns grown under different drought conditions. Front. Plant Sci.6. doi: 10.3389/fpls.2015.00479
62
WangL. K.FengZ. X.WangX.WangX. W.ZhangX. G. (2010). DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics26, 136–138. doi: 10.1093/bioinformatics/btp612
63
WangZ. Y.SunB. J. (2009). Annular multiphase flow behavior during deep water drilling and the effect of hydrate phase transition. Petroleum. Sci.6, 57–63. doi: 10.1007/s12182-009-0010-3
64
WangC.TangR. J.KouS. H.XuX. S.LuY.RauscherK.et al. (2024). Mechanisms of calcium homeostasis orchestrate plant growth and immunity. Nature627, 382–388. doi: 10.1038/s41586-024-07100-0
65
WangS. Q.WangB.HuaW. P.NiuJ. F.DangK. K.QiangY.et al. (2017). De novo assembly and analysis of polygonatum sibiricum transcriptome and identification of genes involved in polysaccharide biosynthesis. Int. J. Mol. Sci.18, 1–17. doi: 10.3390/ijms18091950
66
WuQ. H.WeiD.DongL. L.LiuY. P.RenC. X.LiuQ. Q.et al. (2019). Variation in the microbial community contributes to the improvement of the main active compounds of Magnolia officinalis Rehd. et Wils in the process of sweating. Chin. Med.14, 1–15. doi: 10.1186/s13020-019-0267-4
67
XieQ.WangH. X.GuanH. D.XuN.ZhaoX.ChengX. M.et al. (2023). The in vitro/in vivo metabolic pathways analysis of lobetyol, lobetyolin, and lobetyolinin, three polyacetylenes from Codonopsis Radix, by UHPLC-Q/TOF-MS and UHPLC-MS/MS. J. Pharm. Biomed. Anal.223, 1–16. doi: 10.1016/j.jpba.2022.115140
68
XuL.MaL.WeiR. N.MaY. T.MaT. F.DangJ. H.et al. (2024). Effect of continuous cropping on growth and lobetyolin synthesis of the medicinal plant codonopsis pilosula (Franch.) nannf. Based on the integrated analysis of plant-metabolite-soil factors. J. Agric. Food Chem.72, 19604–19617. doi: 10.1021/acs.jafc.4c03664
69
XuH. H.WuM.WeiW. G.RenW. K.ZhengZ. (2025). Chrysanthemum morifolium Ramat. as a traditional tea material: Unraveling the influence of kill-green process on drying characteristics, phytochemical compounds, and volatile profile. Food Res. Int.200, 1–23. doi: 10.1016/j.foodres.2024.115478
70
YangD. D.ChenY.GuoF. X.HuangB. T.OkyereS. A.WangH. (2019). Comparative analysis of chemical composition, antioxidant and antimicrobial activities of leaves, leaf tea and root from Codonopsis pilosula. Ind. Crops Products142, 1–8. doi: 10.1016/j.indcrop.2019.111844
71
YangZ.WangX. Y.HeS.LuoY. H.ChenS.ShanY.et al. (2021). Heat shock treatment maintains the quality attributes of postharvest jujube fruits and delays their senescence process during cold storage. J. Food Biochem.45, 1–16. doi: 10.1111/jfbc.13937
72
YokoyamaR. (2024). Evolution of aromatic amino acid metabolism in plants: a key driving force behind plant chemical diversity in aromatic natural products. Philos. Trans. R. Soc. B-Biological Sci.379, 1–15. doi: 10.1098/rstb.2023.0352
73
YuanM. H.JiangZ. Y.BiG. Z.NomuraK.LiuM. H.WangY. P.et al. (2021). Pattern-recognition receptors are required for NLR-mediated plant immunity. Nature592, 105. doi: 10.1038/s41586-021-03316-6
74
ZeierJ. (2013). New insights into the regulation of plant immunity by amino acid metabolic pathways. Plant Cell Environ.36, 2085–2103. doi: 10.1111/pce.12122
75
ZhangQ. Q.GuoN. N.ZhangY. H.YuY. B.LiuS. Y. (2022b). Genome-wide characterization and expression analysis of pathogenesis-related 1 (PR-1) gene family in tea plant (Camellia sinensis (L.) O. Kuntze) in response to blister-blight disease stress. Int. J. Mol. Sci.23, 1–16. doi: 10.3390/ijms23031292
76
ZhangX. Y.JiangS. N.SunT. T.ZhiW. B.DingK. R.QiaoZ. Y.et al. (2024). Post-harvest processing methods have critical roles on the contents of active metabolites and pharmacological effects of Astragali Radix. Front. Pharmacol.15. doi: 10.3389/fphar.2024.1489777
77
ZhangB. X.LewisJ. A.KovacsF.SattlerS. E.SarathG.KangC. L. (2023). Activity of Cytosolic Ascorbate Peroxidase (APX) from Panicum virgatum against Ascorbate and Phenylpropanoids. Int. J. Mol. Sci.24, 1–17. doi: 10.3390/ijms24021778
78
ZhangQ.LiM. X.YangG. Y.LiuX. Q.YuZ. D.PengS. B. (2022a). Protocatechuic acid, ferulic acid and relevant defense enzymes correlate closely with walnut resistance to Xanthomonas arboricola pv. juglandis. BMC Plant Biol.22, 1–14. doi: 10.1186/s12870-022-03997-9
79
ZhangS. X.ShiY. Y.HuangL. Q.WangC. K.ZhaoD. R.MaK. L.et al. (2020). Comparative transcriptomic analysis of rhizomes, stems, and leaves of Polygonatum odoratum (Mill.) Druce reveals candidate genes associated with polysaccharide synthesis. Gene744, 1–11. doi: 10.1016/j.gene.2020.144626
80
ZhaoL. J.HuangY. X.KellerA. A. (2018). Comparative Metabolic Response between Cucumber (Cucumis sativus) and Corn (Zea mays) to a Cu(OH)2 Nanopesticide. J. Agric. Food Chem.66, 6628–6636. doi: 10.1021/acs.jafc.7b01306
81
ZhaoM. C.ZhangQ.LiuH.TangS.ShangC. Y.ZhangW.et al. (2023). The osmotic stress-activated receptor-like kinase DPY1 mediates SnRK2 kinase activation and drought tolerance in Setaria. Plant Cell35, 3782–3808. doi: 10.1093/plcell/koad200
82
ZhouS. X.ChenM. T.ZhangY. B.GaoQ.NomanA.WangQ.et al. (2019a). OsMKK3, a Stress-Responsive Protein Kinase, Positively Regulates Rice Resistance to Nilaparvata lugens via Phytohormone Dynamics. Int. J. Mol. Sci.20, 1–15. doi: 10.3390/ijms20123023
83
ZhouW.TorresJ. L. L.BlilouI.ZhangX.ZhaiQ.SmantG.et al. (2019b). A Jasmonate signaling network activates root stem cells and promotes regeneration. Cell.177(4), 942–956.10.1016/j.cell.2019.03.006
84
ZhouZ. W.WuQ. Y.YaoZ. L.DengH. L.LiuB. B.YueC.et al. (2020). Dynamics of ADH and related genes responsible for the transformation of C6-aldehydes to C6-alcohols during the postharvest process of oolong tea. Food Sci. Nutr.8, 104–113. doi: 10.1002/fsn3.1272
85
ZouY. F.ChenX. F.MalterudK. E.RiseF.BarsettH.InngjerdingenK. T.et al. (2014). Structural features and complement fixing activity of polysaccharides from Codonopsis pilosula Nannf. var. modesta L.T.Shen roots. Carbohydr. Polymers113, 420–429. doi: 10.1016/j.carbpol.2014.07.036
Summary
Keywords
Codonopsis Radix, rubbing-sweating dried, shade-dried, postharvest processing, stress response genes, polysaccharide and lobetyolin biosynthesis
Citation
Liang W, Bai G, Sun J, Tao W, Li Q, Dong P, Wang H, Cheng J, Guo F and Chen Y (2025) Effects of post-harvest ‘ rubbing-sweating ‘ drying treatment on the accumulation of bioactive compounds in Codonopsis pilosula: a transcriptomic analysis. Front. Plant Sci. 16:1650787. doi: 10.3389/fpls.2025.1650787
Received
23 June 2025
Accepted
12 September 2025
Published
30 September 2025
Volume
16 - 2025
Edited by
Xuming Li, Hugo Biotechnologies Co., Ltd., China
Reviewed by
Annamalai Muthusamy, Manipal Academy of Higher Education, India
Rathna Silviya Lodi, Shandong Academy of Agricultural Sciences, China
Updates
Copyright
© 2025 Liang, Bai, Sun, Tao, Li, Dong, Wang, Cheng, Guo and Chen.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wei Liang, liangw@gsau.edu.cn; Fengxia Guo, guofx@gsau.edu.cn; Yuan Chen, chenyuan@gsau.edu.cn
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.