De novo transcriptome assembly and characterization of nine tissues of Lonicera japonica to identify potential candidate genes involved in chlorogenic acid, luteolosides, and secoiridoid biosynthesis pathways

Regulation of gene expression in the secondary metabolic synthesis pathway of Lonicera japonica Flos by yeast polysaccharides

Lonicera japonica Flos is a valuable herb in the Lonicerae family. While transcriptomic studies on L. japonica have focused on different tissues (stems, leaves, flowers) or flowering stages, few have investigated the molecular mechanisms underlying chemical composition synthesis influenced by exogenous factors, such as foliar fertilization. Moreover, most transcriptomic studies on L. Japonica have been conducted on chlorogenic acid and luteoloside, and the molecular synthesis mechanism of the overall chemical composition has not been analyzed. Methods: We conducted a single-factor, four-level foliar fertilization experiment using yeast polysaccharides. Different yeast polysaccharides concentrations were sprayed on L. japonica for six consecutive days with dynamic sampling. High-performance liquid chromatography determined the active ingredients in each group. The two groups exhibiting the most significant differences were selected for transcriptomic analysis to identify key synthetic genes responsible for L. japonica's active ingredients. Key results: Principal component analysis conducted on samples collected on September 8 revealed significant differences in the active ingredient amounts between the 0.1 g/L yeast polysaccharides treatment group and the control group. Transcriptome sequencing analysis identified 218 significantly differentially expressed genes, including 60 upregulated and 158 downregulated genes. Twelve differential genes involved in the chemical components synthesis pathway of L. japonica under yeast polysaccharides treatment were identified: PAL1, PAL2, PAL3, 4CL1, 4CL, CHS1, CHS2, CHS, CHI1, CHI2, F3H, and SOH. Conclusions: This study contributes to the theoretical understanding of essential synthetic genes associated with L. japonica's active ingredients. It offers data support for further gene exploration and sheds light on the molecular mechanisms underlying L. japonica quality formation. These findings hold significant implications for enhancing the content of secondary metabolites of L. japonica.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01482-1.

Genome-wide analysis of the MADS-box gene family in Lonicera japonica and a proposed floral organ identity model

BACKGROUND

Lonicera japonica Thunb. is widely used in traditional Chinese medicine. Medicinal L. japonica mainly consists of dried flower buds and partially opened flowers, thus flowers are an important quality indicator. MADS-box genes encode transcription factors that regulate flower development. However, little is known about these genes in L. japonica.

RESULTS

In this study, 48 MADS-box genes were identified in L. japonica, including 20 Type-I genes (8 Mα, 2 Mβ, and 10 Mγ) and 28 Type-II genes (26 MIKCc and 2 MIKC*). The Type-I and Type-II genes differed significantly in gene structure, conserved domains, protein structure, chromosomal distribution, phylogenesis, and expression pattern. Type-I genes had a simpler gene structure, lacked the K domain, had low protein structure conservation, were tandemly distributed on the chromosomes, had more frequent lineage-specific duplications, and were expressed at low levels. In contrast, Type-II genes had a more complex gene structure; contained conserved M, I, K, and C domains; had highly conserved protein structure; and were expressed at high levels throughout the flowering period. Eleven floral homeotic MADS-box genes that are orthologous to the proposed Arabidopsis ABCDE model of floral organ identity determination, were identified in L. japonica. By integrating expression pattern and protein interaction data for these genes, we developed a possible model for floral organ identity determination.

CONCLUSION

This study genome-widely identified and characterized the MADS-box gene family in L. japonica. Eleven floral homeotic MADS-box genes were identified and a possible model for floral organ identity determination was also developed. This study contributes to our understanding of the MADS-box gene family and its possible involvement in floral organ development in L. japonica.

Volatile Profiling and Transcriptome Sequencing Provide Insights into the Biosynthesis of α-Pinene and β-Pinene in Liquidambar formosana Hance Leaves

Liquidambar formosana Hance is a pinene-rich deciduous plant species in the Altingiaceae family that is used as a medicinal plant in China. However, the regulatory mechanisms underlying α-pinene and β-pinene biosynthesis in L. formosana leaves remain unknown. Here, a joint analysis of the volatile compounds and transcriptomes of L. formosana leaves was performed to comprehensively explore the terpene synthase (TPS) that may participate in α-pinene and β-pinene biosynthesis. Headspace solid-phase microextraction (HS-SPME) and gas chromatography-mass spectrometry (GC-MS) jointly detected volatile L. formosana leaves. Trees with high and low levels of both α-pinene and β-pinene were defined as the H group and L group, respectively. RNA sequencing data revealed that DXR (1-deoxy-D-xylulose-5-phosphate reductoisomerase), HDS [(E)-4-hydroxy-3-methylbut-2-eny-l-diphosphate synthase], and TPS may be the major regulators of monoterpenoid biosynthesis. We identified three TPSs (LfTPS1, LfTPS2, and LfTPS3), which are highly homologous to α-pinene and β-pinene synthases of other species in phylogenetic analysis. Four TPS genes (LfTPS1, LfTPS2, LfTPS4, LfTPS5) may be critically involved in the biosynthesis and regulation of α-pinene and β-pinene in L. formosana. Bioinformatic and transcriptomic results were verified using quantitative real-time PCR. We identified LfTPS1, LfTPS2 as candidate genes for α-pinene and β-pinene biosynthesis that significantly improve the yield of beneficial terpenoids.

Medicinal Plants: Guests and Hosts in the Heterologous Expression of High-Value Products

Medicinal plants play an important dual role in the context of the heterologous expression of high-value pharmaceutical products. On the one hand, the classical biochemical and modern omics approaches allowed for the discovery of various genes encoding biosynthetic pathways in medicinal plants. Recombinant DNA technology enabled introducing these genes and regulatory elements into host organisms and enhancing the heterologous production of the corresponding secondary metabolites. On the other hand, the transient expression of foreign DNA in plants facilitated the production of numerous proteins of pharmaceutical importance. This review summarizes several success stories of the engineering of plant metabolic pathways in heterologous hosts. Likewise, a few examples of recombinant protein expression in plants for therapeutic purposes are also highlighted. Therefore, the importance of medicinal plants has grown immensely as sources for valuable products of low and high molecular weight. The next step ahead for bioengineering is to achieve more success stories of industrial-scale production of secondary plant metabolites in microbial systems and to fully exploit plant cell factories' commercial potential for recombinant proteins.

