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

Genome-Wide Identification of the CYP716 Gene Family in Platycodon grandiflorus (Jacq.) A. DC. and Its Role in the Regulation of Triterpenoid Saponin Biosynthesis

The main type of saponins occurring in the root of Platycodon grandiflorus (Jacq.) A. DC. are oleanolic acid glycosides. The CYP716 gene family plays a major role in catalyzing the conversion of β-amyrin into oleanolic acid. However, studies on the CYP716 genes in P. grandiflorus are limited, and its evolutionary history remains poorly understood. In this study, 22 PgCYP716 genes were identified, distributed among seven subfamilies. Cis-acting elements of the PgCYP716 promoters were mainly involved in plant hormone regulation and responses to abiotic stresses. PgCYP716A264, PgCYP716A391, PgCYP716A291, and PgCYP716BWv3 genes were upregulated in the root and during saponin accumulation, as shown by RNA-seq analysis, suggesting that these four genes play an important role in saponin synthesis. The results of subcellular localization indicated that these four genes encoded membrane proteins. Furthermore, the catalytic activity of these four genes was proved in the yeast, which catalyzed the conversion of β-amyrin into oleanolic acid. We found that the content of β-amyrin, platycodin D, platycoside E, platycodin D3, and total saponins increased significantly when either of the four genes was over expressed in the transgenic hair root. In addition, the expression of PgSS, PgGPPS2, PgHMGS, and PgSE was also upregulated while these four genes were overexpressed. These data support that these four PgCYP716 enzymes oxidize β-amyrin to produce oleanolic acid, ultimately promoting saponin accumulation by activating the expression of upstream pathway genes. Our results enhanced the understanding of the functional variation among the PgCYP716 gene family involved in triterpenoid biosynthesis and provided a theoretical foundation for improving saponin content and enriching the saponin biosynthetic pathway in P. grandiflorus.

Genome-wide identification of bZIP transcription factors and their expression analysis in Platycodon grandiflorus under abiotic stress

The Basic Leucine Zipper (bZIP) transcription factors (TFs) family is among of the largest and most diverse gene families found in plant species, and members of the bZIP TFs family perform important functions in plant developmental processes and stress response. To date, bZIP genes in Platycodon grandiflorus have not been characterized. In this work, a number of 47 PgbZIP genes were identified from the genome of P. grandiflorus, divided into 11 subfamilies. The distribution of these PgbZIP genes on the chromosome and gene replication events were analyzed. The motif, gene structure, cis-elements, and collinearity relationships of the PgbZIP genes were simultaneously analyzed. In addition, gene expression pattern analysis identified ten candidate genes involved in the developmental process of different tissue parts of P. grandiflorus. Among them, Four genes (PgbZIP5, PgbZIP21, PgbZIP25 and PgbZIP28) responded to drought and salt stress, which may have potential biological roles in P. grandiflorus development under salt and drought stress. Four hub genes (PgbZIP13, PgbZIP30, PgbZIP32 and PgbZIP45) mined in correlation network analysis, suggesting that these PgbZIP genes may form a regulatory network with other transcription factors to participate in regulating the growth and development of P. grandiflorus. This study provides new insights regarding the understanding of the comprehensive characterization of the PgbZIP TFs for further exploration of the functions of growth and developmental regulation in P. grandiflorus and the mechanisms for coping with abiotic stress response.

Genome-wide methylation, transcriptome and characteristic metabolites reveal the balance between diosgenin and brassinosteroids in Dioscorea zingiberensis

