pmc logo imageJournal ListSearchpmc logo image
Logo of plntphysJournal URL: redirect3.cgi?&&auth=0nOOfJWdapggFMBxYA2_fHyzCL1ia-dZ4HQzLukFP&reftype=publisher&article-id=2330308&issue-id=164620&journal-id=69&FROM=Article|Banner&TO=Publisher|Other|N%2FA&rendering-type=normal&&http://www.plantphysiol.org
Plant Physiol. 2008 May; 147(1): 415–428.
doi: 10.1104/pp.107.115279.
PMCID: PMC2330308
Genotype, Age, Tissue, and Environment Regulate the Structural Outcome of Glucosinolate Activation1[W][OA]
Adam M. Wentzell and Daniel J. Kliebenstein*
Genetics Graduate Group and Department of Plant Sciences, University of California, Davis, Davis, California 95616
*Corresponding author; e-mail kliebenstein/at/ucdavis.edu.
Received December 20, 2007; Accepted March 16, 2008.
Abstract
Glucosinolates are the inert storage form of a two-part phytochemical defense system in which the enzyme myrosinase generates an unstable intermediate that rapidly rearranges into the biologically active product. This rearrangement step generates simple nitriles, epithionitriles, or isothiocyanates, depending on the structure of the parent glucosinolate and the presence of proteins that promote specific structural outcomes. Glucosinolate accumulation and myrosinase activity differ by plant age and tissue type and respond to environmental stimuli such as planting density and herbivory; however, the influence of these factors on the structural outcome of the rearrangement step remains unknown. We show that the structural outcome of glucosinolate activation is controlled by interactions among plant age, planting density, and natural genetic variation in Arabidopsis (Arabidopsis thaliana) rosette leaves using six well-studied accessions. We identified a similarly complex interaction between tissue type and the natural genetic variation present within these accessions. This raises questions about the relative importance of these novel levels of regulation in the evolution of plant defense. Using mutants in the structural specifier and glucosinolate activation genes identified previously in Arabidopsis rosette leaves, we demonstrate the requirement for additional myrosinases and structural specifiers controlling these processes in the roots and seedlings. Finally, we present evidence for a novel EPITHIOSPECIFIER PROTEIN-independent, simple nitrile-specifying activity that promotes the formation of simple nitriles but not epithionitriles from all glucosinolates tested.
 
