Analysis of the transcriptome of Panax notoginseng root uncovers putative triterpene saponin-biosynthetic genes and genetic markers

Panax notoginseng: panoramagram of phytochemical and pharmacological properties, biosynthesis, and regulation and production of ginsenosides

Panax notoginseng is a famous perennial herb widely used as material for medicine and health-care food. Due to its various therapeutic effects, research work on P. notoginseng has rapidly increased in recent years, urging a comprehensive review of research progress on this important medicinal plant. Here, we summarize the latest studies on the representative bioactive constituents of P. notoginseng and their multiple pharmacological effects, like cardiovascular protection, anti-tumor, and immunomodulatory activities. More importantly, we emphasize the biosynthesis and regulation of ginsenosides, which are the main bioactive ingredients of P. notoginseng. Key enzymes and transcription factors (TFs) involved in the biosynthesis of ginsenosides are reviewed, including diverse CYP450s, UGTs, bHLH, and ERF TFs. We also construct a transcriptional regulatory network based on multi-omics data and predicted candidate TFs mediating the biosynthesis of ginsenosides. Finally, the current three major biotechnological approaches for ginsenoside production are highlighted. This review covers advances in the past decades, providing insights into quality evaluation and perspectives for the rational utilization and development of P. notoginseng resources. Modern omics technologies facilitate the exploration of the molecular mechanisms of ginsenoside biosynthesis, which is crucial to the breeding of novel P. notoginseng varieties. The identification of functional enzymes for biosynthesizing ginsenosides will lead to the formulation of potential strategies for the efficient and large-scale production of specific ginsenosides.

Mining genomic regions associated with agronomic and biochemical traits in quinoa through GWAS

Quinoa (Chenopodium quinoa Willd.), an Andean crop, is a facultative halophyte food crop recognized globally for its high nutritional value and plasticity to adapt to harsh conditions. We conducted a genome-wide association study on a diverse set of quinoa germplasm accessions. These accessions were evaluated for the following agronomic and biochemical traits: days to 50% flowering (DTF), plant height (PH), panicle length (PL), stem diameter (SD), seed yield (SY), grain diameter (GD), and thousand-grain weight (TGW). These accessions underwent genotyping-by-sequencing using the DNBSeq-G400R platform. Among all evaluated traits, TGW represented maximum broad-sense heritability. Our study revealed average SNP density of ≈ 3.11 SNPs/10 kb for the whole genome, with the lowest and highest on chromosomes Cq1B and Cq9A, respectively. Principal component analysis clustered the quinoa population in three main clusters, one clearly representing lowland Chilean accessions, whereas the other two groups corresponded to germplasm from the highlands of Peru and Bolivia. In our germplasm set, we estimated linkage disequilibrium decay to be ≈ 118.5 kb. Marker-trait analyses revealed major and consistent effect associations for DTF on chromosomes 3A, 4B, 5B, 6A, 7A, 7B and 8B, with phenotypic variance explained (PVE) as high as 19.15%. Nine associations across eight chromosomes were also found for saponin content with 20% PVE by qSPN5A.1. More QTLs were identified for PL and TGW on multiple chromosomal locations. We identified putative candidate genes in the genomic regions associated with DTF and saponin content. The consistent and major-effect genomic associations can be used in fast-tracking quinoa breeding for wider adaptation across marginal environments.

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.

Genome-wide identification, expression analysis, and transcriptome analysis of the IAA gene family in Zoysia japonica

BACKGROUND

AUX/IAA is an essential signaling molecule and has great physiological importance in various plants, but its function in Zoysia japonica remains unknown.

METHODS AND RESULTS

Genome-wide identification and analysis of AUX/IAA genes used bioinformatics methods to investigate the ZjIAA genes' expression of exogenous IAA hydroponics treatment for 2 h by qRT-PCR, control and exogenous IAA treated zoysia were subjected to transcriptome sequencing. ZjIAAs were distributed across the 13 subfamilies by phylogenetic analysis with Oryza sativa and Arabidopsis thaliana. Multiple sequence alignment revealed that the majority of genes were non-canonical ZjIAAs with incomplete domain. The optimal growth concentration of the IAA hormone was 0.05 mM, and the qRT-PCR analysis revealed that eight ZjIAAs were differentially expressed, with seven genes considerably upregulating and one gene significantly downregulating. The result of transcriptome sequencing revealed that 515 differentially expressed genes (DEGs) were identified, with 344 upregulated genes and 171 downregulated genes. A total of 18 genes were annotated as involved in the plant hormone signal transduction pathway. And 8 ZjIAAs exhibited distinct expressions, 7 upregulated, and only one downregulated, according to the qRT-PCR study.

CONCLUSIONS

Genome-wide identification and analysis increased the understanding of the evolution and function of the IAA family in zoysia. DEGs of control and treatment with 0.05 mM exogenous IAA hormone were investigated by transcriptome sequencing. ZjIAAs had substantial variations in the expression of associated genes, with the majority of genes upregulated and 18 genes implicated in plant hormone signal transduction.

Computational studies reveal co-occurrence of two mutations in IL7R gene of high-grade serous carcinoma patients

Major cause of mortality in ovarian cancer can be attributed to a lack of specific and sensitive biomarkers for diagnosis and prognosis of the disease. Uncovering the mutations in genes involved in crucial oncogenic pathways is a key step in discovery and development of novel biomarkers. Whole exome sequencing (WES) is a powerful method for the detection of cancer driver mutations. The present work focuses on identifying functionally damaging mutations in patients with high-grade serous ovarian carcinoma (HGSC) through computational analysis of WES. In this study, WES data of HGSC patients was retrieved from the genomic literature available in sequence read archive, the variants were identified and comprehensive structural and functional analysis was performed. Interestingly, I66T and V138I mutations were found to be co-occurring in the IL7R gene in four out of five HGSC patient samples investigated in this study. The V138I mutation was located in the fibronectin type-3 domain and computationally assessed to be causing disruptive effects on the structure and dynamics of IL7R protein. This mutation was found to be co-occurring with the neutral I66T mutation in the same domain which compensated the disruptive effects of V138I variant. These comprehensive studies point to a hitherto unexplored significant role of the IL7R gene in ovarian carcinoma. It is envisaged that the work will lay the foundation for the development of a novel biomarker with potential application in molecular profiling and in estimation of the disease prognosis.Communicated by Ramaswamy H. Sarma.

