Transcriptome analysis reveals ginsenosides biosynthetic genes, microRNAs and simple sequence repeats in Panax ginseng C. A. Meyer

Beyond genome: Advanced omics progress of Panax ginseng

Ginseng is a perennial herb of the genus Panax in the family Araliaceae as one of the most important traditional medicine. Genomic studies of ginseng assist in the systematic discovery of genes related to bioactive ginsenosides biosynthesis and resistance to stress, which are of great significance in the conservation of genetic resources and variety improvement. The transcriptome reflects the difference and consistency of gene expression, and transcriptomics studies of ginseng assist in screening ginseng differentially expressed genes to further explore the powerful gene source of ginseng. Protein is the ultimate bearer of ginseng life activities, and proteomic studies of ginseng assist in exploring the biosynthesis and regulation of secondary metabolites like ginsenosides and the molecular mechanism of ginseng adversity adaptation at the overall level. In this review, we summarize the current status of ginseng research in genomics, transcriptomics and proteomics, respectively. We also discuss and look forward to the development of ginseng genome allele mapping, ginseng spatiotemporal, single-cell transcriptome, as well as ginseng post-translational modification proteome. We hope that this review will contribute to the in-depth study of ginseng and provide a reference for future analysis of ginseng from a systems biology perspective.

Transcriptome-Wide Identification and Integrated Analysis of a UGT Gene Involved in Ginsenoside Ro Biosynthesis in Panax ginseng

Panax ginseng as a traditional medicinal plant with a long history of medicinal use. Ginsenoside Ro is the only oleanane-type ginsenoside in ginseng, and has various pharmacological activities, including anti-inflammatory, detoxification, and antithrombotic activities. UDP-dependent glycosyltransferase (UGT) plays a key role in the synthesis of ginsenoside, and the excavation of UGT genes involved in the biosynthesis of ginsenoside Ro has great significance in enriching ginsenoside genetic resources and further revealing the synthesis mechanism of ginsenoside. In this work, ginsenoside-Ro-synthesis-related genes were mined using the P. ginseng reference-free transcriptome database. Fourteen hub transcripts were identified by differential expression analysis and weighted gene co-expression network analysis. Phylogenetic and synteny block analyses of PgUGAT252645, a UGT transcript among the hub transcripts, showed that PgUGAT252645 belonged to the UGT73 subfamily and was relatively conserved in ginseng plants. Functional analysis showed that PgUGAT252645 encodes a glucuronosyltransferase that catalyzes the glucuronide modification of the C3 position of oleanolic acid using uridine diphosphate glucuronide as the substrate. Furthermore, the mutation at 622 bp of its open reading frame resulted in amino acid substitutions that may significantly affect the catalytic activity of the enzyme, and, as a consequence, affect the biosynthesis of ginsenoside Ro. Results of the in vitro enzyme activity assay of the heterologous expression product in E. coli of PgUGAT252645 verified the above analyses. The function of PgUGAT252645 was further verified by the result that its overexpression in ginseng adventitious roots significantly increased the content of ginsenoside Ro. The present work identified a new UGT gene involved in the biosynthesis of ginsenoside Ro, which not only enriches the functional genes in the ginsenoside synthesis pathway, but also provides the technical basis and theoretical basis for the in-depth excavation of ginsenoside-synthesis-related genes.

Carvacrol as a Stimulant of the Expression of Key Genes of the Ginsenoside Biosynthesis Pathway and Its Effect on the Production of Ginseng Saponins in Panax quinquefolium Hairy Root Cultures

The accumulation of ginsenosides (triterpenic saponins) was determined in Panax quinquefolium hairy root cultures subjected to an elicitation process using carvacrol at 5, 10, 25, 50, 100, 250, and 500 μM concentrations during 24 and 72 h exposure. This study was the first one in which carvacrol was applied as an elicitor. The content of eight ginsenosides, Rb1, Rb2, Rb3, Rc, Rd, Rg1, Rg2, and Re, was determined using HPLC analysis. Moreover, the quantitative RT-PCR method was applied to assess the relative expression level of farnesyl diphosphate synthase, squalene synthase, and dammarenediol synthase genes in the studied cultures. The addition of carvacrol (100 μM) was an effective approach to increase the production of ginsenosides. The highest content and productivity of all detected saponins were, respectively, 20.01 mg∙g-1 d.w. and 5.74 mg∙L-1∙day-1 after 72 h elicitation. The production profile of individual metabolites in P. quinquefolium cultures changed under the influence of carvacrol. The biosynthesis of most examined protopanaxadiol derivatives was reduced under carvacrol treatment. In contrast, the levels of ginsenosides belonging to the Rg group increased. The strongest effect of carvacrol was noticed for Re metabolites, achieving a 7.72-fold increase in comparison to the control. Saponin Rg2, not detected in untreated samples, was accumulated after carvacrol stimulation, reaching its maximum concentration after 72 h exposure to 10 μM elicitor.

Genetic and molecular dissection of ginseng (Panax ginseng Mey.) germplasm using high-density genic SNP markers, secondary metabolites, and gene expressions

Genetic and molecular knowledge of a species is crucial to its gene discovery and enhanced breeding. Here, we report the genetic and molecular dissection of ginseng, an important herb for healthy food and medicine. A mini-core collection consisting of 344 cultivars and landraces was developed for ginseng that represents the genetic variation of ginseng existing in its origin and diversity center. We sequenced the transcriptomes of all 344 cultivars and landraces; identified over 1.5 million genic SNPs, thereby revealing the genic diversity of ginseng; and analyzed them with 26,600 high-quality genic SNPs or a selection of them. Ginseng had a wide molecular diversity and was clustered into three subpopulations. Analysis of 16 ginsenosides, the major bioactive components for healthy food and medicine, showed that ginseng had a wide variation in the contents of all 16 ginsenosides and an extensive correlation of their contents, suggesting that they are synthesized through a single or multiple correlated pathways. Furthermore, we pair-wisely examined the relationships between the cultivars and landraces, revealing their relationships in gene expression, gene variation, and ginsenoside biosynthesis. These results provide new knowledge and new genetic and genic resources for advanced research and breeding of ginseng and related species.

Comparative transcriptional analysis of metabolic pathways and mechanisms regulating essential oil biosynthesis in four elite Cymbopogon spp

