DNA methylation regulates the secondary metabolism of saponins to improve the adaptability of Eleutherococcus senticosus during drought stress.

Plant growth and development can be significantly impacted by drought stress. Plants will adjust the synthesis and accumulation of secondary metabolites to improve survival in times of water constraint. Simultaneously, drought stress can lead to modifications in the DNA methylation status of plants, and these modifications can directly impact gene expression and product synthesis by changing the DNA methylation status of functional genes involved in secondary metabolite synthesis. However, further research is needed to fully understand the extent to which DNA methylation modifies the content of secondary metabolites to mediate plants' responses to drought stress, as well as the underlying mechanisms involved. Our study found that in Eleutherococcus senticosus (E. senticosus), moderate water deprivation significantly decreased DNA methylation levels throughout the genome and at the promoters of EsFPS, EsSS, and EsSE. Transcription factors like EsMYB-r1, previously inhibited by DNA methylation, can re-bind to the EsFPS promotor region following DNA demethylation. This process promotes gene expression and, ultimately, saponin synthesis and accumulation. The increased saponin levels in E. senticosus acted as antioxidants, enhancing the plant's adaptability to drought stress.

Identification of Eleutherococcus senticosus NAC transcription factors and their mechanisms in mediating DNA methylation of EsFPS, EsSS, and EsSE promoters to regulate saponin synthesis.

BACKGROUND: The formation of pharmacologically active components in medicinal plants is significantly impacted by DNA methylation. However, the exact mechanisms through which DNA methylation regulates secondary metabolism remain incompletely understood. Research in model species has demonstrated that DNA methylation at the transcription factor binding site within functional gene promoters can impact the binding of transcription factors to target DNA, subsequently influencing gene expression. These findings suggest that the interaction between transcription factors and target DNA could be a significant mechanism through which DNA methylation regulates secondary metabolism in medicinal plants. RESULTS: This research conducted a comprehensive analysis of the NAC family in E. senticosus, encompassing genome-wide characterization and functional analysis. A total of 117 EsNAC genes were identified and phylogenetically divided into 15 subfamilies. Tandem duplications and chromosome segment duplications were found to be the primary replication modes of these genes. Motif 2 was identified as the core conserved motif of the genes, and the cis-acting elements, gene structures, and expression patterns of each EsNAC gene were different. EsJUB1, EsNAC047, EsNAC098, and EsNAC005 were significantly associated with the DNA methylation ratio in E. senticosus. These four genes were located in the nucleus or cytoplasm and exhibited transcriptional self-activation activity. DNA methylation in EsFPS, EsSS, and EsSE promoters significantly reduced their activity. The methyl groups added to cytosine directly hindered the binding of the promoters to EsJUB1, EsNAC047, EsNAC098, and EsNAC005 and altered the expression of EsFPS, EsSS, and EsSE genes, eventually leading to changes in saponin synthesis in E. senticosus. CONCLUSIONS: NAC transcription factors that are hindered from binding by methylated DNA are found in E. senticosus. The incapacity of these NACs to bind to the promoter of the methylated saponin synthase gene leads to subsequent alterations in gene expression and saponin synthesis. This research is the initial evidence showcasing the involvement of EsNAC in governing the impact of DNA methylation on saponin production in E. senticosus.

Genome-wide identification of MBD gene family members in Eleutherococcus senticosus with their expression motifs under drought stress and DNA demethylation.

BACKGROUND: Methyl-binding domain (MBD) is a class of methyl-CpG-binding domain proteins that affects the regulation of gene expression through epigenetic modifications. MBD genes are not only inseparable from DNA methylation but have also been identified and validated in various plants. Although MBD is involved in a group of physiological processes and stress regulation in these plants, MBD genes in Eleutherococcus senticosus remain largely unknown. RESULTS: Twenty EsMBD genes were identified in E. senticosus. Among the 24 chromosomes of E. senticosus, EsMBD genes were unevenly distributed on 12 chromosomes, and only one tandem repeat gene existed. Collinearity analysis showed that the fragment duplication was the main motif for EsMBD gene expansion. As the species of Araliaceae evolved, MBD genes also evolved and gradually exhibited different functional differentiation. Furthermore, cis-acting element analysis showed that there were numerous cis-acting elements in the EsMBD promoter region, among which light response elements and anaerobic induction elements were dominant. The expression motif analysis revealed that 60% of the EsMBDs were up-regulated in the 30% water content group. CONCLUSIONS: By comparing the transcriptome data of different saponin contents of E. senticosus and integrating them with the outcomes of molecular docking analysis, we hypothesized that EsMBD2 and EsMBD5 jointly affect the secondary metabolic processes of E. senticosus saponins by binding to methylated CpG under conditions of drought stress. The results of this study laid the foundation for subsequent research on the E. senticosus and MBD genes.

