Metabolic engineering in Nicotiana benthamiana reveals key enzyme functions in Arabidopsis indole glucosinolate modification

Indole glucosinolates, derived from the amino acid Trp, are plant secondary metabolites that mediate numerous biological interactions between cruciferous plants and their natural enemies, such as herbivorous insects, pathogens, and other pests. While the genes and enzymes involved in the Arabidopsis thaliana core biosynthetic pathway, leading to indol-3-yl-methyl glucosinolate (I3M), have been identified and characterized, the genes and gene products responsible for modification reactions of the indole ring are largely unknown. Here, we combine the analysis of Arabidopsis mutant lines with a bioengineering approach to clarify which genes are involved in the remaining biosynthetic steps in indole glucosinolate modification. We engineered the indole glucosinolate biosynthesis pathway into Nicotiana benthamiana, showing that it is possible to produce indole glucosinolates in a noncruciferous plant. Building upon this setup, we demonstrate that all members of a small gene subfamily of cytochrome P450 monooxygenases, CYP81Fs, are capable of carrying out hydroxylation reactions of the glucosinolate indole ring, leading from I3M to 4-hydroxy-indol-3-yl-methyl and/or 1-hydroxy-indol-3-yl-methyl glucosinolate intermediates, and that these hydroxy intermediates are converted to 4-methoxy-indol-3-yl-methyl and 1-methoxy-indol-3-yl-methyl glucosinolates by either of two family 2 O-methyltransferases, termed indole glucosinolate methyltransferase 1 (IGMT1) and IGMT2.     

De novo indol-3-ylmethyl glucosinolate biosynthesis,and not long-distance transport,contributes to defence of Arabidopsis against powdery mildew

Powdery mildew is a fungal disease that affects a wide range of plants and reduces crop yield worldwide. As obligate biotrophs, powdery mildew fungi manipulate living host cells to suppress defence responses and to obtain nutrients. Members of the plant order Brassicales produce indole glucosinolates that effectively protect them from attack by non-adapted fungi. Indol-3-ylmethyl glucosinolate is constitutively produced in the phloem and transported to epidermal cells for storage. Upon attack, indol-3-ylmethyl glucosinolate is activated by CYP81F2 to provide broad-spectrum defence against fungi. How de novo biosynthesis and transport contribute to defence of powdery mildew-attacked epidermal cells is unknown. Bioassays and glucosinolate analysis demonstrate that GTR glucosinolate transporters are not involved in antifungal defence. Using quantitative live-cell imaging of fluorophore-tagged markers, we show that accumulation of the glucosinolate biosynthetic enzymes CYP83B1 and SUR1 is induced in epidermal cells attacked by the non-adapted barley powdery mildew Blumeria graminis f.sp. hordei. By contrast, glucosinolate biosynthesis is attenuated during interaction with the virulent powdery mildew Golovinomyces orontii. Interestingly, SUR1 induction is delayed during the Golovinomyces orontii interaction. We conclude that epidermal de novo synthesis of indol-3-ylmethyl glucosinolate contributes to CYP81F2-mediated broad-spectrum antifungal resistance and that adapted powdery mildews may target this process.     

植物激素与芥子油苷在生物合成上的相互作用

植物激素在植物的生长发育中起着关键性作用,芥子油苷是一类重要的次生代谢物质。植物激素与芥子油苷之间存在复杂的相互作用。生长素与吲哚类芥子油苷在生物合成上存在着相互作用。植物防卫信号分子与芥子油苷之间也存在相互作用,茉莉酸强烈诱导吲哚类芥子油苷生物合成相关基因CYP7982和CYP7983的表达,从而诱导吲哚-3-甲基芥子油苷和N-甲氧吲哚-3-甲基芥子油苷等吲哚类芥子油苷的生成,水杨酸和乙烯则能轻度诱导4-甲氧吲哚-3-甲基芥子油苷的生成。植物防卫信号转导途径相互作用以精细调节不同种类吲哚类芥子油苷的生成。     

