Metabolomics-centered mining of plant metabolic diversity and function: Past decade and future perspectives.

Plants are natural experts in organic synthesis, being able to generate large numbers of specific metabolites with widely varying structures that help them adapt to variable survival challenges. Metabolomics is a research discipline that integrates the capabilities of several types of research including analytical chemistry, statistics, and biochemistry. Its ongoing development provides strategies for gaining a systematic understanding of quantitative changes in the levels of metabolites. Metabolomics is usually performed by targeting either a specific cell, a specific tissue, or the entire organism. Considerable advances in science and technology over the last three decades have propelled us into the era of multi-omics, in which metabolomics, despite at an earlier developmental stage than genomics, transcriptomics, and proteomics, offers the distinct advantage of studying the cellular entities that have the greatest influence on end phenotype. Here, we summarize the state of the art of metabolite detection and identification, and illustrate these techniques with four case study applications: (i) comparing metabolite composition within and between species, (ii) assessing spatio-temporal metabolic changes during plant development, (iii) mining characteristic metabolites of plants in different ecological environments and upon exposure to various stresses, and (iv) assessing the performance of metabolomics as a means of functional gene identification , metabolic pathway elucidation, and metabolomics-assisted breeding through analyzing plant populations with diverse genetic variations. In addition, we highlight the prominent contributions of joint analyses of plant metabolomics and other omics datasets, including those from genomics, transcriptomics, proteomics, epigenomics, phenomics, microbiomes, and ion-omics studies. Finally, we discuss future directions and challenges exploiting metabolomics-centered approaches in understanding plant metabolic diversity.

Disease resistance conferred by components of essential chrysanthemum oil and the epigenetic regulation of OsTPS1.

The sesquiterpene alpha-bisabolol is the predominant active ingredient in essential oils that are highly valued in the cosmetics industry due to its wound healing, anti-inflammatory, and skin-soothing properties. Alpha-bisabolol was thought to be restricted to Compositae plants. Here we reveal that alpha-bisabolol is also synthesized in rice, a non-Compositae plant, where it acts as a novel sesquiterpene phytoalexin. Overexpressing the gene responsible for the biosynthesis of alpha-bisabolol, OsTPS1, conferred bacterial blight resistance in rice. Phylogenomic analyses revealed that alpha-bisabolol-synthesizing enzymes in rice and Compositae evolved independently. Further experiments demonstrated that the natural variation in the disease resistance level was associated with differential transcription of OsTPS1 due to polymorphisms in its promoter. We demonstrated that OsTPS1 was regulated at the epigenetic level by JMJ705 through the methyl jasmonate pathway. These data reveal the cross-family accumulation and regulatory mechanisms of alpha-bisabolol production.

Integrative metabolomic and transcriptomic analyses reveal the mechanisms of Tibetan hulless barley grain coloration.