Upcoming progress of transcriptomics studies on plants: An overview

Transcriptome sequencing or RNA-Sequencing is a high-resolution, sensitive and high-throughput next-generation sequencing (NGS) approach used to study non-model plants and other organisms. In other words, it is an assembly of RNA transcripts from individual or whole samples of functional and developmental stages. RNA-Seq is a significant technique for identifying gene predictions and mining functional analysis that improves gene ontology understanding mechanisms of biological processes, molecular functions, and cellular components, but there is limited information available on this topic. Transcriptomics research on different types of plants can assist researchers to understand functional genes in better ways and regulatory processes to improve breeding selection and cultivation practices. In recent years, several advancements in RNA-Seq technology have been made for the characterization of the transcriptomes of distinct cell types in biological tissues in an efficient manner. RNA-Seq technologies are briefly introduced and examined in terms of their scientific applications. In a nutshell, it introduces all transcriptome sequencing and analysis techniques, as well as their applications in plant biology research. This review will focus on numerous existing and forthcoming strategies for improving transcriptome sequencing technologies for functional gene mining in various plants using RNA- Seq technology, based on the principles, development, and applications.

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

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

Transcriptomics and Metabolomics Changes Triggered by Inflorescence Removal in Panax notoginseng (Burk.)

The root of Panax notoginseng (Burk.), in which saponins are the major active components, is a famous traditional Chinese medicine used to stop bleeding and to decrease inflammation and heart disease. Inflorescence removal increases the yield and quality of P. notoginseng, but the underlying molecular mechanisms are unknown. Here, the differences between inflorescence-removal treatment and control groups of P. notoginseng were compared using transcriptomics and metabolomics analyses. Illumina sequencing of cDNA libraries prepared from the rhizomes, leaves and roots of the two groups independently identified 6,464, 4,584, and 7,220 differentially expressed genes (DEG), respectively. In total, 345 differentially expressed transcription factors (TFs), including MYB and WRKY family members, were induced by the inflorescence-removal treatment. Additionally, 215 DEGs involved in saponin terpenoid backbone biosynthetic pathways were identified. Most genes involved in the mevalonic acid (MVA) and methylerythritol phosphate (MEP) pathways were activated by inflorescence removal. The co-expression analysis showed that the low expression levels of flavonoid biosynthesis-related genes (e.g., C4H and F3H) decreased the biosynthesis and accumulation of some flavonoids after inflorescence removal. The results not only provide new insights into the fundamental mechanisms underlying the poorly studied inflorescence-removal process in P. notoginseng and other rhizome crops, but they also represent an important resource for future research on gene functions during inflorescence-removal treatments and the reproductive stage.

De novo whole-genome assembly and discovery of genes involved in triterpenoid saponin biosynthesis of Vietnamese ginseng (Panax vietnamensis Ha et Grushv.)

Vietnamese ginseng (Panax vietnamensis Ha et Grushv.), also known as Ngoc Linh ginseng, is a high-value herb in Vietnam. Vietnamese ginseng has been proven to be effective in enhancing the immune system, human memory, anti-stress, anti-inflammatory, anti-cancer, and prevent aging. The present study reports the first draft whole-genome of Vietnamese ginseng and the identification of potential genes involved in the triterpenoid metabolic pathway. De novo whole-genome assembly was performed successfully from a data of approximately 139 Gbps of 394,802,120 high quality reads to generate 9815 scaffolds with an N50 value of 572,722 bp from the leaf of Vietnamese ginseng. The assembled genome of Vietnamese ginseng is 3,001,967,204 bp long containing 79,374 gene models. Among them, there are 55,012 genes (69.30%) were annotated by various public molecular biology databases. The potential genes involved in triterpenoid saponin biosynthesis in Vietnamese ginseng and their metabolic pathway were also predicted." Three genes encoding squalene monooxygenase isozymes in Vietnamese ginseng were cloned, sequenced and characterized. Moreover, expression levels of several key genes involved in terpenoid biosynthesis in different parts of Vietnamese ginseng were also analyzed. The SSR markers were detected by various programs from both of assembly full dataset of Vietnamese ginseng genome and predicted genes. The present work provided important data of the draft whole-genome of Vietnamese ginseng for further studies to understand the role of genes involved in ginsenoside biosynthesis and their metabolic pathway at the molecular level of this rare medicinal species.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01076-1.

Gibberellin Signaling Promotes the Secondary Growth of Storage Roots in Panax ginseng

