Genome-wide association studies identify loci controlling specialized seed metabolites in Arabidopsis.

Plants synthesize specialized metabolites to facilitate environmental and ecological interactions. During evolution, plants diversified in their potential to synthesize these metabolites. Quantitative differences in metabolite levels of natural Arabidopsis (Arabidopsis thaliana) accessions can be employed to unravel the genetic basis for metabolic traits using genome-wide association studies (GWAS). Here, we performed metabolic GWAS on seeds of a panel of 315 A. thaliana natural accessions, including the reference genotypes C24 and Col-0, for polar and semi-polar seed metabolites using untargeted ultra-performance liquid chromatography-mass spectrometry. As a complementary approach, we performed quantitative trait locus (QTL) mapping of near-isogenic introgression lines between C24 and Col-0 for specific seed specialized metabolites. Besides common QTL between seeds and leaves, GWAS revealed seed-specific QTL for specialized metabolites, indicating differences in the genetic architecture of seeds and leaves. In seeds, aliphatic methylsulfinylalkyl and methylthioalkyl glucosinolates associated with the ALKENYL HYDROXYALKYL PRODUCING loci (GS-ALK and GS-OHP) on chromosome 4 containing alkenyl hydroxyalkyl producing 2 (AOP2) and 3 (AOP3) or with the GS-ELONG locus on chromosome 5 containing methylthioalkyl malate synthase (MAM1) and MAM3. We detected two unknown sulfur-containing compounds that were also mapped to these loci. In GWAS, some of the annotated flavonoids (kaempferol 3-O-rhamnoside-7-O-rhamnoside, quercetin 3-O-rhamnoside-7-O-rhamnoside) were mapped to transparent testa 7 (AT5G07990), encoding a cytochrome P450 75B1 monooxygenase. Three additional mass signals corresponding to quercetin-containing flavonols were mapped to UGT78D2 (AT5G17050). The association of the loci and associating metabolic features were functionally verified in knockdown mutant lines. By performing GWAS and QTL mapping, we were able to leverage variation of natural populations and parental lines to study seed specialized metabolism. The GWAS data set generated here is a high-quality resource that can be investigated in further studies.

Seed-coat protective neolignans are produced by the dirigent protein AtDP1 and the laccase AtLAC5 in Arabidopsis.

Lignans/neolignans are generally synthesized from coniferyl alcohol (CA) in the cinnamate/monolignol pathway by oxidation to generate the corresponding radicals with subsequent stereoselective dimerization aided by dirigent proteins (DIRs). Genes encoding oxidases and DIRs for neolignan biosynthesis have not been identified previously. In Arabidopsis thaliana, the DIR AtDP1/AtDIR12 plays an essential role in the 8-O-4' coupling in neolignan biosynthesis by unequivocal structural determination of the compound missing in the atdp1 mutant as a sinapoylcholine (SC)-conjugated neolignan, erythro-3-{4-[2-hydroxy-2-(4-hydroxy-3-methoxyphenyl)-1-hydroxymethylethoxy]-3,5-dimethoxyphenyl}acryloylcholine. Phylogenetic analyses showed that AtDP1/AtDIR12 belongs to the DIR-a subfamily composed of DIRs for 8-8' coupling of monolignol radicals. AtDP1/AtDIR12 is specifically expressed in outer integument 1 cells in developing seeds. As a putative oxidase for neolignan biosynthesis, we focused on AtLAC5, a laccase gene coexpressed with AtDP1/AtDIR12. In lac5 mutants, the abundance of feruloylcholine (FC)-conjugated neolignans decreased to a level comparable to those in the atdp1 mutant. In addition, SC/FC-conjugated neolignans were missing in the seeds of mutants defective in SCT/SCPL19, an enzyme that synthesizes SC. These results strongly suggest that AtDP1/AtDIR12 and AtLAC5 are involved in neolignan biosynthesis via SC/FC. A tetrazolium penetration assay showed that seed coat permeability increased in atdp1 mutants, suggesting a protective role of neolignans in A. thaliana seeds.

Retrograde sulfur flow from glucosinolates to cysteine in Arabidopsis thaliana.

