The chromosome-scale high-quality genome assembly of Panax notoginseng provides insight into dencichine biosynthesis

Genome-Wide Identification and Expression Pattern of the NAC Gene Family in Panax notoginseng

BACKGROUND

The NAC transcription factor family of genes is one of the largest families of transcription factors in plants, playing important functions in plant growth and development, response to adversity stress, disease resistance, and hormone signaling. In this study, we identified the number of members of the Panax notoginseng NAC (PnNAC) gene family and conducted a comprehensive analysis of their physicochemical characteristics, chromosomal location, evolutionary features, and expression patterns both in different parts of the plant at different growth stages and in response to infection by Alternaria panax.

METHODS

The NAC gene family in P. notoginseng was identified using Hidden Markov Model (HMMER) and National Center of Biotechnology Information Conserved Domain Database (NCBI CDD), and their physicochemical properties were analyzed with Perl scripts. Phylogenetic relationships were determined using Clustal Omega and FastTree, and gene structures were visualized with an R script. Promoter regions were analyzed with PlantCARE, motifs with MEME and ggmotif, and transcriptome data were processed using Hical Indexing for Spliced Alignment of Transcripts (HISAT2) and HTseq.

RESULTS

This study identified 98 PnNAC genes in P. notoginseng, analyzed their characteristics (protein lengths 104-882 aa, molecular weights 11.78-100.20 kDa, isoelectric points 4.12-9.75), location (unevenly distributed on 12 chromosomes, no tandem repeats), evolution, and expression patterns (distinct in different parts, growth stages, and after A. panax infection).

CONCLUSIONS

PnNAC plays an important role in the growth and development of P. notoginseng and in its response to A. panax. PnNAC could be a candidate gene for further research on and functional analysis of P. notoginseng disease resistance.

Genome-wide identification and characterization of WRKYs family involved in responses to Cylindrocarpon destructans in Panax notoginseng

BACKGROUND

WRKY transcription factors (TFs) are key regulators of plant responses to biotic and abiotic stresses. Previous studies demonstrated that the role of WRKY TFs play in the disease resistance of Panax notoginseng, the causal agent of root rot disease. However, comprehensive genome-wide analyses of WRKY genes in this species remain scarce.

RESULT

We identified 79 WRKY genes in the P. notoginseng genome, classifing them into three groups based on structural features and phylogenetic relationships: Class I (14 genes), Class II (55 genes), and Class III (10 genes). Of these, 58 PnWRKY genes were mapped to the P. notoginseng chromosomes and showed collinearity with Arabidopsis thaliana, Daucus carota, and three Solanaceae species. Expression analysis revealed that 53 PnWRKY genes were actively transcribed across various tissues, including roots, flowers, stems, rhizomes, and different root parts. Furthermore, PnWRKY genes responded to Cylindrocarpon destructans infection and were induced by jasmonic acid (JA) and salicylic acid (SA). Notably, the ectopic expression of PnWRKY35 in tobacco enhanced resistance to C. destructans, accompanied by increased levels of gibberellins (GA24, GA3, GA8), JA, SA, and brassinolide. KEGG enrichment analysis of 547 differentially expressed genes (DEGs) indicated their involvement pathways related to disease resistance, including the biosynthesis of sesquiterpenes and triterpenes, JA biosynthesis, plant-pathogen interactions and the MAPK signaling pathway.

CONCLUSION

This study demonstrated that the WRKY family in P. notoginseng plays a significant role in resistance to root rot disease. PnWRKY genes are responsive to MeJA and SA induction as well as C. destructans infection. Moreover, ectopic expression of PnWRKY35 activates multiple plant disease resistance pathways, increases phytohormone levels, and enhance resistance to C. destructans. These findings provide a foundation for future exploration of the mechanism underlying P. notoginseng resistance to root rot disease.

Herbal Multiomics Provide Insights into Gene Discovery and Bioproduction of Triterpenoids by Engineered Microbes

Triterpenoids are natural products found in plants that exhibit industrial and agricultural importance. Triterpenoids are typically synthesized through two main pathways: the mevalonate (MVA) and methylerythritol 4-phosphate (MEP) pathways. They then undergo structural diversification with the help of squalene cyclases (OSCs), cytochrome P450 monooxygenases (P450s), UDP glycosyltransferases (UGTs), and acyltransferases (ATs). Advances in multiomics technologies for herbal plants have led to the identification of novel triterpenoid biosynthetic pathways. The application of various analytical techniques facilitates the qualitative and quantitative analysis of triterpenoids. Progress in synthetic biology and metabolic engineering has also facilitated the heterologous production of triterpenoids in microorganisms, such as Escherichia coli and Saccharomyces cerevisiae. This review summarizes recent advances in biotechnological approaches aimed at elucidating the complex pathway of triterpenoid biosynthesis. It also discusses the metabolic engineering strategies employed to increase the level of triterpenoid production in chassis cells.

Transcriptomics and metabolomics revealed the molecular basis of the color formation in the roots of Panax notoginseng

Panax notoginseng is a traditional Chinese medicine rich in many pharmacological components. The root of the 'Miaoxiang Sanqi No. 2' is yellow or greenish yellow, while a novel cultivar-'Wenyuan Ziqi No. 1' shows purple root and is thought to have high medicinal value. Little information is available about the anthocyanin biosynthesis in P. notoginseng root. In this study, we compared the 'Miaoxiang Sanqi No. 2' and 'Wenyuan Ziqi No. 1' in morphological, transcriptional and metabolic levels. The results showed that purple rich in the periderm, rhizome and phloem around cambium of the 'Wenyuan Ziqi No. 1' root and cyanidin 3-O-galactoside was the main anthocyanin causing purple. Moreover, 'Wenyuan Ziqi No. 1' highly accumulated in 155 metabolites, including flavones, phenylpropanoids and lipids. Transcriptome data showed that phenylpropanoid biosynthesis pathway genes are highly expressed in 'Wenyuan Ziqi No. 1'. Conjoint analysis showed that anthocyanin biosynthesis pathway substances were highly accumulated in 'Wenyuan Ziqi No. 1', and the expression level of structural genes involved in anthocyanin biosynthesis pathway was higher in 'Wenyuan Ziqi No. 1'. Meanwhile, eight R2R3-MYB genes that might be involved in anthocyanin biosynthesis were identified. The comprehensive analysis of two cultivars provides new insights into the understanding of root coloration in P. notoginseng.

