The Medicinal Herb Panax notoginseng Genome Provides Insights into Ginsenoside Biosynthesis and Genome Evolution

Spatiotemporal expression analysis of jasmonic acid and saponin-related genes uncovers a potential biosynthetic regulation in Panax notoginseng

BACKGROUND

Sanqi, the root of Panax notoginseng, has long been recognized for its therapeutic effects on cardiovascular diseases. Saponins, including ginsenosides and notoginsenosides, are the main bioactive components of P. notoginseng. The biosynthesis of saponins is closely related to the defense responses orchestrated by endogenous hormones.

RESULTS

To provide new insights into the underlying role of phytohormone jasmonic acid (JA) in the synthesis and regulation of saponins, we performed an ultra-performance liquid chromatography analysis of different tissues of P. notoginseng aged 2-4 years. Moreover, by combined evaluation of saponin content and transcriptome profiling of each tissue, the spatial and temporal distribution of saponins was analyzed. N notoginsenoside R1, ginsenoside Rb1 and ginsenoside Rd accumulated in the underground tissues, including the root, tuqi, fibril and rhizome. In agreement with this data, the corresponding genes of the endogenous hormone JAs, especially coronatine insensitive 1 (COI1) and myelocytomatosis proteins 2 (MYC2), were predominantly expressed in the underground tissues. The tissue- and age-specific distribution of saponins was consistent with the expression of genes involved in JA biosynthetic, metabolic and signaling pathways.

CONCLUSION

The present study has revealed the temporal and spatial effects of endogenous phtohormones in the synthesis and regulation of notoginsenosides, which will provide a significant impact on improving the ecological planting technology, cultivating new high-quality varieties and protecting the rare resources of medicinal P. notoginseng. © 2024 Society of Chemical Industry.

Molecular regulation of the key specialized metabolism pathways in medicinal plants

The basis of modern pharmacology is the human ability to exploit the production of specialized metabolites from medical plants, for example, terpenoids, alkaloids, and phenolic acids. However, in most cases, the availability of these valuable compounds is limited by cellular or organelle barriers or spatio-temporal accumulation patterns within different plant tissues. Transcription factors (TFs) regulate biosynthesis of these specialized metabolites by tightly controlling the expression of biosynthetic genes. Cutting-edge technologies and/or combining multiple strategies and approaches have been applied to elucidate the role of TFs. In this review, we focus on recent progress in the transcription regulation mechanism of representative high-value products and describe the transcriptional regulatory network, and future perspectives are discussed, which will help develop high-yield plant resources.

Cucumber mosaic virus impairs the physiological homeostasis of Panax notoginseng and induces saponin-mediated resistance

As an important medicinal plant, Panax notoginseng often suffers from various abiotic and biotic stresses during its growth, such as drought, heavy metals, fungi, bacteria and viruses. In this study, the symptom and physiological parameters of cucumber mosaic virus (CMV)-infected P. notoginseng were analyzed and the RNA-seq was performed. The results showed that CMV infection affected the photosynthesis of P. notoginseng, caused serious oxidative damage to P. notoginseng and increased the activity of several antioxidant enzymes. Results of transcriptome analysis and corresponding verification showed that CMV infection changed the expression of genes related to plant defense and promoted the synthesis of P. notoginseng saponins to a certain extent, which may be defensive ways of P. notoginseng against CMV infection. Furthermore, pretreatment plants with saponins reduced the accumulation of CMV. Thus, our results provide new insights into the role of saponins in P. notoginseng response to virus infection.

First Contiguous Genome Assembly of Japanese Lady Bell (Adenophora triphylla) and Insights into Development of Different Leaf Types

Adenophora triphylla is an important medicinal and food plant found in East Asia. This plant is rich in secondary metabolites such as triterpenoid saponin, and its leaves can develop into different types, such as round and linear, depending on the origin of germination even within the same species. Despite this, few studies have comprehensively characterized the development processes of different leaf types and triterpenoid saponin pathways in this plant. Herein, we provide the first report of a high-quality genome assembly of A. triphylla based on a combination of Oxford Nanopore Technologies and Illumina sequencing methods. Its genome size was estimated to be 2.6 Gb, and the assembled genome finalized as 2.48 Gb, containing 57,729 protein-coding genes. Genome completeness was assessed as 95.6% using the Benchmarking Universal Single-Copy Orthologs score. The evolutionary divergence of A. triphylla was investigated using the genomes of five plant species, including two other species in the Campanulaceae family. The species A. triphylla diverged approximately 51-118 million years ago from the other four plants, and 579 expanded/contracted gene families were clustered in the Gene Ontology terms. The expansion of the β-amyrin synthase (bAS) gene, a key enzyme in the triterpenoid saponin pathway, was identified in the A. triphylla genome. Furthermore, transcriptome analysis of the two leaf types revealed differences in the activity of starch, sucrose, unsaturated fatty acid pathways, and oxidoreductase enzymes. The heat and endoplasmic reticulum pathways related to plant stress were active in the development of round type leaf, while an enhancement of pyrimidine metabolism related to cell development was confirmed in the development of the linear type leaf. This study provides insight into the evolution of bAS genes and the development of different leaf types in A. triphylla.

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.

