Coupling deep transcriptome analysis with untargeted metabolic profiling in Ophiorrhiza pumila to further the understanding of the biosynthesis of the anti-cancer alkaloid camptothecin and anthraquinones

Transcriptomic divergence of the Rheum palmatum complex derived from top-geoherb and non-geoherb areas provides the insights into geoherbalism properties of rhubarb

Geoherb usually represents high-quality medicinal herbs with better clinical therapeutic effects, and elucidating the geoherbalism is essential for the quality improvement of traditional Chinese Medicine. However, few researches were conducted to clarify the geoherbalism based on a large scale of transcriptomics. In the present study, we compared the transcriptomes of Rheum palmatum complex derived from top-geoherb and non-geoherb areas to show the geoherbalism properties of rhubarb. A total of 412.32 Gb clean reads were obtained with unigene numbers of 100,615 after assembly. Based on the obtained transcriptome datasets, key enzyme-encoding genes involved in the anthraquinones biosynthesis were also obtained. We also found that 21 anthraquinone-related unigenes were differentially expressed between two different groups, and some of these DEGs were correlated to the content accumulation of five free anthraquinones, indicating that the gene expression profiles may promote the geoherbalism formation of rhubarb. In addition, the selective pressure analyses indicated that most paired orthologous genes between these two groups were subject to negative selection, and only a low proportion of orthologs under positive selection were detected. Functional annotation analyses indicated that these positive-selected genes related to the functions such as gene expression, substance transport, stress response and metabolism, indicating that discrepant environment also enhanced the formation of geoherbalism. Our study not only provided insights for the genetic mechanism of geoherbalism of rhubarb, but also laid more genetic cues for the future rhubarb germplasms improvement and utilization.

Transcriptome and metabolome profiling of the medicinal plant Veratrum mengtzeanum reveal key components of the alkaloid biosynthesis

Veratrum mengtzeanum is the main ingredient for Chinese folk medicine known as "Pimacao" due to its unique alkaloids. A diverse class of plant-specific metabolites having key pharmacological activities. There are limited studies on alkaloid synthesis and its metabolic pathways in plants. To elucidate the alkaloid pathway and identify novel biosynthetic enzymes and compounds in V. mengtzeanum, transcriptome and metabolome profiling has been conducted in leaves and roots. The transcriptome of V. mengtzeanum leaves and roots yielded 190,161 unigenes, of which 33,942 genes expressed differentially (DEGs) in both tissues. Three enriched regulatory pathways (isoquinoline alkaloid biosynthesis, indole alkaloid biosynthesis and tropane, piperidine and pyridine alkaloid biosynthesis) and a considerable number of genes such as AED3-like, A4U43, 21 kDa protein-like, 3-O-glycotransferase 2-like, AtDIR19, MST4, CASP-like protein 1D1 were discovered in association with the biosynthesis of alkaloids in leaves and roots. Some transcription factor families, i.e., AP2/ERF, GRAS, NAC, bHLH, MYB-related, C3H, FARI, WRKY, HB-HD-ZIP, C2H2, and bZIP were also found to have a prominent role in regulating the synthesis of alkaloids and steroidal alkaloids in the leaves and roots of V. mengtzeanum. The metabolome analysis revealed 74 significantly accumulated metabolites, with 55 differentially accumulated in leaves compared to root tissues. Out of 74 metabolites, 18 alkaloids were highly accumulated in the roots. A novel alkaloid compound viz; 3-Vanilloylygadenine was discovered in root samples. Conjoint analysis of transcriptome and metabolome studies has also highlighted potential genes involved in regulation and transport of alkaloid compounds. Here, we have presented a comprehensive metabolic and transcriptome profiling of V. mengtzeanum tissues. In earlier reports, only the roots were reported as a rich source of alkaloid biosynthesis, but the current findings revealed both leaves and roots as significant manufacturing factories for alkaloid biosynthesis.

OpNAC1 transcription factor regulates the biosynthesis of the anticancer drug camptothecin by targeting loganic acid O-methyltransferase in Ophiorrhiza pumila

Camptothecin (CPT) is an anticancer pentacyclic quinoline alkaloid widely used to treat cancer patients worldwide. However, the biosynthetic pathway and transcriptional regulation of camptothecin are largely unknown. Ophiorrhiza pumila, the herbaceous plant from the Rubiaceae family, has emerged as a model plant for studying camptothecin biosynthesis and regulation. In this study, a high-quality reference genome of O. pumila with estimated size of ~456.90 Mb was reported, and the accumulation level of camptothecin in roots was higher than that in stems and leaves. Based on its spatial distribution in the plant, we examined gene functions and expression by combining genomics with transcriptomic analysis. Two loganic acid O-methyltransferase (OpLAMTs) were identified in strictosidine-producing plant O. pumila, and enzyme catalysis assays showed that OpLAMT1 and not OpLAMT2 could convert loganic acid into loganin. Further knock-out of OpLAMT1 expression led to the elimination of loganin and camptothecin accumulation in O. pumila hairy roots. Four key residues were identified in OpLAMT1 protein crucial for the catalytic activity of loganic acid to loganin. By co-expression network, we identified a NAC transcription factor, OpNAC1, as a candidate gene for regulating camptothecin biosynthesis. Transgenic hairy roots and biochemical assays demonstrated that OpNAC1 suppressed OpLAMT1 expression. Here, we reported on two camptothecin metabolic engineering strategies paving the road for industrial-scale production of camptothecin in CPT-producing plants.

Comparative transcriptome analysis of different tissues of Rheum tanguticum Maxim. ex Balf. (Polygonaceae) reveals putative genes involved in anthraquinone biosynthesis

Rheum tanguticum is a perennial herb and an important medicinal plant, with anthraquinones as its main bioactive compounds. However, the specific pathway of anthraquinone biosynthesis in rhubarb is still unclear. The accumulation of anthraquinones in different tissues (root, leaf, stem and seed) of R. tanguticum revealed considerable variation, suggesting possible differences in metabolite biosynthetic pathways and accumulation among various tissues. To better illustrate the biosynthetic pathway of anthraquinones, we assembled transcriptome sequences from the root, leaf, stem and seed tissues yielding 157,564 transcripts and 88,142 unigenes. Putative functions could be assigned to 56,911 unigenes (64.57%) based on BLAST searches against annotation databases, including GO, KEGG, Swiss-Prot, NR, and Pfam. In addition, putative genes involved in the biosynthetic pathway of anthraquinone were identified. The expression profiles of nine unigenes involved in anthraquinone biosynthesis were verified in different tissues of R. tanguticum by qRT-PCR. Various transcription factors, including bHLH, MYB_related, and C2H2, were identified by searching unigenes against plantTFDB. This is the first transcriptome analysis of different tissues of R. tanguticum and can be utilized to describe the genes involved in the biosynthetic pathway of anthraquiones, understanding the molecular mechanism of active compounds in R. tanguticum.

Biotechnology for micropropagation and camptothecin production in Ophiorrhiza sp

