Global Identification of the Full-Length Transcripts and Alternative Splicing Related to Phenolic Acid Biosynthetic Genes in Salvia miltiorrhiza

Plant-Derived Caffeic Acid and Its Derivatives: An Overview of Their NMR Data and Biosynthetic Pathways

In recent years, caffeic acid and its derivatives have received increasing attention due to their obvious physiological activities and wide distribution in nature. In this paper, to clarify the status of research on plant-derived caffeic acid and its derivatives, nuclear magnetic resonance spectroscopy data and possible biosynthetic pathways of these compounds were collected from scientific databases (SciFinder, PubMed and China Knowledge). According to different types of substituents, 17 caffeic acid and its derivatives can be divided into the following classes: caffeoyl ester derivatives, caffeyltartaric acid, caffeic acid amide derivatives, caffeoyl shikimic acid, caffeoyl quinic acid, caffeoyl danshens and caffeoyl glycoside. Generalization of their 13C-NMR and 1H-NMR data revealed that acylation with caffeic acid to form esters involves acylation shifts, which increase the chemical shift values of the corresponding carbons and decrease the chemical shift values of the corresponding carbons of caffeoyl. Once the hydroxyl group is ester, the hydrogen signal connected to the same carbon shifts to the low field (1.1~1.6). The biosynthetic pathways were summarized, and it was found that caffeic acid and its derivatives are first synthesized in plants through the shikimic acid pathway, in which phenylalanine is deaminated to cinnamic acid and then transformed into caffeic acid and its derivatives. The purpose of this review is to provide a reference for further research on the rapid structural identification and biofabrication of caffeic acid and its derivatives.

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.

Applications of some advanced sequencing, analytical, and computational approaches in medicinal plant research: a review

The potential active chemicals found in medicinal plants, which have long been employed as natural medicines, are abundant. Exploring the genes responsible for producing these compounds has given new insights into medicinal plant research. Previously, the authentication of medicinal plants was done via DNA marker sequencing. With the advancement of sequencing technology, several new techniques like next-generation sequencing, single molecule sequencing, and fourth-generation sequencing have emerged. These techniques enshrined the role of molecular approaches for medicinal plants because all the genes involved in the biosynthesis of medicinal compound(s) could be identified through RNA-seq analysis. In several research insights, transcriptome data have also been used for the identification of biosynthesis pathways. miRNAs in several medicinal plants and their role in the biosynthesis pathway as well as regulation of the disease-causing genes were also identified. In several research articles, an in silico study was also found to be effective in identifying the inhibitory effect of medicinal plant-based compounds against virus' gene(s). The use of advanced analytical methods like spectroscopy and chromatography in metabolite proofing of secondary metabolites has also been reported in several recent research findings. Furthermore, advancement in molecular and analytic methods will give new insight into studying the traditionally important medicinal plants that are still unexplored.

Application of third-generation sequencing to herbal genomics

There is a long history of traditional medicine use. However, little genetic information is available for the plants used in traditional medicine, which limits the exploitation of these natural resources. Third-generation sequencing (TGS) techniques have made it possible to gather invaluable genetic information and develop herbal genomics. In this review, we introduce two main TGS techniques, PacBio SMRT technology and Oxford Nanopore technology, and compare the two techniques against Illumina, the predominant next-generation sequencing technique. In addition, we summarize the nuclear and organelle genome assemblies of commonly used medicinal plants, choose several examples from genomics, transcriptomics, and molecular identification studies to dissect the specific processes and summarize the advantages and disadvantages of the two TGS techniques when applied to medicinal organisms. Finally, we describe how we expect that TGS techniques will be widely utilized to assemble telomere-to-telomere (T2T) genomes and in epigenomics research involving medicinal plants.

Seven non-differentially expressed 'dark biomarkers' show transcriptional dysregulation in chronic lymphocytic leukemia

Aim: Transcriptional regulation is actively involved in the onset and progression of various diseases. This study used the feature-engineering approach model-based quantitative transcription regulation to quantitatively measure the correlation between mRNA and transcription factors in a reference dataset of chronic lymphocytic leukemia (CLL) transcriptomes. Methods: A comprehensive investigation of transcriptional regulation changes in CLL was conducted using 973 samples in six independent datasets. Results & conclusion: Seven mRNAs were detected to have significantly differential model-based quantitative transcription regulation values but no differential expression between CLL patients and controls. We called these genes 'dark biomarkers' because their original expression levels did not show differential changes in the CLL patients. The overlapping lncRNAs might have contributed their transcripts to the expression miscalculations of these dark biomarkers.

Overexpression of SmLAC25 promotes lignin accumulation and decreases salvianolic acid content in Salvia miltiorrhiza

Laccase (LAC) is a blue multicopper oxidase that contains four copper ions, which is involved in lignin polymerization and flavonoid biosynthesis in plants. Although dozens of LAC genes have been identified in Salvia miltiorrhiza Bunge (a model medicinal plant), most have not been functionally characterized. Here, we explored the expression patterns and the functionality of SmLAC25 in S. miltiorrhiza. SmLAC25 has a higher expression level in roots and responds to methyl jasmonate, auxin, abscisic acid, and gibberellin stimuli. The SmLAC25 protein is localized in the cytoplasm and chloroplasts. Recombinant SmLAC25 protein could oxidize coniferyl alcohol and sinapyl alcohol, two monomers of G-lignin and S-lignin. To investigate its function, we generated SmLAC25-overexpressed S. miltiorrhiza plantlets and hairy roots. The lignin content increased significantly in all SmLAC25-overexpressed plantlets and hairy roots, compared with the controls. However, the concentrations of rosmarinic acid and salvianolic acid B decreased significantly in all the SmLAC25-overexpressed lines. Further studies revealed that the transcription levels of some key enzyme genes in the lignin synthesis pathway (e.g., SmCCR and SmCOMT) were significantly improved in the SmLAC25-overexpressed lines, while the expression levels of multiple enzyme genes in the salvianolic acid biosynthesis pathway were inhibited. We speculated that the overexpression of SmLAC25 promoted the metabolic flux of lignin synthesis, which resulted in a decreased metabolic flux to the salvianolic acid biosynthesis pathway.

A high-resolution single-molecule sequencing-based Arabidopsis transcriptome using novel methods of Iso-seq analysis

BACKGROUND