Elucidation of the melitidin biosynthesis pathway in pummelo

Specialized plant metabolism is a rich resource of compounds for drug discovery. The acylated flavonoid glycoside melitidin is being developed as an anti-cholesterol statin drug candidate, but its biosynthetic route in plants has not yet been fully characterized. Here, we describe the gene discovery and functional characterization of a new flavonoid gene cluster (UDP-glucuronosyltransferases (CgUGTs), 1,2 rhamnosyltransferase (Cg1,2RhaT), acyltransferases (CgATs)) that is responsible for melitidin biosynthesis in pummelo (Citrus grandis (L.) Osbeck). Population variation analysis indicated that the tailoring of acyltransferases, specific for bitter substrates, mainly determine the natural abundance of melitidin. Moreover, 3-hydroxy-3-methylglutaryl-CoA reductase enzyme inhibition assays showed that the product from this metabolic gene cluster, melitidin, may be an effective anti-cholesterol statin drug candidate. Co-expression of these clustered genes in Nicotiana benthamiana resulted in the formation of melitidin, demonstrating the potential for metabolic engineering of melitidin in a heterologous plant system. This study establishes a biosynthetic pathway for melitidin, which provides genetic resources for the breeding and genetic improvement of pummelo aimed at fortifying the content of biologically active metabolites.

Genome-wide analysis of Panax MADS-box genes reveals role of PgMADS41 and PgMADS44 in modulation of root development and ginsenoside synthesis

Panax root is an important material used in food and medicine. Its cultivation and production usually depend on root shape and ginsenoside content. There is limited understanding about the synergistic regulatory mechanisms underlying root development and ginsenoside accumulation in Panax. MADS-box transcription factors possibly play a significant role in regulation of root growth and secondary metabolites. In this study, we identified MADS-box transcription factors of Panax, and found high expression levels of SVP, ANR1 and SOC1-like clade genes in its roots. We confirmed that two SOC1-like genes, PgMADS41 and PgMADS44, bind to expansion gene promoters (PgEXLB5 and PgEXPA13), which contribute to root growth, and to SE-4, CYP716A52v2-4, and β-AS-13 promoters, which participate in ginsenoside Ro biosynthesis. These two genes were found to increase lateral root number and main root length in transgenic Arabidopsis thaliana by improving AtEXLA1, AtEXLA3, AtEXPA5, and AtEXPA6 gene expression. As a non-phytohormone regulatory tool, Ro can stimulate adventitious root growth by influencing their expression and ginsenoside accumulation. Our study provides new insights into the coordinated regulatory function of SOC1-like clade genes in Panax root development and triterpenoid accumulation, paving the way towards understanding root formation and genetic improvement in Panax.

Phylogenomics of Aralia sect. Aralia (Araliaceae): Signals of hybridization and insights into its species delimitations and intercontinental biogeography

Genome-scale data have significantly increased the number of informative characters for phylogenetic analyses and recent studies have also revealed widespread phylogenomic discordance in many plant lineages. Aralia sect. Aralia is a small plant lineage (14 spp.) of the ginseng family Araliaceae with a disjunct distribution between eastern Asia (11 spp.) and North America (3 spp.). We herein employ sequences of hundreds of nuclear loci and the complete plastomes using targeted sequence capture and genome skimming to reconstruct the phylogenetic and biogeographic history of this section. We detected substantial conflicts among nuclear genes, yet different analytical strategies generated largely congruent topologies from the nuclear data. Significant cytonuclear discordance was detected, especially concerning the positions of the three North American species. The phylogenomic results support two intercontinental disjunctions: (1) Aralia californica of western North America is sister to the eastern Asian clade consisting of A. cordata and A. continentalis in the nuclear tree, and (2) the eastern North American A. racemosa forms a clade with A. bicrenata from southwestern North America, and the North American A. racemosa - A. bicrenata clade is then sister to the eastern Asian clade consisting of A. glabra (Japan), A. fargesii (C China), and A. apioides and A. atropurpurea (the Hengduan Mountains). Aralia cordata is supported to be disjunctly distributed in Japan, Taiwan, the Ulleung island of Korea, and in Central, Southwest and South China, and Aralia continentalis is redefined with a narrower distribution in Northeast China, eastern Russia and peninsular Korea.

Chromosome-Scale Genome Assembly and Triterpenoid Saponin Biosynthesis in Korean Bellflower (Platycodon grandiflorum)

Platycodon grandiflorum belongs to the Campanulaceae family and is an important medicinal and food plant in East Asia. However, on the whole, the genome evolution of P. grandiflorum and the molecular basis of its major biochemical pathways are poorly understood. We reported a chromosome-scale genome assembly of P. grandiflorum based on a hybrid method using Oxford Nanopore Technologies, Illumina sequences, and high-throughput chromosome conformation capture (Hi-C) analysis. The assembled genome was finalized as 574 Mb, containing 41,355 protein-coding genes, and the genome completeness was assessed as 97.6% using a Benchmarking Universal Single-Copy Orthologs analysis. The P. grandiflorum genome comprises nine pseudo-chromosomes with 56.9% repeat sequences, and the transcriptome analysis revealed an expansion of the 14 beta-amylin genes related to triterpenoid saponin biosynthesis. Our findings provide an understanding of P. grandiflorum genome evolution and enable genomic-assisted breeding for the mass production of important components such as triterpenoid saponins.