Transcriptomic Analysis Under Drought and Salt Stress Provides Insight into Genes Putatively Involved in Ginsenoside Biosynthesis in Panax japonicus Meyer

Panax japonicus Meyer, a perennial herb of the dicotyledonaceae family Araliaceae, is a rare folk traditional Chinese medicine, known as "the king of herbal medicine" in China. To understand the genes involved in secondary pathways under drought and salt stress, the transcriptomic analysis of P. japonicus is of vital importance. The transcriptome of underground rhizomes, stems, and leaves under drought and salt stress in P. japonicus were performed using the Illumina HiSeq platform. After de novo assembly of transcripts, expression profiling and identified differentially expressed genes (DEGs) were performed. Furthermore, putative functions of identified DEGs correlated with ginsenoside in P. japonicus were explored using Gene Ontology terms and Kyoto Encyclopedia of Genes and Genome (KEGG) pathway enrichment analysis. A total of 221,804 unigenes were obtained from the transcriptome of P. japonicus. The further analysis revealed that 10,839 unigenes were mapped to 91 KEGG pathways. Furthermore, a total of two metabolic pathways of P. japonicus in response to drought and salt stress related to triterpene saponin synthesis were screened. The sesquiterpene and triterpene metabolic pathways were annotated and finally putatively involved in ginsenoside content and correlation analysis of the expression of these genes were analyzed to identify four genes, β-amyrin synthase, isoprene synthase, squalene epoxidase, and 1-deoxy-D-ketose-5-phosphate synthase, respectively. Our results paves the way for screening highly expressed genes and mining genes related to triterpenoid saponin synthesis. It also provides valuable references for the study of genes involved in ginsenoside biosynthesis and signal pathway of P. japonicus.

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.

Integrated metabolome and transcriptome analysis reveals the regulatory mechanism of low nitrogen-driven biosynthesis of saponins and flavonoids in Panax notoginseng

BACKGROUND

Nitrogen (N) is an important macronutrient involved in the biosynthesis of primary and secondary metabolites in plants. However, the metabolic regulatory mechanism of low-N-induced triterpenoid saponin and flavonoid accumulation in rhizomatous medicinal Panax notoginseng (Burk.) F. H. Chen remains unclear.

METHODS

To explore the potential regulatory mechanism and metabolic basis controlling the response of P. notoginseng to N deficiency, the transcriptome and metabolome were analysed in the roots.

RESULTS

The N content was significantly reduced in roots of N0-treated P. notoginseng (0 kg·N·667 m-2). The C/N ratio was enhanced in the N-deficient P. notoginseng. N deficiency promotes the accumulation of amino acids (L-proline, L-leucine, L-isoleucine, L-norleucine, L-arginine, and L-citrulline) and sugar (arabinose, xylose, glucose, fructose, and mannose), thus providing precursor metabolites for the biosynthesis of flavonoids and triterpenoid saponins. Downregulation of key structural genes (PAL, PAL3, ACC1, CHS2, PPO, CHI3, F3H, DFR, and FGT), in particular with the key genes of F3H, involved in the flavonoid biosynthesis pathway possibly induced the decrease in flavonoid content with increased N supply. Notoginsenoside R1, ginsenoside Re, Rg1, Rd, F1, R1 + Rg1 + Rb1 and total triterpenoid saponins were enhanced in the N0 groups than in the N15 (15 kg·N·667 m-2) plants. Higher phosphoenolpyruvate (an intermediate of glycolyticwith pathway metabolism) and serine (an intermediate of photorespiration) levels induced by N deficiency possibly promote saponin biosynthesis through mevalonic acid (MVA) and methylerythritol (MEP) pathways. Genes (MVD2, HMGS, HMGR1, HMGR2, DXR, and HMGR1) encoding the primary enzymes HMGS, HMGR, DXR, and MVD in the MVA and MEP pathways were significantly upregulated in the N0-treated P. notoginseng. The saponin biosynthesis genes DDS, DDS, CYP716A52, CYP716A47, UGT74AE2, and FPS were upregulated in the N-deficient plants. Upregulation of genes involved in saponin biosynthesis promotes the accumulation of triterpenoid saponins in the N0-grown P. notoginseng.

CONCLUSIONS

N deficiency enhances primary metabolisms, such as amino acids and sugar accumulation, laying the foundation for the synthesis of flavonoids and triterpenoid saponins in P. notoginseng. F3H, DDS, FPS, HMGR, HMGS and UGT74AE2 can be considered as candidates for functional characterisation of the N-regulated accumulation of triterpenoid saponins and flavonoids in future.

Identification of candidate genes and residues for improving nitrogen use efficiency in the N-sensitive medicinal plant Panax notoginseng

BACKGROUND

Nitrogen (N) metabolism-related key genes and conserved amino acid sites in key enzymes play a crucial role in improving N use efficiency (NUE) under N stress. However, it is not clearly known about the molecular mechanism of N deficiency-induced improvement of NUE in the N-sensitive rhizomatous medicinal plant Panax notoginseng (Burk.) F. H. Chen. To explore the potential regulatory mechanism, the transcriptome and proteome were analyzed and the three-dimensional (3D) information and molecular docking models of key genes were compared in the roots of P. notoginseng grown under N regimes.

RESULTS

Total N uptake and the proportion of N distribution to roots were significantly reduced, but the NUE, N use efficiency in biomass production (NUEb), the recovery of N fertilizer (RNF) and the proportion of N distribution to shoot were increased in the N0-treated (without N addition) plants. The expression of N uptake- and transport-related genes NPF1.2, NRT2.4, NPF8.1, NPF4.6, AVP, proteins AMT and NRT2 were obviously up-regulated in the N0-grown plants. Meanwhile, the expression of CIPK23, PLC2, NLP6, TCP20, and BT1 related to the nitrate signal-sensing and transduction were up-regulated under the N0 condition. Glutamine synthetase (GS) activity was decreased in the N-deficient plants, while the activity of glutamate dehydrogenase (GDH) increased. The expression of genes GS1-1 and GDH1, and proteins GDH1 and GDH2 were up-regulated in the N0-grown plants, there was a significantly positive correlation between the expression of protein GDH1 and of gene GDH1. Glu192, Glu199 and Glu400 in PnGS1 and PnGDH1were the key amino acid residues that affect the NUE and lead to the differences in GDH enzyme activity. The 3D structure, docking model, and residues of Solanum tuberosum and P. notoginseng was similar.