Cymbopogon is an important aromatic and medicinal grass with several species of ethnopharmaceutical importance. The genus is extremely rich in secondary metabolites, monoterpenes like geraniol and citral being principal constituents, also used as biomarker for classification and identification of Cymbopogon chemotypes. In the light of this, present study involved RNA sequencing and comparison of expression profiles of four contrasting Cymbopogon species namely C. flexuosus var. Chirharit (citral rich and frost resistant), C. martinii var. PRC-1 (geraniol rich), C. pendulus var. Praman (the most stable and citral-rich genotype), and Jamrosa (a hybrid of C. nardus var. confertiflorus × C. jwarancusa (rich in geraniol and geranyl acetate). The transcriptome profiles revealed marked differences in gene expression patterns of 28 differentially expressed genes (DEGs) of terpenoid metabolic pathways between the four Cymbopogon sp. The major DEGs were Carotenoid Cleavage Dioxygenases (CCD), Aspartate aminotransferase (ASP amino), Mevalonate E-4 hydroxy, AKR, GGPS, FDPS, and AAT. In addition, few TFs related to different regulatory pathways were also identified. The gene expression profiles of DEGs were correlated to the EO yield and their monoterpene compositions. Overall, the PRC-1 (C. martinii) shows distinguished gene expression profiles from all other genotypes. Thus, the transcriptome sequence database expanded our understanding of terpenoid metabolism and its molecular regulation in Cymbopogon species. Additionally, this data also serves as an important source of knowledge for enhancing oil yield and quality in Cymbopogon and closely related taxa. KEY MESSAGE: Unfolding the new secretes surrounding EO biosynthesis and regulation in four contrasting Cymbopogon species.

Transcriptome analysis and development of EST-SSR markers in Anoectochilus emeiensis (Orchidaceae)

Anoectochilus emeiensis K. Y. Lang, together with other Anoectochilus species, has long been used as the main source of many traditional Chinese medicines. Owing to the shortcomings of molecular markers, the study of the genetic diversity and medicinal component synthesis mechanism of the endemic Anoectochilus species has been delayed. In this study, I carried out a transcriptome analysis of A. emeiensis. A total of 78,381 unigenes were assembled from 64.2 million reads, and 47,541 (60.65%) unigenes were matched to known proteins in the public databases. Then, 9284 expressed sequence tag-derived simple sequence repeats (EST-SSRs) were identified, and the frequency of SSRs in the A. emeiensis transcriptome was 9.88%. Trinucleotide repeats (3699, 39.84%) were the most common type, followed by dinucleotide (3251, 35.02%) and mononucleotide (1750, 18.85%) repeats. Based on the SSR sequence, 6683 primer pairs were successfully designed, 40 primer pairs were randomly selected, and 10 primer pairs were identified as polymorphic loci from 186 individuals of A. emeiensis. The EST-SSR markers examined in this study will be informative for future population genetic studies of A. emeiensis.

Transcriptome and Phenotype Integrated Analysis Identifies Genes Controlling Ginsenoside Rb1 Biosynthesis and Reveals Their Interactions in the Process in Panax ginseng

Genes are the keys to deciphering the molecular mechanism underlying a biological trait and designing approaches desirable for plant genetic improvement. Ginseng is an important medicinal herb in which ginsenosides have been shown to be the major bioactive component; however, only a few genes involved in ginsenoside biosynthesis have been cloned through orthologue analysis. Here, we report the identification of 21 genes controlling Rb1 biosynthesis by stepwise ginseng transcriptome and Rb1 content integrated analysis. We first identified the candidate genes for Rb1 biosynthesis by integrated analysis of genes with the trait from four aspects, including gene transcript differential expression between highest- and lowest-Rb1 content cultivars, gene transcript expression-Rb1 content correlation, and biological impacts of gene mutations on Rb1 content, followed by the gene transcript co-expression network. Twenty-two candidate genes were identified, of which 21 were functionally validated for Rb1 biosynthesis by gene regulation, genetic transformation, and mutation analysis. These genes were strongly correlated in expression with the previously cloned genes encoding key enzymes for Rb1 biosynthesis. Based on the correlations, a pathway for Rb1 biosynthesis was deduced to indicate the roles of the genes in Rb1 biosynthesis. Moreover, the genes formed a strong co-expression network with the previously cloned Rb1 biosynthesis genes, and the variation in the network was associated with the variation in the Rb1 content. These results indicate that Rb1 biosynthesis is a process of correlative interactions among Rb1 biosynthesis genes. Therefore, this study provides new knowledge, 21 new genes, and 96 biomarkers for Rb1 biosynthesis useful for enhanced research and breeding in ginseng.

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.

Expression Profiles of Microsatellites in Fruit Tissues of Akebia trifoliata and Development of Efficient EST-SSR Markers

Akebia trifoliata, a member of the family Lardizabalaceae, has high exploitation potential for multiple economic purposes, so genetic improvements to meet requirements for commercial demand are needed. However, this progress is largely impeded by a lack of effective selection markers. In this study, we obtained 271.49 Gb of clean transcriptomic data from 12 samples (three tissues at four developmental stages) of A. trifoliata fruit. We identified 175,604, 194,370, and 207,906 SSRs from the de novo assembled 416,363, 463,756, and 491,680 unigene sequences obtained from the flesh, seed, and rind tissues, respectively. The profile and proportion of SSR motifs expressed in each fruit tissue and developmental stage were remarkably similar, but many trinucleotide repeats had differential expression levels among different tissues or at different developmental stages. In addition, we successfully designed 16,869 functional EST-SSR primers according to the annotated unigenes. Finally, 94 and 72 primer pairs out of 100 randomly selected primer pairs produced clear bands and polymorphic bands, respectively. These results were also used to elucidate the expression profiles of different tissues at various stages. Additionally, we provided a set of effective, polymorphic, and reliable EST-SSR markers sufficient for accelerating the discovery of metabolic and pathway-specific functional genes for genetic improvement and increased commercial productivity.

An insight into microRNA biogenesis and its regulatory role in plant secondary metabolism

The present review highlights the regulatory roles of microRNAs in plant secondary metabolism and focuses on different bioengineering strategies to modulate secondary metabolite content in plants. MicroRNAs (miRNAs) are the class of small endogenous, essential, non-coding RNAs that riboregulate the gene expression involved in various biological processes in most eukaryotes. MiRNAs has emerged as important regulators in plants that function by silencing target genes through cleavage or translational inhibition. These miRNAs plays an important role in a wide range of plant biological and metabolic processes, including plant development and various environmental response controls. Several important plant secondary metabolites like alkaloids, terpenoids, and phenolics are well studied for their function in plant defense against different types of pests and herbivores. Due to the presence of a wide range of biological and pharmaceutical properties of plant secondary metabolites, it is important to study the regulation of their biosynthetic pathways. The contribution of miRNAs in regulating plant secondary metabolism is not well explored. Recent advancements in molecular techniques have improved our knowledge in understanding the molecular function of genes, proteins, enzymes, and small RNAs involved in different steps of secondary metabolic pathways. In the present review, we have discussed the recent progress made on miRNA biogenesis, its regulation, and highlighted the current research developed in the field of identification, analysis, and characterizations of various miRNAs that regulate plant secondary metabolism. We have also discussed how different bioengineering strategies such as artificial miRNA (amiRNA), endogenous target mimicry, and CRISPR/Cas9 could be utilized to enhance the secondary metabolite production in plants.

SSR Loci Analysis in Transcriptome and Molecular Marker Development in Polygonatum sibiricum

To study the SSR loci information and develop molecular markers, a total of 435,858 unigenes in transcriptome of Polygonatum sibiricum were used to explore SSR. The distribution frequency of SSR and the basic characteristics of repeat motifs were analyzed using MISA software, and SSR primers were designed by Primer 3.0 software and then validated by PCR. Moreover, the gene function analysis of SSR Unigene was obtained by Blast. The results showed that 112,728 SSR loci were found in the transcriptome of Polygonatum sibiricum, which distributed in 435,858 unigenes with a distribution frequency of 25.86%. Mo-nucleotide and Di-nucleotide repeat were the main types, accounted for 83.83% of all SSRs. The repeat motifs of A/T and AC/GT were the predominant repeat types of Mo-nucleotide and Di-nucleotide, respectively. A total of 113,305 pairs of SSR primers with the potential to produce polymorphism were designed for maker development. One hundred and fifty-four of the 500 randomly selected primers not only produced fragments with expected molecular size but also had high polymorphism, which could accurately separate the tested varieties. The gene function of unigenes containing SSR was mostly related to the molecular function of Polygonatum sibiricum. The SSR markers in transcriptome of Polygonatum sibiricum show rich type, strong specificity, and high potential of polymorphism, which will benefit the candidate gene mining and marker-assisted breeding. The developed markers can also provide technical methods for molecular identification of intraspecific species of Polygonatum Mill. and maker-assisted breeding of superior varieties of Polygonatum Mill.

Ginsenosides in central nervous system diseases: Pharmacological actions, mechanisms, and therapeutics

The nervous system is one of the most complex physiological systems, and central nervous system diseases (CNSDs) are serious diseases that affect human health. Ginseng (Panax L.), the root of Panax species, are famous Chinese herbs that have been used for various diseases in China, Japan, and Korea since ancient times, and remain a popular natural medicine used worldwide in modern times. Ginsenosides are the main active components of ginseng, and increasing evidence has demonstrated that ginsenosides can prevent CNSDs, including neurodegenerative diseases, memory and cognitive impairment, cerebral ischemia injury, depression, brain glioma, multiple sclerosis, which has been confirmed in numerous studies. Therefore, this review summarizes the potential pathways by which ginsenosides affect the pathogenesis of CNSDs mainly including antioxidant effects, anti-inflammatory effects, anti-apoptotic effects, and nerve protection, which provides novel ideas for the treatment of CNSDs.

Chromosome-scale assembly and population diversity analyses provide insights into the evolution of Sapindus mukorossi

Sapindus mukorossi is an environmentally friendly plant and renewable energy source whose fruit has been widely used for biomedicine, biodiesel, and biological chemicals due to its richness in saponin and oil contents. Here, we report the first chromosome-scale genome assembly of S. mukorossi (covering ~391 Mb with a scaffold N50 of 24.66 Mb) and characterize its genetic architecture and evolution by resequencing 104 S. mukorossi accessions. Population genetic analyses showed that genetic diversity in the southwestern distribution area was relatively higher than that in the northeastern distribution area. Gene flow events indicated that southwest species may be the donor population for the distribution areas in China. Genome-wide selective sweep analysis showed that a large number of genes are involved in defense responses, growth and development, including SmRPS2, SmRPS4, SmRPS7, SmNAC2, SmNAC23, SmNAC102, SmWRKY6, SmWRKY26, and SmWRKY33. We also identified several candidate genes controlling six agronomic traits by genome-wide association studies, including SmPCBP2, SmbHLH1, SmCSLD1, SmPP2C, SmLRR-RKs, and SmAHP. Our study not only provides a rich genomic resource for further basic research on Sapindaceae woody trees but also identifies several economically significant genes for genomics-enabled improvements in molecular breeding.

Nucleotide Sequence Variation in Long-Term Tissue Cultures of Chinese Ginseng (Panax ginseng C. A. Mey.)

Plants have the salient biological property of totipotency, i.e., the capacity to regenerate a whole plant from virtually any kind of fully differentiated somatic cells after a process of dedifferentiation. This property has been well-documented by successful plant regeneration from tissue cultures of diverse plant species. However, the accumulation of somaclonal variation, especially karyotype alteration, during the tissue culture process compromises cell totipotency. In this respect, Chinese ginseng (Panax ginseng C. A. Mey.) is an exception in that it shows little decline in cell totipotency accompanied by remarkable chromosomal stability even after prolonged tissue cultures. However, it remains unclear whether chromosomal level stability necessarily couples with molecular genetic stability at the nucleotide sequence level, given that the two types of stabilities are generated by largely distinct mechanisms. Here, we addressed this issue by genome-wide comparisons at the single-base resolution of long-term tissue culture-regenerated P. ginseng plants. We identified abundant single nucleotide polymorphisms (SNPs) that have accumulated in cultured ginseng callus and are retained in the process of plant regeneration. These SNPs did not occur at random but showed differences among chromosomes and biased regional aggregation along a given chromosome. In addition, our results demonstrate that, compared with the overall genes, genes related to processes of cell totipotency and chromosomal stability possess lower mutation rates at both coding and flanking regions. In addition, collectively, the mutated genes exhibited higher expression levels than non-mutated genes and are significantly enriched in fundamental biological processes, including cellular component organization, development, and reproduction. These attributes suggest that the precipitated molecular level genetic variations during the process of regeneration in P. ginseng are likely under selection to fortify normal development. As such, they likely did not undermine chromosomal stability and totipotency of the long-term ginseng cultures.

Transcriptomic and Metabolomic Differences Between Two Saposhnikovia divaricata (Turcz.) Schischk Phenotypes With Single- and Double-Headed Roots

Saposhnikovia divaricata is derived from the dried roots of Saposhnikovia divaricata (Turcz.) Schischk and used as a Chinese herbal medicine for treating respiratory, immune, and nervous system diseases. The continuously increasing market demand for traditional Chinese medicine requires the commercial cultivation of Saposhnikovia divaricata using standardized methods and high yielding genotypes, such as double-headed root plants, for achieving consistent quality and a reliable supply. In this study, we aimed to identify the quantitative differences in chromone, a precursor of flavonoid biosynthesis, between plants with single- and double-headed roots using high-performance liquid chromatography and further explore the two phenotypes at the transcriptomic and metabolomic levels. Our results showed that the chromone content was significantly higher in plants with double-headed roots than in those with single-headed roots. Transcriptomic analysis revealed six significantly differentially expressed genes between the two phenotypes, including five key genes in the flavonoid biosynthesis pathway (4-coumarate-CoA ligase, chalcone synthase 1, vinorine synthase, chalcone-flavonone isomerase 1, and flavanone 3 beta-hydroxylase) and one key gene in the abscisic acid biosynthetic pathway (zeaxanthin epoxidase). Moreover, metabolomic analysis showed that the 126 differentially expressed metabolites were mainly enriched in the biosynthesis of secondary metabolites and phytohormones. Overall, our results suggest that plants with double-headed roots have higher medicinal value than those with single-headed roots, probably due to differences in various biosynthetic pathways. These data might help select the genotypes with superior yield and therapeutic properties.

Bm-miR172c-5p Regulates Lignin Biosynthesis and Secondary Xylem Thickness by Altering the Ferulate 5 Hydroxylase Gene in Bacopa monnieri

MicroRNAs (miRNAs) are small non-coding, endogenous RNAs containing 20-24 nucleotides that regulate the expression of target genes involved in various plant processes. A total of 1,429 conserved miRNAs belonging to 95 conserved miRNA families and 12 novel miRNAs were identified from Bacopa monnieri using small RNA sequencing. The Bm-miRNA target transcripts related to the secondary metabolism were further selected for validation. The Bm-miRNA expression in shoot and root tissues was negatively correlated with their target transcripts. The Bm-miRNA cleavage sites were mapped within the coding or untranslated region as depicted by the modified RLM-RACE. In the present study, we validate three miRNA targets, including asparagine synthetase, cycloartenol synthase and ferulate 5 hydroxylase (F5H) and elucidate the regulatory role of Bm-miR172c-5p, which cleaves the F5H gene involved in the lignin biosynthesis. Overexpression (OE) of Bm-miR172c-5p precursor in B. monnieri suppresses F5H gene, leading to reduced lignification and secondary xylem thickness under control and drought stress. By contrast, OE of endogenous target mimics (eTMs) showed enhanced lignification and secondary xylem thickness leading to better physiological response under drought stress. Taken together, we suggest that Bm-miRNA172c-5p might be a key player in maintaining the native phenotype of B. monnieri under control and different environmental conditions.

MicroRNAs Roles in Plants Secondary Metabolism

Plant growth and development is dependent on the regulation of classes of microRNAs (miRNAs) that have emerged as important gene regulators. These miRNAs can regulate plant gene expression to function. They play an important roles in biological homeostasis and environmental response controls. A wide range of plant biological and metabolic processes, including developmental timing, tissues specific development, and differentiation, depends on miRNAs. They perpetually regulate secondary metabolite functions in different plant family lines. Mapping of molecular phylogenies shows the distribution of secondary metabolism in the plant territory. More importantly, a lot of information related to miRNA regulatory processes in plants is revealed, but the role of miRNAs in secondary metabolism regulation and functions of the metabolites are still unclear. In this review, we pinnacle some potential miRNAs regulating the secondary metabolite biosynthesis activities in plants. This will provide an alternative knowledge for functional studies of secondary metabolism.

Advance in glycosyltransferases, the important bioparts for production of diversified ginsenosides

Ginsenosides are a series of glycosylated triterpenoids predominantly originated from Panax species with multiple pharmacological activities such as anti-aging, mediatory effect on the immune system and the nervous system. During the biosynthesis of ginsenosides, glycosyltransferases play essential roles by transferring various sugar moieties to the sapogenins in contributing to form structure and bioactivity diversified ginsenosides, which makes them important bioparts for synthetic biology-based production of these valuable ginsenosides. In this review, we summarized the functional elucidated glycosyltransferases responsible for ginsenoside biosynthesis, the advance in the protein engineering of UDP-glycosyltransferases (UGTs) and their application with the aim to provide in-depth understanding on ginsenoside-related UGTs for the production of rare ginsenosides applying synthetic biology-based microbial cell factories in the future.

Ginseng Omics for Ginsenoside Biosynthesis

Ginseng, also known as the king of herbs, has been regarded as an important traditional medicine for several millennia. Ginsenosides, a group of triterpenoid saponins, have been characterized as bioactive compounds of ginseng. The complexity of ginsenosides hindered ginseng research and development both in cultivation and clinical research. Therefore, deciphering the ginsenoside biosynthesis pathway has been a focus of interest for researchers worldwide. The new emergence of biological research tools consisting of omics and bioinformatic tools or computational biology tools are the research trend in the new century. Ginseng is one of the main subjects analyzed using these new quantification tools, including tools of genomics, transcriptomics, and proteomics. Here, we review the current progress of ginseng omics research and provide results for the ginsenoside biosynthesis pathway. Organization and expression of the entire pathway, including the upstream MVA pathway, the cyclization of ginsenoside precursors, and the glycosylation process, are illustrated. Regulatory gene families such as transcriptional factors and transporters are also discussed in this review.

Noncoding RNAs in Medicinal Plants and their Regulatory Roles in Bioactive Compound Production

BACKGROUND

Noncoding RNAs (ncRNAs), such as microRNAs (miRNAs), small interfering RNAs (siRNAs) and long noncoding RNAs (lncRNAs), play significant regulatory roles in plant development and secondary metabolism and are involved in plant response to biotic and abiotic stresses. They have been intensively studied in model systems and crops for approximately two decades and massive amount of information have been obtained. However, for medicinal plants, ncRNAs, particularly their regulatory roles in bioactive compound biosynthesis, are just emerging as a hot research field.

OBJECTIVE

This review aims to summarize current knowledge on herbal ncRNAs and their regulatory roles in bioactive compound production.

RESULTS

So far, scientists have identified thousands of miRNA candidates from over 50 medicinal plant species and 11794 lncRNAs from Salvia miltiorrhiza, Panax ginseng, and Digitalis purpurea. Among them, more than 30 miRNAs and five lncRNAs have been predicted to regulate bioactive compound production.

CONCLUSION

The regulation may achieve through various regulatory modules and pathways, such as the miR397-LAC module, the miR12112-PPO module, the miR156-SPL module, the miR828-MYB module, the miR858-MYB module, and other siRNA and lncRNA regulatory pathways. Further functional analysis of herbal ncRNAs will provide useful information for quality and quantity improvement of medicinal plants.

Transcriptomic Analysis, Cloning, Characterization, and Expression Analysis of Triterpene Biosynthetic Genes and Triterpene Accumulation in the Hairy Roots of Platycodon grandiflorum Exposed to Methyl Jasmonate

Platycodon grandiflorum is a perennial plant that has been used for medicinal purposes. Specifically, it is widely used in Northern China and Korea for the treatment of various diseases. Terpenoids belong to a group called secondary metabolites and have attracted a wide range of interest. Here, we determined the expressed sequence tag (EST) library of the methyl jasmonate (MeJA)-treated hairy root of P. grandiflorum. In total, 5760 ESTs were obtained, but after deleting the vector sequences and removing poor-quality sequences, a total of 2536 ESTs were attained. Of these, 811 contigs and 1725 singletons were annotated. The data were further analyzed with a focus on the gene families of the terpenoid biosynthetic pathway (TBP). We identified and characterized four TBP genes; among these were three full-length cDNAs encoding PgHMGS, PgMK, and PgMVD, whereas PgHMGR had a partial sequence based on the EST library database. Phylogenetic analysis and a pairwise identity matrix showed that these identified genes were closely related to those of other higher plants. Moreover, the tertiary structure and multiple alignment analysis showed that several distinct conserved motifs were present. Quantitative reverse transcription-polymerase chain reaction results revealed that TBP genes were constitutively expressed in all organs of P. grandiflorum, while the expression of transcript levels of these genes varied distinctly among the organs. Additionally, the total amount of platycosides was highly detected in the root, accumulating in concentrations 7.8 and 2.6 times higher than in the hairy root and stem, respectively, and 1.4 times higher than in the leaf and flower. The highest amount of total phytosterols was found to accumulate in the leaves at 9.3, 9.1, 1.8, and 1.6 times higher than that of the stem, root, hairy root, and flower, respectively. After the hairy root was exposed to the MeJA treatment, the transcript levels of PgHMGS, PgHMGR, PgMK, and PgMVD had significantly increased. The highest level of transcript accumulation occurred at 3 h after initial exposure for most of the genes. Meanwhile, triterpene saponin accumulation increased with the increase in the time of exposure, and at 48 h after initial exposure, the total saponin content was the highest recorded.

Azadirachta indica MicroRNAs: Genome-Wide Identification, Target Transcript Prediction, and Expression Analyses

MicroRNAs are short, endogenous, non-coding RNAs, liable for essential regulatory function. Numerous miRNAs have been identified and studied in plants with known genomic or small RNA resources. Despite the availability of genomic and transcriptomic resources, the miRNAs have not been reported in the medicinal tree Azadirachta indica (Neem) till date. Here for the first time, we report extensive identification of miRNAs and their possible targets in A. indica which might help to unravel their therapeutic potential. A comprehensive search of miRNAs in the A. indica genome by C-mii tool was performed. Overall, 123 miRNAs classified into 63 families and their stem-loop hairpin structures were predicted. The size of the A. indica (ain)-miRNAs ranged between 19 and 23 nt in length, and their corresponding ain-miRNA precursor sequence MFEI value averaged as -1.147 kcal/mol. The targets of ain-miRNAs were predicted in A. indica as well as Arabidopsis thaliana plant. The gene ontology (GO) annotation revealed the involvement of ain-miRNA targets in developmental processes, transport, stress, and metabolic processes including secondary metabolism. Stem-loop qRT-PCR was carried out for 25 randomly selected ain-miRNAs and differential expression patterns were observed in different A. indica tissues. Expression of miRNAs and its targets shows negative correlation in a dependent manner.

De novo RNA sequencing and analysis reveal the putative genes involved in diterpenoid biosynthesis in Aconitum vilmorinianum roots

In this study, the putative genes involved in diterpenoid alkaloids biosynthesis in A. vilmorinianum roots were revealed by transcriptome sequencing. 59.39 GB of clean bases and 119,660 unigenes were assembled, of which 69,978 unigenes (58.48%) were annotated. We identified 27 classes of genes (139 candidate genes) involved in the synthesis of diterpenoid alkaloids, including the mevalonate (MVA) pathway, the methylerythritol 4-phosphate (MEP) pathway, the farnesyl diphosphate regulatory pathway, and the diterpenoid scaffold synthetic pathway. 12 CYP450 genes were identified. We found that hydroxymethylglutaryl-CoA reductase was the key enzyme in MVA metabolism, which was regulated by miR6300. Transcription factors, such as bHLH, AP2/EREBP, and MYB, used to synthesize the diterpenes were analyzed.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02646-6.

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.

Comparative transcriptome analysis of root, stem, and leaf tissues of Entada phaseoloides reveals potential genes involved in triterpenoid saponin biosynthesis

Background

Entada phaseoloides (L.) Merr. is an important traditional medicinal plant. The stem of Entada phaseoloides is popularly used as traditional medicine because of its significance in dispelling wind and dampness and remarkable anti-inflammatory activities. Triterpenoid saponins are the major bioactive compounds of Entada phaseoloides. However, genomic or transcriptomic technologies have not been used to study the triterpenoid saponin biosynthetic pathway in this plant.

Results

We performed comparative transcriptome analysis of the root, stem, and leaf tissues of Entada phaseoloides with three independent biological replicates and obtained a total of 53.26 Gb clean data and 116,910 unigenes, with an average N50 length of 1218 bp. Putative functions could be annotated to 42,191 unigenes (36.1%) based on BLASTx searches against the Non-redundant, Uniprot, KEGG, Pfam, GO, KEGG and COG databases. Most of the unigenes related to triterpenoid saponin backbone biosynthesis were specifically upregulated in the stem. A total of 26 cytochrome P450 and 17 uridine diphosphate glycosyltransferase candidate genes related to triterpenoid saponin biosynthesis were identified. The differential expressions of selected genes were further verified by qPT-PCR.

Conclusions

The dataset reported here will facilitate the research about the functional genomics of triterpenoid saponin biosynthesis and genetic engineering of Entada phaseoloides.

Discovery and modification of cytochrome P450 for plant natural products biosynthesis

Cytochrome P450s are widespread in nature and play key roles in the diversification and functional modification of plant natural products. Over the last few years, there has been remarkable progress in plant P450s identification with the rapid development of sequencing technology, "omics" analysis and synthetic biology. However, challenges still persist in respect of crystal structure, heterologous expression and enzyme engineering. Here, we reviewed several research hotspots of P450 enzymes involved in the biosynthesis of plant natural products, including P450 databases, gene mining, heterologous expression and protein engineering.

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.

De novo transcriptome of Gymnema sylvestre identified putative lncRNA and genes regulating terpenoid biosynthesis pathway

Gymnema sylvestre is a highly valuable medicinal plant in traditional Indian system of medicine and used in many polyherbal formulations especially in treating diabetes. However, the lack of genomic resources has impeded its research at molecular level. The present study investigated functional gene profile of G. sylvestre via RNA sequencing technology. The de novo assembly of 88.9 million high quality reads yielded 23,126 unigenes, of which 18116 were annotated against databases such as NCBI nr database, gene ontology (GO), KEGG, Pfam, CDD, PlantTFcat, UniProt & GreeNC. Total 808 unigenes mapped to 78 different Transcription Factor families, whereas 39 unigenes assigned to CYP450 and 111 unigenes coding for enzymes involved in the biosynthesis of terpenoids including transcripts for synthesis of important compounds like Vitamin E, beta-amyrin and squalene. Among them, presence of six important enzyme coding transcripts were validated using qRT-PCR, which showed high expression of enzymes involved in methyl-erythritol phosphate (MEP) pathway. This study also revealed 1428 simple sequence repeats (SSRs), which may aid in molecular breeding studies. Besides this, 8 putative long non-coding RNAs (lncRNAs) were predicted from un-annotated sequences, which may hold key role in regulation of essential biological processes in G. sylvestre. The study provides an opportunity for future functional genomic studies and to uncover functions of the lncRNAs in G. sylvestre.

SNP discovery and functional annotation in the Panax japonicus var. major transcriptome

Due to the lack of a Panax japonicus var. major reference genome, we assembled a reference transcriptome from P. japonicus C. A. Mey transcriptome sequencing data, and 203 283 unigenes were obtained. In this study, with the assistance from the Trinity, Bowtie2 and SAMtools softwares, 218 465 single nucleotide polymorphisms (SNPs) were identified by mapping the Illumina sequences to the reference transcriptome. The SNP forms included 126 262 transformations and 92 203 transversions. A large number of SNP loci were associated with triterpenoid saponin synthesis: 54 SNPs were associated with cytochrome P450, one with glycosyl transferase and 94 with the biosynthesis of the triterpenoid saponin backbone.

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.

Triterpenoid-biosynthetic UDP-glycosyltransferases from plants

Triterpenoid saponins are naturally occurring structurally diverse glycosides of triterpenes that are widely distributed among plant species. Great interest has been expressed by pharmaceutical and agriculture industries for the glycosylation of triterpenes. Such modifications alter their taste and bio-absorbability, affect their intra-/extracellular transport and storage in plants, and induce novel biological activities in the human body. Uridine diphosphate (UDP)-glycosyltransferases (UGTs) catalyze glycosylation using UDP sugar donors. These enzymes belong to a multigene family and recognize diverse natural products, including triterpenes, as the acceptor molecules. For this review, we collected and analyzed all of the UGT sequences found in Arabidopsis thaliana as well as 31 other species of triterpene-producing plants. To identify potential UGTs with novel functions in triterpene glycosylation, we screened and classified those candidates based on similarity with UGTs from Panax ginseng, Glycine max, Medicago truncatula, Saponaria vaccaria, and Barbarea vulgaris that are known to function in glycosylate triterpenes. We highlight recent findings on UGT inducibility by methyl jasmonate, tissue-specific expression, and subcellular localization, while also describing their catalytic activity in terms of regioselectivity for potential key UGTs dedicated to triterpene glycosylation in plants. Discovering these new UGTs expands our capacity to manipulate the biological and physicochemical properties of such valuable molecules.

The inhibitory effects of ginsenoside Rd on the human glioma U251 cells and its underlying mechanisms

OBJECTIVE

The current study was designed to investigate the inhibitory effects of ginsenoside Rd (Gs-Rd) on human glioma U251 cells in vitro and its possible underlying mechanisms.

METHODS

The groups included blank control group, low concentration Gs-Rd treatment group (20 μM), mid concentration Gs-Rd treatment group (40 μM), and high concentration Gs-Rd treatment group (80 μM). The proliferative activity of human glioma U251 cells was detected by the MTT assay. Flow cytometry was performed to measure cell apoptosis of human glioma U251 cells. In addition, the ELISA assay was used to measure the telomerase activities in different groups on 24 hours, 48 hours, and 72 hours. Furthermore, real-time quantitative polymerase chain reaction (RT-PCR) and Western blot analysis were performed to measure the expression of Bcl-2, human telomerase catalytic subunit (hTERT), and caspase-3 in different groups on 48 hours at both messenger RNA (mRNA) and protein levels.

RESULTS

The proliferation of U251 cells was inhibited by Gs-Rd with different concentrations in the dose- and time-dependent manners. In addition, Gs-Rd promoted U251 cell apoptosis rate in a dose-dependent manner. Gs-Rd with different concentrations (20 μM, 40 μM, and 80 μM) significantly enhanced the expression of teleomerase on 24 hours and 48 hours. In addition, Gs-Rd with different concentrations significantly increased caspase-3 and decreased Bcl-2 and hTERT expressions at both mRNA and protein levels.

CONCLUSION

The Gs-Rd can remarkably inhibit the proliferation and promote cell apoptosis of human glioma U251 cells. The possible underlying mechanisms could be related to inhibiting telomerase activity, downregulating expression of Bcl-2 and hTERT, and upregulating expression of caspase-3 of human glioma U251 cells.

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.

Identification of candidate UDP-glycosyltransferases involved in protopanaxadiol-type ginsenoside biosynthesis in Panax ginseng

Ginsenosides are dammarane-type or triterpenoidal saponins that contribute to the various pharmacological activities of the medicinal herb Panax ginseng. The putative biosynthetic pathway for ginsenoside biosynthesis is known in P. ginseng, as are some of the transcripts and enzyme-encoding genes. However, few genes related to the UDP-glycosyltransferases (UGTs), enzymes that mediate glycosylation processes in final saponin biosynthesis, have been identified. Here, we generated three replicated Illumina RNA-Seq datasets from the adventitious roots of P. ginseng cultivar Cheongsun (CS) after 0, 12, 24, and 48 h of treatment with methyl jasmonate (MeJA). Using the same CS cultivar, metabolomic data were also generated at 0 h and every 12-24 h thereafter until 120 h of MeJA treatment. Differential gene expression, phylogenetic analysis, and metabolic profiling were used to identify candidate UGTs. Eleven candidate UGTs likely to be involved in ginsenoside glycosylation were identified. Eight of these were considered novel UGTs, newly identified in this study, and three were matched to previously characterized UGTs in P. ginseng. Phylogenetic analysis further asserted their association with ginsenoside biosynthesis. Additionally, metabolomic analysis revealed that the newly identified UGTs might be involved in the elongation of glycosyl chains of ginsenosides, especially of protopanaxadiol (PPD)-type ginsenosides.

Ginsenoside Rb1 administration attenuates focal cerebral ischemic reperfusion injury through inhibition of HMGB1 and inflammation signals

High-mobility group box 1 (HMGB1) is released after focal cerebral ischemia/reperfusion (I/R), and aggravates brain tissue damage. Ginsenoside Rb1 (Rb1), isolated from Panax ginseng, has been reported to inhibit I/R-induced cell death in the brain. The present study aimed to investigate the protective ability of GRb1 on focal cerebral I/R rats and to explore its further mechanisms. A middle cerebral artery occlusion (MCAO) rat model was established and treated with different doses of Rb1. The neurological deficits were examined after reperfusion, and TTC staining was applied to assess the infarct volume. Histology and TUNEL staining were performed to evaluate pathological changes and neuronal cell apoptosis in brain tissues. HMGB1 and levels of inflammatory factors and proteins, were examined by ELISA or western blotting. Rb1 treatment notably improved the neurological deficits in an MCAO model, accompanied by decreased infarct volume in the brain tissues. Histological examination revealed that the necrotic tissue area in MCAO rats was also diminished by Rb1 treatment. Apoptosis induced by cerebral I/R was also attenuated by Rb1 treatment via downregulation of cleaved caspase-3 and caspase-9 levels. HMGB1 release was inhibited by Rb1 treatment in MCAO rats, and the levels of nuclear factor-κB, tumor necrosis factor-α, interleukin-6, inducible nitric oxide synthase and nitric oxide were also decreased. The present study suggests that Rb1 serves a protective role in I/R-induced cerebral-neuron injury, due to the decreased cerebral infarct volume of brain tissue. The mechanisms underlying these effects may be associated with the inhibition of HMGB1 inflammatory signals.

Integrated Chemical and Transcriptomic Analysis Reveals the Distribution of Protopanaxadiol- and Protopanaxatriol-Type Saponins in Panax notoginseng

Panax notoginseng is famous for its important therapeutic effects and commonly used worldwide. The active ingredients saponins have distinct contents in different tissues of P. notoginseng, and they may be related to the expression of key genes in the synthesis pathway. In our study, high-performance liquid chromatography results indicated that the contents of protopanaxadiol-(Rb1, Rc, Rb2, and Rd) and protopanaxatriol-type (R1, Rg1, and Re) saponins in below ground tissues were higher than those in above ground tissues. Clustering dendrogram and PCA analysis suggested that the below and above ground tissues were clustered into two separate groups. A total of 482 and 882 unigenes were shared in the below and above ground tissues, respectively. A total of 75 distinct expressions of CYPs transcripts (RPKM ≥ 10) were detected. Of these transcripts, 38 and 37 were highly expressed in the below ground and above ground tissues, respectively. RT-qPCR analysis showed that CYP716A47 gene was abundantly expressed in the above ground tissues, especially in the flower, whose expression was 31.5-fold higher than that in the root. CYP716A53v2 gene was predominantly expressed in the below ground tissues, especially in the rhizome, whose expression was 20.1-fold higher than that in the flower. Pearson's analysis revealed that the CYP716A47 expression was significantly correlated with the contents of ginsenoside Rc and Rb2. The CYP716A53v2 expression was associated with the saponin contents of protopanaxadiol-type (Rb1 and Rd) and protopanaxatriol-type (R1, Rg1, and Re). Results indicated that the expression patterns of CYP716A47 and CYP716A53v2 were correlated with the distribution of protopanaxadiol-type and protopanaxatriol-type saponins in P. notoginseng. This study identified the pivotal genes regulating saponin distribution and provided valuable information for further research on the mechanisms of saponin synthesis, transportation, and accumulation.

Panax ginseng genome examination for ginsenoside biosynthesis

Ginseng, which contains ginsenosides as bioactive compounds, has been regarded as an important traditional medicine for several millennia. However, the genetic background of ginseng remains poorly understood, partly because of the plant's large and complex genome composition. We report the entire genome sequence of Panax ginseng using next-generation sequencing. The 3.5-Gb nucleotide sequence contains more than 60% repeats and encodes 42 006 predicted genes. Twenty-two transcriptome datasets and mass spectrometry images of ginseng roots were adopted to precisely quantify the functional genes. Thirty-one genes were identified to be involved in the mevalonic acid pathway. Eight of these genes were annotated as 3-hydroxy-3-methylglutaryl-CoA reductases, which displayed diverse structures and expression characteristics. A total of 225 UDP-glycosyltransferases (UGTs) were identified, and these UGTs accounted for one of the largest gene families of ginseng. Tandem repeats contributed to the duplication and divergence of UGTs. Molecular modeling of UGTs in the 71st, 74th, and 94th families revealed a regiospecific conserved motif located at the N-terminus. Molecular docking predicted that this motif captures ginsenoside precursors. The ginseng genome represents a valuable resource for understanding and improving the breeding, cultivation, and synthesis biology of this key herb.

Ginseng Genome Database: an open-access platform for genomics of Panax ginseng

BACKGROUND

The ginseng (Panax ginseng C.A. Meyer) is a perennial herbaceous plant that has been used in traditional oriental medicine for thousands of years. Ginsenosides, which have significant pharmacological effects on human health, are the foremost bioactive constituents in this plant. Having realized the importance of this plant to humans, an integrated omics resource becomes indispensable to facilitate genomic research, molecular breeding and pharmacological study of this herb.

DESCRIPTION

The first draft genome sequences of P. ginseng cultivar "Chunpoong" were reported recently. Here, using the draft genome, transcriptome, and functional annotation datasets of P. ginseng, we have constructed the Ginseng Genome Database http://ginsengdb.snu.ac.kr /, the first open-access platform to provide comprehensive genomic resources of P. ginseng. The current version of this database provides the most up-to-date draft genome sequence (of approximately 3000 Mbp of scaffold sequences) along with the structural and functional annotations for 59,352 genes and digital expression of genes based on transcriptome data from different tissues, growth stages and treatments. In addition, tools for visualization and the genomic data from various analyses are provided. All data in the database were manually curated and integrated within a user-friendly query page.

CONCLUSION

This database provides valuable resources for a range of research fields related to P. ginseng and other species belonging to the Apiales order as well as for plant research communities in general. Ginseng genome database can be accessed at http://ginsengdb.snu.ac.kr /.

The mRNA and miRNA transcriptomic landscape of Panax ginseng under the high ambient temperature

BACKGROUND

Ginseng is a popular traditional herbal medicine in north-eastern Asia. It has been used for human health for over thousands of years. With the rise in global temperature, the production of Korean ginseng (Panax ginseng C.A.Meyer) in Korea have migrated from mid to northern parts of the Korean peninsula to escape from the various higher temperature related stresses. Under the high ambient temperature, vegetative growth was accelerated, which resulted in early flowering. This precocious phase change led to yield loss. Despite of its importance as a traditional medicine, biological mechanisms of ginseng has not been well studied and even the genome sequence of ginseng is yet to be determined due to its complex genome structure. Thus, it is challenging to investigate the molecular biology mechanisms at the transcript level.

RESULTS

To investigate how ginseng responds to the high ambient temperature environment, we performed high throughput RNA sequencing and implemented a bioinformatics pipeline for the integrated analysis of small-RNA and mRNA-seq data without a reference genome. By performing reverse transcriptase (RT) PCR and sanger sequencing of transcripts that were assembled using our pipeline, we validated that their sequences were expressed in our samples. Furthermore, to investigate the interaction between genes and non-coding small RNAs and their regulation status under the high ambient temperature, we identified potential gene regulatory miRNAs. As a result, 100,672 contigs with significant expression level were identified and 6 known, 214 conserved and 60 potential novel miRNAs were predicted to be expressed under the high ambient temperature.

CONCLUSION

Collectively, we have found that development, flowering and temperature responsive genes were induced under high ambient temperature, whereas photosynthesis related genes were repressed. Functional miRNAs were down-regulated under the high ambient temperature. Among them are miR156 and miR396 that target flowering (SPL6/9) and growth regulating genes (GRF) respectively.

De novo transcriptomic analysis of leaf and fruit tissue of Cornus officinalis using Illumina platform

Cornus officinalis is one of the most widely used medicinal plants in China and other East Asian countries to cure diseases such as liver, kidney, cardiovascular diseases and frequent urination for thousands of years. It is a Level 3 protected species, and is one of the 42 national key protected wild species of animals and plants in China. However, the genetics and molecular biology of C. officinalis are poorly understood, which has hindered research on the molecular mechanism of its metabolism and utilization. Hence, enriching its genomic data and information is very important. In recent years, the fast-growing technology of next generation sequencing has provided an effective path to gain genomic information from nonmodel species. This study is the first to explore the leaf and fruit tissue transcriptome of C. officinalis using the Illumina HiSeq 4000 platform. A total of 57,954,134 and 60,971,652 clean reads from leaf and fruit were acquired, respectively (GenBank number SRP115440). The pooled reads from all two libraries were assembled into 56,392 unigenes with an average length 856 bp. Among these, 41,146 unigenes matched with sequences in the NCBI nonredundant protein database. The Gene Ontology database assigned 24,336 unigenes with biological process (83.26%), cellular components (53.58%), and molecular function (83.93%). In addition, 10,808 unigenes were assigned a KOG functional classification by the KOG database. Searching against the KEGG pathway database indicated that 18,435 unigenes were mapped to 371 KEGG pathways. Moreover, the edgeR database identified 4,585 significant differentially expressed genes (DEGs), of which 1,392 were up-regulated and 3,193 were down-regulated in fruit tissue compared with leaf tissue. Finally, we explored 581 transcription factors with 50 transcription factor gene families. Most DEGs and transcription factors were related to terpene biosynthesis and secondary metabolic regulation. This study not only represented the first de novo transcriptomic analysis of C. officinalis but also provided fundamental information on its genes and biosynthetic pathway. These findings will help us explore the molecular metabolism mechanism of terpene biosynthesis in C. officinalis.

Isoform Sequencing Provides a More Comprehensive View of the Panax ginseng Transcriptome

Korean ginseng (Panax ginseng C.A. Meyer) has been widely used for medicinal purposes and contains potent plant secondary metabolites, including ginsenosides. To obtain transcriptomic data that offers a more comprehensive view of functional genomics in P. ginseng, we generated genome-wide transcriptome data from four different P. ginseng tissues using PacBio isoform sequencing (Iso-Seq) technology. A total of 135,317 assembled transcripts were generated with an average length of 3.2 kb and high assembly completeness. Of those unigenes, 67.5% were predicted to be complete full-length (FL) open reading frames (ORFs) and exhibited a high gene annotation rate. Furthermore, we successfully identified unique full-length genes involved in triterpenoid saponin synthesis and plant hormonal signaling pathways, including auxin and cytokinin. Studies on the functional genomics of P. ginseng seedlings have confirmed the rapid upregulation of negative feed-back loops by auxin and cytokinin signaling cues. The conserved evolutionary mechanisms in the auxin and cytokinin canonical signaling pathways of P. ginseng are more complex than those in Arabidopsis thaliana. Our analysis also revealed a more detailed view of transcriptome-wide alternative isoforms for 88 genes. Finally, transposable elements (TEs) were also identified, suggesting transcriptional activity of TEs in P. ginseng. In conclusion, our results suggest that long-read, full-length or partial-unigene data with high-quality assemblies are invaluable resources as transcriptomic references in P. ginseng and can be used for comparative analyses in closely related medicinal plants.

Molecular differentiation of Russian wild ginseng using mitochondrial nad7 intron 3 region

Background

Cultivated ginseng is often introduced as a substitute and adulterant of Russian wild ginseng due to its lower cost or misidentification caused by similarity in appearance with wild ginseng. The aim of this study is to develop a simple and reliable method to differentiate Russian wild ginseng from cultivated ginseng.

Methods

The mitochondrial NADH dehydrogenase subunit 7 (nad7) intron 3 regions of Russian wild ginseng and Chinese cultivated ginseng were analyzed. Based on the multiple sequence alignment result, a specific primer for Russian wild ginseng was designed by introducing additional mismatch and allele-specific polymerase chain reaction (PCR) was performed for identification of wild ginseng. Real-time allele-specific PCR with endpoint analysis was used for validation of the developed Russian wild ginseng single nucleotide polymorphism (SNP) marker.

Results

An SNP site specific to Russian wild ginseng was exploited by multiple alignments of mitochondrial nad7 intron 3 regions of different ginseng samples. With the SNP-based specific primer, Russian wild ginseng was successfully discriminated from Chinese and Korean cultivated ginseng samples by allele-specific PCR. The reliability and specificity of the SNP marker was validated by checking 20 individuals of Russian wild ginseng samples with real-time allele-specific PCR assay.

Conclusion

An effective DNA method for molecular discrimination of Russian wild ginseng from Chinese and Korean cultivated ginseng was developed. The established real-time allele-specific PCR was simple and reliable, and the present method should be a crucial complement of chemical analysis for authentication of Russian wild ginseng.

Comparative transcriptome analysis of shoot and root tissue of Bacopa monnieri identifies potential genes related to triterpenoid saponin biosynthesis

Background

Bacopa monnieri commonly known as Brahmi is utilized in Ayurveda to improve memory and many other human health benefits. Bacosides enriched standardized extract of Bacopa monnieri is being marketed as a memory enhancing agent. In spite of its well known pharmacological properties it is not much studied in terms of transcripts involved in biosynthetic pathway and its regulation that controls the secondary metabolic pathway in this plant. The aim of this study was to identify the potential transcripts and provide a framework of identified transcripts involved in bacosides production through transcriptome assembly.

Results

We performed comparative transcriptome analysis of shoot and root tissue of Bacopa monnieri in two independent biological replicate and obtained 22.48 million and 22.0 million high quality processed reads in shoot and root respectively. After de novo assembly and quantitative assessment total 26,412 genes got annotated in root and 18,500 genes annotated in shoot sample. Quality of raw reads was determined by using SeqQC-V2.2. Assembled sequences were annotated using BLASTX against public database such as NR or UniProt. Searching against the KEGG pathway database indicated that 37,918 unigenes from root and 35,130 unigenes from shoot were mapped to 133 KEGG pathways. Based on the DGE data we found that most of the transcript related to CYP450s and UDP-glucosyltransferases were specifically upregulated in shoot tissue as compared to root tissue. Finally, we have selected 43 transcripts related to secondary metabolism including transcription factor families which are differentially expressed in shoot and root tissues were validated by qRT-PCR and their expression level were monitored after MeJA treatment and wounding for 1, 3 and 5 h.

Conclusions

This study not only represents the first de novo transcriptome analysis of Bacopa monnieri but also provides information about the identification, expression and differential tissues specific distribution of transcripts related to triterpenoid sapogenin which is one of the most important pharmacologically active secondary metabolite present in Bacopa monnieri. The identified transcripts in this study will establish a foundation for future studies related to carrying out the metabolic engineering for increasing the bacosides biosynthesis and its regulation for human health benefits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3865-5) contains supplementary material, which is available to authorized users.

