Methyl Jasmonate and 1-Methylcyclopropene Treatment Effects on Quinone Reductase Inducing Activity and Post-Harvest Quality of Broccoli

Effect of pre-harvest methyl jasmonate (MeJA) and post-harvest 1-methylcyclopropene (1-MCP) treatments on broccoli floret glucosinolate (GS) concentrations and quinone reductase (QR, an in vitro anti-cancer biomarker) inducing activity were evaluated two days prior to harvest, at harvest and at 10, 20, and 30 days of post-harvest storage at 4 °C. MeJA treatments four days prior to harvest of broccoli heads was observed to significantly increase floret ethylene biosynthesis resulting in chlorophyll catabolism during post-harvest storage and reduced product quality. Post-harvest treatment with 1-methylcyclopropene (1-MCP), which competitively binds to protein ethylene receptors, maintained post-harvest floret chlorophyll concentrations and product visual quality in both control and MeJA-treated broccoli. Transcript abundance of BoPPH, a gene which is responsible for the synthesis of pheophytinase, the primary enzyme associated with chlorophyll catabolism in broccoli, was reduced by 1-MCP treatment and showed a significant, negative correlation with floret chlorophyll concentrations. The GS, glucobrassicin, neoglucobrassicin, and gluconasturtiin were significantly increased by MeJA treatments. The products of some of the GS from endogenous myrosinase hydrolysis [sulforaphane (SF), neoascorbigen (NeoASG), N-methoxyindole-3-carbinol (NI3C), and phenethyl isothiocyanate (PEITC)] were also quantified and found to be significantly correlated with QR. Sulforaphane, the isothiocyanate hydrolysis product of the GS glucoraphanin, was found to be the most potent QR induction agent. Increased sulforaphane formation from the hydrolysis of glucoraphanin was associated with up-regulated gene expression of myrosinase (BoMyo) and the myrosinase enzyme co-factor gene, epithiospecifier modifier1 (BoESM1). This study demonstrates the combined treatment of MeJA and 1-MCP increased QR activity without post-harvest quality loss.


