- Research article
- Open Access
Comparative transcriptomic characterization of aluminum, sodium chloride, cadmium and copper rhizotoxicities in Arabidopsis thaliana
© Zhao et al; licensee BioMed Central Ltd. 2009
- Received: 06 October 2008
- Accepted: 23 March 2009
- Published: 23 March 2009
Rhizotoxic ions in problem soils inhibit nutrient and water acquisition by roots, which in turn leads to reduced crop yields. Previous studies on the effects of rhizotoxic ions on root growth and physiological functions suggested that some mechanisms were common to all rhizotoxins, while others were more specific. To understand this complex system, we performed comparative transcriptomic analysis with various rhizotoxic ions, followed by bioinformatics analysis, in the model plant Arabidopsis thaliana.
Roots of Arabidopsis were treated with the major rhizotoxic stressors, aluminum (Al) ions, cadmium (Cd) ions, copper (Cu) ions and sodium (NaCl) chloride, and the gene expression responses were analyzed by DNA array technology. The top 2.5% of genes whose expression was most increased by each stressor were compared with identify common and specific gene expression responses induced by these stressors. A number of genes encoding glutathione-S-transferases, peroxidases, Ca-binding proteins and a trehalose-synthesizing enzyme were induced by all stressors. In contrast, gene ontological categorization identified sets of genes uniquely induced by each stressor, with distinct patterns of biological processes and molecular function. These contained known resistance genes for each stressor, such as AtALMT1 (encoding Al-activated malate transporter) in the Al-specific group and DREB (encoding dehydration responsive element binding protein) in the NaCl-specific group. These gene groups are likely to reflect the common and differential cellular responses and the induction of defense systems in response to each ion. We also identified co-expressed gene groups specific to rhizotoxic ions, which might aid further detailed investigation of the response mechanisms.
In order to understand the complex responses of roots to rhizotoxic ions, we performed comparative transcriptomic analysis followed by bioinformatics characterization. Our analyses revealed that both general and specific genes were induced in Arabidopsis roots exposed to various rhizotoxic ions. Several defense systems, such as the production of reactive oxygen species and disturbance of Ca homeostasis, were triggered by all stressors, while specific defense genes were also induced by individual stressors. Similar studies in different plant species could help to clarify the resistance mechanisms at the molecular level to provide information that can be utilized for marker-assisted selection.
- Gene Ontology
- Reactive Oxygen Species Production
- Gene Group
- Gene Expression Response
- Relative Fold Change
Poor root growth is caused by various rhizotoxic factors present in problem soils, and is linked to susceptibility to other stress factors. For example, aluminum (Al) ions cause severe damage to the roots of plants growing in acid soil, accentuating nutrient deficiency and increasing their sensitivity to drought stress . Other metal rhizotoxins, such as cadmium (Cd) and copper (Cu) ions, also inhibit root growth . The poor development of roots occurs because Al, sodium (Na) and Cu ions have negative impacts on the shoot yield of crop plants in problem soils, while Cd ions decrease the efficiency of phytoremediation in Cd-contaminated soils. Improving the tolerance of roots to rhizotoxic ions is therefore an important target in plant breeding. Understanding of the molecular responses of plants to rhizotoxic ions is a critical step towards molecular breeding of stress tolerant crops using marker-assisted selection or genetic engineering.
Several critical genes regulating tolerance to rhizotoxic ions have been identified in studies using hypersensitive mutants. Studies with salt overly sensitive (SOS) mutants identified genes encoding proteins critical for salt sensitivity, including the Na+/H+ antiporter (SOS1)  and its regulating protein kinase, SOS2 . Using the Cd- and Al-sensitive mutants, cad and als, revealed that genes for phytochelatin synthase (CAD1)  and a putative ATP-binding Al-translocator (ALS3)  were involved in tolerance mechanisms to these ions. The identification of stress-responsive genes is a useful approach, because some stress-inducible genes might also be involved in tolerance mechanisms associated with abiotic rhizotoxins. For example, the cis-element DRE , and its binding protein DREB, were identified from a series of studies on dehydration-inducible genes. Several Al-tolerant genes are also responsive to Al ions, such as ALS3 , GST  and AtALMT1 . Analyses of those genes that are responsive to individual rhizotoxic treatments could also improve our knowledge of the mechanisms of toxicity of the different ions.
Genome-wide transcript analysis can be performed in Arabidopsis and other plant species using commercially available oligo-microarray techniques. These techniques have recently been applied to the identification of rhizotoxin-responsive genes in Arabidopsis (e.g. NaCl  and Al ) and other plant species (e.g. Al in maize [12, 13] and Medicago ). Those studies demonstrated that various genes were induced by each rhizotoxin. In order to understand the functions and impacts of such gene expression responses to each rhizotoxin, it is important to distinguish those genes induced as part of the general stress responses from those specific to individual stressors. The comparison of transcriptomes among different treatments and the application of bioinformatics procedures (e.g. co-expression gene analysis) are potentially useful approaches for determining the characteristics of these different gene groups.
In order to determine the effects of rhizotoxic treatments on gene expression in Arabidopsis using this microarray approach, it is necessary to minimize the effects of other factors on gene expression during the course of the experiment. For example, mechanical damage to the roots triggers the expression of "general" stress-responsive genes , and may lead to false conclusions if such a "general response" is not involved in each stress treatment. We previously developed a hydroponic culture system that enhanced rhizotoxicity while minimizing mechanical damage when changing culture solutions [16, 17]. This method has been applied to quantitative trait locus analysis of rhizotoxicities  and for monitoring root tip viability , suggesting that it would also be suitable for obtaining root samples to determine the direct effects of rhizotoxins using microarray analyses. We have also developed an RNA extraction method for Arabidopsis that allows the isolation of high quality RNA from various tissues, including roots, at different developmental stages . This can be adapted to rhizotoxin-damaged roots, allowing the isolation of RNA of sufficiently high quality to allow the determination of the complex patterns of gene expression in response to rhizotoxins, using DNA microarray technology.
In the present study, we combined these experimental procedures to analyze gene expression responses in roots by microarray analysis, following treatment with Al, Cu and Cd ions, or NaCl. By comparing microarray data, we were able to separate the general (i.e. common to all rhizotoxic ions) and specific (i.e. more specific to each ion) gene expression responses that were induced by each rhizotoxic ion. Analyses of the separated gene groups based on Arabidopsis gene information and bioinformatics tools revealed that both general and individual toxic mechanisms and defense responses were triggered by each rhizotoxic ion.
Identification of genes responsive to all ions and to individual rhizotoxic ions
Classification by GO categories defined by TAIR for whole genome genes and for gene groups upregulated by rhizotoxic ions identified by a comparative microarray approach.
Proportion of genes among GO categories (%)
GO slim category
DNA or RNA binding
nucleic acid binding
other enzyme activity
other molecular functions
receptor binding or activity
structural molecule activity
transcription factor activity
unknown molecular functions
Characteristics of genes induced by all ions
Forty-one genes were identified that responded to all the tested ions (Figure 3; Additional file 3). This group contained a significantly larger percentage of genes with "other binding" activity by GO categorization of molecular function (Table 1), including six Ca-binding proteins, such as calmodulin-like proteins (CML38 and 37/39) and Ca-binding EF hand proteins, which were rare in other gene groups (Additional file 4). Three disease resistance proteins, one belonging to the TIR (Toll-Interleukin-Resistance) class of proteins with molecular transducer activities, were also included in this group, which was previously identified as one of the typical stress responsive genes. The group also contained typical reactive oxygen species (ROS)-responsive genes that encoded ROS-scavenging enzymes (three glutathione transferases and two peroxidases), as well as those involved in the signal transduction pathway for ROS responses, namely MYB15 and tolB-related protein. A putative trehalose-phosphate phosphatase gene belonged to this gene group and might be related to the reduction of cellular damage from ROS via the accumulation of trehalose. Induction of these genes could account for the results of previous physiological studies, which reported that ROS production and Ca-alleviation were common features of various rhizotoxicities.
Characteristics of genes uniquely induced by individual ions
1. Genes uniquely induced by Al ions
The Al-responsive group consisted of 103 genes (Figure 3), and included a significantly higher percentage of genes encoding proteins with transporter (10.7%) and transferase (16.5%) activities, by GO categorization of molecular function. Genes encoding transporters were concentrated (i.e. about 19%) in a gene cluster containing 32 genes (Figure 4A), which were relatively specific to Al ions (Table 1). Major transporters for sulfate (SULTR3;1) and borate (BOR2) were found in this specific cluster, together with AtALMT1 and other organic molecule transporters [e.g. mannitol and the organic cation/carnitine transporter (AtOCT1)]. This specific gene cluster also contained genes encoding an auxin/Al-responsive protein, an auxin carrier protein, and a gene encoding purple acid phosphatase.
Although genes encoding transferases were not concentrated in a specific gene cluster, two S-adenosyl-L-methionine:carboxyl methyltransferase family proteins and three carbohydrate transferases (e.g. glycosyltransferase) belonged to this gene group. A large number of genes involved in carbon and nitrogen metabolism were also identified in this Al-specific group, including glutamate dehydrogenase (GDH2), malic enzymes (AtNADP-ME1 and 2), and some carbohydrate decarboxylases, including a pyruvate decarboxylase (Additional file 3).
2. Genes uniquely induced by NaCl
Venn diagram analysis identified 57 genes that were uniquely induced by NaCl treatment (Figure 3). GO analysis for molecular function suggested that this gene group contained a significantly higher percentage of genes encoding transcription factors (24.6%) (Table 1), while GO analysis for biological process found a higher percentage of genes in the transcription category (Figure 2). This group contained more transcription factors, including some DREB family proteins (three of a total of six DREB families identified in all gene groups), which have been recognized as playing a role in salt tolerance. Cluster analysis revealed that 22 genes, including seven transcription factors, were more specific to NaCl than were the other genes (Figure 4B). Some cold-responsive genes (e.g. COR6.6, COR78), whose signal transduction pathways overlap with NaCl stress, were also identified in this cluster. No genes for major catalytic enzymes involved in carbon or nitrogen metabolism, and only one transporter, were found in the NaCl group.
3. Genes uniquely induced by Cd ions
The Cd ion-induced gene group contained no catalytic enzymes involved in major primary or secondary metabolism, but did include some protein kinases, such as receptor-like protein kinases (CRK6 and 10)  (Additional file 3). This could account for the significantly higher ratio of genes with "kinase activity" (16.7%), when genes were categorized by molecular function (Table 1). Genes belonging to the enriched GO category were not enriched in the specific gene cluster (Figure 4C). The specific cluster, also, contained several stress-responsive genes, whose functions such as heat-shock and defense-response, have not yet been clarified (Figure 4C). One gene categorized by GO as having kinase activity, a leucine-rich repeat family protein (AtRLP38) similar to disease resistant proteins, was also identified in this specific gene cluster.
4. Genes uniquely induced by Cu ions
The Cu ion group contained known Cu-detoxifying and binding molecules, such as metallothionein (MT2A) (Additional file 3). A large number of secondary metabolite-synthesizing enzymes involved in "other metabolic processes" (Figure 2), such as strictosidine synthase 3 (SS3) (involved in alkaloid synthesis), anthranilate synthase and six isoforms of cytochrome P450 were also identified in this group. These could account for the significantly higher percentages of genes encoding proteins with other enzyme activities (23.4%) and transferase activities (19.5%) (Table 1). Two trehalose synthases (ATTPS8 and 11) and a ROS-scavenging protein, namely thioredoxin H-8 (ATH-8), may reflect the relative severity of ROS production induced by Cu ion treatment, compared with the other ions. In the Cu ion-specific gene cluster, an l-aminocyclopropane-1-carboxylate synthase (ACC synthase; ACS8) belonging to the ethylene biosynthesis pathway was identified, together with an enzyme relating to auxin synthesis [i.e. an indoleacetic acid (IAA) amide synthase (AUR3)]. An enzyme synthesizing the precursor of IAA, tryptophan, namely tryptophan synthase alpha chain (TSA1) and beta chain (TSB1), were identified in the Cu ion-responsive gene group.
Root tip viability, cell damage and ROS production following rhizotoxic treatments
Co-expression gene analysis within each group
The group of genes that was responsive to all ions contained a large number of genes encoding ROS-scavenging enzymes, such as glutathione transferase and peroxidases, and an enzyme for producing trehalose, whose accumulation stabilizes cellular structure against ROS damage . Overexpression of these genes conferred abiotic stress tolerance [26–28] and their induction would therefore act as part of the defense responses against ROS damage induced by all stressors. One large cluster of genes in this group, identified by co-expression gene analysis by KAGIANA search (see Methods), contained various Ca-binding proteins, including previously identified calmodulin-like proteins (CML37/39 and 38), which were inducible by various stimuli , suggesting that Ca-mediated signaling pathways could play important roles in the Arabidopsis response to rhizotoxic stressors. A previously identified transcription factor MYB15, which is involved in the cold stress-mediated defense system associated with ICE1 (inducer of CBF expression) , was also included in this cluster. It seems likely that this gene group, which was responsive to all ions, is related to the general stress-responsive system in plants. Although the pattern of staining was different, ROS accumulation occurred in the roots subjected to milder rhizotoxic conditions (i.e. concentrations causing 50% growth inhibition) (Additional file 7). This suggests that ROS production is a general feature of rhizotoxic treatments. It is interesting to note that the induction of ROS-scavenging enzymes was common to all stressors, and occurred even under mild stress conditions.
The gene groups responsive to individual ions included those genes typically upregulated by each stressor. For example,AtALMT1 was highly upregulated by Al ions, but was not responsive to other ions . In addition, the upregulation of this gene was the largest detected among all the genes (Additional file 3), suggesting that it plays a critical role in the active Al ion defense system of this plant species . Interestingly, the bypass pathways of tricarboxylic acid and glutamate metabolism were also relatively upregulated by Al ion treatment, compared with other treatments (Additional file 3). This could be related to malate efflux, because organic acid excretion can be enhanced by transgenic modification of several enzymes involved in tricarboxylic acid metabolism and its bypass (e.g. citrate synthase ), though the regulation of cytosolic pH caused by changes in these bypass pathways is a possible alternative mechanism. These possibilities need to be tested by future research.
Other rhizotoxic ions, namely NaCl, Cd and Cu ions, induced distinct and specific sets of genes (Figure 3, Additional file 3). For example, gene clusters in the Cu ion-responsive group consisted of senescence-responsive genes, including a gene encoding a previously identified senescence related protein (AtSEN1), which enhances mRNA degradation . This may be related to the stimulation of secondary metabolism pathways, such as terpenoid indole alkaloid metabolism, involving strictosidine synthases (SS2; ) and tryptophan synthases (TSB2 ; TSA1 ), which are activated in mature and senescent tissues. The Cu ion-responsive group also contained various defensive genes, such as ATTPS8 and 11 , which are involved in trehalose synthesis, in addition to the well-characterized Cu ion-detoxifying protein metallothionein (MT2), indicating enhancement of ROS-scavenging capacity. Cu treatment also stimulated thioredoxin gene expression (thioredoxin H-8 (ATH-8) ), which is involved in the Cu ion tolerance mechanism of some organisms . The Salmonella thioredoxin homolog possibly acts by reducing free Cu ions through regulating the binding capacity of the reduced form of thioredoxin to Cu ions . Taken together, the Cu ion group contained genes reflecting the toxicity of Cu and defensive genes that produced proteins to alleviate Cu toxicity.
The other uniquely identified gene groups had similar compositions. The NaCl group demonstrated the importance of the DREB system in defense . Although previous studies have reported that the DREB1A family was responsive to cold treatment, but not to Na ions , our data indicate that this family is also involved in the NaCl-responsive system in the root. This discrepancy might be because of differences in strength of the NaCl used, as our treatment was almost five times milder (50 mM) than that used in previous molecular biological studies (e.g. ). On the other hand, the Cd ion-responsive gene group consisted of unidentified stress-responsive proteins, which were categorized as heat shock and pathogen-related proteins. Further research is needed to clarify the role of these proteins in Cd tolerance.
When we applied the same experimental design using the lower 2.5 percentile as the threshold, we are able to characterize the groups of genes downregulated by each stressor (Additional files 8, 9). GO annotation by molecular function (Additional file 10) showed that uniquely identified groups of genes had distinct patterns. For example, genes with "hydrase activity" were increased by NaCl, or Cd and Cu ions, while those for "transporters" were increased by Al ion treatment (Additional file 9). On the other hand, several genes relating to defense responses, such as disease-resistance related protein, were found to be downregulated by all stressors. This suggests that a combination of up-regulation and downregulation of stress responsive genes may be important in optimizing the adaptation of particular biological pathways to stress conditions.
Co-expressed gene clusters may reflect the cellular conditions and activated defense systems induced by each stressor. For example, Al ions induce phosphate deficiency as a secondary effect , while defense systems for abiotic stressors are activated by phytohormones (e.g. ABA in Cd and Na tolerance ). Based on the upregulations recorded by GENEVESTIGATOR , we may infer that the ABA signaling pathway was activated by both Cu and Al treatments, because a large portion of one cluster in both the Cu ion- (6/7 in the upper cluster; Figure 6E) and Al ion- (3/4 in the middle cluster; Figure 6B) responsive groups consisted of ABA-responsive genes. Furthermore, activation of the salicylic acid signaling pathway was involved in the responses to all treatments, because a cluster responsive to salicylic acid was identified in the shared gene group. These results could explain the involvement of these signaling pathways in the tolerance mechanisms for each stressor (e.g. ABA signal in Al  and Cu tolerance ; salicylic acid signal in Al , NaCl , Cd  and Cu tolerance ).
To investigate the changes in gene expression caused by various rhizotoxic ions, we employed a simple experimental design using a limited number of microarrays (i.e. single time point and single treatment for each ion). This could be advantageous in terms of experimental costs when applying a similar approach to other plant species. Accurate information (e.g. GO) provided by recent developments in the functional genomics of Arabidopsis, is critically important for the success of this approach. Similar developments in genomic research are becoming available for other plant species, and we can therefore apply this procedure to other plant species, and can use comparative genomics to compare the resistance (and damage) systems to rhizotoxic ions among different plant species. Integrated analyses with other -omics data (e.g. metabolomics) would also be interesting to further our understanding of tolerance to and toxicity of rhizotoxic stressors.
There are limitations to our current approach, and several questions remain. For example, we focused on the genes upregulated either collectively or specifically by four different ions. This method excluded genes that were upregulated by two or three stressors, though they may also play an important role in defense and stress-response. For example, some genes encoding cell wall-associated proteins and vacuole loading proteins, which are known to be involved in Cd and Al tolerance, were excluded by our approach. On the other hand, we selected upregulated genes using the upper 2.5 percentile as a threshold. This relative threshold value was preferable to using an absolute fold change threshold value, allowing the selection of a similar number of genes from each treatment group, despite variable distributions of fold changes. This allowed comparison among the groups of genes with similar weights of importance. However, our procedure cut-off the genes if their fold change values were just below the upper 2.5 percentile. The impact of these genes would therefore have been underestimated by the present analysis. Further investigation of these genes using a different method of data analysis is required for a complete understanding of the complex nature of rhizotoxicities.
Using genome-wide DNA microarray technology, we analyzed the impact of rhizotoxic ions (Al, Cd and Cu) and NaCl on gene expression in the roots of Arabidopsis. Comparison of the microarray data allowed the induced genes to be grouped into those common to all treatments, and those unique to individual treatments. Each gene group contained reported tolerance genes, such as AtALMT1 in Al treatment, DREB in NaCl treatment and MT2 in Cu treatment. ROS-scavenging enzymes and Ca-binding proteins, however, were in the group of genes that was upregulated by all stressors. These results were consistent with tolerance mechanisms identified in previous physiological studies. In addition, bioinformatics analyses of the genes groups showed that distinct physiological responses were induced by each stressor. Overall, we showed that comparative microarray analysis with a simple experimental design was a useful technique for identifying gene responses that were consistent with the cytotoxic and tolerance mechanisms of the roots to rhizotoxic ions. Further data analysis, such as a comparison of downregulated genes and the integration of other -omics based technologies (e.g. metabolomics) would be useful for further research into the complex nature of the responses of plant roots to rhizotoxic stressors. Recent developments in genomic research in other plant species may allow us to use similar approaches in various plant species, allowing useful comparisons to be made regarding the similarities and differences in the tolerance mechanisms among different plant species.
Accession Col-4 (N933, NASC; Nottingham Arabidopsis Stock Center) seedlings were pre-grown for 10 days in modified MGRL medium (pH 5.0), as described previously  (2% MGRL nutrients, but the with the Ca ion concentration adjusted to 200 μM) The culture apparatus was made from plastic and resin fibers, which were stable in the presence of the chemical rhizotoxins  (Additional file 1). About 200 seedlings were grown in each culture apparatus. At day 10, the seedlings were transferred, together with the culture apparatus, to 500 ml of modified MGRL medium, containing either 25 μM AlCl3 (pH 4.95), 50 mM NaCl (pH 5.0), 15 μM CdCl2 (pH 5.0) or 1.6 μM CuSO4 (pH 5.0). In all cases, appropriate stock solutions were used to minimize precipitation and change of ionic form. Two liters of each toxic solution was used for the treatment of 200 seedlings. Under these conditions, the pH of the solution was stable (ΔpH < 0.03) and the toxic ions remained soluble (i.e. concentration in the supernatant fractions obtained by centrifugation at 20,000 g, for 15 min) after 24 h incubation. Gene expression could be influenced by secondary effects (e.g. apoptosis or necrosis) if the treatment was too severe, whilst too weak treatment might not trigger the expression of some of the rhizotoxin-sensitive genes. As a compromise, rhizotoxic treatments were carried out using concentrations of Al, Cd, and Cu ions, and NaCl that caused approximately 90% growth inhibition during a 1-week growth test (data not shown). Control pre-grown seedlings were transferred to the basal test solution on day 10 (no stress). Room temperature was maintained at 23–25°C and illumination was controlled at 12 h daytime (30 μmol E m-2 s-1)/night time (no illumination) cycles during pre-growth, and continuous illumination during rhizotoxic treatments. After treatment with rhizotoxins for 24 h, seedlings were removed from the apparatus using forceps. Roots were then rinsed in distilled water and excess water was removed by absorption with tissue. Roots were excised with scissors, immediately frozen in liquid N2 and stored at -80°C in plastic sample tubes (5 ml) until use (Additional file 1A–c, d).
Total RNA was isolated using the method described by Suzuki et al.  then quantified at A260, using a NanoDrop spectrophotometer (ND-1000, NanoDrop Technologies, Wilmington, DE, USA). The quality of RNA used for microarray analysis was measured using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA), according to the manufacturer's instructions.
Microarray analyses were carried out using a competitive hybridization method (i.e. dye-flip method) using the Agilent microarray system (Agilent Technologies, Palo Alto, CA, USA). All procedures were carried out according to the manufacturer's protocols. Briefly, 1 μg of total RNA from each sample was used to synthesize cRNA and was labeled with cyanine-5 (Cy5)- or cyanine-3 (Cy3)-labeled CTP (Perkin Elmer/NEN Life Sciences, Tokyo, Japan). The labeled cRNAs (including a treatment and a control sample) were competitively hybridized to the Agilent Arabidopsis 2 Oligo Microarray, and then washed. The hybridized slides were scanned using Agilent DNA Microarray Scanner (Software Version 6.1) and data points were extracted using Agilent Feature Extraction software (Version 8.1). Three comparisons, including one dye-flip, were made between biologically independent samples. All microarray data have been deposited in a public database (we will complete this before final acceptance). Expression profile trends for the selected genes were examined by reverse transcription-polymerase chain reaction (RT-PCR) using specific primers. Expression patterns judged by RT-PCR were similar to those examined by microarray analysis (Additional file 11). All microarray data are available through the ARRAYEXPRESS database http://www.ebi.ac.uk/microarray-as/ae/ with accession code of E-MEXP-1907 (Transcription profiling of Al, Cu, Cd, NaCl, stress).
Statistical analyses of microarray data and drawing of scatter plots were performed using GeneSpring GX 7.3 (Silicon Genetics, Redwood City, CA, USA), while identification of GO and classification were carried out using software (e.g. http://www.arabidopsis.org/tools/bulk/go/index.jsp) available from the Arabidopsis Information Resource (TAIR) database http://www.arabidopsis.org. Data points based on less than three measurements and spots with low fluorescence intensity (i.e. 80 <) were excluded from the analyses. After this quality control, 15,523 genes were selected from those spotted on the Agilent Arabidopsis 2 microarray slides. Genes were considered to be highly upregulated by a stressor if their fold change values were in the upper 2.5% (i.e. upper 388 genes, based on fold change values) in all three replications. Based on these criteria, 233, 181, 221 and 245 genes were identified as being highly upregulated by Al, Cd or Cu ions, or NaCl, respectively. The means and standard errors (SE) of log2 (fold changes) were calculated for each data point using GeneSpring GX 7.3. Genes upregulated by all stressors and those uniquely responding to each stressor (designated as unique gene groups) were identified by Venn diagrams.
Cluster analyses were performed using CLUSTER software  (available at http://rana.lbl.gov/eisen/) to group the genes in the unique gene groups by specificity of their responses to a particular stressor. The mean FC with other ion treatments was divided by the mean FC with a particular stressor (e.g. FCs for NaCl, Cd and Cu treatments in the Al unique group were divided by the FC for Al treatment), and were designated as "relative fold change" (RFC). The RFCs of the unique gene groups were separately introduced to the Cluster software and then grouped with hierarchical clustering using the average linkage clustering method. Data were normalized to the mean RFC for each treatment. The output data were visualized using the TREEVIEW program. Under these conditions, genes that were highly specific to a particular stressor showed as greenish in color. Genes were manually sub-grouped by the formed cluster and its color, and specificity was assessed by comparison of sub-groups using the Scheffe test (P < 0.05).
Co-expression gene analysis was performed using KAGIANA software http://pmnedo.kazusa.or.jp/kagiana/ using a co-expression gene file obtained from ATTED-II , which consists of 58 experiments including 1388 arrays, and was downloaded from TAIR. Co-expression gene clusters were connected by lines if their Pearson's coefficients of correlation for gene expression were > 0.6. Expression profiles of the clustered genes to various treatments were collected from the response viewer of the GENEVESTIGATOR database (https://www.genevestigator.ethz.ch/at/; see Additional file 6).
Cell viability, plasma membrane integrity and accumulation of ROS (H2O2 and O2-) were visualized using fluorescent probes, as described previously [19, 52]. Briefly, growing roots (day 5) were transferred to stress conditions identical to those used for microarray experiments. Cell viability was determined by FDA staining (5 μg/ml for 30 s), which generates green fluorescence in viable cells. Plasma membrane damage was visualized by PI (3 μg/ml for 1 min), while H2O2 and O2- accumulation were visualized using H2DCFDA (10 μM for 10 min) and DHE) (10 μM for 30 min, 37°C), respectively. Fluorescence in the root tip was observed using a fluorescence microscope (IMT-2-21-RFL, Olympus, Tokyo) equipped with appropriate dichroic mirror units (PI, IMT-2-DMG; FDA, H2DCFDA; DHE, IMT-2-DMIB). Images were photographed using a digital camera unit (PMDC α/OL-1, Olympus). GO was searched for on TAIR using a web tool for GO annotations and categorization .
This work was supported by the Japan Society for the Promotion of Science and the Ministry of Economy for HK, Trade and Industry, Japan for HK and DS. We thank for the Nottingham Arabidopsis Stock Center providing Arabidopsis seeds.
- Kochian LV, Hoekenga OA, Pineros MA: How do crop plants tolerate acid soils? Mechanisms of aluminum tolerance and phosphorous efficiency. Annu Rev Plant Biol. 2004, 55: 459-493.PubMedView ArticleGoogle Scholar
- Kinraide TB, Pedler JF, Parker DR: Relative effectiveness of calcium and magnesium in the alleviation of rhizotoxicity in wheat induced by copper, zinc, aluminum, sodium, and low pH. Plant Soil. 2004, 259: 201-208.View ArticleGoogle Scholar
- Wu SJ, Ding L, Zhu JK: SOS1, a genetic locus essential for salt tolerance and potassium acquisition. Plant Cell. 1996, 8: 617-627.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhu JK, Liu J, Xiong L: Genetic analysis of salt tolerance in Arabidopsis. Evidence for a critical role of potassium nutrition. Plant Cell. 1998, 10: 1181-1191.PubMedPubMed CentralView ArticleGoogle Scholar
- Howden R, Cobbett CS: Cadmium-Sensitive Mutants of Arabidopsis thaliana. Plant Physiol. 1992, 100: 100-107.PubMedPubMed CentralView ArticleGoogle Scholar
- Larsen PB, Geisler MJ, Jones CA, Williams KM, Cancel JD: ALS3 encodes a phloem-localized ABC transporter-like protein that is required for aluminum tolerance in Arabidopsis. Plant J. 2005, 41: 353-363.PubMedView ArticleGoogle Scholar
- Yamaguchi-Shinozaki K, Shinozaki K: A novel cis-acting element in an Arabidopsis gene is involved in responsiveness to drought, low-temperature, or high-salt stress. Plant Cell. 1994, 6: 251-264.PubMedPubMed CentralView ArticleGoogle Scholar
- Ezaki B, Gardner RC, Ezaki Y, Matsumoto H: Expression of aluminum-induced genes in transgenic Arabidopsis plants can ameliorate aluminum stress and/or oxidative stress. Plant Physiol. 2000, 122: 657-665.PubMedPubMed CentralView ArticleGoogle Scholar
- Hoekenga OA, Maron LG, Pineros MA, Cancado GM, Shaff J, Kobayashi Y, Ryan PR, Dong B, Delhaize E, Sasaki T, Matsumoto H, Yamamoto Y, Koyama H, Kochian LV: AtALMT1, which encodes a malate transporter, is identified as one of several genes critical for aluminum tolerance in Arabidopsis. Proc Natl Acad Sci USA. 2006, 103: 9738-9743.PubMedPubMed CentralView ArticleGoogle Scholar
- Jiang Y, Deyholos MK: Comprehensive transcriptional profiling of NaCl-stressed Arabidopsis roots reveals novel classes of responsive genes. BMC Plant Biol. 2006, 6: 25PubMedPubMed CentralView ArticleGoogle Scholar
- Kumari M, Taylor GJ, Deyholos MK: Transcriptomic responses to aluminum stress in roots of Arabidopsis thaliana. Mol Genet Genomics. 2008, 279: 339-357.PubMedView ArticleGoogle Scholar
- Guo P, Bai G, Carver B, Li R, Bernardo A, Baum M: Transcriptional analysis between two wheat near-isogenic lines contrasting in aluminum tolerance under aluminum stress. Mol Genet Genomics. 2007, 277: 1-12.PubMedView ArticleGoogle Scholar
- Maron LG, Kirst M, Mao C, Milner MJ, Menossi M, Kochian LV: Transcriptional profiling of aluminum toxicity and tolerance responses in maize roots. New Phytol. 2008, 179: 116-128.PubMedView ArticleGoogle Scholar
- Chandran D, Sharopova N, VandenBosch KA, Garvin DF, Samac DA: Physiological and molecular characterization of aluminum resistance in Medicago truncatula. BMC Plant Biol. 2008, 8: 89PubMedPubMed CentralView ArticleGoogle Scholar
- Reymond P, Weber H, Damond M, Farmer EE: Differential gene expression in response to mechanical wounding and insect feeding in Arabidopsis. Plant Cell. 2000, 12: 707-720.PubMedPubMed CentralView ArticleGoogle Scholar
- Toda T, Koyama H, Hara T: A simple hydroponic culture method for the development of a highly viable root system in Arabidopsis thaliana. Biosci Biotechnol Biochem. 1999, 63: 210-212.PubMedView ArticleGoogle Scholar
- Koyama H, Toda T, Hara T: Brief exposure to low-pH stress causes irreversible damage to the growing root in Arabidopsis thaliana : pectin-Ca interaction may play an important role in proton rhizotoxicity. J Exp Bot. 2001, 52: 361-368.PubMedView ArticleGoogle Scholar
- Kobayashi Y, Ikka T, Kimura K, Yasuda O, Koyama H: Characterisation of lanthanum toxicity for root growth of Arabidopsis thaliana from the aspect of natural genetic variation. Funct Plant Biol. 2007, 34: 984-994.View ArticleGoogle Scholar
- Ikka T, Kobayashi Y, Iuchi S, Sakurai N, Shibata D, Kobayashi M, Koyama H: Natural variation of Arabidopsis thaliana reveals that aluminum resistance and proton resistance are controlled by different genetic factors. Theor Appl Genet. 2007, 115: 709-719.PubMedView ArticleGoogle Scholar
- Suzuki Y, Kawazu T, Koyama H: RNA isolation from siliques, dry seeds, and other tissues of Arabidopsis thaliana. BioTechniques. 2004, 37: 542-544.PubMedGoogle Scholar
- Chen K, Du L, Chen Z: Sensitization of defense responses and activation of programmed cell death by a pathogen-induced receptor-like protein kinase in Arabidopsis. Plant Mol Biol. 2003, 53: 61-74.PubMedView ArticleGoogle Scholar
- Kinraide TB: Interactions among Ca2+, Na+ and K+ in salinity toxicity: quantitative resolution of multiple toxic and ameliorative effects. J Exp Bot. 1999, 50: 1495-1505.View ArticleGoogle Scholar
- Heidenreich B, Mayer K, Sandermann H, Ernst D: Mercury-induced genes in Arabidopsis thaliana : identification of induced genes upon long-term mercuric ion exposure. Plant, Cell Environ. 2001, 24: 1227-1234.View ArticleGoogle Scholar
- Kawaura K, Mochida K, Yamazaki Y, Ogihara Y: Transcriptome analysis of salinity stress responses in common wheat using a 22 k oligo-DNA microarray. Funct Integr Genomics. 2006, 6: 132-142.PubMedView ArticleGoogle Scholar
- Elbein AD, Pan YT, Pastuszak I, Carroll D: New insights on trehalose: a multifunctional molecule. Glycobiology. 2003, 13: 17R-27R.PubMedView ArticleGoogle Scholar
- Roxas V, Smith R, Allen E, Allen R: Overexpression of glutathione S-transferase/glutathione peroxidase enhances the growth of transgenic tobacco seedlings during stress. Nat Biotechnol. 1997, 15: 988-991.PubMedView ArticleGoogle Scholar
- Garg AK, Kim JK, Owens TG, Ranwala AP, Choi YD, Kochian LV, Wu RJ: Trehalose accumulation in rice plants confers high tolerance levels to different abiotic stresses. Proc Natl Acad Sci USA. 2002, 99: 15898-15903.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang W, Vinocur B, Altman A: Plant responses to drought, salinity and extreme temperatures: towards genetic engineering for stress tolerance. Planta. 2003, 218: 1-14.PubMedView ArticleGoogle Scholar
- Vanderbeld B, Snedden WA: Developmental and stimulus-induced expression patterns of Arabidopsis calmodulin-like genes CML37, CML38 and CML39. Plant Mol Biol. 2007, 64: 683-697.PubMedView ArticleGoogle Scholar
- Ciftci-Yilmaz S, Morsy MR, Song L, Coutu A, Krizek BA, Lewis MW, Warren D, Cushman J, Connolly EL, Mittler R: The EAR-motif of the Cys2/His2-type zinc finger protein Zat7 plays a key role in the defense response of Arabidopsis to salinity stress. J Biol Chem. 2007, 282: 9260-9268.PubMedView ArticleGoogle Scholar
- Kobayashi Y, Hoekenga OA, Itoh H, Nakashima M, Saito S, Shaff JE, Maron LG, Pineros MA, Kochian LV, Koyama H: Characterization of AtALMT1 expression in aluminum-inducible malate release and its role for rhizotoxic stress tolerance in Arabidopsis. Plant Physiol. 2007, 145: 843-852.PubMedPubMed CentralView ArticleGoogle Scholar
- de la Fuente JM, Ramirez-Rodriguez V, Cabrera-Ponce JL, Herrera-Estrella L: Aluminum tolerance in transgenic plants by alteration of citrate synthesis. Science. 1997, 276: 1566-1568.PubMedView ArticleGoogle Scholar
- Schenk PM, Kazan K, Rusu AG, Manners JM, Maclean DJ: The SEN1 gene of Arabidopsis is regulated by signals that link plant defence responses and senescence. Plant Physiol Biochem. 2005, 43: 997-1005.PubMedView ArticleGoogle Scholar
- Devoto A, Ellis C, Magusin A, Chang HS, Chilcott C, Zhu T, Turner JG: Expression profiling reveals COI1 to be a key regulator of genes involved in wound- and methyl jasmonate-induced secondary metabolism, defence, and hormone interactions. Plant Mol Biol. 2005, 58: 497-513.PubMedView ArticleGoogle Scholar
- Last RL, Bissinger PH, Mahoney DJ, Radwanski ER, Fink GR: Tryptophan mutants in Arabidopsis : the consequences of duplicated tryptophan synthase beta genes. Plant Cell. 1991, 3: 345-358.PubMedPubMed CentralGoogle Scholar
- Radwanski ER, Zhao J, Last RL: Arabidopsis thaliana tryptophan synthase alpha: gene cloning, expression, and subunit interaction. Mol Gen Genet. 1995, 248: 657-667.PubMedView ArticleGoogle Scholar
- Vogel G, Fiehn O, Jean-Richard-dit-Bressel L, Boller T, Wiemken A, Aeschbacher RA, Wingler A: Trehalose metabolism in Arabidopsis : occurrence of trehalose and molecular cloning and characterization of trehalose-6-phosphate synthase homologues. J Exp Bot. 2001, 52: 1817-1826.PubMedView ArticleGoogle Scholar
- Gelhaye E, Rouhier N, Navrot N, Jacquot JP: The plant thioredoxin system. Cell Mol Life Sci. 2005, 62: 24-35.PubMedView ArticleGoogle Scholar
- Gupta SD, Wu HC, Rick PD: A Salmonella typhimurium genetic locus which confers copper tolerance on copper-sensitive mutants of Escherichia coli. J Bacteriol. 1997, 179: 4977PubMedPubMed CentralGoogle Scholar
- Holmgren A: Thioredoxin structure and mechanism: conformational changes on oxidation of the active-site sulfhydryls to a disulfide. Structure. 1995, 3: 239-243.PubMedView ArticleGoogle Scholar
- Liu Q, Kasuga M, Sakuma Y, Abe H, Miura S, Yamaguchi-Shinozaki K, Shinozaki K: Two transcription factors, DREB1 and DREB2, with an EREBP/AP2 DNA binding domain separate two cellular signal transduction pathways in drought- and low-temperature-responsive gene expression, respectively, in Arabidopsis. Plant Cell. 1998, 10: 1391-1406.PubMedPubMed CentralView ArticleGoogle Scholar
- Barcelo J, Poschenrieder C, Andreu I, Gunse B: Cadmium-induced decrease of water stress resistance in bush bean plants (Phaseolus vulgaris L. cv. Contender). I. Effects of Cd on water potential, relative water content and cell wall elasticity. J Plant Physiol. 125: 17-25.Google Scholar
- Zimmermann P, Hirsch-Hoffmann M, Hennig L, Gruissem W: GENEVESTIGATOR. Arabidopsis microarray database and analysis toolbox. Plant Physiol. 2004, 136: 2621-2632.PubMedPubMed CentralView ArticleGoogle Scholar
- Matsumoto H: Cell biology of aluminum toxicity and tolerance in higher plants. Int Rev Cytol. 2000, 200: 1-46.PubMedView ArticleGoogle Scholar
- Sudo E, Itouga M, Yoshida-Hatanaka K, Ono Y, Sakakibara H: Gene expression and sensitivity in response to copper stress in rice leaves. J Exp Bot. 2008, 59: 3465-3474.PubMedPubMed CentralView ArticleGoogle Scholar
- Yang ZM, Wang J, Wang SH, Xu LL: Salicylic acid-induced aluminum tolerance by modulation of citrate efflux from roots of Cassia tora L. Planta. 2003, 217: 168-174.PubMedGoogle Scholar
- Borsani O, Valpuesta V, Botella MA: Evidence for a role of salicylic acid in the oxidative damage generated by NaCl and osmotic stress in Arabidopsis seedlings. Plant Physiol. 2001, 126: 1024-1030.PubMedPubMed CentralView ArticleGoogle Scholar
- Metwally A, Finkemeier I, Georgi M, Dietz KJ: Salicylic acid alleviates the cadmium toxicity in barley seedlings. Plant Physiol. 2003, 132: 272-281.PubMedPubMed CentralView ArticleGoogle Scholar
- El-Tayeb MA, El-Enany AE, Ahmed NL: Salicylic acid-induced adaptive response to copper stress in sunflower (Helianthus annuus L.). Plant Growth Regulation. 2006, 50: 191-199.View ArticleGoogle Scholar
- Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95: 14863-14868.PubMedPubMed CentralView ArticleGoogle Scholar
- Obayashi T, Kinoshita K, Nakai K, Shibaoka M, Hayashi S, Saeki M, Shibata D, Saito K, Ohta H: ATTED-II: a database of co-expressed genes and cis elements for identifying co-regulated gene groups in Arabidopsis. Nucleic Acids Res. 2007, D863-9. 35 Database.Google Scholar
- Kobayashi Y, Furuta Y, Ohno T, Hara T, Koyama H: Quantitative trait loci controlling aluminium tolerance in two accessions of Arabidopsis thaliana (Landsberg erecta and Cape Verde Islands). Plant Cell Environ. 2005, 28: 1516-1524.View ArticleGoogle Scholar
- Berardini TZ, Mundodi S, Reiser L, Huala E, Garcia-Hernandez M, Zhang P, Mueller LA, Yoon J, Doyle A, Lander G, Moseyko N, Yoo D, Xu I, Zoeckler B, Montoya M, Miller N, Weems D, Rhee SY: Functional annotation of the Arabidopsis genome using controlled vocabularies. Plant Physiol. 2004, 135: 745-755.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.