CONCLUSIONS

N deficiency might promote the expression of key genes for N uptake (genes NPF8.1, NPF4.6, AMT, AVP and NRT2), transport (NPF1.2 and NRT2.4), assimilation (proteins GS1 and GDH1), signaling and transduction (genes CIPK23, PLC2, NLP6, TCP20, and BT1) to enhance NUE in the rhizomatous species. N deficiency might induce Glu192, Glu199 and Glu400 to improve the biological activity of GS1 and GDH, this has been hypothesized to be the main reason for the enhanced ability of N assimilation in N-deficient rhizomatous species. The key genes and residues involved in improving NUE provide excellent candidates for the breeding of medicinal plants.

Comparative Transcriptome Analysis Reveals the Molecular Mechanism of Bacillus velezensis GJ-7 Assisting Panax notoginseng against Meloidogyne hapla

The rhizosphere bacteria Bacillus velezensis GJ-7, as a biological control agent (BCA), has significant biological control effects on Meloidogyne hapla, and has strong colonization ability in the root of Panax notoginseng. In this study, we conducted a comparative transcriptome analysis using P. notoginseng plant roots treated with B. velezensis GJ-7 or sterile water alone and in combination with M. hapla inoculation to explore the interactions involving the P. notoginseng plant, B. velezensis GJ-7, and M. hapla. Four treatments from P. notoginseng roots were sequenced, and twelve high-quality total clean bases were obtained, ranging from 3.57 to 4.74 Gb. The Gene Ontology (GO) classification and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment showed that numerous DEGs are involved in the phenylpropane biosynthesis pathway and the MAPK signaling pathway in the roots of P. notoginseng with B. velezensis GJ-7 treatments. The analysis results of the two signaling pathways indicated that B. velezensis GJ-7 could enhance the expression of lignin- and camalexin-synthesis-related genes in plant roots to resist M. hapla. In addition, B. velezensis GJ-7 could enhance plant resistance to M. hapla by regulating the expression of resistance-related genes and transcription factors (TFs), including ETR, ERF, ChiB, WRKY22, and PR1. The expression of plant disease resistance genes in the roots of P. notoginseng with different treatments was validated by using real-time quantitative PCR (qRT-PCR), and the results were consistent with transcriptome sequencing. Taken together, this study indicated that B. velezensis GJ-7 can trigger a stronger defense response of P. notoginseng against M. hapla.

The mitochondrial genomes of Panax notoginseng reveal recombination mediated by repeats associated with DNA replication

Panax notoginseng is one of the most valuable medicinal species. However, its mitochondrial genome has not been reported yet. We aimed to determine the mitogenome sequence of P. notoginseng. We de novo assembled the mitogenome with Illumina short reads and Nanopore long reads. The mitochondrial genome of P. notoginseng has a multipartite structure consisting of interconversion between a "master circle" and numerous "subgenomic circles" through recombinations mediated by 64 pairs of repetitive sequences. Among the multipartite structure, seven subgenomic circles were best supported. Six of the seven subgenomic circles shared an 852 bp conserved fragment. The complete mitogenome of P. notoginseng was 662,479 bp long including 34 mitochondrial protein-coding genes (PCGs), three rRNA, and 19 tRNA genes. We identified 166 microsatellite repeats and 26 long-tandem repeats. Phylogenetic analysis resolved a tree that was mostly congruent with the phylogeny of Apiales species described in the APG IV system and the tree built with the chloroplast genome sequences. A total of 12 mitochondrial plastid DNA fragments were identified. Lastly, we predicted 591C-to-U RNA editing sites in the coding regions of mitochondrial PCGs. The mitochondrial genome will lay the foundation for understanding the evolution of Panax species.

Characteristics of heavy metal accumulation and risk assessment in understory Panax notoginseng planting system

Yunnan Province is the main planting area of the precious Chinese herbal medicines (CHM) Panax notoginseng; however, it locates the geological area with high soil heavy metals in China. The frequent land replacement due to continuous cropping obstacles and excessive application of chemicals makes P. notoginseng prone to be contaminated by heavy metals under the farmland P. notoginseng (FPn) planting. To overcome farmland shortage, understory P. notoginseng (UPn) was developed as a new ecological planting model featured by no chemicals input. However, this newly developed planting system requires urgently the soil-plant heavy metal characteristics and risk assessment. This study aimed to evaluate the pollution status of eight heavy metals in the tillage layer (0-20 cm), subsoil layer (20-40 cm) and the plants of UPn in Lancang County, Yunnan Province. Pollution index (Pi) showed that the contamination degree of heavy metals in the tillage layer and subsoil layer was Cd > Pb > Ni > Cu > Zn > Cr > Hg > As and Pb > Cd > Cu > Ni > Cr > Hg > Zn > As, respectively. Potential ecological risk index (PERI) for the tillage layer and subsoil layer was slight and middle, respectively. The exceeding standard rate of Cd, As, Pb, Hg, Cu in the UPn roots was 5.33%, 5.33%, 13.33%, 26.67% and 1.33%, respectively, while only Cd and Hg in the UPn leaves exceeded the standard 10% and 14%, respectively. The enrichment abilities of Cd and Hg in the roots and leaves of UPn were the strongest, while that of Pb was the weakest. The Hazard index (HI) and target hazard quotient (THQ) of eight heavy metals in the roots and leaves of UPn were less than 1.Therefore, our results prove that Upn has no human health risk and provide a scientific basis for the safety evaluation and extension of UPn.

Remodeling metabolism of Corynebacterium glutamicum for high-level dencichine production

Dencichine, a sought-after compound in the medical industry, requires a more efficient and sustainable production method than the current plant extraction process. This study successfully remodeled the metabolic pathway of Corynebacterium glutamicum to produce dencichine from the precursors of L-2,3-diaminopropionate (L-DAP) and oxalyl-coenzyme A. Firstly, a synthetic pathway for L-DAP was established by introducing exogenous enzymes ZmaU/ZmaV. This resulted in a production of 628 mg/L by overexpressing key genes and reducing the endogenous competitive pathway. Secondly, an oxalyl-CoA synthetic pathway was created through the enzymatic conversion of glyoxylate by introducing heterologous enzymes. Finally, with the integration of the exogenous enzyme BAHD, de novo synthesis of dencichine in C. glutamicum was achieved, and production reached 31.75 mg/L within 48-hour fermentation. This achievement represents the first successful biosynthesis of dencichine in C. glutamicum, offering a promising approach for natural product through microbial fermentation.