Accurate and comprehensive annotation of transcript sequences is essential for transcript quantification and differential gene and transcript expression analysis. Single-molecule long-read sequencing technologies provide improved integrity of transcript structures including alternative splicing, and transcription start and polyadenylation sites. However, accuracy is significantly affected by sequencing errors, mRNA degradation, or incomplete cDNA synthesis.

RESULTS

We present a new and comprehensive Arabidopsis thaliana Reference Transcript Dataset 3 (AtRTD3). AtRTD3 contains over 169,000 transcripts-twice that of the best current Arabidopsis transcriptome and including over 1500 novel genes. Seventy-eight percent of transcripts are from Iso-seq with accurately defined splice junctions and transcription start and end sites. We develop novel methods to determine splice junctions and transcription start and end sites accurately. Mismatch profiles around splice junctions provide a powerful feature to distinguish correct splice junctions and remove false splice junctions. Stratified approaches identify high-confidence transcription start and end sites and remove fragmentary transcripts due to degradation. AtRTD3 is a major improvement over existing transcriptomes as demonstrated by analysis of an Arabidopsis cold response RNA-seq time-series. AtRTD3 provides higher resolution of transcript expression profiling and identifies cold-induced differential transcription start and polyadenylation site usage.

CONCLUSIONS

AtRTD3 is the most comprehensive Arabidopsis transcriptome currently. It improves the precision of differential gene and transcript expression, differential alternative splicing, and transcription start/end site usage analysis from RNA-seq data. The novel methods for identifying accurate splice junctions and transcription start/end sites are widely applicable and will improve single-molecule sequencing analysis from any species.

Alternative Splicing Dynamics During the Lifecycle of Salvia miltiorrhiza Root Revealed the Fine Tuning in Root Development and Ingredients Biosynthesis

Alternative splicing (AS) is an essential post-transcriptional process that enhances the coding and regulatory potential of the genome, thereby strongly influencing multiple plant physiology processes, such as metabolic biosynthesis. To explore how AS affects the root development and synthesis of tanshinones and phenolic acid pathways in Salvia miltiorrhiza roots, we investigated the dynamic landscape of AS events in S. miltiorrhiza roots during an annual life history. Temporal profiling represented a distinct temporal variation of AS during the entire development stages, showing the most abundant AS events at the early seedling stage (ES stage) and troughs in 45 days after germination (DAG) and 120 DAG. Gene ontology (GO) analysis indicated that physiological and molecular events, such as lateral root formation, gravity response, RNA splicing regulation, and mitogen-activated protein kinase (MAPK) cascade, were greatly affected by AS at the ES stage. AS events were identified in the tanshinones and phenolic acids pathways as well, especially for the genes for the branch points of the pathways as SmRAS and SmKSL1. Fifteen Ser/Arg-rich (SR) proteins and eight phosphokinases (PKs) were identified with high transcription levels at the ES stage, showing their regulatory roles for the high frequency of AS in this stage. Simultaneously, a co-expression network that includes 521 highly expressed AS genes, SRs, and PKs, provides deeper insight into the mechanism for the variable programming of AS.

Metabolomics integrated with transcriptomics reveals the distribution of iridoid and crocin metabolic flux in Gardenia jasminoides Ellis

Gardenia jasminoides Ellis (G. jasminoides) fruits are used as a resource for obtaining natural colorants and in traditional Chinese herbal medicine. However, G. jasminoides presents a relatively long flowering period and different ripening periods, so there are significant differences in the accumulation of metabolites in fruits of different colors. In addition, the complete metabolic pathways of iridoidsand crocins, which are used as medicinal composition of G. jasminoides, are poorly understood at present. In this research, we comprehensively compared the transcriptome and metabolites profiles of the developmental stages and locations of iridoid and crocin biosynthesis. A large number of differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs) were detected in four groups of samples, and clear variation in the pattern of metabolite abundance and gene expression were observed among different fruit colors and parts. Geniposide and gardenoside mainly accumulated in the sarcocarp of green fruit (GFS) and the sarcocarp of red fruit (FS), respectively. Crocin mainly accumulated in the peel and sarcocarp of red fruits. In the iridoid pathway, we hypothesized that there was a transport mechanism from the sarcocarp to the peel of G. jasminoides because of the inconsistent expression of G8O, 10-HGO and IS associated with differences in fruit ripening. UGTs play an important role in the biosynthesis of the active components of G. jasminoides. Combined transcriptome and metabonomics analysis showed a negative correlation between the biosynthesis of geniposide and crocin. The redirection of the metabolic flux and the regulation of key enzymes may be the main reasons for the changes in the biosynthesis of iridoid and crocin in G. jasminoides fruit. Our study expended valuable information for functional genomic library and provided new insights for metabolic engineering of secondary metabolite in G. Jasminoides.

Transcription factor SmSPL7 promotes anthocyanin accumulation and negatively regulates phenolic acid biosynthesis in Salvia miltiorrhiza

Plant-specific SQUAMOSA promoter-binding protein-like (SPL) transcription factors play critical regulatory roles during plant growth and development. However, the functions of SPLs in Salvia miltiorrhiza (SmSPLs; a model medicinal plant) have not been reported. Here, the expression patterns and functions of SmSPL7 were characterized in S. miltiorrhiza. SmSPL7 was expressed in all parts of S. miltiorrhiza, with the highest expression level in the leaves, and could be inhibited by multiple hormones, including methyl jasmonate, auxin, abscisic acid, and gibberellin. SmSPL7 is localized within the nucleus and exhibits robust transcriptional activation activity. Transgenic lines overexpressing SmSPL7 demonstrated pronounced growth inhibition, accompanied by increased anthocyanin accumulation via the genetic activation of the anthocyanin biosynthesis pathway. However, SmSPL7 overexpression significantly decreased salvianolic acid B (SalB) production by inhibiting the transcripts of genes implicated in its biosynthesis pathway. Further analysis indicated that SmSPL7 directly binds to SmTAT1 and Sm4CL9 promoters and blocks their expression to inhibit the biosynthesis of SalB. Taken together, these results indicate that SmSPL7 is a negative regulator of SalB biosynthesis but positively regulates anthocyanin accumulation in S. miltiorrhiza. These findings provide new insights into the functionality of the SPL family while establishing an important foundation for further uncovering the crucial roles of SmSPL7 in the growth of S. miltiorrhiza.

Full-length transcriptome analysis of pecan (Carya illinoinensis) kernels

