|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 280, Issue 23, 21867-21881, June 10, 2005
Detoxification and Transcriptome Response in Arabidopsis Seedlings Exposed to the Allelochemical Benzoxazolin-2(3H)-one*
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Despite the ecological and agronomic importance of this class of natural products, relatively little is known concerning the molecular target sites mediating their toxicity or the adaptive strategies mounted by plants in defense against these compounds. Furthermore, in comparison to other areas of chemical ecology, molecular and genomics based approaches have seldom been employed in the field of allelopathy. One recent exception is the use of DNA microarrays to characterize the transcriptome response of Arabidopsis seedlings exposed to ()-catechin, an allelochemical identified in root exudates of Centaurea maculosa (6). Benzoxazinoids and their benzoxazolinone derivatives represent one of the more intensively studied classes of allelochemicals (7, 8). For example, the genes encoding all of the enzymes required for the biosynthesis of the benzoxazinoid 2,4-dihydroxy-7-methoxy-2H-1,4-benzoxazin-3(4H)-one have been identified in corn via transposon tagging, and they represent the first known example of a plant secondary metabolic pathway organized as a gene cluster (9). The biosynthesis of benzoxazinoids, particularly in young seedlings, is generally associated with cereals such as corn, rye, and wheat but have also been identified in species of Acanthaceae, Ranunculaceae, and Scrophulariaceae (7). Benzoxazinoids and benzoxazolinones act as defense compounds against microbial pathogens as well as insect herbivores, and within the rhizosphere they play an important role as allelochemicals (7, 8). Formation of the benzoxazolinone, benzoxazolin-2(3H)-one (BOA),1 results from a two-step degradation of the glucoside of DIBOA (2,4-dihydroxy-2H-1,4-benzoxazin-3(4H)-one; see Fig. 1). Numerous plant species exhibit tolerance to benzoxazinoids, such as BOA, and can rapidly metabolize them to less phytotoxic glucoside and glucoside carbamate derivatives (Fig. 1), potentially due to having co-evolved in association with allelopathic species within the same communities (4).
A primary response to the presence of xenobiotic compounds in both prokaryotic and eukaryotic organisms involves the induction of detoxifying enzymes and transporters, which facilitate the inactivation and elimination of toxins, and the associated metabolic processes can be divided into discrete phases (1012). In phase I, compounds are typically modified such that a functional group such as a hydroxyl moiety is added or exposed through the action of hydrolases, cytochrome P450s, or peroxidases. The availability of functional groups then facilitates the formation of glucosyl, glutathione, and malonyl conjugates in phase II, through the action of specific glucosyltransferases, glutathione S-transferases, and less frequently, malonyltransferases. Conjugated forms of xenobiotics can then be recognized by specific membrane-associated transporters such as ABC transporters (1315) in phase III detoxification, resulting in their vacuolar sequestration or release into the apoplasmic space via exocytosis. In addition, subsequent less well characterized processes can also occur, resulting in further enzymatic modification in vacuoles or deposition of metabolites as cell wall-associated bound residues (11), sometimes referred to as "phase IV" detoxification.
The signaling pathways associated with the response to xenobiotic stress in mammals have been extensively studied and involve both xenobiotic ligand-activated transcription factors and redox sensing proteins (16, 17). In plant systems, it is also well established that the expression of specific genes involved in processes related to the detoxification of allelochemicals and synthetic herbicides can be induced in the presence of these compounds (18); however, evidence for the existence of integrated, coordinately regulated chemical detoxification gene networks, such as those identified in animal systems, is scarce. Furthermore, the signaling components involved in plant responses to xenobiotic stress have not been identified at present. An important class of agrichemicals, the herbicide safeners, may act as potent inducers of these signaling pathways, thereby rendering crops less susceptible to herbicide-induced injury (19, 20).
Arabidopsis thaliana represents an excellent model for the study of plant responses to allelochemicals and other environmental toxins, as its genome has been fully sequenced and well characterized. The commercial availability of nearly full-genome oligonucleotide microarrays (21) opens the possibility for identifying specific enzyme isoforms within complex gene families, such as cytochrome P450s, glutathione S-transferases, and UDP-glucosyltransferases (2224), representing potential components of a coordinated xenobiotic defense response network. In addition to gaining basic insight into plant xenobiotic stress responses, molecular ecotoxicological studies in plants are also of considerable importance to human health and agriculture, given the incidents of soil and groundwater contamination occurring worldwide from industrial pollutants (11, 25).
In the present work we have focused on the analysis of the physiological response, detoxification pathways, and transcriptome responses in Arabidopsis seedlings exposed to the model benzoxazolinone allelochemical BOA. Feeding studies using exogenously supplied BOA revealed that detoxification in Arabidopsis occurs predominantly through O-glucosylation of the intermediate BOA-6-OH, most likely requiring the sequential action of as-yet-unidentified cytochrome P450 and UDP glucosyltransferase activities. Transcriptional profiling experiments using microarrays representing
24,000 transcripts identified a significant number of genes potentially involved in phase I, II, and III detoxification processes that are induced following exposure to this allelochemical. By using a subset of these genes, we further demonstrate their induction in response to a variety of structurally diverse xenobiotic compounds, suggesting they comprise components of a general xenobiotic response network. These data significantly expand upon previous studies examining plant transcriptional responses to allelochemicals and other environmental toxins and provide a foundation for elucidating both the enzymes and regulatory mechanisms involved.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
The extracts were analyzed for detoxification products by high pressure liquid chromatography using a model 126 chromatograph (Beckman Instruments, Fullerton, CA) equipped with a diode array detector (model 168) and an Ultrasphere ODS RP 18 column. Compounds were eluted with the following gradients: 1 min, 100% eluent A (0.1% trifluoroacetic acid in H2O); 121 min, 20% eluent B (methanol) linear; 2141 min, 80% eluent B linear; 4143 min, 100% eluent B linear, using a flow rate of 1 ml/min. The detection wavelengths used were 280 and 405 nm. Major metabolites were identified by co-chromatography with natural BOA-6-O-glucoside and with synthetic BOA-6-OH and glucoside carbamate prepared as described previously (4, 26). These compounds were also used as external standards for quantification.
BOA Treatments for Growth Inhibition and Microarray StudiesFor all BOA growth inhibition and microarray experiments, aseptically germinated A. thaliana (Col-0) seedlings were maintained in a growth chamber at 21 °C under a 16-h photoperiod and light intensity of 150 µmol m2 s1. Seeds were surface-sterilized in 70% ethanol for 5 min and then rinsed two times in sterile distilled water, followed by treatment with 0.5x bleach (3% sodium hypochlorite) and 0.05% Tween 20 for 10 min, and then finally rinsed four times in sterile distilled water. Following surface sterilization, seeds were placed on top of an
2.0-cm-high stack of 9.0-cm Whatman No. 4 filter disks and allowed to air-dry in a sterile hood. Liquid and semi-solid germination media used for all experiments contained 0.5 x Murashige and Skoog salts, 1x Gamborg's B5 vitamins, and 1.0% sucrose (w/v), adjusted to pH 5.7 with KOH. Semi-solid media also contained 1.0% (w/v) agar.
For root elongation assays, seeds were placed in a horizontal line
2.0 cm from the edge in 9.0 x 9.0 x 1.5-cm square Petri dishes containing semi-solid germination media supplemented with different concentrations of BOA. Control plates contained solvent alone (0.1% ethanol). Plates were cold-treated for 3 days, transferred to a growth chamber, and then maintained in a vertical position for 10 days, at which time root lengths were scored. All BOA and control treatments were performed in duplicate; each replicate consisted of 25 seedlings.
For BOA treatments prior to microarray analyses,
200 seeds were scooped into a microspatula and then distributed evenly over the surface of a sterile 0.3-µm microporous membrane raft supported by a buoyant float (Osmotek Ltd., Rehovat, Israel). Seeds, rafts, and floats were placed on liquid germination media in Lifeguard tissue culture vessels with 4.0-cm vented lids (Osmotek Ltd., Rehovat, Israel), cold-treated for 3 days, and then transferred to a growth chamber. After 10 days, BOA (or 0.1% ethanol for control treatments) was added to the media and then the vessels were returned to the growth chambers until harvest. At the end of the treatment period, seedlings were removed from the vessels, flash-frozen in liquid nitrogen, and then stored at 80 °C prior to analysis.
Chemical Treatments for Real Time PCR AssaysFollow-up chemical treatments for quantitative real time PCR experiments were performed as described above for microarray experiments. Ten-day-old seedlings grown on floating microporous membrane/raft assemblies were exposed to two different concentrations of fenclorim (Toronto Research Chemicals, Inc., Ontario, Canada), benoxacor (Sigma), 2,4-dichlorophenoxyacetic acid (Sigma), phenobarbital (Sigma), p-hydroxybenzoic acid (sodium salt; Sigma), and 4-dimethylaminoantipyrine (free base; Sigma). For 2,4,5-trichlorophenol (Sigma), due to extensive injury observed on seedlings exposed to 10 mM treatments, only 100 µM treatments were used. Stock solutions were prepared in either Me2SO or ethanol and then added to seedlings growing on liquid media as described above in duplicate treatments. Control treatments (0.25% ethanol and 0.5% Me2SO) were also performed in duplicate. At 24 h post-treatment, seedlings were flash-frozen in liquid nitrogen and stored at 80 °C prior to analysis.
RNA IsolationTotal RNAs for use in microarray experiments were isolated from 0.5 g of flash-frozen, pulverized 10-day-old seedling tissues using the Trizol reagent (Invitrogen), with an additional homogenization step of 30 s at 25,000 rpm using a hand held homogenizer. The RNA recovered was then re-purified with an RNeasy plant mini-kit (Qiagen, Valencia, CA) per the manufacturer's instructions. RNA recovery and purity were determined spectrophotometrically, and sample integrity was assessed by agarose gel electrophoresis.
Total RNAs for real time PCR experiments shown in Fig. 5 were isolated from 50 mg of flash-frozen, pulverized 10-day-old seedling tissues using an RNeasy plant mini-kit, with an additional homogenization step of 30 s at 25,000 rpm as described above. The RNA samples were also treated with DNase I "on column" using an RNase-free DNase kit as per the manufacturer's instructions (Qiagen, Inc., Valencia, CA) to remove residual DNA contamination. RNA recovery and purity were also determined spectrophotometrically for these samples, and sample integrity was also assessed by agarose gel electrophoresis.
Microarray Hybridization and AnalysisMicroarray hybridizations were performed for three independent replicates with Affymetrix Arabidopsis ATH1 Genome Arrays, using protocols described by Affymetrix, Inc. (Santa Clara, CA). GeneTraffic software (Iobion Informatics, La Jolla, CA) was used to conduct a two-class comparison analysis on normalized and log-transformed signal values obtained from Affymetrix Microarray Suite software version 5.0. A significance analysis test (unpaired t test with Benjamini-Hochberg false discovery rate correction) was performed to test the equality of the mean signal values between the two classes. Means for each class were then inverse-transformed to provide a geometrical mean as an overall estimate of expression. In this manner, a more robust estimate of overall expression, less impacted by outliers or skewed expression levels (as compared with a simple arithmetic average of the raw signal values from each array), was obtained. Fold change was then calculated as the simple ratio of overall signal values from the two classes. Genes with p values of
0.05 were considered to be significantly differentially expressed. Genes that were common to both the I50 and I80 data were identified. A floor adjustment of 64 was applied to genes with very low signal values to avoid artifactually large fold change calculations. Genes that were induced or repressed by at least 2.5-fold in either the I50 or I80 data were identified and retained for further analysis. Differentially expressed genes, thus identified, were then annotated using the NetAffx data base provided by Affymetrix (www.affymetrix.com). In some cases the annotations of unknown genes were further refined by performing additional BLAST searches or updated as additional literature reports became available.
Quantitative Real Time RT-PCR AssaysFirst strand cDNAs were synthesized from 2 µg of total RNA in a 100-µl reaction volume using the TaqMan reverse transcription reagents kit (Applied Biosystems, Foster City, CA) as per the manufacturer's instructions. Quantitative real time PCRs were performed in triplicate using the GenAmp® 5700 sequence detection system (Applied Biosystems). Independent PCRs were performed using the same cDNA for both the gene of interest and 18 S rRNA, using the SYBR® Green PCR Master Mix (Applied Biosystems). Gene-specific primers were designed for the gene of interest and 18 S rRNA using Primer Express® software (Applied Biosystems) and the Amplify program (27). Closely related sequences within the Arabidopsis (Col-0) genome were identified via BLASTN queries of the AGI transcripts data base using the BLAST server at The Arabidopsis Information Resource (www.arabidopsis.org/Blast/). All sequences thus identified were then aligned using the ClustalW alignment function of MegAlign software (DNAstar, Inc. Madison, WI). Gene-specific primer pairs were then manually selected such that at least one primer per pair contained a minimum of two consecutive mismatches at the 3' end when compared against all related Arabidopsis transcripts. For almost all genes analyzed, both primers within a pair fulfilled the above criteria and contained numerous additional mismatches when compared against related sequences. The PCR conditions consisted of denaturation at 95 °C for 10 min, followed by 40 cycles of denaturation at 95 °C for 15 s and annealing/extension at 60 °C for 1 min. A dissociation curve was generated at the end of each PCR cycle to verify that a single product was amplified using software provided with the GeneAmp® 5700 sequence detection system. A negative control reaction in the absence of template (no template control) was also routinely performed in triplicate for each primer pair. The change in fluorescence of SYBR® Green I dye in every cycle was monitored by the GenAmp® 5700 system software, and the threshold cycle (CT) above background for each reaction was calculated. The CT value of 18 S rRNA was subtracted from that of the gene of interest to obtain a
CT value. The CT value of an arbitrary calibrator (e.g. untreated sample in the case of up-regulated genes) was subtracted from the
CT value to obtain a 
CT value. The fold changes in expression level relative to the calibrator were expressed as 2
CT. For determination of statistical significance, pairwise comparisons were performed between treated and control sample 
CT values using an independent two-tailed t test, assuming common variance. Differences associated with p values
0.05 were considered significant.
Motif Searches5' upstream sequences from 1500 to 1 (relative to potential transcription start sites if available) were retrieved using RSA tools (see Ref. 28; rsat.ulb.ac.be/rsat/), by selecting the "mRNA" feature type, and preventing overlap with upstream open reading frames. The program Motif Sampler (see Ref. 29; www.esat.kuleuven.ac.be/~thijs/Work/MotifSampler.html) was then used for the identification of over-represented motifs within the retrieved sequences. All searches were performed using a precompiled 3rd order Markov background model based on Arabidopsis upstream sequences (29), prior probability of finding 1 motif instance = 0.5, maximum number of motif instances per sequence = 0 (no limit), and maximum allowed overlap between different motifs = 2. Each data set was analyzed 10 times using the same parameters to reduce local optima (30), and only cases where an identical or similar consensus sequence was returned from multiple runs were further considered. In addition, two independent statistical tests were employed for motif validation. For the first, searches were performed on 28,577 Arabidopsis 1500 to 1 (see above) 5' upstream sequences to determine the total number of motif instances for all predicted Arabidopsis genes using the RSA-tools "genome-scale DNA pattern" search function (28). These values were used to estimate the probability (p) of occurrence for a given motif within the total number of nucleotides searched on both DNA strands, and then p values for each experimental result were determined based on calculated binomial probabilities using the "PROBBNML" function within the SAS version 9.1 statistical analysis software package (SAS Institute, Inc., Cary, NC). The second test involved a bootstrap analysis using the POBO program (31). For each motif, 1000 bootstrap pseudoclusters were generated by random sampling with replacement from the input promoter data set and then compared with 1000 pseudoclusters generated in a similar manner from a background data set consisting of all predicted Arabidopsis 5' upstream regions. The number of sequences in each pseudo-cluster generated was equal to the number in the input data set, and all background data set sequences were 1500 bp in length. An independent t test was then performed using pseudo-cluster values to estimate the probability for the number of motifs observed in the input data set to occur by chance. A p value
0.05 was considered significant for both of these tests.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
-glucoside 1,2-carbamate) shown in Fig. 1 (26). An additional metabolite, gentiobioside carbamate (1-(2-hydroxyphenylamino)-1-deoxy-
-gentiobioside 1,2-carbamate) derived from the glucoside carbamate, has also been identified in Zea mays but was undetectable in several other Poaceae species analyzed (26, 32). Sensitivity to BOA is highly species-dependent, although dicot species are in general less tolerant than monocots (33). This difference could be explained, at least in part, by the observation that dicots tend to detoxify BOA via O-glucosylation, whereas in several monocot species examined a significant proportion of exogenously applied BOA is metabolized via N-glucosylation, and for Z. mays, virtually all metabolism occurs via N-glucosylation (4). The glucoside carbamate exhibits dramatically reduced phytotoxicity relative to BOA-6-O-glucoside in bioassays, whereas the O-glucoside precursor, BOA-6-OH, is approximately twice as phytotoxic as the parent compound BOA (4). The pathway leading to the formation of the glucoside carbamate via N-glucosylation of BOA therefore represents a more efficient mechanism for detoxification of the allelochemical and would presumably confer increased tolerance relative to species that favor accumulation of O-glucoside.
|
11:1 for both the 250 and 500 µM treatments. Taken together, these data clearly indicate that BOA detoxification in Arabidopsis occurs mainly via the O-glucosylation route, likely requiring the sequential action of as-yet-unidentified cytochrome P450 and UDP glucosyltransferase activities. The surprisingly high levels of BOA-6-OH, which accumulate in Arabidopsis as compared with other species examined (4), suggest either an inefficient conversion of BOA-6-OH to BOA-6-O-glucoside or alternatively, competing, higher levels of
-glucosidase, which markedly increase steady-state levels of BOA-6-OH. Although these data strongly suggest that BOA detoxification in Arabidopsis is a relatively inefficient process, vacuolar sequestration and/or apoplastic extrusion mechanisms must also be taken into account when considering whole plant tolerance levels against allelochemicals and other xenobiotic compounds.
|
540 µM and 1mM, respectively (Fig. 2B). Similar dose-dependent responses for root elongation have been reported for other species grown in the presence of BOA (33), indicating the suitability of A. thaliana as a model species for examining plant responses to this allelochemical.
Transcriptome responses were next analyzed using commercially prepared oligonucleotide arrays (Affymetrix, Inc., Santa Clara, CA) representing
81% of the 29,454 predicted Arabidopsis genes (34). Approximately 200 seedlings per treatment were germinated aseptically on membrane rafts floated on germination media, stratified for 3 days, and then maintained for 10 days in a growth chamber under a 16-h photoperiod (see "Experimental Procedures"). After 10 days, BOA was added to the media at a final concentration of 540 µM or 1.0 mM, and seedlings were harvested 24 h later. Mock treatments were also performed, which consisted of ethanol alone (final 0.1% v/v), followed by a 24-h incubation. Three experimental replicates were prepared for all treatments and controls and grown in separate culture vessels to account for experimental variation, and then total RNAs were isolated from each replicate and used for hybridizations with Affymetrix ATH1 arrays.
|
0.05), which includes a 2.5-fold cut-off for change in expression observed in at least one of the two treatments. One hundred fifty eight genes were induced, and 30 were repressed in both the I50 and I80 treatments, representing in total
0.8% of all Arabidopsis genes included on the ATH1 array. To confirm further the reliability of the observed changes, quantitative real time PCR assays were developed for 16 differentially expressed genes identified by microarray analysis, and transcript levels were measured using the same RNA samples used for the microarray hybridizations. The gene-specific primer pairs designed for these assays are shown in Table II. As shown in Fig. 3, all of the selected genes indicated as being differentially expressed by transcriptional profiling were confirmed by quantitative real time PCR. Remarkably, even differences between expression levels in the I50 and I80 treatments relative to controls were confirmed by these analyses. The magnitude of the fold changes determined for most of the genes tested were similar for the two methods; however, for four of the genes (At4g12490, At1g15520, At4g12500, and At2g15490) the increases were substantially larger when analyzed by real time PCR. Relative fold changes seen in microarray data are often larger when analyzed by real time PCR, predominantly due to limitations in the dynamic range and sensitivity of microarray experiments and global normalization methods applied to microarray data sets. In addition, cross-hybridization may occur in microarray experiments with probe sets designed for closely related genes, resulting in complex signals that will not occur in real time PCR assays using appropriately designed gene-specific primers. In the case of a differentially expressed putative ABC-type transporter At1g15520 (Fig. 3), the discrepancy is likely attributable to the fact that the Affymetrix software assigned an "absent" call to the corresponding probe set in untreated samples, thus a true measure of fold induction could not be determined.
|
27.9 and 26.2% of all functions assigned, respectively (Fig. 4). The metabolism category included genes potentially involved in chemical detoxification processes; thus there is significant overlap with genes assigned to cell rescue, defense, and virulence, as each gene could be assigned to multiple categories. Also of significant interest are putative transcription factors (4.1%) and genes associated with cellular communication and signal transduction mechanisms (7.8%), which could potentially mediate responses for other genes identified in the data set. A striking number of the most highly up-regulated genes represent enzymes potentially involved in the modification, conjugation, or elimination of xenobiotics (supplemental Table I). Included in this category would be various quinone oxidoreductases, aldo-ketoreductases, peroxidases, cytochrome P450-like proteins, glutathione S-transferases, glucosyltransferases, and representatives of several families of multidrug efflux transporters, including the ABC, major facilitator superfamily (MFS), and multidrug and toxic compound extrusion (MATE) transporters. A list of selected up-regulated genes potentially associated with detoxification processes, as well as putative signal transduction components that could potentially be involved in regulating their expression, is shown in Table III.
|
Induction of Putative Phase I-related ProteinsEnzymes associated with phase I-related detoxification processes primarily perform oxidative modifications, frequently resulting in decreased toxicity and increased solubility of parent molecules (11). A significant body of evidence suggests that cytochrome P450s perform the majority of oxidative phase I reactions in plants, and this may in part explain the complexity of P450-related protein families that exist in most plant species. For example, a recent survey of the Arabidopsis genome suggested that as many as 286 genes comprise a superfamily of P450-related sequences (22). Cytochrome P450s are hemeproteins located in the endoplasmic reticulum and catalyze a wide array of oxidative reactions using diverse substrates. Thus, in addition to their general metabolic role, the large families of plant P450 proteins potentially represent an arsenal of degradative activities for the detoxification of foreign and endogenous compounds. Results from the microarray analysis revealed six cytochrome P450 sequences significantly up-regulated in response to I50 and I80 BOA treatments (Table III), whose physiological roles and substrate specificities are not known at present.
Less well studied in plant species, but well known to participate in xenobiotic detoxification pathways in animals, are members of the aldo-ketoreductase (AKR) superfamily (reviewed in Refs. 37 and 38). The Arabidopsis genome encodes at least 18 genes related to AKRs (www.arabidopsis.org) whose functions remain unknown. AKRs have been identified in all kingdoms and metabolize structurally diverse substrates, including monosaccharides, steroids, aromatic aldehydes and hydrocarbons, and aliphatic aldehydes, utilizing nicotinamide cofactors.s Unlike other detoxifying protein families, relatively little is known in plants or animals regarding the transcriptional responses for AKR genes in tissues exposed to toxins or other inducers. Four Arabidopsis AKR-like sequences were significantly induced by I50 and I80 BOA treatments (Table III). Most interestingly, an AKR-like sequence, In2-2, was isolated from corn roots by differential screening of safener-treated tissues (39). Northern analyses revealed that In2-2 expression increased from undetectable levels within 30 min of exposure to the substituted benzenesulfonamide safener 2-N-(aminocarbonyl)-2-chlorobenzenesulfonamide and that this induction was not observed following various stress treatments. Subsequent experiments using In2-2 promoter:reporter gene fusions in transformed tobacco also demonstrated inducibility by the herbicide chlorsulfuron and exogenously supplied branched-chain amino acids, both known to inhibit the enzyme acetolactate synthase (40). Further experiments will be necessary to definitively assign a role for AKR enzymes in plant xenobiotic detoxification, as has been demonstrated for animals (37); however, the transcriptional responses observed in Arabidopsis (see below) for AKR-like sequences to safeners, herbicides, and other xenobiotics suggest they may in fact play such a role.
|
|
Induction of Putative Phase II and Phase III-related ProteinsA significant portion of genes transcriptionally induced by BOA represent protein families such as glutathione S-transferases (GSTs), UDP-glucosyltransferases (UGTs), quinone oxidoreductases, and various multidrug efflux transporters potentially involved in phase II and phase III detoxification processes (Table III). The induction of a GST at the level of transcription in response to herbicide safener exposure was first reported by Wiegand et al. (47), and numerous reports have subsequently described the regulation of GSTs in response to a variety of agrichemicals and other xenobiotics (reviewed in Ref. 48). GSTs play a central role in phase II reactions, catalyzing the transfer of glutathione to a substrate to form S-glutathione conjugates, which are then substrates for ABC transporters for removal (13). Not unexpectedly, the seven putative GSTs shown to be up-regulated in the present study include some of the most highly induced transcripts in BOA-treated Arabidopsis seedlings (Table III).
Plant UGTs play an equally well established role in phase II metabolism and catalyze the transfer of glucose to diverse xenobiotic substrates (18), producing both N- and O-glucoside derivatives of the parent compounds. A recent survey of the Arabidopsis genome revealed the presence of 107 glucosyltransferase genes (24), and significant progress has been made toward elucidating the biochemical activities of this complex family (24, 49). In contrast to glutathionylated compounds, less is known regarding the transport of glucosides, and it is possible that both ABC transporters and tonoplast-associated H+ antiporters are involved (50).
In BOA-treated Arabidopsis seedlings, UGTs were the most highly represented family of up-regulated genes, with 11 different members being significantly induced (Table III). Most striking is the almost 70-fold induction observed by real time PCR in the I80 treatments for UGT73B4 (At2g15490; Fig. 3), which was below the limit of detection for microarray analysis in untreated samples (i.e. was assigned an absent call; data not shown). The induction of UGT enzymatic activities in plants in response to xenobiotic compounds has been documented previously (51); however, increases in steady-state levels for the corresponding mRNAs have not to our knowledge been directly demonstrated prior to this study. Given that the glycoside of BOA is the predominant metabolite formed in plants, including Arabidopsis (32), the possibility cannot be excluded that plant UGT gene promoter activities are specifically induced following exposure to BOA and related compounds.
Fifteen membrane transporter-like sequences were induced by BOA, including six putative ABC-type transporters, an MFS-type transporter (At5g13750), and a MATE-type transporter (At2g04040; Table III). ABC-, MFS-, and MATE-type transporters represent different multidrug efflux protein superfamilies associated with resistance to antibiotics and other toxins in microorganisms (52). In contrast to MFS and MATE-type transporters, induction of ABC-type transporters at the level of steady-state mRNA accumulation in response to xenobiotic compounds has been reported previously for plants (5356). In studies performed by van den Brûle and Smart (55), the BOA-induced ABC transporter AtPDR12 (At1g15520; Table III) was shown to be highly induced by cycloheximide and the antifungal diterpene sclareol, but not by hormone, cadmium, and cold or salt stress treatments. In addition, a recent study performed by Campbell et al. (15) showed induction of AtPDR12 transcript levels in response to pathogen infection, salicylic acid, ethylene, and methyl jasmonate. By using various signal transduction mutants, it was also demonstrated that induction of AtPDR12 by pathogen infection requires salicylic acid, NPR1 function, and sensitivity to jasmonic acid and ethylene. A T-DNA insertion line for AtPDR12 was more sensitive to sclareol than wild-type seedlings in germination assays, thus collectively these data suggested a potential role for AtPDR12 in transporting antifungal compounds accumulating in response to pathogen infection (15). The data generated in the present study suggest that this transporter may also have a more general role in the elimination of phytotoxins unrelated to pathogen defense (see below).
|
Response of Selected Genes to Xenobiotics and Inducers of Xenobiotic MetabolismThe genome-wide responses revealed by transcriptional profiling of BOA-treated Arabidopsis seedlings suggests that exposure to this allelochemical induces the expression of a battery of defense genes involved in the detoxification of xenobiotic compounds (Table III). Presumably, the signal transduction pathway(s) involved in this response would not be highly specific, given the structural diversity of environmental and endogenous toxins typically encountered during the life cycle of a plant (62). Alternatively, the expression levels for these genes could be coordinately increased in response to a wide variety of chemical agents through the action of broad specificity xenobiotic-sensing mechanisms, as has been observed in other organisms (6365). A more detailed examination of the response of select BOA-responsive genes following exposure to a variety of xenobiotics, and proposed inducers of plant chemical detoxification pathways, could therefore provide additional insights concerning the nature of plant xenobiotic defense responses and potentially identify key target genes and regulatory factors involved as well.
To explore this further, Arabidopsis seedlings were exposed to different concentrations of the herbicide 2,4-D, p-hydroxybenzoic acid (an allelochemical agent in wheat and potentially Sorghum bicolor; see Refs. 66 and 67), the environmental pollutant 2,4,5-trichlorophenol, the herbicide safeners fenclorim and benoxacor, and the drugs aminopyrine and phenobarbital, which have been shown to induce putative detoxification pathways involving ABC transporters, glutathione S-transferases, or P450s in both plants and animals (56, 68, 69). Chemical treatments were performed in an identical manner to the BOA treatments described above, and then plants were harvested at 24 h post-treatment, and RNA levels were monitored for 20 selected BOA-responsive genes by real time PCR (Table III). The results from these experiments are shown in Fig. 5.
Overall, steady-state mRNA levels were increased in response to multiple treatments for a majority of the genes analyzed, although the fold changes observed for a given treatment as well as the relative responsiveness to different treatments varied significantly from gene to gene. More importantly, a subset of the genes analyzed were strongly induced by a majority of the compounds tested, suggestive of a physiological role involving defense against toxins. For example, steady-state mRNA levels of ABC transporter AtPDR12 (At1g15520) exhibited greater than 5-fold increases following exposure to all of the xenobiotics tested, with the exception of aminopyrine, where only a modest (3.1-fold) increase was observed at the highest concentration (Fig. 5). Benoxacor and 2,4,5-trichlorophenol elicited surprisingly dramatic increases in AtPDR12 expression; transcript levels were induced more than 5,000-fold by 1 mM benoxacor treatments and more than 400-fold following exposure to 100 µM 2,4,5-trichlorophenol. Glucosyltransferase UGT73B4 (At2g15490) transcript levels were induced more than 30-fold by five of seven compounds tested and increased
17- and 14-fold at the highest concentrations tested for phenobarbital and aminopyrine, respectively. Dramatic responses were also observed for UGT73B4 following exposure to safeners and 2,4,5-trichlorophenol;
100-fold increases were observed in response to 100 µM fenclorim and 2,4,5-trichlorophenol treatments and
800-fold for 1 mM benoxacor treatments. Most interestingly, 100 µM p-hydroxybenzoate treatments induced UGT73B4 transcript levels
32-fold, and this compound has also been shown to be a substrate for this enzyme in addition to 3,4-dihydroxybenzoate and the hydroxycoumarin allelochemicals esculetin and scopoletin (24, 49). For both of these genes (AtPDR12 and UGT73B4), it is interesting to note that basal expression levels (i.e. expression levels in control-treated seedlings) were below the threshold of detection by microarray analysis (not shown) and yet were induced to high levels in response to the treatments performed. These observations would be consistent with physiological roles for AtPDR12 and UGT73B4 principally involving chemical defense, as opposed to more general housekeeping functions. Similarly, "marginal" or absent calls were assigned by the Affymetrix software for basal expression levels of quinone oxidoredutase At5g16980 which, as discussed below, was also highly induced by a majority of the treatments performed.
In addition to AtPDR12 and UGT73B4, steady-state mRNA levels for nine other genes analyzed were induced
5-fold or greater following exposure to at least five of the seven compounds tested (Fig. 5). These included glutathione S-transferase AtGSTU24 (At1g17170), cytochrome P450 CYP81D11 (At3g28740), an aldo-ketoreductase (At1g60730), a short chain type dehydrogenase/reductase (At3g04000), two different quinone oxidoreductases (At5g16980 and At5g16970), glucosyltransferase UGT73B5 (At2g15480), a malonyltransferase-related sequence (At5g39050), and At4g12490, a sequence of undetermined function related to Arabidopsis pEARLI 1 (70). As observed for AtPDR12 and UGT73B4, steady-state mRNA levels increased dramatically (more than 100-fold) in response to specific treatments for many of these genes, particularly those involving herbicide safeners (Fig. 5). For the class III peroxidase sequences tested, PER62 (At5g39580) and PER49 (At4g36430), neither exhibited a general pattern of chemical responsiveness nor were they significantly up-regulated following exposure to either herbicide safener.
All of the compounds tested in the present work caused some level of induction for the majority of the genes analyzed, although striking differences in their overall effectiveness were observed (Fig. 5). The safener benoxacor was clearly the most effective agent, with 1 mM treatments increasing steady-state mRNA levels 5-fold or greater for 18 of the 20 selected genes. The median level of induction for the 18 benoxacor-responsive genes was
21-fold, nearly three times that observed for the safener fenclorim (Fig. 5). Previous studies using Arabidopsis have demonstrated the ability of benoxacor to induce the expression of specific phi, theta, and tau class GST gene family members (71), and modest increases for several ABC transporter mRNA sequences have also been reported for benoxacor-treated Arabidopsis cell suspensions (54). The data shown in Fig. 5 significantly expand upon existing information concerning the molecular action of this important class of agrichemicals. For example, the data strongly suggest that, at least for benoxacor and fenclorim, herbicide-safening action involves the induction of a battery of genes potentially involving all aspects of xenobiotic metabolism and elimination, representing a wide range of enzymatic and transport-related functions. In addition, the observation that transcript levels for many of these genes also increase following exposure to numerous, structurally unrelated xenobiotic compounds further suggests that these safeners activate broad specificity xenobiotic defense pathways responsive to a wide array of compounds.
Recent studies in mammalian systems have revealed that transcriptional responses of genes involved in phase I, II, and III detoxification pathways are mediated through receptor-xenobiotic ligand interactions as well as redox-sensitive factors associated with the actin cytoskeleton (16, 17). For example, the aryl hydrocarbon nuclear receptor (AhR) stimulates the expression of a battery of detoxifying enzymes by binding promiscuously to structurally diverse xenobiotic compounds (64). Following binding to ligand, complexed AhR is translocated to the nucleus and subsequently dimerizes with the co-activator Arnt, and then the heterodimeric ligand-Ahr-Arnt complex stimulates transcription via xenobiotic-response elements in target gene promoters. A second important mechanism involves members of the steroid family of nuclear receptors, the pregnane X receptor, and the constitutive androstane receptor, which are activated by binding to xenobiotic ligands, and subsequently activate transcription of target promoters by binding as heterodimers with the 9-cis-retinoic acid receptor (63, 72). As is the case observed for the AhR, pregnane X receptor binds promiscuously to structurally diverse xenobiotic ligands, and thus serves as a broad specificity xenobiotic sensor that can induce the expression of a battery of detoxification-related genes (63).
Although ligand-activated transcription factors, such as AhR and pregnane X receptor, have no apparent plant counterparts, a third major mammalian xenobiotic-sensing system involving the basic leucine zipper transcription factor Nrf2 and its repressor Keap1 has been elucidated recently (17, 73). In the absence of chemical inducers, Nrf2 is localized cytoplasmically via interaction with the kelch domain-containing protein Keap1, which is itself tethered to the actin cytoskeleton. In the presence of a wide range of structurally diverse sulfhydryl-reactive electrophilic compounds, two critical highly reactive cysteine residues on the Keap1 protein (Cys-273 and Cys-288) form an intermolecular disulfide bridge between two Keap1 monomers, which then releases Nrf2 for transport into the nucleus. In association with other transcription factors, possibly small Maf proteins, Nrf2 subsequently activates target promoters containing the antioxidant-response element, which include a battery of phase II-detoxifying enzymes such as GSTs, glucosyltransferases, and quinone oxidoreductases (65). Thus the Keap1 sensor-Nrf2 system provides a relatively simple feed-forward mechanism for responding to a large diversity of xenobiotic electrophiles, aspects of which could readily be envisaged to occur in plant cells. For example, binding of the Arabidopsis transcription factor TGA1 with its co-activator NPR1 has been shown recently to be redox-regulated through critical cysteine residues within the protein (74). Most interestingly, mRNA levels for a kelch repeat-containing F-box protein family member with sequence similarity to Keap1 were increased by exposure to BOA in the present study (At2g44130; see supplemental Table I). Members of this plant-specific protein family (75) could conceivably fulfill some of the functions attributed to Keap1, including cytoskeletal anchoring via the kelch repeats, and facilitating the ubiquitin-mediated degradation of its interacting partners via the F-box motif (17, 73, 75). Of further interest is the observation that chemical structures for all known herbicide safeners contain at least one "soft" electrophilic center that could potentially modify protein sulfhydryl groups (76, 77). Given their potency as transcriptional inducers of plant phase I, II, and III detoxifying enzymes as demonstrated in the present study as well as by others, it is conceivable that the effects of these compounds are mediated via as-yet-unidentified protein redox sensors analogous to the mammalian Keap1 sensor-Nrf2 system.
Emerging evidence suggests that broad specificity chemosensory mechanisms such as those described above are commonplace among both prokaryotes and eukaryotes (17, 64). Their ability to coordinately up-regulate batteries of xenobiotic-detoxifying enzymes enables organisms to expeditiously metabolize and eliminate a wide array of xenobiotics and endotoxins, and this capacity is often further enhanced through the induction of enzymes exhibiting broad substrate specificities (63). The identification of specific Arabidopsis genes potentially representing components of a xenobiotic-detoxifying arsenal, as in the present study, provides the tools necessary to elucidate the signal transduction pathway(s) underlying the xenobiotic defense response in Arabidopsis via molecular genetic and bioinformatics based approaches. Potent elicitors for these pathways, such as the herbicide safener benoxacor, will undoubtedly serve as important pharmacological probes for these studies as well.
Analysis of 5' Upstream SequencesIt is possible that many or a subset of the genes responsive to BOA share common cis-regulatory elements in their promoter sequences, which could mediate either general transcriptional responses to stresses associated with xenobiotic exposure or, alternatively, responses specific to benzoxazinoids such as BOA. With the exception of the involvement of TGA transcription factors, basic/leucine zipper protein family members, which mediate responses to the auxenic herbicide 2,4-D via interaction with as-1-like elements present in some GST promoter sequences (78, 79), little is known regarding other transcription factors and cis-regulatory elements involved in responses to xenobiotic stress in plants.
To examine this question, the 5' upstream regions from genes identified by our microarray studies were analyzed for 8-,10-, and 12-mer statistically over-represented motifs using the program Motif Sampler (29), and then identified motifs were compared with known cis-regulatory elements in the PlantCARE data base (80). The 5' upstream sequences analyzed were grouped by those most highly induced in the microarray results, as well as by enzyme family (Table IV). As only the 8-mer searches yielded promising motifs observed in multiple runs of Motif Sampler (29), the 10- and 12-mer results were not considered further. More importantly, all of the motifs listed in Table IV were significantly over-represented in the data sets analyzed (p
0.0043), as compared with the number of motif occurrences observed in searches performed with 28,577 Arabidopsis predicted upstream regions.
|
As mentioned, TGA transcription factors have been implicated in the activation of some xenobiotic-responsive promoters; therefore, we were also interested in examining upstream regions in selected groups of promoters for the presence of as-1-like elements. Characteristic of these elements are two imperfect TGACGTCA palindrome motifs whose centers are separated by 12 bp (87). Upstream regions for all genes analyzed by Motif Sampler were also scanned for the presence of as-1-like elements using the pattern finding algorithm within RSA-Tools (28). Among the GST genes found to be up-regulated (Table III), two of the gene upstream sequences (AtGSTU7; At2g29420 and AtGSTU1; At2g29490), contain TGAC half-sites exhibiting the 12-bp spacing consistent with an as-1-like element (not shown). One of the six P450 gene upstream regions (CYP81D11; At3g28740; Table III), contained an as-1-like element, and corresponding mRNA levels increased following 2,4-D exposure (Fig. 5), suggesting this motif may be active. In addition, the glucosyltransferase UGT73B1 (At4g34138) identified by our microarray studies contains an as-1-like element in its upstream region and was also shown recently to be induced by 2,4-D and salicylic acid (88), suggesting a potential role for these elements in the activation of some Arabidopsis glucosyltransferase promoters during xenobiotic stress. Inspection of the upstream regions for the remaining 10 UGTs shown in Table III revealed only one additional UGT gene (UGT75B1; At1g05560) containing an as-1-like element on the "" strand (not shown); however, steady-state mRNA levels for this UGT were unresponsive to 10 µM 2,4-D treatments and exhibited only modest increases (
3-fold) in response to 100 µM treatments (Fig. 5), suggesting that this motif may possess only minimal activity. as-1-like elements were identified in the upstream regions for two additional genes selected for the experiments shown in Fig. 5, ABC transporter AtPDR12 and a putative FAD-linked oxidoreductase (At4g20860). AtPDR12 transcript levels increased
7- and 5-fold following exposure to 10 and 100 µM 2,4-D, respectively (Fig. 5), potentially due to the presence of the identified as-1-like element. Only modest increases (
2.5-fold) were observed for transcript levels for the putative oxidoreductase following exposure to 10 µM 2,4-D, suggesting minimal or no functionality for the as-1 element identified in its upstream region. The absence of obvious as-1-like elements in the majority of the upstream sequences identified in our results, coupled with their responsiveness to herbicide safeners (Fig. 5), indicates that other mechanisms play a more central role in their activation following exposure to xenobiotics.
An important caveat concerning data generated by promoter motif discovery algorithms is that the functionality of an over-represented motif must be determined experimentally and that motif detection in silico can only provide likely candidates for such follow-up studies (30). In addition, it is possible that motifs identified in the present study are functional but not associated with physiological processes related to phytotoxin exposure. Nevertheless, the list of statistically over-represented motifs shown in Table IV generated from these analyses provides a potentially useful framework for investigating transcriptional control mechanisms influencing a critical biological process.
In summary, in the present study we have identified O-glucosylation as the predominant detoxification mechanism in Arabidopsis seedlings exposed to the allelochemical BOA. Using nearly full-genome transcriptional profiling studies and quantitative real time RT-PCR, we have further identified specific genes that may represent integrated components of a diverse broad specificity chemical defense network. These data significantly expand upon previous studies examining plant transcriptional responses to environmental toxins and provide a foundation for further biochemical and genetic experiments to more fully elucidate plant xenobiotic detoxification pathways, and the chemosensory mechanisms critical to plant survival in the presence of allelochemicals and other environmental toxins.
| FOOTNOTES |
|---|
The on-line version of this article (available at www.jbc.org) contains Table I. ![]()
|| Present address: Dept. of Biochemistry, Max-Planck-Institute for Chemical Ecology, Hans-Knöll-Strasse 8, D-07745 Jena, Germany. ![]()