Integrated analysis of transcriptome, metabolome, and histochemistry reveals the response mechanisms of different ages Panax notoginseng to root-knot nematode infection

Panax notoginseng (P. notoginseng) is an invaluable perennial medicinal herb. However, the roots of P. notoginseng are frequently subjected to severe damage caused by root-knot nematode (RKN) infestation. Although we have observed that P. notoginseng possessed adult-plant resistance (APR) against RKN disease, the defense response mechanisms against RKN disease in different age groups of P. notoginseng remain unexplored. We aimed to elucidate the response mechanisms of P. notoginseng at different stages of development to RKN infection by employing transcriptome, metabolome, and histochemistry analyses. Our findings indicated that distinct age groups of P. notoginseng may activate the phenylpropanoid and flavonoid biosynthesis pathways in varying ways, leading to the synthesis of phenolics, flavonoids, lignin, and anthocyanin pigments as both the response and defense mechanism against RKN attacks. Specifically, one-year-old P. notoginseng exhibited resistance to RKN through the upregulation of 5-O-p-coumaroylquinic acid and key genes involved in monolignol biosynthesis, such as PAL, CCR, CYP73A, CYP98A, POD, and CAD. Moreover, two-year-old P. notoginseng enhanced the resistance by depleting chlorogenic acid and downregulating most genes associated with monolignol biosynthesis, while concurrently increasing cyanidin and ANR in flavonoid biosynthesis. Three-year-old P. notoginseng reinforced its resistance by significantly increasing five phenolic acids related to monolignol biosynthesis, namely p-coumaric acid, chlorogenic acid, 1-O-sinapoyl-D-glucose, coniferyl alcohol, and ferulic acid. Notably, P. notoginseng can establish a lignin barrier that restricted RKN to the infection site. In summary, P. notoginseng exhibited a potential ability to impede the further propagation of RKN through the accumulation or depletion of the compounds relevant to resistance within the phenylpropanoid and flavonoid pathways, as well as the induction of lignification in tissue cells.

Evolutionary roles of polyploidization-derived structural variations in the phenotypic diversification of Panax species

Genomic structural variations (SVs) are widespread in plant and animal genomes and play important roles in phenotypic novelty and species adaptation. Frequent whole genome duplications followed by (re)diploidizations have resulted in high diversity of genome architecture among extant species. In this study, we identified abundant genomic SVs in the Panax genus that are hypothesized to have occurred through during the repeated polyploidizations/(re)diploidizations. Our genome-wide comparisons demonstrated that although these polyploidization-derived SVs have evolved at distinct evolutionary stages, a large number of SV-intersecting genes showed enrichment in functionally important pathways related to secondary metabolites, photosynthesis and basic cellular activities. In line with these observations, our metabolic analyses of these Panax species revealed high diversity of primary and secondary metabolites both at the tissue and interspecific levels. In particular, genomic SVs identified at ginsenoside biosynthesis genes, including copy number variation and large fragment deletion, appear to have played important roles in the evolution and diversification of ginsenosides. A further herbivore deterrence experiment demonstrated that, as major triterpenoidal saponins found exclusively in Panax, ginsenosides provide protection against insect herbivores. Our study provides new insights on how polyploidization-derived SVs have contributed to phenotypic novelty and plant adaptation.

Determination of Dencichine in Panax notoginseng in the Forest and Field Using High Performance Liquid Chromatography

Dencichine is a nonprotein amino acid, an effective ingredient in Panax notoginseng with hemostatic and anti-inflammatory effects. There are few studies on the effects of regions and cultivation models on the accumulation of dencichine. In the current study, the content of dencichine in P. notoginseng collected from its global cultivation and trading center Yunnan, China, (>640 samples) was determined using an optimized high-performance liquid chromatography method coupled with a diode array detector but without derivatization. The recovery rate of this method was 80-110%, the relative standard deviation was <10%, and the limits of detection and quantification were 0.003% (w/w) and 0.01% (w/w), respectively. The content of dencichine in each part of P. notoginseng was as follows: rootlets (39.59%) > main roots (29.91%) > leaves (16.21%) > stems (14.29%). For leaves, P. notoginseng in the forest (5.52 ± 2.26 mg/g) was significantly higher than that in the field (3.93 ± 1.72 mg/g) but opposite for main roots. The origins and altitudes made different contributions to the accumulation of dencichine in P. notoginseng. This study provides an effective analytical method to determine dencichines in various parts of P. notoginseng from different origins and altitudes and supports quality control and product development of P. notoginseng.

Multi-algorithm cooperation research of WRKY genes under nitrogen stress in Panax notoginseng

WRKY transcription factors play an important role in the immune system and the innate defense response of plants. WRKY transcription factors have great feedback on nitrogen stress. In this study, bioinformatics was used to detect the WRKYs of Panax notoginseng (PnWRKYs). The response of PnWRKYs under nitrogen stress was also well studied. PnWRKYs were distributed on 11 chromosomes. According to PnWRKY and Arabidopsis thaliana WRKY (AtWRKY) domains, these PnWRKY proteins were divided into three groups by phylogenetic analysis. MEME analysis showed that almost every member contained motif 1 and motif 2. PlantCARE online predicted the cis-acting elements of the promoter. PnWRKY gene family members obtained 22 pairs of repeat fragments by collinearity analysis. The expression levels of PnWRKYs in different parts (roots, flowers, and leafs) were analyzed by the gene expression pattern. They reflected tissue-specific expressions. The qRT-PCR experiments were used to detect 74 PnWRKYs under nitrogen stress. The results showed that the expression levels of 8 PnWRKYs were significantly induced. The PnWRKY gene family may be involved in biotic/abiotic stresses and hormone induction. This study will not only lay the foundation to explore the functions of PnWRKYs but also provide candidate genes for the future improvement of P. notoginseng.

Reconstruction of engineered yeast factory for high yield production of ginsenosides Rg3 and Rd