Camptothecin (CPT) is a monoterpenoid-alkaloid, an anticancer compound from plant. Ever since its discovery in 1996 from the bark of Camptotheca acuminata, various researches have been conducted for enhancing its production. CPT has also been reported in several other species belonging to the plant families Icacinaceae, Rubiaceae, Apocynaceae, Nyssaceae, Betulaceae, Violaceae, Meliaceae, and Gelseminaceae. Out of these, Ophiorrhiza sp. (Rubiaceae) is the next possible candidate for sustainable CPT production after C. acuminata and Nothapodytes nimoonia. Various biotechnological-studies have been conducted on Ophiorrhiza sp. for searching the elite species and the most optimal strategies for CPT production. The genus Ophiorrhiza has been used as medicines for antiviral, antifungal, antimalarial, and anticancer activities. Phytochemical analysis has revealed the presence of alkaloids, flavonoids, triterpenes, and CPT from the plant. Because of the presence of CPT and its herbaceous habit, Ophiorrhiza sp. has now become a hot topic in research area. Currently, for mass production of the elite spp., tissue culture techniques have been implemented. In the past decades, several researchers have contributed on the diversity assessment, phytochemical analysis, mass production, and in vitro production of CPT in Ophiorrhiza sp. In this paper, we review the on the biotechnological strategies, optimal culture medium, micropropagation of Ophiorrhiza sp., effect of PGR on shoot formation, rhizogenesis, callus formation, and enhanced production of CPT for commercial use. KEY POINTS: • Latest literature on in vitro propagation of Ophiorrhiza sp. • Biotechnological production of camptothecin and related compounds • Optimization, elicitation, and transgenic studies in Ophiorrhiza sp.

De Novo Transcriptome Analysis Reveals Putative Genes Involved in Anthraquinone Biosynthesis in Rubia yunnanensis

Rubia yunnanensis Diels (R. yunnanensis), a Chinese perennial plant, is well-known for its medicinal values such as rheumatism, contusion, and anemia. It is rich in bioactive anthraquinones, but the biosynthetic pathways of anthraquinones in R. yunnanensis remain unknown. To investigate genes involved in anthraquinone biosynthesis in R. yunnanensis, we generated a de novo transcriptome of R. yunnanensis using the Illumina HiSeq 2500 sequencing platform. A total of 636,198 transcripts were obtained, in which 140,078 transcripts were successfully annotated. A differential gene expression analysis identified 15 putative genes involved in anthraquinone biosynthesis. Additionally, the hairy roots of R. yunnanensis were treated with 200 µM Methyl Jasmonate (MeJA). The contents of six bioactive anthraquinones and gene expression levels of 15 putative genes were measured using ultra performance liquid chromatography coupled with mass spectrometry (UPLC-MS/MS) and real-time quantitative polymerase chain reaction (RT-qPCR), respectively. The results showed that the expressions levels for 11 of the 15 genes and the contents of two of six anthraquinones significantly increased by MeJA treatment. Pearson's correlation analyses indicated that the expressions of 4 of the 15 putative genes were positively correlated with the contents of rubiquinone (Q3) and rubiquinone-3-O-β-d-xylopranosyl-(1→6)-β-d-glucopyranoside (Q20). This study reported the first de novo transcriptome of R. yunnanensis and shed light on the anthraquinone biosynthesis and genetic information for R. yunnanensis.

Allelopathy as an evolutionary game

In plants, most competition is resource competition, where one plant simply preempts the resources away from its neighbors. Interference competition, as the name implies, is a form of direct interference to prevent resource access. Interference competition is common among animals that can physically fight, but in plants, one of the main mechanisms of interference competition is allelopathy. Allelopathic plants release cytotoxic chemicals into the environment which can increase their ability to compete with surrounding organisms for limited resources. The circumstances and conditions favoring the development and maintenance of allelochemicals, however, are not well understood. Particularly, despite the obvious benefits of allelopathy, current data suggest it seems to have only rarely evolved. To gain insight into the cost and benefit of allelopathy, we have developed a 2 × 2 matrix game to model the interaction between plants that produce allelochemicals and plants that do not. Production of an allelochemical introduces novel cost associated with both synthesis and detoxifying a toxic chemical but may also convey a competitive advantage. A plant that does not produce an allelochemical will suffer the cost of encountering one. Our model predicts three cases in which the evolutionarily stable strategies are different. In the first, the nonallelopathic plant is a stronger competitor, and not producing allelochemicals is the evolutionarily stable strategy. In the second, the allelopathic plant is the better competitor, and production of allelochemicals is the more beneficial strategy. In the last case, neither is the evolutionarily stable strategy. Instead, there are alternating stable states, depending on whether the allelopathic or nonallelopathic plant arrived first. The generated model reveals circumstances leading to the evolution of allelochemicals and sheds light on utilizing allelochemicals as part of weed management strategies. In particular, the wide region of alternative stable states in most parameterizations, combined with the fact that the absence of allelopathy is likely the ancestral state, provides an elegant answer to the question of why allelopathy seems to rarely evolve despite its obvious benefits. Allelopathic plants can indeed outcompete nonallelopathic plants, but this benefit is simply not great enough to allow them to go to fixation and spread through the population. Thus, most populations would remain purely nonallelopathic.

Identification of WRKY transcription factors involved in regulating the biosynthesis of the anti-cancer drug camptothecin in Ophiorrhiza pumila

Camptothecin is a chemotherapeutic drug widely used to treat various cancers. Ophiorrhiza pumila is an ideal plant model for the study of camptothecin production, with various advantages for studying camptothecin biosynthesis and regulation. The DNA-binding WRKY transcription factors have a key regulatory role in secondary metabolite biosynthesis in plants. However, little is currently known about their involvement in camptothecin biosynthesis in O. pumila. We identified 46 OpWRKY genes unevenly distributed on the 11 chromosomes of O. pumila. Phylogenetic and multiple sequence alignment analyses divided the OpWRKY proteins into three subfamilies. Based on spatial expression and co-expression, we targeted the candidate gene OpWRKY6. Overexpression of OpWRKY6 significantly reduced the accumulation of camptothecin compared with the control. Conversely, camptothecin accumulation increased in OpWRKY6 knockout lines. Further biochemical assays showed that OpWRKY6 negatively regulates camptothecin biosynthesis from both the iridoid and shikimate pathways by directly downregulating the gene expression of OpGES, Op10HGO, Op7DLH, and OpTDC. Our data provide direct evidence for the involvement of WRKYs in the regulation of camptothecin biosynthesis and offer valuable information for enriching the production of camptothecin in plant systems.

Biotechnology of camptothecin production in Nothapodytes nimmoniana, Ophiorrhiza sp. and Camptotheca acuminata

Cancer is probably the deadliest human disease in recent years. In the past few years, rapid clinical progress has been made in the field of anticancer drug development. Plant secondary metabolites have been noted as extremely efficacious as promising natural source for anticancer therapy for many years. Camptothecin (CPT) is one of the popularly used anti-tumor drugs possessing clinically proven properties against a plethora of human malignancies that include ovarian and colorectal cancers. For the first time, CPT was obtained from the extracts of a Chinese medicinal tree, Camptotheca acuminata Decne. from the family Cornaceae. Subsequently, CPT was also isolated from the bark of Nothapodytes foetida (Wight) Sleumer (Icacinaceae). However, the availability of enough natural sources for obtaining CPT is a major constraint. Due to overexploitation and harvesting, loss of habitat, excessive trading, and unfavorable environmental factors, the natural source of CPT has become extinct or extremely limited and hence they are red listed under endangered species. Conventional propagation has also failed to meet the ever-expanding demand for CPT production. With this, biotechnological toolkits have constantly been used as a boon to produce sustainable source, utilization, and ex situ conservation of medicinal plants. The approaches serve as a supplement to traditional agriculture in the mass production of plant metabolites with potent bioactivities. Non-availability of enough anticancer medicine and the requirement to satisfy current demands need a sustainable source of CPT. With this background, we present a comprehensive review on CPT discovery, its occurrence in the plant kingdom, biosynthesis, phytochemistry, pharmacological properties, clinical studies, patterns of CPT accumulation, and biotechnological aspects of CPT production in three plants, viz., N. nimmoniana, Ophiorrhiza species, and C. acuminata.Key points• Biotechnological approaches on production of camptothecin from Nothapodytes nimmoniana, Ophiorrhiza species, and Camptotheca acuminata• In vitro propagation of camptothecin-producing plants• Genetic diversity and transgenic research on camptothecin-producing plants.