Comparative transcriptome and tissue-specific expression analysis of genes reveal tissue-cultured plants as an alternative source for phenylethanoids and phenylpropanoids in Rhodiola imbricata (Edgew.)

Rhodiola imbricata (Crassulaceae) is a traditional trans-Himalayan endangered medicinal herb with immense therapeutic applications. Over the years, over-exploitation, un-managed harvesting, and lack of captive cultivation procedures persuaded threat to its wild habitat. Plant tissue culture and RNA-Seq-based molecular bioprospection of key regulatory genes aid the understanding of molecular dynamics involved in specialized metabolites (phenylethanoids and phenylpropanoids) biosynthesis and its sustainable production. Hence, comparative transcriptomic analysis was performed using leaf and root tissues from the wild and tissue-cultured plants, revealing tissue-specific production of salidroside and rosavin. The transcriptome profiling resulted in 345 million high-quality reads yielding 92,380 unique transcripts with an N₅₀ of 1260 bp. Tissue-specific gene expression analysis revealed that both phenylethanoids and phenylpropanoids biosynthesis are predominantly associated with the shikimate pathway. In addition to RNA-Seq data, the downstream biosynthesis pathways genes viz., phospho-2-dehydro-3-deoxyheptonate aldolase (DAHPS), 3-dehydroquinate synthase (DHQS), shikimate kinase (SK), chorismate mutase (CM), arogenate dehydrogenase (TYRAAT), aromatic-L-amino-acid decarboxylase (TDC), phenylalanine ammonia-lyase (PAL), 4-coumarate-CoA ligase (4-CL), cinnamoyl-CoA reductase (CCR), and cinnamyl alcohol dehydrogenase (CAD) showed higher expression pattern in wild plant tissues compared to tissue-cultured plants. The transcript fold expression determined by RT-qPCR results followed similar patterns as those observed in RNA-seq and targeted metabolite profiling data. Salidroside and rosavin content in wild plants exhibited 2.40 fold and 1.77 fold increase accumulation compared to the tissue-cultured plant. The present investigation explained the tissue and condition-specific significant differences between the expression of proposed biosynthetic pathway genes and salidroside and rosavin content. Additionally, NAC, bHLH, and ARF were the most abundant transcription factor families found in the transcriptomic analysis of R. imbricata. The generated transcriptome dataset provides a valuable gene(s)/transcription factors hub that can be used for the sustainable production of salidroside and rosavin in R. imbricata under tissue culture conditions.

Differences in gene expression and endophytic bacterial diversity in Atractylodes macrocephala Koidz. rhizomes from different growth years

Atractylodes macrocephala Koidz. (AMK) is widely used owing to its pharmacological activity in traditional Chinese medicine (TCM). Here, we aimed to characterize the differentially expressed genes (DEGs) of one- and three-year growth (OYG and TYG) rhizomes of AMK combined with the endophytic bacterial diversity analysis using high-throughput RNA-sequencing. 114,572 unigenes were annotated in six public databases. 3570 DEGs revealed a clear difference, of which 936 and 2634 genes were up- and down-regulated, respectively. The results of KEGG pathway analysis indicated that DEGs corresponding to the terpenoid synthesis gene were downregulated in TYG rhizomes. 414,424 sequences corresponding to the 16S rRNA gene were divided into 1267 operational taxonomic units (OTUs). Moreover, the diversity of endophytic bacteria changed with species in OYG (773) and TYG (1201) rhizomes at OTU level, and Proteobacteria, Actinobacteria, and Bacteroidetes were the dominant phyla. Comparison of species differences among different growth years revealed that some species were significantly different, such as Actinomycetes, Variovorax, Cloacibacterium, etc. Interestingly, the decrease in the function-related metabolism of terpenoids and polyketides was found to be correlated the low expression of terpene synthesis genes in TYG rhizomes assessed using PICRUSt2. These data provide a scientific basis for elucidating the mechanism underlying metabolite accumulation and endophytic bacterial diversity in relation to the growth years in AMK.

Comparative transcriptome analysis reveals variations of bioactive constituents in Lonicera japonica flowers under salt stress

Lonicera japonica flowers (LJF) is a traditional Chinese medicine packed with phenols constituents and widely used in the treatments of various diseases throughout the world. However, there is still very little known on how LJF identifies and resists salt stress. Here in, we systematically investigated the effect of salt on the phenotypic, metabolite, and transcriptomic in LJF. During long term stress (35 days), 1055 differential expression genes (DEGs) involved in the biosynthesis of secondary metabolites were screened through transcriptome analysis, among which the candidate genes and pathways involved in phenols biosynthesis were highlighted; and performed by phylogenetic tree analysis and multiple nucleotide sequence alignment. Ninety compounds were identified and their relative levels were compared between the control and stressed groups based on the LC-MS analysis, Putative biosynthesis networks of phenolic acid and flavonoid were con-structed with structural DEGs. Strikingly, the expression patterns of structural DEGs were mostly consistent with the variations of phenols under salt stress. Notably, the upregulation of UDP-glycosyl transferases under salt stress indicated post-modification of glycosyl transferases may participate in downstream flavonoids synthesis. This study reveals the relationships of the gene regulation and the phenols biosynthesis in LJF under salt stress, paving the way for the use of gene-specific expression to improve the yield of biocomponent.

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.

LjaFGD: Lonicera japonica functional genomics database

