Whole-Genome Sequencing and Analysis of the Chinese Herbal Plant Panax notoginseng

Genomic and metabolomic insights into the antimicrobial compounds and plant growth-promoting potential of Bacillus velezensis B115

The B115 strain, isolated from the inter-root soil of healthy plants in a continuous cropping site of Panax notoginseng, was identified as Bacillus velezensis B115 by 16S rDNA sequence comparison and comparative genomic analysis. B115 is a strain of beneficial microorganisms present in the inter-root zone of plants, with favorable plant growth-promoting properties and antagonistic effects against the plant pathogen Fusarium oxysporum. However, the whole genome of B115 remains unclear, thus restricting its potential applications. To address this gap, the whole genome of B115 has been sequenced and annotated to elucidate the molecular mechanisms underlying its plant growth-promoting and antimicrobial activities. The genome analysis revealed that B115 comprises a single circular chromosome of 4,200,774 bp and a plasmid region 16,878 bp long, possessing a GC content of 45.95%. Moreover, 4349 protein-coding genes were predicted. Notably, the B115 genome contains a substantial number of genes (103) involved in the biosynthesis, transport, and catabolism of secondary metabolites. Through genome mining, 13 BGCs and 540 genes encoding secondary metabolites with predicted roles were identified, including members of the surfactin and fengycin families. Utilizing LC-MS/MS technologies, 2318 metabolites were detected in the fermentation broth of B. velezensis B115, encompassing compounds such as Corynebactin, Gamabufotalin, Pracinostat, Indoleacetic acid, (8)-Gingerol, Luteolin, Liquiritigenin, and other metabolites with antimicrobial, growth-promoting, antioxidant, and antitumor properties. By exploring secondary metabolite-related genes and predicting potential secondary metabolites from the B115 genome based on the whole-genome sequence results, we further elucidate the genomic basis for its ability to promote plant growth and inhibit pathogen activity.

A long road ahead to reliable and complete medicinal plant genomes

Long-read DNA sequencing has propelled medicinal plant genomics forward, with over 400 genomes from 203 plants sequenced by February 2025. However, many genomes still have assembly and annotation flaws, with only 11 gapless telomere-to-telomere assemblies. The core challenge remains identifying genes linked to secondary metabolite biosynthesis, regulation and evolution. High-quality complete genomes are essential for characterizing biosynthetic gene clusters and for enabling robust functional genomics and synthetic biology applications. We propose to focus on achieving more complete genome assemblies in diverse varieties on the basis of refining the currently available ones, leverage lessons from crop genomics research, and apply the cutting-edge genomics technologies in research of medicinal plant genomics.

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.

Unveiling the Contamination Patterns of Neonicotinoid Insecticides: Detection, Distribution, and Risk Assessment in Panax notoginseng across Plant Parts

This study analyzed neonicotinoid insecticides (NEOs) and metabolite (m-NEOs) residues in 136 Panax notoginseng samples via ultra-performance liquid chromatography-tandem mass spectrometry. Imidacloprid was the most detected NEO (88.24% of samples), ranging from 1.50 to 2850 μg/kg. To the best of our knowledge, some novel NEOs were detected in P. notoginseng for the first time. NEO clustering patterns varied among plant parts, with higher contamination in leaves and flowers. Fourteen NEO/m-NEOs, including cycloxaprid and acetamiprid, showed site-specific behavior, indicating the possibility of using multiple NEOs simultaneously during planting, resulting in formation of distinct metabolites in different plant parts. Transfer rates in decoction and infusion ranged from 10.06 to 32.33%, reducing residues postprocessing. Dietary risk assessment showed low hazard quotients (HQa: 7.05 × 10-7 to 2.09 × 10-2; HQc: 3.74 × 10-7 to 2.38 × 10-3), but risk-ranking scores indicated potential hazards with imidacloprid and acetamiprid in flowers and leaves. The findings are expected to promote safety assessment and distribution research of NEOs in plants.

Comparative Metabolome and Transcriptome Analyses Reveal Differential Enrichment of Metabolites with Age in Panax notoginseng Roots

Panax notoginseng is a perennial plant well known for its versatile medicinal properties, including hepatoprotective, antioxidant, anti-inflammatory, anti-tumor, estrogen-like, and antidepressant characteristics. It has been reported that plant age affects the quality of P. notoginseng. This study aimed to explore the differential metabolome and transcriptome of 2-year (PN2) and 3-year-old (PN3) P. notoginseng plant root samples. Principal component analysis of metabolome and transcriptome data revealed major differences between the two groups (PN2 vs. PN3). A total of 1813 metabolites and 28,587 genes were detected in this study, of which 255 metabolites and 3141 genes were found to be differential (p < 0.05) between PN2 vs. PN3, respectively. Among differential metabolites and genes, 155 metabolites and 1217 genes were up-regulated, while 100 metabolites and 1924 genes were down-regulated. The KEGG pathway analysis revealed differentially enriched metabolites belonging to class lipids ("13S-hydroperoxy-9Z, 11E-octadecadionic acid", "9S-hydroxy-10E, 12Z-octadecadionic acid", "9S-oxo-10E, 12Z-octadecadionic acid", and "9,10,13-trihydroxy-11-octadecadionic acid"), nucleotides and derivatives (guanine and cytidine), and phenolic acids (chlorogenic acid) were found to be enriched (p < 0.05) in PN3 compared to PN2. Further, these differentially enriched metabolites were found to be significantly (p < 0.05) regulated via linoleic acid metabolism, nucleotide metabolism, plant hormone signal transduction, and arachidonic acid metabolism pathways. Furthermore, the transcriptome analysis showed the up-regulation of key genes MAT, DMAS, SDH, gallate 1-beta-glucosyltransferase, and beta-D-glucosidase in various plants' secondary metabolic pathways and SAUR, GID1, PP2C, ETR, CTR1, EBF1/2, and ERF1/2 genes observed in phytohormone signal transduction pathway that is involved in plant growth and development, and protection against the various stressors. This study concluded that the roots of a 3-year-old P. notoginseng plant have better metabolome and transcriptome profiles compared to a 2-year-old plant with importantly enriched metabolites and genes in pathways related to metabolism, plant hormone signal transduction, and various biological processes. These findings provide insights into the plant's dynamic biochemical and molecular changes during its growth that have several implications regarding its therapeutic use.

