-
PDF
- Split View
-
Views
-
Cite
Cite
Putri Pratiwi, Genta Tanaka, Tomohiro Takahashi, Xiaonan Xie, Koichi Yoneyama, Hideyuki Matsuura, Kosaku Takahashi, Identification of Jasmonic Acid and Jasmonoyl-Isoleucine, and Characterization of AOS, AOC, OPR and JAR1 in the Model Lycophyte Selaginella moellendorffii, Plant and Cell Physiology, Volume 58, Issue 4, April 2017, Pages 789–801, https://doi.org/10.1093/pcp/pcx031
- Share Icon Share
Jasmonic acid (JA) is involved in a variety of physiological responses in seed plants. However, the detection and role of JA in lycophytes, a group of seedless vascular plants, have remained elusive until recently. This study provides the first evidence of 12-oxo-phytodienoic acid (OPDA), JA and jasmonoyl-isoleucine (JA-Ile) in the model lycophyte Selaginella moellendorffii. Mechanical wounding stimulated the accumulation of OPDA, JA and JA-Ile. These data were corroborated by the detection of enzymatically active allene oxide synthase (AOS), allene oxide cyclase (AOC), 12-oxo-phytodienoic acid reductase 3 (OPR3) and JA-Ile synthase (JAR1) in S. moellendorffii. SmAOS2 is involved in the first committed step of JA biosynthesis. SmAOC1 is a crucial enzyme for generating the basic structure of jasmonates and is actively involved in the formation of OPDA. SmOPR5, a functionally active OPR3-like enzyme, is also vital for the reduction of (+)-cis-OPDA, the only isomer of the JA precursor. The conjugation of JA to Ile by SmJAR1 demonstrates that S. moellendorffii produces JA-Ile. Thus, the four active enzymes have characteristics similar to those in seed plants. Wounding and JA treatment induced the expression of SmAOC1 and SmOPR5. Furthermore, JA inhibited the growth of shoots in S. moellendorffii, which suggests that JA functions as a signaling molecule in S. moellendorffii. This study proposes that JA evolved as a plant hormone for stress adaptation, beginning with the emergence of vascular plants.
Introduction
As sessile organisms, plants must adapt to many environmental stresses that can negatively affect their development, growth, productivity and fertility by altering their physiological, morphological and/or developmental processes. Many studies have found that JA modulates numerous plant physiological processes that are related to development and defense responses (Wasternack and Parthier 1997, Ziegler et al. 2000, Ishiguro et al. 2001, Wasternack 2007). Indeed, JA is a naturally occurring phytohormone that is ubiquitous in seed plant species (Creelman and Mullet 1997).