Lonicera japonica Thunb., a traditional Chinese herb, has been used for treating human diseases for thousands of years. Recently, the genome of L. japonica has been decoded, providing valuable information for research into gene function. However, no comprehensive database for gene functional analysis and mining is available for L. japonica. We therefore constructed LjaFGD (www.gzybioinformatics.cn/LjaFGD and bioinformatics.cau.edu.cn/LjaFGD), a database for analyzing and comparing gene function in L. japonica. We constructed a gene co-expression network based on 77 RNA-seq samples, and then annotated genes of L. japonica by alignment against protein sequences from public databases. We also introduced several tools for gene functional analysis, including Blast, motif analysis, gene set enrichment analysis, heatmap analysis, and JBrowse. Our co-expression network revealed that MYB and WRKY transcription factor family genes were co-expressed with genes encoding key enzymes in the biosynthesis of chlorogenic acid and luteolin in L. japonica. We used flavonol synthase 1 (LjFLS1) as an example to show the reliability and applicability of our database. LjaFGD and its various associated tools will provide researchers with an accessible platform for retrieving functional information on L. japonica genes to further biological discovery.

Transcriptomic and metabolomic analyses provide insights into the biosynthesis of chlorogenic acids in Lonicera macranthoides Hand.-Mazz

Lonicera macranthoides Hand.-Mazz (L. macranthoides) is a medicinal herb that is widely distributed in South China. The developmental stage and corolla dehiscence of the flower are the important factors affecting the quality of medicinal ingredients. However, neither the regulatory mechanism controlling chlorogenic acids biosynthesis in L. macranthoides nor the molecular basis of effect of corolla dehiscence on the quality of medicinal materials is fully understood. In this study, metabolomics and transcriptomics were used to analyze the metabolic and transcriptional differences of two different cultivars closed bud type (Bt), and flowering type (Ft), as well as the effect of jasmonic acid methyl ester (MeJA) on chlorogenic acids (CGAs) biosynthesis. In total, large number of differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs) were filtered among three lines of samples. Gene metabolite correlation analyses revealed a ‘core set’ of 30 genes and 54 genes that were strongly correlated with CGAs biosynthesis and regulating the flowering, respectively. Quantitative real-time polymerase chain reaction results proved the alterations in the expression levels of genes encoding the pathways involved in CGAs biosynthesis. The ion abundances of CGAs were most significantly increased, while some of the CGAs derived and Caffeoyl-CoA-derived substances showed the most largely reduced abundances in the closed bud type (Bt) compared to the flowering type (Ft). MeJA may leads to the activation of downstream genes in CGAs biosynthesis pathway. Overall, there were significant differences in the transcriptional and metabolic levels of CGAs biosynthesis pathway in flower buds of different flowering cultivars. The redirection of metabolic flux may contribute to increased accumulation of CGAs. However, whether MeJA and flowering have direct effects on the accumulation of CGAs needs further studied. These researches effectively expanded the functional genomic library and provide new insights into CGAs biosynthesis in L. macranthoides.

Multiomics-based characterization of specialized metabolites biosynthesis in Cornus Officinalis

Cornus officinalis, an important traditional medicinal plant, is used as major constituents of tonics, analgesics, and diuretics. While several studies have focused on its characteristic bioactive compounds, little is known on their biosynthesis. In this study, we performed LC-QTOF-MS-based metabolome and RNA-seq-based transcriptome profiling for seven tissues of C. officinalis. Untargeted metabolome analysis assigned chemical identities to 1,215 metabolites and showed tissue-specific accumulation for specialized metabolites with medicinal properties. De novo transcriptome assembly established for C. officinalis showed 96% of transcriptome completeness. Co-expression analysis identified candidate genes involved in the biosynthesis of iridoids, triterpenoids, and gallotannins, the major group of bioactive metabolites identified in C. officinalis. Integrative omics analysis identified 45 cytochrome P450s genes correlated with iridoids accumulation in C. officinalis. Network-based integration of genes assigned to iridoids biosynthesis pathways with these candidate CYPs further identified seven promising CYPs associated with iridoids’ metabolism. This study provides a valuable resource for further investigation of specialized metabolites’ biosynthesis in C. officinalis.

Toxicity as prime selection criterion among SARS-active herbal medications

We present here a new selection criterion for prioritizing research on efficacious drugs for the fight against COVID-19: the relative toxicity versus safety of herbal medications, which were effective against SARS in the 2002/2003 epidemic. We rank these medicines according to their toxicity versus safety as basis for preferential rapid research on their potential in the treatment of COVID-19. The data demonstrate that from toxicological information nothing speaks against immediate investigation on, followed by rapid implementation of Lonicera japonica, Morus alba, Forsythia suspensa, and Codonopsis spec. for treatment of COVID-19 patients. Glycyrrhiza spec. and Panax ginseng are ranked in second priority and ephedrine-free Herba Ephedrae extract in third priority (followed by several drugs in lower preferences). Rapid research on their efficacy in the therapy - as well as safety under the specific circumstances of COVID-19 - followed by equally rapid implementation will provide substantial advantages to Public Health including immediate availability, enlargement of medicinal possibilities, in cases where other means are not successful (non-responders), not tolerated (sensitive individuals) or just not available (as is presently the case) and thus minimize sufferings and save lives. Moreover, their moderate costs and convenient oral application are especially advantageous for underprivileged populations in developing countries.

Enhanced Biosynthesis of Chlorogenic Acid and Its Derivatives in Methyl-Jasmonate-Treated Gardenia jasminoides Cells: A Study on Metabolic and Transcriptional Responses of Cells

Chlorogenic acid and its derivatives (CQAs) are considered as important bioactive secondary metabolites in Gardenia jasminoides Ellis (G. jasminoides). However, few studies have investigated the biosynthesis and regulation of CQAs in G. jasminoides. In this study, methyl jasmonate (MeJA) was used to enhance CQAs accumulation in cultured G. jasminoides cells. Moreover, the possible molecular mechanism of MeJA-mediated accumulation of CQAs is also explored. To this end, time-course transcriptional profiles of G. jasminoides cells responding to MeJA were used to investigate the mechanism from different aspects, including jasmonate (JAs) biosynthesis, signal transduction, biosynthesis of precursor, CQAs biosynthesis, transporters, and transcription factors (TFs). A total of 57,069 unigenes were assembled from the clean reads, in which 80.7% unigenes were successfully annotated. Furthermore, comparative transcriptomic results indicated that differentially expressed genes (DEGs) were mainly involved in JAs biosynthesis and signal transduction (25 DEGs), biosynthesis of precursor for CQAs (18 DEGs), CQAs biosynthesis (19 DEGs), and transporters (29 DEGs). Most of these DEGs showed continuously upregulated expressions over time, which might activate the jasmonic acid (JA) signal transduction network, boost precursor supply, and ultimately stimulate CQAs biosynthesis. Additionally, various TFs from different TF families also responded to MeJA elicitation. Interestingly, 38 DEGs from different subgroups of the MYB family might display positive or negative regulations on phenylpropanoids, especially on CQAs biosynthesis. Conclusively, our results provide insight into the possible molecular mechanism of regulation on CQAs biosynthesis, which led to a high CQAs yield in the G. jasminoides cells under MeJA treatment.