Gibberellins (GAs) are an important group of phytohormones associated with diverse growth and developmental processes, including cell elongation, seed germination, and secondary growth. Recent genomic and genetic analyses have advanced our knowledge of GA signaling pathways and related genes in model plant species. However, functional genomics analyses of GA signaling pathways in Panax ginseng, a perennial herb, have rarely been carried out, despite its well-known economical and medicinal importance. Here, we conducted functional characterization of GA receptors and investigated their physiological roles in the secondary growth of P. ginseng storage roots. We found that the physiological and genetic functions of P. ginseng gibberellin-insensitive dwarf1s (PgGID1s) have been evolutionarily conserved. Additionally, the essential domains and residues in the primary protein structure for interaction with active GAs and DELLA proteins are well-conserved. Overexpression of PgGID1s in Arabidopsis completely restored the GA deficient phenotype of the Arabidopsis gid1a gid1c (atgid1a/c) double mutant. Exogenous GA treatment greatly enhanced the secondary growth of tap roots; however, paclobutrazol (PCZ), a GA biosynthetic inhibitor, reduced root growth in P. ginseng. Transcriptome profiling of P. ginseng roots revealed that GA-induced root secondary growth is closely associated with cell wall biogenesis, the cell cycle, the jasmonic acid (JA) response, and nitrate assimilation, suggesting that a transcriptional network regulate root secondary growth in P. ginseng. These results provide novel insights into the mechanism controlling secondary root growth in P. ginseng.

Research Progress and Future Development Trends in Medicinal Plant Transcriptomics

Transcriptomics is one of the most popular topics in biology in recent times. Transcriptome sequencing (RNA-Seq) is a high-throughput, high-sensitivity, and high-resolution technique that can be used to study model and non-model organisms. Transcriptome sequencing is also an important method for studying the genomes of medicinal plants, a topic on which limited information is available. The study of medicinal plants through transcriptomics can help researchers analyze functional genes and regulatory mechanisms of medicinal plants and improve breeding selection and cultivation techniques. This article analyzes and compares the applications of transcriptome sequencing in medicinal plants over the past decade and briefly introduces the methods of transcriptome sequencing and analysis, their applications in medicinal plant research, and potential development trends. We will focus on the research and application progress of transcriptome sequencing in the following four areas: the mining of functional genes in medicinal plants, development of molecular markers, biosynthetic pathways of secondary metabolites, and developmental mechanisms of medicinal plants. Our review will provide ideas for the mining of functional genes of medicinal plants and breeding new varieties.

Phytochemistry of ginsenosides: Recent advancements and emerging roles

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

Transcriptional Factor-Mediated Regulation of Active Component Biosynthesis in Medicinal Plants

Plants produce thousands of chemically diverse secondary metabolites, many of which have valuable pharmaceutical properties. There is much interest in the synthesis of these pharmaceuticallyvaluable compounds, including the key enzymes and the transcription factors involved. The function and regulatory mechanism of transcription factors in biotic and abiotic stresses have been studied in depth. However, their regulatory roles in the biosynthesis of bioactive compounds, especially in medicinal plants, have only begun. Here, we review what is currently known about how transcription factors contribute to the synthesis of bioactive compounds (alkaloids, terpenoids, flavonoids, and phenolic acids) in medicinal plants. Recent progress has been made in the cloning and characterization of transcription factors in medicinal plants on the genome scale. So far, several large transcription factors have been identified in MYB, WRKY, bHLH, ZIP, AP2/ERF transcription factors. These transcription factors have been predicted to regulate bioactive compound production. These transcription factors positively or negatively regulate the expression of multiple genes encoding key enzymes, and thereby control the metabolic flow through the biosynthetic pathway. Although the research addressing this niche topic is in its infancy, significant progress has been made, and advances in high-throughput sequencing technology are expected to accelerate the discovery of key regulatory transcription factors in medicinal plants. This review is likely to be useful for those interested in the synthesis of pharmaceutically- valuable plant compounds, especially those aiming to breed or engineer plants that produce greater yields of these compounds.

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.

Full-Length Transcriptome Analyses of Genes Involved in Triterpenoid Saponin Biosynthesis of Psammosilene tunicoides Hairy Root Cultures With Exogenous Salicylic Acid

Triterpenoid saponins constitute a diverse class of bioactive compounds in medicinal plants. Salicylic acid (SA) is an efficient elicitor for secondary metabolite production, but a transcriptome-wide regulatory network of SA-promoted triterpenoid saponin biosynthesis remains little understood. In the current study, we described the establishment of the hairy root culture system for Psammosilene tunicoides, a triterpenoid saponin-producing medicinal herb in China, using genetic transformation by Agrobacterium rhizogenes. Compared to controls, we found that total saponin content was dramatically increased (up to 2.49-fold) by the addition of 5 mg/L SA in hairy roots for 1 day. A combination of single-molecule real-time (SMRT) and next-generation sequencing (Illumina RNA-seq) was generated to analyze the full-length transcriptome data for P. tunicoides, as well as the transcript profiles in treated (8 and 24 h) and non-treated (0 h) groups with 5 mg/L SA in hairy roots. A total of 430,117 circular consensus sequence (CCS) reads, 16,375 unigenes and 4,678 long non-coding RNAs (lncRNAs) were obtained. The average length of unigenes (2,776 bp) was much higher in full-length transcriptome than that derived from single RNA-seq (1,457 bp). The differentially expressed genes (DEGs) were mainly enriched in the metabolic process. SA up-regulated the unigenes encoding SA-binding proteins and antioxidant enzymes in comparison with controls. Additionally, we identified 89 full-length transcripts encoding enzymes putatively involved in saponin biosynthesis. The candidate transcription factors (WRKY, NAC) and structural genes (AACT, DXS, SE, CYP72A) might be the key regulators in SA-elicited saponin accumulation. Their expression was further validated by quantitative real-time PCR (qRT-PCR). These findings preliminarily elucidate the regulatory mechanisms of SA on triterpenoid saponin biosynthesis in the transcriptomic level, laying a foundation for SA-elicited saponin augmentation in P. tunicoides.

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.

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.

De novo assembly and Transcriptome characterization of an endemic species of Vietnam, Panax vietnamensis Ha et Grushv., including the development of EST-SSR markers for population genetics

BACKGROUND

Understanding the genetic diversity in endangered species that occur inforest remnants is necessary to establish efficient strategies for the species conservation, restoration and management. Panax vietnamensis Ha et Grushv. is medicinally important, endemic and endangered species of Vietnam. However, genetic diversity and structure of population are unknown due to lack of efficient molecular markers.