Effects of different soil moisture on the growth, quality, and root rot disease of organic Panax notoginseng cultivated under pine forests

The under-forest economy in the agroforestry system can improve land use efficiency, protect ecological environment, and promote arable land sustainable development. However, the effects of soil moisture in the forest and irrigation strategies on the healthy growth of intercropping crops are still incomplete. Here, considering the organic Panax notoginseng cultivated under pine forests (PPF) as the research object, we explored the effects of different soil moisture on the physiological state, yield, quality and disease occurrence of PPF. Our results suggested that 80-85% and 95-100% field capacity (FC) treatments were more conducive to increased photosynthetic rate and biomass accumulation of PPF, but 50-55% and 65-70% FC treatments were more conducive to the accumulation of saponins in PPF leaves. Notably, the root rot index of PPF was highest under 95-100% FC (19.51) treatment, significantly higher than that under 65-70% FC (8.44) and 80-85% FC (10.21) treatments. Further, the rhizosphere microorganisms of PPF under different soil moisture treatments were sequenced, and the sequencing data analysis revealed that high soil moisture (95-100% FC) could destroy the microbial diversity balance and cause the accumulation of pathogens (Fusarium oxysporum and Ilyonectria radicicola), leading to a high incidence of root rot. The incidence of PPF root rot was negatively correlated with rhizosphere microbial diversity. Overall, our results highlight that the quantitative irrigation (80-85% FC) is conducive to maintaining the balance between yield, saponin content and disease occurrence of PPF, providing a practical basis for PPF irrigation strategy and promoting the sustainable development of PPF agroforestry system.

New Glycosyltransferases in Panax notoginseng Perfect Main Ginsenosides Biosynthetic Pathways

Ginsenosides, the main bioactive ingredients of the Panax genus, are dammarane or oleanane triterpenoids with glycosylated modifications at C3/C6/C20 hydroxyls or C28 carboxyl, and their diverse glycosylation pattern has attracted great attention. However, the biosynthesis of some important saponins is still unclear. In this study, six UGTs were characterized, two of which were novel. PnUGT71A3 catalyzes not only the C6 hydroxyl glycosylation of protopanaxatriol (PPT) and F1 to form Rh1 and Rg1, respectively, but also the C20 hydroxyl glycosylation of protopanaxadiol (PPD)-type Rg3 to generate Rd. Especially, PnUGT94M1 is UDP-β-l-rhamnose (UDP-Rha)-dependent, regioselectively catalyzing the C2' hydroxyl rhamnosylation of C6 glucose of the PPT-type ginsenosides Rg1 and Rh1 to generate ginsenosides Re and Rg2, respectively. Site-directed mutagenesis showed that His21, Asp120, Ser363, and Pro372 are key residues, and the triple mutant (G344S/G345S/L346T) highly improved the activity toward Rg1 and Rh1. The findings in this study, perfect main ginsenosides biosynthetic pathways in the Panax genus, expand the biocatalyst toolbox for ginsenoside production and show that the PSPG motif is one of the options to modify UGTs to improve their activities.

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.

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.

Using interdisciplinary, phylogeny-guided approaches to understand the evolution of plant metabolism

To cope with relentless environmental pressures, plants produce an arsenal of structurally diverse chemicals, often called specialized metabolites. These lineage-specific compounds are derived from the simple building blocks made by ubiquitous core metabolic pathways. Although the structures of many specialized metabolites are known, the underlying metabolic pathways and the evolutionary events that have shaped the plant chemical diversity landscape are only beginning to be understood. However, with the advent of multi-omics data sets and the relative ease of studying pathways in previously intractable non-model species, plant specialized metabolic pathways are now being systematically identified. These large datasets also provide a foundation for comparative, phylogeny-guided studies of plant metabolism. Comparisons of metabolic traits and features like chemical abundances, enzyme activities, or gene sequences from phylogenetically diverse plants provide insights into how metabolic pathways evolved. This review highlights the power of studying evolution through the lens of comparative biochemistry, particularly how placing metabolism into a phylogenetic context can help a researcher identify the metabolic innovations enabling the evolution of structurally diverse plant metabolites.

Chromosome-level genome of Entada phaseoloides provides insights into genome evolution and triterpenoid saponins biosynthesis