Metabolomics Reveals Distinct Metabolites between Lonicera japonica and Lonicera macranthoides Based on GC-MS

Lonicera japonica Thunb. (LJ) and Lonicera macranthoides Hand. -Mazz. (LM) have been widely used in Chinese medicine for thousands of years. Although the morphological characteristics of LJ and LM are quite similar, there are significant distinctions of medicinal ingredients (mainly the secondary metabolites) and clinical indications between them. However, the in-depth differences of primary metabolites have not thoroughly been studied yet. Therefore, gas chromatography-mass spectrometry- (GC-MS-) based metabolomics method combined with chemometric methods were performed to analyze the distinction in this study. The results showed that LJ and LM were obviously classified into two groups. 10 metabolites were obtained as biomarkers on account of their p values, pcorr values, and differing variable importance in projection (VIP) values. Metabolic pathway analysis showed that the galactose metabolism and starch and sucrose metabolism gathered as potential pathways caused these extraordinary differences of primary metabolites between LJ and LM. Further, we found that the differences of main medicinal ingredients between LJ and LM could be interpreted from these metabolites according to the analysis of mainly related pathways. The metabolites involved in the starch and sucrose metabolism presented upregulated in LJ, while almost all metabolites in the galactose metabolism, the TCA cycle, and the phenolic acid part of phenylpropanoid metabolism were downregulated in LJ. Therefore, the energy stored in the starch and sucrose metabolism may be saved to produce flavonoid, which could be the reason that the level of flavonoid of phenylpropanoid metabolism is higher in LJ compared to LM. Consequently, this study presented an effective tool for quality evaluation of LJ and LM and laid a foundation for further studies of the metabolic mechanisms and high-quality manufacturing of them.

The honeysuckle genome provides insight into the molecular mechanism of carotenoid metabolism underlying dynamic flower coloration

Lonicera japonica is a widespread member of the Caprifoliaceae (honeysuckle) family utilized in traditional medical practices. This twining vine honeysuckle also is a much-sought ornamental, in part due to its dynamic flower coloration, which changes from white to gold during development. The molecular mechanism underlying dynamic flower coloration in L. japonica was elucidated by integrating whole genome sequencing, transcriptomic analysis and biochemical assays. Here, we report a chromosome-level genome assembly of L. japonica, comprising nine pseudochromosomes with a total size of 843.2 Mb. We also provide evidence for a whole-genome duplication event in the lineage leading to L. japonica, which occurred after its divergence from Dipsacales and Asterales. Moreover, gene expression analysis not only revealed correlated expression of the relevant biosynthetic genes with carotenoid accumulation, but also suggested a role for carotenoid degradation in L. japonica's dynamic flower coloration. The variation of flower color is consistent with not only the observed carotenoid accumulation pattern, but also with the release of volatile apocarotenoids that presumably serve as pollinator attractants. Beyond novel insights into the evolution and dynamics of flower coloration, the high-quality L. japonica genome sequence also provides a foundation for molecular breeding to improve desired characteristics.

Comparative transcriptomics analysis revealing flower trichome development during flower development in two Lonicera japonica Thunb. cultivars using RNA-seq

BACKGROUND

Lonicera japonica Thunb. (L. japonica) has the functions of clearing away heat and detoxifying, broad-spectrum antibacterial and anti-virus, etc. More than 70% of anti-inflammatory and cold Chinese patent medicines contain L. japonica. Trichomes comprise specialized multicellular structures that have the capacity to synthesize and secrete secondary metabolites and protect plants from biotic and abiotic stresses. The extraction of trichome secretions has great commercial value. However, little is known about the trichome formation mechanism in L. japonica. Therefore, the study of trichome development between different varieties provides a basis for selecting suitable planting resources.

RESULTS

Here, we present a genome-wide comparative transcriptome analysis between two L. japonica cultivars, toward the identification of biological processes and functional gene activities that occur during flowering stage trichome development. In this study, the density and average lengths of flower trichomes were at their highest during three-green periods (S2). Using the Illumina RNA-Seq method, we obtained 134,304 unigenes, 33,733 of which were differentially expressed. In an analysis of 40 differentially expressed unigenes (DEGs) involved in trichome development, 29 of these were transcription factors. The DEGs analysis of plant hormone signal transduction indicated that plant growth and development may be independent of gibberellin (GA) and cytokinine (CTK) signaling pathways, and plant stress may be independent of jasmonic acid (JA) and ethylene (ET) signaling pathways. We screened several genes involved in the floral biosynthesis of odors, tastes, colors, and plant hormones, and proposed biosynthetic pathways for sesquiterpenoid, triterpenoid, monoterpenoid, flavonoid, and plant hormones. Furthermore, 82 DEGs were assigned to cell cycles and 2616 were predicted as plant resistance genes (PRGs).

CONCLUSIONS

This study provides a comprehensive characterization of the expression profiles of flower development during the seven developmental stages of L. japonica, thereby offering valuable insights into the molecular networks that underly flower development in L. japonica.

Cross-Species Annotation of Expressed Genes and Detection of Different Functional Gene Modules Between 10 Cold- and 10 Hot-Propertied Chinese Herbal Medicines