RESULTS

In this study, we employed Illumina HiSeq™ 4000 sequencing to analyze the transcriptomes of P. vietnamensis (roots, leaves and stems). Raw reads total of 23,741,783 was obtained and then assembled, from which the generated unigenes were 89,271 (average length = 598.3191 nt). The 31,686 unigenes were annotated in different databases i.e. Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, Nucleotide Collection (NR/NT) and Swiss-Prot for functional annotation. Further, 11,343 EST-SSRs were detected. From 7774 primer pairs, 101 were selected for polymorphism validation, in which; 20 primer pairs were successfully amplified to DNA fragments and significant amounts of polymorphism was observed within population. The nine polymorphic microsatellite loci were used for population structure and diversity analyses. The obtained results revealed high levels of genetic diversity in populations, the average observed and expected heterozygosity were HO = 0.422 and HE = 0.479, respectively. During the Bottleneck analysis using TPM and SMM models (p < 0.01) shows that targeted population is significantly heterozygote deficient. This suggests sign of the bottleneck in all populations. Genetic differentiation between populations was moderate (FST = 0.133) and indicating slightly high level of gene flow (Nm = 1.63). Analysis of molecular variance (AMOVA) showed 63.17% of variation within individuals and 12.45% among populations. Our results shows two genetic clusters related to geographical distances.

CONCLUSION

Our study will assist conservators in future conservation management, breeding, production and habitats restoration of the species.

Comparative transcriptome analysis of MeJA-responsive AP2/ERF transcription factors involved in notoginsenosides biosynthesis

Differential transcriptome analysis is an effective method for gene selection of triterpene saponin biosynthetic pathways. MeJA-induced differential transcriptome of Panax notoginseng has not been analyzed yet. In this study, comparative transcriptome analysis of P. notoginseng roots and methyl jasmonate (MeJA)-induced roots revealed 83,532 assembled unigenes and 21,947 differentially expressed unigenes. Sixteen AP2/ERF transcription factors, which were significantly induced by MeJA treatment in the root of P. notoginseng, were selected for further analysis. Real-time quantitative PCR (RT-qPCR) and co-expression network analysis of the 16 AP2/ERF transcription factors showed that PnERF2 and PnERF3 had significant correlation with dammarenediol II synthase gene (DS) and squalene epoxidase gene (SE), which are key genes in notoginsenoside biosynthesis, in different tissues and MeJA-induced roots. A phylogenetic tree was conducted to analyze the 16 candidate AP2/ERF transcription factors and other 38 transcription factors. The phylogenetic tree analysis showed PnERF2, AtERF3, AtERF7, TcERF12 and other seven transcriptional factors are in same branch, while PnERF3 had close evolutionary relationships with AtDREB1A, GhERF38 and TcAP2. The results of comparative transcriptomes and AP2/ERF transcriptional factors analysis laid a solid foundation for further investigations of disease resistance and notoginsenoside biosynthesis in P. notoginseng.

Transcriptomic response to cold of thermophilous medicinal plant Marsdenia tenacissima

Extracts from Marsdenia tenacissima, involving tenacissoside H, I and G, have been used as remedies of cancer, inflammation and asthma. Low temperature serves as one of the main factors constrain the planting expansion and quality of M. tenacissima, but its functional mechanism has been known scarcely for the lack of genomic information and transcriptional profile. Here we investigated the transcriptomic responses of M. tenacissima under cold stress to gain insight into the molecular mechanism of low temperature sensitivity. Total RNAs were collected from samples obtained at 4-time points (after 0, 3, 6 and 48 h cold treatments with 4 °C, respectively), then used for library construction and sequenced on the Illumina Hiseq™ 4000 platform. Passing quality assessments, 500794 transcripts, and 206137 unigenes were de novo assembly out in Trinity v2.4.0, holding contig N50 of 2566 bp and unigene mean length of 754 bp. 44.20% of assembled unigenes were annotated to the well-known public protein database on a basis of sequence similarity. Using statistical comparison of the fragments per kilo base of transcript per million reads mapped (FPKM) values between conditions, 6082 group-specific differentially expressed genes (DEGs) were identified and considered as cold-responsive genes, which contained copious transcription factors and active secondary metabolism. Among them, 43 unigenes were constantly up-regulated expression along with cold time, which mainly implicated in the biosynthesis of secondary metabolites, carbon metabolism, RNA and DNA metabolism. Conversely, 21 unigenes involved in photosynthesis, cell wall, protein degradation, and transporters were downregulated continually with cold timescale. Experimentally, MtEF1α was chosen as the best housekeeping gene. Functional enrichments found that damaging of cold stress on M. tenacissima may be ascribed to inability of photosynthesis, ribsome processing, flavonoid biosynthesis and terpenoids degradation. Correlation analysis between cold induced transcription factors and tenacissoside biosynthesis-related genes indicated that 3β-HSD significant positively correlated with bHLH51, and 4-MSO with NF-YB, GRAS3, Trihelix, FAR1, MYB60, MYBS1, bZIP43. Further promoter clone found MYB-binding site in the promoter of 4-MSO. In view of the reported cold tolerance of MYB60, it is recommended as a potential candidate suitable for future molecular design of exaptation cultivation with high bioactive constituents.

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.

Recent Advances in the Metabolic Engineering of Yeasts for Ginsenoside Biosynthesis