As a medicinal herbal plant, Entada phaseoloides has high levels of secondary metabolites, particularly triterpenoid saponins, which are important resources for scientific research and medical applications. However, the lack of a reference genome for this genus has limited research on its evolution and utilization of its medicinal potential. In this study, we report a chromosome-scale genome assembly for E. phaseoloides using Illumina, Nanopore long reads, and high-throughput chromosome conformation capture technology. The assembled reference genome is 456.18 Mb (scaffold N50=30.9 Mb; contig N50=6.34 Mb) with 95.71 % of the sequences anchored onto 14 pseudochromosomes. E. phaseoloides was estimated to diverge from the Leguminosae lineage at approximately 72.0 million years ago. With the integration of transcriptomic and metabolomic data, gene expression patterns and metabolite profiling of E. phaseoloides were determined in different tissues. The pattern of gene expression and metabolic profile of the kernel were distinct from those of other tissues. Furthermore, the evolution of certain gene families involved in the biosynthesis of triterpenoid saponins and terpenes was analyzed and offer new insights into the formation of these two metabolites. Four CYP genes, one UGT gene and related transcription factors were identified as candidate genes contributing to regulation of triterpenoid saponins biosynthesis. As the first high-quality assembled reference genome in the genus Entada, it will not only provide new information for the evolutionary study of this genus and conservation biology of E. phaseoloides but also lay a foundation for the formation and utilization of secondary metabolites in medicinal plants.

Reshuffling of the ancestral core-eudicot genome shaped chromatin topology and epigenetic modification in Panax

All extant core-eudicot plants share a common ancestral genome that has experienced cyclic polyploidizations and (re)diploidizations. Reshuffling of the ancestral core-eudicot genome generates abundant genomic diversity, but the role of this diversity in shaping the hierarchical genome architecture, such as chromatin topology and gene expression, remains poorly understood. Here, we assemble chromosome-level genomes of one diploid and three tetraploid Panax species and conduct in-depth comparative genomic and epigenomic analyses. We show that chromosomal interactions within each duplicated ancestral chromosome largely maintain in extant Panax species, albeit experiencing ca. 100–150 million years of evolution from a shared ancestor. Biased genetic fractionation and epigenetic regulation divergence during polyploidization/(re)diploidization processes generate remarkable biochemical diversity of secondary metabolites in the Panax genus. Our study provides a paleo-polyploidization perspective of how reshuffling of the ancestral core-eudicot genome leads to a highly dynamic genome and to the metabolic diversification of extant eudicot plants.

The role of polyploidization generated genomic diversity in shaping the hierarchical genome architecture remains unclear. Here, the authors show that repatterning of the ancestral eudicot genome has resulted in multi-dimensional genome plasticity and secondary metabolite diversification via comparisons of Panax genomes.

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.

Comparative Transcriptome and Phytochemical Analysis Provides Insight into Triterpene Saponin Biosynthesis in Seeds and Flowers of the Tea Plant (Camellia sinensis)

Triterpene saponins exhibit various biological and pharmacological activities. However, the knowledge on saponin biosynthesis in tea plants (Camellia sinensis L.) is still limited. In this work, tea flower and seed samples at different developmental stages and leaves were collected and analyzed with UPLC-PDA-MS and RNA sequencing for saponin determination and transcriptome comparison. The saponin content reached around 19% in the freshly mature seeds and 7% in the green flower buds, and decreased with the fruit ripeness and flower blooming. Almost no saponins were detected in leaf samples. PCA and KEGG analysis suggested that the gene expression pattern and secondary metabolism in TF1 and TS2 vs. leaf samples were significantly different. Weighted gene coexpression network analysis (WGCNA) uncovered two modules related to saponin content. The mevalonate (MVA) instead of 2-C-methyl-d-erythritol-4-phospate (MEP) pathway was responsible for saponin accumulation in tea plants, and 3-hydroxy-3-methylglutaryl-CoA synthase (HMGS), diphosphomevalonate decarboxylase (MVD) and isopentenyl diphosphate isomerase (IDI) may be the key enzymes involved in saponin biosynthesis in tea seeds and flowers. Moreover, ten transcription factors (TFs) were predicted to regulate saponin biosynthesis in the tea plant. Taken together, our study provides a global insight into the saponin biosynthesis and accumulation in the tea plant.

Characterization of a Group of UDP-Glycosyltransferases Involved in the Biosynthesis of Triterpenoid Saponins of Panax notoginseng

UDP-glycosyltransferase (UGT)-mediated glycosylation is a common modification in triterpene saponins, which exhibit a wide range of bioactivities and important pharmacological effects. However, few UGTs involved in saponin biosynthesis have been identified, limiting the biosynthesis of saponins. In this study, an efficient heterologous expression system was established for evaluating the UGT-mediated glycosylation process of triterpene saponins. Six UGTs (UGTPn17, UGTPn42, UGTPn35, UGTPn87, UGTPn19, and UGTPn12) from Panax notoginseng were predicted and found to be responsible for efficient and direct enzymatic biotransformation of 21 triterpenoid saponins via 26 various glycosylation reactions. Among them, UGTPn87 exhibited promiscuous sugar-donor specificity of UDP-glucose (UDP-Glc) and UDP-xylose (UDP-Xyl) by catalyzing the elongation of the second sugar chain at the C3 or/and C20 sites of protopanaxadiol-type saponins with a UDP-Glc or UDP-Xyl donor, as well as at the C20 site of protopanaxadiol-type saponins with a UDP-Glc donor. Two new saponins, Fd-Xyl and Fe-Xyl, were generated by catalyzing the C3-O-Glc xylosylations of notoginsenoside Fd and notoginsenoside Fe when incubated with UGTPn87. Moreover, the complete biosynthetic pathways of 17 saponins were elucidated, among which notoginsenoside L, vinaginsenoside R16, gypenoside LXXV, and gypenoside XVII were revealed in Panax for the first time. A yeast cell factory was constructed with a yield of Rh2 at 354.69 mg/L and a glycosylation ratio of 60.40% in flasks. Our results reveal the biosynthetic pathway of a group of saponins in P. notoginseng and provide a theoretical basis for producing rare and valuable saponins, promoting their industrial application in medicine and functional foods.

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.