Comprehensive Characterization for Ginsenosides Biosynthesis in Ginseng Root by Integration Analysis of Chemical and Transcriptome

Herbgenomics provides a global platform to explore the genetics and biology of herbs on the genome level. Panax ginseng C.A. Meyer is an important medicinal plant with numerous pharmaceutical effects. Previous reports mainly discussed the transcriptome of ginseng at the organ level. However, based on mass spectrometry imaging analyses, the ginsenosides varied among different tissues. In this work, ginseng root was separated into three tissues-periderm, cortex and stele-each for five duplicates. The chemical analysis and transcriptome analysis were conducted simultaneously. Gene-encoding enzymes involved in ginsenosides biosynthesis and modification were studied based on gene and molecule data. Eight widely-used ginsenosides were distributed unevenly in ginseng roots. A total of 182,881 unigenes were assembled with an N50 contig size of 1374 bp. About 21,000 of these unigenes were positively correlated with the content of ginsenosides. Additionally, we identified 192 transcripts encoding enzymes involved in two triterpenoid biosynthesis pathways and 290 transcripts encoding UDP-glycosyltransferases (UGTs). Of these UGTs, 195 UGTs (67.2%) were more highly expressed in the periderm, and that seven UGTs and one UGT were specifically expressed in the periderm and stele, respectively. This genetic resource will help to improve the interpretation on complex mechanisms of ginsenosides biosynthesis, accumulation, and transportation.