Ginsenosides are a group of glycosylated triterpenes isolated from Panax species. Ginsenosides are promising candidates for the prevention and treatment of cancer as well as food additives. However, owing to a lack of efficient approaches for ginsenoside production from plants and chemical synthesis, ginsenosides may not yet have reached their full potential as medicinal resources. In recent years, an alternative approach for ginsenoside production has been developed using the model yeast Saccharomyces cerevisiae and non-conventional yeasts such as Yarrowia lipolytica and Pichia pastoris. In this review, various metabolic engineering strategies, including heterologous gene expression, balancing, and increasing metabolic flux, and enzyme engineering, have been described as recent advanced engineering techniques for improving ginsenoside production. Furthermore, the usefulness of a systems approach and fermentation strategy has been presented. Finally, the present challenges and future research direction for industrial cell factories have been discussed.

Mining genes associated with furanocoumarin biosynthesis in an endangered medicinal plant, Glehnia littoralis

The endangered medicinal plant Glehnia littoralis is one of the important natural source of furanocoumarin, which has been used as mucolytic, antitussive, antitumour and antibacterial. However, the genetic information of furanocoumarin biosynthesis in G. littoralis is scarce at present. The objective of this study was to mine the putative candidate genes involved in the biosynthesis pathwayof furanocoumarin and provide references for gene identification, and functional genomics of G. littoralis. We carried out the transcriptome analysis of leaves and roots in G. littoralis, which provided a dataset for gene mining. Psoralen, imperatorin and isoimperatorin were detected in G. littoralis by high performance liquid chromatography analysis. Candidate key genes were mined based on the annotations and local BLAST with homologous sequences using BioEdit software. The relative expression of genes was analysed using quantitative real-time polymerase chain reaction. Further, the CYP450 genes were mined using phylogenetic analyses using MEGA 6.0 software. Atotal of 156,949 unigenes were generated, of which 9021 were differentially-expressed between leaves and roots. A total of 82 unigenes encoding eight enzymes in furanocoumarin biosynthetic pathway were first obtained. Seven genes that encoded key enzymes in the downstream furanocoumarin biosynthetic pathway and expressed more in roots than leaves were screened. Twenty-six candidate CYP450 unigenes expressed abundantly in roots and were chiefly concentrated in CYP71, CYP85 and CYP72 clans. Finally, we filtered 102 differentially expressed transcription factors (TFs) unigenes. The transcriptome of G. littoralis was characterized which would help to elucidate the furanocoumarin biosynthetic pathway in G. littoralis and provide an invaluable resource for further study of furanocoumarin.

Structure and Location Studies on Key Enzymes in Saponins Biosynthesis of Panax notoginseng

Panax notoginseng is one of the most widely used traditional herbs for the treatment of various diseases, in which saponins were the main active components. At present, the research of P. notoginseng mainly focused on the discovery of new compounds and pharmacology. However, there were few studies on the molecular mechanism of the synthesis of secondary metabolites of P. notoginseng. In our study, four coding sequences (CDS) encoding the key enzymes involved in saponin biosynthesis were cloned, namely farnesyl diphosphate synthase (FPS), squalene synthase (SS), squalene epoxidase (SE), and dammarenediol-II synthase (DS), which contained open reading frame (ORF) of 1029 bp, 1248 bp, 1614 bp, and 2310 bp, and coded 342, 415, 537, and 769 amino acids, respectively. At the same time, their domains, secondary structures, three-dimensional structures, and phylogenetics trees were analyzed by kinds of bioinformatics tools. Their phylogenetics relationships were also analyzed. In addition, GFP (Green fluorescent protein) fusion genes were constructed by the plasmid transformation system to determine the subcellular localization. The results of subcellular localization showed that FPS, SE, and DS were mainly located in cytomembrane and its surrounding, while SS was located both in cytoplasm and cytomembrane. Our findings provided data demonstrating the expression patterns of genes involved in saponin biosynthesis and would facilitate efforts to further elucidate the biosynthesis of the bioactive components in P. notoginseng.

Glycosyltransformation of ginsenoside Rh2 into two novel ginsenosides using recombinant glycosyltransferase from Lactobacillus rhamnosus and its in vitro applications

Background

Ginsenoside Rh2 is well known for many pharmacological activities, such as anticancer, antidiabetes, antiinflammatory, and antiobesity properties. Glycosyltransferases (GTs) are ubiquitous enzymes present in nature and are widely used for the synthesis of oligosaccharides, polysaccharides, glycoconjugates, and novel derivatives. We aimed to synthesize new ginsenosides from Rh2 using the recombinant GT enzyme and investigate its cytotoxicity with diverse cell lines.

Methods

We have used a GT gene with 1,224-bp gene sequence cloned from Lactobacillus rhamnosus (LRGT) and then expressed in Escherichia coli BL21 (DE3). The recombinant GT protein was purified and demonstrated to transform Rh2 into two novel ginsenosides, and they were characterized by nuclear magnetic resonance (NMR) techniques and evaluated by 3-(4, 5-dimethylthiazol-2-yl)-2-5-diphenyltetrazolium bromide assay.

Results

Two novel ginsenosides with an additional glucopyranosyl (6→1) and two additional glucopyranosyl (6→1) linked with the C-3 position of the substrate Rh2 were synthesized, respectively. Cell viability assay in the lung cancer (A549) cell line showed that glucosyl ginsenoside Rh2 inhibited cell viability more potently than ginsenoside Rg3 and Rh2 at a concentration of 10 μM. Furthermore, glucosyl ginsenoside Rh2 did not exhibit any cytotoxic effect in murine macrophage cells (RAW264.7), mouse embryo fibroblasts cells (3T3-L1), and skin cells (B16BL6) at a concentration of 10 μM compared with ginsenoside Rh2 and Rg3.

Conclusion

This is the first report on the synthesis of two novel ginsenosides, namely, glucosyl ginsenoside Rh2 and diglucosyl ginsenoside Rh2 from Rh2 by using recombinant GT isolated from L. rhamnosus. Moreover, diglucosyl ginsenoside Rh2 might be a new candidate for treatment of inflammation, obesity, and skin whiting, and especially for anticancer.