Structure and Function of Rhizosphere Soil and Root Endophytic Microbial Communities Associated With Root Rot of Panax notoginseng

Panax notoginseng (Burk.) F. H. Chen is a Chinese medicinal plant of the Araliaceae family used for the treatment of cardiovascular and cerebrovascular diseases in Asia. P. notoginseng is vulnerable to root rot disease, which reduces the yield of P. notoginseng. In this study, we analyzed the rhizosphere soil and root endophyte microbial communities of P. notoginseng from different geographical locations using high-throughput sequencing. Our results revealed that the P. notoginseng rhizosphere soil microbial community was more diverse than the root endophyte community. Rhodopseudomonas, Actinoplanes, Burkholderia, and Variovorax paradoxus can help P. notoginseng resist the invasion of root rot disease. Ilyonectria mors-panacis, Pseudomonas fluorescens, and Pseudopyrenochaeta lycopersici are pathogenic bacteria of P. notoginseng. The upregulation of amino acid transport and metabolism in the soil would help to resist pathogens and improve the resistance of P. notoginseng. The ABC transporter and gene modulating resistance genes can improve the disease resistance of P. notoginseng, and the increase in the number of GTs (glycosyltransferases) and GHs (glycoside hydrolases) families may be a molecular manifestation of P. notoginseng root rot. In addition, the complete genomes of two Flavobacteriaceae species and one Bacteroides species were obtained. This study demonstrated the microbial and functional diversity in the rhizosphere and root microbial community of P. notoginseng and provided useful information for a better understanding of the microbial community in P. notoginseng root rot. Our results provide insights into the molecular mechanism underlying P. notoginseng root rot and other plant rhizosphere microbial communities.

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.

Metabolomics analysis provides new insights into the medicinal value of flavonoids in tobacco leaves

Tobacco is a traditional Chinese medicine containing a variety of biologically active substances. In addition to being used to make cigarettes, tobacco is also a vastly underdeveloped medicinal resource. In order to identify and clarify the biological activities and medicinal value of tobacco leaves, the metabolomes of tobacco leaves were analyzed by liquid chromatography-tandem mass spectrometry (LC-MS/MS) based on multiple reaction monitoring (MRM). In total, 1169 metabolites were identified and quantified. The results showed that the metabolic profiles of the tobacco cultivars K326 and Yun87 are similar to each other but different from that of Hongda. Moreover, the curing process affects the metabolic profiles of tobacco leaves. Flavonoids are the largest class of metabolites in tobacco leaves. Flavonoids have multiple biological functions; for example, they can promote or inhibit inflammation. We found that quercetin provides anti-inflammatory activity by inhibiting the il-1β mRNA expression, while glycitin and neohesperidin can promote il-1β and il-6 production. Our results provide in-depth insights into the medical uses and biological mechanisms of tobacco leaves.

Transcriptome analysis identifies putative genes involved in triterpenoid biosynthesis in Platycodon grandiflorus

Comprehensive transcriptome analysis of different Platycodon grandiflorus tissues discovered genes related to triterpenoid saponin biosynthesis. Platycodon grandiflorus (Jacq.) A. DC. (P. grandiflorus), a traditional Chinese medicine, contains considerable triterpenoid saponins with broad pharmacological activities. Triterpenoid saponins are the major components of P. grandiflorus. Here, single-molecule real-time and next-generation sequencing technologies were combined to comprehensively analyse the transcriptome and identify genes involved in triterpenoid saponin biosynthesis in P. grandiflorus. We quantified four saponins in P. grandiflorus and found that their total content was highest in the roots and lowest in the stems and leaves. A total of 173,354 non-redundant transcripts were generated from the PacBio platform, and three full-length transcripts of β-amyrin synthase, the key synthase of β-amyrin, were identified. A total of 132,610 clean reads obtained from the DNBSEQ platform were utilised to explore key genes related to the triterpenoid saponin biosynthetic pathway in P. grandiflorus, and 96 differentially expressed genes were selected as candidates. The expression levels of these genes were verified by quantitative real-time PCR. Our reliable transcriptome data provide valuable information on the related biosynthesis pathway and may provide insights into the molecular mechanisms of triterpenoid saponin biosynthesis in P. grandiflorus.

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.

A Comparative Analysis on the Structure and Function of the Panax notoginseng Rhizosphere Microbiome