Panax notoginseng is one of the most valuable traditional Chinese herbs. The main active ingredients, dammarane-type ginsenosides, show multiple pharmacological activities. Recently, the key UDP-dependent glycosyltransferases (UGTs) involved in the biosynthesis of common ginsenosides have been widely studied. However, only a few UGTs that catalyze ginsenoside formation have been reported. This study further investigated the new catalytic function of 10 characterized UGTs from the public database. PnUGT31(PnUGT94B2) and PnUGT53 (PnUGT71B8)exhibited promiscuous sugar-donor specificity of UDP-glucose and UDP-xylose, which could catalyze the glycosylation of C20-OH sites and elongation of the sugar chain at the C3 and/or C20 sites. We further analyzed the expression patterns in P. notoginseng and predicted the catalytic mechanisms of PnUGT31 and PnUGT53 using molecular docking simulations. Moreover, different gene modules were built to increase the yield of ginsenosides in engineered yeast. The metabolic flow of the proginsenediol (PPD) synthetic pathway was enhanced by LPPDS gene modules based on the engineered strain. The resulting yeast was constructed to produce 1.72 g/L PPD in a shaking flask, but cell growth was significantly inhibited. EGH and LKG gene modules were constructed to achieve high-level production of dammarane-type ginsenosides. The production of G-Rg3 controlled by LKG modules increased 3.84 times (254.07 mg/ L), whereas the G-Rd titer reached 56.68 mg/L after 96 h in shaking flask culture under the control of all modules, both of which yielded the highest values for known microbes.

Comparative genomics reveals the diversification of triterpenoid biosynthesis and origin of ocotillol-type triterpenes in Panax

Gene duplication is assumed to be the major force driving the evolution of metabolite biosynthesis in plants. Freed from functional burdens, duplicated genes can mutate toward novelties until fixed due to selective fitness. However, the extent to which this mechanism has driven the diversification of metabolite biosynthesis remains to be tested. Here we performed comparative genomics analysis and functional characterization to evaluate the impact of gene duplication on the evolution of triterpenoid biosynthesis using Panax species as models. We found that whole-genome duplications (WGDs) occurred independently in Araliaceae and Apiaceae lineages. Comparative genomics revealed the evolutionary trajectories of triterpenoid biosynthesis in plants, which was mainly promoted by WGDs and tandem duplication. Lanosterol synthase (LAS) was likely derived from a tandem duplicate of cycloartenol synthase that predated the emergence of Nymphaeales. Under episodic diversifying selection, the LAS gene duplicates produced by γ whole-genome triplication have given rise to triterpene biosynthesis in core eudicots through neofunctionalization. Moreover, functional characterization revealed that oxidosqualene cyclases (OSCs) responsible for synthesizing dammarane-type triterpenes in Panax species were also capable of producing ocotillol-type triterpenes. Genomic and biochemical evidence suggested that Panax genes encoding the above OSCs originated from the specialization of one OSC gene duplicate produced from a recent WGD shared by Araliaceae (Pg-β). Our results reveal the crucial role of gene duplication in diversification of triterpenoid biosynthesis in plants and provide insight into the origin of ocotillol-type triterpenes in Panax species.

2,3-Butanediol from the leachates of pine needles induces the resistance of Panax notoginseng to the leaf pathogen Alternaria panax

Compared with the use of monocultures in the field, cultivation of medicinal herbs in forests is an effective strategy to alleviate disease. Chemical interactions between herbs and trees play an important role in disease suppression in forests. We evaluated the ability of leachates from needles of Pinus armandii to induce resistance in Panax notoginseng leaves, identified the components via gas chromatography-mass spectrometry (GC-MS), and then deciphered the mechanism of 2,3-Butanediol as the main component in the leachates responsible for resistance induction via RNA sequencing (RNA-seq). Prespraying leachates and 2,3-Butanediol onto leaves could induce the resistance of P. notoginseng to Alternaria panax. The RNA-seq results showed that prespraying 2,3-Butanediol onto leaves with or without A. panax infection upregulated the expression of large number of genes, many of which are involved in transcription factor activity and the mitogen-activated protein kinase (MAPK) signaling pathway. Specifically, 2,3-Butanediol spraying resulted in jasmonic acid (JA) -mediated induced systemic resistance (ISR) by activating MYC2 and ERF1. Moreover, 2,3-Butanediol induced systemic acquired resistance (SAR) by upregulating pattern-triggered immunity (PTI)- and effector-triggered immunity (ETI)-related genes and activated camalexin biosynthesis through activation of WRKY33. Overall, 2,3-Butanediol from the leachates of pine needles could activate the resistance of P. notoginseng to leaf disease infection through ISR, SAR and camalexin biosynthesis. Thus, 2,3-Butanediol is worth developing as a chemical inducer for agricultural production.

Genome-wide identification and expression analysis of MYB gene family under nitrogen stress in Panax notoginseng

The myeloblastosis (MYB) gene family, involved in regulating many important physiological and biochemical processes, is one of the largest transcript factor superfamilies in plants. Since the identification of genome sequencing of Panax notoginseng has been completed, there was little known about the whole genome of its specific MYB gene family and the response to abiotic stresses, in consideration of the excessive application of nitrogen fertilizers in P. notoginseng. In this study, 123 PnMYB genes (MYB genes of P. notoginseng) have been identified and divided into 3 subfamilies by the phylogenetic analysis. These PnMYB genes were unevenly located on 12 chromosomes. Meanwhile, the gene structure and protein conserved domain were established by MEME Suite. The analysis of collinear relationships reflected that there were 121 homologous genes between P. notoginseng and Arabidopsis and 30 between P. notoginseng and rice. Moreover, cis-acting elements of PnMYB gene promoters were predicted which indicated that PnMYBs are involved in biotic, abiotic stress, and hormone induction. The expressions of PnMYB transcription factors in its roots, flowers, and leaves were detected by qRT-PCR and they had tissue-specific expressions and related to the growth of different tissues. Under nitrogen stress, MYB transcription factors had great feedback. Ten R2R3-MYB subfamily genes were significantly induced and indicated the possible function of protecting P. notoginseng from excess nitrogen. With further knowledge on identification of PnMYB gene related to tissue selectivity and abiotic stresses, this study laid the foundation for the functional development of PnMYB gene family and improved the cultivation of P. notoginseng.