Glucosinolate profiling of Brassica rapa cultivars after infection by Leptosphaeria maculans and Fusarium oxysporum

The glucosinolate contents of two different cultivars of Brassica rapa (Herfstraap and Oleifera) infected with Leptosphaeria maculans and Fusarium oxysporum were determined. Infection triggered the accumulation of aliphatic glucosinolates (gluconapin, progoitrin, glucobrassicanapin and gluconapoleiferin) and indole glucosinolate (4-hydroxy-glucobrassicin) in Herfstraap and of two indole glucosinolates (glucobrassicin and 4-hydroxy-glucobrassicin) in Oleifera. While total and aliphatic glucosinolates decreased significantly in Oleifera, a large increase was observed in Herfstraap after fungal infection. The indole glucosinolate glucobrassicin accumulated in Oleifera at a higher rate than Herfstraap especially after infection with F. oxysporum. Apparently the interaction between fungus and B. rapa is cultivar and fungal species specific.     

Indole Glucosinolate Biosynthesis Limits Phenylpropanoid Accumulation in Arabidopsis thaliana

Plants produce an array of metabolites (including lignin monomers and soluble UV-protective metabolites) from phenylalanine through the phenylpropanoid biosynthetic pathway. A subset of plants, including many related to Arabidopsis thaliana, synthesizes glucosinolates, nitrogen- and sulfur-containing secondary metabolites that serve as components of a plant defense system that deters herbivores and pathogens. Here, we report that the Arabidopsis thaliana reduced epidermal fluorescence5 (ref5-1) mutant, identified in a screen for plants with defects in soluble phenylpropanoid accumulation, has a missense mutation in CYP83B1 and displays defects in glucosinolate biosynthesis and in phenylpropanoid accumulation. CYP79B2 and CYP79B3 are responsible for the production of the CYP83B1 substrate indole-3-acetaldoxime (IAOx), and we found that the phenylpropanoid content of cyp79b2 cyp79b3 and ref5-1 cyp79b2 cyp79b3 plants is increased compared with the wild type. These data suggest that levels of IAOx or a subsequent metabolite negatively influence phenylpropanoid accumulation in ref5 and more importantly that this crosstalk is relevant in the wild type. Additional biochemical and genetic evidence indicates that this inhibition impacts the early steps of the phenylpropanoid biosynthetic pathway and restoration of phenylpropanoid accumulation in a ref5-1 med5a/b triple mutant suggests that the function of the Mediator complex is required for the crosstalk.     

Integration of Biosynthesis and Long-Distance Transport Establish Organ-Specific Glucosinolate Profiles in Vegetative Arabidopsis

Although it is essential for plant survival to synthesize and transport defense compounds, little is known about the coordination of these processes. Here, we investigate the above- and belowground source-sink relationship of the defense compounds glucosinolates in vegetative Arabidopsis thaliana. In vivo feeding experiments demonstrate that the glucosinolate transporters1 and 2 (GTR1 and GTR2), which are essential for accumulation of glucosinolates in seeds, are likely to also be involved in bidirectional distribution of glucosinolates between the roots and rosettes, indicating phloem and xylem as their transport pathways. Grafting of wild-type, biosynthetic, and transport mutants show that both the rosette and roots are able to synthesize aliphatic and indole glucosinolates. While rosettes constitute the major source and storage site for short-chained aliphatic glucosinolates, long-chained aliphatic glucosinolates are synthesized both in roots and rosettes with roots as the major storage site. Our grafting experiments thus indicate that in vegetative Arabidopsis, GTR1 and GTR2 are involved in bidirectional long-distance transport of aliphatic but not indole glucosinolates. Our data further suggest that the distinct rosette and root glucosinolate profiles in Arabidopsis are shaped by long-distance transport and spatially separated biosynthesis, suggesting that integration of these processes is critical for plant fitness in complex natural environments.     

Indole glucosinolate breakdown and its biological effects