Pecan is rich in bioactive components such as fatty acids (FAs) and flavonoids and is an important nut type worldwide. Therefore, the molecular mechanisms of phytochemical biosynthesis in pecan are a focus of research. Recently, a draft genome and several transcriptomes have been published. However, the full-length mRNA transcripts remain unclear, and the regulatory mechanisms behind the quality components biosynthesis and accumulation have not been fully investigated. In this study, single-molecule long-read sequencing technology was used to obtain full-length transcripts of pecan kernels. In total, 37,504 isoforms of 16,702 genes were mapped to the reference genome. The numbers of known isoforms, new isoforms, and novel isoforms were 9013 (24.03%), 26,080 (69.54%), and 2411 (6.51%), respectively. Over 80% of the transcripts (30,751, 81.99%) had functional annotations. A total of 15,465 alternative splicing (AS) events and 65,761 alternative polyadenylation events were detected; wherein, the retained intron was the predominant type (5652, 36.55%) of AS. Furthermore, 1894 long noncoding RNAs and 1643 transcription factors were predicted using bioinformatics methods. Finally, the structural genes associated with FA and flavonoid biosynthesis were characterized. A high frequency of AS accuracy (70.31%) was observed in FA synthesis-associated genes. This study provides a full-length transcriptome data set of pecan kernels, which will significantly enhance the understanding of the regulatory basis of phytochemical biosynthesis during pecan kernel maturation.

SmSPL6 Induces Phenolic Acid Biosynthesis and Affects Root Development in Salvia miltiorrhiza

Salvia miltiorrhiza is a renowned model medicinal plant species for which 15 SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) family genes have been identified; however, the specific functions of SmSPLs have not been well characterized as of yet. For this study, the expression patterns of SmSPL6 were determined through its responses to treatments of exogenous hormones, including indole acetic acid (IAA), gibberellic acid (GA₃), methyl jasmonic acid (MeJA), and abscisic acid (ABA). To characterize its functionality, we obtained SmSPL6-ovexpressed transgenic S. miltiorrhiza plants and found that overexpressed SmSPL6 promoted the accumulation of phenolic acids and repressed the biosynthesis of anthocyanin. Meanwhile, the root lengths of the SmSPL6-overexpressed lines were significantly longer than the control; however, both the fresh weights and lateral root numbers decreased. Further investigations indicated that SmSPL6 regulated the biosynthesis of phenolic acid by directly binding to the promoter regions of the enzyme genes Sm4CL9 and SmCYP98A14 and activated their expression. We concluded that SmSPL6 regulates not only the biosynthesis of phenolic acids, but also the development of roots in S. miltiorrhiza.

Insights into salvianolic acid B biosynthesis from chromosome-scale assembly of the Salvia bowleyana genome

Salvia bowleyana is a traditional Chinese medicinal plant that is a source of nutritional supplements rich in salvianolic acid B and a potential experimental system for the exploration of salvianolic acid B biosynthesis in the Labiatae. Here, we report a high-quality chromosome-scale genome assembly of S. bowleyana covering 462.44 Mb, with a scaffold N50 value of 57.96 Mb and 44,044 annotated protein-coding genes. Evolutionary analysis revealed an estimated divergence time between S. bowleyana and its close relative S. miltiorrhiza of ~3.94 million years. We also observed evidence of a whole-genome duplication in the S. bowleyana genome. Transcriptome analysis showed that SbPAL1 (PHENYLALANINE AMMONIA-LYASE1) is highly expressed in roots relative to stem and leaves, paralleling the location of salvianolic acid B accumulation. The laccase gene family in S. bowleyana outnumbered their counterparts in both S. miltiorrhiza and Arabidopsis thaliana, suggesting that the gene family has undergone expansion in S. bowleyana. Several laccase genes were also highly expressed in roots, where their encoded proteins may catalyze the oxidative reaction from rosmarinic acid to salvianolic acid B. These findings provide an invaluable genomic resource for understanding salvianolic acid B biosynthesis and its regulation, and will be useful for exploring the evolution of the Labiatae.

Full-length transcriptome sequences by a combination of sequencing platforms applied to isoflavonoid and triterpenoid saponin biosynthesis of Astragalus mongholicus Bunge

BACKGROUND

Astragalus mongholicus Bunge is an important medicinal plant used in traditional Chinese medicine. It is rich in isoflavonoids and triterpenoid saponins. Although these active constituents of A. mongholicus have been discovered for a long time, the genetic basis of isoflavonoid and triterpenoid saponin biosynthesis in this plant is virtually unknown because of the lack of a reference genome. Here, we used a combination of next-generation sequencing (NGS) and single-molecule real-time (SMRT) sequencing to identify genes involved in the biosynthetic pathway of secondary metabolites in A. mongholicus.

RESULTS

In this study, NGS, SMRT sequencing, and targeted compound analysis were combined to investigate the association between isoflavonoid and triterpenoid saponin content, and specific gene expression in the root, stem, and leaves of A. mongholicus. Overall, 643,812 CCS reads were generated, yielding 121,107 non-redundant transcript isoforms with an N50 value of 2124 bp. Based on these highly accurate transcripts, 104,756 (86.50%) transcripts were successfully annotated by any of the seven databases (NR, NT, Swissprot, KEGG, KOG, Pfam and GO). Levels of four isoflavonoids and four astragalosides (triterpenoid saponins) were determined. Forty-four differentially expressed genes (DEGs) involved in isoflavonoid biosynthesis and 44 DEGs from 16 gene families that encode enzymes involved in triterpenoid saponin biosynthesis were identified. Transcription factors (TFs) associated with isoflavonoid and triterpenoid saponin biosynthesis, including 72 MYBs, 53 bHLHs, 64 AP2-EREBPs, and 11 bZIPs, were also identified. The above transcripts showed different expression trends in different plant organs.

CONCLUSIONS

This study provides important genetic information on the A. mongholicus genes that are essential for isoflavonoid and triterpenoid saponin biosynthesis, and provides a basis for developing the medicinal value of this plant.

Genome-wide identification and functional characterization of natural antisense transcripts in Salvia miltiorrhiza