Adaptive Evolution of Chalcone Isomerase Superfamily in Fagaceae.

Chalcone Isomerase (CHI) catalyzes the biosynthesis of flavonoids and secondary metabolism in plants. Currently, there is no systematic analysis of CHIs gene family in Fagaceae which is available. In this study, twenty-two CHI proteins were identified in five species of the Fagaceae family. The CHI superfamily in Fagaceae can be classified into three subfamilies and five groups using phylogenetic analysis, analysis of physicochemical properties, and structural prediction. Results indicated that serine (Ser) and isoleucine (Ile) residues determine the substrate preferred by active Type I Fagaceae CHI, and the chalcone isomerase-like (CHIL) of Fagaceae had active site residues. Adaptive analysis of CHIs showed that CHIs are subject to selection pressure. The active CHI gene of Fagaceae was located in the cytoplasm, and it had the typical gene structure of CHI and contains four exons. All the twenty-two identified CHIs had the conserved domain motif 3, and the different groups had their own structural characteristics. In the process of fatty acid binding protein (FAP) evolution to CHIL and CHI, the physical and chemical properties of proteins also had significant differences in addition to changes in protein functions.

Actinophytocola xanthii sp. nov., an actinomycete isolated from rhizosphere soil of the plant Xanthium sibiricum.

A novel actinomycete strain, 11-183T, was isolated from the rhizosphere soil of Xanthium sibiricum, which was collected in Tangshan, Hebei, China. A phylogenetic analysis based on 16S rRNA gene sequences showed that strain 11-183T formed a clade within the genus Actinophytocola, with a maximum similarity of 98.44 % to Actinophytocola xinjiangensis QAIII60T, followed by 97.76 % similarity to Actinophytocola sediminis YIM M13705T. The average nucleotide identity and digital DNA-DNA hybridization values differed by 79.24 and 23.4 %, respectively, between strain 11-183T and Actinophytocolaxinjiangensis QAIII60T. Strain 11-183T grew well on N-Z-amine agar, and it produced a scant, white aerial mycelium. The isolate formed pale yellow to brown-black colonies and a dense, non-fragmented, branched substrate mycelium, and produced aerial hyphae on which nodular spore chains formed. Growth was observed at salinities ranging from 0 to 2 %, at pH values ranging from pH 6.5 to 8.0 and at temperatures ranging from 15 to 37 °C. The cell-wall amino acids included meso-diaminopimelic acid. Whole cell hydrolysates contained galactose and glucose. The principal fatty acids were iso-C16 : 0, iso-C16 : 1 H and C17 : 1ω6c. Diphosphatidylglycerol, phosphatidylmonomethylethanolamine and phosphatidylethanolamine were the diagnostic phospholipids. The isoprenoid quinones included MK-9(H4) and MK-10(H4). The G+C content of the genomic DNA was 71.7 mol%. Based on the genotypic and phenotypic data, we conclude that strain 11-183T belongs to a novel species of the genus Actinophytocola. The name proposed for the novel species is Actinophytocola xanthii sp. nov., with the type strain 11-183T (=KCTC 39690T= MCCC 1K02062T).

[Analysis of transicriptomes and exploring flavonoid biosynthetic pathway genes in Lithocarpus polystachyus].

The sweet taste and health effect of Lithocarpus polystachyus are mainly related flavonoid. To obtain Lithocarpus transcriptome database and flavonoid biosynthesis-related genes, the RNA-Seq techology (Illumina HiSeq 4000) was used to sequence its transcriptome. Six Gb database was assembled after assembly steps, and 41 043 of L. polystachyus unigenes were obtained. With blasting them with 7 data banks, all unigenes were involved in 51 GO-terms and 237 metabolic pathways. And furthermore 28 genes of the flavonoid biosynthesis-related were found. After using the MicroSatallite, 18 161 SSR were obtained, the single-nucleotide-repeated was the richest at 7 346. These data represent abundant messages about transcripts and provide valuable genome data sources in molecular biology of L. polystachyus.