Introduction
Brassica vegetables are recognized as functional foods that have putative cancer preventive effects as shown in epidemiological and animal carcinogenesis studies [1]. The glucosinolates (GS) including glucoraphanin, gluconasturtiin, and sinigrin found in the tissues of accessions of Brassica oleracea have been identified as potent cancer prevention agents because products of their hydrolysis by the endogenous enzyme myrosinase generate sulforaphane, phenethyl isothiocyanate (PEITC), and allyl isothiocyanate (AITC). These isothiocyanate products have been shown to induce phase II detoxification enzymes such as glutathione S-transferases (GSTs) and quinone reductase (QR) in in vitro or in vivo systems that can enhance detoxification and elimination of carcinogens from the human body [2][3][4]. QR activity elevation with in vitro and in vivo systems has been shown to correlate with induction of other protective phase II enzymes such as the GSTs and provides a reasonable biomarker for the potential chemoprotective effect of phytochemicals against initiation of carcinogenesis [5]. In addition, several hydrolysis products from other GS found in B. oleracea cultivars including indole-3carbinol [6] and N-methoxyindole-3-carbinol [7,8] have also been reported to promote antiproliferative bioactivity in human cancer cell lines.
GS biosynthesis genes have been intensively studied in Arabidopsis using biochemical assays. There is high homology of the genes in GS biosynthesis between Arabidopsis and Brassicaceae [9]. CYP79 catalyzes the conversion of amino acids to aldoximes. CYP79F1 and CYP79F2 genes are responsible for aldoxime metabolism leading to aliphatic GS derived from chain-elongated methionine derivatives, whereas CYP79B2 and CYP79B3 have distinct functions in indolyl GS biosynthesis, which is derived from tryptophan [10]. In the biosynthetic pathway of indolyl GS, CYP79B2 gene product catalyzes the conversion of tryptophan to indole-3acetaldoxime, with CYP83A1 and CYP83B1 metabolizing the phenylalanine-and tyrosine-derived aldoximes [10]. It has been reported that indolyl GS biosynthesis is increased by jasmonic (JA) treatment [11]. The biological activity of GS varies with diversity of structure of the side chain that is the last step of GS biosynthesis [12]. Recently, methoxylation genes involved in glucobrassicin biosynthesis such as CYP81F2, CYP81F3, and CYP81F4 were identified by genetic engineering Arabidopsis indolyl GS biosynthesis into Nicotiana benthamiana [13]. CYP81F2 gene product is responsible for the methoxylation of glucobrassicin resulting in 4-methoxyglucobrassicin. The hydrolysis product of 4-methoxyglucobrassicin has been reported to be antibiotic to fungal pathogens and to the green peach aphid (Myzus persicae) [14,15].
Intact GS do not display bioactivity but following hydrolysis by the endogenous enzyme myrosinase, isothiocynates and other products are generated, which have been associated with insect resistance and anti-cancer activity. When the plant tissue is disrupted, myrosinase and substrates (GS) come into contact, resulting in GS hydrolysis. The chemical structure of hydrolysis products depends on the structure of the GS side chain and reaction conditions such as pH, concentration of Fe 2+ and presence of epithiospecifier protein (ESP), a myrosinase co-factor that will favor formation of nitriles [16]. In the absence of ESP, the addition of Fe 2+ ions also promotes nitrile formation, which are essentially without anti-cancer activity compared to the isothiocyanates like sulforaphane, PEITC, and AITC [17]. The epithiospecifier modifier 1 (ESM1) gene in Arabidopsis encodes a protein shown to inhibit function of ESP, leading to increased isothiocyanate production from GS hydrolysis [18].
Methyl jasmonate (MeJA), a plant signal transduction compound associated with herbivore defense, can act as an elicitor to enhance GS biosynthesis [19]. Previous research has shown that MeJA treatments can significantly increase QR inducing activity mediated by enhancement of GS biosynthesis including glucoraphanin, glucobrassicin and neoglucobrassicin in cauliflower [20]. However, 1 mM MeJA treatment was also found to significantly promote ethylene production and increased 1-aminocyclopropane-1-carboxylate acid (ACC) concentrations and ACC oxidase activity associated with senescence and loss of product quality in broccoli [21]. Thus, while MeJA treated broccoli can display enhanced QR activity associated with increased GS concentrations, elevated ethylene production can accelerate post-harvest senescence, phytochemical degradation and visual quality loss.
Inhibition of plant ethylene perception or blocking the ethylene receptor using 1-methylcyclopropene (1-MCP) is an effective way to improve shelf life and quality of fruits and vegetables [22] and application of 1-MCP increases shelf life of broccoli [23]. 1-MCP application has also been found to maintain the phytochemicals in broccoli such as chlorophylls, carotenoids, ascorbic acid and GS after harvest by binding to the ethylene receptors ETR1 and ETR2 [24][25][26]. Visual color is a critical factor in retailer and consumer evaluation of product quality and subsequent purchasing decisions [27]. Chlorophyll content is considered a good indicator of broccoli post-harvest quality. Previous studies reported that pheophytinase (PPH) and pheophorbide a oxygenase (PaO) are key enzymes in post-harvest chlorophyll breakdown [28,29]. It was reported that gene expression of BoPPH and BoPaO is reduced by 1-MCP treatment [30].
In this experiment, we evaluate MeJA treatment on postharvest quality and phytochemical content of broccoli to see if 1-MCP can modulate MeJA initiated postharvest senescence. Thus, the objectives of this research are to evaluate the effect of pre-harvest MeJA and post-harvest 1-MCP treatment on postharvest physicochemistry and quinone reductase bioactivity of broccoli floret extracts.

Plant Cultivation and Sample Preparation with Treatments
'Green Magic' broccoli (Sakata Seed Co., Morgan Hill, CA) was used for this experiment. Broccoli seeds were germinated in 32 cell plant plug trays filled with sunshine® LC1 professional soil mix (Sun Gro Horticulture, Vancouver, British Columbia, Canada). Seedlings were grown in a greenhouse at the University of Illinois at Champaign-Urbana under a 25°C /15 °C and 14 h/10 h: day/night temperature regime and with supplemental lighting. Forty days after seed germination, seedlings were first transferred into 1-liter pots and then after a month 150 broccoli seedlings were repotted into 3.75-L pots. These broccoli seedlings were evenly placed on three greenhouse benches and control and MeJA treatment assigned within each bench to minimize micro-environmental variation. 500 micromoles of MeJA (Sigma-Aldrich, St. Louis, MO, USA) in solution containing 0.1% ethanol was sprayed on aerial tissues of each of the treated plants four days prior to harvest at commercial maturity. Timing of MeJA sprays and concentration of solution was previously determined to optimize up-regulation of indolyl GS [31]. For the control group, only a 0.1% ethanol solution was applied. At commercial market maturity 50 broccoli heads were harvested from both the control and MeJA treated plants, transported to the laboratory, and divided into branchlets of broccoli florets. Branchlets of control and MeJA treated plants were each randomly divided Enhancing Broccoli Quinone Reductase Activity PLOS ONE | www.plosone.org into two groups generating four treatment groups: (1) No MeJA or 1-MCP (Control); (2) No MeJA and 500 ppb treatment with 1-MCP for 24 h (1-MCP); (3) MeJA without 1-MCP (MeJA), and (4) MeJA and 500 ppb treatment with 1-MCP for 24 h (MeJA_1-MCP). Treatments (1,3) and treatments (2,4) were placed in airtight plastic containers at 20 °C. 1-MCP was generated in containers holding treatments (2,4) by adding an activator and a Smartfresh ® tablet (AgroFresh, Inc. a division of Rohm and Hass, Philadelphia, PA, USA) to the activation solution following the instructions provided by the company. After treatment, broccoli branchlets were stored in a walk-in cooler at 4 °C. At each sampling date [-2 (2 days after MeJA treatment), at harvest (4 days after MeJA treatment, day 0 for 1-MCP), and at 10, 20, and 30 days of post-harvest storage], three random subsamples of branchlets (replications) of each treatment group were selected and assayed for ethylene and CO 2 production and visual quality. Pictures of broccoli florets and their relative visual quality from each assay date are presented in Figure 1. After measuring CO 2 , ethylene production, and hue angle, a measurement of floret color change, a subsample of tissue from each replication was taken, frozen in liquid nitrogen, and stored at -80 °C until ground with a mortar and pestle in liquid nitrogen for RNA extraction. Residual sample tissues from each replication were freeze-dried. Freeze-dried broccoli floret tissue of each sample was finely ground with a commercial coffee grinder. The ground freeze-dried broccoli samples were stored at -20 °C prior to GS quantification and quinone reductase bioactivity assay.