Present address: United States Department of Agriculture, Agricultural Research Service, Forage-Animal Production Research Unit, N220 Agricultural Sciences North, Lexington, KY 40546. ![]()
To whom correspondence should be addressed: United States Department of Agriculture, Agricultural Research Service, Natural Products Utilization Research Unit, P. O. Box 8048, University, MS 38677. Tel.: 662-915-7965; Fax: 662-915-1035; E-mail: sbaerson{at}olemiss.edu.
1 The abbreviations used are: BOA, benzoxazolin-2(3H)-one; DIBOA, 2,4-dihydroxy-2H-1,4-benzoxazin-3(4H)-one; ABC, ATP binding cassette; MATE, multidrug and toxic compound extrusion; MFS, major facilitator superfamily; AKR, aldo-keto reductase; GSTs, glutathione S-transferases; UGTs, UDP-glucosyltransferases; AhR, aryl hydrocarbon receptor; RT, reverse transcription; 2,4-D, 2,4-dichlorophenoxyacetic acid. ![]()
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
B. Fode, T. Siemsen, C. Thurow, R. Weigel, and C. Gatz The Arabidopsis GRAS Protein SCL14 Interacts with Class II TGA Transcription Factors and Is Essential for the Activation of Stress-Inducible Promoters PLANT CELL, November 1, 2008; 20(11): 3122 - 3135. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Golisz, M. Sugano, and Y. Fujii Microarray expression profiling of Arabidopsis thaliana L. in response to allelochemicals identified in buckwheat J. Exp. Bot., August 1, 2008; 59(11): 3099 - 3109. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Brazier-Hicks, K. M. Evans, O. D. Cunningham, D. R. W. Hodgson, P. G. Steel, and R. Edwards Catabolism of Glutathione Conjugates in Arabidopsis thaliana: ROLE IN METABOLIC REACTIVATION OF THE HERBICIDE SAFENER FENCLORIM J. Biol. Chem., July 25, 2008; 283(30): 21102 - 21112. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Mueller, B. Hilbert, K. Dueckershoff, T. Roitsch, M. Krischke, M. J. Mueller, and S. Berger General Detoxification and Stress Responses Are Mediated by Oxidized Lipids through TGA Transcription Factors in Arabidopsis PLANT CELL, March 1, 2008; 20(3): 768 - 785. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. R. Baerson, F. E. Dayan, A. M. Rimando, N. P. D. Nanayakkara, C.-J. Liu, J. Schroder, M. Fishbein, Z. Pan, I. A. Kagan, L. H. Pratt, et al. A Functional Genomics Investigation of Allelochemical Biosynthesis in Sorghum bicolor Root Hairs J. Biol. Chem., February 8, 2008; 283(6): 3231 - 3247. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Manabe, N. Tinker, A. Colville, and B. Miki CSR1, the Sole Target of Imidazolinone Herbicide in Arabidopsis thaliana Plant Cell Physiol., September 1, 2007; 48(9): 1340 - 1358. [Abstract] [Full Text] [PDF] |
||||
![]() |
Z. Pan, A. M. Rimando, S. R. Baerson, M. Fishbein, and S. O. Duke Functional Characterization of Desaturases Involved in the Formation of the Terminal Double Bond of an Unusual 16:3{Delta}9, 12, 15 Fatty Acid Isolated from Sorghum bicolor Root Hairs J. Biol. Chem., February 16, 2007; 282(7): 4326 - 4335. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. E. Sattler, L. Mene-Saffrane, E. E. Farmer, M. Krischke, M. J. Mueller, and D. DellaPenna Nonenzymatic Lipid Peroxidation Reprograms Gene Expression and Activates Defense Markers in Arabidopsis Tocopherol-Deficient Mutants PLANT CELL, December 1, 2006; 18(12): 3706 - 3720. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| All ASBMB Journals | Molecular and Cellular Proteomics |
| Journal of Lipid Research | ASBMB Today |