Cereal grains accumulate anthocyanin during developmental process. The anthocyanin content increases at grain filling stages to develop grain coloration in cereals. However, anthocyanin biosynthesis responsible for grain coloring and its regulatory mechanisms controlled by structural and functional genes remain unclear. Therefore, this study aimed to explore the global map of metabolic changes linked to grain coloration of Tibetan hulless barley (qingke) using an integrative metabolome and transcriptome approach. Grains from three colored qingke cultivars at different developmental stages were considered for molecular and metabolic investigations. A total of 120 differentially accumulated metabolites (DAMs) and 8,327 differentially expressed genes (DEGs) were filtered. DEGs were mainly enriched in the phenylpropanoid and flavonoid pathways. The transcript levels of anthocyanin biosynthesis genes (PAL, C4H, 4CL, CHS, FLS, F3H, F3'H, DFR, ANS, GT, OMT, and MAT) significantly upregulate in colored qingke compared to the non-colored variety. During grain development and maturation, the strong correlation of HvMYC2 expression with anthocyanin contents and anthocyanin biosynthesis genes suggested it as a critical gene in anthocyanin accumulation. Further results confirmed that HvMYC2 could be activated by HvMYB and be a positive regulator of UV-B and cold tolerance in qingke. In addition, verification based on enzymatic assays indicated that six key modifier enzymes could catalyze glycosylation, malonylation, and methylation of anthocyanins, thereby dissecting the major anthocyanin modification pathway in colored qingke. Overall, our study provides global insight into anthocyanin accumulation and the mechanism underlying grain coloration in qingke.

Resistance to Powdery Mildew in Qingke Involves the Accumulation of Aromatic Phenolamides Through Jasmonate-Mediated Activation of Defense-Related Genes.

Powdery mildew (PM) leads to severe yield reduction in qingke (Hordeum vulgare L. var. nudum). Although studies have focused on identifying PM-related resistance genes, mechanistic insights into the metabolic regulation networks of resistance against PM have rarely been explored in qingke. Here, we integrated transcriptomic, proteomic and metabolomic data using PM-susceptible (G72) and PM-resistant (K69) accessions to systemically explore the mechanisms of PM resistance. The integrated results show that a rapidly transduction of jasmonic acid (JA) and (+)-7-iso-jasmonoyl-L-isoleucine (JA-Ile), and importantly, a inducing accumulation of aromatic PAs conferred qingke-specific resistance for PM stress. Functional analysis revealed that the four BAHD N-acyltransferase genes were responsible for the synthesis of aliphatic and aromatic PAs. The expression of the four genes are induced by methyl jasmonate (MeJA) and PM treatment. Co-expression network analysis shows that a histone lysine demethylase, JMJ705 gene, also induced by MeJA and PM treatment, had highly correlation with PAs biosynthesis. Chromatin immunoprecipitation (ChIP)-seq assays revealed that the level of trimethylated histone H3 lysine 27 (H3K27me3) of the four genes in MeJA and PM-treated plants was significantly reduced. Overall, our results suggest that a novel strategy for jasmonic acid signal-mediated demethylation controlling the accumulation of aromatic PAs to enhance plant immune resistance through removal of H3K27me3 and activating defense-related gene expression.

Plant metabolic gene clusters in the multi-omics era.

Secondary metabolism in plants gives rise to a vast array of small-molecule natural products. The discovery of operon-like gene clusters in plants has provided a new perspective on the evolution of specialized metabolism and the opportunity to rapidly advance the metabolic engineering of natural product production. Here, we review historical aspects of the study of plant metabolic gene clusters as well as general strategies for identifying plant metabolic gene clusters in the multi-omics era. We also emphasize the exploration of their natural variation and evolution, as well as new strategies for the prospecting of plant metabolic gene clusters and a deeper understanding of how their structure influences their function.

OsRLCK160 contributes to flavonoid accumulation and UV-B tolerance by regulating OsbZIP48 in rice.

Plants produce specialized metabolites to adapt to the ever-changing environments. Flavonoids are antioxidants essential for growth, development, and breeding with increased stress resistance in crops. However, the mechanism of the involvement of flavonoids in ultraviolet-B (UV-B) stress in rice (Oryza sativa) is largely unknown. In this study, we cloned and functionally identified a receptor-like kinase (OsRLCK160) and a bZIP transcription factor (OsbZIP48) positively regulating flavonoid accumulation through metabolite-based genome-wide association study of the flavonoid content in rice. Meanwhile, OsRLCK160 interacted with and phosphorylated OsbZIP48 to regulate the flavonoid accumulation and participate in UV-B tolerance in rice. Our study indicates the importance of applying OsRLCK160 and OsbZIP48 to advance the fundamental understanding of stable rice production and breed UV-B-tolerant rice varieties, which may contribute to breeding high-yield rice varieties.

Selection of a subspecies-specific diterpene gene cluster implicated in rice disease resistance.

Diterpenoids are the major group of antimicrobial phytoalexins in rice1,2. Here, we report the discovery of a rice diterpenoid gene cluster on chromosome 7 (DGC7) encoding the entire biosynthetic pathway to 5,10-diketo-casbene, a member of the monocyclic casbene-derived diterpenoids. We revealed that DGC7 is regulated directly by JMJ705 through methyl jasmonate-mediated epigenetic control3. Functional characterization of pathway genes revealed OsCYP71Z21 to encode a casbene C10 oxidase, sought after for the biosynthesis of an array of medicinally important diterpenoids. We further show that DGC7 arose relatively recently in the Oryza genus, and that it was partly formed in Oryza rufipogon and positively selected for in japonica during domestication. Casbene-synthesizing enzymes that are functionally equivalent to OsTPS28 are present in several species of Euphorbiaceae but gene tree analysis shows that these and other casbene-modifying enzymes have evolved independently. As such, combining casbene-modifying enzymes from these different families of plants may prove effective in producing a diverse array of bioactive diterpenoid natural products.

Cytochrome P450 CYP716A254 catalyzes the formation of oleanolic acid from ß-amyrin during oleanane-type triterpenoid saponins biosynthesis in Anemone flaccida.

Anemone flaccida Fr. Shmidt (Ranunculaceae), known as 'Di Wu' in China, is a perennial herb which has long been used to treat arthritis. The rhizome of A. flaccida contains pharmacologically active components i.e. oleanane-type triterpenoid saponins. Oleanolic acid is natural triterpenoid in plants with diverse biological activities. The biosynthesis of oleanolic acid involves cyclization of 2,3-oxidosqualene to the oleanane-type triterpenoid skeleton, followed by a series of oxidation reactions catalyzed by cytochrome P450 monooxygenase (CYP450). Previously, we identified four possible cytochrome P450 genes belonging to CYP716A subfamily from the transcriptome of A. flaccida. In this study, we identified one of those genes "CYP716A254" encoding a cytochrome P450 monooxygenase from A. flaccida that catalyzes the conversion of the ß-amyrin into oleanolic acid. The heterologous expression of CYP716A254 in yeast resulted in oxidation of ß-amyrin at the C-18 position to oleanolic acid production. These results provide an important basis for further studies of oleanane-type triterpenoid saponins synthesis in A. flaccida.

Comprehensive Analysis of the Triterpenoid Saponins Biosynthetic Pathway in Anemone flaccida by Transcriptome and Proteome Profiling.

BACKGROUND: Anemone flaccida Fr. Shmidt (Ranunculaceae), commonly known as 'Di Wu' in China, is a perennial herb with limited distribution. The rhizome of A. flaccida has long been used to treat arthritis as a tradition in China. Studies disclosed that the plant contains a rich source of triterpenoid saponins. However, little is known about triterpenoid saponins biosynthesis in A. flaccida. RESULTS: In this study, we conducted the tandem transcriptome and proteome profiling of a non-model medicinal plant, A. flaccida. Using Illumina HiSeq 2000 sequencing and iTRAQ technique, a total of 46,962 high-quality unigenes were obtained with an average sequence length of 1,310 bp, along with 1473 unique proteins from A. flaccida. Among the A. flaccida transcripts, 36,617 (77.97%) showed significant similarity (E-value < 1e (-5)) to the known proteins in the public database. Of the total 46,962 unigenes, 36,617 open reading frame (ORFs) were predicted. By the fragments per kilobases per million reads (FPKM) statistics, 14,004 isoforms/unigenes were found to be upregulated, and 14,090 isoforms/unigenes were down-regulated in the rhizomes as compared to those in the leaves. Based on the bioinformatics analysis, all possible enzymes involved in the triterpenoid saponins biosynthetic pathway of A. flaccida were identified, including cytosolic mevalonate pathway (MVA) and the plastidial methylerythritol pathway (MEP). Additionally, a total of 126 putative cytochrome P450 (CYP450) and 32 putative UDP glycosyltransferases were selected as the candidates of triterpenoid saponins modifiers. Among them, four of them were annotated as the gene of CYP716A subfamily, the key enzyme in the oleanane-type triterpenoid saponins biosynthetic pathway. Furthermore, based on RNA-Seq and proteome analysis, as well as quantitative RT-PCR verification, the expression level of gene and protein committed to triterpenoids biosynthesis in the leaf versus the rhizome was compared. CONCLUSION: A combination of the de novo transcriptome and proteome profiling based on the Illumina HiSeq 2000 sequencing platform and iTRAQ technique was shown to be a powerful method for the discovery of candidate genes, which encoded enzymes that were responsible for the biosynthesis of novel secondary metabolites in a non-model plant. The transcriptome data of our study provides a very important resource for the understanding of the triterpenoid saponins biosynthesis of A. flaccida.

The Transcript Profile of a Traditional Chinese Medicine, Atractylodes lancea, Revealing Its Sesquiterpenoid Biosynthesis of the Major Active Components.

Atractylodes lancea (Thunb.) DC., named "Cangzhu" in China, which belongs to the Asteraceae family. In some countries of Southeast Asia (China, Thailand, Korea, Japan etc.) its rhizome, commonly called rhizoma atractylodis, is used to treat many diseases as it contains a variety of sesquiterpenoids and other components of medicinal importance. Despite its medicinal value, the information of the sesquiterpenoid biosynthesis is largely unknown. In this study, we investigated the transcriptome analysis of different tissues of non-model plant A. lancea by using short read sequencing technology (Illumina). We found 62,352 high quality unigenes with an average sequence length of 913 bp in the transcripts of A. Lancea. Among these, 43,049 (69.04%), 30,264 (48.53%), 26,233 (42.07%), 17,881 (28.67%) and 29,057(46.60%) unigenes showed significant similarity (E-value<1e(-5)) to known proteins in Nr, KEGG, SWISS-PROT, GO, and COG databases, respectively. Of the total 62,352 unigenes, 43,049 (Nr Database) open reading frames were predicted. On the basis of different bioinformatics tools we identify all the enzymes that take part in the terpenoid biosynthesis as well as five different known sesquiterpenoids via cytosolic mevalonic acid (MVA) pathway and plastidal methylerythritol phosphate (MEP) pathways. In our study, 6, 864 Simple Sequence Repeats (SSRs) were also found as great potential markers in A. lancea. This transcriptomic resource of A. lancea provides a great contribution in advancement of research for this specific medicinal plant and more specifically for the gene mining of different classes of terpenoids and other chemical compounds that have medicinal as well as economic importance.