Respiration and Ethylene Production Measurement
Respiration was measured as tissue CO 2 production. Three subsamples (300 g each) of broccoli branchlets from each treatment were placed into 3 L jars and enclosed with a silicon rubber cap for 1 h at 20 °C. Sample CO 2 was estimated using 2% CO 2 in a nitrogen gas (v/v) standard for each experiment. The headspace gas in the jar was sampled with a 0.2 mL plastic hypodermic syringe and injected into a GC (model Perkin Elmer AutoSystem Gas Chromatograph) equipped with Enhancing Broccoli Quinone Reductase Activity PLOS ONE | www.plosone.org a Propak® (Waters Co., Milford, MA) column and thermal conductivity detector (TCD). Temperature of the injector, detector and column was 100, 150 and 30 °C, respectively. The results were expressed as mL of CO 2 /kg/h. Ethylene measurement was measured as previously reported [20] using pure ethylene gas as a standard for estimating sample concentrations.

Determination of Total Chlorophyll Content
Frozen floret tissue samples from each sampling date (75 mg) were ground and extracted in 1.5 ml of 80% acetone in a 2 mL tube with vigorous vortexing for 1 h. Total chlorophyll content was determined by using the equations listed below [32].

Hue Angle Measurement
As a reliable measure of color change during post-harvest storage sample, hue angle was measured as by using a LabScan XE colorimeter (Hunter Associates Laboratory, Reston, VA, USA) generating values for a* (redness and greenness), and b* (yellowness and blueness). The instrument was calibrated with a standard white and black tile. The average of four different broccoli branchlets was recorded in each replication. Hue degree (h°) was calculated as h° = tan -1 (b*/a*) when a*>0 and b*>0, or as h° = 180° -tan -1 (b*/a*) when a<0 and b>0 [20].

Determination of Sample GS Concentrations
Freeze-dried broccoli powder (0.2 g) and 2 mL of 70% methanol were added to 10 mL tubes (Nalgene, Rochester, NY, USA) and heated on a heating block at 95 °C for 10 min. After cooling on ice, 0.5 mL benzylglucosinolate (1 mM) was added as an internal standard (POS Pilot Plant Corp, Saskatoon, SK, Canada), mixed, and centrifuged at 3,000 × g for 15 min at 4 °C. The supernatant was saved and the pellet was re-extracted with 2 mL 70% methanol at 95 °C for 10 min and the two extracts combined. A subsample (1 mL) from each pooled extract was transferred into a 2-mL microcentrifuge tube. Protein was precipitated with 0.15 mL of a 1:1 mixture of 1 M lead acetate and 1 M barium acetate. After centrifuging at 12,000 × g for 1 min, each sample was then loaded onto a column containing DEAE Sephadex A-25 resin (Sigma-Aldrich) for desulfation with arylsulfatase (Helix pomatia Type-1, Sigma-Aldrich) for 18 h and the desulfo-GS eluted. One hundred µL of each sample were injected on to a HPLC. Quantification of GS using high-performance liquid chromatography was performed using a previously described protocol [33].

