High-throughput RNA sequencing analysis of Mallotus japonicus revealed novel polerovirus and amalgavirus.

High-throughput RNA sequencing (RNA-seq) analysis of samples from Mallotus japonicus, a traditional medicinal plant, yielded two novel RNA viruses tentatively named Mallotus japonicus virus A (MjVA) and Mallotus japonicus virus B (MjVB). The MjVA and MjVB genomes encode proteins showing amino acid sequence similarities to those of poleroviruses (the genus Polerovirus, the family Solemoviridae) and amalgaviruses (the genus Amalgavirus, the family Amalgaviridae), respectively. The MjVA genome contains seven highly overlapping open reading frames, which are translated to seven proteins through various translational mechanisms, including -1 programmed ribosomal frameshifting (PRF) at the slippery motif GGGAAAC, non-AUG translational initiation, and stop codon readthrough. The MjVB genome encodes two proteins; one of which is translated by +1 PRF mechanism at the slippery motif UUUCGN. The abundance analysis of virus-derived RNA fragments revealed that MjVA is highly concentrated in plant parts with well-developed phloem tissues as previously demonstrated in other poleroviruses, which are transmitted by phloem feeders, such as aphids. MjVB, an amalgavirus generally transmitted by seeds, is distributed in all samples at low concentrations. Thus, this study demonstrates the effectiveness and usefulness of RNA-seq analysis of plant samples for the identification of novel RNA viruses and analysis of their tissue distribution. Keywords: Polerovirus; Amalgavirus; Mallotus japonicus; RNA virus; viral genome; programmed ribosomal frameshifting.

Daidzein Hydroxylation by CYP81E63 Is Involved in the Biosynthesis of Miroestrol in Pueraria mirifica.

White Kwao Krua (Pueraria candollei var. mirifica), a Thai medicinal plant, is a rich source of phytoestrogens, especially isoflavonoids and chromenes. These phytoestrogens are well known; however, their biosynthetic genes remain largely uncharacterized. Cytochrome P450 (P450) is a large protein family that plays a crucial role in the biosynthesis of various compounds in plants, including phytoestrogens. Thus, we focused on P450s involved in the isoflavone hydroxylation that potentially participates in the biosynthesis of miroestrol. Three candidate P450s were isolated from the transcriptome libraries by considering the phylogenetic and expression data of each tissue of P. mirifica. The candidate P450s were functionally characterized both in vitro and in planta. Accordingly, the yeast microsome harboring PmCYP81E63 regiospecifically exhibited either 2' or 3' daidzein hydroxylation and genistein hydroxylation. Based on in silico calculation, PmCYP81E63 had higher binding energy with daidzein than with genistein, which supported the in vitro result of the isoflavone specificity. To confirm in planta function, the candidate P450s were then transiently co-expressed with isoflavone-related genes in Nicotiana benthamiana. Despite no daidzein in the infiltrated N. benthamiana leaves, genistein and hydroxygenistein biosynthesis were detectable by liquid Chromatography with tandem mass spectrometry (LC-MS/MS). Additionally, we demonstrated that PmCYP81E63 interacted with several enzymes related to isoflavone biosynthesis using bimolecular fluorescence complementation studies and a yeast two-hybrid analysis, suggesting a scheme of metabolon formation in the pathway. Our findings provide compelling evidence regarding the involvement of PmCYP81E63 in the early step of the proposed miroestrol biosynthesis in P. mirifica.

Chromosome-scale genome assembly of Glycyrrhiza uralensis revealed metabolic gene cluster centred specialized metabolites biosynthesis.