Diosgenin (DG) is a bioactive metabolite isolated from Dioscorea species, renowned for its medicinal properties. Brassinosteroids (BRs) are a class of crucial plant steroidal hormones. Cholesterol and campesterol are important intermediates of DG and BR biosynthesis, respectively. DG and BRs are structurally similar components; however, the regulatory network and metabolic interplays have not been fully elucidated. In an effort to decode these complex networks, we conducted a comprehensive study integrating genome-wide methylation, transcriptome and characteristic metabolite data from Dioscorea zingiberensis. Leveraging these data, we were able to construct a comprehensive regulatory network linking DG and BRs. Mass spectrometry results enabled us to clarify the alterations in cholesterol, campesterol, diosgenin, and castasterone (one of the major active BRs). The DG content decreased by 27.72% at 6 h after brassinolide treatment, whereas the content increased by 85.34% at 6 h after brassinazole treatment. Moreover, we pinpointed DG/BR-related genes, such as CASs, CYP90s, and B3-ARFs, implicated in the metabolic pathways of DG and BRs. Moreover, CASs and CYP90s exhibit hypomethylation, which is closely related to their high transcription. These findings provide robust evidence for the homeostasis between DG and BRs. In conclusion, our research revealed the existence of a balance between DG and BRs in D. zingiberensis. Furthermore, our work not only provides new insights into the relationship between the two pathways but also offers a fresh perspective on the functions of secondary metabolites.

Whole-genome sequencing in medicinal plants: current progress and prospect

Advancements in genomics have dramatically accelerated the research on medicinal plants, and the development of herbgenomics has promoted the "Project of 1K Medicinal Plant Genome" to decipher their genetic code. However, it is difficult to obtain their high-quality whole genomes because of the prevalence of polyploidy and/or high genomic heterozygosity. Whole genomes of 123 medicinal plants were published until September 2022. These published genome sequences were investigated in this review, covering their classification, research teams, ploidy, medicinal functions, and sequencing strategies. More than 1,000 institutes or universities around the world and 50 countries are conducting research on medicinal plant genomes. Diploid species account for a majority of sequenced medicinal plants. The whole genomes of plants in the Poaceae family are the most studied. Almost 40% of the published papers studied species with tonifying, replenishing, and heat-cleaning medicinal effects. Medicinal plants are still in the process of domestication as compared with crops, thereby resulting in unclear genetic backgrounds and the lack of pure lines, thus making their genomes more difficult to complete. In addition, there is still no clear routine framework for a medicinal plant to obtain a high-quality whole genome. Herein, a clear and complete strategy has been originally proposed for creating a high-quality whole genome of medicinal plants. Moreover, whole genome-based biological studies of medicinal plants, including breeding and biosynthesis, were reviewed. We also advocate that a research platform of model medicinal plants should be established to promote the genomics research of medicinal plants.

Integrative analysis of transcriptome and metabolome reveals the sesquiterpenoids and polyacetylenes biosynthesis regulation in Atractylodes lancea (Thunb.) DC

Atractylodes lancea (Thunb.) is a perennial medicinal herb, with its dry rhizomes are rich in various sesquiterpenoids and polyacetylenes components (including atractylodin, atractylon and β-eudesmol). However, the contents of these compounds are various and germplasms specific, and the mechanisms of biosynthesis in A. lancea are still unknown. In this study, we identified the differentially expressed candidate genes and metabolites involved in the biosynthesis of sesquiterpenoids and polyacetylenes, and speculated the anabolic pathways of these pharmaceutical components by transcriptome and metabolomic analysis. In the sesquiterpenoids biosynthesis, a total of 28 differentially expressed genes (DEGs) and 6 differentially expressed metabolites (DEMs) were identified. The beta-Selinene is likely to play a role in the synthesis of atractylon and β-eudesmol. Additionally, the polyacetylenes biosynthesis showed the presence of 3 DEGs and 4 DEMs. Notably, some fatty acid desaturase (FAB2 and FAD2) significantly down-regulated in polyacetylenes biosynthesis. The gamma-Linolenic acid is likely involved in the biosynthesis of polyacetylenes and thus further synthesis of atractylodin. Overall, these studies have investigated the biosynthetic pathways of atractylodin, atractylon and β-eudesmol in A. lancea for the first time, and present potential new anchor points for further exploration of sesquiterpenoids and polyacetylenes compound biosynthesis pathways in A. lancea.

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

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

Insights into Metabolic Engineering of Bioactive Molecules in Tetrastigma hemsleyanum Diels & Gilg: A Traditional Medicinal Herb