Measurement of Myrosinase Activity
Myrosinase activity was optimized according to previous studies [34][35][36]. Crude extracts were prepared by adding 0.3 g of a finely ground freeze-dried sample in 4 mL of an extraction buffer consisting of 10 mM potassium phosphate, 1 mM ethylenediaminetetraacetic acid (EDTA), 3 mM dithiothreitol (DTT) and 5% glycerol (pH 7.0) for 20 min in an ice bath [35]. The crude extracts were centrifuged at 15,000 × g for 30 min at 4 °C. To remove endogenous GS and glucose, the crude extract was filtered through an Amicon ultrafiltration cell (Millipore, Billerica, MA, USA) with a 10 KDa molecular weight cutoff [34] and washed several times at 4 °C using the same extraction buffer at pH 7.0 [35]. Fifty µL of purified extracts and 450 µL of 0.2 mM sinigrin in 33.3 mM phosphate buffer, pH 6.5 were mixed and incubated for 40 min [36]. To stop the enzyme reaction extracts were heated at 95 °C for 10 min. The release of glucose was determined by the glucose oxidase/peroxidase/ ABTS method [37] using a microplate reader (Biotek Instruments, Winooski, VT, USA). Glucose concentrations were calculated using a linear standard curve. By calculating the glucose amount in aliquots of purified extracts without sinigrin, endogenous glucose levels were subtracted in purified extracts for myrosinase activity measurement.

Quinone Reductase (QR) Inducing Activity
For the QR assay, 75 mg of broccoli floret powder from each sample was suspended in 1.5 mL of water in the absence of light for 4 h at room temperature in a sealed 2 mL microcentrifuge tube (Fisher Scientific, Waltham, MA, USA) to facilitate GS hydrolysis by endogenous myrosinase. Slurries were then centrifuged at 12,000 × g for 10 min and supernatants were diluted to 0.25% final concentration in the QR activity assays. The QR inducing activities of different samples were determined by means of the protocol described by Prochaska and Santamaria [38].

Analysis of Glucosinolate Hydrolysis Products
The extraction and analysis of isothiocyanates and other hydrolysis products was carried out according to previously published methods with some modifications [39]. Broccoli powder (75 mg) was suspended in 1.5 mL of water in the absence of light for 4 h (optimal time for hydrolysis products of indolyl GS) and 24 h (optimal time for sulforaphane and PEITC generation from glucoraphanin and gluconasturtiin) at room temperature in a sealed 2 mL microcentrifuge tube (Fisher Scientific) to facilitate GS hydrolysis by endogenous myrosinase. Slurries were then centrifuged at 12,000 × g for 5 min and supernatants was decanted into a 2 mL microcentrifuge tube. Twenty μL of butyl isothiocyanate (0.5 mg/mL) and 4-methoxyindole (1 mg/mL, Synthonix, Wake Forest, NC, USA) were added as the internal standards for isothiocyanates and hydrolysis products of indolyl GS to quantify indole-3-carbinol (I3C), NI3C, and neoascorbigen (NeoASG), respectively, with 0.5 mL of methylene chloride. NI3C and NeoASG were purified and identified by LC-MS [8,40]. Tubes were shaken vigorously before being centrifuged for 2 min at 9,600 × g. The methylene chloride layer (200 µL) was transferred into a 350 µL flat bottom insert (Fisher Scientific) in a 2 mL HPLC autosampler vial (Agilent, Santa Clara, CA, USA) for mixing with 100 µL of a reagent containing 20 mM triethylamine and 200 mM ß-mercaptoethanol (derivatization reagent) in methylene chloride. For SF and PEITC, unlike other hydrolysis products of GS measurement, 0.5 mL of fresh broccoli extracts were mixed with 0.5 mL of Enhancing Broccoli Quinone Reductase Activity PLOS ONE | www.plosone.org derivatization reagent using an orbital shaker at 220 rpm for 24 hours. Then, internal standards were added as described above. The mixture was incubated at 30 °C for 60 min under constant stirring, and then dried under a stream of nitrogen. The residue containing isothiocyanate derivatives (isothiocyanate-mercaptoethanol derivatives) and other hydrolysis compounds were dissolved in 200 µL of acetonitrile/ water (1:1) (v/v) and 10 µL of this solution injected onto a Agilent 1100 HPLC system (Agilent). Extracts were separated on a Eclipse XDB-C18 column (150 × 4 mm, particle size 5 μm, Agilent) with a Adsorbosphere C18 all-guard™ cartridge precolumn (Grace, Deerfield, IL). Mobile phase A was water and B methanol. Mobile phase B was 0% at injection, increasing to 10% by 10 min, 100% at 35 min, and held 5 min, then decreased to 0% by 50 min. Flow rates were kept at 0.8 mL/ min. The detector wavelength was set at 227 and 271 nm. Response factors for monomeric indolyl derivatives were used from a previous report [40]. Due to a lack of standards for NI3C and NeoASG the standard curve of I3C was applied for quantification of both NI3C and NeoASG. The quantities were expressed as I3C molar equivalent concentrations.