Tissue-specific transcriptome for Rheum tanguticum reveals candidate genes related to the anthraquinones biosynthesis

Rheum tanguticum (Maxim. ex Regel) Maxim. ex Balf. is a herbaceous perennial plant indigenous to China, and its root and rhizomes were usually used as an important traditional Chinese medicine. However, the genomic resources are still scarce for R. tanguticum and even for Rheum genus. Transcriptome datasets from different tissues of R. tanguticum were obtained to screen the genes related to anthraquinones biosynthesis, and five free anthraquinones were also determined. Nine cDNA libraries of roots, stems and leaves were generated, and a total of 272 million high-quality reads were assembled into 257,942 unigenes. Based on the functional annotation, A total of 227 candidate enzyme genes involved in the MVA, MEP, shikimate and polyketide pathways were identified, and several differentially expressed genes found functionally associated with anthraquinones biosynthesis showed distinct tissue-specific expression patterns. Especially, we found that the expression levels of PKS III genes might result in the content differences of free anthraquinones in different tissues of R. tanguticum. Besides, 137,400 SSR loci were identified, and 64,081 SSR primer pairs were successfully designed based on these loci. Our results not only provide cues for the genetic mechanism of anthraquinone content differences in different tissues of R. tanguticum, but also lay genomic foundation for the subsequent genetic engineering and breeding for Rheum species.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01099-8.

Divergent camptothecin biosynthetic pathway in Ophiorrhiza pumila

BACKGROUND

The anticancer drug camptothecin (CPT), first isolated from Camptotheca acuminata, was subsequently discovered in unrelated plants, including Ophiorrhiza pumila. Unlike known monoterpene indole alkaloids, CPT in C. acuminata is biosynthesized via the key intermediate strictosidinic acid, but how O. pumila synthesizes CPT has not been determined.

RESULTS

In this study, we used nontargeted metabolite profiling to show that 3α-(S)-strictosidine and 3-(S), 21-(S)-strictosidinic acid coexist in O. pumila. After identifying the enzymes OpLAMT, OpSLS, and OpSTR as participants in CPT biosynthesis, we compared these enzymes to their homologues from two other representative CPT-producing plants, C. acuminata and Nothapodytes nimmoniana, to elucidate their phylogenetic relationship. Finally, using labelled intermediates to resolve the CPT biosynthesis pathway in O. pumila, we showed that 3α-(S)-strictosidine, not 3-(S), 21-(S)-strictosidinic acid, is the exclusive intermediate in CPT biosynthesis.

CONCLUSIONS

In our study, we found that O. pumila, another representative CPT-producing plant, exhibits metabolite diversity in its central intermediates consisting of both 3-(S), 21-(S)-strictosidinic acid and 3α-(S)-strictosidine and utilizes 3α-(S)-strictosidine as the exclusive intermediate in the CPT biosynthetic pathway, which differs from C. acuminata. Our results show that enzymes likely to be involved in CPT biosynthesis in O. pumila, C. acuminata, and N. nimmoniana have evolved divergently. Overall, our new data regarding CPT biosynthesis in O. pumila suggest evolutionary divergence in CPT-producing plants. These results shed new light on CPT biosynthesis and pave the way towards its industrial production through enzymatic or metabolic engineering approaches.

A high-quality genome assembly of Morinda officinalis, a famous native southern herb in the Lingnan region of southern China

Morinda officinalis is a well-known medicinal and edible plant that is widely cultivated in the Lingnan region of southern China. Its dried roots (called bajitian in traditional Chinese medicine) are broadly used to treat various diseases, such as impotence and rheumatism. Here, we report a high-quality chromosome-scale genome assembly of M. officinalis using Nanopore single-molecule sequencing and Hi-C technology. The assembled genome size was 484.85 Mb with a scaffold N50 of 40.97 Mb, and 90.77% of the assembled sequences were anchored on eleven pseudochromosomes. The genome includes 27,698 protein-coding genes, and most of the assemblies are repetitive sequences. Genome evolution analysis revealed that M. officinalis underwent core eudicot γ genome triplication events but no recent whole-genome duplication (WGD). Likewise, comparative genomic analysis showed no large-scale structural variation after species divergence between M. officinalis and Coffea canephora. Moreover, gene family analysis indicated that gene families associated with plant-pathogen interactions and sugar metabolism were significantly expanded in M. officinalis. Furthermore, we identified many candidate genes involved in the biosynthesis of major active components such as anthraquinones, iridoids and polysaccharides. In addition, we also found that the DHQS, GGPPS, TPS-Clin, TPS04, sacA, and UGDH gene families-which include the critical genes for active component biosynthesis-were expanded in M. officinalis. This study provides a valuable resource for understanding M. officinalis genome evolution and active component biosynthesis. This work will facilitate genetic improvement and molecular breeding of this commercially important plant.

Crosstalk of Multi-Omics Platforms with Plants of Therapeutic Importance

From time immemorial, humans have exploited plants as a source of food and medicines. The World Health Organization (WHO) has recorded 21,000 plants with medicinal value out of 300,000 species available worldwide. The promising modern "multi-omics" platforms and tools have been proven as functional platforms able to endow us with comprehensive knowledge of the proteome, genome, transcriptome, and metabolome of medicinal plant systems so as to reveal the novel connected genetic (gene) pathways, proteins, regulator sequences and secondary metabolite (molecule) biosynthetic pathways of various drug and protein molecules from a variety of plants with therapeutic significance. This review paper endeavors to abridge the contemporary advancements in research areas of multi-omics and the information involved in decoding its prospective relevance to the utilization of plants with medicinal value in the present global scenario. The crosstalk of medicinal plants with genomics, transcriptomics, proteomics, and metabolomics approaches will be discussed.

Convergent evolution of plant specialized 1,4-naphthoquinones: metabolism, trafficking, and resistance to their allelopathic effects

Plant 1,4-naphthoquinones encompass a class of specialized metabolites known to mediate numerous plant-biotic interactions. This class of compounds also presents a remarkable case of convergent evolution. The 1,4-naphthoquinones are synthesized by species belonging to nearly 20 disparate orders spread throughout vascular plants, and their production occurs via one of four known biochemically distinct pathways. Recent developments from large-scale biology and genetic studies corroborate the existence of multiple pathways to synthesize plant 1,4-naphthoquinones and indicate that extraordinary events of metabolic innovation and links to respiratory and photosynthetic quinone metabolism probably contributed to their independent evolution. Moreover, because many 1,4-naphthoquinones are excreted into the rhizosphere and they are highly reactive in biological systems, plants that synthesize these compounds also needed to independently evolve strategies to deploy them and to resist their effects. In this review, we highlight new progress made in understanding specialized 1,4-naphthoquinone biosynthesis and trafficking with a focus on how these discoveries have shed light on the convergent evolution and diversification of this class of compounds in plants. We also discuss how emerging themes in metabolism-based herbicide resistance may provide clues to mechanisms plants employ to tolerate allelopathic 1,4-naphthoquinones.

Chromosome-level genome assembly of Ophiorrhiza pumila reveals the evolution of camptothecin biosynthesis