Most species in the Brassicaceae produce one or more indole glucosinolates. In addition to the parent indol-3-ylmethylglucosinolate (IMG), other commonly encountered indole glucosinolates are 1-methoxyIMG, 4-hydroxyIMG, and 4-methoxyIMG. Upon tissue disruption, enzymatic hydrolysis of IMG produces an unstable aglucone, which reacts rapidly to form indole-3-acetonitrile and indol-3-ylmethyl isothiocyanate. The isothiocyanate, in turn, can react with water, ascorbate, glutathione, amino acids, and other plant metabolites to produce a variety of physiologically active indole compounds. Myrosinase-initiated breakdown of the substituted indole glucosinolates proceeds in a similar manner to that of IMG. Induction of indole glucosinolate production in response to biotic stress, experiments with mutant plants, and artificial diet assays suggest a significant role for indole glucosinolates in plant defense. However, some crucifer-feeding specialist herbivores recognize indole glucosinolates and their breakdown products as oviposition and/or feeding stimulants. In mammalian diets, IMG can have both beneficial and deleterious effects. Most IMG breakdown products induce the synthesis of phase 1 detoxifying enzymes, which may in some cases prevent carcinogenesis, but in other cases promote carcinogenesis. Recent advances in indole glucosinolate research have been fueled by their occurrence in the well-studied model plant Arabidopsis thaliana. Knowledge gained from genetic and biochemical experiments with A. thaliana can be applied to gain new insight into the ecological and nutritional properties of indole glucosinolates in other plant species.     

Glucosinolate changes in developing pods of single and double low varieties of autumn-sown oilseed rape (B. napus)

Single and double low varieties of oilseed rape were grown in the 1987/88 and 1988/89 seasons to study changes in the concentrations of total and individual glucosinolates within pods during development. Total glucosinolate concentration in seeds of all varieties increased during development when expressed on a fresh weight basis. The levels of the major alkenyl glucosinolates present in the seed; 2–hydroxy-3–butenyl, 3–butenyl and 4–pentenyl had been reduced in the transition from single to double low varieties. The major indole glucosinolates in the seed, 4–hydroxy-3–indolylmethyl and 3–indolylmethyl were present in the same amounts in single and double low varieties but in the latter represented a greater proportion of the total seed glucosinolate content. A decline in the total glucosinolate concentration in the pod walls with time together with the analogous profile of individual glucosinolates in the seeds and pod walls suggests that the pod wall is a major site of seed glucosinolate synthesis. Other plant parts may also have an important role to play in provision of intact glucosinolates or precursors to the pod walls for glucosinolate biosynthesis.     

Cytochrome P450 CYP79B2 from Arabidopsis catalyzes the conversion of tryptophan to indole-3-acetaldoxime, a precursor of indole glucosinolates and indole-3-acetic acid

Glucosinolates are natural plant products known as flavor compounds, cancer-preventing agents, and biopesticides. We report cloning and characterization of the cytochrome P450 CYP79B2 from Arabidopsis. Heterologous expression of CYP79B2 in Escherichia coli shows that CYP79B2 catalyzes the conversion of tryptophan to indole-3-acetaldoxime. Recombinant CYP79B2 has a K(m) of 21 microm and a V(max) of 7.78 nmol/h/ml culture. Inhibitor studies show that CYP79B2 is different from a previously described enzyme activity that converts tryptophan to indole-3-acetaldoxime (Ludwig-Müller, J. , and Hilgenberg, W. (1990) Phytochemistry, 29, 1397-1400). CYP79B2 is wound-inducible and expressed in leaves, stem, flowers, and roots, with the highest expression in roots. Arabidopsis overexpressing CYP79B2 has increased levels of indole glucosinolates, which strongly indicates that CYP79B2 is involved in indole glucosinolate biosynthesis. Our data show that oxime production by CYP79s is not restricted to those amino acids that are precursors for cyanogenic glucosides. Our data are consistent with the hypothesis that indole glucosinolates have evolved from cyanogenesis. Indole-3-acetaldoxime is a precursor of the plant hormone indole-3-acetic acid, which suggests that CYP79B2 might function in biosynthesis of indole-3-acetic acid. Identification of CYP79B2 provides an important tool for modification of the indole glucosinolate content to improve nutritional value and pest resistance.     