Cloning of Broccoli Epithiospecifier Modifier 1 (BoESM1)
Using known Arabidopsis (NM_112278.2), Brassica rapa (FJ830451.1) and Brassica napus (FJ830448.1) gene sequence information, PCR primers were designed with the Primer3 software package (http://frodo.wi.mit.edu/primer3) to isolate the broccoli, cabbage, and cauliflower homologous ESM1 gene, known to be associated with GS hydrolysis. PCR amplification was performed using the GoTaq® PCR Core System (Promega, Madison, WI, USA) following the protocol described by the manufacturer. The PCR product was separated on 1% TAE gel and purified by using a Qiagen gel extraction kit (QIAGEN, Valencia, CA, USA) according to the manufacturer's protocols. The amplified PCR products were cloned with pGEM®-T Easy Vector System (Promega), and the clones were sequenced in the W. Carver Biotechnology Center, University of Illinois at Urbana-Champaign. The amino acid sequences deduced from the isolated cDNA sequences were subjected to phylogenetic tree analysis using Clustal W2 (http://www.ebi.ac.uk/Tools/clustalw2/).

RNA Extraction and Quantitative Real Time-PCR
Total RNA was isolated from control and MeJA treated floret tissue samples using the RNeasy Mini Kit (QIAGEN) according to the manufacturer's instructions. The quantity of RNA was measured using a NanoDrop 3300 spectrophotometer (Thermo Scientific, Waltham, MA). One μg of the total RNA was reversetranscribed with Superscript™ III First-Strand Synthesis SuperMix for qRT-PCR (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions. The resulting cDNA samples were diluted to 1/10 their concentrations (v/v) for qRT-PCR. Previously reported primer sets of GS biosynthesis (BoCYP79B2, BoCYP83B1, BoCYP79F1, and BoCYP83A1) genes, hydrolysis (BoMyo, BoESP, and BoESM1) genes, a pathogenesis-related (PR) protein (BoPR) gene known to be responsive to MeJA, chlorophyll catabolism (BoPPH and BoPaO) genes, and the broccoli actin gene (BoACT1) as a normalization standard were used for qRT-PCR [41,42] (Table 1). Transcript abundance of the broccoli actin (BoACT1) gene from each treatment was stable throughout

Statistical Analysis
Statistical analyses were conducted using the JMP 10 software (SAS institute Inc., Cary, NC). Student's T-tests were used for comparing treatment groups. Fisher's Least Significant Difference (LSD) test was conducted for comparing treatment group means at P ≤ 0.05. Pearson correlation was conducted on all pairs of a GS, hydrolysis product and chlorophyll concentrations, gene expression, and QR inducing activity based on the mean values of each treatment across postharvest storage dates. The results are presented as means ± SD.

Ethylene Production and Respiration Rate for Broccoli Florets Subjected to MeJA Treatments
Treatment with 500 µM MeJA significantly increased ethylene production (1.9 fold) in broccoli floret tissues four days after treatment ( Figure S3A). Ethylene production dropped significantly during postharvest storage at 4 °C. Ethylene production between control and MeJA treated groups were not significantly different at 10 and 20 days of storage regardless of 1-MCP treatment. There was also no consistent difference in respiration rates among the different treatments at harvest or during postharvest storage ( Figure S3B).

Product Color Measurement for Visual Quality Change and Chlorophyll Concentrations
Images of broccoli florets from the four different treatments from each assay date are presented in Figure 1. There were differences in visual quality between control and 1-MCP treatment groups over the period of post-harvest storage, regardless of MeJA treatment (Figure 1). In order to determine this objectively, tissue chlorophyll concentrations and floret hue angle were measured to quantify visual quality. There was a significant reduction in total chlorophyll content in broccoli florets two (day-2) and four days (head harvest) after MeJA treatment compared to controls (Figure 2A). It was previously reported that MeJA treatment reduced total chlorophyll content in Arabidopsis thaliana [43]. Chlorophyll b concentrations are much higher than chlorophyll a in broccoli and showed more dramatic losses during post-harvest storage (data not shown). Hue angle measurements indicated that 1-MCP treatments were associated with superior visual quality throughout the period of post-harvest storage compared to controls and the MeJA treatment ( Figure 2B).

Chlorophyll Catabolism Gene Expression by Preharvest MeJA and Post-harvest 1-MCP Treatments
To investigate the mechanism of chlorophyll and visual quality loss during the post-harvest period, relative transcript abundance of two genes associated with chlorophyll catabolism BoPPH and BoPaO were assayed by qRT-PCR. Transcript abundance of these genes was significantly greater in the control or MeJA treated broccoli for BoPPH at 20 and 30 days of post-harvest storage and for BoPaO at 30 days postharvest storage ( Figure 2C, D). There was a significant negative correlation (r = -0.642, P = 0.007) between BoPPH gene expression ( Figure 2C) and total chlorophyll concentrations ( Figure 2A, Table 2). Hue angle measurements of visual quality were negatively correlated with both BoPPH (r = -0.868, P < 0.001) and BoPaO gene expression (r = -0.641, P = 0.014, Table 2). Down-regulation of expression of these genes to maintain visual quality has been previously reported in broccoli [28]. 1-MCP mediated reduction of ethylene binding to receptor proteins has previously been shown to be responsible for reduced expression of BoPPH and BoPaO, chlorophyll degradation and associated visual quality loss during post-harvest storage [30].

Pre-harvest MeJA and Post-harvest 1-MCP Treatments Influence GS and GS Hydrolysis Product Concentrations
Treatment with MeJA significantly increased glucobrassicin, 4-methoxyglucobrassicin, neoglucobrassicin, gluconasturtiin, and total glucosinolate concentrations in broccoli floret samples. At harvest (day 0), the relative increase of these GS was 1.58, 4.75, 4.71, 2.28, and 1.49 fold over controls, respectively ( Figure 3). In the case of 4-methoxyglucobrassicin, there was more dramatic increase in MeJA treated broccoli compared to control during post-harvest storage. Previous reports indicate MeJA treatments increased glucoraphanin, glucobrassicin, and neoglucobrassicin concentrations in cauliflower curds [20]. Both jasmonate (JA) and MeJA induced significant increases (up to 20-fold) in the concentration of specific indolyl GS in Brassica napus (primarily the GS, glucobrassicin), in B. rapa (primarily 4-hydroxy glucobrassicin), and in B. juncea (both increased) [44]. The different responses to MeJA treatment suggest that variation in GS response to MeJA is species specific. Postharvest 1-MCP treatment maintained glucoraphanin ( Figure 3A) and glucobrassicin ( Figure 3B) concentrations regardless of MeJA treatment. The combination of MeJA and 1-MCP treatment showed the highest concentrations for glucobrassicin, neoglucobrassicin, gluconasturtiin, and total glucosinolate during post-harvest storage (Figure 3). MeJA treated broccoli showed significant reduction in glucoraphanin, glucobrassicin, and neoglucobrassicin concentrations during the first 10 days of post-harvest storage ( Figure 3A, B, and D). This implies that wounding or damage to broccoli heads during post-harvest storage should be minimized to maintain concentrations of glucosinolates associated with health promotion.
The glucobrassicin concentrations of all four treatment groups observed rapidly decreased ( Figure 3B) Figure 3D), which has a structure similar to glucobrassicin (methoxylated glucobrassicin). The concentration of 4-methoxyglucobrassicin was increased in all four treatment groups during post-harvest storage with the MeJA treatment showing the largest increment ( Figure 3C). Neoglucobrassicin and 4-methoxyglucobrassicin are products of the glucobrassicin biosynthesis pathway following hydroxylation then methylation, respectively [10]. Considering this biosynthetic pathway, reduction of glucobrassicin may be associated with 4-methoxylation of glucobrassicin. We found that there was a significant correlation between the loss of glucobrassicin and increases of 4-methoxyglucobrassicin (r = -0.578, P = 0.019, Table 2), implying active GS conversion [13] during post-harvest storage at 4° C.
MeJA treatments were also found to significantly increase gluconasturtiin concentrations at harvest ( Figure 3E).
Concentrations of gluconasturtiin were observed to increase during the period of post-harvest storage in all the treatments except for the MeJA treated group. Since indolyl GS (glucobrassicin, 4-methoxyglucobrassicin, and neoglucobrassicin) and aromatic GS (gluconasturtiin) share the same biosynthetic pathway, increases in specific GS concentrations may reduce other GS levels with shared precursors. The increment of 4-methoxyglucobrassicin in MeJA treated broccoli florets may interfere with neoglucobrassicin and gluconasturtiin biosynthesis. Post-harvest 1-MCP treatment may contribute to the accumulation of gluconasturtiin by delaying senescence but the effect was greater in MeJA treated broccoli compared to non-MeJA treated broccoli.
Interestingly, SF formation was significantly increased by MeJA treatment even though there were no significant increases in glucoraphanin concentrations (Figure 3 and Figure 4). This can be explained by the significantly increased levels in gene expression of myrosinase (BoMyo) and BoESM1 compared to the BoESP gene ( Figure 5). Myrosinase activity was significantly (58%) increased by MeJA treatment at four days after treatment ( Figure S4). MeJA treatment not only increased gluconasturtiin concentration but also PEITC formation (Figure 3 and Figure S5). NeoASG concentrations were significantly increased by MeJA treatment and maintained elevated concentrations during post-harvest storage (Figure 4). The major hydrolysis product of neoglucobrassicin was NeoASG ( Figure 4C, D). I3C concentrations were observed to be significantly increased in only the MeJA treatments at harvest ( Figure 4B).