Plants synthesize thousands of secondary metabolites, including glucosinolates, alkaloids, terpenoids, and phenylpropanoids. These diverse chemicals often serve as the primary medium for plant-ecosystem interactions, attracting pollinators and seed dispersers and providing defense against herbivores and pathogens (Wink, 1988). Many secondary metabolites have demonstrated defense functions as growth inhibitors, toxins, or feeding deterrents for myriad herbivores and pathogens, ranging from phloem-feeding aphids to leaf-chewing lepidopteran larvae (Burow et al., 2006b; Kim and Jander, 2007). However, adapted herbivores have evolved strategies to detoxify or avoid plant defense systems and often use these phytochemicals to locate appropriate plant hosts. For example, several insect species specialize on glucosinolate-containing plants and use particular glucosinolate structures as cues for host recognition, feeding, and oviposition (Giamoustaris and Mithen, 1995; Griffiths et al., 2001; Moyes and Raybould, 2001; Miles et al., 2005). These opposing selective pressures imposed by generalist and specialist herbivores likely favor the structural diversification of defensive metabolites, as plants seek to avoid detection by adapted specialists yet maintain effective defenses against polyphagous herbivores (Wink, 2003; Lankau, 2007).
Population density can impact secondary metabolite defenses as a consequence of specialist herbivores using these metabolites to locate their hosts. A single isolated plant may not produce sufficient volatile metabolite to attract specialist herbivores and therefore would benefit from an emphasis on direct defenses. However, the combined volatiles of numerous plants producing the same metabolites will increase the risk of attracting specialist herbivores, thus favoring structures that are weaker attractants despite providing less effective defense against polyphagous insects (Fig. 1C). As such, plastic regulation of defensive metabolite structures in response to conspecific density could be advantageous when the defensive metabolite is also a volatile host recognition cue for adapted herbivores (Lankau, 2007).
Figure 1.Figure 1.
Rearrangement during glucosinolate activation can generate multiple structural outcomes and mediate plant-herbivore interactions. A, Variable structural outcomes across plant development could optimize a plant's defensive chemistry to temporal changes (more ...)
Extrinsic factors controlling herbivore populations may also influence secondary metabolite diversity. For instance, seasonal changes can alter the herbivore community over the course of plant development (Wolda, 1988); therefore, selection might favor defensive metabolite systems capable of regulating structural variation over an individual plant's life (Fig. 1A). Because most insect populations are composed of multiple feeding guilds, they may drive metabolic defenses to respond to differences in insects' preferred feeding tissues (Burow et al., 2006b; Kim and Jander, 2007). If metabolites that provide defense against one guild are ineffective against herbivores with different feeding preferences, defensive metabolite systems capable of changing structural variation across tissues would be favored (Fig. 1B). Here, we focus on the structural outcome of glucosinolate activation in Arabidopsis (Arabidopsis thaliana) as a model secondary metabolite system to test whether structural diversity is modulated in response to development, tissue, and population density, consistent with the above hypotheses.
Glucosinolates are structurally diverse metabolites produced throughout the Brassicaceae, including the genetic model plant Arabidopsis (Fahey et al., 2001; Wittstock and Halkier, 2002). Glucosinolates are the inert storage form of a two-part defense system in which the biologically active structures are catabolically produced by the enzyme myrosinase (Louda and Mole, 1991; Bones and Rossiter, 1996). Structural variation in intact glucosinolates is achieved through the use of different amino acid precursors and sequential modification of the side chains. Natural genetic variation in the biosynthetic enzymes within Arabidopsis leads to diverse chemical mixtures influencing insect herbivory (Kliebenstein et al., 2001a, 2001b, 2001c). The accumulation and structural diversity of intact glucosinolates are regulated by plant development and tissue type (Petersen et al., 2002; Reichelt et al., 2002; Brown et al., 2003) and a suite of biotic and abiotic factors (Leach et al., 1999; Agrawal et al., 2002; Kliebenstein et al., 2002; Bidart-Bouzat et al., 2005; Mewis et al., 2006; Kim and Jander, 2007; Reifenrath and Muller, 2007).
Further structural diversity is introduced into the glucosinolate system during catabolic glucosinolate activation. The enzyme myrosinase generates an unstable intermediate that subsequently rearranges into the different glucosinolate activation products, with divergent biological activities dependent upon the specific structures produced (Fig. 1D; Louda and Mole, 1991; Bones and Rossiter, 1996; Nastruzzi et al., 1996). The principal structural outcomes of glucosinolate activation are simple nitriles, epithionitriles, or the default isothiocyanate (Fig. 1D). Structural outcomes in planta depend on both the structure of the parent glucosinolate and natural genetic variation for the presence of structural specifier proteins, such as EPITHIOSPECIFIER PROTEIN (ESP) and EPITHIOSPECIFIER MODIFIER1 (ESM1; Fig. 1D; Bernardi et al., 2000; Foo et al., 2000; Lambrix et al., 2001; Zabala et al., 2005; Burow et al., 2006a; Zhang et al., 2006). Structural outcomes are often reported as the percentage of a particular structure relative to all of the products formed from the parent glucosinolate (Lambrix et al., 2001; Burow et al., 2006a, 2006b, 2007; Zhang et al., 2006). This allows the analysis of biochemical processes that generate a particular structure with minimal influence from any variation in myrosinase activity or intact glucosinolate accumulation. In spite of the biological importance of the structural outcome of glucosinolate activation, it remains to be determined how this process is regulated during development and in different tissues.
In this study, we test the role that plant age, growth conditions, and tissue types play in modulating metabolite diversity using the structural outcome of glucosinolate activation in Arabidopsis as our model. We demonstrate that there is quantitative age-dependent regulation of structural outcomes that is modulated by planting density and that exhibits natural genetic variation using six well-characterized accessions. We also detected differential regulation of structural outcomes between roots, flowers, and seedlings. As previous research focused on the structural outcome of glucosinolate activation in mature rosettes, we used previously described mutants in rosette myrosinases and structural specifier genes to demonstrate that glucosinolate activation and the control of structural outcomes in roots and seedlings use distinct myrosinases and structural specifiers, which remain to be identified. This analysis of control over structural variation in the glucosinolate system should aid future studies designed to test whether this variation optimizes plant defenses in response to heterogeneous herbivore populations across different individual and community scales.
RESULTS
Does the Control of Structural Outcomes Show Age-Dependent Regulation?
Previous research on the control of structural rearrangement following glucosinolate activation in Arabidopsis has focused primarily on rosette leaves at a single stage in plant development. To determine whether the structural outcome of glucosinolate activation depends on the age of the plant, we determined the proportion of glucosinolate products in each structural class at weekly intervals in rosette leaves in six accessions.
Developmental Trajectory of Simple Nitrile Production in Columbia
Over the course of development in the common reference accession Columbia (Col-0), we observed significant changes in the proportion of simple nitrile structures produced. Both exogenous allyl glucosinolate and endogenous 4-methylsulfinylbutyl glucosinolate showed a 3-fold increase in simple nitrile proportions from 3 weeks after germination to a maximum at 8 weeks (Fig. 2A; Table I). The absolute magnitude of simple nitrile proportions for the two substrates was not identical, but similar differences in structural outcomes using diverse substrates have been observed previously (Lambrix et al., 2001; Zhang et al., 2006) and may indicate distinct rearrangement kinetics depending on parent glucosinolate structure. The simple nitrile proportions observed from weeks 4 through 8 are significantly higher than those in the controls using Sinapis alba myrosinase with allyl glucosinolate as the substrate.
Figure 2.Figure 2.
Structural outcomes of glucosinolate activation in Col-0 rosettes as a function of age and growth conditions. Percentages of simple nitrile (white circles) and isothiocyanate (black triangles) were determined for both the exogenously supplied allyl glucosinolate (more ...)
Table I.Table I.
Simple nitrile proportions in Col-0 depend on age and growth conditions
Genetic Variation in Simple Nitrile Production
To determine whether the developmental trajectory of simple nitrile proportions observed in Col-0 was generalizable, five additional accessions were analyzed (Table II; Fig. 3A). Most of these diverse accessions had simple nitrile proportion trajectories similar to those of Col-0 (Fig. 3A). However, there were significant differences between accessions for this trajectory, highlighted by Bayreuth (Bay-0), which had the largest increase in simple nitrile proportions, and Monte Tosso (Mr-0), with the slowest increase during development (Fig. 3A; Table II; accession [ESP] × week). Previous analysis associated genetic polymorphism at ESP with both simple nitrile and epithionitrile formation; however, we did not find a significant effect of ESP genotype on the proportion of simple nitrile structures produced in this collection of accessions. This suggests that the observed simple nitrile production is distinct from ESP-mediated nitrile production and may be controlled by different structural specifier genes (Table II). The trajectory of simple nitrile proportions for endogenous but-3-enyl glucosinolate in Cape Verde Islands (Cvi) and Shahdara (Sha) is nearly identical to that for the exogenous allyl substrate (Supplemental Fig. S1A; Supplemental Table S1).
Table II.Table II.
GLIMMIX analysis of the developmental trajectories of structural outcomes for allyl glucosinolate
Figure 3.Figure 3.
Developmental trajectories of structural outcomes for six accessions. Structural outcomes were determined in rosette leaves, and the average percentages of simple nitrile (white circles), epithionitrile (black squares), and isothiocyanate (black triangles) (more ...)
The Developmental Trajectory of Epithionitrile Production
In addition to simple nitriles, the accessions Cvi, Landsberg erecta (Ler), and Sha possess ESP activity leading to the formation of epithionitrile structures from allyl glucosinolate. We tested whether the proportion of epithionitriles produced in these three accessions also depends on plant age (Fig. 3A; Table II). From a low at 3 weeks after germination, epithionitrile proportions increased gradually throughout development until it was the predominant structure produced following allyl or but-3-enyl glucosinolate activation (Fig. 3A; Table II, Supplemental Fig. S1A; Supplemental Table S1). In contrast to simple nitriles, the developmental trajectories for epithionitrile proportions did not vary significantly between Cvi, Ler, and Sha (Table II; accession × week). The proportion of isothiocyanates formed was affected by variation in both plant age and natural genetic diversity, with similar patterns between endogenous and exogenous glucosinolates (Table II; Fig. 3A; Supplemental Fig. S1A). The directionality of the changes agrees with the concept that isothiocyanate production is the default outcome of glucosinolate hydrolysis and is reduced as substrate is diverted into simple nitrile and epithionitrile structures (Fig. 3A).
Structural Outcomes and Rosette Leaf Age
The analysis of developmental effects on structural outcomes of glucosinolate activation could be complicated by changes controlled by the ages of individual leaves within a rosette. To test for a relationship between leaf age and structural outcomes, a leaf was sampled from the innermost, middle, and outermost whorls of the rosettes of 8-week-old plants. The leaf sampled from the innermost whorl was still expanding, while the oldest leaf showed signs of senescence. Most accessions showed no significant variation in endogenous or exogenous glucosinolate structural outcomes between these three leaf cohorts (Supplemental Fig. S2). This shows that variation in the structural outcomes of glucosinolate activation between leaf cohorts within the rosette cannot explain the observed developmental trajectories (Fig. 3A). In particular, the youngest leaves of Cvi, Ler, and Sha showed high proportions of epithionitrile structures, while young plants showed low epithionitrile proportions, suggesting that plant age is the main determinant of the structural outcome of glucosinolate activation in our analysis (Fig. 3A; Supplemental Fig. S2). The only accession showing variation between differently aged leaves within a rosette was Cvi (Supplemental Fig. S2; P = 0.034 and 0.022, respectively).
Does the Developmental Trajectory of Structural Outcomes Depend on Planting Density?
As Arabidopsis in the wild can grow as isolated plants or in more dense populations, the structural outcome of glucosinolate activation might be regulated by the presence of neighboring plants, using either direct signals from the neighboring plants or indirect competition for space, light, or nutrients. To test this, structural outcomes were assayed in plants grown at high and low densities (Fig. 3; Table II).
Density Modulates Simple Nitrile Production in Col-0
Within the reference accession, Col-0, the structural outcome of glucosinolate activation depends strongly on planting density and plant age (Table I). When grown at high density, simple nitrile proportions in Col-0 were highest at 2 weeks and declined with age. In contrast, simple nitrile proportions in plants grown at low density started low and increased gradually over time (Table I; Fig. 2B). This suggests that Arabidopsis can modulate structural outcomes in response to planting density. Future experiments are required to identify the critical aspects of planting density that control this plasticity, such as increased competition for nutrients or the physical presence of neighbors.
Genetic Variation in Developmental Trajectories in Response to Planting Density
We next tested whether the observed response to planting density in Col-0 is found in all six of these accessions. In general, all accessions showed higher simple nitrile proportions in young plants at high density compared with the low-density planting (Table II; Fig. 3). Interestingly, there was a statistically significant effect of the accession (ESP) × density term, indicating that there are genetic differences among these accessions in their response to planting density (Table II). This interaction is most obvious in the accession Bay-0, in which strong increases in simple nitrile proportions were observed in older plants at high density in comparison with the other accessions (Fig. 3, Bay-0 versus Col-0).
Epithionitrile proportions were also affected by planting density, with high density accelerating the age-dependent increase in exogenous and endogenous epithionitrile proportions in both Cvi and Sha (Tables I and II; Figs. 2 and 3; Supplemental Fig. S1; Supplemental Tables S1 and S2). In contrast, simple nitrile proportions decline over time in these accessions at elevated density. The observation of distinct trajectories for simple nitrile and epithionitrile proportions also supports the idea that simple nitrile formation following glucosinolate activation is not an exclusively ESP-mediated process.
Structural Outcomes and the Onset of Bolting
We noted flowering times for each accession to investigate a possible connection between the regulation of the structural outcome of glucosinolate activation in rosette leaves and the transition to reproductive growth. Under our growth conditions, bolting in Cvi was accelerated at high density, with bolting at 6 and 7 weeks in the high- and low-density plantings, respectively. In contrast, bolting was delayed by 1 week at high density in Bay-0 and Sha, with bolting at 8 and 7 weeks in high and low density. Mr-0, Ler, and Col-0 bolted after the conclusion of this experiment. We observed no relationship between structural outcomes and the onset of bolting in Cvi, Sha, and Bay-0. For example, in Bay-0, bolting was delayed by high density, but the age-dependent increase in simple nitrile proportion was accelerated (Fig. 3). Furthermore, Cvi and Sha showed similar responses to density throughout development for both simple nitrile and epithionitrile proportions, despite the different effects of density on flowering time (Fig. 3). This suggests that the developmental regulation of structural outcomes is not a secondary effect of the transition to flowering.
Do Structural Outcomes Vary Among Tissues?
The preferences of different insect guilds for specific tissues could lead to tissue-specific regulation of the structural outcome of glucosinolate activation. To determine whether structural outcomes vary between plant tissues, we compared rosettes, flowers, cauline leaves, seedlings, and roots in these accessions.
Structural Outcomes in Cauline Leaves and Florets
Flowers showed a pronounced decrease in epithionitrile proportions relative to the rosettes. All three accessions containing functional ESP went from 80% to 90% epithionitrile products in the rosettes to less than 40% in flowers (Fig. 4). The observed declines in epithionitrile production in flowers were balanced by increased proportions of isothiocyanate structures. In contrast, cauline leaf epithionitrile proportions were statistically indistinguishable from those in concurrently sampled rosette tissue. Cvi was the lone exception, with increased epithionitrile proportions in the cauline leaves (Fig. 5; P = 0.041). Of these six accessions, only Bay-0 showed a significant decrease in simple nitrile proportions between rosette and cauline leaves and floral organs (Fig. 4, rosettes versus cauline leaves [P = 0.002], cauline leaves versus flowers [P = 0.001]). Endogenous glucosinolate activation in Col-0, Cvi, and Sha exhibited similar patterns of structural outcomes as observed for allyl glucosinolate (Supplemental Fig. S3). Thus, the structural outcome of glucosinolate activation is differentially regulated between flowers and leaves in Arabidopsis, which may indicate that simple nitriles and epithionitriles are less effective at providing defense against insects that attack Arabidopsis flowers.
Figure 4.Figure 4.
Changes in structural outcomes in flowers and cauline leaves. Glucosinolate activation was assayed in flowers, cauline leaves, and rosette leaves using allyl glucosinolate in four plants of each accession. Average percentages of simple nitrile (dark shading), (more ...)
Figure 5.Figure 5.
Structural outcomes in seedlings. Structural outcomes were measured using allyl glucosinolate in 3- and 7-d-old seedlings and are reported as average percentages of a given structure. The accessions are Bay-0 (♦), Col-0 (□) Mr-0 ([filled triangle]), (more ...)
Structural Outcomes in Seedlings
Arabidopsis seedlings have a different glucosinolate composition than mature tissues, and we proceeded to test whether they also produce different structures following glucosinolate activation (Petersen et al., 2002; Brown et al., 2003). All accessions had high simple nitrile proportions at 3 and 7 d after germination, with higher levels of simple nitrile structures at 7 d than at any age of rosette leaf (Fig. 3). Simple nitrile proportions had a statistically significant accession (ESP) × day interaction, indicating differences in the temporal control of simple nitrile formation in seedlings, with some accessions showing a rapid increase in simple nitrile proportions and others showing little to no change (Fig. 5A; Table III). Similar to the rosette, ESP did not affect simple nitrile proportions in seedlings (Table III).
Table III.Table III.
GLIMMIX analysis of structural outcomes in seedlings
Epithionitrile proportions were low in 3-d-old seedlings and increased from 3 to 7 d in Ler and Cvi but not in Sha (Fig. 5B; Table III). Endogenous glucosinolate structural outcomes exhibited the same trends observed for allyl glucosinolate in Col-0, Cvi, and Sha seedlings, although the magnitude of simple nitrile proportions for 4-methylsulfinylbutyl glucosinolate was considerably lower (Supplemental Fig. S4). Col-0 and Ler seedlings tested using the exogenous benzyl glucosinolate showed the same trend of high simple nitrile proportions at 3 and 7 d (data not shown). These data suggest that juvenile Arabidopsis seedlings tend to form more simple nitrile structures relative to mature rosettes. It is possible that this may play a role in the turnover of endogenous glucosinolates during seed germination (Petersen et al., 2002), as the less toxic simple nitrile structures are more readily metabolized.
Glucosinolate Activation in Roots
This assay revealed very high proportions of simple nitriles in the roots of 2-week-old plants in all accessions, with no detectable epithionitriles and isothiocyanates forming the remainder of the observed products (Fig. 6). This was not an artifact of the Murashige and Skoog (MS) medium, as the observed simple nitrile proportions in roots were at least 10-fold higher than those observed in control reactions using S. alba myrosinase with added MS medium. The absence of epithionitrile structures further supports the observation that simple nitrile production is distinct from ESP-mediated processes. The structural outcome of glucosinolate activation in Col-0 and Ler root tissue was also tested using benzyl glucosinolate as the substrate (Fig. 6B). No significant difference in simple nitrile proportions was observed for this substrate relative to allyl glucosinolate (P = 0.220 for Ler and 0.072 for Col-0). These data show that Arabidopsis roots have a different structural outcome relative to the rosettes, supporting the concept that each tissue may be optimized for the particular challenges it confronts.
Figure 6.Figure 6.
Structural outcomes in roots. The structural outcome of glucosinolate activation was assayed using allyl glucosinolate in the roots of 2-week-old plants. A, Representative gas chromatography-flame ionization detection chromatogram of glucosinolate activation (more ...)
Do Known Structural Specifier Genes Control Structural Outcomes in Roots and Seedlings?
All previous genetic and biochemical analyses of the structural outcome of glucosinolate activation in Arabidopsis focused on mature rosette tissue. However, the previous data show that each tissue produces a different mixture of glucosinolate activation products. Therefore, we tested the known rosette leaf structural specifier genes and myrosinases for a role in glucosinolate activation in root and seedling tissues.
ESP-Independent Simple Nitrile Formation
To test directly whether simple nitrile formation is independent of ESP, we assayed the seedlings and roots of an ESP knockout line and a 35S::ESP overexpression line (Fig. 7). Simple nitrile proportions in the ESP knockout line did not differ significantly from those of wild-type Col-0 in seedlings or roots for either the endogenous or exogenous glucosinolates, demonstrating that ESP is not required for this activity (Fig. 7; Supplemental Fig. S5). This was further supported by simple nitrile production in seedlings and roots of the accession Mr-0, which contains a large deletion in the ESP open reading frame (Fig. 5; Lambrix et al., 2001). The 35S::ESP line, which increases the production of simple nitriles from endogenous 4-methylsulfinylbutyl glucosinolate in Arabidopsis rosettes, produced high levels of epithionitriles from allyl glucosinolate in seedlings and roots, showing that ESP can function in these tissues. Boiled root controls showed less than 2% simple nitriles, demonstrating that the root-specific simple nitrile structural specifier activity is heat labile and likely attributable to unknown structural specifier proteins.
Figure 7.Figure 7.
Structural outcomes in glucosinolate activation mutants. Structural outcomes were measured using allyl glucosinolate in the roots and seedlings of mutants with alterations in the myrosinases and structural specifier genes known to function in rosette (more ...)
Structural Outcomes in ESM1 Mutants
ESM1 is another glucosinolate structural specifier gene, which inhibits simple nitrile formation in rosettes with functional ESM1 (Zhang et al., 2006). ESM1 knockout seedlings had elevated simple nitriles relative to the parental Col-0 at 7 d after germination (P = 0.004) but not at 3 d, indicating that ESM1 does not affect simple nitrile proportions until after day 3 (Fig. 7A; Supplemental Fig. S5). In contrast, simple nitrile proportions were unchanged in the roots of the ESM1 knockout (Fig. 7C). The ESM1 overexpression line was created in a recombinant inbred line, CS1995, derived from a cross between Col-0 and Ler (Zhang et al., 2006). No significant changes in simple nitrile or epithionitrile proportions were detected between CS1995 and 35S::ESM1 lines in either seedlings or roots (Fig. 7; Supplemental Fig. S5). This contrasts with the rosettes, in which both the ESM1 knockout and 35S lines alter simple nitrile proportions (Zhang et al., 2006). Thus, ESM1 does not inhibit the unidentified simple nitrile specifier activity in Col-0 seedlings or roots.
Glucosinolate Activation and Structural Outcomes in TGG1/2 Mutants
The final known components of glucosinolate activation in Arabidopsis rosettes are the myrosinase isoenzymes TGG1 and TGG2 (Barth and Jander, 2006). To test whether these myrosinases are responsible for glucosinolate activation in seedlings and roots, glucosinolate activation products were assayed in the TGG1/2 double knockout line (Fig. 7). Although the activation assay is not a quantitative measurement of myrosinase activity, nearly wild-type levels of activation products were detectable in seedlings and roots of the double knockout line, whereas no activation products were detectable in 5-week-old rosette leaves of this line (data not shown). Furthermore, the TGG1/2 mutant had no change in structural outcomes of glucosinolate activation in roots or 3-d-old seedlings and had only a slight but significant decline in simple nitrile production at 7 d (Fig. 7; Supplemental Fig. S5; P = 0.031).
DISCUSSION
The variable biological activities of the different glucosinolate activation product structures, and evidence of complex regulatory patterns in glucosinolate accumulation and myrosinase activity (Petersen et al., 2002; Brown et al., 2003; Barth and Jander, 2006), led us to hypothesize that structural outcomes may be subject to diverse developmental and environmental regulatory inputs. To test this, we assayed the structural outcome of glucosinolate activation in rosette leaves throughout development in six accessions grown at two different planting densities and in seedlings, flowers, cauline leaves, and roots. Our results reveal previously unknown levels of regulation in structural outcomes, dependent on plant age, tissue, and planting density. Microarray transcript profiling has shown that genes involved in glucosinolate activation and controlling structural outcomes respond to a variety of stimuli, including herbivory, fertilizer, and light regimen (Reymond et al., 2004; Hirai et al., 2005; Maruyama-Nakashita et al., 2006; Dombrecht et al., 2007). However, it was not tested directly whether this gene regulation results in age-, tissue-, or environment-dependent regulation of structural outcomes. Recently, the leaf myrosinases and ESP were also shown to have distinct patterns of regulation between specific plant tissues that varied among Arabidopsis accessions (Barth and Jander, 2006; Burow et al., 2007). Assuming that the primary role of the glucosinolate system is plant defense, regulation of structural outcomes might allow plants to optimize their defenses in response to changing risks of herbivory.
Optimization of Herbivory Defense
Striking differences in the structural outcome of glucosinolate activation were observed over the course of plant development and among tissues tested. All accessions produced primarily isothiocyanates in the flowers and high levels of simple nitriles in the roots and seedlings (Figs. 5–7). This suggests important differences in the defensive needs of these tissues, possibly as a consequence of the distinct herbivore communities they confront. Isothiocyanates may provide better defense against pollen-feeding insects in flowers, whereas simple nitriles seem to be preferred in the roots. Additionally, the age of the plant modulates structural outcomes in rosette leaves (Figs. 2 and 3), possibly in response to seasonal changes in the herbivore community. This age-dependent regulation was affected by planting density, because crowding appeared to increase simple nitrile and epithionitrile proportions, particularly in younger plants. Since isothiocyanates are potent attractants and oviposition cues for specialized herbivores, the distinct trajectories observed at different planting densities may reflect the changing risk of attack by polyphagous and specialist herbivores. Detection and oviposition by specialists can lead to sustained herbivory by the developing insect larvae, and young plants may be unable to tolerate this damage. As such, a single isolated plant, which is not producing enough total volatiles to be at risk of detection, may prefer the more attractive but also more effective herbivory-deterring isothiocyanate structures. In contrast, younger plants at high densities might emphasize less attractive nitrile structures to avoid detection, despite reduced defensive potential toward generalist herbivores. Now that the patterns of regulation have been described, the stage is set for detailed studies to directly test these ecological hypotheses and investigate the regulatory networks involved.
Interactions of Natural Genetic Variation and Development
Significant natural variation was identified among the six Arabidopsis accessions tested in tissue-, age-, and density-dependent regulation of the structural outcomes of glucosinolate activation. Nearly every tissue or condition tested showed evidence for control by natural genetic variation (Figs. 3–7; Tables II and III). These changes in structural outcomes could be mediated by the differential regulation of a single gene. This is clearly supported by the observed tissue-specific regulation of ESP, which is tightly correlated with the production of epithionitriles (Burow et al., 2007). In contrast, interactions between natural variation and the regulation of simple nitrile proportions suggest that this may be a modular system, such that natural variation can affect the percentage of simple nitriles in one tissue without affecting other tissues or conditions. This is supported by the observation that structural specifier genes discovered in rosette leaves do not play a significant role in simple nitrile formation in roots and seedlings (Fig. 7). Further support for modularity comes from the observation that structural outcomes in CS1995 roots are identical to those in the Col-0 parent, yet in shoots, CS1995 is identical to its Ler parent and distinct from Col-0 (Fig. 7). Such modularity could be mediated by members of a gene family each acting in a limited set of tissues, a result of an evolutionary process called subfunctionalization (Hughes, 1994; Lynch and Force, 2000; Fraser et al., 2002). There are numerous Arabidopsis genes homologous with ESP, ESM1, and the myrosinase binding proteins, all with different tissue-specific expression patterns (Capella et al., 2001; Zhang et al., 2002; Barth and Jander, 2006). However, the homologs of ESP and ESM1 have no demonstrated function, so further experiments are required to test their importance in the tissue-specific regulation of structural outcomes.
Identification of a Novel Simple Nitrile Specifier Activity
To date, ESP is the only known structural specifier protein in plants capable of increasing the production of epithionitriles from alkenyl glucosinolates and simple nitriles from all other substrates (Tookey, 1973; Bernardi et al., 2000; Foo et al., 2000; Lambrix et al., 2001; Zabala et al., 2005; Burow et al., 2006a, 2007). Our study identified one or more novel simple nitrile specifier activities that are genetically and biochemically distinct from ESP. The strongest evidence came from root and seedling simple nitrile production in the absence of a functional ESP in Mr-0, which has a naturally occurring deletion (Lambrix et al., 2001), and the T-DNA knockout of ESP in Col-0 (Figs. 6–8). All accessions containing functional ESP produced simple nitriles in the roots, not epithionitriles, as would be expected if ESP were mediating structural outcomes (Fig. 6). Detection of epithionitriles in roots of the ESP overexpression lines confirmed that epithionitrile production is possible in Arabidopsis roots (Fig. 7C). We also observed high proportions of simple nitriles at various ages in shoot tissues of all accessions, which was statistically independent of ESP functionality (Figs. 3 and 6; Tables II and III). Combined, this evidence strongly supports the existence of a novel simple nitrile structural specifier activity that exclusively promotes the formation of simple nitriles, but not epithionitriles, regardless of the parent glucosinolate structure.
CONCLUSION
We have described several patterns of age- and tissue-dependent regulation of the structural outcome of glucosinolate activation and interactions with natural variation, setting the stage for detailed molecular studies to resolve the basis of this variation. These novel regulatory inputs indicate that Arabidopsis regulates the mixture of activation product structures over tissue, development, and planting density. However, the ecological roles of many of these defensive metabolites remain undefined, requiring further experiments to understand the ecological significance of these results. In particular, the biological role of simple nitrile activation products remains unclear, as simple nitrile structures are generally less toxic and provide less effective defense against polyphagous herbivores when compared directly with the corresponding isothiocyanates (Nastruzzi et al., 2000; Jander et al., 2001; Lambrix et al., 2001; Barth and Jander, 2006; Zhang et al., 2006). We hypothesize that increased proportions of simple nitrile and epithionitrile products reduce the risk of detection by adapted herbivores or provide indirect defense by attracting predators and parasites of adapted herbivores. The use of knockout and overexpression mutants with altered regulation of structural outcomes of glucosinolate activation will clarify the ecological significance of these results.
MATERIALS AND METHODS
Glucosinolate Activation Assay
The structural outcome of glucosinolate activation was assayed using a modified version of the previously published protocol (Lambrix et al., 2001). Briefly, 100 mg of Arabidopsis (Arabidopsis thaliana) tissue was harvested, weighed, and immediately crushed in an 8-mL reaction vial containing 1 mL of MES acid buffer at pH 6.0 and 0.4 μmol of allyl glucosinolate. Since the accessions and tissues used in this study vary in their endogenous glucosinolate content (Kliebenstein et al., 2001c; Petersen et al., 2002; Brown et al., 2003; Wentzell et al., 2007), exogenous allyl glucosinolate enabled comparisons of structural outcomes using a common substrate. For the flowers and roots and some of the earlier time points in the developmental trajectory, tissue was limiting. Analysis of samples ranging from 20 to 200 mg showed that the tissue mass sampled did not significantly affect structural outcomes (data not shown). Upon complete tissue homogenization, the reaction vial was capped and incubated for 5 min. The reaction was stopped and glucosinolate activation products were extracted with 4 mL of dichloromethane. The organic phase was removed, dried, and concentrated to 200 μL for gas chromatography analysis using an Agilent HP 5890 apparatus with a flame ionization detector (Lambrix et al., 2001). Peak identities were confirmed using a gas chromatograph-mass spectral detector (Agilent HP 6890 and Agilent 5973N MSD), by comparison with published mass spectra (Spencer and Daxenbichler, 1980).
Structural outcomes are reported as percentages of simple nitrile, epithionitrile, or isothiocyanate products for a particular glucosinolate. By dividing the absolute amount of a particular structure by the sum of all possible products, the effects of myrosinase activity and differences in the biosynthesis and accumulation of the endogenous substrates are canceled, since they affect both the numerator and denominator equally. This assay is not a quantitative measure of total myrosinase activity, both because three of the accessions used in this study produce allyl glucosinolate and, therefore, the substrate concentration is not constant and because this assay reaches saturation for some samples.
Commercially available myrosinase purified from Sinapis alba (Sigma-Aldrich catalog no. T4528) was assayed for structural outcomes in the absence of plant tissues using allyl glucosinolate as the substrate. To control for abiotic factors influencing nonenzymatic production of simple nitrile and epithionitrile product structures, activation was assayed in the presence of MS medium and soil with and without S. alba myrosinase (Tookey, 1973; Foo et al., 2000). No glucosinolate activation was observed without myrosinase. The soil and MS medium controls with added myrosinase showed 3% and 6% simple nitrile products and no epithionitriles.
Plant Material
Six accessions of Arabidopsis were analyzed to test for natural genetic variation in the regulation of structural outcomes of glucosinolate activation. The accessions Bay-0, Col-0, Cvi, Ler, and Sha were selected because they are the parents of three large recombinant inbred line populations (Lister and Dean, 1993; Alonso-Blanco et al., 1998; Loudet et al., 2002) and have been characterized for their glucosinolate content and glucosinolate activation in rosette leaves (Kliebenstein et al., 2001c; Lambrix et al., 2001; Wentzell et al., 2007). The accessions Ler, Cvi, and Sha all possess functional ESP and produce predominantly epithionitriles and simple nitriles following glucosinolate activation in 4-week-old rosette leaves, whereas Bay-0, Col-0, and Mr-0 lack functional ESP and produce primarily isothiocyanates (Lambrix et al., 2001). Mr-0 has a large deletion in the open reading frame of ESP, creating a natural knockout, whereas the loss of expression and a splice site defect are the causal polymorphisms in Col-0 and Bay-0 (Lambrix et al., 2001). Seeds were obtained from the Arabidopsis Biological Resource Center (ABRC).
Mr-0, Cvi, and Sha produce allyl glucosinolate, while Cvi and Sha also accumulate but-3-enyl glucosinolate in rosette leaves (Kliebenstein et al., 2001c; Wentzell et al., 2007), and all three but-3-enyl activation products were detectable. Col-0 produces 4-methylsulfinylbutyl glucosinolate in the rosettes, and its simple nitrile and isothiocyanate products were detected (Kliebenstein et al., 2001c). Under these conditions, the products from endogenous 3-hydroxypropyl glucosinolate in Bay-0 and Ler could not be quantified reliably (Kliebenstein et al., 2001c; Wentzell et al., 2007). Activation products from endogenous glucosinolates in the roots were below the limit of detection in all accessions.
Growth Conditions for the Analysis of Age- and Planting Density-Dependent Regulation
Rearrangement outcomes were surveyed at weekly intervals in the rosette leaves of plants grown alone and with numerous conspecific neighbors. Plants were maintained at 25°C in a controlled-environment chamber with a 12-h light cycle under 150 μE fluorescent lighting. Seeds were cold stratified at 4°C for 48 h, then germinated and grown in Premier ProMix potting soil. Plants were watered and fertilized by subirrigation at weekly intervals. The low-density treatment consisted of a single plant per 49-cm2 pot. The high-density plants were maintained at approximately 100 plants per 49-cm2 pot, ensuring extensive contact and crowding between adjacent plants throughout the experiment. The high-density plants were maintained in the same chamber and grown at the same time as the low-density plants. Rosette leaves were assayed beginning at 2 and 3 weeks for the high- and low-density treatments, respectively.
Pots were removed from the growth chamber, and rosette tissue was harvested in three independent biological replicates sampled from separate pots for each time point, accession, and planting density. All replicates were grown in independent randomized complete blocks. At weeks 3 and 4 in the low-density treatment and at weeks 2 through 5 at high density, rosette leaves from multiple plants were pooled to obtain sufficient tissue for each replicate. At least three leaves were collected from each individual included in every replicate for all pooled samples. Trajectories were ended after 7 weeks for Ler and after 8 weeks for Col-0 and Mr-0 in the low-density treatment, and the high-density treatment of Ler was not included due to limited germination. Wounded plants were never returned to the growth chamber, and at no point during this experiment were the plants subjected to visible pests.
To test for an association between structural outcomes and leaf age, three leaves were sampled from 8-week-old plants of each accession and assayed individually. These plants were maintained at low density under conditions identical to the developmental trajectory and never subjected to visible pests. The youngest leaf was sampled from the innermost whorl and not yet fully expanded. A healthy, fully expanded leaf was sampled from the middle whorl of the plant. The oldest leaf was sampled from the outermost whorl and was beginning to senesce, with approximately 10% of the leaf area yellow. This was repeated using four plants of each accession grown simultaneously.
Growth Conditions for the Analysis of Structural Outcomes in Different Tissues
Structural outcomes were assayed in the flowers and cauline leaves using four plants of each accession grown at low densities in 49-cm2 pots and maintained in the same chambers and conditions as the developmental trajectory and free of pests throughout their life cycle. When the primary inflorescence reached 20 cm in height, all flowers from a single plant were collected and pooled to obtain sufficient tissue. Cauline and rosette leaves from these plants were assayed simultaneously for these comparisons. This provided four independent biological replicates per tissue per accession.
Root analysis required the growth of each accession on vertically oriented plates of MS medium (Gibco-BRL). Seeds were surface sterilized and cold stratified to ensure synchronous germination and maintained at 20°C under 160 μE fluorescent lighting. Root tissue from approximately 100 plants was pooled for each replicate, and three independent biological replicates taken from separate plates were assayed for each accession. Care was taken to exclude shoot tissue and any MS medium. Structural outcomes were also assayed using exogenous benzyl glucosinolate in Ler and Col-0. To control for the possibility that the presence of root chemical or ion concentrations were altering structural outcomes, boiled roots were incubated for 5 min with S. alba myrosinase and then assayed for structural outcomes.
For seedlings, each accession was germinated and grown on MS medium under conditions identical to those for the roots. At 3 and 7 d after germination, seedlings were assayed for structural outcomes. Approximately 50 seedlings were pooled for each sample, with three independent biological replicates sampled from different plates for each accession and time point. Care was taken to sample the aerial portions of the seedlings. All plates were free of pests and contamination by visual examination.
Analysis of Structural Outcomes in Known Mutants
We obtained previously described mutants in ESP, ESM1, and the myrosinases TGG1 and TGG2, the known glucosinolate activation and structural specifier genes identified in the rosettes, to determine whether they also function in roots and seedlings (Lambrix et al., 2001; Barth and Jander, 2006; Zhang et al., 2006). The confirmed ESP knockout line (Salk_010349) and the ESP overexpression line are in the Col-0 background (Lambrix et al., 2001; Burow et al., 2006a; Zhang et al., 2006). The ESM1 knockout line (Salk_043148) is in the Col-0 background (Zhang et al., 2006). The ESM1 overexpression line was created in a Col-0 × Ler recombinant inbred line (CS1995) that has a predominantly Col-0 genetic background with the functional Ler allele of ESP and the nonfunctional Ler allele of ESM1 (Zhang et al., 2006). To test the principal myrosinases in Arabidopsis rosette tissues, we obtained the TGG1/2 double knockout, also in the Col-0 background (Barth and Jander, 2006). Each mutant and control wild-type plant was assayed in three biologically independent replicates for each tissue as described above.
Developmental Trajectory Data Analysis
Although all three structural outcomes are contingent on one another, they each describe biochemically distinct outcomes, with potentially different regulatory patterns; thus, each was analyzed separately. Simple nitrile and isothiocyanate proportions were analyzed using a general linear mixed model, with density (high or low) and ESP genotype (functional or nonfunctional) as fixed factors, accession as a random factor nested within ESP, and age (in weeks) as a repeated measure to account for sampling these genotypes multiple times throughout development. Since epithionitrile formation was confined exclusively to accessions with functional ESP, analysis of epithionitrile proportions was restricted to Cvi, Ler, and Sha, and ESP genotype was dropped from the model. All interactions between factors were also tested (Table II).
Since wounding can alter structural outcomes, sampled plants were discarded after each assay. Thus, repeated sampling occurred at the level of accession rather than in individual plants. For each trait, we also investigated whether endogenous allyl glucosinolate synthesis in Cvi, Sha, and Mr-0 and the mass of tissue sampled affected structural outcomes by including each as a covariate in preliminary analyses. Neither had significant effects for any trait; thus, both were dropped from the final models.
As a result of excluding the high-density treatment of the accession Ler, and because the high-density treatment had enough tissue to be sampled starting at 2 weeks, as opposed to 3 weeks for the low density treatment, these data sets were unbalanced. Additionally, structural proportions were approximately log-normally distributed. To account for these departures from standard assumptions, the trajectory data were analyzed using the GLIMMIX procedure in SAS. This procedure minimizes the potential biases introduced by unbalanced and nonnormal data by using restricted penalized likelihood to fit the model and estimate random effects (Wolfinger and O'Connell, 1993).
To determine whether age, density, and their interaction have similar effects on exogenous and endogenous glucosinolate structural outcomes in Col-0, Cvi, and Sha, data for each substrate were analyzed within an accession according to the following model: proportion = age + density + age × density. The results obtained for each substrate were compared within each accession to determine whether the significance and directionality of each term were the same. Finally, to determine whether structural outcomes varied with leaf age within individual plants, differently aged leaves within a single rosette were compared within each accession using paired t tests. All post hoc comparisons were performed using Tukey adjustments.
Analysis of Structural Outcomes in Different Tissues and Mutants
Differences in structural outcomes between flowers, cauline leaves, and rosette leaves within an accession were analyzed using paired t tests within an accession. Root structural outcome data were compared between accessions, and within Col-0 and Ler between the benzyl and allyl glucosinolate substrates, using two-tailed heteroscedastic t tests, which allow for unequal sample variances. Post hoc comparisons were performed using Tukey adjustments.
Seedling structural outcome data were also analyzed using the GLIMMIX procedure. Each structural outcome was analyzed separately for the effects of days after germination, accession, ESP genotype (either functional or nonfunctional), and all possible interactions (Table III). Accession was nested within ESP genotype for simple nitrile and isothiocyanate proportions. For epithionitrile proportions, only the data from Cvi, Ler, and Sha were analyzed, and ESP genotype was removed from the model. Days were treated as repeated measures of the accessions. Post hoc comparisons were performed using Tukey adjustments.
Mutants in the known structural specifier and myrosinase genes identified in the rosettes were tested for their effects on structural outcomes using two-tailed heteroscedastic t tests. ESP overexpression and ESP knockout lines were compared with Col-0, as were the ESM1 knockout and the TGG1/2 double knockout lines. The ESM1 overexpression line was compared with the recombinant inbred line CS1995.
Supplemental Data
The following materials are available in the online version of this article.
  • Supplemental Figure S1. Endogenous glucosinolate structural outcomes in Cvi and Sha as a function of age and growth conditions.
  • Supplemental Figure S2. Leaf age-dependent changes in structural outcomes.
  • Supplemental Figure S3. Changes in endogenous glucosinolate structural outcomes in flowers and cauline leaves.
  • Supplemental Figure S4. Endogenous glucosinolate structural outcomes in seedlings of Col-0, Cvi, and Sha.
  • Supplemental Figure S5. Endogenous glucosinolate structural outcomes in seedlings of glucosinolate activation and structural specifier mutants.
  • Supplemental Table S1. Structural outcomes in Cvi depend on age and growth conditions.
  • Supplemental Table S2. Structural outcomes in Sha depend on age and growth conditions.
Supplementary Material
[Supplemental Data]
Acknowledgments
We thank Beverly Ajie (Center for Population Biology, University of California, Davis [UC Davis]) for assistance with the statistical analysis of the developmental data and many helpful discussions. We also thank Heather Rowe (Genetics Graduate Group, UC Davis) for critical review of the manuscript and Dr. Joshua Chang Mell (Genetics Graduate Group, UC Davis) for numerous helpful discussions. We are extremely grateful to Dr. Andrew Waterhouse and Mauri Anderson (Department of Viticulture and Enology, UC Davis) for the use of their gas chromatograph. Seeds of the TGG1/2 double knockout line were kindly provided by Dr. Georg Jander (Cornell University). Seeds of the 35S::ESP line were kindly provided by Dr. Ute Wittstock (Technische Universität Braunschweig). Benzyl glucosinolate was a generous gift of Dr. Steffen Abel (UC Davis).
Notes
1This work was supported by the National Science Foundation (grant nos. DBI 0642481 and MCB 0323759 to D.J.K.; Dissertation Improvement Grant DEB 0608516 to A.M.W.).
The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantphysiol.org) is: Daniel J. Kliebenstein (kliebenstein/at/ucdavis.edu).
[W]The online version of this article contains Web-only data.
[OA]Open Access articles can be viewed online without a subscription.
References
  • Agrawal AA, Conner JK, Johnson MT, Wallsgrove R (2002) Ecological genetics of an induced plant defense against herbivores: additive genetic variance and costs of phenotypic plasticity. Evolution Int J Org Evolution 56: 2206–2213 [PubMed]
  • Alonso-Blanco C, Peeters AJM, Koornneef M, Lister C, Dean C, van den Bosch N, Pot J, Kuiper MTR (1998) Development of an AFLP based linkage map of Ler, Col and Cvi Arabidopsis thaliana ecotypes and construction of a Ler/Cvi recombinant inbred line population. Plant J 14: 259–271 [PubMed]
  • Barth C, Jander G (2006) Arabidopsis myrosinases TGG1 and TGG2 have redundant function in glucosinolate breakdown and insect defense. Plant J 46: 549–562 [PubMed]
  • Bernardi R, Negri A, Ronchi S, Palmieri S (2000) Isolation of the epithiospecifier protein from oil-rape (Brassica napus ssp. oleifera) seed and its characterization. FEBS Lett 467: 296–298 [PubMed]
  • Bidart-Bouzat MG, Mithen R, Berenbaum MR (2005) Elevated CO(2) influences herbivory-induced defense responses of Arabidopsis thaliana. Oecologia 145: 415–424 [PubMed]
  • Bones AM, Rossiter JT (1996) The myrosinase-glucosinolate system, its organisation and biochemistry. Physiol Plant 97: 194–208.
  • Brown PD, Tokuhisa JG, Reichelt M, Gershenzon J (2003) Variation of glucosinolate accumulation among different organs and developmental stages of Arabidopsis thaliana. Phytochemistry 62: 471–481 [PubMed]
  • Burow M, Marekert J, Gershenzon J, Wittstock U (2006a) Comparative biochemical characterization of nitrile-forming proteins from plants and insects that alter myrosinase-catalysed hydrolysis of glucosinolates. FEBS J 273: 2432–2446 [PubMed]
  • Burow M, Muller R, Gershenzon J, Wittstock U (2006b) Altered glucosinolate hydrolysis in genetically engineered Arabidopsis thaliana and its influence on the larval development of Spodoptera littoralis. J Chem Ecol 32: 2333–2349 [PubMed]
  • Burow M, Rice M, Hause B, Gershenzon J, Wittstock U (2007) Cell- and tissue-specific localization and regulation of the epithiospecifier protein in Arabidopsis thaliana. Plant Mol Biol 64: 173–185 [PubMed]
  • Capella AN, Menossi M, Arruda P, Benedetti CE (2001) COI1 affects myrosinase activity and controls the expression of two flower-specific myrosinase-binding protein homologues in Arabidopsis. Planta 213: 691–699 [PubMed]
  • Dombrecht B, Xue GP, Sprague SJ, Kirkegaard JA, Ross JJ, Reid JB, Fitt GP, Sewelam N, Schenk PM, Manners JM, et al (2007) MYC2 differentially modulates diverse jasmonate-dependent functions in Arabidopsis. Plant Cell 19: 2225–2245 [PubMed]
  • Fahey JW, Zalcmann AT, Talalay P (2001) The chemical diversity and distribution of glucosinolates and isothiocyanates among plants. Phytochemistry 56: 5–51 [PubMed]
  • Foo HL, Gronning LM, Goodenough L, Bones AM, Danielsen B-E, Whiting DA, Rossiter JT (2000) Purification and characterisation of epithiospecifier protein from Brassica napus: enzymic intramolecular sulphur addition within alkenyl thiohydroximates derived from alkenyl glucosinolate hydrolysis. FEBS Lett 468: 243–246 [PubMed]
  • Fraser HB, Hirsh AE, Steinmetz LM, Scharfe C, Feldman MW (2002) Evolutionary rate in the protein interaction network. Science 296: 750–752 [PubMed]
  • Giamoustaris A, Mithen R (1995) The effect of modifying the glucosinolate content of leaves of oilseed rape (Brassica napus ssp oleifera) on its interaction with specialist and generalist pests. Ann Appl Biol 126: 347–363.
  • Griffiths DW, Deighton N, Birch AN, Patrian B, Baur R, Stadler E (2001) Identification of glucosinolates on the leaf surface of plants from the Cruciferae and other closely related species. Phytochemistry 57: 693–700 [PubMed]
  • Hirai MY, Klein M, Fujikawa Y, Yano M, Goodenowe DB, Yamazaki Y, Kanaya S, Nakamura Y, Kitayama M, Suzuki H, et al (2005) Elucidation of gene-to-gene and metabolite-to-gene networks in Arabidopsis by integration of metabolomics and transcriptomics. J Biol Chem 280: 25590–25595 [PubMed]
  • Hughes AL (1994) The evolution of functionally novel proteins after gene duplication. Proc Biol Sci 256: 119–124 [PubMed]
  • Jander G, Cui J, Nhan B, Pierce NE, Ausubel FM (2001) The TASTY locus on chromosome 1 of Arabidopsis affects feeding of the insect herbivore Trichoplusia ni. Plant Physiol 126: 890–898 [PubMed]
  • Kim JH, Jander G (2007) Myzus persicae (green peach aphid) feeding on Arabidopsis induces the formation of a deterrent indole glucosinolate. Plant J 49: 1008–1019 [PubMed]
  • Kliebenstein D, Lambrix V, Reichelt M, Gershenzon J, Mitchell-Olds T (2001a) Gene duplication and the diversification of secondary metabolism: side chain modification of glucosinolates in Arabidopsis thaliana. Plant Cell 13: 681–693 [PubMed]
  • Kliebenstein DJ, Figuth A, Mitchell-Olds T (2002) Genetic architecture of plastic methyl jasmonate responses in Arabidopsis thaliana. Genetics 161: 1685–1696 [PubMed]
  • Kliebenstein DJ, Gershenzon J, Mitchell-Olds T (2001b) Comparative quantitative trait loci mapping of aliphatic, indolic and benzylic glucosinolate production in Arabidopsis thaliana leaves and seeds. Genetics 159: 359–370 [PubMed]
  • Kliebenstein DJ, Kroymann J, Brown P, Figuth A, Pedersen D, Gershenzon J, Mitchell-Olds T (2001c) Genetic control of natural variation in Arabidopsis thaliana glucosinolate accumulation. Plant Physiol 126: 811–825 [PubMed]
  • Lambrix V, Reichelt M, Mitchell-Olds T, Kliebenstein D, Gershenzon J (2001) The Arabidopsis epithiospecifier protein promotes the hydrolysis of glucosinolates to nitriles and influences Trichoplusia ni herbivory. Plant Cell 13: 2793–2807 [PubMed]
  • Lankau RA (2007) Specialist and generalist herbivores exert opposing selection on a chemical defense. New Phytol 175: 176–184 [PubMed]
  • Leach J, Stevenson H, Rainbow A, Mullen L (1999) Effects of high plant populations on the growth and yield of winter oilseed rape (Brassica napus). J Agric Sci 132: 173–180.
  • Lister C, Dean C (1993) Recombinant inbred lines for mapping RFLP and phenotypic markers in Arabidopsis thaliana. Plant J 4: 745–750.
  • Louda S, Mole S (1991) Glucosinolates: chemistry and ecology. In GA Rosenthal, MR Baerenbaum, eds, Herbivores: Their Interactions with Secondary Plant Metabolites, Ed 2. Academic Press, San Diego, pp 123–164.
  • Loudet O, Chaillou S, Camilleri C, Bouchez D, Daniel-Vedele F (2002) Bay-0 × Shahdara recombinant inbred line population: a powerful tool for the genetic dissection of complex traits in Arabidopsis. Theor Appl Genet 104: 1173–1184 [PubMed]
  • Lynch M, Force A (2000) The probability of duplicate gene preservation by subfunctionalization. Genetics 154: 459–473 [PubMed]
  • Maruyama-Nakashita A, Nakamura Y, Tohge T, Saito K, Takahashi H (2006) Arabidopsis SLIM1 is a central transcriptional regulator of plant sulfur response and metabolism. Plant Cell 18: 3235–3251 [PubMed]
  • Mewis I, Tokuhisa JG, Schultz JC, Appel HM, Ulrichs C, Gershenzon J (2006) Gene expression and glucosinolate accumulation in Arabidopsis thaliana in response to generalist and specialist herbivores of different feeding guilds and the role of defense signaling pathways. Phytochemistry 67: 2450–2462 [PubMed]
  • Miles CI, del Campo ML, Renwick JA (2005) Behavioral and chemosensory responses to a host recognition cue by larvae of Pieris rapae. J Comp Physiol A Neuroethol Sens Neural Behav Physiol 191: 147–155 [PubMed]
  • Moyes CL, Raybould AF (2001) The role of spatial scale and intraspecific variation in secondary chemistry in host-plant location by Ceutorhynchus assimilis (Coleoptera: Curculionidae). Proc Biol Sci 268: 1567–1573 [PubMed]
  • Nastruzzi C, Cortesi R, Esposito E, Menegatti E, Leoni O, Iori R, Palmieri S (1996) In vitro cytotoxic activity of some glucosinolate-derived products generated by myrosinase hydrolysis. J Agric Food Chem 44: 1014–1021.
  • Nastruzzi C, Cortesi R, Esposito E, Menegatti E, Leoni O, Iori R, Palmieri S (2000) In vitro antiproliferative activity of isothiocyanates and nitriles generated by myrosinase-mediated hydrolysis of glucosinolates from seeds of cruciferous vegetables. J Agric Food Chem 48: 3572–3575 [PubMed]
  • Petersen BL, Chen S, Hansen CH, Olsen CE, Halkier BA (2002) Composition and content of glucosinolates in developing Arabidopsis thaliana. Planta 214: 562–571 [PubMed]
  • Reichelt M, Brown PD, Schneider B, Oldham NJ, Stauber E, Tokuhisa J, Kliebenstein DJ, Mitchell-Olds T, Gershenzon J (2002) Benzoic acid glucosinolate esters and other glucosinolates from Arabidopsis thaliana. Phytochemistry 59: 663–671 [PubMed]
  • Reifenrath K, Muller C (2007) Species-specific and leaf-age dependent effects of ultraviolet radiation on two Brassicaceae. Phytochemistry 68: 875–885 [PubMed]
  • Reymond P, Bodenhausen N, Van Poecke RMP, Krishnamurthy V, Dicke M, Farmer EE (2004) A conserved transcript pattern in response to a specialist and a generalist herbivore. Plant Cell 16: 3132–3147 [PubMed]
  • Spencer GF, Daxenbichler ME (1980) Gas chromatography-mass spectrometry of nitriles, isothiocyanates and oxazolidinethiones derived from cruciferous glucosinolates. J Sci Food Agric 31: 359–367.
  • Tookey HL (1973) Crambe thioglucoside glucohydrolase (EC 3.2.3.1): separation of a protein required for epithiobutane formation. Can J Biochem 51: 1654–1660 [PubMed]
  • Wentzell AM, Rowe HC, Hansen BG, Ticconi C, Halkier BA, Kliebenstein DJ (2007) Linking metabolic QTLs with network and cis-eQTLs controlling biosynthetic pathways. PLoS Genet 3: 1687–1701 [PubMed]
  • Wink M (1988) Plant breeding: importance of plant secondary metabolites for protection against pathogens and herbivores. Theor Appl Genet 75: 225–233.
  • Wink M (2003) Evolution of secondary metabolites from an ecological and molecular phylogenetic perspective. Phytochemistry 64: 3–19 [PubMed]
  • Wittstock U, Halkier BA (2002) Glucosinolate research in the Arabidopsis era. Trends Plant Sci 7: 263–270 [PubMed]
  • Wolda H (1988) Insect seasonality: why? Annu Rev Ecol Syst 19: 1–18.
  • Wolfinger R, O'Connell M (1993) Generalized linear mixed models: a pseudo-likelihood approach. J Statist Comput Simul 48: 233–243.
  • Zabala MD, Grant M, Bones AM, Bennett R, Lim YS, Kissen R, Rossiter JT (2005) Characterisation of recombinant epithiospecifier protein and its over-expression in Arabidopsis thaliana. Phytochemistry 66: 859–867 [PubMed]
  • Zhang J, Pontoppidan B, Xue J, Rask L, Meijer J (2002) The third myrosinase gene TGG3 in Arabidopsis thaliana is a pseudogene specifically expressed in stamen and petal. Physiol Plant 115: 25–34 [PubMed]
  • Zhang Z-Y, Ober JA, Kliebenstein DJ (2006) The gene controlling the quantitative trait locus EPITHIOSPECIFIER MODIFIER1 alters glucosinolate hydrolysis and insect resistance in Arabidopsis. Plant Cell 18: 1524–1536 [PubMed]