A transcriptome analysis uncovers Panax notoginseng resistance to Fusarium solani induced by methyl jasmonate

BACKGROUND

Panax notoginseng is a famous Chinese herbal medicine, but the root rot disease mainly caused by Fusarium solani severely reduces the yield and quality of its medicinal materials.

OBJECTIVE

The defense priming in P. notoginseng through exogenous application of signaling molecule will supply theoretical support for the exogenous regulation of disease resistance in P. notoginseng.

METHODS

In this study, the exogenous application of methyl jasmonate (MeJA) increased P. notoginseng's resistance to F. solani. Furthermore, the P. notoginseng transcriptome during F. solani infection was investigated through next-generation sequencing to uncover the resistance mechanism of P. notogingseng induced by MeJA.

RESULTS

The de novo assembly of transcriptome sequences produced 80,551 unigenes, and 36,771 of these unigenes were annotated by at least one database. A differentially expressed gene analysis revealed that a large number of genes related to terpenoid backbone biosynthesis, phenylalanine metabolism, and plant-pathogen interactions were predominantly up-regulated by MeJA. Moreover, jasmonic acid (JA) biosynthesis-related genes and the JA signaling pathway genes, such as linoleate 13S-lipoxygenase, allene oxide cyclase, allene oxide synthase, TIFY, defensin, and pathogenesis-related proteins, showed increased transcriptional levels after inoculation with F. solani. Notably, according to the gene expression analysis, JA and ethylene signaling pathways may act synergistically to positively regulate the defense responses of P. notoginseng to F. solani.

CONCLUSION

JA signaling appears to play a vital role in P. notoginseng responses to F. solani infection, which will be helpful in improving the disease resistance of P. notoginseng cultivars as well as in developing an environmentally friendly biological control method for root rot disease.

Role of Chinese Herbal Medicines in Regulation of Energy Metabolism in Treating Cardiovascular Diseases

Recently, studying myocardial energy metabolism pathways or improving myocardial metabolism through drugs is another effective strategy for treating ischemic heart disease. Many active components of Chinese herbal medicines (CHMs) have been found to modulate energy metabolism in myocardial cells, cerebral vascular cells, endothelial cells and tumour cells. This paper reviews the advances in studies on the active components of CHMs that modulating energy metabolism in treating cardiovascular diseases over the past five years.

Role of Saponins in Plant Defense Against Specialist Herbivores

The diamondback moth (DBM), Plutella xylostella (Lepidoptera: Plutellidae) is a very destructive crucifer-specialized pest that has resulted in significant crop losses worldwide. DBM is well attracted to glucosinolates (which act as fingerprints and essential for herbivores in host plant recognition) containing crucifers such as wintercress, Barbarea vulgaris (Brassicaceae) despite poor larval survival on it due to high-to-low concentration of saponins and generally to other plants in the genus Barbarea. B. vulgaris build up resistance against DBM and other herbivorous insects using glucosinulates which are used in plant defense. Aside glucosinolates, Barbarea genus also contains triterpenoid saponins, which are toxic to insects and act as feeding deterrents for plant specialist herbivores (such as DBM). Previous studies have found interesting relationship between the host plant and secondary metabolite contents, which indicate that attraction or resistance to specialist herbivore DBM, is due to higher concentrations of glucosinolates and saponins in younger leaves in contrast to the older leaves of Barbarea genus. As a response to this phenomenon, herbivores as DBM has developed a strategy of defense against these plant biochemicals. Because there is a lack of full knowledge in understanding bioactive molecules (such as saponins) role in plant defense against plant herbivores. Thus, in this review, we discuss the role of secondary plant metabolites in plant defense mechanisms against the specialist herbivores. In the future, trials by plant breeders could aim at transferring these bioactive molecules against herbivore to cash crops.

Metabolomes and transcriptomes revealed the saponin distribution in root tissues of Panax quinquefolius and Panax notoginseng

Background

Panax quinquefolius and Panax notoginseng are widely used and well known for their pharmacological effects. As main pharmacological components, saponins have different distribution patterns in the root tissues of Panax plants.

Methods

In this study, the representative ginsenosides were detected and quantified by desorption electrospray ionization mass spectrometry and high-performance liquid chromatography analysis to demonstrate saponin distribution in the root tissues of P. quinquefolius and P. notoginseng, and saponin metabolite profiles were analyzed by metabolomes to obtain the biomarkers of different root tissues. Finally, the transcriptome analysis was performed to demonstrate the molecular mechanisms of saponin distribution by gene profiles.

Results

There was saponin distribution in the root tissues differed between P. quinquefolius and P. notoginseng. Eight-eight and 24 potential biomarkers were detected by metabolome analysis, and a total of 340 and 122 transcripts involved in saponin synthesis that were positively correlated with the saponin contents (R > 0.6, P < 0.05) in the root tissues of P. quinquefolius and P. notoginseng, respectively. Among them, GDPS1, CYP51, CYP64, and UGT11 were significantly correlated with the contents of Rg1, Re, Rc, Rb2, and Rd in P. quinquefolius. UGT255 was markedly related to the content of R1; CYP74, CYP89, CYP100, CYP103, CYP109, and UGT190 were markedly correlated with the Rd content in P. notoginseng.

Conclusions

These results provided the visual and quantitative profiles of and confirmed the pivotal transcripts of CYPs and UGTs regulating the saponin distribution in the root tissues of P. quinquefolius and P. notoginseng.

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.

Oleanane-Type Saponins Biosynthesis in Panax notoginseng via Transformation of β-Amyrin Synthase Gene from Panax japonicus