Salvia miltiorrhiza is one of the most widely used traditional medicines. Natural antisense transcripts (NATs) are a class of long noncoding RNAs that can regulate gene expression. Here, we identified 812 NATs, including 168 cis-NATs and 644 trans-NATs from twelve root, flower, and leaf samples of S. miltiorrhiza using RNA-seq. The expression profiles for 41 of 50 NATs and their sense transcripts (STs) obtained from RNA-Seq were validated using qRT-PCR. The expression profiles of 17 NATs positively correlated with their STs. GO and KEGG pathway analyses mapped the STs for cis-NATs to pathways for biosynthesis of secondary metabolites. We characterized four NATs in detail, including NAT0001, NAT0002, NAT0004, and NAT00023. Their STs are kaurene synthase-like 1 and the homologs of UDP-glucose flavonoid 3-O-glucosyltransferase 6, UDP-glycosyltransferase 90A1, and beta-glucosidase 40, respectively. The first gene is involved in the biosynthesis of bioactive tanshinones, the next two are involved in anthocyanin biosynthesis, whereas the last is involved in phenylpropanoid biosynthesis. Besides, we found seven STs that are potential targets of miRNAs. And we found two miRNAs including miR156a and miR7208, might originate from NATs, NAT0112 and NAT0086. The results suggest that S. miltiorrhiza NATs might interact with STs, produce miRNAs, and be regulated by miRNAs. They potentially play significant regulatory roles in the biosynthesis of bioactive compounds.

Full-length transcriptome analysis and identification of genes involved in asarinin and aristolochic acid biosynthesis in medicinal plant Asarum sieboldii

Asarum sieboldii, a well-known traditional Chinese medicinal herb, is used for curing inflammation and ache. It contains both the bioactive ingredient asarinin and the toxic compound aristolochic acid. To address further breeding demand, genes involved in the biosynthetic pathways of asarinin and aristolochic acid should be explored. Therefore, the full-length transcriptome of A. sieboldii was sequenced using PacBio Iso-Seq to determine the candidate transcripts that encode the biosynthetic enzymes of asarinin and aristolochic acid. In this study, 63 023 full-length transcripts were generated with an average length of 1371 bp from roots, stems, and leaves, of which 49 593 transcripts (78.69%) were annotated against public databases. Furthermore, 555 alternative splicing (AS) events, 10 869 long noncoding RNAs (lncRNAs) as well as their 11 291 target genes, and 17 909 simple sequence repeats (SSRs) were identified. The data also revealed 97 candidate transcripts related to asarinin metabolism, of which six novel genes that encoded enzymes involved in asarinin biosynthesis were initially reported. In addition, 56 transcripts related to aristolochic acid biosynthesis were also identified, including CYP81B. In summary, these transcriptome data provide a useful resource to study gene function and genetic engineering in A. sieboldii.

Transcriptomic analysis reveals the GRAS family genes respond to gibberellin in Salvia miltiorrhiza hairy roots

BACKGROUND

Salvia miltiorrhiza is one of the most important traditional Chinese medicinal plants with high medicinal value. Gibberellins are growth-promoting phytohormones that regulate numerous growth and developmental processes in plants. However, their role on the secondary metabolism regulation has not been investigated.

RESULTS

In this study, we found that gibberellic acid (GA) can promote hairy roots growth and increase the contents of tanshinones and phenolic acids. Transcriptomic sequencing revealed that many genes involved in the secondary metabolism pathway were the GA-responsive. After further analysis of GA signaling pathway genes, which their expression profiles have significantly changed, it was found that the GRAS transcription factor family had a significant response to GA. We identified 35 SmGRAS genes in S. miltiorrhiza, which can be divided into 10 subfamilies. Thereafter, members of the same subfamily showed similar conserved motifs and gene structures, suggesting possible conserved functions.

CONCLUSIONS

Most SmGRAS genes were significantly responsive to GA, indicating that they may play an important role in the GA signaling pathway, also participating in the GA regulation of root growth and secondary metabolism in S. miltiorrhiza.

smi-miR396b targeted SmGRFs, SmHDT1, and SmMYB37/4 synergistically regulates cell growth and active ingredient accumulation in Salvia miltiorrhiza hairy roots

MIR396b had been cloned and overexpressed in Salvia miltiorrhiza hairy roots. MiR396b targets SmGRFs, SmHDT1, and SmMYB37/4 to regulate cell growth and secondary metabolism in S. miltiorrhiza hairy roots. Danshen (Salvia miltiorrhiza Bunge) is a valuable medicinal herb with two kinds of clinically used natural products, salvianolic acids and tanshinones. miR396 is a conserved microRNA and plays extensive roles in plants. However, it is still unclear how miR396 works in S. miltiorrhiza. In this study, an smi-MIR396b has been cloned from S. miltiorrhiza. Overexpression of miR396b in danshen hairy roots inhibited hairy root growth, reduced salvianolic acid concentration, but enhanced tanshinone accumulation, resulting in the biomass and total salvianolic acids respectively reduced to 55.5 and 72.1% of the control and total tanshinones increased up to 1.91-fold of the control. Applied degradome sequencing, 5'RLM-RACE, and qRT-PCR, 13 targets for miR396b were identified including seven conserved SmGRF1-7 and six novel ones. Comparative transcriptomics and microRNomics analysis together with qRT-PCR results confirmed that miR396b targets SmGRFs, SmHDT1, and SmMYB37/4 to mediate the phytohormone, especially gibberellin signaling pathways and consequentially resulted in the phenotype variation of miR396b-OE hairy roots. Furthermore, miR396b could be activated by methyl jasmonate, abscisic acid, gibberellin, salt, and drought stresses. The findings in this study indicated that smi-miR396b acts as an upstream and central regulator in cell growth and the biosynthesis of tanshinones and salvianolic acids, shedding light on the coordinated regulation of plant growth and biosynthesis of active ingredients in S. miltiorrhiza.

SmMYC2b Enhances Tanshinone Accumulation in Salvia miltiorrhiza by Activating Pathway Genes and Promoting Lateral Root Development

