Abstract
Plant specialized metabolism emerged from the land colonization by ancient plants, becoming diversified along with plant evolution. To date, more than 1 million metabolites have been predicted to exist in the plant kingdom, and their metabolic processes have been revealed on the molecular level. Previous studies have reported that rates of evolution are greater for genes involved in plant specialized metabolism than in primary metabolism. This perspective introduces topics on the enigmatic molecular evolution of some plant specialized metabolic processes. Two transferase families, BAHD acyltransferases and aromatic prenyltransferases, which are involved in the biosynthesis of paclitaxel and meroterpenes, respectively, have shown apparent expansion. The latter family has been shown to beinvolved in the biosynthesis of a variety of aromatic substances, including prenylated coumarins in citrus plants and shikonin in Lithospermum erythrorhizon. These genes have evolved in the development of each special subfamily within the plant lineage. The broadness of substrate specificity and the exon-intron structure of their genes may provide hints to explain the evolutionary process underlying chemodiversity in plants.
Introduction
Since land plant colonization 500 million years ago, plant specialized metabolic processes have expanded considerably, resulting in the development of diverse traits within the plant kingdom (Weng et al., 2012). The chemical diversity of those natural products provides various metabolites beneficial for human life, including compounds associated with flavor, color, taste, and medicine. A comparative genome analysis strongly suggested that gene duplications played a major role in the evolution of divergent metabolic pathways (Fani and Fondi, 2009). The increase in the number of gene copies may have allowed promiscuous diversity of the encoded enzymes, resulting in the synthesis of new metabolites and providing organismal fitness that enhances the establishment of biosynthetic pathways in the plant lineage. The expansion of plant specialized metabolism has been observed in the genome of Selaginella moellendorffii, a plant that diverged shortly after the establishment of vascular tissues in plant evolution (Banks et al., 2011). A representative example of these expanded gene families is cytochrome P450-dependent monooxygenases, which constitute 1% of the predicted proteome in Selaginella. The genome of liverwort, Marchantia polymorpha, also encodes many terpenoid biosynthetic enzymes sharing a common isoprenoid pathway, a derivative designated taxadiene for the synthesis of plant hormones like gibberellin (Bowman et al., 2017). In Physcomitrella patens, a diterpene ent-kaurene is converted to gibberellin-type diterpenes, which act as regulators of protonema differentiation (Hayashi et al., 2010).
Species of the gymnosperm Taxus synthesize unique diterpene compounds called “taxoids,” which include an important anticancer drug, paclitaxel, a derivative designated taxadiene (Guerra-Bubb et al., 2012). Over 350 taxoid compounds were identified by 1999, with these compounds having variable side residues at the C1, C2, C4, C5, C7, C9, C10, C13, and C14 positions of the core taxadiene skeleton (Baloglu and Kingston, 1999). Except for a partial biosynthetic route (Croteau et al., 2006), knowledge about the biosynthetic pathway of taxoids that contribute to the chemodiversity in Taxus is limited.
Because of their fine-tuned genome data resources, angiosperm species provide good model systems to study molecular mechanisms underlying the chemodiversity of plant metabolites (Kroymann, 2011). For example, meroterpenes, including furanocoumarin derivatives (Bourgaud et al., 2006) and shikonin derivatives that are lipophilic red naphthoquinone (Yazaki, 2017), are specialized metabolites synthesized through branched routes from a metabolic pathway common to the general phenylpropanoid and isoprenoid biosynthetic pathways (Yazaki et al., 2017). The term “primary metabolism” indicates processes required to sustain life, such as energy acquisition from glucose. These processes include, for example, the biosynthesis of ubiquinone, a component of the respiratory chain in mitochondria. The biosynthesis of shikonin derivatives involves steps common to those involved in ubiquinone biosynthesis. To avoid confusion in distinguishing between primary and specialized (secondary) metabolism, this article uses the term “common metabolism” rather than “primary metabolism” to indicate biosynthetic pathways conserved in a broad variety of organisms.
This perspective focuses on two enzyme families as examples of molecular evolutionary events: the aromatic substrate prenyltransferase family, which plays a key role in the diversification of phenolics, and the BAHD (BEAT-AHCT-HCBT-DAT; initials of representative members) acyltransferase family, which is responsible for the derivatization of a core metabolite.
Evolution of the Citrus Prenyltransferase Gene Family
Among prenyltransferase superfamily including prenyl chain elogation enzymes, aromatic prenyltransferases represent a family responsible for the prenylation of aromatic substances. An aromatic prenyltransferase of Citrus limon, ClPT1, is responsible for the biosynthesis of 8-geranylumbelliferone, a coumarin derivative of a plant specialized metabolite (Munakata et al., 2014). The chemical diversity of coumarin derivatives is greatly increased by the involvement of aromatic prenyltransferases, which have been identified in many plant lineages during the last decade (Karamat et al., 2014; Munakata et al., 2014). Phylogenetic analysis has suggested that the diverse prenyltransferases developed independently in each plant family rather than developing from a common ancestor within the prenyltransferase gene family (Munakata et al., 2016). The plant prenyltransferase gene family contains conserved subfamilies responsible for the ubiquinone, plastoquinone, and vitamin E biosynthesis pathways (Li, 2016).
An outline of the evolutionary development of plant aromatic prenyltransferases in Citrus species was revealed by a phylogenetic analysis of previously characterized prenyltransferases and prenyltransferases of the model species P. patens, S. moellendorffii, Arabidopsis thaliana, Glycine max, and Lithospermum erythrorhizon (see below), in addition to Citrus sinensis (Figure 1A). Phylogenetically, these intrinsic membrane proteins can begrouped into three major subfamilies, i.e., those involved in the biosynthesis of vitamin E, plastoquinone, and ubiquinone (shown as yellow and gray backgrounds and as the black triangle, respectively in Figure 1A, with the black triangle expanded in Figure 1B). The biochemical functions of AtVTE2-1 (Savidge et al., 2002), AtVTE2–2 (Venkatesh et al., 2006), and OsPPT1 (Ohara et al., 2006) have been described. As expected from their fundamental roles, all model plant species had one or more proteins in each subfamily. In contrast, a search of the C. sinensis database revealed nine prenyltransferase-like proteins, forming a Citrus-specific subfamily within the vitamin E clade (shown in red in Figure 1A). A similar result was obtained by searching Citrus clementina genome sequences. These results suggest that Citrus species have developed a unique, expanded gene subfamily for specialized metabolism, with ClPT1 being biochemically characterized. This analysis also identified a similar unique subfamily expansion in G. max (shown in blue in Figure 1A). The first flavonoid-specific prenyltransferase SfN8DT1 from a legume species Sophora flavescens (Sasaki et al., 2008) is in this group, suggesting that flavonoid prenyltransferases in soybeans were derived from a vitamin E biosynthetic enzyme. Other, later detected flavonoid prenyltransferases were all classified in this subgroup (Akashi et al., 2009; Yoneyama et al., 2016). Most of these enzymes involved in specialized metabolism show strict substrate specificity in relation to a particular prenyl diphosphate.
Figure 1
Prenyltransferases involved in common metabolism show broad specificity in relation to substrates of different side chain lengths; i.e., they accept various prenyl diphosphates of different chain lengths (Sadre et al., 2010). For example, the ubiquinone biosynthesis pathway in rice can bemodified by introducing a decaprenyl diphosphate synthase, resulting in the production of non-native UQ10 rather than native UQ9 (Ohara et al., 2006; Takahashi et al., 2006). These expanded gene families and the broad substrate specificity of prenyltransferases may provide the opportunity for neo-functionalization of new enzymes in plant evolutionary history.
Evolution of the p-hydroxybenzoic Acid Geranyltransferase Gene for Shikonin Biosynthesis
A boraginaceaeous medicinal plant, L. erythrorhizon, possesses a unique subfamily of p-hydroxybenzoic acid geranyltransferases (PGTs) (Figure 1B) that are specifically involved in shikonin biosynthesis (Yazaki et al., 2002). An overview of the evolutionary history of PGT was attained by assessing genome sequences and transcriptomes of L. erythrorhizon from the GenBank datasets SRP108575 and SRP141330, respectively, as well as by reassembling our original data (Takanashi et al., 2019). The hypothetical PGT-like proteins were found to beclosely related to the ubiquinone prenyltransferase subfamily involved in common metabolism (magenta in Figure 1B), which was closer to these hypothetical PGT-like proteins than the specialized citrus prenyltransferases (Figure 1A). Most PGT-like proteins are encoded by genes with a single exon, whereas general ubiquinone biosynthetic polyprenyltransferases (PPTs) are encoded by genes containing multiple exons (Figure 1C). It is of interest to determine how the single exon structure was generated during the evolution of plant specialized metabolism.
Missing Ubiquinone Prenyltransferase in L. erythrorhizon
Although ubiquinone is a common metabolite in all eukaryotes, and the genes encoding PPTs are essential for the survival of a broad range of organisms, no orthologous ubiquinone PPT gene was found in the L. erythrorhizon transcriptome. Experiments in yeast showed that LePGT cannot synthesize ubiquinone (Yazaki et al., 2002), and LePGT heterologously expressed in E. coli was found to inhibit ubiquinone biosynthesis (Wu et al., 2015). Genomic sequencing identified a contig fragment that could code for PPT-like proteins (asterisk in Figure 1B) and that its amino acid sequence was moderately similar to that of OsPPT1, which is responsible for ubiquinone biosynthesis in rice. In addition, there were three contigs that wecould not classify, which are labeled “unclassified genes” (“unclassified” in Figure 1B). In contrast to the particular PGT that catalyzes shikonin biosynthesis, an intron insertion was found in the hypothetical gene, at the same position as in the PGTs of Nicotiana tabacum and Oryza sativa (Figure 1C). This conserved exon-intron organization was also observed in the PPT genes from Arabidopsis and rice (Ohara et al., 2006). This gene product is a strong candidate for a ubiquinone prenyltransferase in L. erythrorhizon, and its biochemical characterization is expected in the future.
Evolution of the Taxus Acyltransferase Gene Family
Acyltransferases also substantially contribute to the diversification of specialized metabolites, in which BAHD and SCPL (serine carboxypeptidase-like) are representatives. Taxoids such as paclitaxel present in Taxus species are specialized metabolites and highly acylated compounds. Five known taxoid acyltransferases are closely related to each other, with all grouped in clade V of the BAHD acyltransferase family (D’Auria, 2006). These Taxus proteins differ in substrate specificities for both acyl donors and acceptors; i.e., they can utilize acetyl-CoA, benzoyl-CoA or phenylalanoyl-CoA for O- and N-acylation of various taxoid molecules (D’Auria, 2006).
To understand the evolutionary development of the Taxus BAHD acyltransferase family, BAHD clade V was analyzed phylogenetically in detail (yellow background in Figure 2A). The amino acid sequences of Taxus BAHD members were obtained from the transcriptome data of Taxus x media cultured cells (Yukimune et al., 1996). Phylogenetic analysis showed that the Taxus BAHD proteins form a Taxus-specific clade (red bracket in Figure 2A), containing all five characterized acyltransferases (white background in the Taxus-specific clade), as well as other Taxus proteins of unknown function (asterisk in Figure 2A). Within this clade of the BAHD family, O. sativa and A. thaliana each form a unique clade, suggesting that lineage-specific subfamily expansion of the BAHD acyltransferases plays a major role in plant evolution (Fani and Fondi, 2009). In addition to this Taxus-specific subgroup, other Taxus BAHD proteins have been identified, with these classified with other model plant BAHD members (Supplementary Figure S1), suggesting that Taxus species possess genes encoding general BAHD clade V proteins that are conserved among a broad range of plant species.
Figure 2
It can behypothesized that neo-functionalization is induced by the acquisition of promiscuous enzymatic activity during plant evolution. Wehave examined the enzymatic activity of recombinant proteins prepared from seven isolated cDNAs encoding BAHD members of the Taxus-specific subfamily (dagger in Figure 2A). Each crude recombinant enzyme was prepared using pET22a and OrigamiB as a host-vector system (Novagen), without a periplasmic signal sequence, according to the conventional method. Each enzyme was reacted with acetyl-CoA and 10-deacetyl baccatin III (10-DAB) as substrates, and the reaction products were analyzed using an UPLC–MS/MS system equipped with a BEH C18 column (Waters). The clone encoding 5-hydroxytaxadiene 5-O-acetyltransferase (TAT) had 10-DAB:10-O-acetyltransferase (DBAT) activity (Walker et al., 2000), as well as the canonical enzyme DBAT (Figures 2B,C; Walker and Croteau, 2000). The amount of the product formed by the substrate was 1.4mol% for TAT and 10.4% for DBAT, suggesting that the activity of TAT was 13.2% that of DBAT. This promiscuity of enzymatic activity may represent the evolutionary footprint of a biosynthetic enzyme that acquires a new functionality through the alteration of substrate and product specificities, resulting in the production of a unique specialized metabolite.
Conclusions and Perspectives
Using two transferase subfamilies as examples, wehave shown the “heritage” of expansion of a gene family, which is relevant for the development of plant specialized metabolic pathways. A protein in the specific BAHD subfamily of Taxus species showed promiscuous enzymatic activity for noncanonical substrates containing side chains at a noncanonical carbon position. These observations fit the general context of developmental molecular evolution that explains the development and establishment of new canonical enzymatic activity (Weng et al., 2012). The generation in L. erythrorhizon of a PGT gene subfamily, each containing a single exon and involved in shikonin biosynthesis, suggests the putative involvement of the reverse transcription of mature mRNA. If this surmise is valid for other enzyme families, single exon genes may provide clues to identifying missing proteins responsible for biosynthetic pathways for other valuable plant specialized metabolites.
There are yet many missing links, even in actively studied shikonin and taxoid biosynthetic pathways. The applicable range of the single exon hypothesis may not belimited only to biosynthetic enzymes, but to regulatory factors. The identification of regulatory factors will beessential to understanding the production of plant specialized metabolites, including membrane transporters. Comparative genomics will enable the assessment of the evolutionary footprint of these genes, e.g., the expansion of specific subfamilies and the proliferation of single exon genes. Further biochemical and molecular genetics studies may provide experimental evidence for the involvement of hypothetical proteins in plant specialized metabolism.
Supplementary Material
The Supplementary Material for this article can befound online at: https://www.frontiersin.org/articles/10.3389/fpls.2019.00794/full#supplementary-material
FIGURE S1Expanded phylogenetic tree of Figure 2A. Phylogenetic analysis with hypothetical Taxus BAHD acyltransferase-like proteins and related proteins from model plant species. Asterisks indicate Taxus proteins found in this study. Functionally identified BAHD proteins are highlighted in yellow background.
Statements
Data availability statement
The datasets generated for this study can befound in GenBank.
Author contributions
HK and KY wrote the manuscript and performed the phylogenetic and biochemical analyses. HL was involved in the assembly of genomic contigs and the analysis of the exon-intron structure of Lithospermum erythrorhizon genes. HM, YK, and HT were responsible for transcriptome analysis of Taxus spp.
Funding
This work was supported in part by the New Energy and Industrial Technology Development Organization (NEDO, No. 16100890 to KY). Additional support was provided by the Mission Research of RISH, Kyoto University.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could beconstrued as a potential conflict of interest.
References
1
AkashiT.SasakiK.AokiT.AyabeS.YazakiK. (2009). Molecular cloning and characterization of a cDNA for pterocarpan 4-dimethylallyltransferase catalyzing the key prenylation step in the biosynthesis of glyceollin, a soybean phytoalexin. Plant Physiol.149, 683–693. doi: 10.1104/pp.108.123679
2
BalogluE.KingstonD. G. I. (1999). The taxane diterpenoids. J. Nat. Prod.62, 1448–1472. doi: 10.1021/np990176i
3
BanksJ. A.NishiyamaT.HasebeM.BowmanJ. L.GribskovM.dePamphilisC.et al. (2011). The Selaginella genome identifies genetic changes associated with the evolution of vascular plants. Science332, 960–963. doi: 10.1126/science.1203810
4
BourgaudF.HehnA.LarbatR.DoerperS.GontierE.KellnerS.et al. (2006). Biosynthesis of coumarins in plants: a major pathway still to beunravelled for cytochrome P450 enzymes. Phytochem. Rev.5, 293–308. doi: 10.1007/s11101-006-9040-2
5
BowmanJ. L.KohchiT.YamatoK. T.JenkinsJ.ShuS.IshizakiK.et al. (2017). Insights into land plant evolution garnered from the Marchantia polymorpha genome. Cell171, 287–304. doi: 10.1016/j.cell.2017.09.030
6
CroteauR.KetchumR. E. B.LongR. M.KasperaR.WildungM. R. (2006). Taxol biosynthesis and molecular genetics. Phytochem. Rev.5, 75–97. doi: 10.1007/s11101-005-3748-2
7
D’AuriaJ. C. (2006). Acyltransferases in plants: a good time to beBAHD. Curr. Opin. Plant Biol.9, 331–340. doi: 10.1016/j.pbi.2006.03.016
8
FaniR.FondiM. (2009). Origin and evolution of metabolic pathways. Phys. Life. Rev.6, 23–52. doi: 10.1016/j.plrev.2008.12.003
9
Guerra-BubbJ.CroteauR.WilliamsR. M. (2012). The early stages of taxol biosynthesis: an interim report on the synthesis and identification of early pathway metabolites. Nat. Prod. Rep.29, 683–696. doi: 10.1039/c2np20021j
10
HayashiK.HorieK.HiwatashiY.KawaideH.YamaguchiS.HanadaA.et al. (2010). Endogenous diterpenes derived from ent-kaurene, a common gibberellin precursor, regulate protonema differentiation of the moss Physcomitrella patens. Plant Physiol.153, 1085–1097. doi: 10.1104/pp.110.157909
11
KaramatF.OlryA.MunakataR.KoedukaT.SugiyamaA.ParisC.et al. (2014). A coumarin-specific prenyltransferase catalyzes the crucial biosynthetic reaction for furanocoumarin formation in parsley. Plant J.77, 627–638. doi: 10.1111/tpj.12409
12
KroymannJ. (2011). Natural diversity and adaptation in plant secondary metabolism. Curr. Opin. Plant Biol.14, 246–251. doi: 10.1016/j.pbi.2011.03.021
13
LiW. (2016). Bringing bioactive compounds into membranes: the UbiA superfamily of intramembrane aromatic prenyltransferases. Trends Biochem. Sci.41, 356–370. doi: 10.1016/j.tibs.2016.01.007
14
MunakataR.InoueT.KoedukaT.KaramatF.OlryA.SugiyamaA.et al. (2014). Molecular cloning and characterization of a geranyl diphosphate-specific aromatic prenyltransferase from lemon. Plant Physiol.166, 80–90. doi: 10.1104/pp.114.246892
15
MunakataR.OlryA.KaramatF.CourdavaultV.SugiyamaA.DateY.et al. (2016). Molecular evolution of parsnip (Pastinaca sativa) membrane-bound prenyltransferases for linear and/or angular furanocoumarin biosynthesis. New Phytol.211, 332–344. doi: 10.1111/nph.13899
16
OharaK.YamamotoK.HamamotoM.SasakiK.YazakiK. (2006). Functional characterization of OsPPT1, which encodes p-hydroxybenzoate polyprenyltransferase involved in ubiquinone biosynthesis in Oryza sativa. Plant Cell Physiol.47, 581–590. doi: 10.1093/pcp/pcj025
17
SadreR.FrentzenM.SaeedM.HawkesT. (2010). Catalytic reactions of the homogentisate prenyl transferase involved in plastoquinone-9 biosynthesis. J. Biol. Chem.285, 18191–18198. doi: 10.1074/jbc.M110.117929
18
SasakiK.MitoK.OharaK.YamamotoH.YazakiK. (2008). Cloning and characterization of naringenin 8-prenyltransferase, a flavonoid-specific prenyltransferase of Sophora flavescens. Plant Physiol.146, 1075–1084. doi: 10.1104/pp.107.110544
19
SavidgeB.WeissJ. D.WongY. H.LassnerM. W.MitskyT. A.ShewmakerC. K.et al. (2002). Isolation and characterization of homogentisate phytyltransferase genes from Synechocystis sp. PCC 6803 and Arabidopsis. Plant Physiol.129, 321–332. doi: 10.1104/pp.010747
20
TakahashiS.OgiyamaY.KusanoH.ShimadaH.KawamukaiM.KadowakiK. (2006). Metabolic engineering of coenzyme Q by modification of isoprenoid side chain in plant. FEBS Lett.580, 955–959. doi: 10.1016/j.febslet.2006.01.023
21
TakanashiK.NakagawaY.AburayaS.KaminadeK.AokiW.Saida-MunakataY.et al. (2019). Comparative proteomic analysis of Lithospermum erythrorhizon reveals regulation of a variety of metabolic enzymes leading to comprehensive understanding of the shikonin biosynthetic pathway. Plant Cell Physiol.60, 19–28. doi: 10.1093/pcp/pcy183
22
VenkateshT. V.KarunanandaaB.FreeD. L.RottnekJ. M.BaszisS. R.ValentinH. E. (2006). Identification and characterization of an Arabidopsis homogentisate phytyltransferase paralog. Planta223, 1134–1144. doi: 10.1007/s00425-005-0180-1
23
WalkerK.CroteauR. (2000). Molecular cloning of a 10-deacetylbaccatin III-10-O-acetyl transferase cDNA from Taxus and functional expression in Escherichia coli. Proc. Natl. Acad. Sci. USA97, 583–587. doi: 10.1073/pnas.97.2.583
24
WalkerK.SchoendorfA.CroteauR. (2000). Molecular cloning of a taxa-4(20),11(12)-dien-5α-ol-O-acetyl transferase cDNA from Taxus and functional expression in Escherichia coli. Arch. Biochem. Biophys.374, 371–380. doi: 10.1006/abbi.1999.1609
25
WengJ. K.PhilippeR. N.NoelJ. P. (2012). The rise of chemodiversity in plants. Science336, 1667–1670. doi: 10.1126/science.1217411
26
WuH.BennettG. N.SanK.-Y. (2015). Metabolic control of respiratory levels in coenzyme Q biosynthesis-deficient Escherichia coli strains leading to fine-tune aerobic lactate fermentation. Biotechnol. Bioeng.112, 1720–1726. doi: 10.1002/bit.25585
27
YazakiK. (2017). Lithospermum erythrorhizon cell cultures: present and future aspects. Plant Biotechnol.34, 131–142. doi: 10.5511/plantbiotechnology.17.0823a
28
YazakiK.ArimuraG.OhnishiT. (2017). “Hidden” terpenoids in plants: their biosynthesis, localization and ecological roles. Plant Cell Physiol.58, 1615–1621. doi: 10.1093/pcp/pcx123
29
YazakiK.KunihisaM.FujisakiT.SatoF. (2002). Geranyl diphosphate:4-hydroxybenzoate geranyltransferase from Lithospermum erythrorhizon. Cloning and characterization of a key enzyme in shikonin biosynthesis. J. Biol. Chem.277, 6240–6246. doi: 10.1074/jbc.M106387200
30
YoneyamaK.AkashiT.AokiT. (2016). Molecular characterization of soybean pterocarpan 2-dimethylallyltransferase in glyceollin biosynthesis: local gene and whole-genome duplications of prenyltransferase genes led to the structural diversity of soybean prenylated isoflavonoids. Plant Cell Physiol.57, 2497–2509. doi: 10.1093/pcp/pcw178
31
YukimuneY.TabataH.HigashiY.HaraY. (1996). Methyl jasmonate-induced overproduction of paclitaxel and baccatin III in Taxus cell suspension cultures. Nat. Biotechnol.14, 1129–1132. doi: 10.1038/nbt0996-1129
Summary
Keywords
prenyltransferase, acyltransferase, BAHD, Citrus, gene family, molecular evolution, specialized metabolism, Lithospermum, Taxus
Citation
Kusano H, Li H, Minami H, Kato Y, Tabata H and Yazaki K (2019) Evolutionary Developments in Plant Specialized Metabolism, Exemplified by Two Transferase Families. Front. Plant Sci. 10:794. doi: 10.3389/fpls.2019.00794
Received
01 March 2019
Accepted
31 May 2019
Published
25 June 2019
Volume
10 - 2019
Edited by
Kazuki Saito, RIKEN Center for Sustainable Resource Science (CSRS), Japan
Reviewed by
Toru Nakayama, Tohoku University, Japan; Suvi Tuulikki Häkkinen, VTT Technical Research Centre of Finland Ltd, Finland
Updates
Copyright
© 2019 Kusano, Li, Minami, Kato, Tabata and Yazaki.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kazufumi Yazaki, yazaki@rish.kyoto-u.ac.jp
This article was submitted to Plant Metabolism and Chemodiversity, a section of the journal Frontiers in Plant Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.