<title>Background</title> <p>Arsenic (As) is a toxic metalloid found ubiquitously in the environment [##UREF##0##1##] and is classified as a human carcinogen [##UREF##1##2##]. Currently, the US Environmental Protection Agency declares arsenic as the highest priority hazardous substance found at contaminated sites in the United States (see Availability and requirements section for URL). Naturally high levels of arsenic in drinking water have caused major human health problems in the United States, China, Argentina, Taiwan, and most notably in Bangladesh and India where tens of millions of people have been affected [##REF##12842773##3##,##REF##12165430##4##]. Arsenic is highly toxic at low concentrations, therefore drinking water safety standards were lowered from 50 to 10 μg/L in the U.S. [##UREF##2##5##].</p> <p>Plants typically encounter arsenic in the anionic forms of arsenate [As (V)] and arsenite [As (III)], both of which have different cytotoxic effects [##REF##12857839##6##]. As (III) reacts with the sulfhydryl groups of enzymes and proteins, thereby inhibiting cellular function and resulting in death [##UREF##3##7##]. Alternatively, As (V) is an analog of the macronutrient phosphate, so it competes with phosphate for uptake in the roots, as well as in the cytoplasm where it might disrupt metabolism by replacing phosphate in ATP to form unstable ADP-As [##UREF##4##8##]. Once taken up by the roots, arsenate is reduced to a more highly toxic species, arsenite, which is subsequently detoxified via soluble thiols such as glutathione and/or phytochelatins (PCs) and transported for vacuolar sequestration [##REF##10759512##9##]. PCs are low molecular weight thiolate peptides of the general structure (γ-Glu-Cys)<sub><italic>n</italic></sub>-Gly (<italic>n </italic>= 2–11) and are synthesized from glutathione by the constitutively present phytochelatin synthase [##REF##16594069##10##]. Both arsenate and arsenite efficiently induce the production of PCs in plants [##UREF##5##11##], however it is believed since arsenate has no affinity for the sulfhydryl groups in PCs, As (V) is reduced in the cytoplasm, resulting in As(III)-PC complexes [##REF##12857839##6##]. Glutathione and PCs have been reported to form As(III)-tris-thiolate complexes in <italic>Brassica juncea </italic>upon exposure to As(V) [##REF##10759512##9##]. Therefore, PC synthesis causes a depletion of cellular glutathione, resulting in a decreased capacity to quench reactive oxygen species (ROS) [##REF##11351093##12##].</p> <p>Phytoremediation has emerged as an alternative technology for removing toxic metals from contaminated soils and groundwater. The potential for phytoremediation to be an effective means of removing arsenic from contaminated sites has been demonstrated in hyperaccumulators of the <italic>Pteris </italic>genus [##REF##17305305##13##, ####REF##12371185##14##, ##REF##16297966##15####16297966##15##] and may be enhanced by a better understanding of plant transcriptional responses to arsenic. Many plant studies have demonstrated the direct involvement of thiol-containing molecules (glutathione, phytochelatins, etc.) in arsenic detoxification, however more robust approaches (i.e. microarrays) should help clarify how arsenic affects plant physiological processes on a global scale. The goals of our study were to test our hypothesis that many genes would be differentially expressed in response to arsenate stress and to identify genes as putative players in As (V) detoxification using Arabidopsis as a model. In this paper, we investigate the transcriptional responses to As (V) in <italic>Arabidopsis thaliana </italic>using oligonucleotide microarrays. Our results demonstrate that As (V) stress strongly induces Cu/Zn superoxide dismutase (SOD) activity, but represses the production of Fe SODs. Our microarray data also suggest the involvement of other antioxidant genes, various transcription factors, tonoplast proteins, and proteins associated with cell wall growth. Of particular interest, we report that As (V) stress represses numerous genes induced by Pi starvation. We discuss the physiological implications of these findings, and suggest new avenues for research of arsenic metabolism in plants.</p>
<title>Methods</title> <title>Plants and growth conditions</title> <p>Seeds of <italic>Arabidopsis thaliana </italic>ecotype Columbia plants were surface sterilized and plated on agar-solidified MS culture medium supplemented with B5 vitamins, 10% sucrose, 2% Gelrite<sup>®</sup>, pH 5.8. Phosphate is supplied as 1.25 mM KH<sub>2</sub>PO<sub>4 </sub>in the culture medium. Arsenic-treated plates were supplemented with 100 μM potassium arsenate (Sigma) according to a previously determined sub-lethal growth response curve. Plates were cold stratified at 4°C for 24 hrs and then placed in a growth chamber at 25°C under a 16 hr photoperiod. At each time point (3 d, 10 d), 2 g of whole plant material (shoots + roots) was harvested from each plate, frozen in liquid nitrogen, and subjected to RNA isolation using Trizol<sup>® </sup>reagent (Invitrogen, Carlsbad, CA) according to manufacturer's protocol. A total of three biological replicates were assayed (3 control, 3 treated) where each pooled 2 g sample represented a single biological replicate.</p> <title>Microarray experiments and aRNA labeling</title> <p>Total RNA from six biological replicates were purified using RNeasy MiniElute columns (Qiagen, Valencia, CA). A total of 1.25 μg of purified total RNA was subjected to Aminoallyl Message Amp II kit (Ambion, Austin, TX) first strand cDNA synthesis, second strand synthesis, and <italic>in vitro </italic>transcription for amplified RNA (aRNA) synthesis. aRNA was purified according to manufacturers protocol (Ambion, Austin, TX) and quantified using a Nanodrop spectrophotometer. Two 4 μg samples of aRNA were labeled with Cy3 and Cy5 monoreactive dyes (Amersham Pharmacia, Pittsburgh, PA) in order to conduct a dye swap technical replicate for each biological replicate. Each aRNA sample was brought to dryness in a Speedvac and dissolved in 5 μL of 0.2 M NaHCO<sub>3 </sub>buffer. Five microliters of Cy3 or Cy5 (in DMSO) was added to each sample and incubated for 2 hrs in the dark at RT. Labeled aRNA was purified according to kit instructions (Ambion, Austin, TX) and quantified using the Nanodrop spectrophotometer. One-hundred pmol Cy3- and Cy5-labeled aRNA targets were denatured by incubating at 65°C for 5 min and added to a hybridization mix containing 9 μl 20× SSC, 5.4 μl Liquid Block (Amersham Pharmacia, Pittsburgh, PA), and 3.6 μl 2% SDS for a 90 μl total volume.</p> <title>Hybridization and data analysis</title> <p>Microarrays comprised of 70-mer oligonucleotides obtained from the University of Arizona (see Availability and requirements section for URL) were immobilized by rehydrating the slide over a 50°C waterbath for 10 s and snap drying on a 65°C heating block for 5 s for a total of four times. Slides were UV-crosslinked at 180 mJ in a UV cross-linker (Stratagene, La Jolla, CA). The slides were then washed in 1% SDS, dipped in 100% EtOH five times followed by 3 min shaking. Slides were spun dry at 1000 rpm for 2 minutes and immediately placed in a light-proof box. The 90 μl hybridization mix was pipetted onto a microarray slide underneath a lifterslip (Lifterslip, Portsmouth, NH) and placed in a hybridization chamber (Corning, Corning, NY) overnight at 55°C. After hybridization, slides were washed in 2× SSC, 0.5% SDS for 5 minutes at 55°C, 0.5× SSC for 5 minutes at room temperature, and 0.05× SSC for 5 minutes at room temperature. Slides were then spun dry at 1000 rpm in a Sorvall centrifuge and scanned with a GenePix 4000B scanner (Axon Instruments, Inc., Union City, CA). The intensity variation was removed by fitting a loess regression using SAS 9.1 (SAS, Cary, NC). Data were log-2 transformed and statistically analyzed using rank product statistics as described by [##REF##15327980##43##] to identify differentially expressed genes. Bioconductor Rank Prod package was used to perform the rank product analysis [##REF##16982708##44##,##REF##15461798##45##]. Significantly different genes reported in this study exhibited <italic>P </italic>&lt; 0.001, as designated by the rank product analysis. The false discovery rate (FDR) [##UREF##6##46##] value obtained was based on 10,000 random permutations. Since 10,000 random permutations was very computer intensive, 1000 random permutations were performed 10 different times each time starting with a different random seed number and the average FDR value calculated was used for further analysis. The genes that had FDR values less than or equal to 0.01 were considered as differentially expressed. Data for all microarray experiments were submitted to the NCBI GEO microarray database and can be viewed under the accession GSE10425.</p> <title>Microarray Data Quality Control</title> <p>Global gene expression profiling comparing arsenate-treated <italic>Arabidopsis </italic>plants with control was carried out to better understand the mechanisms of plant response to arsenate stress and to identify genes involved in arsenic metabolism. For microarray data quality control, we examined both dye dependent effects and distribution of the ratio after normalization. We have included an additional file that illustrates the quality of microarray experiments, as well as the overall gene expression pattern [see Additional file ##SUPPL##1##2##]. Additional file ##SUPPL##1##2a## shows the normalized M vs. A plot, which was generated as a scatter plot of log intensity ratios <italic>M = log</italic><sub>2 </sub><italic>(R/G) </italic>versus average log intensities <italic>A = log</italic><sub>2</sub><italic>(R*G)/2</italic>, where R and G represent the fluorescence intensities in the Cy3 and Cy5 channels, respectively [##REF##12154381##47##]. As shown by the figure, Loess normalization effectively removed dye dependent effects in the microarray and rendered evenly distributed ratios across all signal intensities. The histogram suggests a normal distribution of the logarithm 2-based transformed ratio [see Additional file ##SUPPL##1##2b##]. Overall, the microarray experiments generated high quality data without significant dye-dependent effects and skewness of ratio distribution.</p> <title>Gene ontology analysis</title> <p>Gene ontology annotations were translated from microarray data using the GO annotations bioinformatics tool available at The Arabidopsis Information Resource Web site <ext-link ext-link-type="uri" xlink:href="http://www.arabidopsis.org/tools/"/> where results were based on molecular function.</p> <title>RT-PCR amplification</title> <p>Total RNA was extracted from <italic>Arabidopsis thaliana </italic>ecotype Columbia grown for ten days as described for the microarray experiment. Five micrograms of total RNA was reverse-transcribed with oligo(dT)<sub>20 </sub>primers using the Superscript III first-strand cDNA synthesis kit (Invitrogen, Carlsbad, CA). RT PCR was performed using the ABI 7000 Sequence Detection System (Applied Biosystems, Foster City, CA). PCR was performed in a 15 μl reaction volume containing Power Sybr<sup>® </sup>PCR mix (Applied Biosystems, Foster City, CA) and gene-specific primers were designed with PrimerExpress software. Actin was used as the reference gene, and the primer sequences for <italic>Arabidopsis </italic>actin gene were AGTGGTCGTACAACCGGTATTGT (F) and GAGGAAGAGCATTCCCCTCGTA (R). After the RT PCR experiment, Ct number was extracted for both reference gene and target gene with auto baseline and manual threshold.</p> <title>Cluster Analysis</title> <p>The cluster analysis was conducted with MultiExperiment viewer Version 4.0 (TIGR, Rockville, MD) with logarithm 2 transformed ratio of treated vs. control samples from real-time PCR. The complete linkage hierarchical cluster was used to cluster the genes only. The color scheme is as shown in the figure, with repressed genes shown as green and red color indicating induced genes.</p> <title>SOD activity assay</title> <p>Total soluble protein was extracted from whole <italic>Arabidopsis </italic>plants (root + shoot) grown on plates as described above that were harvested at each respective time point. Total soluble protein was quantified by the method of Bradford [##REF##942051##48##] using BSA as a standard and 50 μg samples were loaded. Bovine SOD (Sigma) was used in each gel to serve as a positive control for SOD activity. Following electrophoretic separation on a 10% non-denaturing polyacrylamide gel, SOD activity was determined as described by Beauchamp and Fridovich (1971) and modified by Azevedo <italic>et al</italic>. [##UREF##7##49##]. The gels were rinsed with DDI water and incubated in the dark for 30 min at room temperature in a reaction mixture containing 50 mM potassium phosphate buffer (pH 7.8), 1 mM EDTA, 0.05 mM riboflavin, 0.1 mM nitroblue tetrazolium and 0.3% (v/v) TEMED. Following incubation, gels were rinsed with DDI water and illuminated in water until SOD bands were visible. The gels were then immersed in a 6% (v/v) acetic acid solution to stop the reaction. To confirm specificity of Cu/Zn-SOD activity, H<sub>2</sub>0<sub>2 </sub>and KCN were used as inhibitors as described by Azevedo <italic>et al</italic>. [##UREF##7##49##] and modified by Vitoria <italic>et al</italic>. [##REF##11397437##50##]. Mn-SOD is resistant to both inhibitors, Fe-SOD is resistant to KCN and inhibited by H<sub>2</sub>0<sub>2</sub>, and Cu/Zn-SOD is inhibited by both inhibitors, thus allowing classification of SOD activity. Prior to SOD staining, gels containing lanes in triplicate were cut into three parts; one gel was treated as described above, the second and third parts were incubated for 20 min in 100 mM potassium phosphate buffer (pH 7.8) containing either 2 mM KCN or 5 mM H<sub>2</sub>O<sub>2</sub>, respectively. Following incubation, gels were rinsed with DDI water and then stained for SOD activity.</p>
<title>Results</title> <title>Root growth under As (V) stress</title> <p>Arsenate exposure resulted in reduced Arabidopsis root growth and branching (Figure ##FIG##0##1##). Exposure of 50 μM As (V) resulted in significantly reduced Arabidopsis root growth (Figure ##FIG##0##1C##). In addition to known relevant physiological data, the exposure study was used to determine suitable As (V) exposure for the microarray study. We noted no seed germination effects with regard to arsenate treatments, notably at 100 μM, the arsenate concentration used for transcriptomics experiments.</p> <title>Gene ontology for genes affected by As (V)</title> <p>Forty-six genes were induced by As (V) treatment (i.e., exceeded a fold-change threshold of 1.5 and met a significance criteria of <italic>P </italic>&lt; 0.001; FDR &lt; 1%) as indicated by microarray analysis. The largest functional categories affected included unknown function, hydrolase, and antioxidant activity. Other functional categories affected by As (V) included genes with transferase, kinase, lyase, transporter, and binding activity [see Additional file ##SUPPL##0##1##; Table ##TAB##0##1##]. Alternatively, 113 genes were repressed by As (V) (i.e. exceeded a fold-change threshold of -1.5 and met a significance criteria of <italic>P </italic>&lt; 0.001; FDR &lt; %1), with unknown function, hydrolase, and binding activity representing the largest categories. Genes with transporter, kinase, transferase, and transcriptional regulator activity were also repressed by As (V) [see Additional file ##SUPPL##0##1##; Table ##TAB##1##2##]. Differentially expressed As (V)-induced and -repressed genes are listed below (Table ##TAB##0##1## and Table ##TAB##1##2##, respectively) and complete lists of all genes affected by As (V) stress are also included as additional files [see Additional files ##SUPPL##2##3## and ##SUPPL##3##4##]. Most interestingly, it was discovered that As (V) stress repressed transcription of many genes involved in the phosphate starvation response, and also repressed several transcriptional factors. Several genes involved in oxidative stress were also highly modulated in response to As.</p> <title>Superoxide dismutases</title> <p>SODs represented the highest ranked of both significantly induced as well as repressed genes in response to As (V) stress (Tables ##TAB##0##1## and ##TAB##1##2##), therefore these genes presented logical primary targets for the validation of our microarray data. Results demonstrated 4.57-fold induction of a chloroplast Cu/Zn SOD (at2g28190), 2.41-fold induction of a Cu/Zn SOD (at1g08830), as well as a 3.16-fold induction of an SOD copper chaperone (at1g12520). Alternatively, Fe SOD (at4g25100) transcripts were downregulated in response to arsenic stress (-5.17-fold change). These findings were confirmed with quantitative RT-PCR (qRT-PCR) (Table ##TAB##2##3##).</p> <p>Based upon our observations of transcript-level changes in SOD gene expression, we performed a nondenaturing PAGE enzyme assay of superoxide dismutase activity [##REF##4943714##16##] to assess whether the observed specific changes in SOD transcript correlated with enzyme activity. This method enables the distinction between the three SOD isoenzymes found in Arabidopsis (CuZnSOD, FeSOD, and MnSOD) by using inhibitors of specific SODs. Gels were preincubated with KCN, which inhibits CuZn SOD, as well as H<sub>2</sub>O<sub>2</sub>, which inhibits both CuZn SOD and Fe SOD. MnSOD is resistant to both inhibitors (Figure ##FIG##1##2##). Plants were harvested from control plates containing no arsenate and treated plates containing 100 μM arsenate at seven-, ten-, and thirteen days post-germination. Irrespective of harvest date, CuZnSOD activity was strongly induced by arsenate treatment, whereas FeSOD activity was repressed, and MnSOD showed no change in activity, therefore providing sufficient evidence to confirm our microarray results. Transcription of other antioxidant genes (i.e., peroxidases, glutathione-S transferases, catalase) were indicated by our microarray experiment as affected by arsenate stress, however these genes were not included in our qRT-PCR validation.</p> <title>Transcription factors</title> <p>Our microarray experiment indicated that eight different genes encoding proteins with known transcription factor activity all displayed lower expression levels in As (V)-stressed plants (Table ##TAB##1##2##). One of these transcription factors (at1g12610) encodes a member of the DREB subfamily A-1 of the <italic>ERF/AP2 </italic>transcription factor family (<italic>DDF1</italic>). One other AP2-domain-containing transcription factor (at4g34410) that encodes a member of the ERF (ethylene response factor) subfamily B-3 of the <italic>ERF/AP2 </italic>transcription factor family was also repressed in response to As (V). Two zinc finger (C2H2 type) genes (at3g46090, at3g46080) encoded a ZAT7 and a protein similar to ZAT7, respectively. Also exhibiting lower expression in As (V)-treated plants were three members of the WRKY family of transcription factors (at2g38470, at4g23810, at1g80840), <italic>WRKY33</italic>, <italic>WRKY53</italic>, and <italic>WRKY40</italic>, respectively as well as one gene encoding NAC domain containing protein 81.</p> <title>As (V) represses genes involved in phosphate starvation response</title> <p>A notable transcriptional trend is that As (V) stress results in repression of some genes involved in the Pi starvation response. Of particular interest, a P-type cyclin (at5g61650) that was affected by As (V)-stress shares significant homology to the <italic>PHO80 </italic>gene from yeast. We performed qRT-PCR for this gene and found that its expression was actually strongly repressed at both day 3 and day 10 (Table ##TAB##2##3##; Figure ##FIG##2##3##). Three genes that were repressed by As (V) in this study have also been reported to be repressed by Pi starvation [##REF##16085708##17##]. Interestingly, the three highest ranking differentially expressed genes found to be strongly induced by Pi starvation (at1g73010 &gt; at5g20790 &gt; at1g17710, respectively) [##REF##17177879##18##], were also repressed by As (V) in our study. These genes are of particular interest on account of their unknown function. Quantitative RT-PCR confirmed that transcription of these genes was strongly repressed at both 3 day and 10 day time points. Additionally, the qRT-PCR data indicate that these genes were more repressed at day 3 than at day 10 (Table ##TAB##2##3##; Figure ##FIG##2##3##). These data also corroborated our microarray experiments, which reflect global expression ratios at 10 days post germination (Table ##TAB##2##3##).</p> <title>Sulfate assimilation</title> <p>The role that thiol groups play in arsenic detoxification has been well characterized, therefore we expected to see induction of genes involved in sulfate assimilation and metabolism in response to arsenic stress. Ferredoxin (at1g10960), a key redox protein found in the chloroplast was As (V)-induced. Expression levels for another gene involved in the sulfate reduction pathway, 5'-adenylylsulfate reductase (<italic>APR3</italic>) (at4g21190) were also elevated in response to As (V) stress. This enzyme catalyzes the reduction of APS to sulfite using glutathione as an electron donor. Although not involved in sulfate assimilation, the cysteine-rich metal-binding protein, metallothionein (MT) 1A (at1g07600) was also induced. <italic>Arabidopsis </italic>knockout mutants that were generated for class 1 MTs accumulated significantly less aboveground As, Cd, and Zn, suggesting that class 1 MTs may play a role in metal and metalloid ion translocation [##REF##16240177##19##].</p> <title>Genes involved in cell wall assembly, architecture, and growth</title> <p>A wide range of genes encoding proteins involved in cell wall activities exhibit altered expression levels in response to As (V) (Table ##TAB##0##1##; Table ##TAB##1##2##). Peroxidases, which were indicated by microarray as affected by As (V) stress, are known to strengthen the cell wall in response to biotic stress via formation of lignin, extension cross-links, and dityrosine bonds [##REF##15856234##20##]. Additionally, As (V) affected transcription of numerous xyloglucan endotransglucosylase/hydrolases (XTHs) and glycosyl hydrolase genes (Table ##TAB##0##1##; Table ##TAB##1##2##), with the majority of these exhibiting lower expression in the presence of As (V).</p>
<title>Discussion</title> <title>Arsenic and oxidative stress</title> <title>Superoxide dismutases</title> <p>Increasing evidence from mammalian studies demonstrates that ROS are generated in response to exposure to inorganic forms of arsenic [##REF##9653147##21##, ####REF##11172004##22##, ##REF##12888263##23####12888263##23##]. The reduction of arsenic is linked with <italic>in vivo </italic>and <italic>in vitro </italic>ROS production in mammalian cells [##REF##9653147##21##], but little is known about the mechanisms by which arsenic-induced ROS generation occurs in plants. It is believed that the reduction of As (V) to As (III), which is well documented in plants, results in the production of ROS [##UREF##4##8##,##REF##9741595##24##]. However, this increase in ROS may also be the result of either depletion of glutathione or inhibition of antioxidant enzymes. Plants have evolved both nonenzymatic antioxidants (i.e., glutathione, ascorbate, and carotenoids), as well as antioxidant enzymes (i.e., superoxide dismutases, catalases, and peroxidases) to manage the balance of ROS in the cell.</p> <p>SODs represent a first line of defense by converting superoxide radicals to H<sub>2</sub>0<sub>2</sub>, whereas catalases and peroxidases remove H<sub>2</sub>O<sub>2</sub>. Three classes of SODs have been identified according to the active site metal cofactor: FeSOD, MnSOD, and Cu/ZnSOD. As (V) and As (III) were both shown to induce expression of glutathione S-transferases (GSTs), catalases, and SODs in <italic>Zea mays </italic>[##REF##9741595##24##]. An increase in SOD activity was correlated with an increase in As (V) treatment in <italic>Holcus lanatus </italic>[##REF##11351093##12##]. Higher levels of SOD, catalase, and ascorbate peroxidase were observed in <italic>Pteris vittata</italic>, an arsenic hyperaccumulator, than in arsenic-sensitive fern species <italic>Pteris ensiformis </italic>and <italic>Nephrolepsis exaltata </italic>[##REF##15781440##25##]. These researchers concluded that arsenic-induced increases in antioxidant enzymes levels may represent a secondary defensive mechanism against oxidative stress in <italic>Pteris vitatta </italic>and correspond with its arsenic accumulation and lack of toxicity symptoms. It was shown that <italic>Pteris vittata </italic>SOD, catalase, and peroxidase levels rose sharply in response to low levels of As (V), but leveled off at As (V) levels &gt; 20 mg kg<sup>-1</sup>, which was consistent with changes in biomass in the arsenic hyperaccumulator [##REF##14720474##26##].</p> <p>Although the strong induction of SODs in response to As (V) stress was not surprising, the dramatically lower levels of FeSODs were unexpected. We suggest the involvement of an NAC domain-containing transcription factor to explain the observed decrease in FeSOD transcription based on our microarray results (Table ##TAB##1##2##). One group recently generated transgenic plants to overexpress three different Arabidopsis NAC transcription factors and identified NAC-dependent genes using microarrays [##REF##15319476##27##]. Not only was at4g25100 (FeSOD) expression found to be NAC-dependent, but transcription of other genes we have observed to be repressed by As (V) stress also appear to be dependent on NAC-domain containing transcription factors. We continue this discussion more thoroughly in the following section on transcription factors.</p> <title>Peroxidases</title> <p>Peroxidases are functionally diverse and participate in two major cycles: the hydroxylic cycle where peroxidases regulate H<sub>2</sub>O<sub>2 </sub>levels and release ROS (·OH, HOO·) and the peroxidative cycle where various substrates (e.g. phenolic compounds) are oxidized or polymerized. Their involvement in a broad range of physiological processes allows peroxidase expression in all plant organs from germination to early senescence, however they are predominantly expressed in the roots [##REF##15856234##20##]. It is not surprising that peroxidases seem to be affected by arsenate stress (Table ##TAB##0##1##; Table ##TAB##1##2##), especially in consideration of the elevated SOD activity, which produces H<sub>2</sub>O<sub>2 </sub>as a product of superoxide radical dismutation.</p> <title>Transcription factors</title> <p>Our microarray data corroborate those of Tran <italic>et al</italic>. [##REF##15319476##27##], suggesting the involvement of a different NAC domain-containing transcription factor (at5g08790) in expression of FeSOD, as well as several other genes known to exhibit NAC-dependent expression. NAC proteins comprises a large gene family (&gt; 100 members in Arabidopsis) of plant-specific transcription factors that have roles in wide-ranging processes such as development, defense, and abiotic stress response [##REF##15708345##28##]. Microarray experiments were carried out on NAC-overexpression Arabidopsis mutants to discover genes exhibiting dependence on NAC transcription factors for transcription [##REF##15319476##27##]. We speculate that repression of <italic>NAC81 </italic>(at5g08790) in As (V)-stressed Arabidopsis may be responsible for the observed repression of FeSOD (at4g25100), ferritin 1 (FER 1) (at5g01600), <italic>XTH15 </italic>(at4g14130), <italic>XTH24 </italic>(at4g30270), <italic>erd1 </italic>ATP-dependent Clp protease ATP-binding subunit (at5g51070), and a branched-chain amino acid amino transferase 2 (at1g10070), as these genes were reported as exhibiting NAC-dependent expression [##REF##15319476##27##].</p> <title>As (V) stress represses genes induced by Pi deprivation</title> <p>Although phosphate is undoubtedly one of the most biologically important nutrients, its availability in soils is quite low. Therefore, plants have evolved mechanisms to maximize Pi accessibility/availability, such as increased root hair growth, lateral root branching, and induction of phosphate transporters and phosphatases [##REF##11493657##29##]. Certain phosphate starvation-induced genes have evolved to release phosphate from plasma membranes by hydrolyzing phospholipids under conditions of low Pi availability, as phospholipids comprise a major Pi pool <italic>in planta </italic>[##REF##11906834##30##]. Conversion from phospholipids to galactolipids is one such strategy and can result from the activity of monogalactosyldiacylglycerol (MGDG) synthase or digalactosyl diacylglycerol (DGDG) synthase [##REF##16762032##31##]. Arabidopsis plants expressing <italic>MGD</italic>2 and <italic>MGD</italic>3 promoter-GUS fusion constructs showed that under Pi starvation, <italic>MGD</italic>3::GUS was expressed in apices of serrated edges (hydathodes) and in the lateral root branch [##REF##16762032##31##]. Through investigation of Arabidopsis MGDG synthase gene expression under Pi starvation, these authors showed that global changes in plant membranes under Pi deprivation are tightly regulated by Pi signaling and that signal transduction through a Pi-sensing mechanism is responsible for regulating MGDG synthase gene expression [##REF##16762032##31##]. We report here that the expression of <italic>MGD3 </italic>(at2g11810) is lower in As (V)-treated Arabidopsis at 3 days and 10 days (Table ##TAB##2##3##; Figure ##FIG##2##3##). Therefore, it is conceivable that our observations may either reflect a Pi/As (V) sensing mechanism or simply the lower number of lateral roots in As (V)-stressed plants (Figure ##FIG##0##1##). <italic>SENESCENCE RELATED GENE </italic>3 (SRG3; at3g02040), a glycerophosphoryl diester phosphor-diesterase, is believed to participate in processes similar to those of the MGDG synthase genes [##REF##17921353##32##]<italic>SRG</italic>3 had lower transcript abundance in As (V)-treated plants in our microarray study (along with other senescence-associated proteins) (Table ##TAB##1##2##), as well as, in 3 day and 10 day As (V)-treated plants (Table ##TAB##2##3##; Figure ##FIG##2##3##). A type 5 acid phosphatase (<italic>ACP</italic>5; at3g17790) was also repressed in our As (V)-treated plants as indicated by microarray (Table ##TAB##1##2##) and was strongly repressed in our qRT-PCR validation experiments at both 3 day and 10 day time points (Table ##TAB##2##3##; Figure ##FIG##2##3##). In Arabidopsis, <italic>ACP</italic>5 has been shown to be induced by H<sub>2</sub>O<sub>2</sub>, but not by paraquat or salicylic acid and is thought to be involved in both phosphate mobilization and in the metabolism of reactive oxygen species [##REF##10504579##33##]. In contrast, <italic>ACP</italic>5 was strongly repressed by As (V) despite elevated SOD levels, which generate H<sub>2</sub>O<sub>2</sub>. Therefore, further study is required to determine the specific cause of As (V)-mediated <italic>ACP</italic>5 repression.</p> <p>Recent investigations into the genome-scale transcriptional changes to phosphate deprivation in <italic>Arabidopsis </italic>have elucidated a broad range of genes involved in phosphate metabolism [##REF##16085708##17##,##REF##17177879##18##]. Our microarray data suggested that many genes repressed by As (V) stress have been reported by others [##REF##16085708##17##,##REF##17177879##18##] to be induced in response to Pi deprivation in <italic>Arabidopsis thaliana</italic>. Because As (V) behaves as a phosphate analog, it is likely that this observation can be explained by a saturation effect of the phosphate analog, As (V), thereby misleading metabolic and regulatory perception of the toxic metalloid as an abundant supply of Pi. However, arsenate likely disrupts critical biological processes that involve reversible phosphorylation, as well as pathways for phosphate signaling, but even under arsenate stress, Arabidopsis accumulates much higher concentrations of As in the root than is translocated to the shoot. In another study, when wild-type (Columbia ecotype) Arabidopsis plants were grown on 100 μM sodium arsenate for 3 weeks, low concentrations of arsenic were accumulated in the shoot, whereas high concentrations of arsenic were observed in roots [##REF##16567632##34##]. However, when the arsenate reductase homolog (ACR2) was silenced, arsenate was translocated to the shoot at concentrations that classified as hyperaccumulation [##REF##16567632##34##]. Nevertheless, the signaling mechanisms by which plants distinguish between As (V) and phosphate are unknown and other mechanisms of As detoxification and storage besides the well documented phytochelatin response [##REF##10759512##9##, ####REF##16594069##10##, ##UREF##5##11##, ##REF##11351093##12####11351093##12##] may exist.</p> <p>In order to confirm the observation that As (V) stress represses genes involved in phosphate starvation/acquisition, we performed qRT-PCR on some of the more interesting candidates (Table ##TAB##2##3##; Figure ##FIG##2##3##). We are particularly interested in elucidating pathways involved in As (V) signaling in plants. The P-type cyclin (at5g61650) that was affected by As (V) (Table ##TAB##2##3##; Figure ##FIG##2##3##) shares significant homology to the <italic>PHO80 </italic>gene from yeast. Cyclins bind and activate cyclin-dependent kinases, which play key roles in cell division via phosphorylation of critical substrates, such as the retinoblastoma protein, transcription factors, nuclear laminar proteins, and histones [##REF##9442875##35##]. Interestingly, it was demonstrated that expression of this cyclin from Arabidopsis restored the phosphate signaling pathway in a <italic>PHO80</italic>-deficient yeast mutant, suggesting a putative key Pi signaling role [##REF##15197472##36##].</p> <p>Protein kinases play crucial roles in signal transduction pathways in all eukaryotes [##REF##7834742##37##]. At3g08720 (<italic>ATPK</italic>19) is one of two nearly identical kinase genes in Arabidopsis that encode for proteins that share high sequence homology with the mammalian 40S ribosomal protein kinases S6K1 and S6K2 [##REF##11008486##38##]. <italic>ATPK</italic>19 was demonstrated to be the functional plant homolog of mammalian p70s6k when ectopic expression of this gene specifically phosphorylated ribosomal protein S6 derived from either plant or animal [##REF##9528776##39##]. <italic>ATPK</italic>19 has recently been implicated as a crucial nodal point in a network evolved for integrating stress signals with plant growth regulation [##REF##16377759##40##]. Lower expression levels observed for <italic>ATPK</italic>19 in As (V)-treated plants, which was most severe at day 3 (Table ##TAB##2##3##; Figure ##FIG##2##3##), lends us to conclude that As (V) stress may suppress plant growth through the downregulation of this growth-regulating kinase, possibly as a result of the chemical similarity between As (V) and phosphate. Alternatively, the downregulation of <italic>ATPK</italic>19 may result from the more general stress responses imposed by the toxic metalloid (e.g. oxidative stress, sulfhydryl group binding, etc.).</p> <p>Our results are in agreement with the recently proposed ideas of Catarecha <italic>et al</italic>. [##REF##17400898##41##] who studied an Arabidopsis mutant that displayed enhanced arsenic accumulation. These authors identified a Pi transporter (PHT1;1) mutant with a decreased rate of As (V) uptake and increased As (V) accumulation. By comparing gene expression of the mutant with wild-type plants, it was shown that in Arabidopsis, As (V) rapidly repressed genes involved in the Pi starvation response and induced the expression of other As (V)-responsive genes [##REF##17400898##41##]. Interestingly, the repression of Pi starvation genes was shown to be specific for As (V), whereas the As (V)-induced genes were also induced by As (III). A model resulted that suggests arsenic acts via two separate signaling pathways [##REF##17400898##41##]. Because of the chemical similarity of As (V) and Pi, As (V) fools the Pi sensor, thus initiating the repression of the Pi starvation response. Although our microarray experiments did not detect differential expression of any high-affinity Pi transporter, which may be due to differences in experimental approach, Catarecha <italic>et al</italic>. [##REF##17400898##41##] illustrated the high sensitivity of the Pi transporter, PHT1;1, to As (V) and suggested that plants have evolved an As (V) sensing system whereby As (V) and Pi signaling pathways oppose each other to protect the plant from arsenic toxicity. Based on our results, it is conceivable that the P-type cyclin (at5g61650) and <italic>ATPK</italic>19 (at3g08720) may be involved in As (V) sensing, but further study is required to confirm this finding.</p> <p>Our comparison of As (V)-repressed genes that have also been shown to be induced by Pi deprivation elucidate some promising candidates for future studies. For example, we are particularly interested in genes with unknown function that are strongly induced in both roots and leaves by Pi starvation (i.e. at1g73010; at1g17710; at2g04460; at5g20790; [##REF##16085708##17##], supplemental data; [##REF##17177879##18##]). Both at1g73010 and at1g17710 are described as phosphoric monoester hydrolases (see Availability and requirements section for URL), but to our knowledge, these have not been studied in this regard. Most recently, a study described gene networks for the Arabidopsis transcriptome based on the graphical Gaussian model of global-scale transcriptional studies [##REF##17921353##32##]. In a constructed subnetwork of genes involved in phosphate starvation, SRG3, at1g73010, at3g17790 (ACP5), at2g11810, and at5g20790 were all closely linked, suggesting their critical roles in phosphate metabolism in Arabidopsis [##REF##17921353##32##]. Based on their strong induction in response to Pi starvation [##REF##16085708##17##,##REF##17177879##18##], it is reasonable to conceive that at1g73010 and at1g17710 have evolved as Pi scavengers for increasing Pi availability. We have confirmed the downregulation of these genes in response to As (V) at both 3 and 10 day time points (Table ##TAB##2##3##; Figure ##FIG##2##3##). In this study, at2g04460 transcript levels were strongly repressed at 3 and 10 day time points, whereas at5g20790 was repressed at day 3 and day 10 (Table ##TAB##2##3##; Figure ##FIG##2##3##). Interestingly, at2g04460 encodes for a putative retroelement <italic>pol </italic>polyprotein that has been reported as highly expressed in salt overly sensitive (<italic>sos</italic>) Arabidopsis mutants [##REF##11351099##42##]. Because the function of these two Pi starvation-induced genes is unknown, these putative gene candidates may provide opportunities for gaining insight into As (V)/Pi dynamics in <italic>Arabidopsis thaliana</italic>.</p> <p>The data presented here have led to the development of new hypotheses for future research. The potential antagonistic effects of various arsenate and Pi concentrations on the expression of the aforementioned genes in Arabidopsis are poorly understood. Additionally, the efficiency of arsenate reduction and subsequent detoxification via phytochelatins or glutathione is poorly understood. Under conditions of arsenate stress (i.e., 100 μM), perhaps the cellular concentrations of arsenate surpass those that may be efficiently reduced by glutathione or arsenate reductase, thus allowing free arsenate to interfere with biological reactions that involve phosphate. Additionally, the selected arsenate concentration to employ in this study could have resulted in free arsenite after reduction <italic>in vivo </italic>that would likely have deleterious consequences. Therefore, it is reasonable to conceive that the observed transcriptional responses, as well as the impaired phenotype seen in this study, may be reflective of either arsenate, arsenite, or both.</p>
<title>Conclusion</title> <p>Our data show that in Arabidopsis, Cu/Zn SODs are strongly induced in response to As (V) stress, while Fe SOD expression is repressed. We also demonstrate that As (V) stress results in the repression of genes involved in phosphate acquisition, redistribution, and phosphorylation, which supports a recent study [##REF##17400898##41##] that suggests As (V) and Pi signaling pathways act in opposition to protect plant health. Although this study identifies some interesting targets for exploring As (V) metabolism, further studies using Arabidopsis mutants with altered expression of these genes are necessary to elucidate their biological significance, as well as to clarify new pathways involved in arsenic signaling in plants.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Arsenic is toxic to plants and a common environmental pollutant. There is a strong chemical similarity between arsenate [As (V)] and phosphate (Pi). Whole genome oligonucleotide microarrays were employed to investigate the transcriptional responses of <italic>Arabidopsis thaliana </italic>plants to As (V) stress.</p> <title>Results</title> <p>Antioxidant-related genes (i.e. coding for superoxide dismutases and peroxidases) play prominent roles in response to arsenate. The microarray experiment revealed induction of chloroplast Cu/Zn superoxide dismutase (SOD) (at2g28190), Cu/Zn SOD (at1g08830), as well as an SOD copper chaperone (at1g12520). On the other hand, Fe SODs were strongly repressed in response to As (V) stress. Non-parametric rank product statistics were used to detect differentially expressed genes. Arsenate stress resulted in the repression of numerous genes known to be induced by phosphate starvation. These observations were confirmed with qRT-PCR and SOD activity assays.</p> <title>Conclusion</title> <p>Microarray data suggest that As (V) induces genes involved in response to oxidative stress and represses transcription of genes induced by phosphate starvation. This study implicates As (V) as a phosphate mimic in the cell by repressing genes normally induced when available phosphate is scarce. Most importantly, these data reveal that arsenate stress affects the expression of several genes with little or unknown biological functions, thereby providing new putative gene targets for future research.</p>
<title>Availability and requirements</title> <p>US Environmental Protection Agency: <ext-link ext-link-type="uri" xlink:href="http://www.atsdr.cdc.gov/cercla/05list.html"/></p> <p>University of Arizona: <ext-link ext-link-type="uri" xlink:href="http://ag.arizona.edu/microarray/"/></p> <p>The Arabidopsis Information Resource:  <ext-link ext-link-type="uri" xlink:href="http://www.arabidopsis.org/tools/"/></p> <title>Authors' contributions</title> <p>Jason Abercrombie conceived of the study, its experimental design, and overall data interpretation and drafting of the manuscript. Matt Halfhill contributed to the experimental design, data interpretation, and technical assistance. Priya Ranjan, Murali Rao, and Arnold Saxton all conducted the non-parametric statistical analyses of the microarray data and reprogrammed the RankProd package to carry out 10,000 random permutations of the data and provided intellectual input for the statistical interpretation of the data. C. Neal Stewart, Jr. supervised organization of the manuscript, provided critical analyses of the data, and gave final approval of its readiness for submission. All authors read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We appreciate our collaborators at Edenspace Systems Corporation for their valuable cooperation, especially, Mark Elless, David Lee, and Bruce Ferguson. We thank Laura Abercrombie and Reggie Millwood for technical assistance. Funding was provided by NIH and NSF grants as well as the Tennessee Agriculture Experiment Station and the Ivan Racheff Endowment. We would also like to thank the BMC Plant Biology editors and 4 anonymous reviewers for their thorough reviews and constructive criticism of this manuscript.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Phenotype of arsenate stress in Arabidopsis.</bold> Representative <italic>Arabidopsis thaliana </italic>"Columbia" plants grown for 10 days on MS medium containing either (Panel A) 0 μM As (V) or (Panel B) 100 μM arsenate. C. Histogram of root length vs. As (V) concentration.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>SOD activity in <italic>Arabidopsis thaliana </italic>'Col' grown on medium containing 100 μM potassium arsenate.</bold> A, superoxide dismutase activity without inhibitors, B, Gels were preincubated with KCN (which inhibits CuZn SOD), C, H<sub>2</sub>O<sub>2 </sub>added as an inhibitor (which inhibits both CuZn SOD and Fe SOD). Lane 1, purified bovine SOD positive control. Lane 2, control plants harvested at 7 days, Lane 3, arsenate-treated plants harvested at 7 days. Lane 4, control plants harvested at 10 days, Lane 5, arsenate-treated plants harvested at 10 days. Lane 6, control plants harvested at 13 days, Lane 7, arsenate-treated plants harvested at 13 days.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Cluster diagram of ΔΔCt values in response to 100 μM arsenate at 3 and 10 day time points as determined by qRT-PCR.</bold> Bright green rectangles represent strongly repressed gene expression (-4-fold and lower), whereas genes strongly induced (4-fold and higher) are shown in bright red. Arabidopsis actin was used as a reference gene to calculate ΔΔCt.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Gene ontology based on molecular function for induced genes of arsenic-treated <italic>Arabidopsis thaliana </italic>Columbia plants.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Molecular function</td><td align="left">Gene ID</td><td align="left">Locus</td><td align="left">FC</td><td align="center">RP</td></tr></thead><tbody><tr><td align="left">Antioxidant activity</td><td align="left">Peroxiredoxin Q</td><td align="left">at3g26060</td><td align="left">1.53</td><td align="right">151.4</td></tr><tr><td/><td align="left">peroxidase</td><td align="left">at5g64100</td><td align="left">2.50</td><td align="right">7.3</td></tr><tr><td/><td align="left">peroxidase</td><td align="left">at1g05250</td><td align="left">1.90</td><td align="right">49.0</td></tr><tr><td/><td align="left">peroxidase 57 (PER57) (P57) (PRXR10)</td><td align="left">at5g17820</td><td align="left">1.68</td><td align="right">56.6</td></tr><tr><td/><td align="left">peroxidase</td><td align="left">at1g05240</td><td align="left">2.05</td><td align="right">16.5</td></tr><tr><td/><td align="left">superoxide dismutase [Cu-Zn], chloroplast</td><td align="left">at2g28190</td><td align="left">4.57</td><td align="right">2.3</td></tr><tr><td/><td align="left">superoxide dismutase [Cu-Zn], (SODCC) (CSD1)</td><td align="left">at1g08830</td><td align="left">2.41</td><td align="right">8.6</td></tr><tr><td/><td align="left">superoxide dismutase copper chaperone</td><td align="left">at1g12520</td><td align="left">3.16</td><td align="right">5.2</td></tr><tr><td align="left">Metal ion binding</td><td align="left">metallothionein-like protein 1A, (MT-1A)</td><td align="left">at1g07600</td><td align="left">1.67</td><td align="right">41.9</td></tr><tr><td/><td align="left">ferredoxin, chloroplast</td><td align="left">at1g10960</td><td align="left">1.53</td><td align="right">95.6</td></tr><tr><td align="left">Kinase activity</td><td align="left">leucine-rich repeat transmembrane protein kinase</td><td align="left">at3g24240</td><td align="left">1.59</td><td align="right">63.2</td></tr><tr><td/><td align="left">Cyclin-dependent protein kinase</td><td align="left">at5g61650</td><td align="left">1.64</td><td align="right">66.9</td></tr><tr><td align="left">Oxygen binding</td><td align="left">non-symbiotic hemoglobin 1 (HB1) (GLB1)</td><td align="left">at2g16060</td><td align="left">1.59</td><td align="right">94.2</td></tr><tr><td align="left">Hydrolase activity</td><td align="left">ATPase, BadF/BadG/BcrA/BcrD-type family</td><td align="left">at1g30540</td><td align="left">1.62</td><td align="right">75.5</td></tr><tr><td/><td align="left">myrosinase-associated protein</td><td align="left">at1g54010</td><td align="left">1.54</td><td align="right">90.5</td></tr><tr><td/><td align="left">myrosinase-associated protein</td><td align="left">at1g54000</td><td align="left">1.68</td><td align="right">47.6</td></tr><tr><td/><td align="left">xyloglucan:xyloglucosyl transferase</td><td align="left">at4g37800</td><td align="left">1.61</td><td align="right">72.0</td></tr><tr><td/><td align="left">glycosyl hydrolase family 1 protein</td><td align="left">at3g09260</td><td align="left">1.67</td><td align="right">50.8</td></tr><tr><td align="left">Isomerase activity</td><td align="left">peptidyl prolyl cis-trans isomerase</td><td align="left">at3g62030</td><td align="left">1.64</td><td align="right">53.1</td></tr><tr><td align="left">Lyase activity</td><td align="left">ribulose bisphosphate carboxylase small chain 2B</td><td/><td/><td/></tr><tr><td/><td align="left">ribulose bisphosphate carboxylase small chain 3B</td><td/><td/><td/></tr><tr><td align="left">Alcohol dehydrogenase activity</td><td align="left">alcohol dehydrogenase (ADH)</td><td align="left">at1g77120</td><td align="left">1.74</td><td align="right">46.5</td></tr><tr><td align="left">Nitrate reductase activity</td><td align="left">nitrate reductase 1 (NR1)</td><td align="left">at1g77760</td><td align="left">1.77</td><td align="right">41.8</td></tr><tr><td align="left">Sulfate reduction</td><td align="left">5'-adenylylsulfate reductase (APR3)</td><td align="left">at4g21990</td><td align="left">1.53</td><td align="right">112.1</td></tr><tr><td align="left">Molecular function unknown</td><td align="left">Photoassimilate-responsive protein</td><td align="left">at3g54040</td><td align="left">1.59</td><td align="right">112.1</td></tr><tr><td/><td align="left">Expressed protein</td><td align="left">at1g09310</td><td align="left">1.67</td><td align="right">63.1</td></tr><tr><td/><td align="left">Replication protein</td><td align="left">At5g35260</td><td align="left">1.62</td><td align="right">975.3</td></tr><tr><td/><td align="left">Hypothetical protein related to GB:AAD15331</td><td align="left">at2g06480</td><td align="left">1.58</td><td align="right">81.6</td></tr><tr><td/><td align="left">DREPP plasma membrane polypeptide-related</td><td align="left">at5g44610</td><td align="left">1.54</td><td align="right">106.5</td></tr><tr><td/><td align="left">Pentatricopeptide repeat-containing protein</td><td align="left">at1g07590</td><td align="left">1.62</td><td align="right">74.6</td></tr><tr><td/><td align="left">Meprin and TRAF domain-containing protein</td><td align="left">at5g26280</td><td align="left">1.60</td><td align="right">70.8</td></tr><tr><td/><td align="left">Expressed protein</td><td align="left">at4g39675</td><td align="left">1.72</td><td align="right">86.0</td></tr><tr><td/><td align="left">C2-domain-containing protein</td><td align="left">at4g15740</td><td align="left">1.76</td><td align="right">42.7</td></tr><tr><td/><td align="left">Expressed protein</td><td align="left">at1g09340</td><td align="left">1.62</td><td align="right">78.4</td></tr><tr><td/><td align="left">Late embryogenesis abundant 3 family protein</td><td align="left">at4g02380</td><td align="left">1.67</td><td align="right">52.1</td></tr><tr><td/><td align="left">Bet v 1 allergen family protein</td><td align="left">at1g24020</td><td align="left">1.54</td><td align="right">98.4</td></tr><tr><td/><td align="left">Universal stress protein</td><td align="left">at3g03270</td><td align="left">1.98</td><td align="right">19.7</td></tr><tr><td align="left">Transporter activity</td><td align="left">Plasma membrane intrinsic protein 2B (PIP2B)</td><td align="left">at2g37170</td><td align="left">1.70</td><td align="right">43.3</td></tr><tr><td/><td align="left">Tonoplast intrinsic protein gamma</td><td align="left">at2g36830</td><td align="left">1.53</td><td align="right">104.2</td></tr><tr><td align="left">Glutathione transferase activity</td><td align="left">Glutathione S-transferase GST20; Tau class</td><td align="left">at1g78370</td><td align="left">1.68</td><td align="right">53.5</td></tr><tr><td align="left">Cell wall structure</td><td align="left">Glycine-rich protein</td><td align="left">at2g05510</td><td align="left">4.31</td><td align="right">2.4</td></tr><tr><td align="left">RNA binding</td><td align="left">Pumilio/Puf RNA-binding domain-containing protein</td><td align="left">At1g78160</td><td align="left">1.60</td><td align="right">67.7</td></tr><tr><td align="left">Stress response</td><td align="left">Drought-responsive protein (Di21)</td><td align="left">at4g15910</td><td align="left">1.67</td><td align="right">69.6</td></tr><tr><td align="left">Electron transport</td><td align="left">Cytochrome B561 family protein</td><td align="left">at5g38360</td><td align="left">1.56</td><td align="right">97.9</td></tr><tr><td align="left">Amino acid biosynthesis</td><td align="left">asparagine synthetase 2</td><td align="left">at5g65010</td><td align="left">1.74</td><td align="right">65.1</td></tr><tr><td align="left">Carbonic anhydrase activity</td><td align="left">Carbonic anhydrase 1, chloroplast</td><td align="left">at3g01500</td><td align="left">1.67</td><td align="right">51.0</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Gene ontology based on molecular function for selected repressed genes of arsenic-treated <italic>Arabidopsis thaliana </italic>Columbia plants.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Molecular function</td><td align="left">Gene ID</td><td align="left">Locus</td><td align="left">FC</td><td align="right">RP</td></tr></thead><tbody><tr><td align="left">Catalase activity</td><td align="left">catalase 3 (SEN2)</td><td align="left">at1g20620</td><td align="left">-1.59</td><td align="right">191.7</td></tr><tr><td align="left">Peroxidase activity</td><td align="left">peroxidase</td><td align="left">at3g49120</td><td align="left">-1.77</td><td align="right">165.3</td></tr><tr><td/><td align="left">peroxidase</td><td align="left">at5g64120</td><td align="left">-1.84</td><td align="right">123.2</td></tr><tr><td/><td align="left">cationic peroxidase</td><td align="left">at4g25980</td><td align="left">-1.52</td><td align="right">333.1</td></tr><tr><td align="left">Oxidoreductase activity</td><td align="left">superoxide dismutase [Fe], chloroplast</td><td align="left">at4g25100</td><td align="left">-5.17</td><td align="right">1.7</td></tr><tr><td/><td align="left">lipoxygenase</td><td align="left">at1g72520</td><td align="left">-2.41</td><td align="right">242.9</td></tr><tr><td/><td align="left">FAD-binding domain-containing protein</td><td align="left">at1g26380</td><td align="left">-1.50</td><td align="right">293.0</td></tr><tr><td/><td align="left">cytochrome p450 83B1</td><td align="left">at4g31500</td><td align="left">-1.71</td><td align="right">142.3</td></tr><tr><td/><td align="left">auxin-responsive family protein</td><td align="left">at5g35735</td><td align="left">-1.59</td><td align="right">156.4</td></tr><tr><td align="left">Metal ion binding</td><td align="left">germin-like protein</td><td align="left">at5g39160</td><td align="left">-1.51</td><td align="right">319.7</td></tr><tr><td/><td align="left">germin-like protein</td><td align="left">at5g39190</td><td align="left">-2.13</td><td align="right">27.3</td></tr><tr><td/><td align="left">calcium-binding EF hand family protein</td><td align="left">at1g76650</td><td align="left">-2.00</td><td align="right">50.4</td></tr><tr><td/><td align="left">C2-domain containing protein</td><td align="left">at4g34150</td><td align="left">-1.53</td><td align="right">313.4</td></tr><tr><td/><td align="left">touch-responsive protein/calmodulin-related</td><td align="left">at2g41100</td><td align="left">-1.64</td><td align="right">137.9</td></tr><tr><td/><td align="left">ferritin 1 (FER 1)</td><td align="left">at5g01600</td><td align="left">-1.78</td><td align="right">83.0</td></tr><tr><td/><td align="left">ferritin 4</td><td align="left">at3g56090</td><td align="left">-1.52</td><td align="right">201.0</td></tr><tr><td/><td align="left">zinc finger (C2H2 type) protein</td><td align="left">at3g46090</td><td align="left">-1.51</td><td align="right">191.2</td></tr><tr><td/><td align="left">zinc finger (C2H2 type) protein</td><td align="left">at3g46080</td><td align="left">-1.59</td><td align="right">178.8</td></tr><tr><td/><td align="left">zinc finger (C3HC4 type) protein</td><td align="left">at5g27420</td><td align="left">-1.75</td><td align="right">82.2</td></tr><tr><td align="left">Hydrolase activity</td><td align="left">lipase class 3 family protein</td><td align="left">at1g02660</td><td align="left">-1.56</td><td align="right">218.8</td></tr><tr><td/><td align="left">invertase/pectin methylesterase family protein</td><td align="left">at5g62360</td><td align="left">-1.75</td><td align="right">85.0</td></tr><tr><td/><td align="left">protein phosphatase 2C</td><td align="left">at2g30020</td><td align="left">-1.52</td><td align="right">237.4</td></tr><tr><td/><td align="left">phosphoric monoester hydrolase</td><td align="left">at1g73010</td><td align="left">-3.01</td><td align="right">7.0</td></tr><tr><td/><td align="left">acid phosphatase type 5 (ACP5)</td><td align="left">at3g17790</td><td align="left">-1.62</td><td align="right">290.7</td></tr><tr><td/><td align="left">phosphoric monoester hydrolase</td><td align="left">at1g17710</td><td align="left">-1.88</td><td align="right">118.6</td></tr><tr><td/><td align="left">glycosyl hydrolase family 17 protein</td><td align="left">at3g55430</td><td align="left">-1.53</td><td align="right">201.2</td></tr><tr><td/><td align="left">glycosyl hydrolase family 17 protein</td><td align="left">at4g31140</td><td align="left">-1.71</td><td align="right">248.5</td></tr><tr><td/><td align="left">glycosyl hydrolase family 17 protein</td><td align="left">at4g19810</td><td align="left">-1.96</td><td align="right">99.2</td></tr><tr><td/><td align="left">glycosyl hydrolase family 36 protein</td><td align="left">at5g20250</td><td align="left">-1.52</td><td align="right">201.8</td></tr><tr><td/><td align="left">xyloglucan endotransglucosylase/hydrolase</td><td align="left">at4g30280</td><td align="left">-1.63</td><td align="right">172.2</td></tr><tr><td/><td align="left">xyloglucan endotransglucosylase/hydrolase</td><td align="left">at4g14130</td><td align="left">-2.00</td><td align="right">47.9</td></tr><tr><td/><td align="left">xyloglucan endotransglucosylase/hydrolase</td><td align="left">at5g57560</td><td align="left">-1.68</td><td align="right">97.3</td></tr><tr><td/><td align="left">nudix hydrolase homolog 4</td><td align="left">at1g18300</td><td align="left">-1.54</td><td align="right">215.9</td></tr><tr><td/><td align="left">MERI-5 endo-xyloglucan transferase</td><td align="left">at4g30270</td><td align="left">-1.96</td><td align="right">45.6</td></tr><tr><td align="left">Protein binding</td><td align="left">calmodulin-binding family protein</td><td align="left">at4g33050</td><td align="left">-1.78</td><td align="right">88.6</td></tr><tr><td/><td align="left">ankyrin repeat family protein</td><td align="left">at5g45110</td><td align="left">-1.58</td><td align="right">252.0</td></tr><tr><td/><td align="left">mitochondrial substrate carrier family protein</td><td align="left">at4g24570</td><td align="left">-1.50</td><td align="right">260.2</td></tr><tr><td/><td align="left">polygalacturonase inhibitory protein</td><td align="left">at5g06860</td><td align="left">-1.91</td><td align="right">65.0</td></tr><tr><td align="left">Chitin binding</td><td align="left">hevein-like protein (HEL)</td><td align="left">at3g04720</td><td align="left">-1.53</td><td align="right">499.9</td></tr><tr><td align="left">Carbohydrate binding</td><td align="left">legume lectin family protein</td><td align="left">at3g16530</td><td align="left">-2.03</td><td align="right">37.8</td></tr><tr><td align="left">Sugar binding</td><td align="left">curculin-like lectin family protein</td><td align="left">at1g78830</td><td align="left">-1.64</td><td align="right">138.3</td></tr><tr><td align="left">ATP binding</td><td align="left">ATP-dependent Clp protease ATP-binding subunit</td><td align="left">at5g51070</td><td align="left">-1.51</td><td align="right">188.7</td></tr><tr><td align="left">Jasmonic acid synthesis</td><td align="left">allene oxide cyclase</td><td align="left">at3g25760</td><td align="left">-1.54</td><td align="right">279.6</td></tr><tr><td align="left">Peptidase activity</td><td align="left">vacuolar processing enzyme gamma</td><td align="left">at4g32940</td><td align="left">-1.65</td><td align="right">117.1</td></tr><tr><td/><td align="left">subtilase family protein</td><td align="left">at1g32970</td><td align="left">-1.91</td><td align="right">81.1</td></tr><tr><td align="left">Ligase activity</td><td align="left">v-box domain-containing protein</td><td align="left">at2g35930</td><td align="left">-1.63</td><td align="right">205.5</td></tr><tr><td/><td align="left">asparagine synthetase 1</td><td align="left">at3g47340</td><td align="left">-1.74</td><td align="right">75.5</td></tr><tr><td align="left">Transferase activity</td><td align="left">glutathione S-transferase (GSTF6); phi class</td><td align="left">at1g02930</td><td align="left">-2.10</td><td align="right">62.9</td></tr><tr><td/><td align="left">glutathione S-transferase (GSTF7); phi class</td><td align="left">at1g02920</td><td align="left">-2.88</td><td align="right">7.6</td></tr><tr><td/><td align="left">branched-chain amino acid amino transferase 2</td><td align="left">at1g10070</td><td align="left">-1.60</td><td align="right">113.8</td></tr><tr><td align="left">Nutrient reservoir activity</td><td align="left">patatin</td><td align="left">at2g26560</td><td align="left">-1.81</td><td align="right">84.8</td></tr><tr><td align="left">Kinase activity</td><td align="left">serine/threonine protein kinase 19</td><td align="left">at3g08720</td><td align="left">-1.55</td><td align="right">249.8</td></tr><tr><td align="left">Molecular function unknown</td><td align="left">hypothetical protein no ATG start</td><td align="left">at3g09922</td><td align="left">-2.15</td><td align="right">33.6</td></tr><tr><td/><td align="left">expressed protein</td><td align="left">at2g25510</td><td align="left">-1.99</td><td align="right">96.9</td></tr><tr><td/><td align="left">expressed protein no ATG start</td><td align="left">at5g03545</td><td align="left">-2.72</td><td align="right">16.8</td></tr><tr><td/><td align="left">expressed protein</td><td align="left">at5g42530</td><td align="left">-2.01</td><td align="right">147.4</td></tr><tr><td/><td align="left">expressed protein</td><td align="left">at4g31570</td><td align="left">-2.21</td><td align="right">45.6</td></tr><tr><td/><td align="left">expressed protein</td><td align="left">at1g69890</td><td align="left">-1.73</td><td align="right">79.9</td></tr><tr><td/><td align="left">expressed protein</td><td align="left">at5g20790</td><td align="left">-2.33</td><td align="right">32.7</td></tr><tr><td/><td align="left">expressed protein</td><td align="left">at2g04460</td><td align="left">-1.73</td><td align="right">190.5</td></tr><tr><td/><td align="left">VQ motif-containing protein</td><td align="left">at2g22880</td><td align="left">-1.83</td><td align="right">57.3</td></tr><tr><td/><td align="left">glycine-rich protein</td><td align="left">at1g07135</td><td align="left">-1.53</td><td align="right">167.4</td></tr><tr><td/><td align="left">glycine-rich protein</td><td align="left">at3g04640</td><td align="left">-1.81</td><td align="right">73.9</td></tr><tr><td/><td align="left">glycine-rich protein</td><td align="left">at2g05540</td><td align="left">-1.85</td><td align="right">47.9</td></tr><tr><td/><td align="left">glycine-rich protein</td><td align="left">at2g05380</td><td align="left">-1.59</td><td align="right">146.7</td></tr><tr><td/><td align="left">integral membrane family protein</td><td align="left">at4g15610</td><td align="left">-1.58</td><td align="right">139.7</td></tr><tr><td/><td align="left">gibberellin-responsive protein</td><td align="left">at1g22690</td><td align="left">-1.88</td><td align="right">51.5</td></tr><tr><td/><td align="left">gibberellin-regulated protein (GASA1)</td><td align="left">at1g75750</td><td align="left">-1.70</td><td align="right">93.7</td></tr><tr><td/><td align="left">dehydrin (RAB18)</td><td align="left">at5g66400</td><td align="left">-1.70</td><td align="right">81.1</td></tr><tr><td/><td align="left">unknown protein – similar to glycosyltransferase</td><td align="left">at2g41640</td><td align="left">-1.54</td><td align="right">185.8</td></tr><tr><td/><td align="left">patatin-like protein 8</td><td align="left">at4g29800</td><td align="left">-1.52</td><td align="right">198.7</td></tr><tr><td/><td align="left">phosphate-responsive protein</td><td align="left">at5g64260</td><td align="left">-1.52</td><td align="right">221.0</td></tr><tr><td/><td align="left">phosphate-responsive protein</td><td align="left">at1g35140</td><td align="left">-1.58</td><td align="right">131.5</td></tr><tr><td/><td align="left">similar to LITAF-domain containing protein</td><td align="left">at5g13190</td><td align="left">-1.54</td><td align="right">293.4</td></tr><tr><td align="left">Transcription factor activity</td><td align="left">DRE-binding protein</td><td align="left">at1g12610</td><td align="left">-1.74</td><td align="right">116.8</td></tr><tr><td/><td align="left">AP2 domain-containing transcription factor</td><td align="left">at4g34410</td><td align="left">-2.01</td><td align="right">50.4</td></tr><tr><td/><td align="left">zinc finger (C2H2 type) protein</td><td align="left">at3g46090</td><td align="left">-1.51</td><td align="right">191.2</td></tr><tr><td/><td align="left">zinc finger (C2H2 type) protein</td><td align="left">at3g46080</td><td align="left">-1.59</td><td align="right">178.8</td></tr><tr><td/><td align="left">WRKY family transcription factor 33</td><td align="left">at2g38470</td><td align="left">-1.63</td><td align="right">157.9</td></tr><tr><td/><td align="left">WRKY family transcription factor 53</td><td align="left">at4g23810</td><td align="left">-1.55</td><td align="right">278.3</td></tr><tr><td/><td align="left">WRKY family transcription factor 40</td><td align="left">at1g80840</td><td align="left">-1.88</td><td align="right">70.9</td></tr><tr><td/><td align="left">NAC domain-containing protein</td><td align="left">at5g08790</td><td align="left">-1.53</td><td align="right">238.7</td></tr><tr><td align="left">Senescence-related</td><td align="left">senescence-associated family protein</td><td align="left">at5g66040</td><td align="left">-1.55</td><td align="right">164.0</td></tr><tr><td/><td align="left">senescence/dehydration-associated protein</td><td align="left">at2g17840</td><td align="left">-1.60</td><td align="right">256.9</td></tr><tr><td/><td align="left">senescence-associated protein (SEN1)</td><td align="left">at4g35770</td><td align="left">-1.59</td><td align="right">123.9</td></tr><tr><td/><td align="left">SRG3 (senescence-related gene 3)</td><td align="left">at3g02040</td><td align="left">-2.65</td><td align="right">11.0</td></tr><tr><td align="left">Transporter activity</td><td align="left">MATE efflux family protein</td><td align="left">at1g61890</td><td align="left">-1.80</td><td align="right">108.3</td></tr><tr><td align="left">Galactolipid biosynthesis</td><td align="left">monogalactosyldiacylglycerol synthase type C</td><td align="left">at2g11810</td><td align="left">-1.78</td><td align="right">130.5</td></tr><tr><td align="left">Electron transport</td><td align="left">cytochrome p450 family 94 subfamily B</td><td align="left">at3g48520</td><td align="left">-1.56</td><td align="right">249.8</td></tr><tr><td align="left">Guanosine tetraphosphate metabolism</td><td align="left">RSH 2 (RELA-SPOT HOMOLOG)</td><td align="left">at3g14050</td><td align="left">-1.52</td><td align="right">188.7</td></tr><tr><td align="left">N-terminal protein myristoylation</td><td align="left">band 7 family protein</td><td align="left">at3g01290</td><td align="left">-1.54</td><td align="right">253.4</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Comparison of microarray expression data (significance criteria of P &lt; 0.001 and FDR of 1%) with RT-PCR data from arsenate-treated <italic>Arabidopsis thaliana</italic>.</p></caption><table frame="hsides" rules="groups"><tbody><tr><td align="left">Gene id</td><td align="left">gene name</td><td align="right">ratio<sup>a</sup></td><td align="center" colspan="2">ΔΔCt<sup>b</sup></td></tr><tr><td/><td/><td/><td align="right"><bold>Day 3</bold></td><td align="right"><bold>Day 10</bold></td></tr><tr><td colspan="1"><hr/></td><td colspan="1"><hr/></td><td colspan="1"><hr/></td><td colspan="2"><hr/></td></tr><tr><td align="left">at2g28190</td><td align="left">Cu Zn SOD (CSD2)</td><td align="right">4.57</td><td align="right">3.06</td><td align="right">3.77</td></tr><tr><td align="left">at1g12520</td><td align="left">Cu Zn SOD Cu chaperone</td><td align="right">3.16</td><td align="right">3.14</td><td align="right">2.44</td></tr><tr><td align="left">at1g08830</td><td align="left">Cu Zn SOD (CSD1)</td><td align="right">2.41</td><td align="right">4.54</td><td align="right">2.96</td></tr><tr><td align="left">at4g25100</td><td align="left">Fe SOD (SODB)</td><td align="right">-5.17</td><td align="right">-2.49</td><td align="right">-1.80</td></tr><tr><td align="left">at1g73010</td><td align="left">phosphoric monoester hydrolase</td><td align="right">-3.01</td><td align="right">-6.70</td><td align="right">-2.97</td></tr><tr><td align="left">at5g20790</td><td align="left">expressed protein</td><td align="right">-2.33</td><td align="right">-4.85</td><td align="right">-2.32</td></tr><tr><td align="left">at2g11810</td><td align="left">MGDG synthase type C</td><td align="right">-1.78</td><td align="right">-1.50</td><td align="right">-2.44</td></tr><tr><td align="left">at3g17790</td><td align="left">acid phosphatase type 5 (ACP5)</td><td align="right">-1.62</td><td align="right">-4.65</td><td align="right">-4.44</td></tr><tr><td align="left">at1g17710</td><td align="left">phosphoric monoester hydrolase</td><td align="right">-1.88</td><td align="right">-9.57</td><td align="right">-3.00</td></tr><tr><td align="left">at2g04460</td><td align="left">expressed protein</td><td align="right">-1.73</td><td align="right">-6.04</td><td align="right">-2.80</td></tr><tr><td align="left">at3g08720</td><td align="left">serine/threonine protein kinase 19</td><td align="right">-1.55</td><td align="right">-5.58</td><td align="right">-1.82</td></tr><tr><td align="left">at3g02040</td><td align="left">senescence-related gene 3 (SRG3)</td><td align="right">-2.65</td><td align="right">-5.12</td><td align="right">-1.85</td></tr><tr><td align="left">at5g61650</td><td align="left">P-type cyclin</td><td align="right">1.64</td><td align="right">-5.60</td><td align="right">-3.39</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>Functional characterization of differentially expressed <italic>A. thaliana </italic>genes in response to As (V) stress.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>Microarray quality control for chips used in this study.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p>Complete microarray dataset for genes induced by As (V).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional file 4</title><p>Complete microarray dataset for genes repressed by As (V).</p></caption></supplementary-material>
<table-wrap-foot><p>RP, rank products; FC, fold change. Only genes upregulated above 1.5-fold and meeting a significance criteria of P &lt; 0.001 and FDR of 1% are shown. Rank products analysis reveals that most significantly induced genes display the lowest RP value.</p></table-wrap-foot> <table-wrap-foot><p>RP, rank products; FC, fold change. Only genes upregulated above -1.5-fold and meeting a significance criteria of P &lt; 0.001 and FDR of 1% are shown. Rank products analysis reveals that most significantly repressed genes display the lowest RP value.</p></table-wrap-foot> <table-wrap-foot><p><sup>a </sup>Microarray data (linear fold-change)</p><p><sup>b </sup>RT-PCR data (ΔΔCt values) Arabidopsis Actin was used as a reference gene (see Materials and Methods)</p><p><italic>P </italic>values for RT-PCR data were &lt; 0.05.</p></table-wrap-foot>
<graphic xlink:href="1471-2229-8-87-1"/> <graphic xlink:href="1471-2229-8-87-2"/> <graphic xlink:href="1471-2229-8-87-3"/>
<media xlink:href="1471-2229-8-87-S1.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2229-8-87-S2.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2229-8-87-S3.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2229-8-87-S4.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Moore", "Fleishcher", "Woolson", "Grossblatt N"], "given-names": ["L", "M", "E"], "article-title": ["Distribution of arsenic in the Environment"], "source": ["Medical and Biologic Effects of Environmental Pollutants: Arsenic"], "year": ["1977"], "publisher-name": ["Washington, D.C.: National Academy of Sciences"], "fpage": ["16"], "lpage": ["26"]}, {"collab": ["International Agency for Research on Cancer"], "article-title": ["Monograph of the evaluation of carcinogenic risk to humans \u2013 Overall evaluation of carcinogenicity"], "source": ["An update of IARC monographs 1 to 42 Lyon"]}, {"collab": ["National Research Council"], "source": ["Arsenic in drinking water"], "year": ["1999"], "publisher-name": ["National Academy Press: Washington, D.C"]}, {"surname": ["Ullrich-Eberius", "Sanz", "Novacky"], "given-names": ["C", "A", "A"], "article-title": ["Evaluation of arsenate- and vanadate-associated changes of electrical membrane potential and phosphate transport in "], "italic": ["Lemna gibba"], "source": ["J Exp Bot"], "year": ["1989"], "volume": ["40"], "fpage": ["119"], "lpage": ["128"]}, {"surname": ["Meharg", "Hartley-Whitaker"], "given-names": ["A", "J"], "article-title": ["Arsenic uptake and metabolism in arsenic-resistant and non-resistant plant species"], "source": ["New Phytologist"], "year": ["2002"], "volume": ["154"], "fpage": ["29"], "lpage": ["43"]}, {"surname": ["Schmoger", "Oven", "Grill"], "given-names": ["M", "M", "E"], "article-title": ["Detoxification of arsenic by phytochelatins in plants"], "source": ["Plant Phys"], "year": ["2000"], "volume": ["122"], "fpage": ["793"], "lpage": ["801"]}, {"surname": ["Benjamini", "Hochberg"], "given-names": ["Y", "Y"], "article-title": ["Controlling the false discovery rate: A practical and powerful approach to multiple testing"], "source": ["J Royal Stat Soc"], "year": ["1995"], "volume": ["57B"], "fpage": ["289"], "lpage": ["300"]}, {"surname": ["Azevedo", "Alas", "Smith", "Lea"], "given-names": ["RA", "RM", "RJ", "PJ"], "article-title": ["Response of antioxidant enzymes to transfer from elevated carbon dioxide to air and ozone fumigation, in the leaves and roots of wild-type and a catalase-deficient mutant of barley"], "source": ["Physiol Plant"], "year": ["1998"], "volume": ["104"], "fpage": ["280"], "lpage": ["292"]}]
BMC Plant Biol. 2008 Aug 6; 8:87