According to the traditional Chinese medicine (TCM) system, Chinese herbal medicines (HMs) can be divided into four categories: hot, warm, cold, and cool. A cool nature usually is categorized as a cold nature, and a warm nature is classified as a hot nature. However, the detectable characteristics of the gene expression profile associated with the cold and hot properties have not been studied. To address this question, a strategy for the cross-species annotation of conserved genes was established in the present study by using transcriptome data of 20 HMs with cold and hot properties. Functional enrichment analysis was performed on group-specific expressed genes inferred from the functional genome of the reference species (i.e., Arabidopsis). Results showed that metabolic pathways relevant to chrysoeriol, luteolin, paniculatin, and wogonin were enriched for cold-specific genes, and pathways of inositol, heptadecane, lauric acid, octanoic acid, hexadecanoic acid, and pentadecanoic acid were enriched for hot-specific genes. Six functional modules were identified in the HMs with the cold property: nucleotide biosynthetic process, peptidy-L-cysteine S-palmitoylation, lipid modification, base-excision repair, dipeptide transport, and response to endoplasmic reticulum stress. For the hot HMs, another six functional modules were identified: embryonic meristem development, embryonic pattern specification, axis specification, regulation of RNA polymerase II transcriptional preinitiation complex assembly, mitochondrial RNA modification, and cell redox homeostasis. The research provided a new insight into HMs’ cold and hot properties from the perspective of the gene expression profile of plants.

Comparative transcriptomic analyses of chlorogenic acid and luteolosides biosynthesis pathways at different flowering stages of diploid and tetraploid Lonicera japonica

The Flos Lonicerae Japonicae (FLJ), Lonicera japonica Thunb, belonging to the Caprifoliaceae family, is an economically important plant that is highly utilized in traditional Chinese medicine as well as in Japanese medicine. The flowers of these plants are rich in chlorogenic acid (CGA) and luteoloside. Our previous study revealed that tetraploid L. japonica has higher fresh/dry weight, phenolic acids and flavonoids contents than those of diploid plants. However, why tetraploid L. japonica can yield higher CGA and luteolosides than that in diploid and what is the difference in the molecular regulatory mechanism of these pathways between diploid and tetraploids remained unclear. Therefore, in the present study, we performed comprehensive transcriptome analyses of different flowering stages of diploid and tetraploid L. japonica. The CGA content of tetraploid was found higher than that of diploid at all the growth stages. While the luteoloside content of diploid was higher than that of tetraploid at S4 and S6 growth stages. We obtained a high-quality transcriptome assembly (N50 = 2,055 bp; Average length = 1,331 bp) compared to earlier studies. Differential expression analysis revealed that several important genes involving in plant hormone signal transduction, carbon metabolism, starch/sucrose metabolism and plant-pathogen interaction were upregulated in tetraploid compared with the diploid L. japonica, reflecting the higher adaptability and resistance of tetraploid species. Furthermore, by associating the phenotypic data and gene expression profiles, we were able to characterize the potential molecular regulatory mechanism of important biosynthetic pathways at different flowering stages. Overall, our work provides a foundation for further research on these important secondary metabolite pathways and their implications in traditional Chinese/Japanese medicine.

De novo transcriptomic analysis of light-induced flavonoid pathway, transcription factors in the flower buds of Lonicera japonica

KEY MESSAGE

Transcriptomic analysis of the relationship between gene expression patterns and flavonoid contents in the flower buds of Lonicera japonica under light-induced conditions, especially the flavonoid pathway genes and transcription factors.

ABSTRACT

Flos Lonicerae Japonicae (FLJ), the flower buds of Lonicera japonica Thunb., has been used to treat some human diseases including severe respiratory syndromes and hand-foot-and-mouth diseases owing to its putative antibacterial, and antiviral effects. Luteoloside is a flavonoid that is used by the Chinese Pharmacopoeia to evaluate the quality of FLJ. Light is an important environmental factor that affects flavonoid biosynthesis in the flower buds of L. japonica. However, how light triggers increases in flavonoid production remains unclear. To enhance our understanding of the mechanism involved in light-regulated flavonoid biosynthesis, we sequenced the transcriptomes of L. japonica exposed to three different light conditions: 100% light intensity (CK), 50% light intensity (LI50), and 25% light intensity (LI25) using an Illumina HiSeq 4000 System. A total of 77,297 unigenes with an average length of 809 bp were obtained. Among them, 43,334 unigenes (56.06%) could be matched to at least one biomolecular database. Additionally, 4188, 1545 and 1023 differentially expressed genes (DEGs) were identified by comparative transcriptomics LI25-vs-CK, LI50-vs-CK, and LI25-vs-LI50, respectively. Of note, genes known to be involved in flavonoid biosynthesis, such as 4-coumarate coenzyme A ligase (4CL), and chalcone synthase (CHS) were up-regulated. In addition, a total of 1649 transcription factors (TFs) were identified and divided into 58 TF families; 98 TFs exhibited highly dynamic changes in response to light intensity. Quantitative real-time PCR (qRT-PCR) was used to test the expression profiles of the RNA sequencing (RNA-Seq) data. This study offers insight into how transcriptional expression pattern is influenced by light in the flower buds of L. japonica, and will enhance the understanding of molecular mechanisms of flavonoid biosynthesis in response to light in L. japonica.

Ameliorating effects of exogenous calcium on the photosynthetic physiology of honeysuckle (Lonicera japonica) under salt stress

Calcium (Ca2+) plays pivotal roles in modulating plant growth, development and stress responses. This work was conducted to study the effects of 20 mM calcium on the biomass, malondialdehyde content, chlorophyll content, ion ratio, chlorophyll a fluorescence and gas-exchange parameters, gene expression of annual honeysuckle under 50, 100 and 200 mM NaCl. At the end of treatment, Na+ concentration was increased with the mounting salinity, but a higher ratio of K+/Na2+, Ca2+/Na+, Mg2+/Na+ were obtained after calcium addition. Salinity exerted an adverse effect on the dry weights and chlorophyll content, whereas CaCl2 played a positive role. Consistent with biomass reduction, the photosynthetic rate and stomatal conductance declined in leaves of honeysuckle exposed to elevated salinity. However, the extent of reduction was much less under CaCl2 combination treatments than one caused by NaCl treatments. Exogenous calcium also protects the photochemical activity of PSII by protecting reaction centre from inactivation and maintaining electron transport from QA- to QB-. Further, exogenous calcium promoted the overexpression of LHCB coding gene Cab and Rubisco large subunit coding gene rbcL under short-term stress. In conclusion, exogenous calcium was effective in improving the salt tolerance of honeysuckle in the photosynthetic base, thereby improving the growth of plants.