Gene Duplication in the Diversification of Secondary Metabolism: Tandem 2-Oxoglutarate–Dependent Dioxygenases Control Glucosinolate Biosynthesis in Arabidopsis

Secondary metabolites are a diverse set of plant compounds believed to have numerous functions in plant–environment interactions. The large chemical diversity of secondary metabolites undoubtedly arises from an equally diverse set of enzymes responsible for their biosynthesis. However, little is known about the evolution of enzymes involved in secondary metabolism. We are studying the biosynthesis of glucosinolates, a large group of secondary metabolites, in Arabidopsis to investigate the evolution of enzymes involved in secondary metabolism. Arabidopsis contains natural variations in the presence of methylsulfinylalkyl, alkenyl, and hydroxyalkyl glucosinolates. In this article, we report the identification of genes encoding two 2-oxoglutarate–dependent dioxygenases that are responsible for this variation. These genes, AOP2 and AOP3, which map to the same position on chromosome IV, result from an apparent gene duplication and control the conversion of methylsulfinylalkyl glucosinolate to either the alkenyl or the hydroxyalkyl form. By heterologous expression in Escherichia and the correlation of gene expression patterns to the glucosinolate phenotype, we show that AOP2 catalyzes the conversion of methylsulfinylalkyl glucosinolates to alkenyl glucosinolates. Conversely, AOP3 directs the formation of hydroxyalkyl glucosinolates from methylsulfinylalkyl glucosinolates. No ecotype coexpressed both genes. Furthermore, the absence of functional AOP2 and AOP3 leads to the accumulation of the precursor methylsulfinylalkyl glucosinolates. A third member of this gene family, AOP1, is present in at least two forms and found in all ecotypes examined. However, its catalytic role is still uncertain.     

Indole Glucosinolates and Camalexin do not Influence the Development of the Clubroot Disease in Arabidopsis thaliana

Root galls of Brassicaceae caused by Plasmodiophora brassicae are dependent on increased auxin and cytokinin formation. In this study we investigated whether indole glucosinolates are involved in indole‐3‐acetic acid (IAA) biosynthesis in root galls, by using a genetic approach. The cytochrome P450 enzymes, CYP79B2 and CYP79B3, convert tryptophan to indole‐3‐acetaldoxime (IAOx), which is a precursor for indole glucosinolates and the phytoalexin camalexin in Arabidopsis thaliana. Root galls of the Arabidopsis ecotypes Wassilewskija (WS) and Columbia (Col) accumulated camalexin, WS at levels up to 320 μg/g dry weight. By contrast, camalexin was absent in root galls of cyp79b2/b3 double mutants. Infection rate and disease index as a measure of club development in mutant and wild‐type plants of the two ecotypes were investigated and no differences were found in gall formation. This demonstrates that camalexin is an ineffective inhibitor of P. brassicae and indole glucosinolates are not the source of elevated levels of IAA in galls, because free IAA levels in mutant galls were comparable with those in wild type.     

Contribution of Glucosinolate Transport to Arabidopsis Defense Responses

Accumulation of glucosinolates, a class of defense-related secondary metabolites found almost exclusively in the Capparales, is induced in response to a variety of biological stresses. It is often assumed that elevated glucosinolate levels result from de novo biosynthesis, but glucosinolate transport from other parts of the plant to the site of herbivory or pathogen infection can also contribute to the defense response. Several studies with Arabidopsis and other crucifers have demonstrated that glucosinolates from vegetative tissue are transported to developing seeds. Here we discuss evidence that long-chain aliphatic glucosinolates are transported to the site of herbivory in response to Myzus persicae (green peach aphid) feeding on Arabidopsis.Key Words: glucosinolate, transport, graft, Arabidopsis, Myzus persicae, aphid     