Oleanane-type saponins considered as the main medicinal ingredients in Panax japonicus are not found in Panax notoginseng. β-Amyrin synthase (βAS) was recognized as the first key enzyme in the biosynthetic branch of oleanane-type saponins. In this study, βAS gene from P. japonicus ( PjβAS) was transferred into P. notoginseng cells. Along with PjβAS expression in the transgenic cells, the expression levels of several key enzyme genes related to triterpenoid saponins biosynthesis and the content of P. notoginseng saponins were also increased. Two oleanane-type saponins, chikusetsusaponin IV and chikusetsusaponin IVa, contained in P. japonicus were first discovered in transgenic P. notoginseng cells. This study successfully constructed a biosynthetic pathway of oleanane-type saponins in P. notoginseng by introducing just one gene into the species. On the basis of this discovery and previous studies, the common biosynthetic pathway of triterpenoid saponins in Panax genus may be unified to some extent.

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

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

Transcriptome analysis of Panax ginseng response to high light stress

Background

Ginseng (Panax ginseng Meyer) is an essential source of pharmaceuticals and functional foods. Ginseng productivity has been compromised by high light (HL) stress, which is one of the major abiotic stresses during the ginseng cultivation period. The genetic improvement for HL tolerance in ginseng could be facilitated by analyzing its genetic and molecular characteristics associated with HL stress.

Methods

Genome-wide analysis of gene expression was performed under HL and recovery conditions in 1-year-old Korean ginseng (P. ginseng cv. Chunpoong) using the Illumina HiSeq platform. After de novo assembly of transcripts, we performed expression profiling and identified differentially expressed genes (DEGs). Furthermore, putative functions of identified DEGs were explored using Gene Ontology terms and Kyoto Encyclopedia of Genes and Genome pathway enrichment analysis.

Results

A total of 438 highly expressed DEGs in response to HL stress were identified and selected from 29,184 representative transcripts. Among the DEGs, 326 and 114 transcripts were upregulated and downregulated, respectively. Based on the functional analysis, most upregulated and a significant number of downregulated transcripts were related to stress responses and cellular metabolic processes, respectively.

Conclusion

Transcriptome profiling could be a strategy to comprehensively elucidate the genetic and molecular mechanisms of HL tolerance and susceptibility. This study would provide a foundation for developing breeding and metabolic engineering strategies to improve the environmental stress tolerance of ginseng.

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.

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.

Integrated metabolomic and transcriptomic analyses revealed the distribution of saponins in Panax notoginseng

Panax notoginseng is famous for its important therapeutic effects. Saponins are bioactive compounds found in different parts and developmental stages of P. notoginseng plants. Thus, it is urgently to study saponins distribution in different parts and growth ages of P. notoginseng plants. In this study, potential biomarkers were found, and their chemical characteristic differences were revealed through metabolomic analysis. High-performance liquid chromatography data indicated the higher content of saponins (i.e., Rg1, Re, Rd, and Rb1) in the underground parts than that in the aerial parts. 20(S)-Protopanaxadiol saponins were mainly distributed in the aerial parts. Additionally, the total saponin content in the 3-year-old P. notoginseng plant (188.0 mg/g) was 1.4-fold higher than that in 2-year-old plant (130.5 mg/g). The transcriptomic analysis indicated the tissue-specific transcription expression of genes, namely, PnFPS, PnSS, PnSE1, PnSE2, and PnDS, which encoded critical synthases in saponin biosyntheses. These genes showed similar expression patterns among the parts of P. notoginseng plants. The expression levels of these genes in the flowers and leaves were 5.2fold higher than that in the roots and fibrils. These results suggested that saponins might be actively synthesized in the aerial parts and transformed to the underground parts. This study provides insights into the chemical and genetic characteristics of P. notoginseng to facilitate the synthesis of its secondary metabolites and a scientific basis for appropriate collection and rational use of this plant.

Progress on the Studies of the Key Enzymes of Ginsenoside Biosynthesis

As the main bioactive constituents of Panax species, ginsenosides possess a wide range of notable medicinal effects such as anti-cancer, anti-oxidative, antiaging, anti-inflammatory, anti-apoptotic and neuroprotective activities. However, the increasing medical demand for ginsenosides cannot be met due to the limited resource of Panax species and the low contents of ginsenosides. In recent years, biotechnological approaches have been utilized to increase the production of ginsenosides by regulating the key enzymes of ginsenoside biosynthesis, while synthetic biology strategies have been adopted to produce ginsenosides by introducing these genes into yeast. This review summarizes the latest research progress on cloning and functional characterization of key genes dedicated to the production of ginsenosides, which not only lays the foundation for their application in plant engineering, but also provides the building blocks for the production of ginsenosides by synthetic biology.

Genome-Wide Variation Patterns Uncover the Origin and Selection in Cultivated Ginseng (Panax ginseng Meyer)

Chinese ginseng (Panax ginseng Meyer) is a medicinally important herb and plays crucial roles in traditional Chinese medicine. Pharmacological analyses identified diverse bioactive components from Chinese ginseng. However, basic biological attributes including domestication and selection of the ginseng plant remain under-investigated. Here, we presented a genome-wide view of the domestication and selection of cultivated ginseng based on the whole genome data. A total of 8,660 protein-coding genes were selected for genome-wide scanning of the 30 wild and cultivated ginseng accessions. In complement, the 45s rDNA, chloroplast and mitochondrial genomes were included to perform phylogenetic and population genetic analyses. The observed spatial genetic structure between northern cultivated ginseng (NCG) and southern cultivated ginseng (SCG) accessions suggested multiple independent origins of cultivated ginseng. Genome-wide scanning further demonstrated that NCG and SCG have undergone distinct selection pressures during the domestication process, with more genes identified in the NCG (97 genes) than in the SCG group (5 genes). Functional analyses revealed that these genes are involved in diverse pathways, including DNA methylation, lignin biosynthesis, and cell differentiation. These findings suggested that the SCG and NCG groups have distinct demographic histories. Candidate genes identified are useful for future molecular breeding of cultivated ginseng.