Vascular plants have survived on earth for the past 420 million years and have diverged into several descendants. Only two descendants have survived: the euphyllophytes and lycophytes (Kenrick and Crane 1997, Banks et al. 2011). Fossil evidence shows that lycophytes diverged shortly after land plants evolved vascular tissues, which generated all other vascular plants. There are currently only three extant lycophyte groups: the Lycopodiales (club mosses), Isoetales (quillworts) and Selaginellales (spike mosses) (Weng et al. 2005, Banks 2009, Hecht et al. 2011). Recently, the approximately 100 Mbp genome of Selaginella moellendorffii was fully sequenced (Banks et al. 2011). This model lycophyte, as a representative of the earliest and still surviving vascular plant lineage, can potentially provide information for deciphering the evolution of the physiological, biochemical and developmental processes that are unique to land plants.
The model plants that transitioned from water to land, such as the moss Physcomitrella patens and the liverwort Marchantia polymorpha, produce OPDA but not JA (Stumpe et al. 2010, Yamamoto et al. 2015). These findings strongly suggest that the first half of the octadecanoid pathway in chloroplasts remains in bryophytes. Thereafter, the difference in JA production between seed plants and bryophytes is one of the significant observations in the study of plant evolution. The lycophytes are intermediate plants between bryophytes and seed plants, and this ancient lineage diverged shortly after land plants evolved vascular tissues. Thus, the origin of JA is of interest for studying plant evolution. In this study, we analyzed the lycophyte model S. moellendorffii to determine the presence of OPDA, JA and JA-Ile, as well as the enzymatic activities of AOS, AOC, OPR3 and JAR1, which are involved in the biosynthesis of jasmonates. Moreover, we showed the growth-inhibitory activity of OPDA, JA and JA-Ile in S. moellendorffii. These data indicate a possible role for JA as a signaling molecule that regulates growth and response to wounding in S. moellendorffii.
Results
Identification of OPDA, JA and JA-Ile in S. moellendorffii
![UPLC-MS/MS analysis of OPDA, JA and JA-Ile in S. moellendorffii. (A−C) MRM mode analysis of a specific daughter peak at m/z 164.99 that was derived from the peak at m/z 291.37 [M−H]−. (D−F) MRM mode analysis of a specific daughter peak at m/z 58.71 that was derived from the peak at m/z 209.09 [M−H]−. (G−I) MRM mode analysis of a specific daughter peak at m/z 129.68 that was derived from the peak at m/z 322.03 [M−H]−. (A, D, G) Standard. (B, E, H) S. moellendorffii extract. (C, F, I) Co-injection of the standard and S. moellendorffii extract.](https://oup.silverchair-cdn.com/oup/backfile/Content_public/Journal/pcp/58/4/10.1093_pcp_pcx031/1/m_pcx031f2.jpeg?Expires=1748011145&Signature=0HbPsG7f2HSddijWqT3E99lPkJgWDVhcG1-s-WjEFMHNfONz412-ZU3T9Uj-A5XY4jPdpBb1LyQ7Lf20Y-ihMu17JA4uH9ZsBBheKBr9oNUO4iN8NQkffKw4oppGYqDPZpXzSpOI3pg2by1Jvy8Lij8YZWDZc6wkeX25N0ijrcMFD4827-WItx0UxKtqSdn-X6VPhyJYWYpHckPt6~xA99FrKqd512Uf2LraNxHInx2uc8B1cNe7ls4PaDdE3Z-iMS58pdlrR9d2rg5Y-NOw9eDPGRtFp4Zf90xO-7K38XwipEOc0xPVn9728dTT6kU6bv4YqBk7n0bMN7IVsIyuSg__&Key-Pair-Id=APKAIE5G5CRDK6RD3PGA)
UPLC-MS/MS analysis of OPDA, JA and JA-Ile in S. moellendorffii. (A−C) MRM mode analysis of a specific daughter peak at m/z 164.99 that was derived from the peak at m/z 291.37 [M−H]−. (D−F) MRM mode analysis of a specific daughter peak at m/z 58.71 that was derived from the peak at m/z 209.09 [M−H]−. (G−I) MRM mode analysis of a specific daughter peak at m/z 129.68 that was derived from the peak at m/z 322.03 [M−H]−. (A, D, G) Standard. (B, E, H) S. moellendorffii extract. (C, F, I) Co-injection of the standard and S. moellendorffii extract.
Accumulation of OPDA, JA and JA-Ile upon wounding stress

Accumulation of OPDA, JA and OPDA in S. moellendorffii after wounding. S. moellendorffii was subjected to mechanical wounding and then collected at the indicated time after wounding. The concentrations of OPDA, JA and JA-Ile were analyzed by UPLC-MS/MS. The data are presented as the mean ± SD (n = 3). The asterisks represent a significant difference between the treated and control plants (Student’s t-test, *P < 0.05, **P < 0.01).
Putative AOS, AOC, OPR3 and JAR1 genes in S. moellendorffii
To confirm further the presence of OPDA, JA and JA-Ile in S. moellendorffii, we investigated the characteristics of important JA biosynthetic enzymes in S. moellendorffii. A genomic database analysis using Phytozome v 11.0 (https://phytozome.jgi.doe.gov/pz/portal.html) was performed to screen AOS, AOC, OPR3 and JAR1 homologous genes in S. moellendorffii using Arabidopsis amino acid sequences corresponding to these proteins as queries. As a result, the presence of putative AOS (SmAOS1, Sm_271334; SmAOS2, Sm_177201; SmAOS3, Sm_228572), AOC (SmAOC1, Sm_91887), OPR3 (SmOPR1, Sm_270843; SmOPR5, Sm_111662) and JAR1 (SmJAR1, Sm_110439) genes was predicted in S. moellendorffii (Supplementary Table S1).
Functional analysis of SmAOS2
The first committed step of JA biosynthesis is catalyzed by AOS, which converts 13(S)-HPOT into unstable 12,13-EOT. A genomic database analysis of Phytozome v 11.0 revealed the presence of nine putative AOS genes in S. moellendorffii. Amino acid sequence analysis of the SmAOS candidates showed that SmAOS1, SmAOS2 and SmAOS3 were promising candidates containing the important sequence motifs: the I-helix GXXX (F/L), EXLR motif and heme-binding PXVXNKQCPG, which are characteristic of members of the CYP74 family (Supplementary Fig. S1) (Koeduka et al. 2015). Additionally, the amino acid sequences of SmAOS1, SmAOS2 and SmAOS3 are quite similar (Supplementary Fig. S1). A phylogenetic tree of the AOSs showed that the SmAOSs were separated from the seed plant AOSs (Supplementary Fig. S2). Computational analysis using ChloroP v 1.1 (http://www.cbs.dtu.dk/services/ChloroP/) and iPSORT (http://ipsort.hgc.jp/) for subcellular localization predicted the absence of a transit peptide in SmAOS1, SmAOS2 and SmAOS3, whereas most known AOSs are localized in chloroplasts (Schaller and Stintzi 2009).

Chiral GC-MS analysis of the OPDA methyl ester. Recombinant SmAOS2 was reacted with the substrate 13(S)-HPOT, directly treated with alkaline solution and then methylated by ethereal diazomethane. The methylated OPDA was analyzed by GC-MS using a chiral capillary column. The data are presented as the total ion current chromatogram. The former and latter peaks represent the trans-(+)- and (−)-OPDA methyl esters, respectively (Laudert et al. 1997).
Functional analysis of SmAOC1
AOC has been successfully cloned from seed plants, a moss and a liverwort (Ziegler et al. 2000, Agrawal et al. 2003, Stenzel et al. 2003, Maucher et al. 2004, Farmaki et al. 2007, Pi et al. 2008, Stumpe et al. 2010, Hashimoto et al. 2011, Yamamoto et al. 2015). The 5′-untranslated region (UTR) of SmAOC1 was examined by 5′-rapid amplification of cDNA ends (RACE) to obtain the full-length sequence. A computational analysis using ChloroP v 1.1 and iPSORT predicted that SmAOC1 has a chloroplast transit peptide at its N-terminus. The amino acid sequence alignment of SmAOC1 and AOCs of other plants showed that SmAOC1 had high similarity to other AOCs (Supplementary Fig. S7). The phylogenetic tree of AOCs demonstrated that SmAOC1 was related to MpAOC from the liverwort M. polymorpha and to AOCs in seed plants, while it was separated from PpAOC1 and PpAOC2 of the moss P. patens (Supplementary Fig. S8).
To determine whether SmAOC1 is involved in the production of OPDA, SmAOC1 was cloned and overexpressed in E. coli. The recombinant SmAOC1 was fused with a His-tag at the N-terminus to replace the chloroplast signal peptide. After purifying SmAOC1 using Ni-NTA affinity column chromatography, a protein with the expected size of 22 kDa was clearly detected as a single band by SDS–PAGE analysis (Supplementary Fig. S3). The purified recombinant SmAOC1 was used for enzymatic analysis.

Chiral GC-MS analysis of the OPDA methyl ester. Recombinant SmAOC1 was reacted with the substrate 13-HPOT and an allene oxide synthase (PpAOS1), directly treated with alkaline solution and then methylated by ethereal diazomethane. The methylated OPDA was analyzed by GC-MS using a chiral capillary column. The molecular ion peak of the OPDA methyl ester at m/z 306 was monitored. The former and latter peaks represent the trans-(+)- and (−)-OPDA methyl ester, respectively (Laudert et al. 1997). (A) Racemic OPDA standard. (B) Product of the SmAOC1 enzymatic reaction.
The first of three enzymes involved in JA biosynthesis (LOX, AOS and AOC) have been reported to reside in chloroplasts (Schaller and Stintzi 2009). To examine whether SmAOC1 is localized in chloroplasts, a 35S::SmAOC1-GFP plasmid was constructed and introduced into P. patens protoplasts and was then finally observed by confocal laser scanning microscopy. The green fluorescent signal of SmAOC1 fused with green fluorescent protein (GFP) was clearly identified by the red autofluorescence of Chl. The overlay data also strongly supported that SmAOC1 is located in chloroplasts, which is similar to the other characterized AOCs (Supplementary Fig. S9). Accordingly, the chloroplast could be the functional organelle in the biosynthesis of OPDA in S. moellendorffii.
Functional analysis of SmOPRs
OPRs are classified into two groups based on their ability to reduce OPDA: OPR3-like enzymes (e.g. Arabidopsis AtOPR3, tomato SlOPR3 and rice OsOPR7) and OPR1-like enzymes (e.g. AtOPR1 and SlOPR1) (Breithaupt et al. 2009). OPR3-like enzymes catalyze the reduction of cis-(+)-OPDA, a natural JA precursor, and of cis-(−)-OPDA, which is not involved in JA biosynthesis. The Arabidopsis opr3 mutant shows male sterility and is desensitized to responses involving JA signaling. In contrast, OPR1-like enzymes preferentially reduce cis-(−)-OPDA rather than cis-(+)-OPDA; therefore, OPR1-like enzymes might function in pathways other than JA biosynthesis.
A genomic database search using Phytozome v 11.0 predicted the presence of six putative OPR genes in S. moellendorffii. All of the amino acid sequences of candidate SmOPRs were aligned with those of other known OPRs (Supplementary Fig. S10). Based on a comparison of two critical amino acids for substrate binding and a phylogenetic analysis of SmOPRs (Supplementary Fig. S10), two genes designated SmOPR1 and SmOPR5 were highly promising candidate genes that encoded OPR3-like enzymes involved in JA biosynthesis. SmOPR1 harbors two important active-site residues, phenylalanine and histidine, as an OPR3 motif in the active site, which is considered to be necessary to reduce the JA precursor cis-(+)-OPDA. However, SmOPR1 was grouped with the other SmOPRs in an independent OPR cluster (Supplementary Fig. S11). In contrast, SmOPR5 was a member of the cluster of OPR3-like enzymes, such as AtOPR3, OsOPR7, SlOPR3, ZmOPR7 and ZmOPR8, which are active in JA biosynthesis (Supplementary Fig. S11). Moreover, SmOPR5 shared 53% amino acid sequence identity with AtOPR3. Nonetheless, instead of phenylalanine and histidine, which are considered unique to the OPR3 substrate filter, the active site residues of SmOPR5 were tryptophan and histidine. Therefore, these two genes were selected for further analysis to verify their enzymatic activity.

Chiral GC-MS analysis of the reaction products catalyzed by SmOPR1 and SmOPR5. Recombinant SmOPR1 and SmOPR5 were each reacted with racemic OPDA, directly treated with alkaline solution and then methylated by ethereal diazomethane. The methylated reaction products were analyzed by GC-MS using a chiral capillary column. The data are presented as the total ion current chromatogram. Peaks 1–4 represent trans-isomers of (+)-OPC-8:0 (natural type), (−)-OPC-8:0, (+)-OPDA (natural type) and (−)-OPDA, respectively, that were derived from cis-isomers of these compounds.
Enzymatic production of JA-Ile by SmJAR1
JA-Ile is a versatile signaling molecule for JA-mediated defense and developmental events and is produced by the conjugation of JA and Ile. JAR1s are GH3 proteins that have the ability to synthesize JA-Ile (Staswick and Tiryaki 2004). An Arabidopsis jar1 mutant, which is unable to produce JA-Ile, fails to trigger JA-mediated responses (Suza and Staswick 2008). Thus, JAR1 activity is necessary for JA signal transduction in seed plants. Here, we found 14 JAR1 homologous genes in the Phytozome database using AtJAR1 as a query (Supplementary Table S1). Multiple sequence alignment and phylogenetic analysis using Clustal Omega software were performed to show the relationship between these JAR1 homologous proteins in S. moellendorffii and the other existing JAR1 proteins (Supplementary Figs. S14, S15). As a result, the candidate Sm_110439, which was designated SmJAR1, was most closely related to AtJAR1, OsJAR1 and OsJAR2, which can synthesize JA-Ile (Staswick and Tiryaki 2004, Wakuta et al. 2011). Moreover, SmJAR1 also has three short conserved motifs involved in ATP/AMP binding (Chang et al. 1997): SSGTSQGRPK (motif 1), YGSSE (motif 2) and YRLGD (motif 3) (Supplementary Fig. S14). To determine whether SmJAR1 encodes functional JA-Ile synthase in S. moellendorffii, SmJAR1 was isolated to check its enzymatic activity.
![UPLC-MS/MS analysis of JA-Ile in E. coli overexpressing SmJAR1. MRM mode analysis of a specific daughter peak of JA-Ile-d6 at m/z 129.68 that was derived from the peak at m/z 328.03 [M−H]− and of the peak at m/z 129.68 that was derived from the peak of JA-Ile at m/z 322.03 [M−H]−. (A) Internal standard JA-Ile-d6. (B) Culture supernatant of E. coli overexpressing SmJAR1. (C) Culture supernatant of E. coli containing pET23a as a control. *The value of 100% is equal to 9 × 105 counts, which was converted by the algorithm in Mass Lynx software (Waters Corporation).](https://oup.silverchair-cdn.com/oup/backfile/Content_public/Journal/pcp/58/4/10.1093_pcp_pcx031/1/m_pcx031f7.jpeg?Expires=1748011146&Signature=WMNfCZJ1IbNf3RKk2twAld8Rttmvoex~r4Oa6TcN9NQqFVfDcCVrmFYZqZgvE~1IHhD6fJaeDHkDVUBPN~PvKUINLvhbfRvWeMuyxOUgOeHEEaOx0IdIoMEGO2fvyhPZcWpKd8VurwQzq-Wp3sKZhxPfj7MCXDXQqKKNti6TWQIAGc94AZtr7XL7ADqXOB3IN-F-oYerKabHf4qL16E1gHiDEm5RB6V1UHdlYMqDmQGkSrSplUNvvzdkfsyNP1ct993hXCBi2BtfliHRPV18V5LE0ddDStX71RvLCzaXSwpJ6bsCjIUySdB7onzZ7oWCy2Tw3Gj4Yd3v~xDy4-la6w__&Key-Pair-Id=APKAIE5G5CRDK6RD3PGA)
UPLC-MS/MS analysis of JA-Ile in E. coli overexpressing SmJAR1. MRM mode analysis of a specific daughter peak of JA-Ile-d6 at m/z 129.68 that was derived from the peak at m/z 328.03 [M−H]− and of the peak at m/z 129.68 that was derived from the peak of JA-Ile at m/z 322.03 [M−H]−. (A) Internal standard JA-Ile-d6. (B) Culture supernatant of E. coli overexpressing SmJAR1. (C) Culture supernatant of E. coli containing pET23a as a control. *The value of 100% is equal to 9 × 105 counts, which was converted by the algorithm in Mass Lynx software (Waters Corporation).
Expression of genes encoding key enzymes in JA-Ile biosynthesis under wounding stress and JA treatment

Quantitative RT-PCR analysis of SmAOC1, SmOPR5 and SmJAR1 in wounded S. moellendorffii. Plants were subjected to mechanical wounding and then collected at the indicated times after wounding. To evaluate the transcript levels of SmAOC1, SmOPR5 and SmJAR1, qRT-PCR was performed using RNA isolated from control and wounded plants. The data are presented as the mean ± SD (n = 3). The asterisks represent a significant difference between the treated and control plants (Student’s t-test, *P < 0.05; **P < 0.01).

Quantitative RT-PCR analysis of SmAOC1, SmOPR5 and SmJAR1 in S. moellendorffii treated with JA. Plants were treated with 50 µM JA and then collected at the indicated times after JA treatment. To evaluate the transcript levels of SmAOC1, SmOPR5 and SmJAR1, qRT-PCR was performed using RNA isolated from control and JA-treated plants. The data are presented as the mean ± SD (n = 3). The asterisks represent a significant difference between the treated and control plants (Student’s t-test, *P < 0.05; **P < 0.01; ***P < 0.001).
Effect of exogenous OPDA, JA and JA-Ile on the growth of S. moellendorffii shoots

Growth-inhibitory activities of OPDA, JA and JA-Ile in S. moellendorffii. S. moellendorffii was treated with OPDA, JA and JA-Ile for 4 weeks, and then the shoots were measured. The data are presented as the mean ± SD (n = 5). The asterisks represent a significant difference between the treated and control plants (Student’s t-test, **P < 0.01; ***P < 0.001).
Discussion
Identification and accumulation of OPDA, JA and JA-Ile
In this study, the presence of OPDA, JA and JA-Ile in S. moellendorffii was successfully demonstrated (Fig. 2). Moreover, wounding transiently increased the endogenous concentrations of OPDA and JA in S. moellendorffii within 10 min (Fig. 3), which is probably due to the role of jasmonates in regulating responses to environmental stress in S. moellendorffii. Given that wounding also stimulated OPDA and JA in the fern Pteridium aquilinum (Radhika et al. 2012), the accumulation of JA in response to wounding is probably a common physiological response among all vascular plant species (Maucher et al. 2004). However, JA-Ile accumulated more slowly than OPDA and JA after wounding. It is possible that the early rise in JA levels within 30 min represents an early response, whereas the slow JA-Ile accumulation may be attributable to a late stress response in this plant.
The transient increases in JA and OPDA after wounding observed in S. moellendorffii are similar to those observed in seed plants. It has been suggested that another metabolic pathway is required to inactivate JA responses in S. moellendorffii. In Arabidopsis, the oxidation of JA-Ile, which decreases the JA-Ile concentration, is important for regulating JA activity (Kitaoka et al. 2011, Koo et al. 2011). A phylogenetic analysis of Cyt P450s in Arabidopsis thaliana, P. patens and S. moellendorffii revealed that SmCYP94Js (SmCYP94J1, SmCYP94J2, SmCYP94J3, SmCYP94J5Pv1 and SmCYP94J5Pv2) are members of a cluster that includes AtCYP94B3 and AtCYP94C1 as enzymes responsible for the oxidation of JA-Ile and 12OH-JA-Ile (Banks et al. 2011). SmCYP94Js are likely to be related to the oxidation of JA-Ile, leading to the inactivation of JA activity.
In studies of P. patens and M. polymorpha, only OPDA has been detected in both plants (Stumpe et al. 2010, Yamamoto et al. 2015). In contrast to vascular plants, JA is not important for the physiology of these model bryophytes. Hence, our data suggest that JA and JA-Ile biosynthesis first appeared after bryophytes in plant evolution (Figs. 2, 3). Additionally, there is a significant difference between bryophytes and lycophytes. Lycophytes have a vascular system for the efficient transport of water, nutrients and molecules, as well as for enhancing plant height and size. It is likely that the emergence of the JA biosynthetic pathway after OPDA is related to the acquisition by plants of a vascular system.
Enzymatic activity
There have been no previous studies on enzymes in the JA biosynthetic pathway of non-seed vascular plants. To support the analytical data showing the presence of OPDA, JA and JA-Ile, the enzymatic activities of AOS, AOC, OPR3 and JAR1 were investigated. Here, genes encoding functional AOS, AOC, OPR3 and JAR1 proteins (SmAOS2, SmAOC1, SmOPR1 and SmOPR5 and SmJAR1, respectively) were identified in S. moellendorffii.
SmAOS2. The AOS reaction is the first committed step of JA biosynthesis in plants. AOS is classified as a member of the CYP74 family and requires an oxygenated fatty acid hydroperoxide substrate instead of oxygen or a redox partner (Howe and Schilmiller 2002, Werck-Reichhart et al. 2002). This study revealed three AOS candidate genes (SmAOS1, SmAOS2 and SmAOS3) in S. moellendorffii (Supplementary Figs. S1, S2). Among them, recombinant SmAOS2, which was successfully produced in E. coli, exhibited AOS activity similar to that of previously reported AOSs (Fig. 4; Supplementary Figs. S4−S6). Computer programs analyzing protein subcellular localization did not predict a transit peptide in the N-terminus of SmAOS2. Because most known AOSs are localized in chloroplasts, SmAOS2 is probably located in chloroplasts for JA biosynthesis. Since the amino acid sequences of SmAOS1, SmAOS2 and SmAOS3 are similar, further investigation is needed to clarify the physiological functions of each of these SmAOSs.
SmAOC1. AOC is a critical enzyme that establishes the enantiomeric structure of OPDA, which contributes to the basic structure of jasmonates. The enzymatic activity of SmAOC1 is required for the synthesis of naturally occurring OPDA (Fig. 5). The stereochemistry of OPDA in S. moellendorffii is the same as that in other studied plants. Our microscopic observations using SmAOC1 fused to GFP showed that SmAOC1 localized to chloroplasts (Supplementary Fig. S9). These results suggest that the first half of the octadecanoid pathway for synthesizing OPDA is present in the chloroplasts of S. moellendorffii. The presence of LOX, AOS and AOC in chloroplasts may have been conserved during the process of land plant evolution.
SmOPR1 and SmOPR5. The enzymatic analysis of SmOPR1 and SmOPR5 revealed that S. moellendorffii possesses two types of OPRs that are classified by substrate preference: OPR1-like enzymes and OPR3-like enzymes. SmOPR1 and SmOPR5 are similar in their molecular mass and isoelectric point (pI) (Table 1) but differ in their substrate specificity for cis-OPDA isomers. Of these two SmOPRs, only SmOPR5 reduced cis-(±)-OPDA, including the endogenous substrate cis-(+)-OPDA, which is the natural JA precursor (Fig. 6; Supplementary Figs. S12, S13). In contrast to SmOPR5, SmOPR1 reduced only cis-(−)-OPDA, which is the unnatural type (Fig. 6). SmOPR1 might instead be required for the oxylipin metabolic pathway.
Gene name . | ORF length (bp) . | Protein length (No. of amino acids) . | Mol. wt. (Da) . | Theoretical pI . |
---|---|---|---|---|
SmAOS1 | 1,458 | 485 | 54,107 | 6.01 |
SmAOS2 | 1,404 | 467 | 52,079 | 6.35 |
SmAOS3 | 1,458 | 485 | 54,071 | 5.89 |
SmAOC1 | 738 | 245 | 26,364 | 9.15 |
SmOPR1a | 1131 | 376 | 41,588 | 6.28 |
SmOPR5a | 1209 | 402 | 44,046 | 6.05 |
SmJAR1 | 1749 | 582 | 64,751 | 5.83 |
Gene name . | ORF length (bp) . | Protein length (No. of amino acids) . | Mol. wt. (Da) . | Theoretical pI . |
---|---|---|---|---|
SmAOS1 | 1,458 | 485 | 54,107 | 6.01 |
SmAOS2 | 1,404 | 467 | 52,079 | 6.35 |
SmAOS3 | 1,458 | 485 | 54,071 | 5.89 |
SmAOC1 | 738 | 245 | 26,364 | 9.15 |
SmOPR1a | 1131 | 376 | 41,588 | 6.28 |
SmOPR5a | 1209 | 402 | 44,046 | 6.05 |
SmJAR1 | 1749 | 582 | 64,751 | 5.83 |
aThe numbering of SmOPR genes was previously described by Li et al. (2009).
Gene name . | ORF length (bp) . | Protein length (No. of amino acids) . | Mol. wt. (Da) . | Theoretical pI . |
---|---|---|---|---|
SmAOS1 | 1,458 | 485 | 54,107 | 6.01 |
SmAOS2 | 1,404 | 467 | 52,079 | 6.35 |
SmAOS3 | 1,458 | 485 | 54,071 | 5.89 |
SmAOC1 | 738 | 245 | 26,364 | 9.15 |
SmOPR1a | 1131 | 376 | 41,588 | 6.28 |
SmOPR5a | 1209 | 402 | 44,046 | 6.05 |
SmJAR1 | 1749 | 582 | 64,751 | 5.83 |
Gene name . | ORF length (bp) . | Protein length (No. of amino acids) . | Mol. wt. (Da) . | Theoretical pI . |
---|---|---|---|---|
SmAOS1 | 1,458 | 485 | 54,107 | 6.01 |
SmAOS2 | 1,404 | 467 | 52,079 | 6.35 |
SmAOS3 | 1,458 | 485 | 54,071 | 5.89 |
SmAOC1 | 738 | 245 | 26,364 | 9.15 |
SmOPR1a | 1131 | 376 | 41,588 | 6.28 |
SmOPR5a | 1209 | 402 | 44,046 | 6.05 |
SmJAR1 | 1749 | 582 | 64,751 | 5.83 |
aThe numbering of SmOPR genes was previously described by Li et al. (2009).
A structural comparison of tomato SlOPR1 and SlOPR3 indicated that two active site residues, Tyr78 and Tyr246 in SlOPR1 and Phe74 and His244 in SlOPR3, are critical for substrate specificity (Breithaupt et al. 2009). SlOPR3 is less enantioselective due to its two relatively smaller amino acid residues that form a larger substrate-binding pocket. In contrast, the relatively larger amino acids (two tyrosine residues) in SlOPR1 permit access to only cis-(−)-OPDA, the unnatural type, at the substrate-binding site. This study showed that SmOPR5, but not SmOPR1, actively converts the endogenous substrate cis-(+)-OPDA to cis-(+)-OPC-8:0. This result is in contrast to the prediction from the substrate preference analysis of their active site residues but agrees with the phylogenetic analysis of OPRs (Supplementary Figs. S10, S11). Mutation of the important amino acids phenylalanine and histidine into tyrosine in SlOPR3 did not completely change the substrate specificity (Breithaupt et al. 2009). There is an unidentified factor in the structure of OPRs that is important for substrate specificity, and this factor must be identified by further crystal structure analysis.
Li et al. (2009) showed that OPRs from the lower land plants P. patens (PpOPR1, PpOPR2, PpOPR4 and PpOPR5) and S. moellendorffii (SmOPR1, SmOPR2, SmOPR3, SmOPR4 and SmOPR6) clustered together in subgroup VI. Because SmOPR1 showed OPR1 activity, the other OPRs in subgroup VI may be unrelated to JA biosynthesis. In contrast, the other known OPR3-like enzymes are positioned in the same cluster in subgroup II. It is most likely that OPR3-like enzymes independently evolved to expand the substrate-binding pocket to accept cis-(+)-OPDA.
SmJAR1. JAR1 is responsible for conjugating JA to amino acids. Because JA-Ile is a versatile molecule in JA signaling, JAR1 is necessary for JA-mediated physiological events in seed plants (Wasternack and Hause 2013). JAR1 is a member of the GH3 proteins, which are classified into three groups based on their enzymatic activity. Group I proteins, such as Arabidopsis AtJAR1, can synthesize the JA–amino acid conjugates JA-Ile, JA-Leu and JA-Val (Staswick and Tiryaki 2004). Group II proteins conjugate IAA and salicylic acid to various amino acids (Staswick et al. 2005). Group III proteins accept benzoates as substrates for amino acid conjugation (Okrent et al. 2009). A group I GH3 protein from seed plants showed JA-Ile synthase activity, which functions in defense (Kang et al. 2006, Wang et al. 2008, Suza et al. 2010).
Selaginella moellendorffii produces JA-Ile based on evidence of a functionally active SmJAR1 (Fig. 7). The identification of JA-Ile and an active SmJAR1 strongly suggests that JA-Ile-mediated signaling is functional in S. moellendorffii. Moreover, SmJAR1 was hypothesized to prefer isoleucine as a substrate for conjugation with JA (Fig. 7; Supplementary Fig. S16). The same characteristic was also found in AtJAR1 (Staswick and Tiryaki 2004, Suza and Staswick 2008). Considering the substrate specificity of SmJAR1, JA-Ile seems to be more important than the other JA–amino acid conjugates. Additionally, endogenous conjugates of JA with another amino acid, such as JA-Trp, JA-Phe, JA-Gly, JA-Leu and JA-Val, have not been found in S. moellendorffii (Supplementary Fig. S19).
In the case of rice (Oryza sativa), OsJAR2 in the cluster of group II GH3 proteins can conjugate JA to Ile, similar to OsJAR1 in the cluster of group I GH3 proteins (Wakuta et al. 2011). The other 13 JAR1 homologs found in S. moellendorffii might also show JA-Ile synthetic activity. Given that GH3 proteins synthesize conjugates of JA and amino acids in P. patens (Ludwig-Müller et al. 2009), land plants are suggested to have enzymes with ubiquitous JAR1 activity. JA synthesis and JAR1 activity might have evolved independently in land plants.
Expression of SmAOC1, SmOPR5 and SmJAR1
Wounding and JA induce the expression of various genes involved in defense responses in seed plants, i.e. genes encoding enzymes in JA biosynthesis. In S. moellendorffii, the expression of SmAOC1 and SmOPR5 under wounding treatment was increased (Fig. 8). These results were consistent with the transient increases in OPDA and JA in S. moellendorffii (Fig. 3). Based on the SmOPR5 expression level, the accumulation of SmOPR5 appears to precede the production of OPDA. Additionally, the late accumulation of JA-Ile coincided with the late expression level of SmJAR1 after wounding (Fig. 8). The expression profiles of SmAOC1, SmOPR5 and SmJAR1 support the accumulation patterns of OPDA, JA and JA-Ile, probably due to the role of jasmonates in regulating responses to environmental stress in S. moellendorffii.
The present study also demonstrated that JA induced the expression of SmAOC1 and SmOPR5 in S. moellendorffii. JA treatment led to the accumulation of genes encoding enzymes in the octadecanoid pathway, indicating positive feedback regulation of JA biosynthesis. The expression of SmAOC1 and SmOPR5 was enhanced upon JA application, and these expression patterns were similar to those observed under wounding treatment (Fig. 9). The JA-induced expression of SmAOC1 and SmOPR5 allows us to infer the presence of positive feedback regulation of JA biosynthesis in S. moellendorffii. In contrast to the up-regulation of SmJAR1 by wounding, JA inhibited SmJAR1 expression at the transcriptional level (Fig. 9). Although the detailed JA signaling pathway has not been elucidated in S. moellendorffii, the suppression of SmJAR1 expression by JA appears to function as negative feedback regulation of JA-Ile-mediated signaling and/or an increase in JA signaling that is independent from JA-Ile. Wounding probably activates various signaling pathways that are JA dependent and JA independent. In the case of SmJAR1 expression, integrated signaling caused by wounding might induce SmJAR1 expression against stresses.
OPDA, JA and JA-Ile inhibit the growth of S. moellendorffii
Physcomitrella patens and M. polymorpha cannot produce JA; instead, OPDA inhibits their growth, and JA does not induce any significant physiological responses (Stumpe et al. 2010, Yamamoto et al. 2015). However, P. patens responds to methyl jasmonate, which has not been detected in P. patens, by reducing its colony size (Ponce de León et al. 2012). Many studies have reported the growth-inhibitory activity of JA in seed plants. OPDA, JA and JA-Ile treatment retarded the growth of S. moellendorffii in a manner similar to that of seed plants (Fig. 10; Supplementary Fig. S18). This report is the first evidence that OPDA, JA and JA-Ile inhibit the growth of non-seed vascular plants. The growth-inhibitory activity of JA and JA-Ile is probably a characteristic physiological response that is conserved in vascular plants. During plant evolution, JA biosynthesis and signaling in vascular plants might have developed due to the need to adapt to life on land and/or to organize more complicated physiological processes.
It has recently been reported that the genes encoding proteins related to JA signaling perception, such as COI1 and JAZ, are present in bryophytes and lycophytes (Wang et al. 2015). The JA responsiveness of S. moellendorffii shown in this study indicates that JA might function through a signal transduction mechanism similar to that observed in seed plants: the COI1–JAZ system. Given that JA does not have any significant effect in bryophytes, vascular plants are suggested to have developed the COI1–JAZ system, which is suitable for JA signaling. Investigating the functional difference in the COI1–JAZ system between bryophytes and lycophytes is important for understanding the functional evolution of JA signal transduction in plants. Selaginella moellendorffii will be a key plant for understanding the function of JA in vascular plants.
In conclusion, we provide the first evidence of the presence of OPDA, JA and JA-Ile in the model lycophyte S. moellendorffii, which is taxonomically positioned between bryophytes and euphyllophytes. Moreover, active homologous enzymes involved in the octadecanoid pathway, including an AOS, an AOC, an OPR3-like enzyme and a JAR1, were also identified in this plant. The application of exogenous OPDA, JA and JA-Ile significantly inhibited the growth of S. moellendorffii. Thus, OPDA, JA and JA-Ile are probably important signaling molecules in stress responses and growth regulators in S. moellendorffii.
Materials and Methods
Plant material
Selaginella moellendorffii was grown at 25°C under white fluorescent light with a photoperiod of 14 h light/10 h dark.
Analysis of OPDA, JA and JA-Ile
The microphylls and stems of S. moellendorffii (∼500 mg) were harvested 10, 20, 30, 60 or 180 min after wounding with tweezers and were homogenized in liquid nitrogen, followed by ethanol extraction for 24 h. UPLC-MS/MS analysis was performed according to Sato et al. (2009) and Yamamoto et al. (2015).
JA growth inhibitory activity assay in S. moellendorffii
Selaginella moellendorffii bulbils (Supplementary Fig. S17) were sprinkled onto a Jiffy-7 pellet (Sakata Seed Corporation). The bulbils were grown at 25°C under white fluorescent light with a photoperiod of 14 h light/10 h dark. After germination and 1 week of growth, the plants were treated with OPDA, JA or JA-Ile at concentrations of 25, 50 and 100 µM every other day for 4 weeks. The shoot growth of the plants was evaluated by morphometric analysis.
Bioinformatics analysis
All of the putative genes in S. moellendorffii were screened in the Phytozome v 11.0 database (https://phytozome.jgi.doe.gov/pz/portal.html) using the Arabidopsis amino acid sequence as a query (Supplementary Table S1). Multiple sequence analysis was performed using Clustal Omega and viewed in the GeneDoc software program. A phylogenetic tree was generated using Mega version 5.2 (Neighbor–Joining method). Bootstrap analysis was performed by resampling the data sets 1,000 times. The subcellular localization of SmAOC1 was predicted using ChloroP v 1.1 (http://www.cbs.dtu.dk/services/ChloroP/) and iPSORT (http://ipsort.hgc.jp/). The theoretical pI and molecular weight were analyzed using the Compute pI/Mw tool (http://web.expasy.org/compute_pi/).
Molecular cloning of SmAOS2, SmAOC1, SmOPR1, SmOPR5 and SmJAR1
The SmAOS2, SmAOC1, SmOPR1, SmOPR5 and SmJAR1 genes were cloned according to a conventional method. The details are described in the Supplementary Information.
Synthesis of recombinant SmAOS2, SmAOC1, SmOPR1 and SmOPR5
Recombinant SmAOS2, SmAOC1, SmOPR1 and SmOPR5 were produced according to a conventional method. The details are described in the Supplementary Information.
Enzyme assay for AOS
To measure AOS activity, SmAOS2 was reacted with 13(S)-HPOT according to Bandara et al. (2009) with some modifications. A 5 µg aliquot of recombinant SmAOS2 was reacted with 6 µg of 13(S)-HPOT as a substrate in a total volume of 0.1 ml of 50 mM phosphate buffer (pH 7.0). The reaction mixture was incubated at 25°C for 1 h. To terminate the enzymatic reaction, the mixture was treated with 1 M HCl until it reached pH 3. To convert cis-OPDA to trans-OPDA, the reaction mixture was subjected to alkaline treatment by stirring with 1 ml of 1 M NaOH for 30 min. After acidification, the resulting solution was extracted three times with an identical volume of ethyl acetate, evaporated, and methylated using ethereal diazomethane. After incubation for 30 min on ice, the reaction was terminated with the addition of one or two drops of acetic acid until the color changed from yellow to colorless. The methylated residue was dissolved in 100 µl of CHCl3 for chiral GC-MS analysis (1,200 liter GC-MS/MS system, Varian) equipped with a β-DEX-fused silica capillary column (0.25 mm × 30 m × 0.25 µm, Supelco). Helium was used as the carrier gas at a constant flow rate of 1.0 ml min–1. The ion source and vaporizing chamber were heated to 200 and 220°C, respectively, and the ionization voltage was 70 eV. The column was initially heated to 50°C and held isothermally for 1 min, and then the temperature was subsequently increased at a rate of 10°C min–1 to a final temperature of 190°C, which was held for 80 min.
Enzyme assay for AOC
The enzymatic reaction was performed in a mixture of 5 µg of 13-HPOT and 5 µg of PpAOS1 in 1 ml of 10 mM sodium phosphate buffer (pH 7.8) with or without 10 µg of SmAOC1. The reaction mixture was incubated at 25°C for 1 h. The products of the SmAOC1 reaction were analyzed according to the method of OPDA analysis described in the previous section.
Enzyme assays for OPR
The enzymatic reaction was performed according to the method of Schaller et al. (1998) with some modifications. A 5 µg aliquot of affinity-purified recombinant SmOPR1 or SmOPR5 was incubated in 50 mM potassium phosphate buffer (pH 7.4) containing 7.5 µg of (±)-cis-OPDA as a substrate and 1.0 mM NADPH in a total volume of 0.1 ml for 1 h at 25°C. At the end of the incubation, (±)-cis-OPDA and (±)-cis-OPC-8:0 were converted to methylated trans-isomers according to the method described above. OPC-8:0 was evaluated by chiral GC-MS analysis (same conditions as the enzymatic analysis for SmAOS2).
Subcellular localization of SmAOC1
Subcellular localization of SmAOC1 was analyzed according to the method of Yamamoto et al. (2015). The details are described in the Supplementary Information.
Enzyme assay for JAR1
SmJAR1 activity was assayed according to Wakuta et al. (2011) with some modifications. Escherichia coli BL21 (DE3) cells carrying pET23a-SmJAR1 were cultured overnight in 5 ml of LB medium supplemented with 100 µg ml–1 ampicillin, at 37°C. A portion of the culture was inoculated into 2 ml of fresh selective LB medium, incubated until the OD600 reached 0.2, and treated with 1 mM isopropyl-β-d-1-thiogalactopyranoside (IPTG). The culture was incubated at 37°C for 5 h to induce protein production. Subsequently, 0.2 mM JA was added to the culture. The incubation was continued at 37°C for 3 h with shaking. JA-Ile in the culture was analyzed using the UPLC-MS/MS system as previously described (Sato et al. 2009).
qRT-PCR for SmAOC1, SmOPR5 and SmJAR in S. moellendorffii exposed to wounding stress and JA treatments
The coding sequences (CDSs) for SmAOC1, SmOPR5 and SmJAR1 were used to design specific primers for qRT-PCR. The primers used to amplify these genes are listed in Supplementary Table S2. Total RNA of S. moellendorffii was extracted 0, 10, 20, 30, 60 and 180 min after wounding or 50 µM JA treatment. The first-strand cDNA was used as a template. In the present study, KOD SYBR qPCR Mix (Toyobo) was used according to the manufacturer’s protocol. Each reaction mixture contained 12.5 µl of KOD SYBR® qPCR Mix, 1 µl of each primer (10 mM), 1 µl of cDNA and 9.5 µl of MilliQ water. qRT-PCR was performed on a Thermal Cycler Dice Real Time system (TAKARA TP800). The qRT-PCR conditions were as follows: pre-incubation at 95°C for 30 s, followed by 40 cycles of 95°C for 5 s and 60°C for 30 s. The specificity of each PCR amplicon was assessed with a dissociation curve (95°C for 15 s, 60°C for 30 s and 95°C at 15 s). Each target gene was calculated and expressed as fold regulation in comparison with the housekeeping genes actin (XM_002976966) for the wounding treatment and ubiquitin EST (DN839032, Kirkbride et al. 2013) for the JA treatment.
Supplementary data
Supplementary data are available at PCP online.
Funding
This research was supported by grants from Hokkaido University.
Abbreviation
- AOC
allene oxide cyclase
- AOS
allene oxide synthase
- COI1
coronatine insensitive 1
- 12
13-EOT, (9Z,15Z)-(13S)-12,13-epoxyoctadeca-9,11,15-trienoic acid
- GC-MS
gas chromatography–mass spectrometry
- GFP
green fluorecent protein
- 13-HPOT
13-hydroperoxyoctadecatrienoic acid
- IPTG
isopropyl-β-d-1-thiogalactopyranoside
- JA
jasmonic acid
- JA-Ile
jasmonoyl-isoleucine
- JAR1
jasmonic acid resistant 1
- JAZ
jasmonate-ZIM domain
- LOX
lipoxygenase
- Ni-NTA
nickel–nitrilotriacetic acid
- OPC-8:0
3-oxo-2-(cis-2′-pentenyl)cyclopentane-1-octanoic acid
- OPDA
12-oxo-phytodienoic acid
- OPR
12-oxo-phytodienoic acid reductase
- qRT-PCR
quantitative real-time PCR
- SCF
Skp–Cullin–F box
- UPLC-MS/MS
ultraperformance liquid chromatography–tandem mass spectrometry
Acknowledgments
The authors are grateful to Dr. Hasabe at the National Institute for Basic Biology for kindly providing S. moellendorffii and the plasmid, pUGWnew5. Methyl jasmonate, which was needed to synthesize deuterium-labeled compounds, was kindly provided by Zeon corporation (Tokyo, Japan).
Disclosures
The authors have no conflicts of interest to declare.