Analysis and Identification of Genes Associated with the Desiccation Sensitivity of Panax notoginseng Seeds

Panax notoginseng (Burk.) F. H. Chen, a species of the genus Panax, radix has been traditionally used to deal with various hematological diseases and cardiovascular diseases since ancient times in East Asia. P. notoginseng produces recalcitrant seeds which are sensitive to desiccation and difficult to store for a long time. However, few data are available on the mechanism of the desiccation sensitivity of P. notoginseng seeds. To gain a comprehensive perspective of the genes associated with desiccation sensitivity, cDNA libraries from seeds under control and desiccation processes were prepared independently for Illumina sequencing. The data generated a total of 70,189,896 reads that were integrated and assembled into 55,097 unigenes with a mean length of 783 bp. In total, 12,025 differentially expressed genes (DEGs) were identified during the desiccation process. Among these DEGs, a number of central metabolism, hormonal network-, fatty acid-, and ascorbate-glutathione-related genes were included. Our data provide a comprehensive resource for identifying the genes associated with the desiccation sensitivity of P. notoginseng seeds.

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.

Genome-wide identification and characterization of members of the LEA gene family in Panax notoginseng and their transcriptional responses to dehydration of recalcitrant seeds

BACKGROUND

Late embryogenesis abundant (LEA) proteins play an important role in dehydration process of seed maturation. The seeds of Panax notoginseng (Burkill) F. H. Chen are typically characterized with the recalcitrance and are highly sensitive to dehydration. However, it is not very well known about the role of LEA proteins in response to dehydration stress in P. notoginseng seeds. We will perform a genome-wide analysis of the LEA gene family and their transcriptional responses to dehydration stress in recalcitrant P. notoginseng seeds.

RESULTS

In this study, 61 LEA genes were identified from the P. notoginseng genome, and they were renamed as PnoLEA. The PnoLEA genes were classified into seven subfamilies based on the phylogenetic relationships, gene structure and conserved domains. The PnoLEA genes family showed relatively few introns and was highly conserved. Unexpectedly, the LEA_6 subfamily was not found, and the LEA_2 subfamily contained 46 (75.4%) members. Within 19 pairs of fragment duplication events, among them 17 pairs were LEA_2 subfamily. In addition, the expression of the PnoLEA genes was obviously induced under dehydration stress, but the germination rate of P. notoginseng seeds decreased as the dehydration time prolonged.

CONCLUSIONS

We found that the lack of the LEA_6 subfamily, the expansion of the LEA_2 subfamily and low transcriptional levels of most PnoLEA genes might be implicated in the recalcitrant formation of P. notoginseng seeds. LEA proteins are essential in the response to dehydration stress in recalcitrant seeds, but the protective effect of LEA protein is not efficient. These results could improve our understanding of the function of LEA proteins in the response of dehydration stress and their contributions to the formation of seed recalcitrance.

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.

PlantNLRatlas: a comprehensive dataset of full- and partial-length NLR resistance genes across 100 chromosome-level plant genomes

Plants have evolved two layers of protection against biotic stress: PAMP-triggered immunity (PTI) and effector-triggered immunity (ETI). The primary mechanism of ETI involves nucleotide-binding leucine-rich repeat immune receptors (NLRs). Although NLR genes have been studied in several plant species, a comprehensive database of NLRs across a diverse array of species is still lacking. Here, we present a thorough analysis of NLR genes across 100 high-quality plant genomes (PlantNLRatlas). The PlantNLRatlas includes a total of 68,452 NLRs, of which 3,689 are full-length and 64,763 are partial-length NLRs. The majority of NLR groups were phyletically clustered. In addition, the domain sequences were found to be highly conserved within each NLR group. Our PlantNLRatlas dataset is complementary to RefPlantNLR, a collection of NLR genes which have been experimentally confirmed. The PlantNLRatlas should prove helpful for comparative investigations of NLRs across a range of plant groups, including understudied taxa. Finally, the PlantNLRatlas resource is intended to help the field move past a monolithic understanding of NLR structure and function.

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.

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.

Panax notoginseng WRKY Transcription Factor 9 Is a Positive Regulator in Responding to Root Rot Pathogen Fusarium solani

Panax notoginseng (Burk) F.H. Chen is a rare and valuable Chinese herb, but root rot mainly caused by Fusarium solani severely affects the yield and quality of P. notoginseng herbal materials. In this study, we isolated 30 P. notoginseng WRKY transcription factors (TFs), which were divided into three groups (I, II, and III) on the basis of a phylogenetic analysis. The expression levels of 10 WRKY genes, including PnWRKY9, in P. notoginseng roots increased in response to a methyl jasmonate (MeJA) treatment and the following F. solani infection. Additionally, PnWRKY9 was functionally characterized. The PnWRKY9 protein was localized to the nucleus. The overexpression of PnWRKY9 in tobacco (Nicotiana tabacum) considerably increased the resistance to F. solani, whereas an RNAi-mediated decrease in the PnWRKY9 expression level in P. notoginseng leaves increased the susceptibility to F. solani. The RNA sequencing and hormone content analyses of PnWRKY9-overexpression tobacco revealed that PnWRKY9 and the jasmonic acid (JA) signaling pathway synergistically enhance disease resistance. The PnWRKY9 recombinant protein was observed to bind specifically to the W-box sequence in the promoter of a JA-responsive and F. solani resistance-related defensin gene (PnDEFL1). A yeast one-hybrid assay indicated that PnWRKY9 can activate the transcription of PnDEFL1. Furthermore, a co-expression assay in tobacco using β-glucuronidase (GUS) as a reporter further verified that PnWRKY9 positively regulates PnDEFL1 expression. Overall, in this study, we identified P. notoginseng WRKY TFs and demonstrated that PnWRKY9 positively affects plant defenses against the root rot pathogen. The data presented herein provide researchers with fundamental information regarding the regulatory mechanism mediating the coordinated activities of WRKY TFs and the JA signaling pathway in P. notoginseng responses to the root rot pathogen.

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.

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.