[Expression profiles analysis of two member of squaleneepoxidase gene family from Eleutherococcus senticosus].

In order to find the characteristics of two members of gene family of squaleneexpoxidase (SE) , a quantitative real time PCR method was developed to analyze the expression of Eleutherococcus senticosus SE1 and SE2 gene from different growth periods and in different organs. The result indicated that all the expression of SE2 more than SE1 in the whole growth period and organs of E. senticosus. And in the whole growth period, expression of SE1 showed a low-high-low characteristic. Both expression of SE2 and growth period showed the same trend. The lowest content of the expression was in the roots. SE1 expression have been improved more than SE2 when treated with MeJA. The expression of E. senticosus SE1 and saponins content had significantly positive correlation (P < 0.05) and the correlation coefficients was 0. 858, while the correlation was not significant for SE2. That indicated that SE1 played a key enzyme gene in the biosynthesis of triterpenoidsaponins

[Establishment of prokaryotic expression and optimization ox expression conditions of Eleutherococcus senticosus P450 gene].

According to the sequence of P450 cDNA of Eleutherococcus senticosus, specific primers were designed. Frokaryotic ex pression vector pET30a-P450 was constructed and the prokaryotic expression conditions were optimized. Results showed that the BL21 after being transformed with the recombinant expression vector accumulated the high amount of recombinant protein. SDS-PAGE analysis showed that the recombinant protein was about 53 kDa. The recombinant accumulated the highest amount of recombinant protein af ter IPTG (1 mmol x L(-1)) at 27-37 degrees C for 24 h. Consequently P450 gene of E. senticosus could be expressed successfully by prokaryotic expression vector pET30a-P450. Induction temperature, IPTG concentration, medium type and amount of induction time could all influence the expression of target protein, but the impact strength was different.

bZIP transcription factor responds to changes in light quality and affects saponins synthesis in Eleutherococcus senticosus.

Light quality considerably influences plant secondary metabolism, yet the precise mechanism underlying its impact on Eleutherococcus senticosus remains elusive. Comprehensive metabolomic and transcriptomic analyses revealed that varying light quality alters the biosynthesis of triterpene saponins by modulating the expression of genes involved in the process in E. senticosus. Through correlation analysis of gene expression and saponin biosynthesis, we identified four light-responsive transcription factors, namely EsbZIP1, EsbZIP2, EsbZIP4, and EsbZIP5. EsbZIP transcription factors function in the nucleus, with light quality-dependent promoter activity. Except for EsbZIP2, the other EsbZIP transcription factors exhibit transcriptional self-activation. Furthermore, EsbZIP can bind to the promoter areas of genes that encode important enzymes (EsFPS, EsSS, and EsSE) involved in triterpene saponin biosynthesis, thereby regulating their expression. Overexpression of EsbZIP resultes in significant down-regulation of most downstream target genes，which leads to a decrease in saponin content. Overall, varying light quality enhances the content of triterpene saponins by suppressing the expression of EsbZIP. This study thus elucidates the molecular mechanism by which E. senticosus adjusts triterpene saponin levels in response to changes in light quality.

[Analysis on codon usage of chloroplast genome of Eleutherococcus senticosus].

OBJECTIVE: To analyze the codon usage of chloroplast genome and the influencing factor in Eleutherococcus senticosus. METHOD: Codon of 52 genes, which were selected from the chloroplast genome sequence of E. senticosus, was multivariate statistical and correspondence analyzed using CodonW and SPSS software. RESULT: GC content at the three position of codons by turns was 46.46%, 38.26%, 29.88%, whereas GC1 and GC2 had a significant correlation coefficient (P < 0.01). The correlation coefficient with GC12, and GC3 was 0.205 and was not significant correlated. There were 30 codons which relative synonymous codon usage was greater than 1 and 29 codons end with A and T. In the corresponding analysis, the first axis shows 10.35% variation. And there was significant correlation coefficient between ENC and GC3. The correlation coefficients with GC3 and ENC were -0.288 and 0.353, respectively. We defined 16 codons from 16 amino acids as the major preference codons in chloroplast genome of E. senticosus. CONCLUSION: The third positions for all codon are preferred to ending with A and T. The codon usage bias is formed under effect of mutation and selection, as well as other factors. But the selection will have a far greater impact than others.