Panax notoginseng, an important Chinese medicinal herb, can be mainly cultivated in two planting patterns, cropland planting (DT) and understory planting (LX). We speculate that the rhizosphere microbiome may vary in DT and LX and may play an important role in promoting the growth and health of P. notoginseng. In the present study, culture-independent Illumina HiSeq was employed to investigate the rhizosphere bacteria and fungi under DT and LX planting patterns. Predominant phyla include Proteobacteria, Acidobacteria, Actinobacteria, Gemmatimonadetes, and Ascomycota in the two planting patterns. DT has higher alpha diversity index than LX. The predominant LX-core genera include Bradyrhizobium, Streptomyces, and Actinomadura, and the predominant DT-core genera include Sphingomonas, Variovorax, and Novosphingobium. Total relative abundance of the disease-suppression phylum (Proteobacteria, Firmicutes, and Actinobacteria) and the potential plant growth-promoting rhizobacteria (PGPR) were both significantly higher in LX than in DT. We also identified over-presented microbial functional traits mediating plant-microbe and microbe-microbe interactions, nutrition acquisition, and plant growth promotion in P. notoginseng rhizosphere. Our findings provide a valuable reference for studying beneficial microbes and pathogens of P. notoginseng planted in DT and LX.

De novo chromosome-level assembly of the Centella asiatica genome

Centella asiatica is a herbaceous, perennial species indigenous to India and Southeast Asia. C. asiatica possesses several medicinal properties: anti-aging, anti-inflammatory, wound healing and memory enhancing. The lack of available genomics resources significantly impedes the improvement of C. asiatica varieties through molecular breeding. Here, we combined the 10× Genomics linked-read technology and the long-range HiC technique to obtain the genome assembly. The final assembly contained nine pseudomolecules, corresponding to the haploid chromosome number in C. asiatica. These nine chromosomes covered 402,536,584 bases or 93.6% of the 430-Mb assembly. Comparative genomics analyses based on single-copy orthologous genes showed that C. asiatica and the common ancestor of Coriandrum sativum (coriander) and Daucus carota (carrot) diverged about 48 million years ago. This assembly provides a valuable reference genome for future molecular studies, varietal development through marker-assisted breeding and comparative genomics studies in C. asiatica.

Asterid Phylogenomics/Phylotranscriptomics Uncover Morphological Evolutionary Histories and Support Phylogenetic Placement for Numerous Whole-Genome Duplications

Asterids are one of the most successful angiosperm lineages, exhibiting extensive morphological diversity and including a number of important crops. Despite their biological prominence and value to humans, the deep asterid phylogeny has not been fully resolved, and the evolutionary landscape underlying their radiation remains unknown. To resolve the asterid phylogeny, we sequenced 213 transcriptomes/genomes and combined them with other data sets, representing all accepted orders and nearly all families of asterids. We show fully supported monophyly of asterids, Berberidopsidales as sister to asterids, monophyly of all orders except Icacinales, Aquifoliales, and Bruniales, and monophyly of all families except Icacinaceae and Ehretiaceae. Novel taxon placements benefited from the expanded sampling with living collections from botanical gardens, resolving hitherto uncertain relationships. The remaining ambiguous placements here are likely due to limited sampling and could be addressed in the future with relevant additional taxa. Using our well-resolved phylogeny as reference, divergence time estimates support an Aptian (Early Cretaceous) origin of asterids and the origin of all orders before the Cretaceous-Paleogene boundary. Ancestral state reconstruction at the family level suggests that the asterid ancestor was a woody terrestrial plant with simple leaves, bisexual, and actinomorphic flowers with free petals and free anthers, a superior ovary with a style, and drupaceous fruits. Whole-genome duplication (WGD) analyses provide strong evidence for 33 WGDs in asterids and one in Berberidopsidales, including four suprafamilial and seven familial/subfamilial WGDs. Our results advance the understanding of asterid phylogeny and provide numerous novel evolutionary insights into their diversification and morphological evolution.

Evolutionary Contribution of Duplicated Genes to Genome Evolution in the Ginseng Species Complex

Genes duplicated by whole genome duplication (WGD) and small-scale duplication (SSD) have played important roles in adaptive evolution of all flowering plants. However, it still remains underinvestigated how the distinct models of duplication events and their contending evolutionary patterns have shaped the genome and epigenomes of extant plant species. In this study, we investigated the contribution of the WGD- and SSD-derived duplicate genes to the genome evolution of one diploid and three closely related allotetraploid Panax species based on genome, methylome, and proteome data sets. Our genome-wide comparative analyses revealed that although the ginseng species complex was recently diverged, they have evolved distinct overall patterns of nucleotide variation, cytosine methylation, and protein-level expression. In particular, genetic and epigenetic asymmetries observed in the recent WGD-derived genes are largely consistent across the ginseng species complex. In addition, our results revealed that gene duplicates generated by ancient WGD and SSD mechanisms exhibited distinct evolutionary patterns. We found the ancient WGD-derived genes (i.e., ancient collinear gene) are genetically more conserved and hypomethylated at the cytosine sites. In contrast, some of the SSD-derived genes (i.e., dispersal duplicated gene) showed hypermethylation and high variance in nucleotide variation pattern. Functional enrichment analyses of the duplicated genes indicated that adaptation-related traits (i.e., photosynthesis) created during the distant ancient WGDs are further strengthened by both the more recent WGD and SSD. Together, our findings suggest that different types of duplicated genes may have played distinct but relaying evolutionary roles in the polyploidization and speciation processes in the ginseng species complex.

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.