Plants are a vital source of bioactive molecules for various drug development processes. Tetrastigma hemsleyanum is one of the endangered medicinal plant species well known to the world due to its wide range of therapeutic effects. Many bioactive molecules have been identified from this plant, including many classes of secondary metabolites such as flavonoids, phenols, terpenoids, steroids, alkaloids, etc. Due to its slow growth, it usually takes 3-5 years to meet commercial medicinal materials for this plant. Also, T. hemsleyanum contains low amounts of specific bioactive compounds, which are challenging to isolate easily. Currently, scientists are attempting to increase bioactive molecules' production from medicinal plants in different ways or to synthesize them chemically. The genomic tools helped to understand medicinal plants' genome organization and led to manipulating genes responsible for various biosynthesis pathways. Metabolic engineering has made it possible to enhance the production of secondary metabolites by introducing manipulated biosynthetic pathways to attain high levels of desirable bioactive molecules. Metabolic engineering is a promising approach for improving the production of secondary metabolites over a short time period. In this review, we have highlighted the scope of various biotechnological approaches for metabolic engineering to enhance the production of secondary metabolites for pharmaceutical applications in T. hemsleyanum. Also, we summarized the progress made in metabolic engineering for bioactive molecule enhancement in T. hemsleyanum. It may lead to reducing the destruction of the natural habitat of T. hemsleyanum and conserving them through the cost-effective production of bioactive molecules in the future.

Heterologous Expression of Platycodon grandiflorus PgF3'5'H Modifies Flower Color Pigmentation in Tobacco

Flavonoid-3',5'-hydroxylase (F3'5'H) is the key enzyme for the biosynthesis of delphinidin-based anthocyanins, which are generally required for purple or blue flowers. Previously, we isolated a full-length cDNA of PgF3'5'H from Platycodon grandiflorus, which shared the highest homology with Campanula medium F3'5'H. In this study, PgF3'5'H was subcloned into a plant over-expression vector and transformed into tobacco via Agrobacterium tumefaciens to investigate its catalytic function. Positive transgenic tobacco T0 plants were obtained by hygromycin resistance screening and PCR detection. PgF3'5'H showed a higher expression level in all PgF3'5'H transgenic tobacco plants than in control plants. Under the drive of the cauliflower mosaic virus (CaMV) 35S promoter, the over-expressed PgF3'5'H produced dihydromyricetin (DHM) and some new anthocyanin pigments (including delphinidin, petunidin, peonidin, and malvidin derivatives), and increased dihydrokaempferol (DHK), taxifolin, tridactyl, cyanidin derivatives, and pelargonidin derivatives in PgF3'5'H transgenic tobacco plants by ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) analysis, resulting in a dramatic color alteration from light pink to magenta. These results indicate that PgF3'5'H products have F3'5'H enzyme activity. In addition, PgF3'5'H transfer alters flavonoid pigment synthesis and accumulation in tobacco. Thus, PgF3'5'H may be considered a candidate gene for gene engineering to enhance anthocyanin accumulation and the molecular breeding project for blue flowers.

Histone Deacetylase Inhibitor, Sodium Butyrate-Induced Metabolic Modulation in Platycodon grandiflorus Roots Enhances Anti-Melanogenic Properties

While the status of histone acetylation is a critical regulator of chromatin's structure with a significant impact on plant physiology, our understanding of epigenetic regulation in the biosynthesis of active compounds in plants is limited. In this study, Platycodon grandiflorus was treated with sodium butyrate (NaB), a histone deacetylase inhibitor, to investigate the influence of histone acetylation on secondary metabolism. Its treatment with NaB increased the acetylation of histone H3 at lysine 9, 14, and 27 and enhanced the anti-melanogenic properties of P. grandiflorus roots. Through transcriptome and differentially expressed gene analyses, we found that NaB influenced the expression of genes that were involved in both primary and secondary metabolic pathways. In addition, NaB treatment caused the accumulation of polyphenolic compounds, including dihydroquercetin, gallic acid, and 2,4-dihydroxybenzoic acid. The NaB-induced transcriptional activation of genes in the phenylpropanoid biosynthetic pathway influenced the anti-melanogenic properties of P. grandiflorus roots. Overall, these findings suggest the potential of an epigenomic approach to enhance the medicinal qualities of medicinal plants.

Multilayered regulation of secondary metabolism in medicinal plants