GS Biosynthetic, Hydrolytic, and Pathogenesis Related-Protein (PR) Gene Expression Changes Mediated by MeJA and/or 1-MCP treatments
The amino acid sequences deduced from the isolated B. oleracea ESM1 gene sequences corresponded to the B. napus (96%), B. rapa (95%) and Arabidopsis thaliana protein, ESM1 (79%), respectively, which suggests similar gene product function in broccoli by sharing homologous protein motifs ( Figure S1). In order to evaluate the effects of MeJA and 1-MCP treatments at harvest and during post-harvest storage, gene expression of GS biosynthetic (BoCYP79F1, BoCYP79B2, BoCYP83A1, and BoCYP83B1), hydrolytic (BoMyo, BoESP, and BoESM1), and pathogenesis related (BoPR) genes were measured by qRT-PCR. While gene expression of BoCYP79F1 (0.7 fold) was significantly decreased by MeJA treatment compared to control at two days after treatment, transcript abundance of BoCYP79B2 (10.0 fold), BoCYP83A1 (1.3 fold, non-significant increase), BoCYP83B1 (2.4 fold), BoMyo (2.9 fold), BoESM1 (1.9 fold), BoESP (1.7 fold), and BoPR (1.5 fold) were non-significantly or significantly increased compared to controls at two days after treatment ( Figure 5). The elevated transcript abundance observed at 2 days after MeJA treatment was dramatically reduced with post-harvest storage at 4 °C ( Figure 5). While gene expression levels during post-harvest storage were low and not dramatically different among treatments for many of the genes, transcript abundance for BoCYP79B2, BoCYP83B1, and BoPR gradually increased over the duration of storage. The mRNA expression levels of these genes may explain the increase in indolyl GS during post-harvest storage as was previously reported [45]. Decreased gene expression of BoCYP79B2 and BoCYP83B1 associated with 1-MCP treatment is in agreement with previous research in Arabidopsis showing that that these GS biosynthetic genes are stimulated by elevated MeJA-mediated ethylene production [11,46,47]. Table 2. Correlations among phytochemicals, QR inducing activity and gene expression of broccoli florets during postharvest storage at 4 °C. 1  2  3  4  5  6  7  8  9  10  11  12  13  14  15  16  17  18  19
The observed reductions in glucobrassicin with complimentary increases in 4-methoxyglucobrassicin concentrations are likely a result of the 4-methoxylation of glucobrassicin. The Arabidopsis CYP81F2 gene is involved in accumulation of 4-methoxyglucobrassicin synthesized from glucobrassicin in response to pathogen infection [13,14]. The loss of glucobrassicin in MeJA treated broccoli during the postharvest storage is possibly due to the higher levels of BoCYP81F subfamily transcripts compared to BoCYP79B2 and BoCYP83B1 gene expression. Mikkelsen et al. (2003) [11] reported that enzyme activities responsible for the Nmethoxylation of glucobrassicin are strongly induced by MeJA treatment and this induction is suppressed by ACC. Recently, it has been shown that the CYP81F4 gene product is involved in N-hydroxylation of glucobrassicin to synthesize neoglucobrassicin [13]. As Mikkelsen et al. (2003) [11] observed, it has been also reported that CYP81F4 was upregulated by MeJA and down-regulated by ethylene in Arabidopsis (5.81 fold) [48,49]. Consequently, ethylene accumulation in MeJA treated broccoli favors the pathway of methoxylation from glucobrassicin to 4-methoxyglucobrassicin rather than formation of neoglucobrassicin during post-harvest storage. Favoring the methoxylation pathway by accumulation of ethylene or the ethylene precursor, ACC would facilitate defense against postharvest pathogens since 4methoxyglucobrassicin has been shown to be antibiotic to fungi [14,50]. PEITC, the hydrolysis product of gluconasturtiin has also been reported to possess antifungal activity [51]. The increased levels of gluconasturtiin may be associated with upregulation of BoCYP83B1, which is involved in both indolyl and aromatic GS biosynthesis. Enhanced levels of PEITC during the post-harvest storage could also be associated with antifungal defense.