A high-quality genome assembly is imperative to explore the evolutionary basis of characteristic attributes that define chemotype and provide essential resources for a molecular breeding strategy for enhanced production of medicinal metabolites. Here, using single-molecule high-fidelity (HiFi) sequencing reads, we report chromosome-scale genome assembly for Chinese licorice (Glycyrrhiza uralensis), a widely used herbal and natural medicine. The entire genome assembly was achieved in eight chromosomes, with contig and scaffold N50 as 36.02 and 60.2 Mb, respectively. With only 17 assembly gaps and half of the chromosomes having no or one assembly gap, the presented genome assembly is among the best plant genomes to date. Our results showed an advantage of using highly accurate long-read HiFi sequencing data for assembling a highly heterozygous genome including its complexed repeat content. Additionally, our analysis revealed that G. uralensis experienced a recent whole-genome duplication at approximately 59.02 million years ago post a gamma (γ) whole-genome triplication event, which contributed to its present chemotype features. The metabolic gene cluster analysis identified 355 gene clusters, which included the entire biosynthesis pathway of glycyrrhizin. The genome assembly and its annotations provide an essential resource for licorice improvement through molecular breeding and the discovery of valuable genes for engineering bioactive components and understanding the evolution of specialized metabolites biosynthesis.

Identification of a regiospecific S-oxygenase for the production of marasmin in traditional medicinal plant Tulbaghia violacea.

Marasmin [S-(methylthiomethyl)-L-cysteine-4-oxide] is a pharmaceutically valuable sulfur-containing compound produced by the traditional medicinal plant, Tulbaghia violacea. Here, we report the identification of an S-oxygenase, TvMAS1, that produces marasmin from its corresponding sulfide, S-(methylthiomethyl)-L-cysteine. The amino acid sequence of TvMAS1 showed high sequence similarity to known flavin-containing S-oxygenating monooxygenases in plants. Recombinant TvMAS1 catalyzed regiospecific S-oxygenation at S4 of S-(methylthiomethyl)-L-cysteine to yield marasmin, with an apparent K m value of 0.55 mM. TvMAS1 mRNA accumulated with S-(methylthiomethyl)-L-cysteine and marasmin in various organs of T. violacea. Our findings suggest that TvMAS1 catalyzes the S-oxygenation reaction during the last step of marasmin biosynthesis in T. violacea.

The ability of callus tissues induced from three Allium plants to accumulate health-beneficial natural products, S-alk(en)ylcysteine sulfoxides.

S-Alk(en)ylcysteine sulfoxides (CSOs), such as methiin, alliin, and isoalliin, are health-beneficial natural products biosynthesized in the genus Allium. Here, we report the induction of multiple callus tissue lines from three Allium vegetables, onion (A. cepa), Welsh onion (A. fistulosum), and Chinese chive (A. tuberosum), and their ability to accumulate CSOs. Callus tissues were initiated and maintained in the presence of picloram and 2-isopentenyladenine as auxin and cytokinin, respectively. For all plant species tested, the callus tissues almost exclusively accumulated methiin as CSO, while the intact plants contained a substantial amount of isoalliin together with methiin. These results suggest that the cellular developmental conditions and the regulatory mechanisms required for the biosynthesis of methiin are different from those of alliin and isoalliin. The methiin content in the callus tissues of onion and Welsh onion was much higher compared to that in the intact plants, and its cellular concentration could be estimated as 1.9-21.7 mM. The activity of alliinase that degrades CSOs in the callus tissues was much lower than that of the intact plants for onion and Welsh onion, but at similar levels as in the intact plants for Chinese chive. Our findings that the callus tissues of onion and Welsh onion showed high methiin content and low alliinase activity highlighted their potential as a plant-based system for methiin production.

Discovery of DS-3801b, a non-macrolide GPR38 agonist with N-methylanilide structure.

Motilin is a 22-amino-acid gastrointestinal (GI) hormone and is involved in the regulation of GI motility through binding to GPR38, the motilin receptor which is expressed on smooth muscle cells in the GI tract. Therefore, GPR38 agonists are expected to be novel gastrointestinal prokinetic agents for the treatment of functional gastrointestinal disorders such as gastroparesis and chronic constipation. We identified a series of N-methylanilide derivatives as novel non-macrolide GPR38 agonists. Among them, 12 di-l-tartrate (DS-3801b) was selected as a clinical candidate for further evaluation.