Plant genomes remain highly fragmented and are often characterized by hundreds to thousands of assembly gaps. Here, we report chromosome-level reference and phased genome assembly of Ophiorrhiza pumila, a camptothecin-producing medicinal plant, through an ordered multi-scaffolding and experimental validation approach. With 21 assembly gaps and a contig N50 of 18.49 Mb, Ophiorrhiza genome is one of the most complete plant genomes assembled to date. We also report 273 nitrogen-containing metabolites, including diverse monoterpene indole alkaloids (MIAs). A comparative genomics approach identifies strictosidine biogenesis as the origin of MIA evolution. The emergence of strictosidine biosynthesis-catalyzing enzymes precede downstream enzymes' evolution post γ whole-genome triplication, which occurred approximately 110 Mya in O. pumila, and before the whole-genome duplication in Camptotheca acuminata identified here. Combining comparative genome analysis, multi-omics analysis, and metabolic gene-cluster analysis, we propose a working model for MIA evolution, and a pangenome for MIA biosynthesis, which will help in establishing a sustainable supply of camptothecin.

Hairy root culture technology: applications, constraints and prospect

Hairy root (HR) culture, a successful biotechnology combining in vitro tissue culture with recombinant DNA machinery, is intended for the genetic improvement of plants. This technology has been put to use since the last three decades for genetic advancement of medicinal and aromatic plants and also to harvest the economical products in the form of secondary metabolites that are significantly important for their ethnobotanical and pharmacological properties. It also provides an efficient way out for the quicker extraction and quantification of the valuable phytochemicals. The current review provides an account of the in vitro HR culture technology and its wide-scale applications in the field of research as well as in pharmaceutical industries. Different facets of HR with respect to the culture establishment, phytochemical production as well as research investigations concerning the areas of gene manipulation, biotransformation of the secondary metabolites, phytoremediation, their industrial utilisations and different problems encountered during the application of this technology have been covered in this appraisal. Eventually, an idea has been provided on HR about the recent trends on the progress of this technology that may open up newer prospects in near future and calls for further research and explorations in this field. KEY POINTS: • Genetic engineering-based HR culture aims towards enhanced secondary metabolite production. • This review explores an insight in the HR technology and its multi-faceted approaches. • Up-to-date ground-breaking research applications and constraints of HR culture are discussed.

Comparative transcriptome analysis to identify candidate genes involved in 2-methoxy-1,4-naphthoquinone (MNQ) biosynthesis in Impatiens balsamina L

Impatiens balsamina L. is a tropical ornamental and traditional medicinal herb rich in natural compounds, especially 2-methoxy-1,4-naphthoquinone (MNQ) which is a bioactive compound with tested anticancer activities. Characterization of key genes involved in the shikimate and 1,4-dihydroxy-2-naphthoate (DHNA) pathways responsible for MNQ biosynthesis and their expression profiles in I. balsamina will facilitate adoption of genetic/metabolic engineering or synthetic biology approaches to further increase production for pre-commercialization. In this study, HPLC analysis showed that MNQ was present in significantly higher quantities in the capsule pericarps throughout three developmental stages (early-, mature- and postbreaker stages) whilst its immediate precursor, 2-hydroxy-1,4-naphthoquinone (lawsone) was mainly detected in mature leaves. Transcriptomes of I. balsamina derived from leaf, flower, and three capsule developmental stages were generated, totalling 59.643 Gb of raw reads that were assembled into 94,659 unigenes (595,828 transcripts). A total of 73.96% of unigenes were functionally annotated against seven public databases and 50,786 differentially expressed genes (DEGs) were identified. Expression profiles of 20 selected genes from four major secondary metabolism pathways were studied and validated using qRT-PCR method. Majority of the DHNA pathway genes were found to be significantly upregulated in early stage capsule compared to flower and leaf, suggesting tissue-specific synthesis of MNQ. Correlation analysis identified 11 candidate unigenes related to three enzymes (NADH-quinone oxidoreductase, UDP-glycosyltransferases and S-adenosylmethionine-dependent O-methyltransferase) important in the final steps of MNQ biosynthesis based on genes expression profiles consistent with MNQ content. This study provides the first molecular insight into the dynamics of MNQ biosynthesis and accumulation across different tissues of I. balsamina and serves as a valuable resource to facilitate further manipulation to increase production of MNQ.

Transcriptome-wide identification, characterization, and phylogenomic analysis of cytochrome P450s from Nothapodytes nimmoniana reveal candidate genes involved in the camptothecin biosynthetic pathway

The plant Nothapodytes nimmoniana is an important source of camptothecin (CPT), an anticancer compound widely used in the treatment of colorectal, lung, and ovarian cancers. CPT is biosynthesized by the combination of the seco-iridoid and indole pathways in plants. The majority of the biosynthetic steps and associated genes still remain unknown. Certain reactions in the seco-iridoid pathway are catalyzed by cytochrome P450 enzymes. Hence, identifying transcriptionally active cytochrome P450 genes becomes essential in the elucidation of the CPT biosynthetic pathway. Here, we report the identification of 94 cytochrome P450s from the assembled transcriptomic data from leaf and root tissues of N. nimmoniana. The identified cytochrome P450 genes were full length and possessed all four conserved characteristic signature motifs of cytochrome P450 genes. Phylogenetic analysis of the protein sequences revealed their evolution and diversification and further categorized them into A-type (52.12%) and non-A-type (47.87%) cytochrome P450s. These 94 sequences represent 38 families and 63 subfamilies of cytochrome P450s. We also compared the transcriptional activity of identified cytochrome P450s with the expression of their homologs in the CPT-producing plant Ophiorrhiza pumila. Based on expression profiles and quantitative PCR validation, we propose NnCYP81CB1 and NnCYP89R1 as candidate cytochrome P450 genes involved in camptothecin biosynthesis in N. nimmoniana.

Molecular characterization and overexpression analyses of secologanin synthase to understand the regulation of camptothecin biosynthesis in Nothapodytes nimmoniana (Graham.) Mabb

Camptothecin is a high-value anti-cancerous compound produced in many taxonomically unrelated species. Its biosynthesis involves a complex network of pathways and a diverse array of intermediates. Here, we report the functional characterization and regulation of secologanin synthase (NnCYP72A1), a cytochrome P450 involved in camptothecin biosynthesis from Nothapodytes nimmoniana. It comprises an open reading frame of 1566 bp in length. Heterologous expression in Saccharomyces cerevisiae and in vitro enzymatic assays using loganin as substrate confirmed the formation of secologanin. In planta transient overexpression analysis of NnCYP72A1 resulted in 4.21- and 2.73-fold increase in transcript levels of NnCYP72A1 on days 3 and 6 respectively. Phytochemical analysis of transformed tissues revealed ~ 1.13-1.43- and 2.02-2.86-fold increase in secologanin and CPT accumulation, respectively. Furthermore, promoter analysis of NnCYP72A1 resulted in the identification of several potential cis-regulatory elements corresponding to different stress-related components. Methyl jasmonate, salicylic acid, and wounding treatments resulted in considerable modulation of mRNA transcripts of NnCYP72A1 gene. Chemical analysis of elicitor-treated samples showed a significant increase in CPT content which was concordant with the mRNA transcript levels. Overall, the functional characterization and overexpression of NnCYP72A1 may plausibly enhance the pathway intermediates and serve as prognostic tool for enhancing CPT accumulation.

Ethanol extract of Ophiorrhiza pumila suppresses liver cancer cell proliferation and migration

Background

Ophiorrhiza pumila, belonging to the genus Ophiorrhiza (Rubiaceae), is distributed throughout tropical and subtropical Asia. In this study, we evaluated for the first time the anti-proliferation and anti-migration effects of ethanol extract of O. pumila (OPE) on HepG2 and SMMC-7721 cells, and explored the related mechanism.