Medicinal plants represent a huge reservoir of secondary metabolites (SMs), substances with significant pharmaceutical and industrial potential. However, obtaining secondary metabolites remains a challenge due to their low-yield accumulation in medicinal plants; moreover, these secondary metabolites are produced through tightly coordinated pathways involving many spatiotemporally and environmentally regulated steps. The first regulatory layer involves a complex network of transcription factors; a second, more recently discovered layer of complexity in the regulation of SMs is epigenetic modification, such as DNA methylation, histone modification and small RNA-based mechanisms, which can jointly or separately influence secondary metabolites by regulating gene expression. Here, we summarize the findings in the fields of genetic and epigenetic regulation with a special emphasis on SMs in medicinal plants, providing a new perspective on the multiple layers of regulation of gene expression.

Elucidation of Geniposide and Crocin Accumulation and Their Biosysnthsis-Related Key Enzymes during Gardenia jasminoides Fruit Growth

Gardenia jasminoides fruits are extensively grown worldwide, with a large harvest, and its major medicinal ingredients are geniposide and crocins. Research on their accumulation and biosynthsis-related enzymes is rare. In this study, the accumulation of geniposide and crocin of G. jasminoides fruits at different developmental stages were clarified by HPLC. The highest cumulative amount of geniposide was 2.035% during the unripe-fruit period, and the highest content of crocin was 1.098% during the mature-fruit period. Furthermore, transcriptome sequencing was performed. A total of 50 unigenes encoding 4 key enzymes related in geniposide biosynthsis pathways were screened, and 41 unigenes encoding 7 key enzymes in the pathways of crocin were elucidated. It was found that the expression levels of differentially expressed genes of DN67890_c0_g1_i2-encoding GGPS, which is highly related to geniposide biosynthesis, and DN81253_c0_g1_i1-encoding lcyB, DN79477_c0_g1_i2-encoding lcyE, and DN84975_c1_g7_i11-encoding CCD, which are highly related to crocin biosynthesis, were consistent with the accumulation of geniposide and crocin content, respectively. The qRT-PCR results showed that the trends of relative expression were consistent with transcribed genes. This study provides insights for understanding the geniposide and crocin accumulation and biosynthsis during fruit development in G. jasminoides.

Physiological and sucrose metabolic responses to waterlogging stress in balloon flower (Platycodon grandiflorus (Jacq.) A. DC)

Waterlogging stress is a major limiting factor resulting in stunted growth and loss of crop productivity, especially for root crops. However, physiological responses to waterlogging have been studied in only a few plant models. To gain insight into how balloon flower (Platycodon grandiflorus (Jacq.) A. DC) responds to waterlogging stress, we investigate changes to sucrose metabolism combined with a physiological analysis. Although waterlogging stress decreased the photosynthetic rate in balloon flower, leaves exhibited an increase in glucose (ninefold), fructose (4.7-fold), and sucrose (2.1-fold), indicating inhibition of sugar transport via the phloem. In addition, roots showed a typical response to hypoxia, such as the accumulation of proline (4.5-fold higher than in control roots) and soluble sugars (2.1-fold higher than in control roots). The activities and expression patterns of sucrose catabolizing enzymes suggest that waterlogging stress leads to a shift in the pathway of sucrose degradation from invertase to sucrose synthase (Susy), which consumes less ATP. Furthermore, we suggest that the waterlogging-stress-induced gene PlgSusy1 encodes the functional Susy enzyme, which may contribute to improving tolerance in balloon flower to waterlogging. As a first step toward understanding the waterlogging-induced regulatory mechanisms in balloon flower, we provide a solid foundation for further understanding waterlogging-induced alterations in source-sink relationships.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01310-y.

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

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

De Novo Transcriptome Sequencing of Codonopsis lanceolata for Identification of Triterpene Synthase and Triterpene Acetyltransferase