Chromosome-Scale Assembly of the Dendrobium nobile Genome Provides Insights Into the Molecular Mechanism of the Biosynthesis of the Medicinal Active Ingredient of Dendrobium

Orchids constitute approximately 10% of flowering plant species. However, only about 10 orchid genomes have been published. Metabolites are the main way through which orchids respond to their environment. Dendrobium nobile, belonging to Dendrobium, the second largest genus in Orchidaceae, has high ornamental, medicinal, and ecological value. D. nobile is the source of many popular horticultural varieties. Among the Dendrobium species, D. nobile has the highest amount of dendrobine, which is regarded as one of the criteria for evaluating medicinal quality. Due to lack of data and analysis at the genomic level, the biosynthesis pathways of dendrobine and other related medicinal ingredients in D. nobile are unknown. In this paper, we report a chromosome-scale reference genome of D. nobile to facilitate the investigation of its genomic characteristics for comparison with other Dendrobium species. The assembled genome size of D. nobile was 1.19 Gb. Of the sequences, 99.45% were anchored to 19 chromosomes. Furthermore, we identified differences in gene number and gene expression patterns compared with two other Dendrobium species by integrating whole-genome sequencing and transcriptomic analysis [e.g., genes in the polysaccharide biosynthesis pathway and upstream of the alkaloid (dendrobine) biosynthesis pathway]. Differences in the TPS and CYP450 gene families were also found among orchid species. All the above differences might contribute to the species-specific medicinal ingredient biosynthesis pathways. The metabolic pathway-related analysis will provide further insight into orchid responses to the environment. Additionally, the reference genome will provide important insights for further molecular elucidation of the medicinal active ingredients of Dendrobium and enhance the understanding of orchid evolution.

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.

Solid-State Fermentation With Aspergillus cristatus Enhances the Protopanaxadiol- and Protopanaxatriol-Associated Skin Anti-aging Activity of Panax notoginseng

A metabolomics approach was used to profile metabolites of Panax notoginseng fermented with Aspergillus cristatus in two ways, liquid-state fermentation (LF-P) and solid-state fermentation (SSF-P) and examine metabolite markers representing antioxidant activity and skin anti-aging. Protopanaxadiol (PPD) and protopanaxatriol (PPT) contents were higher in SSF-P than in LF-P and showed a multiplicative increase over the fermentation period of four days. PPD and PPT levels also correlated with antioxidant and anti-aging effects in skin, based on the mRNA expression of dermal extracellular matrix components. In the bioactivity validation assays, PPD and PPT significantly improved the expression of type-I collagen, fibrillin-1, and elastin in human dermal fibroblasts from both young and old subjects; these were comparable with the effects of the SSF-P extracts. Overall, our results suggest that changes in the metabolites of P. notoginseng fermented with A. cristatus enhance the quality and availability of bioactive compounds associated with skin anti-aging.

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.

Genome and transcriptome analysis to understand the role diversification of cytochrome P450 gene under excess nitrogen treatment

BACKGROUND

Panax notoginseng (Burk.) F. H. Chen (P. notoginseng) is a medicinal plant. Cytochrome P450 (CYP450) monooxygenase superfamily is involved in the synthesis of a variety of plant hormones. Studies have shown that CYP450 is involved in the synthesis of saponins, which are the main medicinal component of P. notoginseng. To date, the P. notoginseng CYP450 family has not been systematically studied, and its gene functions remain unclear.

RESULTS

In this study, a total of 188 PnCYP genes were identified, these genes were divided into 41 subfamilies and clustered into 9 clans. Moreover, we identified 40 paralogous pairs, of which only two had Ka/Ks ratio greater than 1, demonstrating that most PnCYPs underwent purification selection during evolution. In chromosome mapping and gene replication analysis, 8 tandem duplication and 11 segmental duplication events demonstrated that PnCYP genes were continuously replicating during their evolution. Gene ontology (GO) analysis annotated the functions of 188 PnCYPs into 21 functional subclasses, suggesting the functional diversity of these gene families. Functional divergence analyzed the members of the three primitive branches of CYP51, CYP74 and CYP97 at the amino acid level, and found some critical amino acid sites. The expression pattern of PnCYP450 related to nitrogen treatment was studied using transcriptome sequencing data, 10 genes were significantly up-regulated and 37 genes were significantly down-regulated. Combined with transcriptome sequencing analysis, five potential functional genes were screened. Quantitative real-time PCR (qRT-PCR) indicated that these five genes were responded to methyl jasmonate (MEJA) and abscisic acid (ABA) treatment.

CONCLUSIONS

These results provide a valuable basis for comprehending the classification and biological functions of PnCYPs, and offer clues to study their biological functions in response to nitrogen treatment.

Phytochemistry of ginsenosides: Recent advancements and emerging roles