Salvia miltiorrhiza Bunge (Lamiaceae) is an economically important medicinal plant as well as an emerging model plant. Our previous studies indicate that SmMYC2b is a positive transcription factor that can affect the biosynthesis of phenolic acids and tanshinones in S. miltiorrhiza. Moreover, MYC2s are well known to induce the development of lateral roots. As tanshinones are mainly distributed in the periderm, the promotion of lateral root development probably leads to increased accumulation of tanshinones. In this paper, we firstly discovered that SmMYC2b played a dual regulatory role in effectively enhancing the tanshinone accumulation by activating tanshinone biosynthetic pathway and promoting lateral root development. The expression levels of the previously studied pathway genes SmCPS1, SmKSL1, SmCYP76AH1, SmCYP76AH3, and SmCYP76AK1 dramatically increased. In addition, SmMYC2b was proved to exhibit a similar function as other homologs in promoting lateral root development, which increased the tanshinone produced tissue and further enhanced the biosynthesis of tanshinones. RNA-seq assays revealed that SmMYC2b-regulated genes comprised 30.6% (1,901 of 6,210) of JA-responsive genes, confirming that SmMYC2b played a crucial role in transcriptional regulation of JA-regulated genes. Overall, we concluded that SmMYC2b could enhance tanshinone accumulation by activating the tanshinone biosynthetic pathway and promoting lateral root development. Our study provides an effective approach to enhance the production of desired tanshinones and enriches our knowledge of the related regulatory network.

SmGRAS1 and SmGRAS2 Regulate the Biosynthesis of Tanshinones and Phenolic Acids in Salvia miltiorrhiza

Salvia miltiorrhiza is one of the most widely used traditional Chinese medicinal plants because of its excellent performance in treating heart diseases. Tanshinones and phenolic acids are two important classes of effective metabolites, and their biosynthesis has attracted widespread interest. Here, we functionally characterized SmGRAS1 and SmGRAS2, two GRAS family transcription factors from S. miltiorrhiza. SmGRAS1/2 were highly expressed in the root periderm, where tanshinones mainly accumulated in S. miltiorrhiza. Overexpression of SmGRAS1/2 upregulated tanshinones accumulation and downregulated GA, phenolic acids contents, and root biomass. However, antisense expression of SmGRAS1/2 reduced the tanshinones accumulation and increased the GA, phenolic acids contents, and root biomass. The expression patterns of biosynthesis genes were consistent with the changes in compounds accumulation. GA treatment increased tanshinones, phenolic acids, and GA contents in the overexpression lines, and restored the root growth inhibited by overexpressing SmGRAS1/2. Subsequently, yeast one-hybrid, dual-luciferase, and electrophoretic mobility shift assays (EMSA) showed SmGRAS1 promoted tanshinones biosynthesis by directly binding to the GARE motif in the SmKSL1 promoter and activating its expression. Yeast two-hybrid assays showed SmGRAS1 interacted physically with SmGRAS2. Taken together, the results revealed that SmGRAS1/2 acted as repressors in root growth and phenolic acids biosynthesis but as positive regulators in tanshinones biosynthesis. Overall, our findings revealed the potential value of SmGRAS1/2 in genetically engineering changes in secondary metabolism.

The application of combined 1H NMR-based metabolomics and transcriptomics techniques to explore phenolic acid biosynthesis in Salvia miltiorrhiza Bunge

Salvia miltiorrhiza Bunge is a traditional Chinese medicine, and its water-soluble phenolic acid active compounds have very important medicinal value; however, the synthesis pathways of the main active ingredients remain unknown. Here, we employed nuclear magnetic resonance (NMR)-based metabolomics and transcriptomics techniques to study the biosynthesis mechanism of salvianolic acids. High-performance liquid chromatography (HPLC) combined with NMR showed an improvement over traditional techniques, and 54 metabolites were detected. The results of the multivariate statistical analysis showed that salvianolic acid B (SAB), rosmarinic acid (RA), caffeic acid, succinate, and citrate were among the multiple compounds that were increased in the methyl jasmonate (MeJA)-elicited group; the levels of sucrose, fructose, glutamine, and tyrosine were decreased. Combined with the differentially expressed genes (DEGs) found by transcriptome sequencing, we speculate that the synthesis of RA after MeJA treatment mostly occurred through caffeic acid and bypassed 4-hydroxyphenyllactic acid. This provides useful information for the study of salvianolic acids synthesis.

Hybrid-Transcriptome Sequencing and Associated Metabolite Analysis Reveal Putative Genes Involved in Flower Color Difference in Rose Mutants

Gene mutation is a common phenomenon in nature that often leads to phenotype differences, such as the variations in flower color that frequently occur in roses. With the aim of revealing the genomic information and inner mechanisms, the differences in the levels of both transcription and secondary metabolism between a pair of natural rose mutants were investigated by using hybrid RNA-sequencing and metabolite analysis. Metabolite analysis showed that glycosylated derivatives of pelargonidin, e.g., pelargonidin 3,5 diglucoside and pelargonidin 3-glucoside, which were not detected in white flowers (Rosa 'Whilte Mrago Koster'), constituted the major pigments in pink flowers. Conversely, the flavonol contents of petal, such as kaempferol-3-glucoside, quercetin 3-glucoside, and rutin, were higher in white flowers. Hybrid RNA-sequencing obtained a total of 107,280 full-length transcripts in rose petal which were annotated in major databases. Differentially expressed gene (DEG) analysis showed that the expression of genes involved in the flavonoid biosynthesis pathway was significantly different, e.g., CHS, FLS, DFR, LDOX, which was verified by qRT-PCR during flowering. Additionally, two MYB transcription factors were found and named RmMYBAN2 and RmMYBPA1, and their expression patterns during flowering were also analyzed. These findings indicate that these genes may be involved in the flower color difference in the rose mutants, and competition between anthocyanin and flavonol biosynthesis is a primary cause of flower color variation, with its regulation reflected by transcriptional and secondary metabolite levels.

Identification of putative genes for polyphenol biosynthesis in olive fruits and leaves using full-length transcriptome sequencing

Olive (Olea europaea) is a rich source of valuable bioactive polyphenols, which has attracted widespread interest. In this study, we combined targeted metabolome, Pacbio ISOseq transcriptome, and Illumina RNA-seq transcriptome to investigate the association between polyphenols and gene expression in the developing olive fruits and leaves. A total of 12 main polyphenols were measured, and 122 transcripts of 17 gene families, 101 transcripts of 9 gene families, and 106 transcripts of 6 gene families that encode for enzymes involved in flavonoid, oleuropein, and hydroxytyrosol biosynthesis were separately identified. Additionally, 232 alternative splicing events of 18 genes related to polyphenol synthesis were analyzed. This is the first time that the third generations of full-length transcriptome technology were used to study the gene expression pattern of olive fruits and leaves. The results of transcriptome combined with targeted metabolome can help us better understand the polyphenol biosynthesis pathways in the olive.

Genome-Wide Identification and Characterization of Salvia miltiorrhiza Laccases Reveal Potential Targets for Salvianolic Acid B Biosynthesis