Specialized (secondary) metabolic pathways in plants have long been considered one-way routes of leading primary metabolite precursors to bioactive end products. Conversely, endogenous degradation of such "end" products in plant tissues has been observed following environmental stimuli, including nutrition stress. Therefore, it is of general interest whether specialized metabolites can be reintegrated into primary metabolism to recover the invested resources, especially in the case of nitrogen- or sulfur-rich compounds. Here, we demonstrate that endogenous glucosinolates (GLs), a class of sulfur-rich plant metabolites, are exploited as a sulfur source by the reallocation of sulfur atoms to primary metabolites such as cysteine in Arabidopsis thaliana Tracer experiments using 34S- or deuterium-labeled GLs depicted the catabolic processing of GL breakdown products in which sulfur is mobilized from the thioglucoside group in GL molecules, potentially accompanied by the release of the sulfate group. Moreover, we reveal that beta-glucosidases BGLU28 and BGLU30 are the major myrosinases that initiate sulfur reallocation by hydrolyzing particular GL species, conferring sulfur deficiency tolerance in A. thaliana, especially during early development. The results delineate the physiological function of GL as a sulfur reservoir, in addition to their well-known functions as defense chemicals. Overall, our findings demonstrate the bidirectional interaction between primary and specialized metabolism, which enhances our understanding of the underlying metabolic mechanisms via which plants adapt to their environments.

Defective cytokinin signaling reprograms lipid and flavonoid gene-to-metabolite networks to mitigate high salinity in Arabidopsis.

Cytokinin (CK) in plants regulates both developmental processes and adaptation to environmental stresses. Arabidopsis histidine phosphotransfer ahp2,3,5 and type-B Arabidopsis response regulator arr1,10,12 triple mutants are almost completely defective in CK signaling, and the ahp2,3,5 mutant was reported to be salt tolerant. Here, we demonstrate that the arr1,10,12 mutant is also more tolerant to salt stress than wild-type (WT) plants. A comprehensive metabolite profiling coupled with transcriptome analysis of the ahp2,3,5 and arr1,10,12 mutants was conducted to elucidate the salt tolerance mechanisms mediated by CK signaling. Numerous primary (e.g., sugars, amino acids, and lipids) and secondary (e.g., flavonoids and sterols) metabolites accumulated in these mutants under nonsaline and saline conditions, suggesting that both prestress and poststress accumulations of stress-related metabolites contribute to improved salt tolerance in CK-signaling mutants. Specifically, the levels of sugars (e.g., trehalose and galactinol), amino acids (e.g., branched-chain amino acids and γ-aminobutyric acid), anthocyanins, sterols, and unsaturated triacylglycerols were higher in the mutant plants than in WT plants. Notably, the reprograming of flavonoid and lipid pools was highly coordinated and concomitant with the changes in transcriptional levels, indicating that these metabolic pathways are transcriptionally regulated by CK signaling. The discovery of the regulatory role of CK signaling on membrane lipid reprogramming provides a greater understanding of CK-mediated salt tolerance in plants. This knowledge will contribute to the development of salt-tolerant crops with the ability to withstand salinity as a key driver to ensure global food security in the era of climate crisis.

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.

Metabolomics and complementary techniques to investigate the plant phytochemical cosmos.

Covering: up to 2021Plants and their associated microbial communities are known to produce millions of metabolites, a majority of which are still not characterized and are speculated to possess novel bioactive properties. In addition to their role in plant physiology, these metabolites are also relevant as existing and next-generation medicine candidates. Elucidation of the plant metabolite diversity is thus valuable for the successful exploitation of natural resources for humankind. Herein, we present a comprehensive review on recent metabolomics approaches to illuminate molecular networks in plants, including chemical isolation and enzymatic production as well as the modern metabolomics approaches such as stable isotope labeling, ultrahigh-resolution mass spectrometry, metabolome imaging (spatial metabolomics), single-cell analysis, cheminformatics, and computational mass spectrometry. Mass spectrometry-based strategies to characterize plant metabolomes through metabolite identification and annotation are described in detail. We also highlight the use of phytochemical genomics to mine genes associated with specialized metabolites' biosynthesis. Understanding the metabolic diversity through biotechnological advances is fundamental to elucidate the functions of the plant-derived specialized metabolome.

HiC-Hiker: a probabilistic model to determine contig orientation in chromosome-length scaffolds with Hi-C.

