De novo Sequencing and Transcriptome Analysis Reveal Key Genes Regulating Steroid Metabolism in Leaves, Roots, Adventitious Roots and Calli of Periploca sepium Bunge

Integrated transcriptomic and metabolomic analysis reveals the molecular basis of tissue-specific accumulation of bioactive steroidal alkaloids in Fritillaria unibracteata

Fritillaria unibracteata is an endangered medicinal plant whose bulb has been used as a Chinese herb to suppress cough, asthma and excessive phlegm for centuries. Steroidal alkaloids, which are synthesized via the steroid synthesis pathways, are their significant bioactive constituents. However, few studies on genes involved in steroidal alkaloid biosynthesis in F. unibracteata have been reported, mainly due to the lack of the F. unibracteata genome. In this paper, comparative transcriptomic and metabolomic analyses of four different tissues of F. unibracteata (leaves, flowers, stems, and bulbs) were performed. Imperialine, peiminine, and peimisine were among the significant bioactive compounds that were considerably abundant in bulb tissue, according to the metabolomic findings. Then, 83.60 Gb transcriptome sequencing of four different tissues was performed, of which one gene encoding phosphomevalonate kinase was directly functionally characterized to verify the accuracy of sequences obtained from the transcriptome. A total of 9217 differentially expressed unigenes (DEGs) were identified in four different tissues of F. unibracteata. GO and KEGG enrichments revealed that phenylpropanoid biosynthesis, MVA-mediated terpenoid backbone biosynthesis, and steroid biosynthesis were enriched in bulb tissue, whereas enrichment of MEP-mediated terpenoid backbone biosynthesis, photosynthesis, photosynthesis-antenna protein and carotenoid biosynthesis was observed in aerial tissues. Moreover, clustering analysis indicated that the downstream steroid biosynthesis pathway was more important in steroidal alkaloid biosynthesis compared to the upstream terpenoid backbone biosynthesis pathway. Hence, MVA-mediated biosynthesis of steroidal alkaloids was proposed, in which 15 bulb-clustered DEGs were positively correlated with a high accumulation of bioactive steroid alkaloids, further validating our proposal. In addition, 36 CYP450s showing a positive correlation with bioactive steroidal alkaloids provided candidate enzymes to catalyze the subsequent steps of steroidal alkaloid biosynthesis. In addition, the transcription factors and ABC transporters clustered in bulb tissue might be responsible for the regulation and transportation of steroidal alkaloid biosynthesis. Protein-protein interaction analysis implied a highly complex steroid alkaloid biosynthesis network in which delta (24)-sterol reductase might be one of the central catalysts. Based on the integrated transcriptome and metabolome, this current study is a first step in understanding the tissue-specific biosynthesis of steroidal alkaloids in F. unibracteata. Furthermore, key genes and regulators identified herein could facilitate metabolic engineering to improve steroidal alkaloids in F. unibracteata.

Identification of key genes involved in secondary metabolite biosynthesis in Digitalis purpurea

The medicinal plant Digitalis purpurea produces cardiac glycosides that are useful in the pharmaceutical industry. These bioactive compounds are in high demand due to ethnobotany's application to therapeutic procedures. Recent studies have investigated the role of integrative analysis of multi-omics data in understanding cellular metabolic status through systems metabolic engineering approach, as well as its application to genetically engineering metabolic pathways. In spite of numerous omics experiments, most molecular mechanisms involved in metabolic pathways biosynthesis in D. purpurea remain unclear. Using R Package Weighted Gene Co-expression Network Analysis, co-expression analysis was performed on the transcriptome and metabolome data. As a result of our study, we identified transcription factors, transcriptional regulators, protein kinases, transporters, non-coding RNAs, and hub genes that are involved in the production of secondary metabolites. Since jasmonates are involved in the biosynthesis of cardiac glycosides, the candidate genes for Scarecrow-Like Protein 14 (SCL14), Delta24-sterol reductase (DWF1), HYDRA1 (HYD1), and Jasmonate-ZIM domain3 (JAZ3) were validated under methyl jasmonate treatment (MeJA, 100 μM). Despite early induction of JAZ3, which affected downstream genes, it was dramatically suppressed after 48 hours. SCL14, which targets DWF1, and HYD1, which induces cholesterol and cardiac glycoside biosynthesis, were both promoted. The correlation between key genes and main metabolites and validation of expression patterns provide a unique insight into the biosynthesis mechanisms of cardiac glycosides in D. purpurea.