Gene-Metabolite Network Analysis Revealed Tissue-Specific Accumulation of Therapeutic Metabolites in Mallotus japonicus.

Mallotus japonicus is a valuable traditional medicinal plant in East Asia for applications as a gastrointestinal drug. However, the molecular components involved in the biosynthesis of bioactive metabolites have not yet been explored, primarily due to a lack of omics resources. In this study, we established metabolome and transcriptome resources for M. japonicus to capture the diverse metabolite constituents and active transcripts involved in its biosynthesis and regulation. A combination of untargeted metabolite profiling with data-dependent metabolite fragmentation and metabolite annotation through manual curation and feature-based molecular networking established an overall metabospace of M. japonicus represented by 2129 metabolite features. M. japonicus de novo transcriptome assembly showed 96.9% transcriptome completeness, representing 226,250 active transcripts across seven tissues. We identified specialized metabolites biosynthesis in a tissue-specific manner, with a strong correlation between transcripts expression and metabolite accumulations in M. japonicus. The correlation- and network-based integration of metabolome and transcriptome datasets identified candidate genes involved in the biosynthesis of key specialized metabolites of M. japonicus. We further used phylogenetic analysis to identify 13 C-glycosyltransferases and 11 methyltransferases coding candidate genes involved in the biosynthesis of medicinally important bergenin. This study provides comprehensive, high-quality multi-omics resources to further investigate biological properties of specialized metabolites biosynthesis in M. japonicus.

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.

Multiomics-based characterization of specialized metabolites biosynthesis in Cornus Officinalis.

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

DFT Study on the Biosynthesis of Preasperterpenoid A: Role of Secondary Carbocations in the Carbocation Cascade.

Preasperterpenoid A, featuring a 5/7/(3)6/5 pentacyclic structure, is a C25 sesterterpenoid produced by Penicillium verruculosum. The results of density functional calculations on putative biosynthetic carbocation cyclization/rearrangements leading to preasperterpenoid A revealed a highly concerted four-step cyclization mechanism. Interestingly, two secondary carbocation structures were obtained as minima, but appeared almost as shoulders in the energy profile, and may represent essentially transient structures during the highly concerted reaction.

Transcriptome analysis of Pueraria candollei var. mirifica for gene discovery in the biosyntheses of isoflavones and miroestrol.

BACKGROUND: Pueraria candollei var. mirifica, a Thai medicinal plant used traditionally as a rejuvenating herb, is known as a rich source of phytoestrogens, including isoflavonoids and the highly estrogenic miroestrol and deoxymiroestrol. Although these active constituents in P. candollei var. mirifica have been known for some time, actual knowledge regarding their biosynthetic genes remains unknown. RESULTS: Miroestrol biosynthesis was reconsidered and the most plausible mechanism starting from the isoflavonoid daidzein was proposed. A de novo transcriptome analysis was conducted using combined P. candollei var. mirifica tissues of young leaves, mature leaves, tuberous cortices, and cortex-excised tubers. A total of 166,923 contigs was assembled for functional annotation using protein databases and as a library for identification of genes that are potentially involved in the biosynthesis of isoflavonoids and miroestrol. Twenty-one differentially expressed genes from four separate libraries were identified as candidates involved in these biosynthetic pathways, and their respective expressions were validated by quantitative real-time reverse transcription polymerase chain reaction. Notably, isoflavonoid and miroestrol profiling generated by LC-MS/MS was positively correlated with expression levels of isoflavonoid biosynthetic genes across the four types of tissues. Moreover, we identified R2R3 MYB transcription factors that may be involved in the regulation of isoflavonoid biosynthesis in P. candollei var. mirifica. To confirm the function of a key-isoflavone biosynthetic gene, P. candollei var. mirifica isoflavone synthase identified in our library was transiently co-expressed with an Arabidopsis MYB12 transcription factor (AtMYB12) in Nicotiana benthamiana leaves. Remarkably, the combined expression of these proteins led to the production of the isoflavone genistein. CONCLUSIONS: Our results provide compelling evidence regarding the integration of transcriptome and metabolome as a powerful tool for identifying biosynthetic genes and transcription factors possibly involved in the isoflavonoid and miroestrol biosyntheses in P. candollei var. mirifica.