Laccases are widely distributed in plant kingdom catalyzing the polymerization of lignin monolignols. Rosmarinic acid (RA) has a lignin monolignol-like structure and is converted into salvianolic acid B (SAB), which is a representatively effective hydrophilic compound of a well-known medicinal plant Salvia miltiorrhiza and also the final compound of phenolic acids metabolism pathway in the plant. But the roles of laccases in the biosynthesis of SAB are poorly understood. This work systematically characterizes S. miltiorrhiza laccase (SmLAC) gene family and identifies the SAB-specific candidates. Totally, 29 laccase candidates (SmLAC1-SmLAC29) are found to contain three signature Cu-oxidase domains. They present relatively low sequence identity and diverse intron-exon patterns. The phylogenetic clustering of laccases from S. miltiorrhiza and other ten plants indicates that the 29 SmLACs can be divided into seven groups, revealing potential distinct functions. Existence of diverse cis regulatory elements in the SmLACs promoters suggests putative interactions with transcription factors. Seven SmLACs are found to be potential targets of miR397. Putative glycosylation sites and phosphorylation sites are identified in SmLAC amino acid sequences. Moreover, the expression profile of SmLACs in different organs and tissues deciphers that 5 SmLACs (SmLAC7/8/20/27/28) are expressed preferentially in roots, adding the evidence that they may be involved in the phenylpropanoid metabolic pathway. Besides, silencing of SmLAC7, SmLAC20 and SmLAC28, and overexpression of SmLAC7 and SmLAC20 in the hairy roots of S. miltiorrhiza result in diversification of SAB, signifying that SmLAC7 and SmLAC20 take roles in SAB biosynthesis. The results of this study lay a foundation for further elucidation of laccase functions in S. miltiorrhiza, and add to the knowledge for SAB biosynthesis in S. miltiorrhiza.

Apollo: Democratizing genome annotation

Genome annotation is the process of identifying the location and function of a genome's encoded features. Improving the biological accuracy of annotation is a complex and iterative process requiring researchers to review and incorporate multiple sources of information such as transcriptome alignments, predictive models based on sequence profiles, and comparisons to features found in related organisms. Because rapidly decreasing costs are enabling an ever-growing number of scientists to incorporate sequencing as a routine laboratory technique, there is widespread demand for tools that can assist in the deliberative analytical review of genomic information. To this end, we present Apollo, an open source software package that enables researchers to efficiently inspect and refine the precise structure and role of genomic features in a graphical browser-based platform. Some of Apollo’s newer user interface features include support for real-time collaboration, allowing distributed users to simultaneously edit the same encoded features while also instantly seeing the updates made by other researchers on the same region in a manner similar to Google Docs. Its technical architecture enables Apollo to be integrated into multiple existing genomic analysis pipelines and heterogeneous laboratory workflow platforms. Finally, we consider the implications that Apollo and related applications may have on how the results of genome research are published and made accessible.

Biosynthesis, Chemistry, and Pharmacology of Polyphenols from Chinese Salvia Species: A Review

Salvia species find widespread application in food and pharmaceutical products owing to their large polyphenol content. The main polyphenols in Chinese Salvia species are phenolic acids and flavonoids, which exhibit anti-oxygenation, anti-ischemia-reperfusion injury, anti-thrombosis, anti-tumour, and other therapeutic effects. However, there are few peer-reviewed studies on polyphenols in Chinese Salvia species, especially flavonoids. This review is a systematic, comprehensive collation of available information on the biosynthesis, chemistry, and pharmacology of Chinese Salvia species. We believe that our study makes a significant contribution to the literature because this review provides a detailed literary resource on the currently available information on various polyphenolic components of Chinese Salvia species, including their bioactivities and structures. In addition, the study provides information that would encourage further investigation of this plant material as a natural resource with potential for a broad range of applications in various industries, such as the food and pharmaceutical industries.

Expanding Alternative Splicing Identification by Integrating Multiple Sources of Transcription Data in Tomato