Correlations between QR inducing Activity of Broccoli and Enhanced GS and Hydrolysis Products
QR inducing activity of MeJA treated broccoli floret extracts was significantly increased compared to controls at 2 and 4 days after MeJA treatment and throughout the course of postharvest storage ( Figure 4F). In order to elucidate the major QR inducers in MeJA treated broccoli, correlation analysis was conducted between QR inducing activity and other variables (Table 2). Significant positive correlations were observed between QR inducing activity and the GS concentrations of gluconasturtiin (r = 0.712, P = 0.004), 4-methoxyglucobrassicin (r = 0.716, P = 0.002), and neoglucobrassicin (r = 0.804, P < 0.001) in broccoli floret sample extracts during post-harvest storage at 4 °C (Table 2). Previously, 4-methoxyindole-3carbinol has been reported to provide in vitro antiproliferation activity in two different human colon cancer cells lines [52].
This increased isothiocyanate formation induced by MeJA treatment also suggests that plant defense against herbivores may be involved in enhancing isothiocyanate formation as previous research reported that increased allyl isothiocyanate concentrations reduced survival and growth, and increased development time of Pieris rapae [53]. The detailed mechanism of how isothiocyanate formation is increased by MeJA treatment or insect invasion is not yet fully understood. However, a possible explanation could be related to the proteins and/or co-factors associated with hydrolysis. MeJA treatment increased gene expression of broccoli myrosinase as previous research reported [54]. The protein co-factor ESP when bound to myrosinase favors hydrolysis and conversion of GS products into nitriles [17] while elevated levels of EMS1 protein in Arabidopsis has been associated with enhanced hydrolysis of GS into isothiocyanate products [18]. Compared to broccoli myrosinase, transcript abundance of BoESP was not dramatically increased. Increased levels of unbound myrosinase, free of the co-factor, ESP favors the generation of isothiocyanates instead of QR inactive nitrile forms [17]. If cofactors such as ESP and ESM1 act competitively to bind with myrosinase, an increase of ESM1 gene expression should lead to enhanced isothiocyanate formation as was observed [55]. Compared to gene expression of myrosinase and its co-factors at two days after MeJA treatment, transcript abundance decreased at four days after treatment. This indicates how rapidly the regulatory response of herbivore defense mechanisms can act. Since isothiocyanates are toxic to both plants and herbivores, Brassica plants will also produce enhanced levels of detoxifying enzymes such as glutathione transferase for isothiocyanate neutralization [56]. To minimize metabolic costs under non-invasive conditions and reduce isothiocyanate autotoxicity the plant can lower glucosinolates and myrosinase biosynthesis to allocate resources for growth and development while maintaining the capacity to rapidly respond to biotic stress by selectively modifying expression of GS biosynthetic genes and BoMyo and its associated protein co-factors.
A standard format for comparing the efficacy of elicitors of QR induction is determining concentrations required for a twofold increase in activity (CD value). The CD value of NI3C, NeoASG, PEITC and SF are 35, 38.5, 5.0, and 0.2 µM, respectively [3]. Even though SF concentrations observed in our study are smaller than those of NeoASG (Figure 4), the relative SF bioactivity suggests it is the major contributor toward enhanced QR inducing activity in our extracts, although MeJA-mediated increases in other GS and their hydrolysis products are also likely contributors. More research is needed to determine if the hydrolysis products of 4methoxyglucobrassicin induce QR activity or not ( Figure 6).