Inherent atomic mobility changes in carbocation intermediates during the sesterterpene cyclization cascade.

We previously showed that the regio- and stereoselectivity in terpene-forming reactions are determined by the conformations of the carbocation intermediates, which reflect the initial conformation of the substrate, geranylfarnesyl diphosphate (GFPP). However, it remains unclear how the initial conformation of GFPP is controlled, and which part(s) of the GFPP molecule are important for its fixation inside the substrate-binding pocket. Here, we present the first detailed analysis of the inherent atomic mobility in carbocation intermediates during sesterterpene biosynthesis. We identified two methyl groups as the least mobile of all the carbons of the carbocation intermediates in the first half of the cyclization cascade. Our analysis suggests that these two methyl groups are critical for the preorganization of GFPP in the biosynthetic pathways leading to sesterfisherol and quiannulatene.

The complexity of intercellular localisation of alkaloids revealed by single-cell metabolomics.

Catharanthus roseus is a medicinal plant well known for producing bioactive compounds such as vinblastine and vincristine, which are classified as terpenoid indole alkaloids (TIAs). Although the leaves of this plant are the main source of these antitumour drugs, much remains unknown on how TIAs are biosynthesised from a central precursor, strictosidine, to various TIAs in planta. Here, we have succeeded in showing, for the first time in leaf tissue of C. roseus, cell-specific TIAs localisation and accumulation with 10 µm spatial resolution Imaging mass spectrometry (Imaging MS) and live single-cell mass spectrometry (single-cell MS). These metabolomic studies revealed that most TIA precursors (iridoids) are localised in the epidermal cells, but major TIAs including serpentine and vindoline are localised instead in idioblast cells. Interestingly, the central TIA intermediate strictosidine also accumulates in both epidermal and idioblast cells of C. roseus. Moreover, we also found that vindoline accumulation increases in laticifer cells as the leaf expands. These discoveries highlight the complexity of intercellular localisation in plant specialised metabolism.

Metabolic diversification of nitrogen-containing metabolites by the expression of a heterologous lysine decarboxylase gene in Arabidopsis.

Lysine decarboxylase converts l-lysine to cadaverine as a branching point for the biosynthesis of plant Lys-derived alkaloids. Although cadaverine contributes towards the biosynthesis of Lys-derived alkaloids, its catabolism, including metabolic intermediates and the enzymes involved, is not known. Here, we generated transgenic Arabidopsis lines by expressing an exogenous lysine/ornithine decarboxylase gene from Lupinus angustifolius (La-L/ODC) and identified cadaverine-derived metabolites as the products of the emerged biosynthetic pathway. Through untargeted metabolic profiling, we observed the upregulation of polyamine metabolism, phenylpropanoid biosynthesis and the biosynthesis of several Lys-derived alkaloids in the transgenic lines. Moreover, we found several cadaverine-derived metabolites specifically detected in the transgenic lines compared with the non-transformed control. Among these, three specific metabolites were identified and confirmed as 5-aminopentanal, 5-aminopentanoate and δ-valerolactam. Cadaverine catabolism in a representative transgenic line (DC29) was traced by feeding stable isotope-labeled [α-15 N]- or [ε-15 N]-l-lysine. Our results show similar 15 N incorporation ratios from both isotopomers for the specific metabolite features identified, indicating that these metabolites were synthesized via the symmetric structure of cadaverine. We propose biosynthetic pathways for the metabolites on the basis of metabolite chemistry and enzymes known or identified through catalyzing specific biochemical reactions in this study. Our study shows that this pool of enzymes with promiscuous activities is the driving force for metabolite diversification in plants. Thus, this study not only provides valuable information for understanding the catabolic mechanism of cadaverine but also demonstrates that cadaverine accumulation is one of the factors to expand plant chemodiversity, which may lead to the emergence of Lys-derived alkaloid biosynthesis.