Codonopsis lanceolata (Campanulaceae) is a perennial plant commonly known as the bonnet bellflower. This species is widely used in traditional medicine and is considered to have multiple medicinal properties. In this study, we found that shoots and roots of C. lanceolata contained various types of free triterpenes (taraxerol, β-amyrin, α-amyrin, and friedelin) and triterpene acetates (taraxerol acetate, β-amyrin acetate, and α-amyrin acetate). The content of triterpenes and triterpene acetates by GC analysis was higher in the shoot than in the roots. To investigate the transcriptional activity of genes involved in triterpenes and triterpene acetate biosynthesis, we performed de novo transcriptome analysis of shoots and roots of C. lanceolata by sequencing using the Illumina platform. A total of 39,523 representative transcripts were obtained. After functional annotation of the transcripts, the differential expression of genes involved in triterpene biosynthetic pathways was investigated. Generally, the transcriptional activity of unigenes in the upstream region (MVA and MEP pathway) of triterpene biosynthetic pathways was higher in shoots than in roots. Various triterpene synthases (2,3-oxidosqualene cyclase, OSC) participate to produce triterpene skeletons by the cyclization of 2,3-oxidosqualene. A total of fifteen contigs were obtained in annotated OSCs in the representative transcripts. Functional characterization of four OSC sequences by heterologous expression in yeast revealed that ClOSC1 was determined as taraxerol synthase, and ClOSC2 was a mixed-amyrin synthase producing α-amyrin and β-amyrin. Five putative contigs of triterpene acetyltransferases showed high homology to the lettuce triterpene acetyltransferases. Conclusively, this study provides the basis of molecular information, particularly for the biosynthesis of triterpenes and triterpene acetates in C. lanceolata.

Genome-Wide Identification and Expression Profiling Analysis of the Trihelix Gene Family and response of PgGT1 under Abiotic Stresses in Platycodon grandiflorus

The trihelix gene family plays an important role in plant growth and abiotic stress responses. Through the analysis of genomic and transcriptome data, 35 trihelix family members were identified for the first time in Platycodon grandiflorus; they were classified into five subfamilies: GT-1, GT-2, SH4, GTγ, and SIP1. The gene structure, conserved motifs and evolutionary relationships were analyzed. Prediction of physicochemical properties of the 35 trihelix proteins founded, the number of amino acid molecules is between 93 and 960, theoretical isoelectric point is between 4.24 and 9.94, molecular weight is between 9829.77 and 107435.38, 4 proteins among them were stable, and all GRAVY is negative. The full-length cDNA sequence of the PgGT1 gene of the GT-1 subfamily was cloned by PCR. It is a 1165 bp ORF encoding a 387 amino acid protein, with a molecular weight of 43.54 kDa. The predicted subcellular localization of the protein in the nucleus was experimentally verified. After being treated with NaCl, PEG6000, MeJA, ABA, IAA, SA, and ethephon, the expression of PgGT1 gene showed an up-regulated trend except for the roots treated with NaCl and ABA. This study laid a bioinformatics foundation for the research of trihelix gene family and the cultivation of excellent germplasm of P. grandiflorus.

PlgMYBR1, an R2R3-MYB transcription factor, plays as a negative regulator of anthocyanin biosynthesis in Platycodon grandiflorus

Floral color plays a major role in pollinator specificity, and changes in color may result in pollinator shifts and pollinator-mediated speciation. In the purple flowers of Platycodon grandiflorus, anthocyanins are the major pigment metabolites, whereas white flowers result due to the absence of anthocyanins. The lack of anthocyanins may be due to the inhibition of the anthocyanin biosynthesis pathway. However, the molecular mechanism of anthocyanin biosynthesis in P. grandiflorus is not fully understood. Hence, we identified R2R3-MYB transcription factor, PlgMYBR1, as a negative regulator for anthocyanin biosynthesis using sequence homology and tissue-specific expression pattern analyses. A heterologous co-expression assay suggested that PlgMYBR1 inhibited the function of AtPAP1 (Arabidopsis thaliana production of anthocyanin pigment 1), indicating that PlgMYBR1 plays as a repressor of anthocyanin biosynthesis in P. grandiflorus. Our results provide a foundation for future efforts to understand the anthocyanin biosynthesis in P. grandiflorus and, thereby, to improve flower color through genetic engineering.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03490-6.

Creating novel ornamentals via new strategies in the era of genome editing