QTL Analysis Using SNP Markers Developed by Next-Generation Sequencing for Identification of Candidate Genes Controlling 4-Methylthio-3-Butenyl Glucosinolate Contents in Roots of Radish,Raphanus sativus L

SNP markers for QTL analysis of 4-MTB-GSL contents in radish roots were developed by determining nucleotide sequences of bulked PCR products using a next-generation sequencer. DNA fragments were amplified from two radish lines by multiplex PCR with six primer pairs, and those amplified by 2,880 primer pairs were mixed and sequenced. By assembling sequence data, 1,953 SNPs in 750 DNA fragments, 437 of which have been previously mapped in a linkage map, were identified. A linkage map of nine linkage groups was constructed with 188 markers, and five QTLs were detected in two F₂ populations, three of them accounting for more than 50% of the total phenotypic variance being repeatedly detected. In the identified QTL regions, nine SNP markers were newly produced. By synteny analysis of the QTLs regions with Arabidopsis thaliana and Brassica rapa genome sequences, three candidate genes were selected, i.e., RsMAM3 for production of aliphatic glucosinolates linked to GSL-QTL-4, RsIPMDH1 for leucine biosynthesis showing strong co-expression with glucosinolate biosynthesis genes linked to GSL-QTL-2, and RsBCAT4 for branched-chain amino acid aminotransferase linked to GSL-QTL-1. Nucleotide sequences and expression of these genes suggested their possible function in 4MTB-GSL biosynthesis in radish roots.     

CYP79F1 and CYP79F2 have distinct functions in the biosynthesis of aliphatic glucosinolates in Arabidopsis

Cytochromes P450 of the CYP79 family catalyze the conversion of amino acids to oximes in the biosynthesis of glucosinolates, a group of natural plant products known to be involved in plant defense and as a source of flavor compounds, cancer-preventing agents and bioherbicides. We report a detailed biochemical analysis of the substrate specificity and kinetics of CYP79F1 and CYP79F2, two cytochromes P450 involved in the biosynthesis of aliphatic glucosinolates in Arabidopsis thaliana. Using recombinant CYP79F1 and CYP79F2 expressed in Escherichia coli and Saccharomyces cerevisiae, respectively, we show that CYP79F1 metabolizes mono- to hexahomomethionine, resulting in both short- and long-chain aliphatic glucosinolates. In contrast, CYP79F2 exclusively metabolizes long-chain elongated penta- and hexahomomethionines. CYP79F1 and CYP79F2 are spatially and developmentally regulated, with different gene expression patterns. CYP79F2 is highly expressed in hypocotyl and roots, whereas CYP79F1 is strongly expressed in cotyledons, rosette leaves, stems, and siliques. A transposon-tagged CYP79F1 knockout mutant completely lacks short-chain aliphatic glucosinolates, but has an increased level of long-chain aliphatic glucosinolates, especially in leaves and seeds. The level of long-chain aliphatic glucosinolates in a transposon-tagged CYP79F2 knockout mutant is substantially reduced, whereas the level of short-chain aliphatic glucosinolates is not affected. Biochemical characterization of CYP79F1 and CYP79F2, and gene expression analysis, combined with glucosinolate profiling of knockout mutants demonstrate the functional role of these enzymes. This provides valuable insights into the metabolic network leading to the biosynthesis of aliphatic glucosinolates, and into metabolic engineering of altered aliphatic glucosinolate profiles to improve nutritional value and pest resistance.     

High glucosinolate broccoli: a delivery system for sulforaphane