De novo Sequencing and Transcriptome Analysis Reveal Key Genes Regulating Steroid Metabolism in Leaves, Roots, Adventitious Roots and Calli of Periploca sepium Bunge

Periploca sepium Bunge is a traditional medicinal plant, whose root bark is important for Chinese herbal medicine. Its major bioactive compounds are C21 steroids and periplocin, a kind of cardiac glycoside, which are derived from the steroid synthesis pathway. However, research on P. sepium genome or transcriptomes and their related genes has been lacking for a long time. In this study we estimated this species nuclear genome size at 170 Mb (using flow cytometry). Then, RNA sequencing of four different tissue samples of P. sepium (leaves, roots, adventitious roots, and calli) was done using the sequencing platform Illumina/Solexa Hiseq 2,500. After de novo assembly and quantitative assessment, 90,375 all-transcripts and 71,629 all-unigenes were finally generated. Annotation efforts that used a number of public databases resulted in detailed annotation information for the transcripts. In addition, differentially expressed genes (DEGs) were identified by using digital gene profiling based on the reads per kilobase of transcript per million reads mapped (RPKM) values. Compared with the leaf samples (L), up-regulated genes and down-regulated genes were eventually obtained. To deepen our understanding of these DEGs, we performed two enrichment analyses: gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). Here, the analysis focused upon the expression characteristics of those genes involved in the terpene metabolic pathway and the steroid biosynthesis pathway, to better elucidate the molecular mechanism of bioactive steroid synthesis in P. sepium. The bioinformatics analysis enabled us to find many genes that are involved in bioactive steroid biosynthesis. These genes encoded acetyl-CoA acetyltransferase (ACAT), HMG-CoA synthase (HMGS), HMG-CoA reductase (HMGR), mevalonate kinase (MK), phosphomevalonate kinase (PMK), mevalonate diphosphate decarboxylase (MDD), isopentenylpyrophosphate isomerase (IPPI), farnesyl pyrophosphate synthase (FPS), squalene synthase (SS), squalene epoxidase (SE), cycloartenol synthase (CAS), sterol C-24 methyltransferase (SMT1), sterol-4alpha-methyl oxidase 1 (SMO1), sterol 14alpha-demethylase (CYP51/14-SDM), delta(14)-sterol reductase (FK/14SR), C-8,7 sterol isomerase (HYD1), sterol-4alpha-methyl oxidase 2 (SMO2), delta(7)-sterol-C5(6)-desaturase (STE1/SC5DL), 7-dehydrocholesterol reductase (DWF5/DHCR7), delta (24)-sterol reductase (DWF1/DHCR24), sterol 22-desaturase (CYP710A), progesterone 5beta-reductase (5β-POR), 3-beta-hydroxysteroid dehydrogenase (3β-HSD). This research will be helpful to further understand the mechanism of bioactive steroid biosynthesis in P. sepium, namely C21 steroid and periplocin biosynthesis.