Transcriptional regulation of Lonicera japonica Thunb. during flower development as revealed by comprehensive analysis of transcription factors

BACKGROUND

Lonicera japonica Thunb. flower has been used for the treatment of various diseases for a long time and attracted many studies on its potential effects. Transcription factors (TFs) regulate extensive biological processes during plant development. As the restricted reports of L. japonica on TFs, our work was carried out to better understand the TFs' regulatory roles under different developmental stages in L. japonica.

RESULTS

In this study, 1316 TFs belonging to 52 families were identified from the transcriptomic data, and corresponding expression profiles during the L. japonica flower development were comprehensively analyzed. 917 (69.68%) TFs were differentially expressed. TFs in bHLH, ERF, MYB, bZIP, and NAC families exhibited obviously altered expression during flower growth. Based on the analysis of differentially expressed TFs (DETFs), TFs in MYB, WRKY, NAC and LSD families that involved in phenylpropanoids biosynthesis, senescence processes and antioxidant activity were detected. The expression of MYB114 exhibited a positive correlation with the contents of luteoloside; Positive correlation was observed among the expression of MYC12, chalcone synthase (CHS) and flavonol synthase (FLS), while negative correlation was observed between the expression of MYB44 and the synthases; The expression of LSD1 was highly correlated with the expression of SOD and the total antioxidant capacity, while the expression of LOL1 and LOL2 exhibited a negative correlation with them; Many TFs in NAC and WRKY families may be potentially involved in the senescence process regulated by hormones and reactive oxygen species (ROS). The expression of NAC19, NAC29, and NAC53 exhibited a positive correlation with the contents of ABA and H₂O₂, while the expression of WRKY53, WRKY54, and WRKY70 exhibited a negative correlation with the contents of JA, SA and ABA.

CONCLUSIONS

Our study provided a comprehensive characterization of the expression profiles of TFs during the developmental stages of L. japonica. In addition, we detected the key TFs that may play significant roles in controlling active components biosynthesis, antioxidant activity and flower senescence in L. japonica, thereby providing valuable insights into the molecular networks underlying L. japonica flower development.

Computational Methods for Mapping, Assembly and Quantification for Coding and Non-coding Transcripts

The measurement of gene expression has long provided significant insight into biological functions. The development of high-throughput short-read sequencing technology has revealed transcriptional complexity at an unprecedented scale, and informed almost all areas of biology. However, as researchers have sought to gather more insights from the data, these new technologies have also increased the computational analysis burden. In this review, we describe typical computational pipelines for RNA-Seq analysis and discuss their strengths and weaknesses for the assembly, quantification and analysis of coding and non-coding RNAs. We also discuss the assembly of transposable elements into transcripts, and the difficulty these repetitive elements pose. In summary, RNA-Seq is a powerful technology that is likely to remain a key asset in the biologist's toolkit.

Evolutionary Analysis of Plastid Genomes of Seven Lonicera L. Species: Implications for Sequence Divergence and Phylogenetic Relationships

Plant plastomes play crucial roles in species evolution and phylogenetic reconstruction studies due to being maternally inherited and due to the moderate evolutionary rate of genomes. However, patterns of sequence divergence and molecular evolution of the plastid genomes in the horticulturally- and economically-important Lonicera L. species are poorly understood. In this study, we collected the complete plastomes of seven Lonicera species and determined the various repeat sequence variations and protein sequence evolution by comparative genomic analysis. A total of 498 repeats were identified in plastid genomes, which included tandem (130), dispersed (277), and palindromic (91) types of repeat variations. Simple sequence repeat (SSR) elements analysis indicated the enriched SSRs in seven genomes to be mononucleotides, followed by tetra-nucleotides, dinucleotides, tri-nucleotides, hex-nucleotides, and penta-nucleotides. We identified 18 divergence hotspot regions (rps15, rps16, rps18, rpl23, psaJ, infA, ycf1, trnN-GUU-ndhF, rpoC2-rpoC1, rbcL-psaI, trnI-CAU-ycf2, psbZ-trnG-UCC, trnK-UUU-rps16, infA-rps8, rpl14-rpl16, trnV-GAC-rrn16, trnL-UAA intron, and rps12-clpP) that could be used as the potential molecular genetic markers for the further study of population genetics and phylogenetic evolution of Lonicera species. We found that a large number of repeat sequences were distributed in the divergence hotspots of plastid genomes. Interestingly, 16 genes were determined under positive selection, which included four genes for the subunits of ribosome proteins (rps7, rpl2, rpl16, and rpl22), three genes for the subunits of photosystem proteins (psaJ, psbC, and ycf4), three NADH oxidoreductase genes (ndhB, ndhH, and ndhK), two subunits of ATP genes (atpA and atpB), and four other genes (infA, rbcL, ycf1, and ycf2). Phylogenetic analysis based on the whole plastome demonstrated that the seven Lonicera species form a highly-supported monophyletic clade. The availability of these plastid genomes provides important genetic information for further species identification and biological research on Lonicera.

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.

Chemical and genomic diversity of six Lonicera species occurring in Korea

Lonicera spp. (Caprifoliaceae) are important not only as a common medicinal herb in East Asia but also as one of the most problematic invasive species in North America. In the present study, we performed a systemic analysis of genomic and chemical diversity among six Lonicera species occurring in Korea, L. japonica, L. maackii, L. insularis, L. sachalinensis, L. praeflorens, and L. vesicaria, using chloroplast DNA whole genome shotgun (WGS) sequencing and LC-MS analyses. The phylogenetic and phylochemical relationships did not coincide with each other, but partial consistency could be found among them. InDel-based cDNA marker for authentication was developed based on the genome sequences. Flavonoids, iridoids, and organic acids were identified in the LC-MS analyses, and their inter-species distribution and localization were also revealed.

Trends in herbgenomics