Ginsenosides, a group of tetracyclic saponins, accounts for the nutraceutical and pharmaceutical relevance of the ginseng (Panax sp.) herb. Owing to the associated therapeutic potential of ginsenosides, their demand has been increased significantly in the last two decades. However, a slow growth cycle, low seed production, and long generation time of ginseng have created a gap between the demand and supply of ginsenosides. The biosynthesis of ginsenosides involves an intricate network of pathways with multiple oxidation and glycosylation reactions. However, the exact functions of some of the associated genes/proteins are still not completely deciphered. Moreover, ginsenoside estimation and extraction using analytical techniques are not feasible with high efficiency. The present review is a step forward in recapitulating the comprehensive aspects of ginsenosides including their distribution, structural diversity, biotransformation, and functional attributes in both plants and animals including humans. Moreover, ginsenoside biosynthesis in the potential plant sources and their metabolism in the human body along with major regulators and stimulators affecting ginsenoside biosynthesis have also been discussed. Furthermore, this review consolidates biotechnological interventions to enhance the biosynthesis of ginsenosides in their potential sources and advancements in the development of synthetic biosystems for efficient ginsenoside biosynthesis to meet their rising industrial demands.

TCM-Blast for traditional Chinese medicine genome alignment with integrated resources

The traditional Chinese medicine (TCM) genome project aims to reveal the genetic information and regulatory network of herbal medicines, and to clarify their molecular mechanisms in the prevention and treatment of human diseases. Moreover, the TCM genome could provide the basis for the discovery of the functional genes of active ingredients in TCM, and for the breeding and improvement of TCM. The traditional Chinese Medicine Basic Local Alignment Search Tool (TCM-Blast) is a web interface for TCM protein and DNA sequence similarity searches. It contains approximately 40G of genome data on TCMs, including protein and DNA sequence for 36 TCMs with high medical value.The development of a publicly accessible TCM genome alignment database hosted on the TCM-Blast website ( http://viroblast.pungentdb.org.cn/TCM-Blast/viroblast.php ) has expanded to query multiple sequence databases to obtain TCM genome data, and provide user-friendly output for easy analysis and browsing of BLAST results. The genome sequencing of TCMs helps to elucidate the biosynthetic pathways of important secondary metabolites and provides an essential resource for gene discovery studies and molecular breeding. The TCMs genome provides a valuable resource for the investigation of novel bioactive compounds and drugs from these TCMs under the guidance of TCM clinical practice. Our database could be expanded to other TCMs after the determination of their genome data.

Ginsenosides in Panax genus and their biosynthesis

Ginsenosides are a series of glycosylated triterpenoids which belong to protopanaxadiol (PPD)-, protopanaxatriol (PPT)-, ocotillol (OCT)- and oleanane (OA)-type saponins known as active compounds of Panax genus. They are accumulated in plant roots, stems, leaves, and flowers. The content and composition of ginsenosides are varied in different ginseng species, and in different parts of a certain plant. In this review, we summarized the representative saponins structures, their distributions and the contents in nearly 20 Panax species, and updated the biosynthetic pathways of ginsenosides focusing on enzymes responsible for structural diversified ginsenoside biosynthesis. We also emphasized the transcription factors in ginsenoside biosynthesis and non-coding RNAs in the growth of Panax genus plants, and highlighted the current three major biotechnological applications for ginsenosides production. This review covered advances in the past four decades, providing more clues for chemical discrimination and assessment on certain ginseng plants, new perspectives for rational evaluation and utilization of ginseng resource, and potential strategies for production of specific ginsenosides.

Planting Density Affects Panax notoginseng Growth and Ginsenoside Accumulation by Balancing Primary and Secondary Metabolism

Adjusting planting density is a common agricultural practice used to achieve maximum yields. However, whether the quality of medicinal herbs can be improved by implementing appropriate planting densities is still uncertain. The medicinal crop Panax notoginseng was used to analyze the effects of planting density on growth and ginsenoside accumulation, and the possible mechanisms of these effects were revealed through metabonomics. The results showed that P. notoginseng achieved high ginsenoside accumulation at high planting densities (8 × 8 and 10 × 10 cm), while simultaneously achieved high biomass and ginsenoside accumulation at moderate planting density of 15 × 15 cm. At the moderate planting density, the primary metabolism (starch and sucrose metabolism) and secondary metabolism (the biosynthesis of phytohormone IAA and ginsenoside) of the plants were significantly enhanced. However, the strong intraspecific competition at the high planting densities resulted in stress as well as the accumulation of phytohormones (SA and JA), antioxidants (gentiobiose, oxalic acid, dehydroascorbic acid) and other stress resistance-related metabolites. Interestingly, the dry biomass and ginsenoside content were significantly lower at low densities (20 × 20 and 30 × 30 cm) with low intraspecific competition, which disturbed normal carbohydrate metabolism by upregulating galactose metabolism. In summary, an appropriate planting density was benefit for the growth and accumulation of ginsenosides in P. notoginseng by balancing primary metabolism and secondary metabolism.

Genome survey of Zanthoxylum bungeanum and development of genomic-SSR markers in congeneric species

Zanthoxylum bungeanum, a spice and medicinal plant, is cultivated in many parts of China and some countries in Southeast Asia; however, data on its genome are lacking. In the present study, we performed a whole-genome survey and developed novel genomic-SSR markers of Z. bungeanum. Clean data (∼197.16 Gb) were obtained and assembled into 11185221 scaffolds with an N50 of 183 bp. K-mer analysis revealed that Z. bungeanum has an estimated genome size of 3971.92 Mb, and the GC content, heterozygous rate, and repeat sequence rate are 37.21%, 1.73%, and 86.04%, respectively. These results indicate that the genome of Z. bungeanum is complex. Furthermore, 27153 simple sequence repeat (SSR) loci were identified from 57288 scaffolds with a minimum length > 1 kb. Mononucleotide repeats (19706) were the most abundant type, followed by dinucleotide repeats (5154). The most common motifs were A/T, followed by AT/AT; these SSRs accounted for 71.42% and 11.84% of all repeats, respectively. A total of 21243 non-repeating primer pairs were designed, and 100 were randomly selected and validated by PCR analysis using DNA from 10 Z. bungeanum individuals and 5 Zanthoxylum armatum individuals. Finally, 36 polymorphic SSR markers were developed with polymorphism information content (PIC) values ranging from 0.16 to 0.75. Cluster analysis revealed that Z. bungeanum and Z. armatum could be divided into two major clusters, suggesting that these newly developed SSR markers are useful for genetic diversity and germplasm resource identification in Z. bungeanum and Z. armatum.

Dynamic evolution of Panax species

Background

Panax ginseng is one of the most valuable medicinal plants in Korea. However, deciphering its full genome sequence information for crop improvement has been hampered due to its complex genomic, genetic, and growth characteristics. Many efforts have been made in the past decade to overcome these limitations and understand the genome structure and the evolutionary history of P. ginseng.

Methods

This review aims to discuss the current status of genomic studies on P. ginseng and related species, and the experimental clues suggesting phylogenetic classification and evolutionary history of the genus Panax.

Conclusion

The development of sequencing technologies made genome sequencing of the large P. ginseng genome possible, providing fundamental information to deciphering the evolutionary history of P. ginseng and related species. P. ginseng went through two rounds of whole genome duplication events after diverging from the closest family Apiaceae, which was unveiled from complete whole genome sequences. Further in-depth comparative genome analysis with other related species and genera will uncover the evolutionary history as well as important morphological and ecological characteristics of Panax species.

The chromosome-level reference genome assembly for Panax notoginseng and insights into ginsenoside biosynthesis

Panax notoginseng, a perennial herb of the genus Panax in the family Araliaceae, has played an important role in clinical treatment in China for thousands of years because of its extensive pharmacological effects. Here, we report a high-quality reference genome of P. notoginseng, with a genome size up to 2.66 Gb and a contig N50 of 1.12 Mb, produced with third-generation PacBio sequencing technology. This is the first chromosome-level genome assembly for the genus Panax. Through genome evolution analysis, we explored phylogenetic and whole-genome duplication events and examined their impact on saponin biosynthesis. We performed a detailed transcriptional analysis of P. notoginseng and explored gene-level mechanisms that regulate the formation of characteristic tubercles. Next, we studied the biosynthesis and regulation of saponins at temporal and spatial levels. We combined multi-omics data to identify genes that encode key enzymes in the P. notoginseng terpenoid biosynthetic pathway. Finally, we identified five glycosyltransferase genes whose products catalyzed the formation of different ginsenosides in P. notoginseng. The genetic information obtained in this study provides a resource for further exploration of the growth characteristics, cultivation, breeding, and saponin biosynthesis of P. notoginseng.

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.

Advances in biosynthesis of triterpenoid saponins in medicinal plants

In recent years, biosynthesis of triterpenoid saponins in medicinal plants has been widely studied because of their active ingredients with diverse pharmacological activities. Various oxidosqualene cyclases, cytochrome P450 monooxygenases, uridine diphosphate glucuronosyltransferases, and transcription factors related to triterpenoid saponins biosynthesis have been explored and identified. In the biosynthesis of triterpenoid saponins, the progress of gene mining by omics-based sequencing, gene screening, gene function verification, catalyzing mechanism of key enzymes and gene regulation are summarized and discussed. By the progress of the biosynthesis pathway of triterpenoid saponins, the large-scale production of some triterpenoid saponins and aglycones has been achieved through plant tissue culture, transgenic plants and engineered yeast cells. However, the complex biosynthetic pathway and structural diversity limit the biosynthesis of triterpenoid saponins in different system. Special focus can further be placed on the systematic botany information of medicinal plants obtained from omics large dataset, and triterpenoid saponins produced by synthetic biology strategies, gene mutations and gene editing technology.

The Chromosome Level Genome and Genome-wide Association Study for the Agronomic Traits of Panax Notoginseng

The Chinese ginseng Panax notoginseng is a domesticated herb with significant medicinal and economic value. Here we report a chromosome-level P. notoginseng genome assembly with a high (∼79%) repetitive sequence content. The juxtaposition with the widely distributed, closely related Korean ginseng (Panax ginseng) genome revealed contraction of plant defense genes (in particular R-genes) in the P. notoginseng genome. We also investigated the reasons for the larger genome size of Panax species, revealing contributions from two Panax-specific whole-genome duplication events and transposable element expansion. Transcriptome data and comparative genome analysis revealed the candidate genes involved in the ginsenoside synthesis pathway. We also performed a genome-wide association study on 240 cultivated P. notoginseng individuals and identified the associated genes with dry root weight (63 genes) and stem thickness (168 genes). The P. notoginseng genome represents a critical step toward harnessing the full potential of an economically important and enigmatic plant.

The unprecedented diversity of UGT94-family UDP-glycosyltransferases in Panax plants and their contribution to ginsenoside biosynthesis

More than 150 ginsenosides have been isolated and identified from Panax plants. Ginsenosides with different glycosylation degrees have demonstrated different chemical properties and bioactivity. In this study, we systematically cloned and characterized 46 UGT94 family UDP-glycosyltransferases (UGT94s) from a mixed Panax ginseng/callus cDNA sample with high amino acid identity. These UGT94s were found to catalyze sugar chain elongation at C3-O-Glc and/or C20-O-Glc of protopanaxadiol (PPD)-type, C20-O-Glc or C6-O-Glc of protopanaxatriol (PPT)-type or both C3-O-Glc of PPD-type and C6-O-Glc of PPT-type or C20-O-Glc of PPD-type and PPT-type ginsenosides with different efficiencies. We also cloned 26 and 51 UGT94s from individual P. ginseng and P. notoginseng plants, respectively; our characterization results suggest that there is a group of UGT94s with high amino acid identity but diverse functions or catalyzing activities even within individual plants. These UGT94s were classified into three clades of the phylogenetic tree and consistent with their catalytic function. Based on these UGT94s, we elucidated the biosynthetic pathway of a group of ginsenosides. Our present results reveal a series of UGTs involved in second sugar chain elongation of saponins in Panax plants, and provide a scientific basis for understanding the diverse evolution mechanisms of UGT94s among plants.

Chromosome Level Genome Assembly of Andrographis paniculata

Andrographis paniculata (Chinese name: Chuanxinlian) is an annual dicotyledonous medicinal plant widely grown in China and Southeast Asia. The dried plant has a highly acclaimed usage in the traditional Chinese medicine for its antipyretic, anti-inflammatory, and analgesic effects. In order to help delineate the biosynthetic pathways of various secondary metabolites, we report in this study a high-quality reference genome for A. paniculata. With the help of both PacBio single molecule real time sequencing and Illumina sequencing reads for error correction, the A. paniculata genome was assembled into a total size of 284 Mb with a contig N50 size of 5.14 Mb. The contigs were further assembled into 24 pseudo-chromosomes by the Hi-C technique. We also analyzed the gene families (e.g., KSL, and CYP450) whose protein products are essential for synthesizing bioactive compounds in A. paniculata. In conclusion, the high-quality A. paniculata genome assembly builds the foundation for decoding the biosynthetic pathways of various medicinal compounds.

Whole-genome, transcriptome, and methylome analyses provide insights into the evolution of platycoside biosynthesis in Platycodon grandiflorus, a medicinal plant

Triterpenoid saponins (TSs) are common plant defense phytochemicals with potential pharmaceutical properties. Platycodon grandiflorus (Campanulaceae) has been traditionally used to treat bronchitis and asthma in East Asia. The oleanane-type TSs, platycosides, are a major component of the P. grandiflorus root extract. Recent studies show that platycosides exhibit anti-inflammatory, antiobesity, anticancer, antiviral, and antiallergy properties. However, the evolutionary history of platycoside biosynthesis genes remains unknown. In this study, we sequenced the genome of P. grandiflorus and investigated the genes involved in platycoside biosynthesis. The draft genome of P. grandiflorus is 680.1 Mb long and contains 40,017 protein-coding genes. Genomic analysis revealed that the CYP716 family genes play a major role in platycoside oxidation. The CYP716 gene family of P. grandiflorus was much larger than that of other Asterid species. Orthologous gene annotation also revealed the expansion of β-amyrin synthases (bASs) in P. grandiflorus, which was confirmed by tissue-specific gene expression. In these expanded gene families, we identified key genes showing preferential expression in roots and association with platycoside biosynthesis. In addition, whole-genome bisulfite sequencing showed that CYP716 and bAS genes are hypomethylated in P. grandiflorus, suggesting that epigenetic modification of these two gene families affects platycoside biosynthesis. Thus whole-genome, transcriptome, and methylome data of P. grandiflorus provide novel insights into the regulation of platycoside biosynthesis by CYP716 and bAS gene families.

Tandem gene duplications drive divergent evolution of caffeine and crocin biosynthetic pathways in plants

Background

Plants have evolved a panoply of specialized metabolites that increase their environmental fitness. Two examples are caffeine, a purine psychotropic alkaloid, and crocins, a group of glycosylated apocarotenoid pigments. Both classes of compounds are found in a handful of distantly related plant genera (Coffea, Camellia, Paullinia, and Ilex for caffeine; Crocus, Buddleja, and Gardenia for crocins) wherein they presumably evolved through convergent evolution. The closely related Coffea and Gardenia genera belong to the Rubiaceae family and synthesize, respectively, caffeine and crocins in their fruits.

Results

Here, we report a chromosomal-level genome assembly of Gardenia jasminoides, a crocin-producing species, obtained using Oxford Nanopore sequencing and Hi-C technology. Through genomic and functional assays, we completely deciphered for the first time in any plant the dedicated pathway of crocin biosynthesis. Through comparative analyses with Coffea canephora and other eudicot genomes, we show that Coffea caffeine synthases and the first dedicated gene in the Gardenia crocin pathway, GjCCD4a, evolved through recent tandem gene duplications in the two different genera, respectively. In contrast, genes encoding later steps of the Gardenia crocin pathway, ALDH and UGT, evolved through more ancient gene duplications and were presumably recruited into the crocin biosynthetic pathway only after the evolution of the GjCCD4a gene.

Conclusions

This study shows duplication-based divergent evolution within the coffee family (Rubiaceae) of two characteristic secondary metabolic pathways, caffeine and crocin biosynthesis, from a common ancestor that possessed neither complete pathway. These findings provide significant insights on the role of tandem duplications in the evolution of plant specialized metabolism.

Application of Identification and Evaluation Techniques for Ethnobotanical Medicinal Plant of Genus Panax: A Review

Genus Panax, as worldwide medicinal plants, has a medical history for thousands of years. Most of the entire genus are traditional ethnobotanical medicine in China, Myanmar, Thailand, Vietnam and Laos, which have given rise to international attention and use. This paper reviewed more than 210 articles and related books on the research of Panax medicinal plants and their Chinese patent medicines published in the last 30 years. The purpose was to review and summarize the species classification, geographical distribution, and ethnic minorities medicinal records of the genus Panax, and further to review the analytical tools and data analysis methods for the authentication and quality assessment of Panax medicinal materials and Chinese patent medicines. Five main technologies applied in the identification and evaluation of Panax have been introduced and summarized. Chromatography was the most widely used one. Further research and development of molecular identification technology had the potential to become a mainstream identification technology. In addition, some novel, controversial, and worthy methods including electronic noses, electronic eyes, and DNA barcoding were also introduced. At the same time, more than 80% of the researches were carried out by a combination of chemometric pattern-recognition technologies and multi-analysis technologies. All the technologies and methods applied can provide strong support and guarantee for the identification and evaluation of genus Panax, and also conduce to excellent reference value for the development and in-depth research of new technologies in Panax.

Till 2018: a survey of biomolecular sequences in genus Panax

Ginseng is popularly known to be the king of ancient medicines and is used widely in most of the traditional medicinal compositions due to its various pharmaceutical properties. Numerous studies are being focused on this plant's curative effects to discover their potential health benefits in most human diseases, including cancer- the most life-threatening disease worldwide. Modern pharmacological research has focused mainly on ginsenosides, the major bioactive compounds of ginseng, because of their multiple therapeutic applications. Various issues on ginseng plant development, physiological processes, and agricultural issues have also been studied widely through state-of-the-art, high-throughput sequencing technologies. Since the beginning of the 21st century, the number of publications on ginseng has rapidly increased, with a recent count of more than 6,000 articles and reviews focusing notably on ginseng. Owing to the implementation of various technologies and continuous efforts, the ginseng plant genomes have been decoded effectively in recent years. Therefore, this review focuses mainly on the cellular biomolecular sequences in ginseng plants from the perspective of the central molecular dogma, with an emphasis on genomes, transcriptomes, and proteomes, together with a few other related studies.

The transcriptome variations of Panaxnotoginseng roots treated with different forms of nitrogen fertilizers

BACKGROUND

The sensitivity of plants to ammonia is a worldwide problem that limits crop production. Excessive use of ammonium as the sole nitrogen source results in morphological and physiological disorders, and retarded plant growth.

RESULTS

In this study we found that the root growth of Panax notoginseng was inhibited when only adding ammonium nitrogen fertilizer, but the supplement of nitrate fertilizer recovered the integrity, activity and growth of root. Twelve RNA-seq profiles in four sample groups were produced and analyzed to identify deregulated genes in samples with different treatments. In comparisons to NH[Formula: see text] treated samples, ACLA-3 gene is up-regulated in samples treated with NO[Formula: see text] and with both NH[Formula: see text] and NO[Formula: see text], which is further validated by qRT-PCR in another set of samples. Subsequently, we show that the some key metabolites in the TCA cycle are also significantly enhanced when introducing NO[Formula: see text]. These potentially enhance the integrity and recover the growth of Panax notoginseng roots.

CONCLUSION

These results suggest that the activated TCA cycle, as demonstrated by up-regulation of ACLA-3 and several key metabolites in this cycle, contributes to the increased Panax notoginseng root yield when applying both ammonium and nitrate fertilizer.

Genome-Wide Identification of WRKY Transcription Factors in the Asteranae

The WRKY transcription factors family, which participates in many physiological processes in plants, constitutes one of the largest transcription factor families. The Asterales and the Apiales are two orders of flowering plants in the superorder Asteranae. Among the members of the Asterales, globe artichoke (Cynara cardunculus var. scolymus L.), sunflower (Helianthus annuus L.), and lettuce (Lactuca sativa L.) are important economic crops worldwide. Within the Apiales, ginseng (Panax ginseng C. A. Meyer) and Panax notoginseng (Burk.) F.H. Chen are important medicinal plants, while carrot (Daucus carota subsp. carota L.) has significant economic value. Research involving genome-wide identification of WRKY transcription factors in the Asterales and the Apiales has been limited. In this study, 490 WRKY genes, 244 from three species of the Apiales and 246 from three species of the Asterales, were identified and categorized into three groups. Within each group, WRKY motif characteristics and gene structures were similar. WRKY gene promoter sequences contained light responsive elements, core regulatory elements, and 12 abiotic stress cis-acting elements. WRKY genes were evenly distributed on each chromosome. Evidence of segmental and tandem duplication events was found in all six species in the Asterales and the Apiales, with segmental duplication inferred to play a major role in WRKY gene evolution. Among the six species, we uncovered 54 syntenic gene pairs between globe artichoke and lettuce. The six species are thus relatively closely related, consistent with their traditional taxonomic placement in the Asterales. This study, based on traditional species classifications, was the first to identify WRKY transcription factors in six species from the Asteranae. Our results lay a foundation for further understanding of the role of WRKY transcription factors in species evolution and functional differentiation.

Transcriptomic profiling reveals MEP pathway contributing to ginsenoside biosynthesis in Panax ginseng

Background

Panax ginseng C. A. Mey is one of famous medicinal herb plant species. Its major bioactive compounds are various ginsenosides in roots and rhizomes. It is commonly accepted that ginsenosides are synthesized from terpene precursors, IPP and DMAPP, through the cytoplasmic mevalonate (MVA) pathway. Another plastic 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway was proved also contributing to ginsenoside generation in the roots of P. ginseng by using specific chemical inhibitors recently. But their gene expression characteristics are still under reveal in P. ginseng. With the development of the high-throughput next generation sequencing (NGS) technologies, we have opportunities to discover more about the complex ginsenoside biosynthesis pathways in P. ginseng.

Results

We carried out deep RNA sequencing and comprehensive analyses on the ginseng root samples of 1–5 years old and five different tissues of 5 years old ginseng plants. The de novo assembly totally generated 48,165 unigenes, including 380 genes related to ginsenoside biosynthesis and all the genes encoding the enzymes of the MEP pathway and the MVA pathway. We further illustrated the gene expression profiles related to ginsenoside biosynthesis among 1–5 year-old roots and different tissues of 5 year-old ginseng plants. Particularly for the first time, we revealed that the gene transcript abundances of the MEP pathway were similar to those of the MVA pathway in ginseng roots but higher in ginseng leaves. The IspD was predicated to be the rate-limiting enzyme in the MEP pathway through both co-expression network and gene expression profile analyses.

Conclusions

At the transcriptional level, the MEP pathway has similar contribution to ginsenoside biosynthesis in ginseng roots, but much higher in ginseng leaves, compared with the MVA pathway. The IspD might be the key enzyme for ginsenoside generation through the MEP pathway. These results provide new information for further synthetic biology study on ginsenoside metabolic regulation.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5718-x) contains supplementary material, which is available to authorized users.