Functional regulation of ginsenoside biosynthesis by RNA interferences of a UDP-glycosyltransferase gene in Panax ginseng and Panax quinquefolius

Panax ginseng (Asian ginseng) and Panax quinquefolius (American ginseng) have been used as medicinal and functional herbal remedies worldwide. Different properties of P. ginseng and P. quinquefolius were confirmed not only in clinical findings, but also at cellular and molecular levels. The major pharmacological ingredients of P. ginseng and P. quinquefolius are the triterpene saponins known as ginsenosides. The P. ginseng roots contain a higher ratio of ginsenoside Rg1:Rb1 than that in P. quinquefolius. In ginseng plants, various ginsenosides are synthesized via three key reactions: cyclization, hydroxylation and glycosylation. To date, several genes including dammarenediol synthase (DS), protopanaxadiol synthase and protopanaxatriol synthase have been isolated in P. ginseng and P. quinquefolius. Although some glycosyltransferase genes have been isolated and identified association with ginsenoside synthesis in P. ginseng, little is known about the glycosylation mechanism in P. quinquefolius. In this paper, we cloned and identified a UDP-glycosyltransferase gene named Pq3-O-UGT2 from P. quinquefolius (GenBank accession No. KR106207). In vitro enzymatic activity experiments biochemically confirmed that Pq3-O-UGT2 catalyzed the glycosylation of Rh2 and F2 to produce Rg3 and Rd, and the chemical structure of the products were confirmed susing high performance liquid chromatography electrospray ionization mass spectrometry (HPLC/ESI-MS). High sequence similarity between Pq3-O-UGT2 and PgUGT94Q2 indicated a close evolutionary relationship between P. ginseng and P. quinquefolius. Moreover, we established both P. ginseng and P. quinquefolius RNAi transgenic roots lines. RNA interference of Pq3-O-UGT2 and PgUGT94Q2 led to reduce levels of ginsenoside Rd, protopanaxadiol-type and total ginsenosides. Expression of key genes including protopanaxadiol and protopanaxatriol synthases was up-regulated in RNAi lines, while expression of dammarenediol synthase gene was not obviously increased. These results revealed that P. quinquefolius was more sensitive to the RNAi of Pq3-O-UGT2 and PgUGT94Q2 when compared with P. ginseng.