Identification of the AP2/ERF transcription factor family of Eleutherococcus senticosus and its expression correlation with drought stress.

In this study, through analysis of the genome of Eleutherococcus senticosus (ES). 228 AP2/ERF genes were identified and classified into 5 groups AP2 (47 genes), ERF (108 genes), RAV (6 genes), DREB (64 genes), and soloist (3 genes). According to the AP2/ERF classification of Arabidopsis thaliana, the ES AP2/ERF proteins were subdivided into 15 groups. The gene structure and motifs of each group of AP2/ERF in ES were highly similar, which confirmed the conservation of AP2/ERF genes. The ES AP2/ERF genes were unevenly distributed on chromosomes, and a total of four pairs of tandem repeats, and 84 co-linear gene pairs were found, so the AP2/ERF genes expanded in a fragment replication manner, and dominated by pure selection during evolution. By analyzing the transcriptome data of ES under different drought stress conditions, 87 AP2/ERF genes with differential expression were obtained, of which 10 genes with highly significant differences were further analyzed and screened for qRT-PCR validation. To the best of our knowledge, this is the first report on the AP2/ERF gene of Eleutherococcus senticosus, and the bioinformatics analysis and experimental validation provided valuable information about them, which is of great significance for further research on the molecular mechanisms of ES in response to drought stress.

[Molecular cloning of farnesyl diphosphate synthase from Eleutherococcus senticosus and its bioinformatics and expression analysis].

OBJECTIVE: To clone farnesyl diphosphate synthase (FPS) gene from Eleutherococcus senticosus and analyze the bioinformatics and expression of the gene. METHOD: The FPS full length cDNA was cloned by rapid amplification of cDNA ends (RACE). The data was analyzed by bioinformatics method, the structure and function of FPS was deduced. The expression of FPS in different organ of E. senticosus was detected by RT-PCR. RESULT: The full length of FPS cDNA was 1 499 bp containing a 1 029 bp ORF that encoded 342 amino acids. The deduced protein sequence exhibited two Asp riches conserved motifs (DDXXD). Without transmembrane domain, FPS was located in cytoplasm. RT-PCR result showed that FPS gene expressed in different organs of E. senticosus. The expression amounts of FPS in different organs were different significantly (P < 0.05). CONCLUSION: The FPS gene of E. senticosus was successfully cloned for the first time, and provided a stable foundation for studying on its effect and expression control on E. senticosus saponins biosynthesis.

[Cloning of Eleutherococcus senticosus calmodulin gene and effect of endophytic fungus on expression amount of gene].

OBJECTIVE: To clone calmodulin (CaM) gene in Eleutherococcus senticosus, and study the effect of endophytic fungi on expression amount of CaM gene. METHOD: The CaM full length cDNA sequence was cloned by rapid amplification of cDNA ends (RACE). The gene was analyzed and corresponding structure and functions were predicted by the bioinformatics methods. The expression amount of CaM gene affected of endophytic fungus P116-1a, P116-1b, P1094 and P312-1 was detected by RT-PCR. RESULT: The full length of CaM cDNA was 856 bp containing an ORF of 450 bp that encoded a protein of 149 amino acids. The homologous of predicted protein was almost 100% with plants like Panax ginseng and Daucus carota. RT-PCR results showed that endophytic fungus improved CaM expression amount significantly (P<0.05). The highest expression amount of CaM occurred 90 d after reinoculated with endophytic fungi P1094, up to 2.96 times of the control. CONCLUSION: The CaM gene of E. senticosus was successfully cloned for the first time. The results demonstrated that endophytic fungus of E. senticosus improved CaM expression amount significantly.

[Effect of endophytic fungi on expression amount of key enzyme genes in saponins biosynthesis and Eleutherococcus senticosus saponins content].