Transcriptome analysis reveals that jasmonic acid biosynthesis and signaling is associated with the biosynthesis of asperosaponin VI in Dipsacus asperoides

Dipsacus asperoides is a perennial herb, the roots of which are abundant in asperosaponin VI, which has important medicinal value. However, the molecular mechanism underlying the biosynthesis of asperosaponin VI in D. asperoides remains unclear. In present study, a comprehensive investigation of asperosaponin VI biosynthesis was conducted at the levels of metabolite and transcript during root development. The content of asperosaponin VI was significantly accumulated in two-leaf stage roots, and the spatial distribution of asperosaponin VI was localized in the xylem. The concentration of asperosaponin VI gradually increased in the root with the development process. Transcriptome analysis revealed 3916 unique differentially expressed genes (DEGs) including 146 transcription factors (TFs) during root development in D. asperoides. In addition, α-linolenic acid metabolism, jasmonic acid (JA) biosynthesis, JA signal transduction, sesquiterpenoid and triterpenoid biosynthesis, and terpenoid backbone biosynthesis were prominently enriched. Furthermore, the concentration of JA gradually increased, and genes involved in α-linolenic acid metabolism, JA biosynthesis, and triterpenoid biosynthesis were up-regulated during root development. Moreover, the concentration of asperosaponin VI was increased following methyl jasmonate (MeJA) treatment by activating the expression of genes in the triterpenoid biosynthesis pathway, including acetyl-CoA acetyltransferase (DaAACT), 3-hydroxy-3-methylglutaryl coenzyme A synthase (DaHMGCS), 3-hydroxy-3-methylglutaryl coenzyme-A reductase (DaHMGCR). We speculate that JA biosynthesis and signaling regulates the expression of triterpenoid biosynthetic genes and facilitate the biosynthesis of asperosaponin VI. The results suggest a regulatory network wherein triterpenoids, JA, and TFs co-modulate the biosynthesis of asperosaponin VI in D. asperoides.

From single- to multi-omics: future research trends in medicinal plants

Medicinal plants are the main source of natural metabolites with specialised pharmacological activities and have been widely examined by plant researchers. Numerous omics studies of medicinal plants have been performed to identify molecular markers of species and functional genes controlling key biological traits, as well as to understand biosynthetic pathways of bioactive metabolites and the regulatory mechanisms of environmental responses. Omics technologies have been widely applied to medicinal plants, including as taxonomics, transcriptomics, metabolomics, proteomics, genomics, pangenomics, epigenomics and mutagenomics. However, because of the complex biological regulation network, single omics usually fail to explain the specific biological phenomena. In recent years, reports of integrated multi-omics studies of medicinal plants have increased. Until now, there have few assessments of recent developments and upcoming trends in omics studies of medicinal plants. We highlight recent developments in omics research of medicinal plants, summarise the typical bioinformatics resources available for analysing omics datasets, and discuss related future directions and challenges. This information facilitates further studies of medicinal plants, refinement of current approaches and leads to new ideas.

MicroRNAs in Medicinal Plants

Medicinal plant microRNAs (miRNAs) are an endogenous class of small RNA central to the posttranscriptional regulation of gene expression. Biosynthetic research has shown that the mature miRNAs in medicinal plants can be produced from either the standard messenger RNA splicing mechanism or the pre-ribosomal RNA splicing process. The medicinal plant miRNA function is separated into two levels: (1) the cross-kingdom level, which is the regulation of disease-related genes in animal cells by oral intake, and (2) the intra-kingdom level, which is the participation of metabolism, development, and stress adaptation in homologous or heterologous plants. Increasing research continues to enrich the biosynthesis and function of medicinal plant miRNAs. In this review, peer-reviewed papers on medicinal plant miRNAs published on the Web of Science were discussed, covering a total of 78 species. The feasibility of the emerging role of medicinal plant miRNAs in regulating animal gene function was critically evaluated. Staged progress in intra-kingdom miRNA research has only been found in a few medicinal plants, which may be mainly inhibited by their long growth cycle, high demand for growth environment, immature genetic transformation, and difficult RNA extraction. The present review clarifies the research significance, opportunities, and challenges of medicinal plant miRNAs in drug development and agricultural production. The discussion of the latest results furthers the understanding of medicinal plant miRNAs and helps the rational design of the corresponding miRNA/target genes functional modules.

Multi-algorithm cooperation comprehensive research of bZIP genes under Nitrogen stress in Panax notoginseng

Basic leucine zipper (bZIP) transcription factors play an irreplaceable position in the regulation of plant secondary metabolism, growth and development, and resistance to abiotic stress. Panax notoginseng is a traditional medicinal plant in China, but the systematic identification and the resistance of Panax notoginseng bZIP (PnbZIP) family under nitrogen stress have not been reported before, considering the excessive application of N fertilizers. In this study, we conducted a genome-wide identification of the PnbZIP family and analyzed its phylogeny, tissue selectivity, and abiotic resistence. 74 PnbZIPs were distributed on 12 chromosomes and 8 were not successfully located. Through phylogenetic analysis of Arabidopsis and Panax notoginseng, we divided them into 14 subgroups. In the same subgroup, bZIPs had similiar intron/exon structure and conserved motifs. In the analysis of chromosome structure, two PnbZIP genes were duplicated in tandem on chromosome 3. Intraspecific collinearity analysis showed that 28 PnbZIPs participated in segmental replication. Each PnbZIP promoter contained at least one stress response element or stress-related hormone response element. RNA-seq and qRT-PCR methods were used to analyze the expression patterns of the PnbZIP gene in different tissues (roots, flowers, and leaves) and under different nitrogen stresses. The results showed that the PnbZIP gene had the highest expression level in flowers and reflected tissue-specific expressions. Meanwhile, under the stress of ammonium nitrogen fertilizer and nitrate nitrogen fertilizer, PnbZIPs in roots were differently expressed. 10 PnbZIP stress-responsive genes were screened for significant expression, among which PnbZIP46 was significantly up-regulated, which could be a candidate gene for resistance to Nitrogen stress. This study laid the foundation for functional identification of PnbZIPs and improved the cultivation of Panax notoginseng.

Key Glycosyltransferase Genes of Panax notoginseng: Identification and Engineering Yeast Construction of Rare Ginsenosides