Methods

OPE was prepared by percolation with 95% ethanol and its main compounds were analyzed by HPLC-MS². The anti-proliferation effect of OPE was evaluated by the CCK-8 assay and colony formation assay. Cell cycle distribution, apoptosis, and reactive oxygen species (ROS) level were detected by flow cytometry. Migration and invasion abilities were detected by Transwell migration/invasion assays. The expression of correlated proteins was determined using western blotting.

Results

A total of 5 tentative compounds were identified from OPE, including pumiloside, deoxypumiloside, camptothecin, aknadinine, and β-stigmasterol. OPE displayed strong cytostatic effects on HepG2 and SMMC-7721 cells. OPE induced G2/M phase cell cycle arrest, increased apoptosis, and augmented ROS production in these cell lines. In addition, OPE possessed a significant inhibition on cell migration and invasion by reduction of MMP-9 and MMP-2 expression. Moreover, OPE significantly suppressed the phosphorylation of p65.

Conclusions

Our data showed that OPE suppresses liver cancer cell proliferation and migration, which is possibly involved with the inhibition of the NF-κB pathway.

Transcriptome Analysis of Different Tissues Reveals Key Genes Associated With Galanthamine Biosynthesis in Lycoris longituba

L ycoris longituba is a traditional medicinal plant containing the bioactive compound galanthamine (Gal), a type of Amaryllidaceae alkaloid and can be used to treat Alzheimer's disease. However, research on its genome or transcriptome and associated genes in the biosynthetic pathway is incomplete. In this study, we estimated the nuclear genome size of this species to be 29.33 Gb by flow cytometry. Then, RNA sequencing of the leaves, roots, and bulbs of L. longituba was carried out. After de novo assembly, 474,589 all-transcripts and 333,440 all-unigenes were finally generated. In addition, the differentially expressed genes (DEGs) were identified, and genes involved in the galanthamine metabolic pathway encoding tyrosine decarboxylase (TYDC), phenylalanine ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H), p-coumarate 3-hydroxylase (C3H), norbelladine synthase (NBS), norbelladine 4'-O-methyltransferase (OMT), and noroxomaritidine synthase (CYP96T1) were detected and validated by real-time quantitative PCR analysis. One candidate gene, Lycoris longituba O-Methyltransferase (LlOMT), was identified in the proposed galanthamine biosynthetic pathway. Sequence analysis showed that LlOMT is a class I OMT. LlOMT is localized in the cytoplasm, and biochemical analysis indicated that the recombinant LlOMT catalyzes norbelladine to generate 4'-O-methylnorbelladine. The protoplast transformation result showed that the overexpression of LlOMT could increase the Gal content. Our results indicate that LlOMT may play a role in galanthamine biosynthesis in L. longituba. This work provides a useful resource for the metabolic engineering of Amaryllidaceae alkaloids.

Possible clues for camptothecin biosynthesis from the metabolites in camptothecin-producing plants

The plant derived camptothecin (CPT) is a pentacyclic pyrroloquinoline alkaloid with unique antitumor activity. Successive discoveries of new CPT-producing plants occurred in recent years due to market demands. The scattered distribution among angiosperms drew researchers' attention. The aim of this review is to appraise the literature available to date for CPT distribution and the phytochemistry of these CPT-producing plants. Metabolite comparative analyses between the plants were also conducted for tracking of possible clues for CPT biosynthesis. Forty-three plant species in total were reported to possess CPT-producing capability, and one hundred twenty-five alkaloids classified into three major categories are summarized herein. Metabolite comparative analysis between these plants suggests the probability that the formation of the central intermediate for CPT biosynthesis has multiple origins. A more complete biogenetic reasoning for CPT and its structural homolog was delineated based on this fragmentary phytochemical evidence from a chemical point of view. Furthermore, an in-house compound database was constructed for further metabolomic analysis.

De novo genome assembly of Camptotheca acuminata, a natural source of the anti-cancer compound camptothecin

Camptotheca acuminata is 1 of a limited number of species that produce camptothecin, a pentacyclic quinoline alkaloid with anti-cancer activity due to its ability to inhibit DNA topoisomerase. While transcriptome studies have been performed previously with various camptothecin-producing species, no genome sequence for a camptothecin-producing species is available to date. We generated a high-quality de novo genome assembly for C. acuminata representing 403 174 860 bp on 1394 scaffolds with an N50 scaffold size of 1752 kbp. Quality assessments of the assembly revealed robust representation of the genome sequence including genic regions. Using a novel genome annotation method, we annotated 31 825 genes encoding 40 332 gene models. Based on sequence identity and orthology with validated genes from Catharanthus roseus as well as Pfam searches, we identified candidate orthologs for genes potentially involved in camptothecin biosynthesis. Extensive gene duplication including tandem duplication was widespread in the C. acuminata genome, with 2571 genes belonging to 997 tandem duplicated gene clusters. To our knowledge, this is the first genome sequence for a camptothecin-producing species, and access to the C. acuminata genome will permit not only discovery of genes encoding the camptothecin biosynthetic pathway but also reagents that can be used for heterologous expression of camptothecin and camptothecin analogs with novel pharmaceutical applications.