The development of hybrid broccoli genotypes with enhanced levels of 4-methylsulphinylbutyl glucosinolate, the precursor of anticarcinogenic isothiocyanate sulforaphane (SF), by introgressing genomic segments from the wild ancestor Brassica villosa is described. We demonstrate that to obtain enhanced levels of either 3-methylsulphinylpropyl or 4-methylsulphinylbutyl glucosinolate it is necessary to have B. villosa alleles in either a homozygous or heterozygous state at a single quantitative trait locus (QTL) on O2. The ratio of these two glucosinolates, and thus whether iberin or SF is generated upon hydrolysis, is determined by the presence or absence of B. villosa alleles at this QTL, but also at an additional QTL2 on O5. We further demonstrate that following mild cooking high glucosinolate broccoli lines generate about three fold higher levels of SF than conventional varieties. Commercial freezing processes and storage of high glucosinolate broccoli maintains the high level of glucosinolates compared to standard cultivars, although the blanching process denatures the endogenous myrosinase activity.     

Quantitative trait loci analysis of non-enzymatic glucosinolate degradation rates in Brassica oleracea during food processing

Epidemiological and mechanistic studies show health-promoting effects of glucosinolates and their breakdown products. In literature, differences in non-enzymatic glucosinolate degradation rates during food processing between different vegetables are described, which provide the basis for studying the genetic effects of this trait and breeding vegetables with high glucosinolate retention during food processing. Non-enzymatic glucosinolate degradation, induced by heat, was studied in a publicly available Brassica oleracea doubled haploid population. Data were modeled to obtain degradation rate constants that were used as phenotypic traits to perform quantitative trait loci (QTL) mapping. Glucosinolate degradation rate constants were determined for five aliphatic and two indolic glucosinolates. Degradation rates were independent of the initial glucosinolate concentration. Two QTL were identified for the degradation rate of the indolic glucobrassicin and one QTL for the degradation of the aliphatic glucoraphanin, which co-localized with one of the QTL for glucobrassicin. Factors within the plant matrix might influence the degradation of different glucosinolates in different genotypes. In addition to genotypic effects, we demonstrated that growing conditions influenced glucosinolate degradation as well. The study identified QTL for glucosinolate degradation, giving the opportunity to breed vegetables with a high retention of glucosinolates during food processing, although the underlying mechanisms remain unknown.     

Modulation of CYP79 genes and glucosinolate profiles in Arabidopsis by defense signaling pathways

Glucosinolates are natural plant products that function in the defense toward herbivores and pathogens. Plant defense is regulated by multiple signal transduction pathways in which salicylic acid (SA), jasmonic acid, and ethylene function as signaling molecules. Glucosinolate content was analyzed in Arabidopsis wild-type plants in response to single or combinatorial treatments with methyljasmonate (MeJA), 2,6-dichloro-isonicotinic acid, ethylene, and 2,4-dichloro-phenoxyacetic acid, or by wounding. In addition, several signal transduction mutants and the SA-depleted transgenic NahG line were analyzed. In parallel, expression of glucosinolate biosynthetic genes of the CYP79 gene family and the UDPG:thiohydroximate glucosyltransferase was monitored. After MeJA treatment, the amount of indole glucosinolates increased 3- to 4-fold, and the corresponding Trp-metabolizing genes CYP79B2 and CYP79B3 were both highly induced. Specifically, the indole glucosinolate N-methoxy-indol-3-ylmethylglucosinolate accumulated 10-fold in response to MeJA treatment, whereas 4-methoxy-indol-3-ylmethylglucosinolate accumulated 1.5-fold in response to 2,6-dichloro-isonicotinic acid. In general, few changes were seen for the levels of aliphatic glucosinolates, although increases in the levels of 8-methylthiooctyl glucosinolate and 8-methylsulfinyloctyl glucosinolate were observed, particularly after MeJA treatments. The findings were supported by the composition of glucosinolates in the coronatine-insensitive mutant coi1, the ctr1 mutant displaying constitutive triple response, and the SA-overproducing mpk4 and cpr1 mutants. The present data indicate that different indole glucosinolate methoxylating enzymes are induced by the jasmonate and the SA signal transduction pathways, whereas the aliphatic glucosinolates appear to be primarily genetically and not environmentally controlled. Thus, different defense pathways activate subsets of biosynthetic enzymes, leading to the accumulation of specific glucosinolates.     