Panax notoginseng is one of the most famous valuable medical plants in China, and its broad application in clinical treatment has an inseparable relationship with the active molecules, ginsenosides. Ginsenosides are glycoside compounds that have varied structures for the diverse sugar chain. Although extensive work has been done, there are still unknown steps in the biosynthetic pathway of ginsenosides. Here, we screened candidate glycosyltransferase genes based on the previous genome and transcriptome data of P. notoginseng and cloned the full length of 27 UGT genes successfully. Among them, we found that PnUGT33 could catalyze different ginsenoside substrates to produce higher polarity rare ginsenosides by extending the sugar chain. We further analyzed the enzymatic kinetics and predicted the catalytic mechanism of PnUGT33 by simulating molecular docking. After that, we reconstructed the biosynthetic pathway of rare ginsenoside Rg₃ and gypenoside LXXV in yeast. By combining the Golden Gate method and overexpressing the UDPG biosynthetic genes, we further improved the yield of engineering yeast strain. Finally, the shake-flask culture yield of Rg₃ reached 51 mg/L and the fed-batch fermentation yield of gypenoside LXXV reached 94.5 mg/L, which was the first and highest record.

Increasing Expression of PnGAP and PnEXPA4 Provides Insights Into the Enlargement of Panax notoginseng Root Size From Qing Dynasty to Cultivation Era

Root size is a key trait in plant cultivation and can be influenced by the cultivation environment. However, physical evidence of root size change in a secular context is scarce due to the difficulty in preserving ancient root samples, and how they were modified during the domestication and cultivation stays unclear. About 100 ancient root samples of Panax notoginseng, preserved as tribute in the Palace Museum (A.D. 1636 to 1912, Qing dynasty), provided an opportunity to investigate the root size changes during the last 100 years of cultivation. The dry weight of ancient root samples (~120 tou samples, tou represents number of roots per 500 g dry weight) is 0.22-fold of the modern samples with the biggest size (20 tou samples). Transcriptome analysis revealed that PnGAP and PnEXPA4 were highly expressed in 20 tou samples, compared with the 120 tou samples, which might contribute to the thicker cell wall and a higher content of lignin, cellulose, and callose in 20 tou samples. A relatively lower content of dencichine and higher content of ginsenoside Rb1 in 20 tou samples are also consistent with higher expression of ginsenoside biosynthesis-related genes. PnPHL8 was filtrated through transcriptome analysis, which could specifically bind the promoters of PnGAP, PnCYP716A47, and PnGGPPS3, respectively. The results in this study represent the first physical evidence of root size changes in P. notoginseng in the last 100 years of cultivation and contribute to a comprehensive understanding of how the cultivation environment affected root size, chemical composition, and clinical application.

Deletion and tandem duplications of biosynthetic genes drive the diversity of triterpenoids in Aralia elata

Araliaceae species produce various classes of triterpene and triterpenoid saponins, such as the oleanane-type triterpenoids in Aralia species and dammarane-type saponins in Panax, valued for their medicinal properties. The lack of genome sequences of Panax relatives has hindered mechanistic insight into the divergence of triterpene saponins in Araliaceae. Here, we report a chromosome-level genome of Aralia elata with a total length of 1.05 Gb. The loss of 12 exons in the dammarenediol synthase (DDS)-encoding gene in A. elata after divergence from Panax might have caused the lack of dammarane-type saponin production, and a complementation assay shows that overexpression of the PgDDS gene from Panax ginseng in callus of A. elata recovers the accumulation of dammarane-type saponins. Tandem duplication events of triterpene biosynthetic genes are common in the A. elata genome, especially for AeCYP72As, AeCSLMs, and AeUGT73s, which function as tailoring enzymes of oleanane-type saponins and aralosides. More than 13 aralosides are de novo synthesized in Saccharomyces cerevisiae by overexpression of these genes in combination. This study sheds light on the diversity of saponins biosynthetic pathway in Araliaceae and will facilitate heterologous bioproduction of aralosides.

The chromosome-level genome assembly of Gentiana dahurica (Gentianaceae) provides insights into gentiopicroside biosynthesis

Gentiana dahurica Fisch. is a perennial herb of the family Gentianaceae. This species is used as a traditional Tibetan medicine because of its rich gentiopicroside constituents. Here, we generate a high-quality, chromosome-level genome of G. dahurica with a total length of 1,416.54 Mb. Comparative genomic analyses showed that G. dahurica shared one whole-genome duplication (WGD) event with Gelsemium sempervirens of the family Gelsemiaceaei and had one additional species-specific WGD after the ancient whole-genome triplication with other eudicots. Further transcriptome analyses identified numerous enzyme coding genes and the transcription factors related to gentiopicroside biosynthesis. A set of candidate cytochrome P450 genes were identified for being involved in biosynthetic shifts from swertiamarin to gentiopicroside. Both gene expressions and the contents measured by high-performance liquid chromatography indicated that the gentiopicrosides were mainly synthesized in the rhizomes with the highest contents. In addition, we found that two above-mentioned WGDs, contributed greatly to the identified candidate genes involving in gentiopicroside biosynthesis. The first reference genome of Gentianaceae we generated here will definitely accelerate evolutionary, ecological, and pharmaceutical studies of this family.

De novo and comparative transcriptomic analysis explain morphological differences in Panax notoginseng taproots

Background

Panax notoginseng (Burk.) F. H. Chen (PN) belonging to the genus Panax of family Araliaceae is widely used in traditional Chinese medicine to treat various diseases. PN taproot, as the most vital organ for the accumulation of bioactive components, presents a variable morphology (oval or long), even within the same environment. However, no related studies have yet explained the molecular mechanism of phenotypic differences. To investigate the cause of differences in the taproot phenotype, de novo and comparative transcriptomic analysis on PN taproot was performed.

Results

A total of 133,730,886 and 114,761,595 paired-end clean reads were obtained based on high-throughput sequencing from oval and long taproot samples, respectively. 121,955 unigenes with contig N50 = 1,774 bp were generated by using the de novo assembly transcriptome, 63,133 annotations were obtained with the BLAST. And then, 42 genes belong to class III peroxidase (PRX) gene family, 8 genes belong to L-Ascorbate peroxidase (APX) gene family, and 55 genes belong to a series of mitogen-activated protein kinase (MAPK) gene family were identified based on integrated annotation results. Differentially expressed genes analysis indicated substantial up-regulation of PnAPX3 and PnPRX45, which are related to reactive oxygen species metabolism, and the PnMPK3 gene, which is related to cell proliferation and plant root development, in long taproots compared with that in oval taproots. Furthermore, the determination results of real-time quantitative PCR, enzyme activity, and H₂O₂ content verified transcriptomic analysis results.