[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.

De novo transcriptome analyses reveals putative pathway genes involved in biosynthesis and regulation of camptothecin in Nothapodytes nimmoniana (Graham) Mabb

Comprehensive transcriptome analysis of leaf and root tissues of Nothapodytes nimmoniana unravels several putative pathway genes, transcription factors and CYPs related to camptothecin (CPT) biosynthesis. Additionally, post-transcriptional suppression by artificial microRNA (aMIR) of NnCYP76B6 (geraniol 10-hydroxylase) suggests its role in CPT biosynthesis. Tissue-specific LC-MS/MS analysis revealed the presence of secologanin as the central intermediate of MIA pathway in N. nimmoniana. Nothapodytes nimmoniana is a rich source of potent anticancer drug camptothecin (CPT) whose biosynthetic pathway is unresolved due to the lack of genomic and transcriptomic information. Present investigation entails deep transcriptome analysis of N. nimmoniana which led to identification of putative pathway genes and regulatory components involved in CPT biosynthesis. Using Illumina HiSeq 2500 sequencing platform a total of 31,172,889 (6.23 Gb) and 31,218,626 (6.24 Gb) raw reads were generated from leaf and root wood, respectively. These were assembled de novo into 138,183 unique contigs. Additionally, 16 cytochrome P450 transcripts related to secondary metabolism were also identified. Further, transcriptome data pool presented 1683 putative transcription factors of which transcripts corresponding to WRKY TFs were the most abundant (14.14%). A total of 2741 transcripts were differentially expressed out of which 478 contigs showed downregulation in root wood and 2263 contigs were up-regulated. Further, comparative analyses of 17 genes involved in CPT biosynthetic pathway were validated by qRT-PCR. On basis of intermediates, two distinct seco-iridoid pathways are involved in the biosynthesis of monoterpene indole alkaloids either through multiple isomers of strictosidinic acid or strictosidine. Tissue-specific LC-MS/MS analysis revealed the presence of secologanin as the central intermediate of MIA pathway in N. nimmoniana. Geraniol-10 hydroxylase (NnCYP76B6) an important enzyme in CPT biosynthesis which specifically shunts geraniol into the secologanin pathway was also cloned from the trancriptome resource. In planta transient expression of NnCYP76B6 showed a significant enhancement in mRNA transcript levels coincident with enhanced CPT accumulation. Further, artificial microRNA (aMIR) mediated downregulation of NnCYP76B6 resulted in reduction of mRNA transcript levels as well as CPT content in comparison to control. These empirical results suggest a plausible regulatory role for NnCYP76B6 in CPT biosynthesis and also establish a valuable repository for deciphering various structural, rate limiting and regulatory genes of CPT biosynthetic pathway.

Transcriptome and metabolome profiling of Narcissus pseudonarcissus 'King Alfred' reveal components of Amaryllidaceae alkaloid metabolism

Amaryllidaceae alkaloids (AAs) represent a diverse class of plant specialized metabolites and many display potent pharmacological activities. The AA metabolic pathway is poorly understood and resources are minimal. To enable AA pathway elucidation and novel biosynthetic enzymes discovery, we generated comprehensive metabolomic and corresponding transcriptomic datasets from different tissues of Narcissus pseudonarcissus ‘King Alfred’. In this study, we performed untargeted UPLC-QTOF-MS metabolite analysis from different tissues, which generated exhaustive list of compounds, including several AAs, most predominant and diverse in bulbs. RNA sequencing of N. pseudonarcissus ‘King Alfred’ bulbs yielded 195,347 transcripts, after assembly. Top expressed genes belong to process like metabolism, survival, and defense including alkaloid biosynthetic genes. The transcriptome contained complete sequences for all proposed genes encoding AA-biosynthetic enzymes such as tyrosine decarboxylase (TYDC1 and TYDC2), phenylalanine ammonia-lyase (PAL1 and PAL2) and phenolic acids hydroxylases (C4H and C3H) to name a few. Furthermore, transcriptome data were validated using RT-qPCR analysis and expression study in different tissues of N. pseudonarcissus ‘King Alfred’ was performed. Here, we present the first comprehensive metabolome and transcriptome study from N. pseudonarcissus ‘King Alfred’ providing invaluable resources for metabolic engineering and biotechnological applications.

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.

Metabolite Diversity in Alkaloid Biosynthesis: A Multilane (Diastereomer) Highway for Camptothecin Synthesis in Camptotheca acuminata

Camptothecin is a monoterpene indole alkaloid (MIA) used to produce semisynthetic antitumor drugs. We investigated camptothecin synthesis in Camptotheca acuminata by combining transcriptome and expression data with reverse genetics, biochemistry, and metabolite profiling. RNAi silencing of enzymes required for the indole and seco-iridoid (monoterpene) components identified transcriptional crosstalk coordinating their synthesis in roots. Metabolite profiling and labeling studies of wild-type and RNAi lines identified plausible intermediates for missing pathway steps and demonstrated nearly all camptothecin pathway intermediates are present as multiple isomers. Unlike previously characterized MIA-producing plants, C. acuminata does not synthesize 3-α(S)-strictosidine as its central MIA intermediate and instead uses an alternative seco-iridoid pathway that produces multiple isomers of strictosidinic acid. NMR analysis demonstrated that the two major strictosidinic acid isomers are (R) and (S) diastereomers at their glucosylated C21 positions. The presence of multiple diastereomers throughout the pathway is consistent with their use in synthesis before finally being resolved to a single camptothecin isomer after deglucosylation, much as a multilane highway allows parallel tracks to converge at a common destination. A model "diastereomer" pathway for camptothecin biosynthesis in C. acuminata is proposed that fundamentally differs from previously studied MIA pathways.

Renaissance in phytomedicines: promising implications of NGS technologies

Medicinal plant research is growing significantly in faith to discover new and more biologically compatible phytomedicines. Deposition of huge genome/trancriptome sequence data assisted by NGS technologies has revealed the new possibilities for producing upgraded bioactive molecules in medicinal plants. Growing interest of investors and consumers in the herbal drugs raises the need for extensive research to open the facts and details of every inch of life canvas of medicinal plants to produce improved quality of phytomedicines. As in agriculture crops, knowledge emergence from medicinal plant's genome/transcriptome, can be used to assure their amended quality and these improved varieties are then transported to the fields for cultivation. Genome studies generate huge sequence data which can be exploited further for obtaining information regarding genes/gene clusters involved in biosynthesis as well as regulation. This can be achieved rapidly at a very large scale with NGS platforms. Identification of new RNA molecules has become possible, which can lead to the discovery of novel compounds. Sequence information can be combined with advanced phytochemical and bioinformatics tools to discover functional herbal drugs. Qualitative and quantitative analysis of small RNA species put a light on the regulatory aspect of biosynthetic pathways for phytomedicines. Inter or intra genomic as well as transcriptomic interactive processes for biosynthetic pathways can be elucidated in depth. Quality management of herbal material will also become rapid and high throughput. Enrichment of sequence information will be used to engineer the plants to get more efficient phytopharmaceuticals. The present review comprises of role of NGS technologies to boost genomic studies of pharmaceutically important plants and further, applications of sequence information aiming to produce enriched phytomedicines. Emerging knowledge from the medicinal plants genome/transcriptome can give birth to deep understanding of the processes responsible for biosynthesis of medicinally important compounds.

Transcriptome and metabolome analysis in shoot and root of Valeriana fauriei

BACKGROUND

Valeriana fauriei is commonly used in the treatment of cardiovascular diseases in many countries. Several constituents with various pharmacological properties are present in the roots of Valeriana species. Although many researches on V. fauriei have been done since a long time, further studies in the discipline make a limit due to inadequate genomic information. Hence, Illumina HiSeq 2500 system was conducted to obtain the transcriptome data from shoot and root of V. fauriei.

RESULTS

A total of 97,595 unigenes were noticed from 346,771,454 raw reads after preprocessing and assembly. Of these, 47,760 unigens were annotated with Uniprot BLAST hits and mapped to COG, GO and KEGG pathway. Also, 70,013 and 88,827 transcripts were expressed in root and shoot of V. fauriei, respectively. Among the secondary metabolite biosynthesis, terpenoid backbone and phenylpropanoid biosynthesis were large groups, where transcripts was involved. To characterize the molecular basis of terpenoid, carotenoid, and phenylpropanoid biosynthesis, the levels of transcription were determined by qRT-PCR. Also, secondary metabolites content were measured using GC/MS and HPLC analysis for that gene expression correlated with its accumulation respectively between shoot and root of V. fauriei.

CONCLUSIONS

We have identified the transcriptome using Illumina HiSeq system in shoot and root of V. fauriei. Also, we have demonstrated gene expressions associated with secondary metabolism such as terpenoid, carotenoid, and phenylpropanoid.

Function of AP2/ERF Transcription Factors Involved in the Regulation of Specialized Metabolism in Ophiorrhiza pumila Revealed by Transcriptomics and Metabolomics

The hairy roots (HR) of Ophiorrhiza pumila produce camptothecin (CPT), a monoterpenoid indole alkaloid used as a precursor in the synthesis of chemotherapeutic drugs. O. pumila HR culture is considered as a promising alternative source of CPT, however, the knowledge about the biosynthetic pathway and regulatory mechanism is still limited. In this study, five genes that encode AP2/ERF transcription factors, namely OpERF1-OpERF5, were isolated from HR of O. pumila. Phylogenetic analysis of AP2/ERF protein sequences suggested the close evolutionary relationship of OpERF1 with stress-responsive ERF factors in Arabidopsis and of OpERF2 with ERF factors reported to regulate alkaloid production, such as ORCA3 in Catharanthus roseus, NIC2 locus ERF in tobacco, and JRE4 in tomato. We generated the transgenic HR lines of O. pumila, ERF1i and ERF2i, in which the expression of OpERF1 and OpERF2, respectively, was suppressed using RNA interference technique. The transcriptome and metabolome of these suppressed HR were analyzed for functional characterization of OpERF1 and OpERF2. Although significant changes were not observed in the metabolome, including CPT and related compounds, the suppression of OpERF2 resulted in reduced expression of genes in the 2-C-methyl-d-erythritol 4-phosphate and secologanin-strictosidine pathways, which supply a precursor, strictosidine, for CPT biosynthesis. Furthermore, while it was not conclusive for OpERF1, enrichment analysis of differentially expressed genes in the suppressed HR showed that the gene ontology terms for oxidation-reduction, presumably involved in secondary metabolite pathways, were enriched in the ERF2i downregulated gene set. These results suggest a positive role of OpERF2 in regulating specialized metabolism in O. pumila.

Modern plant metabolomics: advanced natural product gene discoveries, improved technologies, and future prospects

Plant metabolomics has matured and modern plant metabolomics has accelerated gene discoveries and the elucidation of a variety of plant natural product biosynthetic pathways. This review covers the approximate period of 2000 to 2014, and highlights specific examples of the discovery and characterization of novel genes and enzymes associated with the biosynthesis of natural products such as flavonoids, glucosinolates, terpenoids, and alkaloids. Additional examples of the integration of metabolomics with genome-based functional characterizations of plant natural products that are important to modern pharmaceutical technology are also reviewed. This article also provides a substantial review of recent technical advances in mass spectrometry imaging, nuclear magnetic resonance imaging, integrated LC-MS-SPE-NMR for metabolite identifications, and X-ray crystallography of microgram quantities for structural determinations. The review closes with a discussion on the future prospects of metabolomics related to crop species and herbal medicine.

Secondary metabolites from Rubiaceae species

This study describes some characteristics of the Rubiaceae family pertaining to the occurrence and distribution of secondary metabolites in the main genera of this family. It reports the review of phytochemical studies addressing all species of Rubiaceae, published between 1990 and 2014. Iridoids, anthraquinones, triterpenes, indole alkaloids as well as other varying alkaloid subclasses, have shown to be the most common. These compounds have been mostly isolated from the genera Uncaria, Psychotria, Hedyotis, Ophiorrhiza and Morinda. The occurrence and distribution of iridoids, alkaloids and anthraquinones point out their chemotaxonomic correlation among tribes and subfamilies. From an evolutionary point of view, Rubioideae is the most ancient subfamily, followed by Ixoroideae and finally Cinchonoideae. The chemical biosynthetic pathway, which is not so specific in Rubioideae, can explain this and large amounts of both iridoids and indole alkaloids are produced. In Ixoroideae, the most active biosysthetic pathway is the one that produces iridoids; while in Cinchonoideae, it produces indole alkaloids together with other alkaloids. The chemical biosynthetic pathway now supports this botanical conclusion.

Next Generation Sequencing and Transcriptome Analysis Predicts Biosynthetic Pathway of Sennosides from Senna (Cassia angustifolia Vahl.), a Non-Model Plant with Potent Laxative Properties

Senna (Cassia angustifolia Vahl.) is a world's natural laxative medicinal plant. Laxative properties are due to sennosides (anthraquinone glycosides) natural products. However, little genetic information is available for this species, especially concerning the biosynthetic pathways of sennosides. We present here the transcriptome sequencing of young and mature leaf tissue of Cassia angustifolia using Illumina MiSeq platform that resulted in a total of 6.34 Gb of raw nucleotide sequence. The sequence assembly resulted in 42230 and 37174 transcripts with an average length of 1119 bp and 1467 bp for young and mature leaf, respectively. The transcripts were annotated using NCBI BLAST with 'green plant database (txid 33090)', Swiss Prot, Kyoto Encylcopedia of Genes & Genomes (KEGG), Cluster of Orthologous Gene (COG) and Gene Ontology (GO). Out of the total transcripts, 40138 (95.0%) and 36349 (97.7%) from young and mature leaf, respectively, were annotated by BLASTX against green plant database of NCBI. We used InterProscan to see protein similarity at domain level, a total of 34031 (young leaf) and 32077 (mature leaf) transcripts were annotated against the Pfam domains. All transcripts from young and mature leaf were assigned to 191 KEGG pathways. There were 166 and 159 CDS, respectively, from young and mature leaf involved in metabolism of terpenoids and polyketides. Many CDS encoding enzymes leading to biosynthesis of sennosides were identified. A total of 10,763 CDS differentially expressing in both young and mature leaf libraries of which 2,343 (21.7%) CDS were up-regulated in young compared to mature leaf. Several differentially expressed genes found functionally associated with sennoside biosynthesis. CDS encoding for many CYPs and TF families were identified having probable roles in metabolism of primary as well as secondary metabolites. We developed SSR markers for molecular breeding of senna. We have identified a set of putative genes involved in various secondary metabolite pathways, especially those related to the synthesis of sennosides which will serve as an important platform for public information about gene expression, genomics, and functional genomics in senna.

Transcriptome analysis of thermogenic Arum concinnatum reveals the molecular components of floral scent production

Several plant species can generate enough heat to increase their internal floral temperature above ambient temperature. Among thermogenic plants, Arum concinnatum shows the highest respiration activity during thermogenesis. However, an overall understanding of the genes related to plant thermogenesis has not yet been achieved. In this study, we performed de novo transcriptome analysis of flower organs in A. concinnatum. The de novo transcriptome assembly represented, in total, 158,490 non-redundant transcripts, and 53,315 of those showed significant homology with known genes. To explore genes associated with thermogenesis, we filtered 1266 transcripts that showed a significant correlation between expression pattern and the temperature trend of each sample. We confirmed five putative alternative oxidase transcripts were included in filtered transcripts as expected. An enrichment analysis of the Gene Ontology terms for the filtered transcripts suggested over-representation of genes involved in 1-deoxy-d-xylulose-5-phosphate synthase (DXS) activity. The expression profiles of DXS transcripts in the methyl-d-erythritol 4-phosphate (MEP) pathway were significantly correlated with thermogenic levels. Our results suggest that the MEP pathway is the main biosynthesis route for producing scent monoterpenes. To our knowledge, this is the first report describing the candidate pathway and the key enzyme for floral scent production in thermogenic plants.

Integration of a Decrescent Transcriptome and Metabolomics Dataset of Peucedanum praeruptorum to Investigate the CYP450 and MDR Genes Involved in Coumarins Biosynthesis and Transport

Peucedanum praeruptorum Dunn is well-known traditional Chinese medicine. However, little is known in the biosynthesis and the transport mechanisms of its coumarin compounds at the molecular level. Although transcriptomic sequence is playing an increasingly significant role in gene discovery, it is not sufficient in predicting the specific function of target gene. Furthermore, there is also a huge database to be analyzed. In this study, RNA sequencing assisted transcriptome dataset and high-performance liquid chromatography (HPLC) coupled with electrospray-ionization quadrupole time-of-flight mass spectrometry (Q-TOF MS)-based metabolomics dataset of P. praeruptorum were firstly constructed for gene discovery and compound identification. Subsequently, methyl jasmonate (MeJA)-induced gene expression analysis and metabolomics analysis were conducted to narrow-down the dataset for selecting the candidate genes and the potential marker metabolites. Finally, the genes involved in coumarins biosynthesis and transport were predicted with parallel analysis of transcript and metabolic profiles. As a result, a total of 40,952 unigenes and 19 coumarin compounds were obtained. Based on the results of gene expression and metabolomics analysis, 7 cytochrome-P450 and 8 multidrug resistance transporter unigenes were selected as candidate genes and 8 marker compounds were selected as biomarkers, respectively. The parallel analysis of gene expression and metabolites accumulation indicated that the gene labeled as 23,746, 228, and 30,922 were related to the formation of the coumarin core compounds whereas 36,276 and 9533 participated in the prenylation, hydroxylation, cyclization or structural modification. Similarly, 1462, 20,815, and 15,318 participated in the transport of coumarin core compounds while 124,029 and 324,293 participated in the transport of the modified compounds. This finding suggested that integration of a decrescent transcriptome and metabolomics dataset could largely narrow down the number of gene to be investigated and significantly improve the efficiency of functional gene predication. In addition, the large amount of transcriptomic data produced from P. praeruptorum and the genes discovered in this study would provide useful information in investigating the biosynthesis and transport mechanism of coumarins.

Cloning and characterization of a norbelladine 4'-O-methyltransferase involved in the biosynthesis of the Alzheimer's drug galanthamine in Narcissus sp. aff. pseudonarcissus

Galanthamine is an Amaryllidaceae alkaloid used to treat the symptoms of Alzheimer’s disease. This compound is primarily isolated from daffodil (Narcissus spp.), snowdrop (Galanthus spp.), and summer snowflake (Leucojum aestivum). Despite its importance as a medicine, no genes involved in the biosynthetic pathway of galanthamine have been identified. This absence of genetic information on biosynthetic pathways is a limiting factor in the development of synthetic biology platforms for many important botanical medicines. The paucity of information is largely due to the limitations of traditional methods for finding biochemical pathway enzymes and genes in non-model organisms. A new bioinformatic approach using several recent technological improvements was applied to search for genes in the proposed galanthamine biosynthetic pathway, first targeting methyltransferases due to strong signature amino acid sequences in the proteins. Using Illumina sequencing, a de novo transcriptome assembly was constructed for daffodil. BLAST was used to identify sequences that contain signatures for plant O-methyltransferases in this transcriptome. The program HAYSTACK was then used to identify methyltransferases that fit a model for galanthamine biosynthesis in leaf, bulb and inflorescence tissues. One candidate gene for the methylation of norbelladine to 4′-O-methylnorbelladine in the proposed galanthamine biosynthetic pathway was identified. This methyltransferase cDNA was expressed in E. coli and the protein purified by affinity chromatography. The resulting protein was found to be a norbelladine 4′-O-methyltransferase (NpN4OMT) of the proposed galanthamine biosynthetic pathway.

De novo assembly and analysis of Cassia obtusifolia seed transcriptome to identify genes involved in the biosynthesis of active metabolites

A cDNA library generated from seeds of Cassia obtusifolia was sequenced using Illumina/Solexa platform. More than 12,968,231 high quality reads were generated, and have been deposited in NCBI SRA (SRR 1012912). A total of 40,102 unigenes (>200 bp) were obtained with an average sequence length of 681 bp by de novo assembly. About 34,089 (85%) unique sequences were annotated and 8694 of the unique sequences were assigned to specific metabolic pathways by the Kyoto Encyclopedia of Genes and Genomes. Among them, 131 unigenes, which are involved in the biosynthesis and (or) regulation of anthraquinone, carotenoid, flavonoid, and lipid, the 4 best known active metabolites, were identified from cDNA library. In addition, three lipid transfer proteins were obtained, which may contribute to the lipid molecules transporting between biological membranes. Meanwhile, 30 cytochrome P450, 12 SAM-dependent methyltransferases, and 12 UDP-glucosyltransferase unigenes were identified, which could also be responsible for the biosynthesis of active metabolites.

Genome-enabled plant metabolomics

The grand challenge currently facing metabolomics is that of comprehensitivity whilst next generation sequencing and advanced proteomics methods now allow almost complete and at least 50% coverage of their respective target molecules, metabolomics platforms at best offer coverage of just 10% of the small molecule complement of the cell. Here we discuss the use of genome sequence information as an enabling tool for peak identity and for translational metabolomics. Whilst we argue that genome information is not sufficient to compute the size of a species metabolome it is highly useful in predicting the occurrence of a wide range of common metabolites. Furthermore, we describe how via gene functional analysis in model species the identity of unknown metabolite peaks can be resolved. Taken together these examples suggest that genome sequence information is current (and likely will remain), a highly effective tool in peak elucidation in mass spectral metabolomics strategies.

Unlocking Triticeae genomics to sustainably feed the future

The tribe Triticeae includes the major crops wheat and barley. Within the last few years, the whole genomes of four Triticeae species-barley, wheat, Tausch's goatgrass (Aegilops tauschii) and wild einkorn wheat (Triticum urartu)-have been sequenced. The availability of these genomic resources for Triticeae plants and innovative analytical applications using next-generation sequencing technologies are helping to revitalize our approaches in genetic work and to accelerate improvement of the Triticeae crops. Comparative genomics and integration of genomic resources from Triticeae plants and the model grass Brachypodium distachyon are aiding the discovery of new genes and functional analyses of genes in Triticeae crops. Innovative approaches and tools such as analysis of next-generation populations, evolutionary genomics and systems approaches with mathematical modeling are new strategies that will help us discover alleles for adaptive traits to future agronomic environments. In this review, we provide an update on genomic tools for use with Triticeae plants and Brachypodium and describe emerging approaches toward crop improvements in Triticeae.

Transcriptome analysis of Japanese pear (Pyrus pyrifolia Nakai) flower buds transitioning through endodormancy

The transcriptomes of endodormant and ecodormant Japanese pear (Pyrus pyrifolia Nakai 'Kosui') flower buds were analyzed using RNA-seq technology and compared. Among de novo assembly of 114,191 unigenes, 76,995 unigenes were successfully annotated by BLAST searches against various databases. Gene Ontology (GO) enrichment analysis revealed that oxidoreductases were enriched in the molecular function category, a result consistent with previous observations of notable changes in hydrogen peroxide concentration during endodormancy release. In the GO categories related to biological process, the abundance of DNA methylation-related gene transcripts also significantly changed during endodormancy release, indicating the involvement of epigenetic regulation. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis also showed the changes in transcript abundance of genes involved in the metabolism of various phytohormones. Genes for both ABA and gibberellin biosynthesis were down-regulated, whereas the genes encoding their degradation enzymes were up-regulated during endodormancy release. In the ethylene pathway, 1-aminocyclopropane-1-carboxylate synthase (ACS), a gene encoding the rate-limiting enzyme for ethylene biosynthesis, was induced towards endodormancy release. All of these results indicated the involvement of phytohormones in endodormancy release. Furthermore, the expression of dormancy-associated MADS-box (DAM) genes was down-regulated concomitant with endodormancy release, although changes in the abundance of these gene transcripts were not as significant as those identified by transcriptome analysis. Consequently, characterization of the Japanese pear transcriptome during the transition from endormancy to ecodormancy will provide researchers with useful information for data mining and will facilitate further experiments on endodormancy especially in rosaceae fruit trees.

Phytochemical genomics--a new trend

Phytochemical genomics is a recently emerging field, which investigates the genomic basis of the synthesis and function of phytochemicals (plant metabolites), particularly based on advanced metabolomics. The chemical diversity of the model plant Arabidopsis thaliana is larger than previously expected, and the gene-to-metabolite correlations have been elucidated mostly by an integrated analysis of transcriptomes and metabolomes. For example, most genes involved in the biosynthesis of flavonoids in Arabidopsis have been characterized by this method. A similar approach has been applied to the functional genomics for production of phytochemicals in crops and medicinal plants. Great promise is seen in metabolic quantitative loci analysis in major crops such as rice and tomato, and identification of novel genes involved in the biosynthesis of bioactive specialized metabolites in medicinal plants.