MOTIVATION: De novo assembly of reference-quality genomes used to require enormously laborious tasks. In particular, it is extremely time-consuming to build genome markers for ordering assembled contigs along chromosomes; thus, they are only available for well-established model organisms. To resolve this issue, recent studies demonstrated that Hi-C could be a powerful and cost-effective means to output chromosome-length scaffolds for non-model species with no genome marker resources, because the Hi-C contact frequency between a pair of two loci can be a good estimator of their genomic distance, even if there is a large gap between them. Indeed, state-of-the-art methods such as 3D-DNA are now widely used for locating contigs in chromosomes. However, it remains challenging to reduce errors in contig orientation because shorter contigs have fewer contacts with their neighboring contigs. These orientation errors lower the accuracy of gene prediction, read alignment, and synteny block estimation in comparative genomics. RESULTS: To reduce these contig orientation errors, we propose a new algorithm, named HiC-Hiker, which has a firm grounding in probabilistic theory, rigorously models Hi-C contacts across contigs, and effectively infers the most probable orientations via the Viterbi algorithm. We compared HiC-Hiker and 3D-DNA using human and worm genome contigs generated from short reads, evaluated their performances, and observed a remarkable reduction in the contig orientation error rate from 4.3% (3D-DNA) to 1.7% (HiC-Hiker). Our algorithm can consider long-range information between distal contigs and precisely estimates Hi-C read contact probabilities among contigs, which may also be useful for determining the ordering of contigs. AVAILABILITY AND IMPLEMENTATION: HiC-Hiker is freely available at: https://github.com/ryought/hic_hiker.

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.

Comparative Metabolome and Transcriptome Analyses of Susceptible Asparagus officinalis and Resistant Wild A. kiusianus Reveal Insights into Stem Blight Disease Resistance.

Phomopsis asparagi is one of the most serious fungal pathogens, which causes stem blight disease in Asparagus officinalis (AO), adversely affecting its production worldwide. Recently, the development of novel asparagus varieties using wild Asparagus genetic resources with natural P. asparagi resistance has become a priority in Japan due to the lack of resistant commercial AO cultivars. In this study, comparative metabolome and transcriptome analyses of susceptible AO and resistant wild Asparagus kiusianus (AK) 24 and 48 h postinoculated (AOI_24 hpi, AOI_48 hpi, AKI_24 hpi and AKI_48 hpi, respectively) with P. asparagi were conducted to gain insights into metabolic and expression changes associated with AK species. Following infection, the resistant wild AK showed rapid metabolic changes with increased levels of flavonoids and steroidal saponins and decreased asparagusic acid glucose ester content, compared with the susceptible AO plants. Transcriptome data revealed a total of 21 differentially expressed genes (DEGs) as the core gene set that displayed upregulation in the resistant AK versus susceptible AO after infection with P. asparagi. Kyoto Encyclopedia of Genes and Genomes pathway analysis of these DEGs identified 11 significantly enriched pathways, including flavonoid biosynthesis and primary metabolite metabolism, in addition to plant signaling and defense-related pathways. In addition, comparative single-nucleotide polymorphism and Indel distributions in susceptible AO and resistant AK plants were evaluated using the latest AO reference genome Aspof.V1. The data generated in this study are important resources for advancing Asparagus breeding programs and for investigations of genetic linkage mapping, phylogenetic diversity and plant defense-related genes.

Metabolomics with 15N Labeling for Characterizing Missing Monoterpene Indole Alkaloids in Plants.

Monoterpene indole alkaloids (MIAs) in medicinal plants remain uncharacterized owing to their complicated structure by metabolomics using liquid chromatography-tandem mass spectrometry (LC-MS/MS) despite their pharmaceutical importance. We demonstrate an untargeted metabolome analysis with 15nitrogen (N) labeling to characterize MIAs having an indolic skeleton in the flowers, leaves, petioles, stems, and roots of Catharanthus roseus. Principal component analysis using 15N- and nonlabeled metabolome data showed that N-containing metabolites (N-metabolites) are labeled with 15N. Paring of the 15N- and nonlabeled precursor ions were performed using the criteria of retention time, difference of m/z value, and a nonlabeled product ion at m/z 144.08 that indicates an indolic skeleton. The mass shift of the m/z value of the product and precursor ions to their 15N-labeled ions identified the number of N of their ions. Finally, molecular formula of 45 MIAs was unambiguously identified using the identified N number. The alkaloid network analysis using the MS/MS similarity showed the structural commonness and uniqueness among the MIAs. Of them, antirhine was identified using an authentic standard compound. Multimetabolomics using LC-MS/MS and imaging mass spectrometry showed that antirhine accumulates considerably in the epidermis and vascular cylinder of the roots. The developed approach showed the existence of the missing MIAs. The modification of this approach will identify other MIAs that contain a hydroxylated or methoxylated indolic skeleton.

The Structural Integrity of Lignin Is Crucial for Resistance against Striga hermonthica Parasitism in Rice.