Next-generation sequencing: recent applications to the analysis of colorectal cancer

Since the establishment of the Sanger sequencing method, scientists around the world focused their efforts to progress in the field to produce the utmost technology. The introduction of next-generation sequencing (NGS) represents a revolutionary step and promises to lead to massive improvements in our understanding on the role of nucleic acids functions. Cancer research began to use this innovative and highly performing method, and interesting results started to appear in colorectal cancer (CRC) analysis. Several studies produced high-quality data in terms of mutation discovery, especially about actionable or less frequently mutated genes, epigenetics, transcriptomics. Analysis of results is unveiling relevant perspectives aiding to evaluate the response to therapies. Novel evidences have been presented also in other directions such as gut microbiota or CRC circulating tumor cells. However, despite its unquestioned potential, NGS poses some issues calling for additional studies. This review intends to offer a view of the state of the art of NGS applications to CRC through examination of the most important technologies and discussion of recent published results.

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.

Construct a gene-to-metabolite network to screen the key genes of triterpene saponin biosynthetic pathway in Panax notoginseng

Triterpene saponins are main active constituents of Panax notoginseng. Metabolites profiling of 12 triterpene saponins was analyzed by high-performance liquid chromatography-mass spectrometry in leaf, petiole, and root extracts of P. notoginseng. Most of the 20(S)-protopanaxatriol (PPT) type saponins, except ginsenoside Re, were mainly distributed in roots, while 20(S)-protopanaxadiol (PPD) type saponins were detected among various tissues. The total content of PPD-type saponins decreased in the order of leaf, petiole, and root. The expression patterns of four key genes (PnFPS, PnSQS, PnDS, and PnSE) in the triterpene saponin biosynthetic pathway were measured by real-time quantitative PCR (RT-qPCR). All the four investigated genes showed high expression levels in leaf. A gene-to-metabolite network was constructed through canonical correlation analysis. The results indicated that the expression levels of PnFPS, PnSQS, PnDS, and PnSE had high correlation with PPD-type saponins ginsenoside Rb₂ , Rb₃ , and Rc, while PnSQS was also highly correlated with Rb₁ . Combining metabolic profiling, RT-qPCR, and gene-to-metabolite network, we inferred that the leaf of P. notoginseng was the main biosynthesis site of PPD-type saponins Rb₂ , Rb₃ , and Rc. The contribution to the biosynthesis of ginsenosides Rb₂ , Rb₃ , and Rc was in the order of PnSE > PnDS > PnSQS > PnFPS. PnSE and PnDS should be the preferred targets to regulate the production of PPD-type saponins Rb₂ , Rb₃ , and Rc in P. notoginseng by plant metabolic engineering.

RNA-sequencing-based transcriptome and biochemical analyses of steroidal saponin pathway in a complete set of Allium fistulosum-A. cepa monosomic addition lines

The genus Allium is a rich source of steroidal saponins, and its medicinal properties have been attributed to these bioactive compounds. The saponin compounds with diverse structures play a pivotal role in Allium’s defense mechanism. Despite numerous studies on the occurrence and chemical structure of steroidal saponins, their biosynthetic pathway in Allium species is poorly understood. The monosomic addition lines (MALs) of the Japanese bunching onion (A. fistulosum, FF) with an extra chromosome from the shallot (A. cepa Aggregatum group, AA) are powerful genetic resources that enable us to understand many physiological traits of Allium. In the present study, we were able to isolate and identify Alliospiroside A saponin compound in A. fistulosum with extra chromosome 2A from shallot (FF2A) and its role in the defense mechanism against Fusarium pathogens. Furthermore, to gain molecular insight into the Allium saponin biosynthesis pathway, high-throughput RNA-Seq of the root, bulb, and leaf of AA, MALs, and FF was carried out using Illumina's HiSeq 2500 platform. An open access Allium Transcript Database (Allium TDB, http://alliumtdb.kazusa.or.jp) was generated based on RNA-Seq data. The resulting assembled transcripts were functionally annotated, revealing 50 unigenes involved in saponin biosynthesis. Differential gene expression (DGE) analyses of AA and MALs as compared with FF (as a control) revealed a strong up-regulation of the saponin downstream pathway, including cytochrome P450, glycosyltransferase, and beta-glucosidase in chromosome 2A. An understanding of the saponin compounds and biosynthesis-related genes would facilitate the development of plants with unique saponin content and, subsequently, improved disease resistance.

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.

De novo assembly and annotation of the Zhe-Maidong (Ophiopogon japonicus (L.f.) Ker-Gawl) transcriptome in different growth stages

Zhe-Maidong (Ophiopogon japonicus (L.f.) Ker-Gawl) is a traditional medicinal herb in the family Liliaceae that has significant pharmacological effects on immunity and cardiovascular disease. In this study, three different growth stages of Zhe-Maidong were investigated using RNA-seq, and a total of 16.4 Gb of raw data was obtained. After filtering and assembling, 96,738 unigenes with an average length of 605.3 bp were ultimately generated. A total of 77,300 unigenes were annotated using information from five databases, including the NT, NR, SwissProt, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) databases. Additionally, the mechanisms of flavonoid, saponin and polysaccharide biosynthesis and of accumulation at different stages of tuber development were also characterized. From the first to third years, the contents of flavonoids, saponins and polysaccharides all increased, whereas the expression levels of related genes decreased in the flavonoid and saponin pathways and first increased and then decreased in the polysaccharide pathway. The results of this study provide the most comprehensive expressed sequence resource for Zhe-Maidong and will expand the available O. japonicus gene library and facilitate further genome-wide research and analyses of this species.