Ornamental breeding has traditionally focused on improving novelty, yield, quality, and resistance to biotic or abiotic stress. However, achieving these goals has often required laborious crossbreeding, while precise breeding techniques have been underutilized. Fortunately, recent advancements in plant genome sequencing and editing technology have opened up exciting new frontiers for revolutionizing ornamental breeding. In this review, we provide an overview of the current state of ornamental transgenic breeding and propose four promising breeding strategies that have already proven successful in crop breeding and could be adapted for ornamental breeding with the help of genome editing. These strategies include recombination manipulation, haploid inducer creation, clonal seed production, and reverse breeding. We also discuss in detail the research progress, application status, and feasibility of each of these tactics.

The Current Developments in Medicinal Plant Genomics Enabled the Diversification of Secondary Metabolites' Biosynthesis

Medicinal plants produce important substrates for their adaptation and defenses against environmental factors and, at the same time, are used for traditional medicine and industrial additives. Plants have relatively little in the way of secondary metabolites via biosynthesis. Recently, the whole-genome sequencing of medicinal plants and the identification of secondary metabolite production were revolutionized by the rapid development and cheap cost of sequencing technology. Advances in functional genomics, such as transcriptomics, proteomics, and metabolomics, pave the way for discoveries in secondary metabolites and related key genes. The multi-omics approaches can offer tremendous insight into the variety, distribution, and development of biosynthetic gene clusters (BGCs). Although many reviews have reported on the plant and medicinal plant genome, chemistry, and pharmacology, there is no review giving a comprehensive report about the medicinal plant genome and multi-omics approaches to study the biosynthesis pathway of secondary metabolites. Here, we introduce the medicinal plant genome and the application of multi-omics tools for identifying genes related to the biosynthesis pathway of secondary metabolites. Moreover, we explore comparative genomics and polyploidy for gene family analysis in medicinal plants. This study promotes medicinal plant genomics, which contributes to the biosynthesis and screening of plant substrates and plant-based drugs and prompts the research efficiency of traditional medicine.

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

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

Integrated Metabolome and Lipidome Strategy to Reveal the Action Pattern of Paclobutrazol, a Plant Growth Retardant, in Varying the Chemical Constituents of Platycodon Root

Platycodon root, a medicinal food homology species which has been used in Asian countries for hundreds of years, is now widely cultivated in China. Treatment with paclobutrazol, a typical plant growth retardant, has raised uncertainties regarding the quality of Platycodon root, which have been rarely investigated. In the present study, metabolomic and lipidomic differences were revealed by ultra-high performance liquid chromatography coupled to ion mobility-quadrupole time of flight mass spectrometry (UPLC-IM-QTOF-MS). A significant decrease of platycodigenin-type saponins was observed in the paclobutrazol-treated sample. Carrying out a comprehensive quantitative analysis, the contents of total saponins and saccharides were determined to illustrate the mode of action of paclobutrazol on Platycodon root. This study demonstrated an exemplary research model in explaining how the exogenous matter influences the chemical properties of medicinal plants, and therefore might provide insights into the reasonable application of plant growth regulators.

Molecular Cloning and Functional Characterization of a β-Glucosidase Gene to Produce Platycodin D in Platycodon grandiflorus

Platycodin D (PD) is a deglycosylated triterpene saponin with much higher pharmacological activity than glycosylated platycoside E (PE). Extensive studies in vitro showed that the transformation of platycoside E to platycodin D can be achieved using β-glucosidase extracted from several bacteria. However, whether similar enzymes in Platycodon grandiflorus could convert platycoside E to platycodin D, as well as the molecular mechanism underlying the deglycosylation process of platycodon E, remain unclear. Here, we identified a β-glucosidase in P. grandiflorus from our previous RNA-seq analysis, with a full-length cDNA of 1,488 bp encoding 495 amino acids. Bioinformatics and phylogenetic analyses showed that β-glucosidases in P. grandiflorus have high homology with other plant β-glucosidases. Subcellular localization showed that there is no subcellular preference for its encoding gene. β-glucosidase was successfully expressed as 6 × His-tagged fusion protein in Escherichia coli BL21 (DE3). Western blot analysis yielded a recombinant protein of approximately 68 kDa. In vitro enzymatic reactions determined that β-glucosidase was functional and could convert PE to PD. RT-qPCR analysis showed that the expression level of β-glucosidase was higher at night than during the day, with the highest expression level between 9:00 and 12:00 at night. Analysis of the promoter sequence showed many light-responsive cis-acting elements, suggesting that the light might regulate the gene. The results will contribute to the further study of the biosynthesis and metabolism regulation of triterpenoid saponins in P. grandiflorus.