Striga species are parasitic weeds that seriously constrain the productivity of food staples, including cereals and legumes, in Sub-Saharan Africa and Asia. In eastern and central Africa, Striga spp. infest as much as 40 million hectares of smallholder farmland causing total crop failure during severe infestation. As the molecular mechanisms underlying resistance are yet to be elucidated, we undertook a comparative metabolome study using the Striga-resistant rice (Oryza sativa) cultivar 'Nipponbare' and the susceptible cultivar 'Koshihikari'. We found that a number of metabolites accumulated preferentially in the Striga-resistant cultivar upon Striga hermonthica infection. Most apparent was increased deposition of lignin, a phenylpropanoid polymer mainly composed of p-hydroxyphenyl (H), guaiacyl (G), and syringyl (S) aromatic units, around the site of interaction in Nipponbare. The increased deposition of lignin was accompanied by induction of the expression of corresponding enzyme-encoding genes in the phenylpropanoid pathway. In addition, perturbing normal lignin composition by knocking down or overexpressing the genes that regulate lignin composition, i.e. p-COUMARATE 3-HYDROXYLASE or FERULATE 5-HYDROXYLASE, enhanced susceptibility of Nipponbare to S hermonthica infection. These results demonstrate that enhanced lignin deposition and maintenance of the structural integrity of lignin polymers deposited at the infection site are crucial for postattachment resistance against S hermonthica.

Keeping the shape of plant tissue for visualizing metabolite features in segmentation and correlation analysis of imaging mass spectrometry in Asparagus officinalis.

INTRODUCTION: Matrix-assisted laser desorption/ionization-imaging mass spectrometry (MALDI-IMS) is a powerful approach for visualizing the localization of metabolites. OBJECTIVES: A method to keep the shape of plant tissue needs to be developed for MALDI-IMS. METHODS: The method was developed using transfer tape and double-sided conductive tape. RESULTS: MALDI-IMS analysis using the developed method enabled to perform segmentation and correlation analysis of mass features. CONCLUSION: This proof-of-concept study showed that rutin localizes in the epidermis, developing tissue, and protoxylem in Asparagus officinalis.

A Highly Specific Genome-Wide Association Study Integrated with Transcriptome Data Reveals the Contribution of Copy Number Variations to Specialized Metabolites in Arabidopsis thaliana Accessions.

Lineage-specific gene duplications contribute to a large variation in specialized metabolites among different plant species. There is also considerable variability in the specialized metabolites within a single plant species. However, it is unclear whether copy number variations (CNVs) derived from gene duplication events contribute to the diversity of specialized metabolites within species. We identified metabolome quantitative trait genes (mQTGs) associated with quantitative metabolite variations and examined the relationship between mQTGs and CNVs. We obtained 1,335 specialized metabolite signals from 53 worldwide A. thaliana accessions using liquid chromatography-quadrupole time-of-flight mass spectrometry. In this study, genes associated with specialized metabolites were inferred by either a generally authorized genome-wide association study (GWAS) approach or a novel analysis of the association between gene expression and metabolite accumulation. Genes qualified by both analyses are defined to be mQTGs. The integrated method enabled us to detect mQTGs with a low false positive rate (=5.71 × 10-4). We also identified 5,654 genes associated with 1,335 specialized metabolites. Of these genes, 4.4% were affected by CNVs, which was more than expected (χ2 test: P < 0.01). This result suggests that CNVs contribute to variations in specialized metabolites within a species. To assess the contribution of CNVs to adaptive evolution in A. thaliana, we examined the selective sweeps around the mQTGs. We observed that the mQTGs with CNVs tended to undergo selective sweeps. These observations imply that variations in specialized metabolites caused by CNVs contribute to the adaptive evolution of A. thaliana.

UGT79B31 is responsible for the final modification step of pollen-specific flavonoid biosynthesis in Petunia hybrida.

