Abstract
As a major staple food, maize (Zea mays) is critical to food security. Shifting environmental pressures increasingly hamper crop defense capacities, causing expanded harvest loss. Specialized labdane-type diterpenoids are key components of maize chemical defense and ecological adaptation. Labdane diterpenoid biosynthesis most commonly requires the pairwise activity of class II and class I diterpene synthases (diTPSs) that convert the central precursor geranylgeranyl diphosphate into distinct diterpenoid scaffolds. Two maize class II diTPSs, ANTHER EAR 1 and 2 (ZmAN1/2), have been previously identified as catalytically redundant ent-copalyl diphosphate (CPP) synthases. ZmAN1 is essential for gibberellin phytohormone biosynthesis, whereas ZmAN2 is stress-inducible and governs the formation of defensive kauralexin and dolabralexin diterpenoids. Here, we report the biochemical characterization of the two remaining class II diTPSs present in the maize genome, COPALYL DIPHOSPHATE SYNTHASE 3 (ZmCPS3) and COPALYL DIPHOSPHATE SYNTHASE 4 (ZmCPS4). Functional analysis via microbial co-expression assays identified ZmCPS3 as a (+)-CPP synthase, with functionally conserved orthologs occurring in wheat (Triticum aestivum) and numerous dicot species. ZmCPS4 formed the unusual prenyl diphosphate, 8,13-CPP (labda-8,13-dien-15-yl diphosphate), as verified by mass spectrometry and nuclear magnetic resonance. As a minor product, ZmCPS4 also produced labda-13-en-8-ol diphosphate (LPP). Root gene expression profiles did not indicate an inducible role of ZmCPS3 in maize stress responses. By contrast, ZmCPS4 showed a pattern of inducible gene expression in roots exposed to oxidative stress, supporting a possible role in abiotic stress responses. Identification of the catalytic activities of ZmCPS3 and ZmCPS4 clarifies the first committed reactions controlling the diversity of defensive diterpenoids in maize, and suggests the existence of additional yet undiscovered diterpenoid pathways.
Introduction
Plant labdane-related diterpenoids represent a diverse group of more than 7,000 metabolites with broad physiological functions in plant development, defense, and ecological adaptation (; Zerbe and Bohlmann, 2015). A few widely conserved diterpenoid metabolites and the corresponding metabolic enzymes are essential for gibberellin biosynthesis of general metabolism (; Zi et al., 2014). From these ancestral diterpenoid pathways, the vast chemical space of specialized diterpenoids has evolved (Zi et al., 2014), which comprises common and species-specific diterpenoid blends with distinct biological functions that range from pest and pathogen defense to allelopathic and signaling activities (; ; ; Tholl, 2015).
Two of the most agroeconomically important grain crops, maize (Zea mays) and rice (Oryza sativa), deploy distinct networks of diterpenoids to mediate the plant response to biotic and abiotic stress (for review, see ; ). Rice produces a complex suite of diterpenoid phytoalexins, including momilactones, oryzalexins, and phytocassanes, that serve as major components of disease resistance and exhibit allelopathic properties to suppress the growth of competing weeds (; ; Toyomasu et al., 2014; ). Similarly, maize produces related yet distinct arsenals of bioactive diterpenoids that, to current knowledge, include the kauralexin and dolabralexin groups with demonstrated and predicted functions in maize chemical defense (Figure 1; ; Vaughan et al., 2015; ; ). For example, pathogen-elicited maize kauralexin and dolabralexin diterpenoids exhibit potent antimicrobial efficacies in vitro and in vivo against several pathogens (; ; ), including species of Fusarium as a major causal agent of maize crop losses and mycotoxin contamination (). In addition, kauralexins exhibit antifeedant activity in response to herbivore attack by the European corn borer (Ostrinia nubilalis; ; ). More recently, both kauralexins and dolabralexins have been shown to accumulate under drought and below-ground oxidative stress, consistent with a possible protective role of diterpenoids in abiotic stress responses (Vaughan et al., 2015; ). Furthermore, a kauralexin- and dolabralexin-deficient maize mutant an2 (anther ear 2) showed increased susceptibility to both drought stress and pathogen attack (Vaughan et al., 2015; ). These findings highlight the biological importance of diterpenoids in conferring stress resilience, and underscore the need to better understand the biosynthesis and chemical ecology of diterpenoids in cereal crops.
FIGURE 1
In angiosperms, the structural diversity of labdane-type diterpenoids is determined by the activity of class II and class I diTPS enzymes that act sequentially to transform the central precursor GGPP into different scaffolds (; Zerbe et al., 2015). Class II diTPSs control the first committed reaction, catalyzing the protonation-initiated cyclization of GGPP into bicyclic labdadienyl/CPP intermediates of distinct stereochemistry and/or regio-specific oxygenations (; Zerbe and Bohlmann, 2015). Class I diTPSs then convert these intermediates via ionization of the diphosphate group and a variety of possible downstream cyclization and rearrangement reactions of the intermediary carbocation (; Zerbe and Bohlmann, 2015). Multi-gene diTPS families of 9–31 members were identified in the genomes of maize, rice, wheat (Triticum aestivum), and switchgrass (Panicum virgatum; ; Xu et al., 2007a; Wu et al., 2012; Zhou et al., 2012; ; ). Functional characterization of class II diTPSs in these species demonstrated the formation of both common and distinct CPP stereoisomers. All species produce ent-CPP as a precursor for gibberellin phytohormones, as well as specialized diterpenoids at least in rice and maize (). In addition, wheat forms the enantiomeric (+)-CPP en route to pimarane and abietane diterpenoids, whereas rice and switchgrass form syn-CPP en route to, for example, rice oryzalexins and momilactones (Xu et al., 2004, 2007a; Wu et al., 2012; Zhou et al., 2012; ). The switchgrass class II diTPS family appears to have functionally diverged more extensively to also produce 8,13-CPP (labda-8,13-dien-15-yl diphosphate) and the clerodane diterpenoid precursor clerodienyl diphosphate (). In maize, only two of the four class II diTPSs present in the genome (B73 RefGen_v4), namely, ANTHER EAR 1 (ZmAN1, Zm00001d032961) and ANTHER EAR 2 (ZmAN2, Zm00001d029648), have been functionally analyzed and demonstrated to produce ent-CPP (; ). Despite their catalytic redundancy, ZmAN1 and ZmAN2 serve different physiological functions in maize. Knock-out mutants of Zman1 display characteristic gibberellin-deficient phenotypes including dwarfism and anther formation in ears that demonstrate a role of ZmAN1 in general metabolism (). By contrast, ZmAN2 is stress-inducible and an2 mutants exhibit normal growth and reproductive phenotypes, but display enhanced susceptibility to fungal disease and environmental stress associated with the lack of kauralexins and/or dolabralexins (; Vaughan et al., 2015; ; ). While these genetic studies clarified the dependency of kauralexin and dolabralexin metabolism on ZmAN2 activity, knowledge of the downstream class I diTPS reactions remains incomplete. Only KAURENE SYNTHASE-LIKE 4 (ZmKSL4, Zm00001d032858) and the P450s CYP71Z16/18 (Zm00001d014136/Zm00001d014134) have been recently identified to catalyze the conversion of ent-CPP into dolabradiene and downstream dolabralexins (). These enzymes expand our knowledge of previously reported maize ent-kaurene synthases (ZmKSL3, ZmKSL5, and ZmTPS1) that convert ent-CPP into ent-kaurene en route to gibberellin biosynthesis and possibly specialized diterpenoid metabolism ().
Herein reported is the biochemical characterization of the two remaining class II diTPSs present in the maize genome, COPALYL DIPHOSPHATE SYNTHASE 3 (ZmCPS3, Zm00001d024512) and COPALYL DIPHOSPHATE SYNTHASE 4 (ZmCPS4, Zm00001d048874), to delineate the scope of yet unknown specialized diterpenoid pathways in maize. Unlike ZmAN2, interrogation of transcriptomic and proteomic datasets did not support a role of ZmCPS3 and ZmCPS4 in highly inducible pathogen defenses. A pattern of largely constitutive gene expression of ZmCPS3 and moderately inducible expression of ZmCPS4 under root exposure to abiotic stress may suggest more constitutive functions in maize ecological adaptation.
Materials and Methods
Gene Constructs
Full-length genes of maize ZmCPS3 (Zm00001d024512) and ZmCPS4 (Zm00001d048874) were synthesized with support of a Department of Energy Joint Genome Institute Community Science Program grant (CSP#2568). Genes were inserted into the second multiple cloning site of a pACYC-Duet plasmid also carrying the GGPP synthase from Abies grandis () to form the constructs pACYC-Duet:AgGGPPS-ZmCPS3 and pACYC-Duet:AgGGPPS-ZmCPS4 (Supplementary Table S1). Additional constructs used in co-expression assays were described previously, including pACYC-Duet:AgGGPPS-ZmAN2 (Zm00001d029648) (Z. mays ent-CPP synthase), pET28b:ZmKSL3 (Zm00001d002349) (Z. mays ent-kaurene synthase), pET28b:ZmKSL4 (Zm00001d032858) (Z. mays dolabradiene synthase), pET28b:GrTPS1 (Grindelia robusta LPP synthase), and pET15b:MvELS (Marrubium vulgare 9,13-epoxy labd-14-ene synthase; Zerbe et al., 2013, 2014; ).
Enzyme Functional Analysis
Functional co-expression of enzymes was carried out using an engineered Escherichia coli platform for enhanced diterpenoid production () as previously described (). The pACYC-Duet:AgGGPPS-ZmCPS3, pACYC-Duet:AgGGPPS-ZmCPS4, or pACYC-Duet:AgGGPPS-ZmAN2 constructs were expressed alone or in combination with ZmKSL3, or ZmKSL4. GrTPS1 was expressed in combination with pACYC-Duet:AgGGPPS to form LPP and an additional combination with pET15b:MvELS to form manoyl oxide for use as authentic standards. In brief, cultures were grown in 50 mL Terrific Broth (TB) medium to an OD600 of ∼0.6 at 37°C. After cooling to 16°C, cultures were induced with 1 mM isopropyl-thio-galactoside (IPTG) and 25 mM sodium pyruvate and incubated for a further 72 h. Enzyme products were extracted with 50 mL of hexane, concentrated under air, and resuspended in 1 mL hexane for GC-MS analysis.
GC-MS Analysis
Gas chromatography mass spectrometry analysis of enzyme products was performed on an Agilent 7890B GC interfaced with a 5977 Extractor XL MS Detector at 70 eV and 1.2 mL min-1 He flow, using an Agilent HP5-MS column (30 m, 250 μm i.d., 0.25 μm film) with a sample volume of 1 μL and the following GC parameters: pulsed splitless injection at 250 and 50°C oven temperature; hold at 50°C for 3 min, 20°C min-1 to 300°C, hold 3 min. MS data from 90 to 600 mass-to-charge ratio (m/z) were collected after a 8 min solvent delay. Products were identified using comparison to authentic standards or, where these were not available, comparison to published mass spectra and the National Institute of Standards and Technology (NIST version 2.0) mass spectral library (Agilent).
NMR Analysis
Diterpenoids were produced via large-scale (12 L) enzyme co-expression cultures as described above. Hexane extracts were dried using rotary evaporation, resuspended in hexane, and purified by silica column chromatography (230–400 mesh, grade 60) using a hexane:ethyl acetate gradient as the mobile phase. Fractions were further purified on an Agilent 1100 series HPLC with diode array UV detector and an Agilent ZORBAX Eclipse Plus-C8 column (4.6 mm × 150 mm, 5 microns) at a 0.5 mL min-1 flow rate and H2O/acetonitrile gradient as mobile phase. Product purity was verified using GC-MS analysis as outlined above. Purified products were dissolved in 0.6 mL deuterated chloroform (CDCl3; Sigma-Aldrich) containing tetramethylsilane (TMS). NMR spectra were acquired at room temperature on a Bruker Avance III 800 spectrometer equipped with a 5 mm CPTCI. Chemical shifts were calculated by reference to known CDCl3 (13C 77.23 ppm, 1H 7.24 ppm) signals offset from TMS. All spectra were acquired using standard experiments on a Bruker TopSpin 3.2 software, including 1D 1H and 1D 13C spectra (201 MHz).
Transcriptomics and Proteomics Analysis
For analysis of class II diTPS gene and protein abundance, publicly available transcriptome and proteome inventories were investigated that represent a range of organs and tissues at different developmental stages of healthy maize plants (Walley et al., 2016). All samples derive from B73, with the exception of 2 cm tassels, 1–2 mm anthers, and mature pollen (W23 inbred), and 5-days-old primary root (Mo17 inbred) (Walley et al., 2016). Transcript abundance (as fragments per kilobase of transcript per million mapped reads, FPKM) and protein expression levels were retrieved directly from this public resource (Walley et al., 2016), and were scaled by color either individually by gene or absolute across all four genes of interest.
Abiotic Elicitation of Maize Roots With CuSO4
Plant samples used for gene expression analysis were prepared previously (). Briefly, maize (var. Golden Queen) seed was germinated in the dark on wetted paper for 4 days at 23°C. Seedlings were transferred and grown hydroponically () for 12 days under 16/8 h light/darkness at 28°C, light intensity of 180 μmol photons m-2 s-1, and ∼60% relative humidity. A total of 1 mM CuSO4 or the corresponding water control was added to the hydroponic medium. Root samples were collected at the time points indicated, with three biological replicates per time point, and immediately frozen in liquid nitrogen for further analysis.
Fungal Elicitation of Mature Maize Roots
Plants samples used for gene expression analysis were derived from a previous study (). Maize (var. Mo17) plants were grown in the greenhouse for 53 days in individual, 10 L pots and supplemented with 14-14-14 Scotts Miracle Grow fertilizer. Large nodal roots (≥2 mm dia) were punctured with a 0.6 mm dia steel pin at 1 cm intervals and inoculated with 10 μL of 1 × 107 conidia mL-1 of Fusarium verticillioides (F.v.), Fusarium graminearum (F.g.), or water control at each wound site. To avoid damage to other tissues not undergoing treatment, sampling was limited to roots on the outer edge of the soil, in contact with the vertical plastic pot wall. Root samples were collected after 7 days and immediately frozen in liquid nitrogen before further processing.
Quantitative Real-Time PCR (qPCR)
Gene expression analysis of ZmCPS3 and ZmCPS4 was performed on the CuSO4-treated or pathogen-treated root samples described above. Total RNA was isolated as described elsewhere () and cDNA was synthesized using SuperScript III First-Strand Synthesis Kit (Invitrogen) according to manufacturer’s instructions. Transcript abundance was measured using a BioRad C1000 Touch Thermo Cycler interfaced with a CFX96 Real-Time System, and iTaq Universal SYBR Green Supermix (BioRad) according to manufacturer’s protocols. Mean cycle threshold (Ct) values of at least two technical and three biological replicates were normalized using elongation factor EF1α as a validated reference gene (), and fold change values were calculated using the 2-ΔΔCt method. Gene-specific oligonucleotides used for qPCR analysis were ZmCPS3 [forward: 5′-TGACGTGTGGATTGGGAAGG-3′, reverse: 5′-CGCTCTGTCTGGCTCAAAGA-3′], ZmCPS4 [forward: 5′-CTCAGGCCAGCTTAACGAC-3′, reverse: 5′-CCTTGCCGATCCATACGTC-3′], and EF1α [forward: 5′-TGGGCCTACTGGTCTTACTACTGA-3′, reverse: 5′-ACA TACCCACGCTTCAGATCCT-3′].
Sequence and Phylogenetic Analysis
Protein sequence alignments were generated using the CLCBio software package (Qiagen), followed by manual curation. Maximum likelihood phylogenetic analysis was performed using PhyML-aBayes version 3.0.1 beta () with four rate substitution categories, LG substitution model, BIONJ starting tree, and 500 bootstrap repetitions.
Results
Active Site Determinants Suggest Distinct Functions of ZmCPS3 and ZmCPS4
The maize genome (B73 RefGen_v4) contains four class II diTPSs (), which comprise the known ent-CPP synthases ZmAN1 and ZmAN2 located on chromosome 1 (; ), and the previously uncharacterized ZmCPS3 and ZmCPS4 positioned on chromosomes 4 and 10, respectively. ZmCPS3 and ZmCPS4 share a protein sequence identity of 55% with each other and 46–56% with ZmAN1 and ZmAN2. Phylogenetic analysis placed ZmCPS3 and ZmCPS4 separate from most ent-CPP synthases on a branch primarily consisting of class II diTPSs that produce prenyl diphosphates of specialized metabolism (Supplementary Figure S1). Although this suggested a role of both enzymes in specialized metabolism, the functional diversity of diTPSs in this group did not allow inference of possible ZmCPS3 and ZmCPS4 functions. Therefore, to inform biochemical analyses, we interrogated key active site residues with known impact on class II diTPS product specificity (Xu et al., 2007b; ; , ; ; ; ; Schulte et al., 2018). Previous studies identified a His-Asn catalytic dyad that is widely conserved among ent-CPP synthases, including ZmAN1 and ZmAN2, and was shown to direct class II diTPS catalysis toward ent-CPP formation (, ; ; ). Neither ZmCPS3 nor ZmCPS4 possess a His-Asn catalytic dyad, but instead feature a Leu-Phe residue pair (Figure 2). These residues are consistent with a recently characterized 8,13-CPP synthase from switchgrass (P. virgatum; ). In addition, presence of a Tyr residue in position 497 and 458 of ZmCPS3 and ZmCPS4, respectively, is consistent with known monocot ent-CPP and (+)-CPP synthases, but contrasts known syn-CPP synthases from rice and switchgrass that feature a His residue in this position (; ). Although insufficient to allow an unambiguous functional annotation, these active site characteristics disfavored an ent-CPP or syn-CPP synthase activity, and supported a possible function of ZmCPS3 and ZmCPS4 as 8,13-CPP or (+)-CPP synthases or related specialized class II diTPSs.
FIGURE 2
ZmCPS3 Functions as a (+)-CPP Synthase
To test the predicted enzyme activity of ZmCPS3, a synthetic full-length gene was co-expressed with a GGPP synthase from A. grandis using an in vivo E. coli expression platform engineered for diterpenoid production (). In vivo expression of class II diTPS enzymes using this system readily yields dephosphorylated products, presumably due to the activity of E. coli endogenous phosphatases and thus enables direct hexane extraction and analysis of the corresponding diterpene alcohols (Zerbe et al., 2013, 2015; ; ). For clarity, structures depicted below represent the native prenyl diphosphate products, but were detected via GC-MS and NMR analysis as the corresponding alcohols. Expression of ZmCPS3 resulted in a major product (compound 1; Figure 3) with a retention time of 11.25 min and a fragmentation pattern showing dominant mass ions of m/z 137, 257, and 275 that closely matched the mass spectrum of copalol (i.e., dephosphorylated CPP) produced by ZmAN2 (compound 8; Figure 3). Two additional byproducts detected in the ZmCPS3 product profile represented unconverted GGPP substrate (compound 2) and an unidentified non-diterpenoid contaminant (compound 3; Supplementary Figure S2).
FIGURE 3
Next we defined the stereochemistry of the ZmCPS3 product by co-expressing ZmCPS3 with characterized maize class I diTPSs that display catalytic specificity toward CPP substrates of different stereochemistries. ZmKSL3 converts ent-CPP into ent-kaurene (), whereas ZmKSL4 forms dolabradiene from ent-CPP and pimarane-type diterpene olefins with (+)-CPP as a substrate (). As controls, we co-expressed ZmAN2 with either ZmKSL3 or ZmKSL4 to generate ent-kaurene (compound 10) and dolabradiene (compound 11), respectively (Figure 3 and Supplementary Figure S2). Co-expression of ZmCPS3 and ZmKSL3 did not result in any detectable class I diTPS product (Figure 3). The combined activity of ZmCPS3 with ZmKSL4 did not yield dolabradiene, but resulted in several pimarane-related products (compounds 4–7), the most abundant of which was identified as pimara-8,14-diene (compound 5) by comparison to reference mass spectra (Supplementary Figures S2, S3). These products are consistent with the previously reported activity of ZmKSL4 with the established (+)-CPP synthase, A. grandis abietadiene synthase variant D621A (). On the basis of these results, ZmCPS3 was designated as a (+)-CPP synthase.
ZmCPS4 Produces 8,13-CPP and LPP
Escherichia coli co-expression of the full-length, synthetic gene encoding ZmCPS4 with the A. grandis GGPP synthase yielded a major product (compound 12) with a retention time of 11.26 min, indicating a related but distinct compound as compared to (+)-CPP and ent-CPP formed by ZmCPS3 and ZmAN2 in vitro, respectively (Figure 4). Presence of signature mass ions of m/z 275 and 257 in the fragmentation pattern of this product indicated the expected labdane structure, but additional major mass ions of m/z 205 and 149 suggested a structure distinct from the common ent-CPP, (+)-CPP, or syn-CPP products observed in monocot crops (Xu et al., 2004; ; Wu et al., 2012). Indeed, the fragmentation pattern of compound 12 matched the product of a recently identified class II diTPS from switchgrass that forms 8,13-CPP (Supplementary Figure S2; ). To further verify the identity of this ZmCPS4 product, nuclear magnetic resonance (NMR) spectroscopy was performed. For this purpose, large-scale E. coli co-expression cultures were used to produce an excess of 1 mg of the product, which was then purified using silica column chromatography and semi-preparative high-pressure liquid chromatography (HPLC). For the purified product, 1D 1H and 13C NMR spectra were acquired and validated the ZmCPS4 product as 8,13-CPP in comparison to published spectra (Supplementary Figure S4; ).
FIGURE 4
In addition to the primary 8,13-CPP product, a minor, more polar, ZmCPS4 product with a retention time of 11.79 min was also observed (compound 13; Figure 4). This product was identified as labda-13-en-8-ol diphosphate (LPP) based on comparison to the retention time and mass spectrum of LPP produced by a known LPP synthase from G. robusta (GrTPS1; Figure 4; Zerbe et al., 2013). Low abundance of this product prevented further stereochemical analysis via NMR. Additional byproducts (compounds 14–16) present in the ZmCPS4 product profile represent manoyl oxide (compound 15) as based on comparison to a manoyl oxide standard produced from the coupled reaction of G. robusta GrTPS1 and MvELS from M. vulgare (Zerbe et al., 2013; ), and a closely related unidentified diterpenoid (compound 14; Figure 4 and Supplementary Figure S2). These products are likely resulting from rearrangement of the ZmCPS4 products 8,13-CPP and LPP after dephosphorylation by endogenous E. coli phosphatases, as has been described for various related class II diTPSs (Zerbe et al., 2015; ; ). To test for possible downstream products of 8,13-CPP and LPP, ZmCPS4 was co-expressed with the currently known maize class I diTPS functions, including the ent-kaurene synthase ZmKSL3 and the dolabradiene synthase ZmKSL4 (; ). No new products were detected when co-expressing ZmCPS4 and ZmKSL3 as compared to the expression of ZmCPS4 alone (Figure 4). When combining ZmCPS4 with ZmKSL4, compound 15 identified as manoyl oxide was significantly increased and, albeit at low abundance, an additional product (compound 17) was formed that predictably represents a labdane diterpene olefin as based on characteristic mass fragments of m/z 272 and 257 (Supplementary Figure S2).
ZmCPS3 and ZmCPS4 Show Distinct Patterns of Tissue-Specific Gene Expression
To further investigate the functions of ZmCPS3 and ZmCPS4, we next examined their gene expression patterns using publicly available transcriptome and proteome data across various organs, tissues, and developmental stages of unstressed maize plants (Walley et al., 2016). With the exception of 5-days-old primary root (Mo17 inbred) and 2 cm tassels, 1–2 mm anthers, and mature pollen (W23 inbred), all samples were from the B73 genotype (Walley et al., 2016). Detectable transcripts of ZmCPS3 were distributed across all organs and tissues with highest abundance included germinating kernels (Figure 5A). Notably, when specifically analyzing gene expression in roots, transcript of ZmCPS3 was observed to be significantly more abundant with expression levels 5–481 times higher as compared to ZmAN1, ZmAN2, and ZmCPS4 in the same tissue (Figure 5B). Transcript of ZmCPS4 was detected at only trace levels and present exclusively in primary and secondary root and root cortex tissues (Figure 5). Consistent with observed transcript abundance, query of public proteome data (Walley et al., 2016) showed that protein levels of ZmCPS3 were highest in root tissues (Supplementary Figure S5). ZmCPS4 protein was not detected in the analyzed proteome.
FIGURE 5
ZmCPS3 and ZmCPS4 Show Distinct Gene Expression Patterns in Response to Biotic and Abiotic Stress
Previous studies demonstrated that ZmAN2 gene expression is induced under both pathogen (F.g. and F.v.) and oxidative stress in above- (; ) and below-ground tissues (). In the context of these findings, the relatively higher transcript abundance of ZmCPS3 in roots and the low but possibly root-specific expression of ZmCPS4 suggested a role of these enzymes in belowground stress responses. To investigate this hypothesis, we analyzed transcript abundance of ZmCPS3 and ZmCPS4 using a previously reported set of samples of 53-days-old maize Mo17 roots incision-inoculated with fungal spores of F.v. and F.g. and harvested 7 days later (). Plants treated by incision-wounding only were used as controls. QPCR analysis showed no significant fold change in ZmCPS3 transcript abundance in roots exposed to F.v. or F.g. as compared to wounded controls (Figure 6). On average, gene expression of ZmCPS4 was significantly decreased in response to F.v. and F.g. elicitation as compared to control plants.
FIGURE 6
To examine ZmCPS3 and ZmCPS4 gene expression in response to abiotic stress, qPCR analyses were performed on two-week-old roots of maize (var. Golden Queen) plants that were exposed to 1 mM CuSO4 treatment (as a proxy for oxidative stress), previously shown to induce diterpenoid biosynthesis (). All treatments were compared to a water-treated control at each time point using the 2-ΔΔCt method in which the control has a fold change of 1. Transcript levels of ZmCPS3 in CuSO4-treated roots did not differ significantly from those observed in the roots of water-treated control plants at 2 and 4 h post treatment, but were significantly reduced at the 24 h post treatment time point (Figure 6) Conversely, ZmCPS4 showed increased transcript abundance in CuSO4-treated roots as early as 2 h post treatment and with an up to sixfold change at 4 h, before decreasing again after 12 h, with a peak of sixfold increase in transcript abundance (Figure 6).
Discussion
Rapidly increasing biotic and abiotic pressures can overcome the natural defense systems of plants, leading to substantial harvest losses in major food crops (; ). Given the proven protective properties of crop-specific diterpenoid arsenals, a deeper knowledge of their biosynthesis and biological functions may aid new solutions to optimize crop stress resilience traits and mitigate yield loss (). Elucidation of the enzyme activities of ZmCPS3 and ZmCPS4 completes the functional range of the maize class II diTPS family, which controls the early committed reactions responsible for the diversity of maize diterpenoid pathways.
Advances in the discovery and mechanistic analysis of diTPSs across the plant kingdom increasingly enable the prediction of diTPS functions (Zerbe and Bohlmann, 2015), as exemplified here for ZmCPS3 and ZmCPS4. However, accurate computational annotation of diTPS activities remains to be constrained by the vast sequence and functional space of the enzyme family, thus necessitating biochemical characterization. Modular co-expression assays with an expanding catalog of diTPSs of known substrate/product-specificity can be leveraged for efficient diTPS functional analysis, and were applied in this study for the identification of ZmCPS3 as a (+)-CPP synthase and ZmCPS4 as an 8,13-CPP synthase (Figures 3, 4). While absent in rice (Xu et al., 2007a), a (+)-CPP synthase has been identified in wheat, where (+)-CPP can be further converted by class I diTPSs to form labdane structures, such as pimara-8,15-diene, abietadiene, and isopimara-7,15-diene (Wu et al., 2012; Zhou et al., 2012). While corresponding end-products and physiological functions have yet to be discovered in maize and wheat, the wide distribution of (+)-CPP-derived diterpenoids in both angiosperm and gymnosperm species suggests roles in stress defense (; Zhou et al., 2012; Zerbe and Bohlmann, 2015). Similar to the recently demonstrated accumulation and predicted defensive functions of dolabralexins in maize roots (), a possible function of ZmCPS3 in belowground stress responses could be hypothesized based on the higher transcript abundance in roots as compared to ZmCPS4, ZmAN1, and ZmAN2. However, the observed lack of elicited ZmCPS3 expression in maize roots in response to Fusarium infection or oxidative stress did not support an inducible defense role for ZmCPS3-derived pathways and metabolites. Although ZmCPS3 and its protein product were more highly expressed in roots as compared to other healthy tissues, its broad tissue- and development-wide presence may suggest a possibly more constitutive function. Tissue-wide expression levels of ZmCPS3 suggest a conceptual parallel to the high constitutive levels of benzoxazinoid pathway enzymes and metabolites in maize seedlings, which have been shown to be important for the maize biotic and abiotic stress responses (; Wouters et al., 2016). The largely constitutive and moderately inducible role is supported by a separate yet related study, where ZmCPS3 transcript levels displayed an approximately twofold accumulation in the leaves of a resistant maize recombinant inbred line two weeks after challenge with gray leaf spot (Cercospora zeina; ).
In addition to the (+)-CPP synthase ZmCPS3, characterization of ZmCPS4 as an 8,13-CPP synthase adds an uncommon scaffold to the diterpenoid network of maize (; ; ). Mechanistically, the proposed ZmCPS4-catalyzed reaction will proceed through terminal deprotonation of the common (+)-labda-13E-en-8-yl+ diphosphate intermediate at C-9 to yield 8,13-CPP, as opposed to the more typical deprotonation of the carbocation at the exocyclic C-17 methyl to form the (+)-CPP or ent-CPP isomers (). ZmCPS4-catalyzed formation of LPP as a minor product will require hydroxylation at C-8 of the (+)-labda-13E-en-8-yl+ diphosphate intermediate (; ; Zerbe et al., 2012). Dual product activity is rarely observed in class II diTPSs (). Lacking available ZmCPS4 maize mutants to enable in planta gene function studies, it can only be speculated if LPP represents a native ZmCPS4 product or results from a possibly reduced enzyme activity in vitro that could cause a slower conversion of the intermediary carbocation and thus facilitate water quenching toward formation of LPP. Notably, an 8,13-CPP synthase was also recently discovered in switchgrass, where the enzyme was characterized as a single product class II diTPS (). By contrast, an 8,13-CPP synthase function is absent in rice and wheat (). Although a gene loss in these species cannot be excluded, this selective presence of 8,13-CPP synthases indicates the independent evolution of this function in maize and switchgrass. Similar to ZmCPS3, no pathogen-elicited gene expression was observed for ZmCPS4; in fact, transcript abundance was significantly decreased in response to Fusarium infection (Figure 6). This contrasts the well-established pathogen-elicited gene expression of ZmAN2 (; ; ; ), and may suggest that ZmAN2-derived pathways are the predominant inducible mediators of pathogen defense responses (; ). Conversely, CuSO4-induced expression of ZmCPS4 in roots points to a possible role in belowground plant-environment interactions as also proposed for kauralexins and dolabralexins (Vaughan et al., 2015; ) (Figure 6).
The identity of predominant pathway end-products derived from (+)-CPP, 8,13-CPP, and possibly LPP, as well as their corresponding roles in plant-environmental adaptation will require further investigation and ultimately the generation and analysis of defined genetic mutants in future studies. Currently, the biochemical characterization of ZmCPS3 and ZmCPS4 expands the known chemical landscape surrounding maize diterpenoid metabolism and completes the characterization of predicted class II diTPSs in the maize genome. Modular pathway networks composed of class II and class I diTPSs are likely operating in maize to convert the ZmCPS3 and/or ZmCPS4 products into a broader spectrum of specialized diterpenoids. In vitro formation of several labdane-type diterpenoids through the sequential activity of ZmKSL4, but not the ent-kaurene synthase ZmKSL3, with ZmCPS3 and ZmCPS4 support this hypothesis. However, no (+)-CPP, 8,13-CPP or LPP derivatives have yet been demonstrated in maize and the biological role of ZmCPS3/4 pathway branches has to be demonstrated in planta. Nevertheless, functional knowledge of all maize class II diTPSs will now enable a detailed investigation of modular diterpenoid-metabolic pathway branches in maize that are formed by class II diTPSs and known or thus far uncharacterized class I diTPSs as the biochemical foundation for maize diterpenoid diversity. Such insight can be of substantial value to elucidate and ultimately harness the genetic basis of crop stress resilience (; ).
Statements
Author contributions
PZ conceived the original research and oversaw data analysis. KMM performed most experiments. L-TM performed sub-cloning experiments and assisted with enzyme characterization. YD and ES performed stress application to plants. KMM and PZ wrote the article with contributions from all authors.
Funding
This work was supported by a grant from the NSF Plant-Biotic Interactions Program (1758976 to ES and PZ), by a DOE Joint Genome Institute Community Science Program (JGI-CSP) grant (CSP2568 to PZ and ES), by a UC Davis Academic Senate Faculty Research Grant (to PZ), by a fellowship through the Hellman Foundation (to PZ), and by a fellowship provided by the NSF Graduate Research Fellowship Program (to KMM).
Acknowledgments
We gratefully acknowledge Dr. Reuben Peters (Iowa State University, United States) for providing the pIRS and pGGxC constructs.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.01542/full#supplementary-material
Abbreviations
- CPP
copalyl diphosphate
- diTPS
diterpene synthase
- GC-MS
gas chromatography mass spectrometry
- GGPP
geranylgeranyl diphosphate
- KSL
kaurene-synthase-like
- P450
cytochrome P450.
References
1
AhmadS.VeyratN.Gordon-WeeksR.ZhangY.MartinJ.SmartL.et al (2011). Benzoxazinoid metabolites regulate innate immunity against aphids and fungi in maize.Plant Physiol.157317–327. 10.1104/pp.111.180224
2
BensenR. J.JohalG. S.CraneV. C.TossbergJ. T.SchnableP. S.MeeleyR. B.et al (1995). Cloning and characterization of the maize An1 gene.Plant Cell775–84. 10.1105/tpc.7.1.75
3
CaniardA.ZerbeP.LegrandS.CohadeA.ValotN.MagnardJ. L.et al (2012). Discovery and functional characterization of two diterpene synthases for sclareol biosynthesis in Salvia sclarea (L.) and their relevance for perfume manufacture.BMC Plant Biol.12:119. 10.1186/1471-2229-12-119
4
ChakrabortyS.NewtonA. C. (2011). Climate change, plant diseases and food security: an overview.Plant Pathol.602–14. 10.1111/j.1365-3059.2010.02411.x
5
ChaturvediR.VenablesB.PetrosR. A.NalamV.LiM.WangX.et al (2012). An abietane diterpenoid is a potent activator of systemic acquired resistance.Plant J.71161–172.
6
ChristensenS. A.SimsJ.VaughanM. M.HunterC.BlockA.WillettD.et al (2018). Commercial hybrids and mutant genotypes reveal complex protective roles for inducible terpenoid defenses in maize.J. Exp. Bot.691693–1705. 10.1093/jxb/erx495
7
CuiG.DuanL.JinB.QianJ.XueZ.ShenG.et al (2015). Functional divergence of diterpene syntheses in the medicinal plant Salvia miltiorrhiza.Plant Physiol.1691607–1618. 10.1104/pp.15.00695
8
DafoeN. J.HuffakerA.VaughanM. M.DuehlA. J.TealP. E.SchmelzE. A. (2011). Rapidly induced chemical defenses in maize stems and their effects on short-term growth of Ostrinia nubilalis.J. Chem. Ecol.37984–991. 10.1007/s10886-011-0002-9
9
de SassiC.TylianakisJ. M. (2012). Climate change disproportionately increases herbivore over plant or parasitoid biomass.PLoS One7:e40557. 10.1371/journal.pone.0040557
10
FalaraV.PicherskyE.KanellisA. K. (2010). A copal-8-ol diphosphate synthase from the angiosperm Cistus creticus subsp. creticus is a putative key enzyme for the formation of pharmacologically active, oxygen-containing labdane-type diterpenes.Plant Physiol.154301–310. 10.1104/pp.110.159566
11
FuJ.RenF.LuX.MaoH.XuM.DegenhardtJ.et al (2016). A Tandem array of ent-Kaurene Synthases in maize with roles in gibberellin and more specialized metabolism.Plant Physiol.170742–751. 10.1104/pp.15.01727
12
GuindonS.DufayardJ. F.LefortV.AnisimovaM.HordijkW.GascuelO. (2010). New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0.Syst. Biol.59307–321. 10.1093/sysbio/syq010
13
HarrisL. J.SaparnoA.JohnstonA.PrisicS.XuM.AllardS.et al (2005). The maize An2 gene is induced by Fusarium attack and encodes an ent-copalyl diphosphate synthase.Plant Mol. Biol.59881–894. 10.1007/s11103-005-1674-8
14
JezJ. M.LeeS. G.SherpA. M. (2016). The next green movement: plant biology for the environment and sustainability.Science3531241–1244. 10.1126/science.aag1698
15
JiaM.PotterK. C.PetersR. J. (2016). Extreme promiscuity of a bacterial and a plant diterpene synthase enables combinatorial biosynthesis.Metab. Eng.3724–34. 10.1016/j.ymben.2016.04.001
16
JiaM.ZhouK.TuftsS.SchulteS.PetersR. J. (2017). A pair of residues that interactively affect diterpene synthase product outcome.ACS Chem. Biol.12862–867. 10.1021/acschembio.6b01075
17
Kato-NoguchiH.PetersR. J. (2013). The role of momilactones in rice allelopathy.J. Chem. Ecol.39175–185. 10.1007/s10886-013-0236-9
18
KeelingC. I.BohlmannJ. (2006a). Diterpene resin acids in conifers.Phytochemistry672415–2423. 10.1016/j.phytochem.2006.08.019
19
KeelingC. I.BohlmannJ. (2006b). Genes, enzymes and chemicals of terpenoid diversity in the constitutive and induced defence of conifers against insects and pathogens.New Phytol.170657–675. 10.1111/j.1469-8137.2006.01716.x
20
KolosovaN.MillerB.RalphS.EllisB. E.DouglasC.RitlandK.et al (2004). Isolation of high-quality RNA from gymnosperm and angiosperm trees.Biotechniques36821–824. 10.2144/04365ST06
21
LuX.ZhangJ.BrownB.LiR.Rodriguez-RomeroJ.BerasateguiA.et al (2018). Inferring roles in defense from metabolic allocation of rice diterpenoids.Plant Cell301119–1131. 10.1105/tpc.18.00205
22
MafuS.DingY.MurphyK. M.YaacoobiO.AddisonJ. B.WangQ.et al (2018). Discovery, biosynthesis and stress-related accumulation of dolabradiene-derived defenses in maize.Plant Physiol.1762677–2690. 10.1104/pp.17.01351
23
MafuS.FischerE.AddissonB.Riberio BarbosanaI.ZerbeP. (2016). Substitution of two active site residues alters C9-hydroxylation in a class II diterpene synthase.ChemBioChem172304–2307. 10.1002/cbic.201600419
24
MannF. M.PrisicS.DavenportE. K.DetermanM. K.CoatesR. M.PetersR. J. (2010). A single residue switch for Mg(2 +)-dependent inhibition characterizes plant class II diterpene cyclases from primary and secondary metabolism.J. Biol. Chem.28520558–20563. 10.1074/jbc.M110.123307
25
MeyerJ.BergerD. K.ChristensenS. A.MurrayS. L. (2017). RNA-Seq analysis of resistant and susceptible sub-tropical maize lines reveals a role for kauralexins in resistance to grey leaf spot disease, caused by Cercospora zeina.BMC Plant Biol.17:197. 10.1186/s12870-017-1137-9
26
MorroneD.LowryL.DetermanM. K.HersheyD. M.XuM.PetersR. J. (2010). Increasing diterpene yield with a modular metabolic engineering system in E. coli: comparison of MEV and MEP isoprenoid precursor pathway engineering.Appl. Microbiol. Biotechnol.851893–1906. 10.1007/s00253-009-2219-x
27
MuellerD. W. (2012). Corn Disease Loss Estimates From the United States and Ontario, Canada—. (2012). Available at: http://fyi.uwex.edu/fieldcroppathology/files/2014/01/BP-96-12-W.pdf
28
PaterakiI.Andersen-RanbergJ.HambergerB.HeskesA. M.MartensH. J.ZerbeP.et al (2014). Manoyl oxide (13R), the biosynthetic precursor of forskolin, is synthesized in specialized root cork cells in Coleus forskohlii.Plant Physiol.1641222–1236. 10.1104/pp.113.228429
29
PelotK. A.ChenR.HagelthornD. M.YoungC. A.AddisonJ. B.MuchlinskiA.et al (2018). Functional diversity of diterpene synthases in the biofuel crop switchgrass.Plant Physiol.10.1104/pp.18.00590Epub ahead of print.
30
PelotK. A.MitchellR.KwonM.HagelthornD. M.WardmanJ. F.ChiangA.et al (2016). Biosynthesis of the psychotropic plant diterpene salvinorin A: discovery and characterization of the Salvia divinorum clerodienyl diphosphate synthase.Plant J.89885–897. 10.1111/tpj.13427
31
PetersR. J. (2006). Uncovering the complex metabolic network underlying diterpenoid phytoalexin biosynthesis in rice and other cereal crop plants.Phytochemistry672307–2317. 10.1016/j.phytochem.2006.08.009
32
PetersR. J. (2010). Two rings in them all: the labdane-related diterpenoids.Nat. Prod. Rep.271521–1530. 10.1039/c0np00019a
33
PotterK.CriswellJ.ZiJ. C.StubbsA.PetersR. J. (2014). Novel product chemistry from mechanistic analysis of ent-copalyl diphosphate synthases from plant hormone biosynthesis.Angew. Chem. Int. Ed. Engl.537198–7202. 10.1002/anie.201402911
34
PotterK. C.JiaM.HongY. J.TantilloD.PetersR. J. (2016a). Product rearrangement from altering a single residue in the rice syn-copalyl diphosphate synthase.Org. Lett.181060–1063. 10.1021/acs.orglett.6b00181
35
PotterK. C.ZiJ.HongY. J.SchulteS.MalchowB.TantilloD. J.et al (2016b). Blocking deprotonation with retention of aromaticity in a plant ent-Copalyl diphosphate synthase leads to product rearrangement.Angew. Chem. Int. Ed. Engl.55634–638. 10.1002/anie.201509060
36
RichardsD. E.KingK. E.Ait-AliT.HarberdN. P. (2001). How gibberellin regulates plant growth and development: a molecular genetic analysis of gibberellin signaling.Annu. Rev. Plant Physiol. Plant Mol. Biol.5267–88. 10.1146/annurev.arplant.52.1.67
37
SchmelzE. A.AlbornH. T.TumlinsonJ. H. (2001). The influence of intact-plant and excised-leaf bioassay designs on volicitin- and jasmonic acid-induced sesquiterpene volatile release in Zea mays.Planta214171–179. 10.1007/s004250100603
38
SchmelzE. A.HuffakerA.SimsJ. W.ChristensenS. A.LuX.OkadaK.et al (2014). Biosynthesis, elicitation and roles of monocot terpenoid phytoalexins.Plant J.79659–678. 10.1111/tpj.12436
39
SchmelzE. A.KaplanF.HuffakerA.DafoeN. J.VaughanM. M.NiX.et al (2011). Identity, regulation, and activity of inducible diterpenoid phytoalexins in maize.Proc. Natl. Acad. Sci. U.S.A.1085455–5460. 10.1073/pnas.1014714108
40
SchulteS.PotterK. C.LemkeC.PetersR. J. (2018). Catalytic bases and stereocontrol in lamiaceae class II diterpene cyclases.Biochemistry573473–3479. 10.1021/acs.biochem.8b00193
41
ThollD. (2015). Biosynthesis and biological functions of terpenoids in plants.Adv. Biochem. Eng. Biotechnol.14863–106. 10.1007/10_2014_295
42
ToyomasuT.UsuiM.SugawaraC.OtomoK.HiroseY.MiyaoA.et al (2014). Reverse-genetic approach to verify physiological roles of rice phytoalexins: characterization of a knockdown mutant of OsCPS4 phytoalexin biosynthetic gene in rice.Physiol. Plant15055–62. 10.1111/ppl.12066
43
VaughanM. M.ChristensenS.SchmelzE. A.HuffakerA.McauslaneH. J.AlbornH. T.et al (2015). Accumulation of terpenoid phytoalexins in maize roots is associated with drought tolerance.Plant Cell. Environ.382195–2207. 10.1111/pce.12482
44
WalleyJ. W.SartorR. C.ShenZ. X.SchmitzR. J.WuK. J.UrichM. A.et al (2016). Integration of omic networks in a developmental atlas of maize.Science353814–818. 10.1126/science.aag1125
45
WoutersF. C.BlanchetteB.GershenzonJ.VassaoD. G. (2016). Plant defense and herbivore counter-defense: benzoxazinoids and insect herbivores.Phytochem. Rev.151127–1151. 10.1007/s11101-016-9481-1
46
WuY.ZhouK.ToyomasuT.SugawaraC.OkuM.AbeS.et al (2012). Functional characterization of wheat copalyl diphosphate synthases sheds light on the early evolution of labdane-related diterpenoid metabolism in the cereals.Phytochemistry8440–46. 10.1016/j.phytochem.2012.08.022
47
XuM.HillwigM. L.PrisicS.CoatesR. M.PetersR. J. (2004). Functional identification of rice syn-copalyl diphosphate synthase and its role in initiating biosynthesis of diterpenoid phytoalexin/allelopathic natural products.Plant J.39309–318. 10.1111/j.1365-313X.2004.02137.x
48
XuM.WildermanP. R.MorroneD.XuJ.RoyA.Margis-PinheiroM.et al (2007a). Functional characterization of the rice kaurene synthase-like gene family.Phytochemistry68312–326. 10.1016/j.phytochem.2006.10.016
49
XuM.WildermanP. R.PetersR. J. (2007b). Following evolution’s lead to a single residue switch for diterpene synthase product outcome.Proc. Natl. Acad. Sci. U.S.A.1047397–7401. 10.1073/pnas.0611454104
50
ZerbeP.BohlmannJ. (2015). Plant diterpene synthases: exploring modularity and metabolic diversity for bioengineering.Trends Biotechnol.33419–428. 10.1016/j.tibtech.2015.04.006
51
ZerbeP.ChiangA.DullatH.O’neil-JohnsonM.StarksC.HambergerB.et al (2014). Diterpene synthases of the biosynthetic system of medicinally active diterpenoids in Marrubium vulgare.Plant J.79914–927. 10.1111/tpj.12589
52
ZerbeP.ChiangA.YuenM.HambergerB.HambergerB.DraperJ. A.et al (2012). Bifunctional cis-abienol synthase from Abies balsamea discovered by transcriptome sequencing and its implications for diterpenoid fragrance production.J. Biol. Chem.28712121–12131. 10.1074/jbc.M111.317669
53
ZerbeP.HambergerB.YuenM. M.ChiangA.SandhuH. K.MadilaoL. L.et al (2013). Gene discovery of modular diterpene metabolism in nonmodel systems.Plant Physiol.1621073–1091. 10.1104/pp.113.218347
54
ZerbeP.RodriguezS. M.MafuS.ChiangA.SandhuH. K.O’neil-JohnsonM.et al (2015). Exploring diterpene metabolism in non-model species: transcriptome-enabled discovery and functional characterization of labda-7,13E-dienyl diphosphate synthase from Grindelia robusta.Plant J.83783–793. 10.1111/tpj.12925
55
ZhouK.XuM.TiernanM.XieQ.ToyomasuT.SugawaraC.et al (2012). Functional characterization of wheat ent-kaurene(-like) synthases indicates continuing evolution of labdane-related diterpenoid metabolism in the cereals.Phytochemistry8447–55. 10.1016/j.phytochem.2012.08.021
56
ZiJ.MafuS.PetersR. J. (2014). To gibberellins and beyond! Surveying the evolution of (di)terpenoid metabolism.Annu. Rev. Plant. Biol.65259–286. 10.1146/annurev-arplant-050213-035705
Summary
Keywords
diterpene synthase, diterpenoid biosynthesis, plant stress response, plant specialized metabolism, Zea mays
Citation
Murphy KM, Ma L-T, Ding Y, Schmelz EA and Zerbe P (2018) Functional Characterization of Two Class II Diterpene Synthases Indicates Additional Specialized Diterpenoid Pathways in Maize (Zea mays). Front. Plant Sci. 9:1542. doi: 10.3389/fpls.2018.01542
Received
21 August 2018
Accepted
02 October 2018
Published
23 October 2018
Volume
9 - 2018
Edited by
Reuben J. Peters, Iowa State University, United States
Reviewed by
Irini Pateraki, University of Copenhagen, Denmark; Courtney M. Starks, Sequoia Sciences (United States), United States
Updates
Copyright
© 2018 Murphy, Ma, Ding, Schmelz and Zerbe.
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: Philipp Zerbe, pzerbe@ucdavis.edu
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.