High-throughput metabolomic and transcriptomic analyses vet the potential route of cerpegin biosynthesis in two varieties of Ceropegia bulbosa Roxb

Exploration with high-throughput transcriptomics and metabolomics of two varieties of Ceropegia bulbosa identifies candidate genes, crucial metabolites and a potential cerpegin biosynthetic pathway. Ceropegia bulbosa is an important medicinal plant, used in the treatment of various ailments including diarrhea, dysentery, and syphilis. This is primarily attributed to the presence of pharmaceutically active secondary metabolites, especially cerpegin. As this plant belongs to an endemic threatened category, genomic resources are not available hampering exploration on the molecular basis of cerpegin accumulation till now. Therefore, we undertook high-throughput metabolomic and transcriptomic analyses using different tissues from two varieties namely, C. bulbosa var. bulbosa and C. bulbosa var. lushii. Metabolomic analysis revealed spatial and differential accumulation of various metabolites. We chemically synthesized and characterized the cerpegin and its derivatives by liquid chromatography tandem-mass spectrometry (LC-MS/MS). Importantly, these comparisons suggested the presence of cerpegin and 5-allyl cerpegin in all C. bulbosa tissues. Further, de novo transcriptome analysis indicated the presence of significant transcripts for secondary metabolic pathways through the Kyoto encyclopedia of genes and genomes database. Tissue-specific profiling of transcripts and metabolites showed a significant correlation, suggesting the intricate mechanism of cerpegin biosynthesis. The expression of potential candidate genes from the proposed cerpegin biosynthetic pathway was further validated by qRT-PCR and NanoString nCounter. Overall, our findings propose a potential route of cerpegin biosynthesis. Identified transcripts and metabolites have built a foundation as new molecular resources that could facilitate future research on biosynthesis, regulation, and engineering of cerpegin or other important metabolites in such non-model plants.

RNA-Seq analysis of differential gene expression in Betula luminifera xylem during the early stages of tension wood formation

Background

Betula luminifera H. Winkler, which is widely distributed in southern China, is an economically important broadleaf tree species. However, little genomic information of B. luminifera is available, and little is known about the molecular mechanisms of wood formation in this species. Meanwhile, few efforts have focused on investigating the early transcriptional changes during tension wood formation in woody plants.

Results

A reference transcriptome dataset was first generated containing 45,700 Unigenes, and 35,135 (76.9%) Unigenes were annotated by a BLAST similarity search against four public databases. Then, based on an anatomical investigation, the global gene expression changes during the early stages of tension wood formation were analyzed. Gene expression profiling showed that a total of 13,273 Unigenes were differentially regulated during the early stages of tension wood formation. Most genes involved in cellulose and lignin biosynthesis were highlighted to reveal their biological importance in tension wood formation. In addition, the transcription levels of many genes involved in the auxin response pathway were significantly changed during the early stages of tension wood formation. Furthermore, 18 TFs co-expressed with key enzymes of cellulose synthesis were identified.

Conclusions

Our results revealed the transcriptional changes associated with TW formation and identified potential key genes in the regulation of this process. These results will help to dissect the molecular mechanism of wood formation and provide key candidate genes for marker-assisted selection in B. luminifera.

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.

Comparative Transcriptome Analysis Identifies Putative Genes Involved in Steroid Biosynthesis in Euphorbia tirucalli

Phytochemical analysis of different Euphorbia tirucalli tissues revealed a contrasting tissue-specificity for the biosynthesis of euphol and β-sitosterol, which represent the two pharmaceutically active steroids in E. tirucalli. To uncover the molecular mechanism underlying this tissue-specificity for phytochemicals, a comprehensive E. tirucalli transcriptome derived from its root, stem, leaf and latex was constructed, and a total of 91,619 unigenes were generated with 51.08% being successfully annotated against the non-redundant (Nr) protein database. A comparison of the transcriptome from different tissues discovered members of unigenes in the upstream steps of sterol backbone biosynthesis leading to this tissue-specific sterol biosynthesis. Among them, the putative oxidosqualene cyclase (OSC) encoding genes involved in euphol synthesis were notably identified, and their expressions were significantly up-regulated in the latex. In addition, genome-wide differentially expressed genes (DEGs) in the different E. tirucalli tissues were identified. The cluster analysis of those DEGs showed a unique expression pattern in the latex compared with other tissues. The DEGs identified in this study would enrich the insights of sterol biosynthesis and the regulation mechanism of this latex-specificity.