MAIN CONCLUSION: UGT79B31 encodes flavonol 3- O -glycoside: 2″- O -glucosyltransferase, an enzyme responsible for the terminal modification of pollen-specific flavonols in Petunia hybrida. Flavonoids are known to be involved in pollen fertility in petunia (P. hybrida) and maize (Zea mays). As a first step toward elucidating the role of flavonoids in pollen, we have identified a glycosyltransferase that is responsible for the terminal modification of petunia pollen-specific flavonoids. An in silico search of the petunia transcriptome database revealed four candidate UDP-glycosyltransferase (UGT) genes. UGT79B31 was selected for further analyses based on a correlation between the accumulation pattern of flavonol glycosides in various tissues and organs and the expression profiles of the candidate genes. Arabidopsis ugt79b6 mutants that lacked kaempferol/quercetin 3-O-glucosyl(1 â†’ 2)glucosides, were complemented by transformation with UGT79B31 cDNA under the control of Arabidopsis UGT79B6 promoter, showing that UGT79B31 functions as a flavonol 3-O-glucoside: 2″-O-glucosyltransferase in planta. Recombinant UGT79B31 protein can convert kaempferol 3-O-galactoside/glucoside to kaempferol 3-O-glucosyl(1 â†’ 2)galactoside/glucoside. UGT79B31 prefers flavonol 3-O-galactosides to the 3-O-glucosides and rarely accepted the 3-O-diglycosides as sugar acceptors. UDP-glucose was the preferred sugar donor for UGT79B31. These results indicated that UGT79B31 encodes a flavonoid 3-O-glycoside: 2″-O-glucosyltransferase. Transient expression of UGT79B31 fused to green fluorescent protein (GFP) in Nicotiana benthamiana showed that UGT79B31 protein was localized in the cytosol.

Biosynthesis of riccionidins and marchantins is regulated by R2R3-MYB transcription factors in Marchantia polymorpha.

R2R3-MYB transcription factors constitute the largest gene family among plant transcription factor families. They became largely divergent during the evolution of land plants and regulate various biological processes. The functions of R2R3-MYBs are mostly characterized in seed plants but are poorly understood in non-seed plants. Here, we examined the function of two R2R3-MYB genes of Marchantia polymorpha (Mapoly0073s0038 and Mapoly0006s0226) that are closely related to subgroup 4 of the R2R3-MYB family. We performed LC/MS/MS metabolomics, RNA-seq analysis and expression analysis in overexpressors and knockout mutants of MpMYB14 and MpMYB02. Overexpression of MpMYB14 remarkably increased the amount of riccionidins, which are specific anthocyanins in liverworts and a few flowering plants. In contrast, overexpression of MpMYB02 increased the amount of several marchantins, which are characteristic cyclic bis (bibenzyl ether) compounds in M. polymorpha and related liverworts. Knockouts of MpMYB14 and MpMYB02 abolished the accumulation of riccionidins and marchantins, respectively. The expression of MpMYB14 was up-regulated by UV-B irradiation, N deficiency, and NaCl treatment, whereas the expression of MpMYB02 was down-regulated by NaCl treatment. Our results suggest that the regulatory framework of phenolic metabolism by R2R3-MYB was already established in early land plants.

Mutations in jasmonoyl-L-isoleucine-12-hydroxylases suppress multiple JA-dependent wound responses in Arabidopsis thaliana.

Plants rapidly perceive tissue damage, such as that inflicted by insects, and activate several key defense responses. The importance of the fatty acid-derived hormone jasmonates (JA) in dictating these wound responses has been recognized for many years. However, important features pertaining to the regulation of the JA pathway are still not well understood. One key unknown is the inactivation mechanism of the JA pathway and its relationship with plant response to wounding. Arabidopsis cytochrome P450 enzymes in the CYP94 clade metabolize jasmonoyl-L-isoleucine (JA-Ile), a major metabolite of JA responsible for many biological effects attributed to the JA signaling pathway; thus, CYP94s are expected to contribute to the attenuation of JA-Ile-dependent wound responses. To directly test this, we created the double and triple knock-out mutants of three CYP94 genes, CYP94B1, CYP94B3, and CYP94C1. The mutations blocked the oxidation steps and caused JA-Ile to accumulate 3-4-fold the WT levels in the wounded leaves. Surprisingly, over accumulation of JA-Ile did not lead to a stronger wound response. On the contrary, the mutants displayed a series of symptoms reminiscent of JA-Ile deficiency, including resistance to wound-induced growth inhibition, decreased anthocyanin and trichomes, and increased susceptibility to insects. The mutants, however, responded normally to exogenous JA treatments, indicating that JA perception or signaling pathways were intact. Untargeted metabolite analyses revealed >40% reduction in wound-inducible metabolites in the mutants. These observations raise questions about the current JA signaling model and point toward a more complex model perhaps involving JA derivatives and/or feedback mechanisms. This article is part of a Special Issue entitled: Plant Lipid Biology edited by Kent D. Chapman and Ivo Feussner.

Top-down Metabolomic Approaches for Nitrogen-Containing Metabolites.