Proposed Model of Pre-harvest MeJA and Post-harvest 1-MCP Treatment Effects on Postharvest Physiology and Health Promoting Bioactivity of Broccoli Florets
Our study suggests that pre-harvest MeJA and postharvest 1-MCP treatments have influence on GS de novo synthesis, conversion, and hydrolysis ( Figure 6). In most previous studies typically only MeJA mediated GS concentration changes were reported without information about the hydrolysis products, which are the active agents of QR induction. This present study suggests that MeJA treatment not only increases concentrations of certain GS (Figure 3), but can also increase Enhancing Broccoli Quinone Reductase Activity sulforaphane and PEITC conversion rates ( Figure S5) from precursor GS by modulating synthesis of myrosinase protein and its cofactors including BoESP and BoEMS1 ( Figure 5). This study confirmed that de novo GS biosynthesis in broccoli florets occurs during post-harvest storage. As Figure 5B and D show, BoCYP79B2 and BoCYP83B1 are responsible for the de novo indolyl and aromatic GS biosynthesis during post-harvest storage. GS biosynthesis genes are regulated by ethylenemediated signals. Thus, 1-MCP treatment reduced gene expression of GS biosynthesis genes including BoCYP79B2 and BoCYP83B1 by blocking ethylene and ethylene receptor protein binding during postharvest storage ( Figure 5B and D). Ethylene accumulation induces 4-methoxylation of glucobrassicin rather than N-methoxylation of glucobrassicin during post-harvest but 1-MCP treatment helped to maintain glucobrassicin and neoglucobrassicin concentrations by inhibiting 4-methoxylation. 1-MCP treatment maintained postharvest quality (Figure 1) by reducing chlorophyll catabolism gene expressions including BoPPH and BoPaO ( Figure 2). This study demonstrates the combined treatment of MeJA and 1-MCP can improve the delivery of anti-cancer compounds to consumers while maintaining postharvest visual quality ( Figure 6). Figure S1. Phylogenetic tree of epithiospecifier modifier 1 (ESM1) co-factor associated with glucosinolate hydrolysis based on the amino acid sequences deduced from the isolated cDNA sequences. Brassica oleracea consensus (cabbage, broccoli, and cauliflower), Brassica rapa ssp. perkinesis (ACO57702.1), Brassica napus (ACO57703.1), and Arabidopsis thaliana ESM1 (ABB90255.1) used to construct phylogenetic tree. The values in parenthesis are amino acid sequence similarity with B. oleracea consensus by using NCBI BLAST search. The tree was constructed using Clustal W2 (http://www.ebi.ac.uk/Tools/clustalw2/).