Tomato (Solanum lycopersicum) is an important vegetable and fruit crop. Its genome was completely sequenced and there are also a large amount of available expressed sequence tags (ESTs) and short reads generated by RNA sequencing (RNA-seq) technologies. Mapping transcripts including mRNA sequences, ESTs, and RNA-seq reads to the genome allows identifying pre-mRNA alternative splicing (AS), a post-transcriptional process generating two or more RNA isoforms from one pre-mRNA transcript. We comprehensively analyzed the AS landscape in tomato by integrating genome mapping information of all available mRNA and ESTs with mapping information of RNA-seq reads which were collected from 27 published projects. A total of 369,911 AS events were identified from 34,419 genomic loci involving 161,913 transcripts. Within the basic AS events, intron retention is the prevalent type (18.9%), followed by alternative acceptor site (12.9%) and alternative donor site (7.3%), with exon skipping as the least type (6.0%). Complex AS types having two or more basic event accounted for 54.9% of total AS events. Within 35,768 annotated protein-coding gene models, 23,233 gene models were found having pre-mRNAs generating AS isoform transcripts. Thus the estimated AS rate was 65.0% in tomato. The list of identified AS genes with their corresponding transcript isoforms serves as a catalog for further detailed examination of gene functions in tomato biology. The post-transcriptional information is also expected to be useful in improving the predicted gene models in tomato. The sequence and annotation information can be accessed at plant alternative splicing database (http://proteomics.ysu.edu/altsplice).

Molecular Regulation of Catalpol and Acteoside Accumulation in Radial Striation and non-Radial Striation of Rehmannia glutinosa Tuberous Root

Rehmannia glutinosa L., a perennial plant of Scrophulariaceae, is one of the most commonly used herbs in traditional Chinese medicine (TCM) that have been widely cultivated in China. However, to date, the biosynthetic pathway of its two quality-control components, catalpol and acteoside, are only partially elucidated and the mechanism for their tissue-specific accumulation remains unknown. To facilitate the basic understanding of the key genes and transcriptional regulators involved in the biosynthesis of catalpol and acteoside, transcriptome sequencing of radial striation (RS) and non-radial striation (nRS) from four R. glutinosa cultivars was performed. A total of 715,158,202 (~107.27 Gb) high quality reads obtained using paired-end Illumina sequencing were de novo assembled into 150,405 transcripts. Functional annotation with multiple public databases identified 155 and 223 unigenes involved in catalpol and acteoside biosynthesis, together with 325 UGTs, and important transcription factor (TF) families. Comparative analysis of the transcriptomes identified 362 unigenes, found to be differentially expressed in all RS vs. nRS comparisons, with 143 upregulated unigenes, including those encoding enzymes of the catalpol and acteoside biosynthetic pathway, such as geranyl diphosphate synthase (RgGPPS), geraniol 8-hydroxylase (RgG10H), and phenylalanine ammonia-lyase (RgPAL). Other differentially expressed unigenes predicted to be related to catalpol and acteoside biosynthesis fall into UDP-dependent glycosyltransferases (UGTs), as well as transcription factors. In addition, 16 differentially expressed genes were selectively confirmed by real-time PCR. In conclusion, a large unigene dataset of R. glutinosa generated in the current study will serve as a resource for the identification of potential candidate genes for investigation of the tuberous root development and biosynthesis of active components.

Full-length transcriptome sequences and the identification of putative genes for flavonoid biosynthesis in safflower

Background

The flower of the safflower (Carthamus tinctorius L.) has been widely used in traditional Chinese medicine for the ability to improve cerebral blood flow. Flavonoids are the primary bioactive components in safflower, and their biosynthesis has attracted widespread interest. Previous studies mostly used second-generation sequencing platforms to survey the putative flavonoid biosynthesis genes. For a better understanding of transcription data and the putative genes involved in flavonoid biosynthesis in safflower, we carry our study.

Results

High-quality RNA was extracted from six types of safflower tissue. The RNAs of different tissues were mixed equally and used for multiple size-fractionated libraries (1–2, 2–3 and 3-6 k) library construction. Five cells were carried (2 cells for 1–2 and for 2-3 k libraries and 1 cell for 3-6 k libraries). 10.43Gb clean data and 38,302 de-redundant sequences were captured. 44 unique isoforms were annotated as encoding enzymes involved in flavonoid biosynthesis. The full length flavonoid genes were characterized and their evolutional relationship and expressional pattern were analyzed. They can be divided into eight families, with a large differences in the tissue expression. The temporal expressions under MeJA treatment were also measured, 9 genes are significantly up-regulated and 2 genes are significantly down-regulated. The genes involved in flavonoid synthesis in safflower were predicted in our study. Besides, the SSR and lncRNA are also analyzed in our study.

Conclusions

Full-length transcriptome sequences were used in our study. The genes involved in flavonoid synthesis in safflower were predicted in our study. Combined the determination of flavonoids, CtC4H2, CtCHS3, CtCHI3, CtF3H3, CtF3H1 are mainly participated in MeJA promoting the synthesis of flavonoids. Our results also provide a valuable resource for further study on safflower.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4946-9) contains supplementary material, which is available to authorized users.

Genomic survey of bZIP transcription factor genes related to tanshinone biosynthesis in Salvia miltiorrhiza

Tanshinones are a class of bioactive components in the traditional Chinese medicine Salvia miltiorrhiza, and their biosynthesis and regulation have been widely studied. Current studies show that basic leucine zipper (bZIP) proteins regulate plant secondary metabolism, growth and developmental processes. However, the bZIP transcription factors involved in tanshinone biosynthesis are unknown. Here, we conducted the first genome-wide survey of the bZIP gene family and analyzed the phylogeny, gene structure, additional conserved motifs and alternative splicing events in S. miltiorrhiza. A total of 70 SmbZIP transcription factors were identified and categorized into 11 subgroups based on their phylogenetic relationships with those in Arabidopsis. Moreover, seventeen SmbZIP genes underwent alternative splicing events. According to the transcriptomic data, the SmbZIP genes that were highly expressed in the Danshen root and periderm were selected. Based on the prediction of bZIP binding sites in the promoters and the co-expression analysis and co-induction patterns in response to Ag⁺ treatment via quantitative real-time polymerase chain reaction (qRT-PCR), we concluded that SmbZIP7 and SmbZIP20 potentially participate in the regulation of tanshinone biosynthesis. These results provide a foundation for further functional characterization of the candidate SmbZIP genes, which have the potential to increase tanshinone production.

Overexpression of SmMYC2 Increases the Production of Phenolic Acids in Salvia miltiorrhiza

MYC2 is a core transcription factor in the plant response to jasmonates. It also functions in secondary metabolism and various processes for growth and development. However, the knowledge about its role in Salvia miltiorrhiza is still very limited. We determined that the biosynthesis of salvianolic acid B (Sal B) was strongly induced in 2-month-old transgenic plants that over-expressed SmMYC2. In the roots of transgenic line 12 that over-expressed SmMYC2 (OEM-12), the Sal B concentration was as high as 5.95 ± 0.07 mg g-1, a level that was 1.88-fold higher than that in control plants that had been transformed with an empty vector. Neither tanshinone IIA nor cryptotanshinone was detected by high-performance liquid chromatography in any of the genotypes. Global transcriptomic analysis using RNA sequencing revealed that most enzyme-encoding genes for the phenylpropanoid biosynthesis pathway were up-regulated in the overexpression lines. Furthermore, both the phenylalanine and tyrosine biosynthesis pathways were activated in those transgenics. Our data demonstrate that overexpression of SmMYC2 promotes the production of phenolic acids by simultaneously activating both primary and secondary pathways for metabolism in S. miltiorrhiza.

Researches on Transcriptome Sequencing in the Study of Traditional Chinese Medicine

Due to its incomparable advantages, the application of transcriptome sequencing in the study of traditional Chinese medicine attracts more and more attention of researchers, which greatly promote the development of traditional Chinese medicine. In this paper, the applications of transcriptome sequencing in traditional Chinese medicine were summarized by reviewing recent related papers.

Overexpression of AtEDT1 promotes root elongation and affects medicinal secondary metabolite biosynthesis in roots of transgenic Salvia miltiorrhiza