Transcriptome analysis of Paris polyphylla var. yunnanensis illuminates the biosynthesis and accumulation of steroidal saponins in rhizomes and leaves

Paris polyphylla var. yunnanensis can synthesize Paris saponins with multiple effective therapies, and its rhizome has become an indispensable ingredient in many patented drugs. However, how Paris saponin content changes in tissues at different stages and the molecular mechanisms underlying the production and accumulation of the bioactive compounds are unclear. This study aimed to uncover the mechanisms underlying the biosynthesis and accumulation by integrating transcriptome sequencing and phytochemical investigation of the leaves and rhizomes at different growth stages. Paris saponin content in leaves was lower during the fruiting stage than the vegetative stage, whereas the content in rhizomes increased during the fruiting stage. The candidate genes related to Paris saponin biosynthesis were determined by transcriptome analyses. Most biosynthetic genes were found to be abundantly expressed in the leaves during the vegetative stage in the light of expression profiles and functional enrichment results. The expression patterns of the differentially expressed genes related to the biosynthesis were positively correlated with the accumulation of saponins in tissues. These findings suggest that both leaves and rhizomes are capable of biosynthesizing Paris saponins, and that aerial plant parts can be used to extract them. The different patterns of biosynthesis and accumulation in the leaves and rhizomes were also determined here. This study will help improve our understanding of the mechanisms underlying the biosynthesis and accumulation of Paris saponins, and aid in the comprehensive development and utilization of this medicinal plant.

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.

SMRT- and Illumina-based RNA-seq analyses unveil the ginsinoside biosynthesis and transcriptomic complexity in Panax notoginseng

Panax notoginseng is one of the most widely used traditional Chinese herbs with particularly valued roots. Triterpenoid saponins are mainly specialized secondary metabolites, which medically act as bioactive components. Knowledge of the ginsenoside biosynthesis in P. notoginseng, which is of great importance in the industrial biosynthesis and genetic breeding program, remains largely undetermined. Here we combined single molecular real time (SMRT) and Second-Generation Sequencing (SGS) technologies to generate a widespread transcriptome atlas of P. notoginseng. We mapped 2,383 full-length non-chimeric (FLNC) reads to adjacently annotated genes, corrected 1,925 mis-annotated genes and merged into 927 new genes. We identified 8,111 novel transcript isoforms that have improved the annotation of the current genome assembly, of which we found 2,664 novel lncRNAs. We characterized more alternative splicing (AS) events from SMRT reads (20,015 AS in 6,324 genes) than Illumina reads (18,498 AS in 9,550 genes), which contained a number of AS events associated with the ginsenoside biosynthesis. The comprehensive transcriptome landscape reveals that the ginsenoside biosynthesis predominantly occurs in flowers compared to leaves and roots, substantiated by levels of gene expression, which is supported by tissue-specific abundance of isoforms in flowers compared to roots and rhizomes. Comparative metabolic analyses further show that a total of 17 characteristic ginsenosides increasingly accumulated, and roots contained the most ginsenosides with variable contents, which are extraordinarily abundant in roots of the three-year old plants. We observed that roots were rich in protopanaxatriol- and protopanaxadiol-type saponins, whereas protopanaxadiol-type saponins predominated in aerial parts (leaves, stems and flowers). The obtained results will greatly enhance our understanding about the ginsenoside biosynthetic machinery in the genus Panax.

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.

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 Chromosome-Based Rubber Tree Genome Provides New Insights into Spurge Genome Evolution and Rubber Biosynthesis

The rubber tree, Hevea brasiliensis, produces natural rubber that serves as an essential industrial raw material. Here, we present a high-quality reference genome for a rubber tree cultivar GT1 using single-molecule real-time sequencing (SMRT) and Hi-C technologies to anchor the ∼1.47-Gb genome assembly into 18 pseudochromosomes. The chromosome-based genome analysis enabled us to establish a model of spurge chromosome evolution, since the common paleopolyploid event occurred before the split of Hevea and Manihot. We show recent and rapid bursts of the three Hevea-specific LTR-retrotransposon families during the last 10 million years, leading to the massive expansion by ∼65.88% (∼970 Mbp) of the whole rubber tree genome since the divergence from Manihot. We identify large-scale expansion of genes associated with whole rubber biosynthesis processes, such as basal metabolic processes, ethylene biosynthesis, and the activation of polysaccharide and glycoprotein lectin, which are important properties for latex production. A map of genomic variation between the cultivated and wild rubber trees was obtained, which contains ∼15.7 million high-quality single-nucleotide polymorphisms. We identified hundreds of candidate domestication genes with drastically lowered genomic diversity in the cultivated but not wild rubber trees despite a relatively short domestication history of rubber tree, some of which are involved in rubber biosynthesis. This genome assembly represents key resources for future rubber tree research and breeding, providing novel targets for improving plant biotic and abiotic tolerance and rubber production.