Conclusion

These results collectively demonstrate that reactive oxygen species (ROS) metabolism and the PnMPK3 gene may play vital roles in regulating the taproot phenotype of PN. This study provides further insights into the genetic mechanisms of phenotypic differences in other species of the genus Panax.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08283-w.

PanaxGDB: A Comprehensive Platform for Panax

The genus Panax is a valuable natural medicinal source used worldwide that contains high levels of triterpenoid saponins with extensive pharmacological activities. In past decades, molecular biotechnology and breeding techniques have been respectively used to generate omics data and information on cultivars primarily from Panax ginseng (ginseng), Panax quinquefolium (American ginseng), and Panax notoginseng (Sanqi) to biosynthesize valuable saponins, improve product quality, and conduct cost-controlled cultivation. Although much data have been produced, there are concerns that redundant data might be generated and that relatively scattered data might be overlooked. Therefore, many scientists desire a reliable, comprehensive omics database of the Panax genus that could save time and promote integrated analysis. Therefore, to provide all-inclusive, reliable, and valuable information on the Panax genus, PanaxGDB, an open comprehensive database that integrates data on omics and information on varieties, was established. The database contains information on nearly 600 compounds from 12 Panax species, draft genomic sequences with annotations and gene expression levels, single nucleotide polymorphisms, genome-wide association analysis based on agronomic traits, globally collected germplasm information, summaries, omics data of the Panax genus, and online versatile analytic tools. The Panax genus database will be updated when new data are released to continue serving as a central portal to boost research on the biology and functions of Panax. PanaxGDB is available at: http://panaxGDB.ynau.edu.cn.

Application of High-Throughput Sequencing on the Chinese Herbal Medicine for the Data-Mining of the Bioactive Compounds

The Chinese Herbal Medicine (CHM) has been used worldwide in clinic to treat the vast majority of human diseases, and the healing effect is remarkable. However, the functional components and the corresponding pharmacological mechanism of the herbs are unclear. As one of the main means, the high-throughput sequencing (HTS) technologies have been employed to discover and parse the active ingredients of CHM. Moreover, a tremendous amount of effort is made to uncover the pharmacodynamic genes associated with the synthesis of active substances. Here, based on the genome-assembly and the downstream bioinformatics analysis, we present a comprehensive summary of the application of HTS on CHM for the synthesis pathways of active ingredients from two aspects: active ingredient properties and disease classification, which are important for pharmacological, herb molecular breeding, and synthetic biology studies.

Genomic and Transcriptomic Analysis Provide Insights Into Root Rot Resistance in Panax notoginseng

Panax notoginseng (Panax notoginseng (Burk.) F.H. Chen), a plant of high medicinal value, is severely affected by root rot during cultivation. Here, we generated a reference genome of P. notoginseng, with a contig N50 size of 241.268 kb, and identified 66 disease-resistance genes (R-genes) as candidate genes for breeding disease-resistant varieties. We then investigated the molecular mechanism underlying the responses of resistant and susceptible P. notoginseng genotypes to Fusarium oxysporum infection at six time points by RNA-seq. Functional analysis of the genes differentially expressed between the two genotypes indicated that genes involved in the defense response biological process like hormone transduction and plant-pathogen interaction are continuously and highly expressed in resistant genotype during infection. Moreover, salicylic acid and jasmonic acid levels gradually increased during infection in the resistant genotype. Coexpression analysis showed that PnWRKY22 acts as a hub gene in the defense response of the resistant genotype. Finally, transiently overexpressing PnWRKY22 increased salicylic acid levels in P. notoginseng leaves. Our findings provide a theoretical basis for studying root rot resistance in P. notoginseng.

Chromosomal-scale genome assembly of Eleutherococcus senticosus provides insights into chromosome evolution in Araliaceae

Siberian ginseng (Eleutherococcus senticosus, also known as ciwujia) belongs to the Araliaceae family, which contains more than 1,500 species in 41 genera with diverse chromosome numbers and genome sizes. General consensus posits that ancient whole-genome duplication events and rapid evolutionary radiation are the driving forces for this variation in genome properties. In an attempt to generate more genomic information for the Araliaceae family, we report a 1.30 Gb high-quality draft genome assembly (contig N50 = 309.43 kb) of E. senticosus via PacBio long reads and Hi-C chromatin interaction maps. We found that transposable elements accounted for 72.82% of the genome and a total of 36,372 protein-coding genes were predicted. Comparative analyses of the E. senticosus, Panax notoginseng and Daucus carota genomes revealed a burst expansion of Tekay chromoviral elements in Araliaceae after its divergence with Apiaceae. We also found that E. senticosus underwent a lineage-specific whole-genome duplication event Es-α and a whole-genome duplication event Araliaceae-β that was probably shared by all Araliaceae species. Even though the rediploidization of the E. senticosus genome is evident, pathway analyses show that these two whole-genome duplication events may have contributed to the adaptation of E. senticosus to a cold environment. Taken together, the high-quality genome assembly of E. senticosus provides a valuable genomic resource for future research on the evolution of Araliaceae.

Review on the Development and Applications of Medicinal Plant Genomes

With the development of sequencing technology, the research on medicinal plants is no longer limited to the aspects of chemistry, pharmacology, and pharmacodynamics, but reveals them from the genetic level. As the price of next-generation sequencing technology becomes affordable, and the long-read sequencing technology is established, the medicinal plant genomes with large sizes have been sequenced and assembled more easily. Although the review of plant genomes has been reported several times, there is no review giving a systematic and comprehensive introduction about the development and application of medicinal plant genomes that have been reported until now. Here, we provide a historical perspective on the current situation of genomes in medicinal plant biology, highlight the use of the rapidly developing sequencing technologies, and conduct a comprehensive summary on how the genomes apply to solve the practical problems in medicinal plants, like genomics-assisted herb breeding, evolution history revelation, herbal synthetic biology study, and geoherbal research, which are important for effective utilization, rational use and sustainable protection of medicinal plants.