Medicinal secondary metabolites (salvianolic acids and tanshinones) are valuable natural bioactive compounds in Salvia miltiorrhiza and have widespread applications. Improvement of medicinal secondary metabolite accumulation through biotechnology is necessary and urgent to satisfy their increasing demand. Herein, it was demonstrated that the overexpression of the transcription factor Arabidopsis thaliana-enhanced drought tolerance 1 (AtEDT1) could affect medicinal secondary metabolite accumulation. In this study, we observed that the transgenic lines significantly conferred drought tolerance phenotype. Meanwhile, we found that the overexpression of AtEDT1 promoted root elongation in S. miltiorrhiza. Interestingly, we also found that the overexpression of AtEDT1 determined the accumulation of salvianolic acids, such as rosmarinic acid, lithospermic acid, salvianolic acid B, and total salvianolic acids due to the induction of the expression levels of salvianolic acid biosynthetic genes. Conversely, S. miltiorrhiza plants overexpressing the AtEDT1 transgene showed a decrease in tanshinone synthesis. Our results demonstrated that the overexpression of AtEDT1 significantly increased the accumulation of salvianolic acids in S. miltiorrhiza. Further studies are required to better elucidate the functional role of AtEDT1 in the regulation of phytochemical compound synthesis.

Transcriptomic analyses reveal biosynthetic genes related to rosmarinic acid in Dracocephalum tanguticum

Dracocephalum tanguticum Maxim, a Lamiaceae species endemic to the Qinghai-Tibetan Plateau and adjacent regions, is an important ornamental, medicinal and aromatic herb. In this study, a comprehensive transcriptome of 18 libraries from six organs namely, roots, stems, leaves, sepals, flowers and seeds of D. tanguticum were generated. More than 100 Gb of sequence data were obtained and assembled de novo into 187,447 transcripts, including 151,463 unigenes, among which the six organs shared 17.7% (26,841). In addition, all unigenes were assigned to 362 pathways, in which 'biosynthesis of secondary metabolites' is the second enriched pathway. Furthermore, rosmarinic acid (RA) is one of the multifunctional phenolic bioactive compounds produced in some Lamiaceae species. The six organs of D. tanguticum were confirmed to produce RA. A total of 22 predicted biosynthetic genes related to RA from the transcriptome were further isolated. Two of these genes were identified as candidates by evaluating the correlation coefficient between the RA contents and the expression of the predicted biosynthetic genes in the six organs. The new sequence information will improve the knowledge of D. tanguticum, as well as provide a reference tool for future studies of biosynthetic genes related to RA in this species.

Transcriptome Analysis of Salicylic Acid Treatment in Rehmannia glutinosa Hairy Roots Using RNA-seq Technique for Identification of Genes Involved in Acteoside Biosynthesis

Rehmannia glutinosa is a common bulk medicinal material that has been widely used in China due to its active ingredients. Acteoside, one of the ingredients, has antioxidant, antinephritic, anti-inflammatory, hepatoprotective, immunomodulatory, and neuroprotective effects, is usually selected as a quality-control component for R. glutinosa herb in the Chinese Pharmacopeia. The acteoside biosynthesis pathway in R. glutinosa has not yet been clearly established. Herein, we describe the establishment of a genetic transformation system for R. glutinosa mediated by Agrobacterium rhizogenes. We screened the optimal elicitors that markedly increased acteoside accumulation in R. glutinosa hairy roots. We found that acteoside accumulation dramatically increased with the addition of salicylic acid (SA); the optimal SA dose was 25 μmol/L for hairy roots. RNA-seq was applied to analyze the transcriptomic changes in hairy roots treated with SA for 24 h in comparison with an untreated control. A total of 3,716, 4,018, and 2,715 differentially expressed transcripts (DETs) were identified in 0 h-vs.-12 h, 0 h-vs.-24 h, and 12 h-vs.-24 h libraries, respectively. KEGG pathway-based analysis revealed that 127 DETs were enriched in "phenylpropanoid biosynthesis." Of 219 putative unigenes involved in acteoside biosynthesis, 54 were found to be up-regulated at at least one of the time points after SA treatment. Selected candidate genes were analyzed by quantitative real-time PCR (qRT-PCR) in hairy roots with SA, methyl jasmonate (MeJA), AgNO₃ (Ag⁺), and putrescine (Put) treatment. All genes investigated were up-regulated by SA treatment, and most candidate genes were weakly increased by MeJA to some degree. Furthermore, transcription abundance of eight candidate genes in tuberous roots of the high-acteoside-content (HA) cultivar QH were higher than those of the low-acteoside-content (LA) cultivar Wen 85-5. These results will pave the way for understanding the molecular basis of acteoside biosynthesis in R. glutinosa, and can serve as a basis for future validation studies.

Transcriptome-Guided Mining of Genes Involved in Crocin Biosynthesis

Gardenia jasminoides is used in traditional Chinese medicine and has drawn attention as a rich source of crocin, a compound with reported activity against various cancers, depression and cardiovascular disease. However, genetic information on the crocin biosynthetic pathway of G. jasminoides is scarce. In this study, we performed a transcriptome analysis of the leaves, green fruits, and red fruits of G. jasminoides to identify and predict the genes that encode key enzymes responsible for crocin production, compared with Crocus sativus. Twenty-seven putative pathway genes were specifically expressed in the fruits, consistent with the distribution of crocin in G. jasminoides. Twenty-four of these genes were reported for the first time, and a novel CCD4a gene was predicted that encodes carotenoid cleavage dioxygenase leading to crocin synthesis, in contrast to CCD2 of C. sativus. In addition, 6 other candidate genes (ALDH12, ALDH14, UGT94U1, UGT86D1, UGT71H4, and UGT85K18) were predicted to be involved in crocin biosynthesis following phylogenetic analysis and different gene expression profiles. Identifying the genes that encode key enzymes should help elucidate the crocin biosynthesis pathway.

Transcriptome Profiling Using Single-Molecule Direct RNA Sequencing Approach for In-depth Understanding of Genes in Secondary Metabolism Pathways of Camellia sinensis

Characteristic secondary metabolites, including flavonoids, theanine and caffeine, are important components of Camellia sinensis, and their biosynthesis has attracted widespread interest. Previous studies on the biosynthesis of these major secondary metabolites using next-generation sequencing technologies limited the accurately prediction of full-length (FL) splice isoforms. Herein, we applied single-molecule sequencing to pooled tea plant tissues, to provide a more complete transcriptome of C. sinensis. Moreover, we identified 94 FL transcripts and four alternative splicing events for enzyme-coding genes involved in the biosynthesis of flavonoids, theanine and caffeine. According to the comparison between long-read isoforms and assemble transcripts, we improved the quality and accuracy of genes sequenced by short-read next-generation sequencing technology. The resulting FL transcripts, together with the improved assembled transcripts and identified alternative splicing events, enhance our understanding of genes involved in the biosynthesis of characteristic secondary metabolites in C. sinensis.