Bioinformatics-assisted, integrated omics studies on medicinal plants

The immense therapeutic and economic values of medicinal plants have attracted increasing attention from the worldwide researchers. It has been recognized that production of the authentic and high-quality herbal drugs became the prerequisite for maintaining the healthy development of the traditional medicine industry. To this end, intensive research efforts have been devoted to the basic studies, in order to pave a way for standardized authentication of the plant materials, and bioengineering of the metabolic pathways in the medicinal plants. In this paper, the recent advances of omics studies on the medicinal plants were summarized from several aspects, including phenomics and taxonomics, genomics, transcriptomics, proteomics and metabolomics. We proposed a multi-omics data-based workflow for medicinal plant research. It was emphasized that integration of the omics data was important for plant authentication and mechanistic studies on plant metabolism. Additionally, the computational tools for proper storage, efficient processing and high-throughput analyses of the omics data have been introduced into the workflow. According to the workflow, authentication of the medicinal plant materials should not only be performed at the phenomics level but also be implemented by genomic and metabolomic marker-based examination. On the other hand, functional genomics studies, transcriptional regulatory networks and protein-protein interactions will contribute greatly for deciphering the secondary metabolic pathways. Finally, we hope that our work could inspire further efforts on the bioinformatics-assisted, integrated omics studies on the medicinal plants.

Assembly and Annotation of a Draft Genome of the Medicinal Plant Polygonum cuspidatum

Polygonum cuspidatum (Japanese knotweed, also known as Huzhang in Chinese), a plant that produces bioactive components such as stilbenes and quinones, has long been recognized as important in traditional Chinese herbal medicine. To better understand the biological features of this plant and to gain genetic insight into the biosynthesis of its natural products, we assembled a draft genome of P. cuspidatum using Illumina sequencing technology. The draft genome is ca. 2.56 Gb long, with 71.54% of the genome annotated as transposable elements. Integrated gene prediction suggested that the P. cuspidatum genome encodes 55,075 functional genes, including 6,776 gene families that are conserved in the five eudicot species examined and 2,386 that are unique to P. cuspidatum. Among the functional genes identified, 4,753 are predicted to encode transcription factors. We traced the gene duplication history of P. cuspidatum and determined that it has undergone two whole-genome duplication events about 65 and 6.6 million years ago. Roots are considered the primary medicinal tissue, and transcriptome analysis identified 2,173 genes that were expressed at higher levels in roots compared to aboveground tissues. Detailed phylogenetic analysis demonstrated expansion of the gene family encoding stilbene synthase and chalcone synthase enzymes in the phenylpropanoid metabolic pathway, which is associated with the biosynthesis of resveratrol, a pharmacologically important stilbene. Analysis of the draft genome identified 7 abscisic acid and water deficit stress-induced protein-coding genes and 14 cysteine-rich transmembrane module genes predicted to be involved in stress responses. The draft de novo genome assembly produced in this study represents a valuable resource for the molecular characterization of medicinal compounds in P. cuspidatum, the improvement of this important medicinal plant, and the exploration of its abiotic stress resistance.

Transcriptome analysis of 1- and 3-year-old Panax notoginseng roots and functional characterization of saponin biosynthetic genes DS and CYP716A47-like

MAIN CONCLUSION

Transcriptome analysis revealed high expression of saponin biosynthetic genes may account for highly accumulated saponins in 3-year-old Panax notoginseng roots and DS and CYP716A47 - like were functionally verified by transgenic tobacco. Panax notoginseng is a well-known traditional medical herb that contains bioactive compounds known as saponins. Three major dammarene-type triterpene saponins including R1, Rb1, and Rg1 were found to be highly accumulated in the roots of 3-year-old plants when compared to those of 1-year-old plants. However, the underlying cellular mechanism is poorly understood. In this study, transcriptome analysis revealed that most genes involved in saponin biosynthesis in P. notoginseng roots augmented during their growth periods. The analysis of the KEGG pathway indicated that the primary metabolism, cell growth, and differentiation were less active in the roots of 3-year-old plant; however, secondary metabolisms were enhanced, thus providing molecular evidence for the harvesting of P. notoginseng roots in the 3rd year of growth. Furthermore, the functional role of DS and CYP716A47-like, two of the candidate genes involved in saponin biosynthesis isolated from P. notoginseng, were verified via overexpression in cultivated tobacco. Approximately, 0.325 µg g-1 of dammarenediol-II and 0.320 µg g-1 of protopanaxadiol were recorded in the dry leaves of transgenic tobacco overexpressed with DS and both DS and CYP716A47-like, respectively. This study provides insights into the molecular mechanisms for saponin accumulation in P. notoginseng roots during its growth period and paves a promising way to produce dammarenediol-II and protopanaxadiol via transgenic techniques.