Transcriptome analysis of Panax zingiberensis identifies genes encoding oleanolic acid glucuronosyltransferase involved in the biosynthesis of oleanane-type ginsenosides

Oleanolic acid glucuronosyltransferase (OAGT) genes synthesizing the direct precursor of oleanane-type ginsenosides were discovered. The four recombinant proteins of OAGT were able to transfer glucuronic acid at C-3 of oleanolic acid that yields oleanolic acid 3-O-β-glucuronide. Ginsenosides are the primary active components in the genus Panax, and great efforts have been made to elucidate the mechanisms underlying dammarane-type ginsenoside biosynthesis. However, there is limited information on oleanane-type ginsenosides. Here, high-performance liquid chromatography analysis demonstrated that oleanane-type ginsenosides (particularly ginsenoside Ro and chikusetsusaponin IV and IVa) are the abundant ginsenosides in Panax zingiberensis, an extremely endangered Panax species in southwest China. These ginsenosides are derived from oleanolic acid 3-O-β-glucuronide, which may be formed from oleanolic acid catalyzed by an unknown oleanolic acid glucuronosyltransferase (OAGT). Transcriptomic analysis of leaves, stems, main roots, and fibrous roots of P. zingiberensis was performed, and a total of 46,098 unigenes were obtained, including all the identified homologous genes involved in ginsenoside biosynthesis. The most upstream genes were highly expressed in the leaves, and the UDP-glucosyltransferase genes were highly expressed in the roots. This finding indicated that the precursors of ginsenosides are mainly synthesized in the leaves and transported to different parts for the formation of particular ginsenosides. For the first time, enzyme activity assay characterized four genes (three from P. zingiberensis and one from P. japonicus var. major, another Panax species with oleanane-type ginsenosides) encoding OAGT, which particularly transfer glucuronic acid at C-3 of oleanolic acid to form oleanolic acid 3-O-β-glucuronide. Taken together, our study provides valuable genetic information for P. zingiberensis and the genes responsible for synthesizing the direct precursor of oleanane-type ginsenosides.