Elucidating the Role of Transport Processes in Leaf Glucosinolate Distribution

In Arabidopsis (Arabidopsis thaliana), a strategy to defend its leaves against herbivores is to accumulate glucosinolates along the midrib and at the margin. Although it is generally assumed that glucosinolates are synthesized along the vasculature in an Arabidopsis leaf, thereby suggesting that the margin accumulation is established through transport, little is known about these transport processes. Here, we show through leaf apoplastic fluid analysis and glucosinolate feeding experiments that two glucosinolate transporters, GTR1 and GTR2, essential for long-distance transport of glucosinolates in Arabidopsis, also play key roles in glucosinolate allocation within a mature leaf by effectively importing apoplastically localized glucosinolates into appropriate cells. Detection of glucosinolates in root xylem sap unambiguously shows that this transport route is involved in root-to-shoot glucosinolate allocation. Detailed leaf dissections show that in the absence of GTR1 and GTR2 transport activity, glucosinolates accumulate predominantly in leaf margins and leaf tips. Furthermore, we show that glucosinolates accumulate in the leaf abaxial epidermis in a GTR-independent manner. Based on our results, we propose a model for how glucosinolates accumulate in the leaf margin and epidermis, which includes symplasmic movement through plasmodesmata, coupled with the activity of putative vacuolar glucosinolate importers in these peripheral cell layers.Feeding behavior of herbivorous insects and distribution of defense compounds in plants have been suggested to be a result of an arms race between plants and insects that has spanned millions of years (Ehrlich and Raven, 1964). Whether insects adapted first to plants or the other way around is an ongoing debate in this research field (Schoonhoven et al., 2005; Ali and Agrawal, 2012). Leaf margin accumulation of defense compounds has been demonstrated in various plant species (Gutterman and Chauser-Volfson, 2000; Chauser-Volfson et al., 2002; Kester et al., 2002; Cooney et al., 2012). In the model plant Arabidopsis (Arabidopsis thaliana), higher concentration of glucosinolates, which constitute a major part of the chemical defense system in this plant (Kliebenstein et al., 2001a; Halkier and Gershenzon, 2006), was found at the leaf midrib and margins compared with the leaf lamina (Shroff et al., 2008; Sønderby et al., 2010). This nonuniform leaf distribution of glucosinolates appeared to explain the feeding pattern of a generalist herbivore (Helicoverpa armigera), as it avoided feeding at the leaf margin and midrib (Shroff et al., 2008). A similar feeding pattern on Arabidopsis was observed for a different generalist herbivore, Spodoptera littoralis (Schweizer et al., 2013). Interestingly, S. littoralis was shown to favor feeding from Arabidopsis leaf margins in glucosinolate-deficient mutants (Schweizer et al., 2013), which could indicate an inherent preference for margin feeding and that Arabidopsis adapted to such behavior by accumulating defense compounds here. A damaged leaf margin may be more critical for leaf stability than damage to inner leaf parts (Shroff et al., 2008), further motivating protection of this tissue. The margin-feeding preference of S. littoralis might be explained by better nutritional value of the leaf margin cells (Schweizer et al., 2013), which has been shown to consist of specialized elongated cell files (Koroleva et al., 2010; Nakata and Okada, 2013).Other distribution patterns have been reported for glucosinolates in an Arabidopsis leaf. A study investigating spatiotemporal metabolic shifts during senescence in Arabidopsis reported that fully expanded mature leaves exhibited a glucosinolate gradient from base to tip, with highest level of glucosinolates at the leaf base (Watanabe et al., 2013). In contrast to the horizontal plane, less has been reported on distribution of glucosinolates in the vertical plane of a leaf. A localization study of cyanogenic glucosides, defense molecules related to glucosinolates (Halkier and Gershenzon, 2006), determined that these compounds primarily were located in the epidermis of sorghum (Sorghum bicolor; Kojima et al., 1979). Whereas epidermis-derived trichomes in Arabidopsis were recently demonstrated to contain glucosinolates and to express glucosinolate biosynthetic genes (Frerigmann et al., 2012), no studies have investigated glucosinolates in the epidermal cell layer.Based on promoter-GUS studies, biosynthesis of glucosinolates in leaves of Arabidopsis has been associated with primarily major and minor veins in leaves and silique walls (Mikkelsen et al., 2000; Reintanz et al., 2001; Tantikanjana et al., 2001; Chen et al., 2003; Grubb et al., 2004; Schuster et al., 2006; Gigolashvili et al., 2007; Li et al., 2011; Redovniković et al., 2012). The discrepancy between vasculature-associated glucosinolate biosynthesis and margin accumulation of glucosinolates suggests that transport processes must be involved in establishing the distribution pattern of glucosinolates within a leaf.Plant transport systems include the apoplastic xylem, the symplastic phloem, and plasmodesmata. Xylem transport is mainly driven by an upward pull generated by transpiration from aerial plant organs, thereby directing transport to sites of rapid evaporation (such as leaf margins; Sattelmacher, 2001). Phloem flow is facilitated by an osmosis-regulated hydrostatic pressure difference between source and sink tissue, primarily generated by Suc bulk flow (Lucas et al., 2013). Plasmodesmata are intercellular channels that establish symplasmic pathways between neighboring cells, and most cell types in a plant are symplastically connected via plasmodesmata (Roberts and Oparka, 2003). Translocation of small molecules in these channels is driven by diffusion and is regulated developmentally as well as spatially to form symplastically connected domains (Roberts and Oparka, 2003; Christensen et al., 2009). To what extent any of these transport processes are involved in establishing specific distribution patterns of glucosinolates within leaves is not known.Recently, two plasma membrane-localized, glucosinolate-specific importers, GLUCOSINOLATE TRANSPORTER1 (GTR1) and GTR2, were identified in Arabidopsis (Nour-Eldin et al., 2012). In leaf, their expression patterns were shown to be in leaf veins (GTR1 and GTR2) and surrounding mesophyll cells (GTR1; Nour-Eldin et al., 2012). Absence of aliphatic and indole glucosinolates in seeds of the gtr1gtr2 double knockout (dKO) mutant (gtr1gtr2 dKO) demonstrated that these transporters are essential for long-distance glucosinolate transport to the seeds and indicates a role in phloem loading (Nour-Eldin et al., 2012). Another study investigating long-distance transport of glucosinolates in the 3-week-old wild type and gtr1gtr2 dKO indicated that GTR1 and GTR2 were involved in bidirectional transport of aliphatic glucosinolates between root and shoot via both phloem and xylem pathways (Andersen et al., 2013). The authors suggested a role for GTR1 and GTR2 in the retention of long-chained aliphatic glucosinolates in roots by removing the compounds from the xylem (Andersen et al., 2013).Identification of the glucosinolate transporters GTR1 and GTR2 has provided a molecular tool to investigate the role of transport processes in establishing leaf glucosinolate distribution. In this study, we have performed a detailed spatial investigation of the distribution of an exogenously fed glucosinolate (sinigrin) and endogenous glucosinolates within mature wild-type and gtr1gtr2 dKO Arabidopsis leaves, achieved by collecting and analyzing leaf parts at the horizontal (y axis: petiole, base, and tip; x axis: midrib, lamina, and margin) as well as at the vertical leaf plane (z axis: abaxial epidermis). Furthermore, we analyze wild-type and gtr1gtr2 dKO root xylem sap and leaf apoplastic fluids for glucosinolates. Based on our results, we propose a model where GTR1 and GTR2 import glucosinolates from the apoplast to the symplast and where the glucosinolate distribution pattern within an Arabidopsis leaf is established via symplasmic movement of glucosinolates through plasmodesmata, coupled with the activity of putative vacuolar glucosinolate importers in peripheral cell layers.