Streamlining the processes that reveal heteroatom-containing metabolites and their biosynthetic genes is essential in integrated metabolomics studies. These metabolites are especially targeted for their potential pharmaceutical activities. By using a Fourier-transform ion cyclotron resonance-mass spectrometry (FTICR-MS) instrument, we provide top-down targeted metabolomic analyses using ultrahigh-resolution liquid chromatography-mass spectrometry (LC-MS), high-resolution matrix-assisted laser desorption/ionization (MALDI), and high-resolution imaging mass spectrometry (IMS) with 15N labeling of nitrogen-containing metabolites. In this study, we efficiently extract known and unknown chemicals and spatial information from the medicinal plant Catharanthus roseus, which sources several cancer drugs. The ultrahigh-resolution LC-MS analysis showed that the molecular formula of 65 N-metabolites were identified using the petals, peduncles, leaves, petioles, stems, and roots of the non- and 15N-labeled Catharanthus plants. The high resolution MALDI analysis showed the molecular formula of 64 N-metabolites using the petals, leaves, and stems of the non- and 15N-labeled Catharanthus. The chemical assignments using molecular formulas stored in databases identified known and unknown metabolites. The comparative analyses using the assigned metabolites revealed that most of the organ-specific ions are derived from unknown N-metabolites. The high-resolution IMS analysis characterized the spatial accumulation patterns of 32 N-metabolites using the buds, leaves, stems, and roots in Catharanthus. The comparative analysis using the non- and 15N-labeled IMS data showed the same spatial accumulation patterns of a non- and 15N-labeled metabolite in the organs, showing that top-down analysis can be performed even in IMS analysis.

Identification of a flavin-containing S-oxygenating monooxygenase involved in alliin biosynthesis in garlic.

S-Alk(en)yl-l-cysteine sulfoxides are cysteine-derived secondary metabolites highly accumulated in the genus Allium. Despite pharmaceutical importance, the enzymes that contribute to the biosynthesis of S-alk-(en)yl-l-cysteine sulfoxides in Allium plants remain largely unknown. Here, we report the identification of a flavin-containing monooxygenase, AsFMO1, in garlic (Allium sativum), which is responsible for the S-oxygenation reaction in the biosynthesis of S-allyl-l-cysteine sulfoxide (alliin). Recombinant AsFMO1 protein catalyzed the stereoselective S-oxygenation of S-allyl-l-cysteine to nearly exclusively yield (RC SS )-S-allylcysteine sulfoxide, which has identical stereochemistry to the major natural form of alliin in garlic. The S-oxygenation reaction catalyzed by AsFMO1 was dependent on the presence of nicotinamide adenine dinucleotide phosphate (NADPH) and flavin adenine dinucleotide (FAD), consistent with other known flavin-containing monooxygenases. AsFMO1 preferred S-allyl-l-cysteine to γ-glutamyl-S-allyl-l-cysteine as the S-oxygenation substrate, suggesting that in garlic, the S-oxygenation of alliin biosynthetic intermediates primarily occurs after deglutamylation. The transient expression of green fluorescent protein (GFP) fusion proteins indicated that AsFMO1 is localized in the cytosol. AsFMO1 mRNA was accumulated in storage leaves of pre-emergent nearly sprouting bulbs, and in various tissues of sprouted bulbs with green foliage leaves. Taken together, our results suggest that AsFMO1 functions as an S-allyl-l-cysteine S-oxygenase, and contributes to the production of alliin both through the conversion of stored γ-glutamyl-S-allyl-l-cysteine to alliin in storage leaves during sprouting and through the de novo biosynthesis of alliin in green foliage leaves.

Metabolome-genome-wide association study dissects genetic architecture for generating natural variation in rice secondary metabolism.

Plants produce structurally diverse secondary (specialized) metabolites to increase their fitness for survival under adverse environments. Several bioactive compounds for new drugs have been identified through screening of plant extracts. In this study, genome-wide association studies (GWAS) were conducted to investigate the genetic architecture behind the natural variation of rice secondary metabolites. GWAS using the metabolome data of 175 rice accessions successfully identified 323 associations among 143 single nucleotide polymorphisms (SNPs) and 89 metabolites. The data analysis highlighted that levels of many metabolites are tightly associated with a small number of strong quantitative trait loci (QTLs). The tight association may be a mechanism generating strains with distinct metabolic composition through the crossing of two different strains. The results indicate that one plant species produces more diverse phytochemicals than previously expected, and plants still contain many useful compounds for human applications.