From Shen Nong's Herbal Classic (Shennong Bencao Jing) to the Compendium of Materia Medica (Bencao Gangmu) and the first scientific Nobel Prize for the mainland of China, each milestone in the historical process of the development of traditional Chinese medicine (TCM) involves screening, testing and integrating. After thousands of years of inheritance and development, herbgenomics (bencaogenomics) has bridged the gap between TCM and international advanced omics studies, promoting the application of frontier technologies in TCM. It is a discipline that uncovers the genetic information and regulatory networks of herbs to clarify their molecular mechanism in the prevention and treatment of human diseases. The main theoretical system includes genomics, functional genomics, proteomics, transcriptomics, metabolomics, epigenomics, metagenomics, synthetic biology, pharmacogenomics of TCM, and bioinformatics, among other fields. Herbgenomics is mainly applicable to the study of medicinal model plants, genomic-assisted breeding, herbal synthetic biology, protection and utilization of gene resources, TCM quality evaluation and control, and TCM drug development. Such studies will accelerate the application of cutting-edge technologies, revitalize herbal research, and strongly promote the development and modernization of TCM.

De novo characterization of the Baphicacanthus cusia(Nees) Bremek transcriptome and analysis of candidate genes involved in indican biosynthesis and metabolism

Baphicacanthus cusia (Nees) Bremek is an herb widely used for the clinical treatment of colds, fever, and influenza in Traditional Chinese Medicine. The roots, stems and leaves can be used as natural medicine, in which indigo and indirubin are two main active ingredients. In this study, quantification of indigo, indirubin, indican and adenosine among various tissues of B. cusia was conducted using HPLC-DAD. Leaves have significantly higher contents than stems and roots (380.66, 315.15, 20,978.26, 4323.15 μg/g in leaves, 306.36, 71.71, 3,056.78, 139.45 μg/g in stems, and 9.31, 7.82, 170.45, 197.48 μg/g in roots, respectively). De novo transcriptome sequencing of B. cusia was performed for the first time. The sequencing yielded 137,216,248, 122,837,394 and 140,240,688 clean reads from leaves, stems and roots respectively, which were assembled into 51,381 unique sequences. A total of 33,317 unigenes could be annotated using the databases of Nr, Swiss-Prot, KEGG and KOG. These analyses provided a detailed view of the enzymes involved in indican backbone biosynthesis, such as cytochrome P450, UDP-glycosyltransferase, glucosidase and tryptophan synthase. Analysis results showed that tryptophan synthase was the candidate gene involved in the tissue-specific biosynthesis of indican. We also detected sixteen types of simple sequence repeats in RNA-Seq data for use in future molecular mark assisted breeding studies. The results will be helpful in further analysis of B. cusia functional genomics, especially in increasing biosynthesis of indican through biotechnological approaches and metabolic regulation.

Transcriptome Analysis Reveals Molecular Signatures of Luteoloside Accumulation in Senescing Leaves of Lonicera macranthoides

Lonicera macranthoides is an important medicinal plant widely used in traditional Chinese medicine. Luteoloside is a critical bioactive compound in L. macranthoides. To date, the molecular mechanisms underlying luteoloside biosynthesis are still largely unknown. In this work, high performance liquid chromatography (HPLC) was employed to determine the luteoloside contents in leaves, stems, and flowers at different developmental stages. Results showed that senescing leaves can accumulate large amounts of luteoloside, extremely higher than that in young and semi-lignified leaves and other tissues. RNA-Seq analysis identified that twenty-four differentially expressed unigenes (DEGs) associated with luteoloside biosynthesis were significantly up-regulated in senescing leaves, which are positively correlated with luteoloside accumulation. These DEGs include phenylalanine ammonia lyase 2, cinnamate 4-hydroxylase 2, thirteen 4-coumarate-CoA ligases, chalcone synthase 2, six flavonoid 3'-monooxygenase (F3'H) and two flavone 7-O-β-glucosyltransferase (UFGT) genes. Further analysis demonstrated that two F3'Hs (CL11828.Contig1 and CL11828.Contig2) and two UFGTs (Unigene2918 and Unigene97915) might play vital roles in luteoloside generation. Furthermore, several transcription factors (TFs) related to flavonoid biosynthesis including MYB, bHLH and WD40, were differentially expressed during leaf senescence. Among these TFs, MYB12, MYB75, bHLH113 and TTG1 were considered to be key factors involved in the regulation of luteoloside biosynthesis. These findings provide insights for elucidating the molecular signatures of luteoloside accumulation in L. macranthoides.

[Development of Plant Metabolomics and Medicinal Plant Genomics]

 A variety of chemicals produced by plants, often referred to as 'phytochemicals', have been used as medicines, food, fuels and industrial raw materials. Recent advances in the study of genomics and metabolomics in plant science have accelerated our understanding of the mechanisms, regulation and evolution of the biosynthesis of specialized plant products. We can now address such questions as how the metabolomic diversity of plants is originated at the levels of genome, and how we should apply this knowledge to drug discovery, industry and agriculture. Our research group has focused on metabolomics-based functional genomics over the last 15 years and we have developed a new research area called 'Phytochemical Genomics'. In this review, the development of a research platform for plant metabolomics is discussed first, to provide a better understanding of the chemical diversity of plants. Then, representative applications of metabolomics to functional genomics in a model plant, Arabidopsis thaliana, are described. The extension of integrated multi-omics analyses to non-model specialized plants, e.g., medicinal plants, is presented, including the identification of novel genes, metabolites and networks for the biosynthesis of flavonoids, alkaloids, sulfur-containing metabolites and terpenoids. Further, functional genomics studies on a variety of medicinal plants is presented. I also discuss future trends in pharmacognosy and related sciences.

Proteomic Contributions to Medicinal Plant Research: From Plant Metabolism to Pharmacological Action