OBJECTIVE: To analyze the effect of endophytic fungi on expression amount of key enzyme genes SS (squalene synthase gene), SE (squalene epoxidase gene) and bAS (beta-amyrin synthase gene) in saponin biosynthesis and saponins content in Eleutherococcus senticosus. METHOD: Wound method was used for back meeting the endophytic fungi to E. senticosus. With GAPDH as internal control gene, the expression of key enzyme genes was detected by real time PCR method. E. senticosus saponins content was measured by spectrophotometry method. RESULT: When wound method back meeting P116-1a and P116-1b after 30 d, the expression content of SS improved significantly (P < 0.05), however the back meeting of P109-4 and P312-1 didnt change the expression of SS. After that SS expression showed reduction-equality-reduction varying trend. Thirty days after back meeting P312-1, the expression content of SE improved significantly (P < 0.05). Ninty days after back meeting P116-1b and P312-1, the expression content of SE improved significantly to 130%,161%, respectively (P < 0.05). After 120 d, back meeting four endophytic fungi, the expression of SE were significantly higher than the control (P < 0.05). Back meeting four endophytic fungi form 60 d to 120 d, the expression of bAS was significantly higher than the control (P < 0.05). The back meeting four endophytic fungi improved E. senticosus saponins content significantly (P < 0.05). CONCLUSION: Endophytic fungi P116-1a, P116-1b, P1094 and P312-1 significantly effected the expression of key enzyme genes SS, SE and bAS and then affected E. senticosus saponins content. Among the genes, bAS was key target gene.

[Cloning and sequence analysis on cDNA of squalene epoxidase gene in Eleutherococcus senticosus].

OBJECTIVE: To clone and sequence the cDNA of squalene epoxidase gene in Eleutherococcus senticosus. METHOD: Total RNA of E. senticosus was extracted by the improved isothiocyanate method and reverse transcripted into cDNA. The primers were designed depending on the reported SE cDNA sequences of Panax ginseng. The SE cDNAs in E. senticosus was amplified using RT-PCR strategy. RESULT: Sequencing results showed two different cDNA fragments (SE1, SE2) with 1665, 1629 bp each ORF which encoded 554,542 amino acids, respectively. The identities of nucleotides and amino acids between SE1, SE2 were 91.49%, 92.55%. SE1, SE2 had the highest amino acids similarity to the SE1 of P. notoginseng, 93.45%, 94.87% respectively. SE1, SE2 both had a FAD binding domain. The deduced speculated amino acids of SE1, SE2 each had 2,4 membrane-spanning helices. CONCLUSION: The two SE sequences in E. senticosus were firstly separated and reported, which has made foundation for E. senticosus secondary metabolite engineering researches.

[Analyze on the expression differences of key genes in saponins biosynthesize in Eleutherococcus senticosus and its callus].

OBJECTIVE: To analyze the expression differences of SS, SE and bAS in leaf, petiole and callus from leaf and petiole of Eleutherococcus senticosus. METHODS: Calluses were induced from explants of E. senticosus leaf and petiole. The expression amount of SS, SE and bAS gene was detected by real time PCR. RESULTS: The expression amount of SS, SE and bAS genes in callus from leaf were 89.3%, 73.8% and 83.4% of that in leaf, respectively. The expression amount of SS, SE and bAS genes in callus from petiole were 80.2%, 87.7% and 70.9% of that in petiole, respectively. CONCLUSION: The expression amount of SS, SE and bAS gene in callus from E. senticosus are significantly lower than that in plant.

Integrative Analysis of Metabolome and Transcriptome Reveals the Mechanism of Flavonoid Biosynthesis in Lithocarpus polystachyus Rehd.

Lithocarpus polystachyus Rehd has received great attention because of its pharmacological activities, such as inhibiting oxidation and lowering blood glucose and blood pressure, and flavonoids are one of its main pharmacodynamic components. It is important to understand the mechanisms of the flavonoid biosynthetic pathway of L. polystachyus, but the regulation of flavonoid biosynthesis is still unclear. In this study, differentially expressed genes and differentially accumulated metabolites in L. polystachyus were studied by integrating transcriptomics and metabolomics technologies. We confirmed the key genes involved in the flavonoid biosynthesis of L. polystachyus, including LpPAL3, LpCHS1, LpCHS2, LpCHI2, and LpF3H, which had consistent expression patterns with their upstream and downstream metabolites, and there is a significantly positive correlation between them. Compared to mature leaves, stems and young leaves are higher in the expression levels of key structural genes. We deduced that the MYB and bHLH transcription factors regulated the biosynthesis of different flavonoid metabolites and their regulatory patterns. Among them, LpMYB2, LpMYB20, LpMYB54, LpMYB12, and LpWD40-113 positively regulated the biosynthesis of flavones and flavanones. This discovery preliminarily revealed the pathways and key genes of flavonoid biosynthesis in L. polystachyus, which provided a reference for further study on flavonoid biosynthesis.

Transcriptome and Metabolome Analysis of the Synthesis Pathways of Allelochemicals in Eupatorium adenophorum.