Whole-Genome Identification and Analysis of Multiple Gene Families Reveal Candidate Genes for Theasaponin Biosynthesis in Camellia oleifera

Camellia oleifera is an economically important oilseed tree. Seed meals of C. oleifera have a long history of use as biocontrol agents in shrimp farming and as cleaning agents in peoples’ daily lives due to the presence of theasaponins, the triterpene saponins from the genus Camellia. To characterize the biosynthetic pathway of theasaponins in C. oleifera, members of gene families involved in triterpenoid biosynthetic pathways were identified and subjected to phylogenetic analysis with corresponding members in Arabidopsis thaliana, Camellia sinensis, Actinidia chinensis, Panax ginseng, and Medicago truncatula. In total, 143 triterpenoid backbone biosynthetic genes, 1169 CYP450s, and 1019 UGTs were identified in C. oleifera. The expression profiles of triterpenoid backbone biosynthetic genes were analyzed in different tissue and seed developmental stages of C. oleifera. The results suggested that MVA is the main pathway for triterpenoid backbone biosynthesis. Moreover, the candidate genes for theasaponin biosynthesis were identified by WGCNA and qRT-PCR analysis; these included 11 CYP450s, 14 UGTs, and eight transcription factors. Our results provide valuable information for further research investigating the biosynthetic and regulatory network of theasaponins.

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

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

Identification and Functional Characterization of FLOWERING LOCUS T in Platycodon grandiflorus

Platycodon grandiflorus roots have been used as a foodstuff and traditional medicine for thousands of years in East Asia. In order to increase the root development of P. grandiflorus, cultivators removed the inflorescences, suggesting the possible negative effect of flowering on root development. This indicates that the genetic improvement of P. grandiflorus by late flowering is a potential approach to increase productivity. However, nothing is known about key genes integrating multiple flowering pathways in P. grandiflorus. In order to fill this gap, we identified potential homologs of the FLOWERING LOCUS T (FT) gene in P. grandiflorus. The alignment with other FT members and phylogenetic analysis revealed that the P. grandiflorus FT (PlgFT) protein contains highly conserved functional domains and belongs to the FT-like clade. The expression analysis revealed spatial variations in the transcription of PlgFT in different organs. In addition, the expression level of PlgFT was increased by high temperature but not by photoperiodic light input signals, presumably due to lacking the CONSTANS binding motif in its promoter region. Furthermore, PlgFT induced early flowering upon its overexpression in P. grandiflorus, suggesting the functional role of PlgFT in flowering. Taken together, we functionally characterized PlgFT as a master regulator of P. grandiflorus flowering under inductive high temperature, which will serve as an important target gene for improving the root productivity.

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

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

The Euscaphis japonica genome and the evolution of malvids

Malvids is one of the largest clades of rosids, includes 58 families and exhibits remarkable morphological and ecological diversity. Here, we report a high-quality chromosome-level genome assembly for Euscaphis japonica, an early-diverging species within malvids. Genome-based phylogenetic analysis suggests that the unstable phylogenetic position of E. japonica may result from incomplete lineage sorting and hybridization event during the diversification of the ancestral population of malvids. Euscaphis japonica experienced two polyploidization events: the ancient whole genome triplication event shared with most eudicots (commonly known as the γ event) and a more recent whole genome duplication event, unique to E. japonica. By resequencing 101 samples from 11 populations, we speculate that the temperature has led to the differentiation of the evergreen and deciduous of E. japonica and the completely different population histories of these two groups. In total, 1012 candidate positively selected genes in the evergreen were detected, some of which are involved in flower and fruit development. We found that reddening and dehiscence of the E. japonica pericarp and long fruit-hanging time promoted the reproduction of E. japonica populations, and revealed the expression patterns of genes related to fruit reddening, dehiscence and abscission. The key genes involved in pentacyclic triterpene synthesis in E. japonica were identified, and different expression patterns of these genes may contribute to pentacyclic triterpene diversification. Our work sheds light on the evolution of E. japonica and malvids, particularly on the diversification of E. japonica and the genetic basis for their fruit dehiscence and abscission.