Herbal medicine is a clinical practice of utilizing medicinal plant derivatives for therapeutic purposes. It has an enduring history worldwide and plays a significant role in the fight against various diseases. Herbal drug combinations often exhibit synergistic therapeutic action compared with single-constituent dosage, and can also enhance the cytotoxicity induced by chemotherapeutic drugs. To explore the mechanism underlying the pharmacological action of herbs, proteomic approaches have been applied to the physiology of medicinal plants and its effects on animals. This review article focuses on the existing proteomics-based medicinal plant research and discusses the following topics: (i) plant metabolic pathways that synthesize an array of bioactive compounds; (ii) pharmacological action of plants tested using in vivo and in vitro studies; and (iii) the application of proteomic approaches to indigenous plants with scarce sequence information. The accumulation of proteomic information in a biological or medicinal context may help in formulating the effective use of medicinal plants.

Characterization of the Lycium barbarum fruit transcriptome and development of EST-SSR markers

Lycium barbarum, commonly known as goji, is important in Chinese herbal medicine and its fruit is a very important agricultural and biological product. However, the molecular mechanism of formation of its fruit and associated medicinal and nutritional components is unexplored. Moreover, this species lacks SSR markers due to lack of genomic and transcriptomic information. In this study, a total of 139,333 unigenes with average length of 1049 bp and N50 of 1579 bp are obtained by trinity assembly from Illumina sequencing reads. A total of 92,498 (66.38%) unigenes showed similarities in at least one database including Nr (46.15%), Nt (56.56%), KO (15.56%), Swiss-prot (33.34%), Pfam (33.43%), GO (33.62%) and KOG/COG (17.55%). Genes in flavonoid and taurine biosynthesis pathways were found and validated by RT-qPCR. A total of 50,093 EST-SSRs were identified from 38,922 unigenes, and 22,537 EST-SSR primer pairs were designed. Four hundred pairs of SSR markers were randomly selected to validate assembly quality, of which 352 (88%) were successful in PCR amplification of genomic DNA from 11 Lycium accessions and 210 produced polymorphisms. The polymorphic loci showed that the genetic similarity of the 11 Lycium accessions ranged from 0.50 to 0.99 and the accessions could be divided into 4 groups. These results will facilitate investigations of the molecular mechanism of formation of L. barbarum fruit and associated medicinal and nutritional components, and will be of value to novel gene discovery and functional genomic studies. The EST-SSR markers will be useful for genetic diversity evaluation, genetic mapping and marker-assisted breeding.

The unexpected depths of genome-skimming data: A case study examining Goodeniaceae floral symmetry genes

PREMISE OF THE STUDY

The use of genome skimming allows systematists to quickly generate large data sets, particularly of sequences in high abundance (e.g., plastomes); however, researchers may be overlooking data in low abundance that could be used for phylogenetic or evo-devo studies. Here, we present a bioinformatics approach that explores the low-abundance portion of genome-skimming next-generation sequencing libraries in the fan-flowered Goodeniaceae.

METHODS

Twenty-four previously constructed Goodeniaceae genome-skimming Illumina libraries were examined for their utility in mining low-copy nuclear genes involved in floral symmetry, specifically the CYCLOIDEA (CYC)-like genes. De novo assemblies were generated using multiple assemblers, and BLAST searches were performed for CYC1, CYC2, and CYC3 genes.

RESULTS

Overall Trinity, SOAPdenovo-Trans, and SOAPdenovo implementing lower k-mer values uncovered the most data, although no assembler consistently outperformed the others. Using SOAPdenovo-Trans across all 24 data sets, we recovered four CYC-like gene groups (CYC1, CYC2, CYC3A, and CYC3B) from a majority of the species. Alignments of the fragments included the entire coding sequence as well as upstream and downstream regions.

DISCUSSION

Genome-skimming data sets can provide a significant source of low-copy nuclear gene sequence data that may be used for multiple downstream applications.

The unexpected depths of genome-skimming data: A case study examining Goodeniaceae floral symmetry genes

PREMISE OF THE STUDY

The use of genome skimming allows systematists to quickly generate large data sets, particularly of sequences in high abundance (e.g., plastomes); however, researchers may be overlooking data in low abundance that could be used for phylogenetic or evo-devo studies. Here, we present a bioinformatics approach that explores the low-abundance portion of genome-skimming next-generation sequencing libraries in the fan-flowered Goodeniaceae.

METHODS

Twenty-four previously constructed Goodeniaceae genome-skimming Illumina libraries were examined for their utility in mining low-copy nuclear genes involved in floral symmetry, specifically the CYCLOIDEA (CYC)-like genes. De novo assemblies were generated using multiple assemblers, and BLAST searches were performed for CYC1, CYC2, and CYC3 genes.

RESULTS

Overall Trinity, SOAPdenovo-Trans, and SOAPdenovo implementing lower k-mer values uncovered the most data, although no assembler consistently outperformed the others. Using SOAPdenovo-Trans across all 24 data sets, we recovered four CYC-like gene groups (CYC1, CYC2, CYC3A, and CYC3B) from a majority of the species. Alignments of the fragments included the entire coding sequence as well as upstream and downstream regions.

DISCUSSION

Genome-skimming data sets can provide a significant source of low-copy nuclear gene sequence data that may be used for multiple downstream applications.

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.

Integrated omics analysis of specialized metabolism in medicinal plants

Medicinal plants are a rich source of highly diverse specialized metabolites with important pharmacological properties. Until recently, plant biologists were limited in their ability to explore the biosynthetic pathways of these metabolites, mainly due to the scarcity of plant genomics resources. However, recent advances in high-throughput large-scale analytical methods have enabled plant biologists to discover biosynthetic pathways for important plant-based medicinal metabolites. The reduced cost of generating omics datasets and the development of computational tools for their analysis and integration have led to the elucidation of biosynthetic pathways of several bioactive metabolites of plant origin. These discoveries have inspired synthetic biology approaches to develop microbial systems to produce bioactive metabolites originating from plants, an alternative sustainable source of medicinally important chemicals. Since the demand for medicinal compounds are increasing with the world's population, understanding the complete biosynthesis of specialized metabolites becomes important to identify or develop reliable sources in the future. Here, we review the contributions of major omics approaches and their integration to our understanding of the biosynthetic pathways of bioactive metabolites. We briefly discuss different approaches for integrating omics datasets to extract biologically relevant knowledge and the application of omics datasets in the construction and reconstruction of metabolic models.