Eupatorium adenophorum (Crofton weed) is an invasive weed in more than 30 countries. It inhibits the growth of surrounding plants by releasing allelochemicals during its invasion. However, the synthetic pathways and molecular mechanisms of its allelochemicals have been rarely reported. In this study, the related genes and pathways of allelochemicals in E. adenophorum were analyzed. Transcriptome analysis showed that differentially expressed genes (DEGs) were mainly enriched in the phenylpropanoid biosynthetic pathway and flavonoid biosynthetic pathway. Thirty-three DEGs involved in the synthesis of allelochemicals were identified, and 30 DEGs showed significant differences in blades and stems. Six allelochemicals were identified from blades and stems by ultraperformance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). Correlation analysis of genes and metabolites showed a strong correlation between the five genes and allelochemicals. In addition, this study supplemented the biosynthetic pathway of Eupatorium adenophorum B (HHO). It was found that acyclic sesquiterpene synthase (NES), δ-cadinene synthase (TPS), and cytochrome P450 (P450) were involved in the synthesis of HHO. These findings provide a dynamic spectrum consisting of allelochemical metabolism and a coexpression network of allelochemical synthesis genes in E. adenophorum.

Integrative analysis of transcriptome and metabolome reveals the effect of DNA methylation of chalcone isomerase gene in promoter region on Lithocarpus polystachyus Rehd flavonoids.

Metabolite biosynthesis is regulated by gene expression, which is altered by DNA methylation in the promoter region. Chalcone isomerase (CHI) gene encodes a key enzyme in the Lithocarpus polystachyus Rehd flavonoid pathway, and the expression of L. polystachyus CHI (LpCHI) is closely related to the synthesis of flavonoid metabolites. In this study, we analyzed the DNA methylation site of the LpCHI promoter and its effect on gene expression and metabolite accumulation. The proportions of three types of LpCHI promoter DNA methylation are 7.5%, 68.75%, 18.75%, determined by bisulfite sequencing. Transcriptome sequencing shows that LpCHI is strongly up-regulated in LpCHI promoter methylation Type A but down-regulated in LpCHI promoter methylation Type B and Type C. The expression of LpCHI shows no significant difference between Type B and Type C. Moreover, nine kinds of differentially expressed transcription factors (DETFs) bind to seven CpG-sites of the LpCHI promoter region to regulate LpCHI expression. The results of metabolomics show that differentially accumulated flavonoids are higher in LpCHI promoter methylation Type A than in LpCHI promoter methylation Type B and Type C. Additionally, a positive correlation was found between the LpCHI expression and flavonoids accumulation. These results show that the effect of CpG site-specificity on gene transcription is great than that of overall promoter DNA methylation on gene transcription. The mechanisms of flavonoid genes regulating metabolite accumulation are further revealed.

Metabolome and transcriptome analysis of eleutheroside B biosynthesis pathway in Eleutherococcus senticosus.

Eleutheroside B (syringin) is a medicinal active ingredient extracted from Eleutherococcus senticosus (Ruper. et Maxim.) Maxim with high clinical application value. However, its synthesis pathway remains unknown. Here, we analyzed the eleutheroside B biosynthesis pathway in E. senticosus. Consequently, metabolomic and transcriptomic analyses identified 461 differentially expressed genes (DEGs) and 425 metabolites. Further, we identified 7 DEGs and 67 metabolites involved in the eleutheroside B biosynthetic pathway in the eleutheroside B high and low plants. The correlation between the gene and metabolites was explored using the pearson correlation coefficient (PCC) analysis. Caffeoyl-CoA O-methyltransferase, caffeic acid-O-methyltransferase, ß-amyrin synthase (ß-AS) genes, NAC5, and HB5 transcription factors were identified as candidate genes and transcription factors related to the eleutheroside B synthesis. Eleutheroside B content was the highest at the young stage of the leaves both in the high and low eleutheroside B plants. Quantitative real-time polymerase chain reaction revealed that phenylalanine ammonia-lyase1, cinnamate 4-hydroxylase, ß-AS, and leucoanthocyanidin reductase gene had higher expression levels at the young stage of the leaves in the low eleutheroside B plants but lower expression levels in the high eleutheroside B plants. In the present study, we complemented the eleutheroside B biosynthetic pathway by analyzing the expression levels of relevant genes and metabolite accumulation patterns.