Genome survey sequencing for the characterization of genetic background of Dracaena cambodiana and its defense response during dragon's blood formation

Dragon's blood collected from the genus Dracaena is used as a renowned traditional medicine in various cultures worldwide. However, the genetics of the genus Dracaena and the formation mechanism of dragon's blood remain poorly understood. Here, we generate the first draft genome reference assembly of an elite Chinese Dracaena species, Dracaena cambodiana, from next-generation sequencing data with 89.46× coverage. The reads were assembled into 2,640,704 contigs with an N50 length of 1.87 kb, and a 1.05 Gb assembly was finally assembled with 2,379,659 scaffolds. Furthermore, 97.75% of the 267,243 simple sequence repeats identified from these scaffolds were mononucleotide, dinucleotide, and trinucleotide repeats. Among all 53,700 predicted genes, 158 genes involved in cell wall and plant hormone synthesis and reactive oxygen species scavenging showed altered regulation during the formation of dragon's blood. This study provides a genomic characterization of D. cambodiana and improves understanding of the molecular mechanism of dragon's blood formation. This report represents the first genome-wide characterization of a Dracaena species in the Asparagaceae.

Identification of potential genes involved in triterpenoid saponins biosynthesis in Gleditsia sinensis by transcriptome and metabolome analyses

Gleditsia sinensis is widely used as a medicinal plant in Asia, especially in China. Triterpenes, alkaloids, and sterols were isolated from Gleditsia species. Among them, triterpenoid saponins are very important metabolites owing to their various pharmacological activities. However, the triterpenoid saponin biosynthesis pathway has not been well characterized. In the present study, we performed de novo transcriptome assembly for 14.3 Gbps of clean reads sequenced from nine tissues of G. sinensis. The results showed that 81,511 unique transcripts (unitranscripts) (47,855 unigenes) were constructed, of which 31,717 unigenes were annotated with Gene Ontology and EC numbers by Blast2GO against the NCBI-nr protein database. We also analyzed the metabolite contents in the same nine tissues by LS–MS/MS, and saponins including gleditsioside I were found in fruit at higher levels. Many of the genes with tissue-specific expression in fruit are involved in the flavonoid biosynthesis pathway, and many of those have UDP-glucosyltransferase (UGT) activity. We constructed a saponin biosynthesis pathway and identified two key enzyme families in the triterpenoid saponin biosynthesis pathway, cytochrome P450 and UDP-glucosyltransferase, that are encoded by 37 unigenes and 77 unigenes, respectively. CYP72A, CYP716A, and CYP88D, which are known as key enzymes for saponin biosynthesis, were also identified among the P450s. Our results provide insight into the secondary metabolite biosynthesis and serve as important resources for future research and cultivation of G. sinensis.

Trends in herbgenomics

From Shen Nong's Herbal Classic (Shennong Bencao Jing) to the Compendium of Materia Medica (Bencao Gangmu) and the first scientific Nobel Prize for the mainland of China, each milestone in the historical process of the development of traditional Chinese medicine (TCM) involves screening, testing and integrating. After thousands of years of inheritance and development, herbgenomics (bencaogenomics) has bridged the gap between TCM and international advanced omics studies, promoting the application of frontier technologies in TCM. It is a discipline that uncovers the genetic information and regulatory networks of herbs to clarify their molecular mechanism in the prevention and treatment of human diseases. The main theoretical system includes genomics, functional genomics, proteomics, transcriptomics, metabolomics, epigenomics, metagenomics, synthetic biology, pharmacogenomics of TCM, and bioinformatics, among other fields. Herbgenomics is mainly applicable to the study of medicinal model plants, genomic-assisted breeding, herbal synthetic biology, protection and utilization of gene resources, TCM quality evaluation and control, and TCM drug development. Such studies will accelerate the application of cutting-edge technologies, revitalize herbal research, and strongly promote the development and modernization of TCM.