Acceleration of Mechanistic Investigation of Plant Secondary Metabolism Based on Computational Chemistry.

This review describes the application of computational chemistry to plant secondary metabolism, focusing on the biosynthetic mechanisms of terpene/terpenoid, alkaloid, flavonoid, and lignin as representative examples. Through these biosynthetic studies, we exhibit several computational methods, including density functional theory (DFT) calculations, theozyme calculation, docking simulation, molecular dynamics (MD) simulation, and quantum mechanics/molecular mechanics (QM/MM) calculation. This review demonstrates how modern computational chemistry can be employed as an effective tool for revealing biosynthetic mechanisms and the potential of computational chemistry-for example, elucidating how enzymes regulate regio- and stereoselectivity, finding the key catalytic residue of an enzyme, and assessing the viability of hypothetical pathways. Furthermore, insights for the next research objective involving application of computational chemistry to plant secondary metabolism are provided herein. This review will be helpful for plant scientists who are not well versed with computational chemistry.

Publisher Correction: A cheminformatics approach to characterize metabolomes in stable-isotope-labeled organisms.

In the originally published Supplementary Information for this paper, the files presented as Supplementary Tables 3, 4, and 7 were duplicates of Supplementary Tables 5, 6, and 9, respectively. All Supplementary Table files are now correct online.

A cheminformatics approach to characterize metabolomes in stable-isotope-labeled organisms.

We report a computational approach (implemented in MS-DIAL 3.0; http://prime.psc.riken.jp/) for metabolite structure characterization using fully 13C-labeled and non-labeled plants and LC-MS/MS. Our approach facilitates carbon number determination and metabolite classification for unknown molecules. Applying our method to 31 tissues from 12 plant species, we assigned 1,092 structures and 344 formulae to 3,604 carbon-determined metabolite ions, 69 of which were found to represent structures currently not listed in metabolome databases.

Metastatic thymic tumor arising from spermatic cord leiomyosarcoma.

We report the first known case of a metastatic thymic tumor arising from spermatic cord leiomyosarcoma, occurring in a 35-year-old man. He had undergone an orchiectomy 24 months previously and a surgical excision of a subcutaneous metastasis 4 months prior to his current presentation. Computed tomography revealed a 1.5-cm, round-shaped anterior mediastinal mass. A thymectomy was performed and the diagnosis of metastatic leiomyosarcoma was made.

Identification of potential genes involved in triterpenoid saponins biosynthesis in Gleditsia sinensis by transcriptome and metabolome analyses.

Gleditsia sinensis is widely used as a medicinal plant in Asia, especially in China. Triterpenes, alkaloids, and sterols were isolated from Gleditsia species. Among them, triterpenoid saponins are very important metabolites owing to their various pharmacological activities. However, the triterpenoid saponin biosynthesis pathway has not been well characterized. In the present study, we performed de novo transcriptome assembly for 14.3 Gbps of clean reads sequenced from nine tissues of G. sinensis. The results showed that 81,511 unique transcripts (unitranscripts) (47,855 unigenes) were constructed, of which 31,717 unigenes were annotated with Gene Ontology and EC numbers by Blast2GO against the NCBI-nr protein database. We also analyzed the metabolite contents in the same nine tissues by LS-MS/MS, and saponins including gleditsioside I were found in fruit at higher levels. Many of the genes with tissue-specific expression in fruit are involved in the flavonoid biosynthesis pathway, and many of those have UDP-glucosyltransferase (UGT) activity. We constructed a saponin biosynthesis pathway and identified two key enzyme families in the triterpenoid saponin biosynthesis pathway, cytochrome P450 and UDP-glucosyltransferase, that are encoded by 37 unigenes and 77 unigenes, respectively. CYP72A, CYP716A, and CYP88D, which are known as key enzymes for saponin biosynthesis, were also identified among the P450s. Our results provide insight into the secondary metabolite biosynthesis and serve as important resources for future research and cultivation of G. sinensis.