The first high-quality chromosomal genome assembly of a medicinal and edible plant Arctium lappa

Arctium lappa has a long medicinal and edible history with great economic importance. Here, the first high-quality chromosome-level draft genome of A. lappa was presented by the Illumina and PacBio sequencing data. The assembled genome was approximately 1.79 Gb with a N50 contig size of 6.88 Mb. Approximately 1.70 Gb (95.4%) of the contig sequences were anchored onto 18 chromosomes using Hi-C data; the scaffold N50 was improved to be 91.64 Mb. Furthermore, we obtained 1.12 Gb (68.46%) of repetitive sequences and 32,771 protein-coding genes; 616 positively selected candidate genes were identified. Among candidate genes related to lignan biosynthesis, the following were found to be highly correlated with the accumulation of arctiin: 4-coumarate-CoA ligase (4CL), dirigent protein (DIR), and hydroxycinnamoyl transferase (HCT). Additionally, we compared the transcriptomes of A. lappa roots at three different developmental stages and identified 8,943 differentially expressed genes (DEGs) in these tissues. These data can be utilized to identify genes related to A. lappa quality or provide a basis for molecular identification and comparative genomics among related species.

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

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

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.

Machine learning, transcriptome, and genotyping chip analyses provide insights into SNP markers identifying flower color in Platycodon grandiflorus

Bellflower is an edible ornamental gardening plant in Asia. For predicting the flower color in bellflower plants, a transcriptome-wide approach based on machine learning, transcriptome, and genotyping chip analyses was used to identify SNP markers. Six machine learning methods were deployed to explore the classification potential of the selected SNPs as features in two datasets, namely training (60 RNA-Seq samples) and validation (480 Fluidigm chip samples). SNP selection was performed in sequential order. Firstly, 96 SNPs were selected from the transcriptome-wide SNPs using the principal compound analysis (PCA). Then, 9 among 96 SNPs were later identified using the Random forest based feature selection method from the Fluidigm chip dataset. Among six machines, the random forest (RF) model produced higher classification performance than the other models. The 9 SNP marker candidates selected for classifying the flower color classification were verified using the genomic DNA PCR with Sanger sequencing. Our results suggest that this methodology could be used for future selection of breeding traits even though the plant accessions are highly heterogeneous.

Review on the Development and Applications of Medicinal Plant Genomes

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

Developmental Methylome of the Medicinal Plant Catharanthus roseus Unravels the Tissue-Specific Control of the Monoterpene Indole Alkaloid Pathway by DNA Methylation

Catharanthus roseus produces a wide spectrum of monoterpene indole alkaloids (MIAs). MIA biosynthesis requires a tightly coordinated pathway involving more than 30 enzymatic steps that are spatio-temporally and environmentally regulated so that some MIAs specifically accumulate in restricted plant parts. The first regulatory layer involves a complex network of transcription factors from the basic Helix Loop Helix (bHLH) or AP2 families. In the present manuscript, we investigated whether an additional epigenetic layer could control the organ-, developmental- and environmental-specificity of MIA accumulation. We used Whole-Genome Bisulfite Sequencing (WGBS) together with RNA-seq to identify differentially methylated and expressed genes among nine samples reflecting different plant organs and experimental conditions. Tissue specific gene expression was associated with specific methylation signatures depending on cytosine contexts and gene parts. Some genes encoding key enzymatic steps from the MIA pathway were found to be simultaneously differentially expressed and methylated in agreement with the corresponding MIA accumulation. In addition, we found that transcription factors were strikingly concerned by DNA methylation variations. Altogether, our integrative analysis supports an epigenetic regulation of specialized metabolisms in plants and more likely targeting transcription factors which in turn may control the expression of enzyme-encoding genes.