- Research article
- Open Access
Identification of genes and pathways associated with aluminum stress and tolerance using transcriptome profiling of wheat near-isogenic lines
BMC Genomics volume 9, Article number: 400 (2008)
Aluminum is considered the most limiting factor for plant productivity in acidic soils, which cover large areas of the world's potential arable lands. The inhibition of root growth is recognized as the primary effect of Al toxicity. To identify genes associated with Al stress and tolerance, transcriptome analyses of four different wheat lines (2 Al-tolerant and 2 Al sensitive) that differ in their response to Al were performed.
Microarray expression profiling revealed that 83 candidate genes are associated with Al stress and 25 are associated with tolerance. The stress-associated genes include important enzymes such as pyruvate dehydrogenase, alternative oxidase, and galactonolactone oxidase, ABC transporter and ascorbate oxido-reducatase. The Al tolerance-associated genes include the ALMT-1 malate transporter, glutathione S-transferase, germin/oxalate oxidase, fructose 1,6-bisphosphatase, cysteine-rich proteins, cytochrome P450 monooxygenase, cellulose synthase, zinc finger transcription factor, disease resistance response protein and F-box containing domain protein.
In this survey, we identified stress- and tolerance-associated genes that may be involved in the detoxification of Al and reactive oxygen species. Alternative pathways could help maintain the supply of important metabolites (H2O2, ascorbate, NADH, and phosphate) needed for Al tolerance and root growth. The Al tolerance-associated genes may be key factors that regulate these pathways.
Aluminum is considered as the most limiting factor for plant productivity in acidic soils. It is estimated that over 50% of the world's potential arable land surface is composed of acid soils mostly distributed in developing countries [1, 2]. Al tolerance is second to drought tolerance for its importance as agronomic trait for worldwide crop production. The root apex is considered the first target of Al toxicity and the reduction in root biomass leads to poor uptake of water and nutrients . At the cellular level, Al toxicity results from a broad spectrum of deleterious effects caused by the Al3+ ion, which is considered the most toxic species of Al for both plant and animal cells under low pH conditions . It was proposed that the toxic effect of Al can be reduced by chelating Al in the rhizosphere with organic anions. Exudation of a variety of organic anions such as malate, citrate, or oxalate upon exposure to Al have been reported [2, 5, 6] and Al tolerance of several species can be enhanced by increasing organic acid biosynthesis [7–10]. Overexpression of the Al-inducible malate transporter improves Al tolerance in barley  while overexpression of the Sb MATE protein, a putative citrate transporter, improves Al tolerance in Arabidopsis and wheat . The release of phosphate was also associated with Al tolerance, indicating that multiple mechanisms are involved in wheat . In Arabidopsis, the major Al tolerance quantitative trait loci (QTL) was proposed to be proximal to the AtALMT1 malate transporter and to involve additional genes . Fine mapping indicated that Al tolerance is more closely associated with a chromosome location that is distal to AtALMT by 1400 to 2100 kbp indicating that other gene products may regulate or complement its activity . This regulation may be controlled at the transcriptional, translational, post-translational and enzymatic levels. The association of organic acid release and Al tolerance is not universal as several Al-sensitive plants were found to secrete large amounts of organic acids [15, 16]. These results suggest that organic acid excretion, per se, is not a tolerance (or resistance) mechanism, but a consequence of biochemical reactions required for Al tolerance. In Arabidopsis, there was no induction of AtALMT1 by Al3+ in the stop1 Arabidopsis mutant suggesting that the STOP1 zinc finger protein is required to activate ALMT1 transcription . The ALS3 gene encodes an ABC transporter-like protein that is required for Al tolerance/resistance and may function to redistribute accumulated Al away from sensitive tissues to protect growing roots from Al toxicity . Several genes up-regulated during Al exposure have been identified using various molecular approaches [18–22]. Overexpression of some of these genes in transgenic plants has resulted in modest improvement of Al tolerance, suggesting that they alleviate part of the toxicity caused by Al [23–27]. In rice, Al tolerance appears to be a complex multigenic trait that involves all twelve chromosomes. However, fewer loci were reported to be involved in other grasses [2, 28, 29]. Two major loci on chromosome 4DL and 6A, and additional loci with additive effects are involved in Al tolerance of the wheat cultivar Atlas66 [30, 31]. The availability of near-isogenic lines derived from the cultivar Atlas66 having similar QTLs for Al tolerance provides a useful tool to identify genes associated with Al stress and tolerance [30, 31].
In this study, several genes associated with Al stress and tolerance were identified using transcriptome analyses. The putative functions of identified genes in several biochemical pathways are discussed in relation to stress responses and the maintenance of root growth.
To identify genes associated with Al stress and tolerance in wheat, a large scale expression profiling study was initiated using the Affymetrix GeneChip® Wheat Genome Array which allows the screening of 55052 transcripts. To better discriminate between genes that are associated with Al stress and tolerance, four different wheat lines (two Al-tolerant and two Al-sensitive) were analyzed. In addition to the well studied tolerant wheat Atlas66 and the sensitive wheat Bounty used in our previous study , we used two near isogenic lines (NILs) derived from a cross between Atlas66 and the sensitive cultivar Century . The sensitive NIL OK91G108 (named Century-S thereafter) has a high degree of genetic similarity (96.9%) with the tolerant NIL OK91G106 (named Century-T thereafter). The major difference between tolerant and sensitive plants is the ability to maintain growth under high Al concentrations. The tolerant cultivars used in this study are able to grow for several days (50% rate of control plants or a root growth inhibition (RGI) = 50%) in the presence of 50 μM Al (result not shown). The sensitive cultivars are unable to grow at this high Al concentration and an RGI of 50% is obtained in the presence of 5 μM Al. We have previously shown that the stress-associated genes are expressed at comparable levels in the sensitive and tolerant cultivars after 24 hours when they are exposed to Al concentrations resulting in a similar RGI . The microarray experiments were thus designed to identify genes that are differentially expressed after 24 hours of Al exposure at concentrations resulting in 50% RGI for all lines used.
The number of Al-regulated genes (differential expression between Al-treated and untreated plants) is over 1000 when cultivars are analyzed separately. However, only 263 genes are differentially expressed in all four cultivars. ANOVA analyses with stringent p values retained 83 of these genes as highly significantly regulated (Table 1; additional file 1). Among these genes, 20 were previously identified as up-regulated by Al (references in Table 1). We randomly selected 4 of the stress candidate genes for qRT-PCR validation (see bold probeset IDs in Table 1). Two of them (genes # 2 and 3) are strongly over-expressed in response to Al, one (gene # 5, oxalate oxidase) was previously shown to be up-regulated by Al and the fourth one (gene # 49) was represented by 5 different genes (genes # 49, 61, 73, 74, and 76).
As found for stress-associated genes, microarray analyses allowed the identification of more than 1000 genes potentially associated with Al tolerance when selection was performed using only one pair of wheat lines (differential expression between Atlas66/Bounty or Century-T/Century-S exposed to Al) but this number was reduced to 69 common genes when both pairs were considered. This result suggests that a large number of genes are associated with cultivar specific responses and only a limited number of them are associated with Al tolerance. This list was reduced to 25 genes using ANOVA analyses with stringent p values (Table 2; additional file 2). These genes were subdivided into constitutively or differentially expressed genes based on the difference in expression between Al-treated and non-treated control plants. Table 2A lists the 6 constitutively expressed genes classified in decreasing order of expression (average fold expression of the two tolerant compared to the two sensitive cultivars). The most strongly expressed gene is the previously identified Al malate transporter ALMT-1 . We selected two of these genes for qRT-PCR validation (gene # 84 encodes ALMT-1 and is more expressed in tolerant cultivars while gene # 89 encodes a protein containing an F-box domain and is more expressed in sensitive cultivars; Table 2A). Table 2B lists the 19 genes that are regulated by Al and differentially expressed between the two tolerant and the two sensitive cultivars. Four of these genes were selected for validation by qRT-PCR (see bold probeset IDs in Table 2B).
Quantitative RT-PCR validation of microarray data
We selected ten genes for validation of the microarray data by qRT-PCR (4 among the stress-associated candidates and 6 among the tolerance-associated candidates, as indicated above). Four different Al concentrations (0, 5, 50 and 250 μM Al) were used to determine the effect on gene expression. Four new replicates of each Al exposure conditions using all wheat lines were prepared. The values indicated in each histogram column represent the expression levels relative to the first column (cultivar Atlas66 not exposed to Al) (Fig. 1A–J). The level of expression of the tested genes under control conditions is generally similar in the tolerant and sensitive cultivars.
Genes were classified as associated with Al tolerance when there was a significant difference in expression in both tolerant cultivars (Atlas66 and Century-T exposed to 50 μM Al) compared to their respective sensitive counterparts (Bounty and Century-S exposed to 5 μM Al). These conditions allow the detection of either 1) Al-regulated genes that may be expressed differentially in tolerant and sensitive lines, or 2) genes not regulated by Al (constitutively expressed) but showing constitutive difference in expression between tolerant and sensitive lines. Once classified as associated with stress or tolerance, gene regulation by Al was confirmed by comparing the Al-treated and non-treated controls. qRT-PCR, analyses confirmed that three (Ta.192.1.S1_at, Ta.12921.1.S1_x_at and Ta.392.2.S1_at, Fig. 1A–C) of the four stress candidate genes are stress-associated while the 4th gene (Ta.5557.1.S1_x_at, Fig. 1D) is more closely associated with Al tolerance. Among the six selected tolerance candidates genes, four were confirmed as associated with Al tolerance (Fig. 1E–H: Ta.30659.1.S1_at, TaAffx.16664.1.A1_at, Ta23271.1.S1_s_at, Ta.8545.1.S1_x_at) while the last two genes (Ta.21314.1.S1_x_at and Ta.24632.1.S1_at, Fig. 1I–J), were not clearly linked to either stress or tolerance.
The correlation between gene expression data obtained by microarray and qRT-PCR is presented in Fig. 2. Al-regulated genes are shown in Fig. 2A for the cultivars Atlas66 and Bounty. The fold change in gene expression (Al-treated VS control) obtained using microarray experiments compared to the fold change obtained using qRT-PCR gives a correlation coefficient of 0.61. A similar analysis was performed to compare the change in gene expression between the tolerant and sensitive cultivars (Atlas66 50 μM Al/Bounty 5 μM Al and Century-T 50 μM Al/Century-S 5 μM Al) calculated from the microarray or qRT-PCR data (Fig. 2B). In this case, a correlation coefficient of 0.77 was obtained between the two experimental methods. These correlation coefficients are very good considering that microarray data are semi-quantitative and subject to error for multigene families where different transcripts could hybridize to similar probes on the array. qRT-PCR data are more specific since the amplification of single transcripts are confirmed by melting curves and gel analyses. The Ta.5557.1.S1_x_at was initially selected as a stress-associated gene by microarray profiling (gene no 5 in Table 1) while it was reclassified as associated with Al tolerance using qRT-PCR (Fig. 1D). This may be related to the presence of multiple oxalate oxidase genes in wheat. Based on these results, we conclude that most of the genes identified in Table 1 and Table 2 are associated with Al stress and tolerance. However, real-time PCR analyses will be required to validate the other genes that were identified using microarrays.
Several Al-associated QTLs with major and minor contributions were reported, indicating that Al tolerance is a multigenic trait. A large scale transcriptome analysis was initiated to identify some gene components of this multigenic trait. The genes that are differentially expressed between tolerant and sensitive cultivars are of particular interest since the tolerant cultivars are exposed to Al concentrations that are ten times higher, yet they are still growing at the same rate as the sensitive cultivars. This design allows the distinction of genes that are regulated in response to stress from those that are associated with Al tolerance. Transcriptome analyses allowed the identification of several genes with known function that could improve our understanding of the biochemical processes and pathways involved in Al stress and tolerance. As indicated in the result section, 20 of the 83 genes associated with the Al stress response were identified in previous studies indicating that the genes associated with Al stress are appropriately identified using this approach.
We have previously shown that Al inhibits a redox reaction associated with root growth and Al tolerance . The perturbation of this redox reaction by Al could lead to the accumulation of different reactive oxygen species which can stimulate a redox signalling pathway and increase the expression of antioxidant enzymes . The accumulation of reactive oxygen species during Al exposure was observed in maize, soybean and wheat [36–38]. Al was also shown to affect mitochondrial functions leading to ROS production . Table 1 shows that the genes encoding WCI-5 (gene # 2), WAS-2 (genes # 10 and 23) and the EF-hand Ca2+-binding protein (genes # 40, 41, 43, 54 and 55) are up-regulated during Al exposure. In Arabidopsis, the WCI-5 and WAS-2 genes were shown to be induced via a MAP kinase pathway associated with the pathogen response . In wheat seedlings, the WCI-5 gene was shown to be induced by pathogen infection  and mediated by the EF-hand Ca2+-binding protein CCD1 . These results suggest that a MAP kinase pathway is stimulated by Al in wheat roots as observed in cell suspension cultures of Coffea arabica L. . Three different protein kinase-like genes (genes # 37, 60, 63) are up-regulated during Al exposure suggesting that different signal transduction pathways may be activated. We identified two other potential signalling genes associated with Al tolerance. The first one encodes a putative C2H2 zinc finger transcription factor (gene # 107) that has a 33 aa region sharing 50% identity with the STOP1 C2H2 zinc finger protein required for Al3+ induction of ALMT1 in Arabidopsis . The other gene encodes an F-box containing protein (gene # 89), which is member of a family known for its involvement in the controlled degradation of target proteins. An Arabidopsis F-box protein homologous to a bacterial redox sensor was proposed to be involved in the cell wall receptor-like associated kinase (WAK1) signalling pathway . Since WAK1 was shown to be transiently overexpressed during Al exposure and associated with Al tolerance , the reduced expression of an F-Box protein may stimulate this pathway by allowing a stronger expression of WAK1 in tolerant cultivars. The TIR1 gene is another F-Box protein in Arabidopsis that interacts with different proteins in the SCFTIR1 complex to mediate responses to auxin [45, 46]. These two genes (C2H2 zinc finger and F-Box) require further characterization and functional analysis to understand their roles in the response to Al.
Management of Al-associated stresses and maintenance of energy supply
The identification of stress-regulated genes provide new tools to reduce Al stress, as shown by the ectopic overexpression of some Al stress-regulated genes [23–26]. Among the 83 candidate genes regulated by Al stress (Table 1), several could play a role in alleviating phosphate deficiency and provide energy to fight oxidative stress. Nutrient deficiency, and especially phosphate, occurs in the presence of Al due to the precipitation of Al-phosphate . Under phosphate deficiency, several pathways are hindered due to the reduced availability of ATP and related nucleoside phosphates. These imbalances lead to the induction of alternative pathways of glycolysis to maintain energy and carbon skeletons for key metabolic processes . The fructose-1,6-bisphosphatase (gene # 101) associated with Al tolerance (Table 2) is normally inhibited by phosphate but is stimulated under low phosphate availability (Fig. 3, Box 3). This activity produces fructose-6-P and liberates Pi. The glycolytic alternative enzyme pyrophosphate phosphofructokinase is stimulated by low phosphate and can resynthesise fructose-1,6-bisphosphate and liberate more Pi. This apparently futile cycle is useful to conserve ATP and recycle Pi that is readily available in the pyrophosphate molecule . Furthermore, since sucrose is the most abundant photosynthate available in the phloem sap, fructose (from sucrose) is readily available and will enter the glycolytic pathway at the level of fructose-6-P thereby maintaining a stable supply of this metabolite. The high availability of fructose-6-P (and other hexose-P) can provide sufficient metabolites to support both the pentose phosphate pathway and organic acid synthesis (Fig. 3, Box 3). Al-tolerant rye was found to maintain a higher content of glucose-6-phosphate than Al-sensitive wheat in response to Al exposure . This metabolite is an important precursor of the pentose phosphate pathway (with glucose-6-phosphate dehydrogenase (G6PDH) as the first enzyme) and is an important source of NADPH production. The importance of energy metabolism in modulating Al tolerance is exemplified in Pseudomonas fluorescens where G6PDH was shown to play a pivotal role in enhancing NADPH production and protecting against ROS toxicity . The up-regulation of pyruvate dehydrogenase (gene # 65) and alternative oxidase (gene # 72) supports the involvement of an alternative pathway to maintain energy production (Fig. 3, Box 3) . Members of the alternative oxidase gene family were found to be highly responsive to oxidative stress and to reduce mitochondrial ROS production while maintaining NADH supply [51, 52]. In Al tolerant plants, a higher NADH supply (and NADPH through the pentose phosphate pathway) would provide reducing power to regenerate ascorbate and glutathione to fight oxidative stress (Fig. 3, Box 3).
Exudation of organic acids has received much attention as an Al tolerance (resistance) mechanism. The secretion of malate or other organic acids (citrate, oxalate) can chelate Al in the rhizosphere and has been associated with Al tolerance in several species [7–10, 33]. Enzymes involved in malate synthesis (PEP carboxylase and malate dehydrogenase, Fig. 3 Box 3) were shown to be stimulated during Al stress [53, 54]. These enzymes are also stimulated during Pi deficiency where malate is known to mobilize external Pi . These two actions are complementary since the chelation of Al will also free Pi from precipitated Al-phosphate in the rhizosphere. Some of the genes identified in this study may complement the ALMT1 malate transporter to improve Al exclusion. An interesting possibility is that organic acids, such as malate or oxalate, could form intracellular Al-bound complexes that are transported outside the cell (Fig. 3, Box 3). It is possible that the ALMT1 malate transporter (gene # 84) associated with Al tolerance can transport free malate but that it would require a co-factor to transport the larger Al-bound malate complex. Oxalate is another organic acid known to be a strong Al chelator associated with Al tolerance in Taro . This organic acid is another potential molecule that could form an Al complex transported out of the cell by an anion transporter such as ALMT1 or the ABC transporter (gene # 22) that is up-regulated during Al exposure.
Several proteins encoded by genes that are up-regulated during Al exposure have a high content of cysteine (Fig. 3, Box 2). The high amount of thiol groups in these proteins may help protect enzymes containing thiol sensitive sites by providing new targets of oxidation. The presence of 11 cysteine residues in the small protein (128 aa) encoded by gene # 90 may explain why this gene is associated with Al tolerance (Fig. 3, Box 2). The up-regulation of four different GSTs (genes # 1, 70, 91, 96), containing cysteine residues in their active sites, indicate that these enzyme may be inactivated by oxidative stress. The synthesis of new enzyme molecules during Al exposure may partially compensate for GST inactivation. On the other hand, two of the four GSTs are associated with Al tolerance (genes # 91, 96). These GSTs may help to reduce oxidative stress as suggested for parB, a tobacco GST isoform .
Maintenance of root growth
The ability to maintain root elongation under Al stress may require complementary mechanisms such as Al sequestering (outside the cell, or Al exclusion) and increased expression of genes involved in detoxification or root growth. We identified several genes that could play a role in maintaining root growth. Different enzymes such as sugar hydrolases and transferases may participate in cell wall remodelling/loosening in the elongation zone (genes # 32, 36, 46, 77). During the cell wall loosening process, ascorbate provides an electron to reduce Cu2+ to Cu+ (Fig. 3, Box 1) [56, 57], which participates in the Fenton reaction with H2O2 to ultimately produce •OH radicals. The non-enzymatic cleavage of cell wall polysaccharides thus requires the presence of three basic molecules: ascorbate, Cu+2 and H2O2 (grey boxes in Fig. 3) [56, 57]. The enzymatic reactions required to maintain/regenerate ascorbate and provide H2O2 are essential components to maintain root growth and are likely associated with Al tolerance. The redox activity previously shown to be associated with Al tolerance  may represent the transfer of electrons required to regenerate the apoplasmic ascorbate. The great avidity of Al3+ for electrons (Fig. 3, Box 1) could prevent the regeneration of apoplasmic ascorbate and inhibit root growth. Oxalate oxidase may help maintain the production of H2O2 for the Fenton reaction (Fig. 3). However, H2O2 could become a substrate for cell wall peroxidases which are not inhibited when ascorbate levels are reduced . The inhibition of peroxidase activity was associated with Al tolerance [59, 60] and may help maintain H2O2 levels needed for the non-enzymatic wall loosening process. In our study, three peroxidases (genes # 13, 21, 58) are up-regulated and could be related to the reduced growth rate under our experimental conditions. Transcripts encoding blue copper binding proteins (BCBP) (gene # 25) is possibly up-regulated to stimulate the Cu2+-dependent Fenton reaction needed for cell wall loosening. When apoplasmic ascorbate loses an electron, it is transformed into monodehydroascorbate/dehydroascorbate (MDHA/DHA) and new electrons are needed to regenerate ascorbate and avoid degradation into oxalate (Fig. 3, Box 1). The up-regulation of an ascorbate dependent oxidoreductase (gene # 24) could be associated with the regeneration of apoplasmic ascorbate (Fig. 3, Box 3). Impaired ascorbate metabolism was proposed to be involved in the reduction of root growth in squash roots exposed to Al  while a higher level of ascorbate and gluthathione was shown to be associated with Al tolerance in tobacco . Furthermore, feeding with D-glucose or L-galactono-γ-lactone to enhance ascorbate levels was able to improve Al tolerance in rice . These results indicate that maintaining a high ascorbate level is an essential aspect of Al tolerance. Gene # 30 annotated as an arabinonolactone oxidase is homologous to galactonolactone oxidase involved in ascorbate synthesis (Fig. 3, Box 3). This last gene participates in the regeneration of cytoplasmic ascorbate while providing additional NADH (from glucose, Fig. 3, Box 3). The arabinonolactone oxidase is known to synthesize erythroascorbate in yeast. Interestingly, the overexpression of the yeast arabinonolactone oxidase enzyme in rice improves Al tolerance .
Genome wide expression profiling and qRT-PCR using different wheat cultivars subjected to controlled stress treatments allowed the identification of several new genes associated with Al stress and tolerance. Several genes associated with Al tolerance could play an important role in maintaining the energy balance and Al exclusion. The maintenance of ascorbate homeostasis is proposed to be a key element to sustain elongation growth. The availability of genes associated with Al tolerance provides new tools for QTL analyses and for breeding programs aimed at improving Al tolerance of cultivated crops.
Plant material, growth and Al exposure conditions
Two wheat cultivars with a high tolerance to Al (Triticum aestivum L. cv. Atlas66 and the near isogenic line OK91G106, named Century-T in this work; ) and two wheat cultivars with low tolerance to Al (T. aestivum L. cv. Bounty and the near isogenic line OK91G108, named Century-S in this work; ) were grown as previously described, and treated under conditions where Al remains mostly in the Al3+ form . To reduce pH variations and ensure that Al speciation was stable throughout the experiment, at least 100 ml of solution was used for each plant. The root growth inhibition (RGI) is expressed as 100 × [1-(root growth of Al-treated seedling divided by the root growth of control seedlings)]. Replicate experiments were performed on different days with one series of Al concentration per day per cultivar in order to rapidly collect root tips.
RNA isolation and microarray profiling
Root tips (5–10 mm) were isolated after 24 hours of exposure to Al and frozen on dry ice. Total RNA was isolated from the root tips (5 mm long) of 50 plants collected from the various genotypes using the RNeasy Plant Mini Kit (Qiagen). The RNA quality was assessed on agarose gels and with the Bioanalyzer 2100 (Agilent). Microarray profiling was performed according to Affymetrix protocols at the Functional Genomics Platform of McGill University and Génome Québec Innovation Centre using the Affymetrix GeneChip® Wheat Genome Array. The microarray results were deposited in the database of ArrayExpress under the accession number E-TABM-454. Three biological replicates of the different lines were treated with Al concentrations resulting in 50% RGI: 50 μM for the tolerant cultivars Atlas66 and Century-T and 5 μM for the sensitive cultivars Bounty and Century-S. The wheat NILs Century-T and Century-S were used only to prepare the Al-treated plants due to the limited number of seeds available. In this case, the average value obtained for the untreated (not exposed to Al) Atlas66 and Bounty was used as the reference.
Quantitative Real-Time PCR
The level of expression of selected genes was validated by Real-Time PCR. Probes were designed using the primer3 software and primer specificity was verified on agarose gels and by Tm measurements at the end of the Real-Time PCR reactions. Total RNA was prepared (RNeasy Plant Mini Kit; Qiagen) from the root tips (5 mm long) of 50 plants collected from the various genotypes and reverse transcribed using SuperScript II reverse transcriptase and random hexamers (SuperSript™ First-Strand; Invitrogen). Quantitative Real-Time was performed on an ABI 7000 Real Time Cycler using the specific primers described in additional file 3. One-tenth dilutions of the cDNAs were used as template for the qRT-PCR in a total volume of 25 μL as follows: 12.5 μL SYBR Green (Platinum® SYBR® Green qPCR SuperMix-UDG, Invitrogen) 1.0 μL primer mix (50:50 mix of forward and reverse primers at 10 pmol/μL each) and 2 μL template. The reaction conditions were: 10 min at 95°C followed by 40 cycles of 1 min at 95°C and 30 s at 60°C.
The microarray data were analyzed using the robust multi-array average (RMA) software (RMA version 0.2) of background adjusted, normalized, and log transformed perfect match values . A two-fold cut-off value (log2 ≥ 1) was arbitrarily set to indicate differential gene expression between two samples. Genes showing a differential expression greater than two-fold (RMA differential expression of log2 ≥ 1) between the Atlas66 and Bounty controls were excluded from the reference set as these may represent inherited cultivar basal expression levels that could bias the analysis. Overall, this excluded 1.26% of the genes on the microarray (773 of 55,052) demonstrating that most genes are expressed at a similar level between the two cultivars. Genes that are differentially expressed to the same extent in all four cultivars exposed to Al (giving 50% RGI) compared to the non-treated cultivars were selected and classified as candidates for stress-associated genes. The differentially expressed genes between the two tolerant cultivars exposed to Al (Atlas66 and Century-T) and their respective sensitive counterparts (Bounty and Century-S) exposed to Al concentrations resulting in 50% RGI were classified as candidate genes associated with Al tolerance. These genes were subdivided in two groups (constitutively expressed or Al-regulated) based on the average level of regulation in the two tolerant cultivars. Genes were considered constitutively expressed when the differential signal between Al-treated and non-treated samples was less than two-fold (log2 < 1). An analysis of variance was performed using GraphPad InStat 3 to select genes that are differentially expressed under the conditions specified for each analysis.
For Quantitative Real-Time PCR analyses, the amplification efficiency (90% to 100%) for the different primer sets was determined by amplification of cDNA dilution series using 80, 20, 10, 5, 2.5, and 1.25 ng per reaction (data not shown). The variance (standard error) was very small between the PCR replicates for a same biological sample (Real-Time experimental replicate) compared to the variance between the different biological replicates. Calculations and statistical analyses were performed by ANOVA on the mean ΔCT (CT of each gene – CT of 18S RNA used as load control) of different biological replicates, as described in the Results section.
von Uexküll HR, Mutert E: Global extent, development and economic impact of acid soils. Plant Soil. 1995, 171: 1-15. 10.1007/BF00009558.
Kochian LV, Pineros MA, Hoekenga OA: The physiology, genetics and molecular biology of plant aluminum resistance and toxicity. Plant Soil. 2005, 274: 175-195. 10.1007/s11104-004-1158-7.
Kochian LV: Cellular mechanisms of aluminum toxicity and resistance in plants. Annu Rev Plant Physiol Plant Mol Biol. 1995, 46: 237-260. 10.1146/annurev.pp.46.060195.001321.
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.
Delhaize E, Ryan PR, Randall PJ: Aluminum tolerance in wheat (Triticum aestivum L.) (II. Aluminum-stimulated excretion of malic acid from root apices). Plant Physiol. 1993, 103: 695-702.
Kidd PS, Llugany M, Poschenrieder C, Gunse B, Barcelo J: The role of root exudates in aluminium resistance and silicon-induced amelioration of aluminium toxicity in three varieties of maize (Zea mays L.). J Exp Bot. 2001, 52: 1339-1352. 10.1093/jexbot/52.359.1339.
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. 10.1126/science.276.5318.1566.
Koyama H, Kawamura A, Kihara T, Hara T, Takita E, Shibata D: Overexpression of mitochondrial citrate synthase in Arabidopsis thaliana improved growth on a phosphorus-limited soil. Plant Cell Physiol. 2000, 41: 1030-1037. 10.1093/pcp/pcd029.
Tesfaye M, Temple SJ, Allan DL, Vance CP, Samac DA: Overexpression of malate dehydrogenase in transgenic alfalfa enhances organic acid synthesis and confers tolerance to aluminum. Plant Physiol. 2001, 127: 1836-1844. 10.1104/pp.127.4.1836.
Anoop VM, Basu U, McCammon MT, McAlister-Henn L, Taylor GJ: Modulation of citrate metabolism alters aluminum tolerance in yeast and transgenic canola overexpressing a mitochondrial citrate synthase. Plant Physiol. 2003, 132: 2205-2217. 10.1104/pp.103.023903.
Delhaize E, Ryan PR, Hebb DM, Yamamoto Y, Sasaki T, Matsumoto H: Engineering high-level aluminum tolerance in barley with the ALMT1 gene. Proc Natl Acad Sci U S A. 2004, 101: 15249-15254. 10.1073/pnas.0406258101.
Magalhaes JV, Liu J, Guimaraes CT, Lana UG, Alves VM, Wang YH, Schaffert RE, Hoekenga OA, Pineros MA, Shaff JE, Klein PE, Carneiro NP, Coelho CM, Trick HN, Kochian LV: A gene in the multidrug and toxic compound extrusion (MATE) family confers aluminum tolerance in sorghum. Nat Genet. 2007, 39: 1156-1161. 10.1038/ng2074.
Pellet DM, Papernik LA, Kochian LV: Multiple Aluminum-Resistance Mechanisms in Wheat (Roles of Root Apical Phosphate and Malate Exudation). Plant Physiol. 1996, 112: 591-597.
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 U S A. 2006, 103: 9738-9743. 10.1073/pnas.0602868103.
Ishikawa S, Wagatsuma T, Sasaki R, Ofei-Manu P: Comparison of the amount of citric and malic acids in Al media of seven plant species and two cultivars each in five plant species. Soil Sci Plant Nutr. 2000, 46: 751-758.
Pineros MA, Shaff JE, Manslank HS, Alves VM, Kochian LV: Aluminum resistance in maize cannot be solely explained by root organic acid exudation. A comparative physiological study. Plant Physiol. 2005, 137: 231-241. 10.1104/pp.104.047357.
Iuchi S, Koyama H, Iuchi A, Kobayashi Y, Kitabayashi S, Kobayashi Y, Ikka T, Hirayama T, Shinozaki K, Kobayashi M: Zinc finger protein STOP1 is critical for proton tolerance in Arabidopsis and coregulates a key gene in aluminum tolerance. Proc Natl Acad Sci U S A. 2007, 104: 9900-9905. 10.1073/pnas.0700117104.
Snowden KC, Gardner RC: Five genes induced by aluminum in wheat (Triticum aestivum L.) roots. Plant Physiol. 1993, 103: 855-861. 10.1104/pp.103.3.855.
Hamel F, Breton C, Houde M: Isolation and characterization of wheat aluminum-regulated genes: possible involvement of aluminum as a pathogenesis response elicitor. Planta. 1998, 205: 531-538. 10.1007/s004250050352.
Richards KD, Schott EJ, Sharma YK, Davis KR, Gardner RC: Aluminum induces oxidative stress genes in Arabidopsis thaliana. Plant Physiol. 1998, 116: 409-418. 10.1104/pp.116.1.409.
Sasaki T, Ezaki B, Matsumoto H: A gene encoding multidrug resistance (MDR)-like protein is induced by aluminum and inhibitors of calcium flux in wheat. Plant Cell Physiol. 2002, 43: 177-185. 10.1093/pcp/pcf025.
Watt DA: Aluminium-responsive genes in sugarcane: identification and analysis of expression under oxidative stress. J Exp Bot. 2003, 54: 1163-1174. 10.1093/jxb/erg128.
Basu U, Good AG, Taylor GJ: Transgenic Brassica napus plants overexpressing aluminium-induced mitochondrial manganese superoxide dismutase cDNA are resistant to aluminium. Plant Cell Environ. 2001, 24: 1269-1278. 10.1046/j.0016-8025.2001.00783.x.
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. 10.1104/pp.122.3.657.
Ezaki B, Katsuhara M, Kawamura M, Matsumoto H: Different mechanisms of four aluminum (Al)-resistant transgenes for Al toxicity in Arabidopsis. Plant Physiol. 2001, 127: 918-927. 10.1104/pp.127.3.918.
Ezaki B, Sasaki K, Matsumoto H, Nakashima S: Functions of two genes in aluminium (Al) stress resistance: repression of oxidative damage by the AtBCB gene and promotion of efflux of Al ions by the NtGDI1 gene. J Exp Bot. 2005, 56: 2661-2671. 10.1093/jxb/eri259.
Sivaguru M, Ezaki B, He ZH, Tong H, Osawa H, Baluska F, Volkmann D, Matsumoto H: Aluminum-induced gene expression and protein localization of a cell wall-associated receptor kinase in Arabidopsis. Plant Physiol. 2003, 132: 2256-2266. 10.1104/pp.103.022129.
Matos M, Camacho MV, Perez-Flores V, Pernaute B, Pinto-Carnide O, Benito C: A new aluminum tolerance gene located on rye chromosome arm 7RS. Theor Appl Genet. 2005, 111: 360-369. 10.1007/s00122-005-2029-1.
Wight CP, Kibite S, Tinker NA, Molnar SJ: Identification of molecular markers for aluminium tolerance in diploid oat through comparative mapping and QTL analysis. Theor Appl Genet. 2006, 112: 222-231. 10.1007/s00122-005-0114-0.
Berzonsky WA: The genomic inheritance of aluminum tolerance in Atlas-66 wheat. Genome. 1992, 35: 689-693.
Ma HX, Bai GH, Carver BF, Zhou LL: Molecular mapping of a quantitative trait locus for aluminum tolerance in wheat cultivar Atlas-66. Theor Appl Genet. 2005, 112: 51-57. 10.1007/s00122-005-0101-5.
Carver BF, Whitmore WE, Smith EL, Bona L: Registration of four aluminum-tolerant winter wheat germplasms and two susceptible near-isolines. Crop Sci. 1993, 33: 1113-1114.
Sasaki T, Yamamoto Y, Ezaki B, Katsuhara M, Ahn SJ, Ryan PR, Delhaize E, Matsumoto H: A wheat gene encoding an aluminum-activated malate transporter. Plant J. 2004, 37: 645-653. 10.1111/j.1365-313X.2003.01991.x.
Maltais K, Houde M: A new biochemical marker for aluminium tolerance in plants. Physiol Plant. 2002, 115: 81-86. 10.1034/j.1399-3054.2002.1150109.x.
Foyer CH, Noctor G: Redox homeostasis and antioxidant signaling: a metabolic interface between stress perception and physiological responses. Plant Cell. 2005, 17: 1866-1875. 10.1105/tpc.105.033589.
Boscolo PR, Menossi M, Jorge RA: Aluminum-induced oxidative stress in maize. Phytochemistry. 2003, 62: 181-189. 10.1016/S0031-9422(02)00491-0.
Cakmak I, Horst W J: Effect of aluminium on lipid peroxidation, superoxide dismutase, catalase, and peroxidase activities in root tips of soybean (Glycinemax). Physiol Plant. 1991, 83: 463-468. 10.1111/j.1399-3054.1991.tb00121.x.
Lukaszewski KM, Blevins DG: Root growth inhibition in boron-deficient or aluminum-stressed Squash may be a result of impaired ascorbate metabolism. Plant Physiol. 1996, 112: 1135-1140.
Yamamoto Y, Kobayashi Y, Devi SR, Rikiishi S, Matsumoto H: Aluminum toxicity is associated with mitochondrial dysfunction and the production of reactive oxygen species in plant cells. Plant Physiol. 2002, 128: 63-72. 10.1104/pp.128.1.63.
Takezawa D: Elicitor- and A23187-induced expression of WCK-1, a gene encoding mitogen-activated protein kinase in wheat. Plant Mol Biol. 1999, 40: 921-933. 10.1023/A:1006263607135.
Gorlach J, Volrath S, Knauf-Beiter G, Hengy G, Beckhove U, Kogel KH, Oostendorp M, Staub T, Ward E, Kessmann H, Ryals J: Benzothiadiazole, a novel class of inducers of systemic acquired resistance, activates gene expression and disease resistance in wheat. Plant Cell. 1996, 8: 629-643. 10.1105/tpc.8.4.629.
Takezawa D: A rapid induction by elicitors of the mRNA encoding CCD-1, a 14kDa Ca2+ -binding protein in wheat cultured cells. Plant Mol Biol. 2000, 42: 807-817. 10.1023/A:1006431724090.
Arroyo-Serralta GA, Ku-Gonzalez A, Hernandez-Sotomayor SM, Zuniga Aguilar JJ: Exposure to toxic concentrations of aluminum activates a MAPK-like protein in cell suspension cultures of Coffea arabica. Plant Physiol Biochem. 2005, 43: 27-35. 10.1016/j.plaphy.2004.12.003.
Vandenabeele S, Van Der KK, Dat J, Gadjev I, Boonefaes T, Morsa S, Rottiers P, Slooten L, Van Montagu M, Zabeau M, Inze D, Van Breusegem F: A comprehensive analysis of hydrogen peroxide-induced gene expression in tobacco. Proc Natl Acad Sci U S A. 2003, 100: 16113-16118. 10.1073/pnas.2136610100.
Gray WM, del Pozo JC, Walker L, Hobbie L, Risseeuw E, Banks T, Crosby WL, Yang M, Ma H, Estelle M: Identification of an SCF ubiquitin-ligase complex required for auxin response in Arabidopsis thaliana. Genes Dev. 1999, 13: 1678-1691. 10.1101/gad.13.13.1678.
Quint M, Ito H, Zhang W, Gray WM: Characterization of a novel temperature-sensitive allele of the CUL1/AXR6 subunit of SCF ubiquitin-ligases. Plant J. 2005, 43: 371-383. 10.1111/j.1365-313X.2005.02449.x.
Jemo M, Abaidoo RC, Nolte C, Horst WJ: Aluminum resistance of cowpea as affected by phosphorus-deficiency stress. J Plant Physiol. 2007, 164: 442-451. 10.1016/j.jplph.2005.12.010.
Plaxton WC: Plant response to stress: biochemical adaptations to phosphate deficiency. Encyclopedia of plant and crop science. Edited by: Goodman RM. 2004, New-York NY USA, Marcel Dekker Inc., 976-980.
Slaski JJ: Differences in the metabolic responses of root tips of wheat and rye to aluminum stress. Plant Soil. 1994, 167: 165-171. 10.1007/BF01587612.
Singh R, Beriault R, Middaugh J, Hamel R, Chenier D, Appanna VD, Kalyuzhnyi S: Aluminum-tolerant Pseudomonas fluorescens: ROS toxicity and enhanced NADPH production. Extremophiles. 2005, 9: 367-373. 10.1007/s00792-005-0450-7.
Maxwell DP, Wang Y, McIntosh L: The alternative oxidase lowers mitochondrial reactive oxygen production in plant cells. Proc Natl Acad Sci U S A. 1999, 96: 8271-8276. 10.1073/pnas.96.14.8271.
Clifton R, Millar AH, Whelan J: Alternative oxidases in Arabidopsis: a comparative analysis of differential expression in the gene family provides new insights into function of non-phosphorylating bypasses. Biochim Biophys Acta. 2006, 1757: 730-741. 10.1016/j.bbabio.2006.03.009.
Gaume A, Mächler F, Frossard E: Aluminum resistance in two cultivars of Zea mays L.: Root exudation of organic acids and influence of phosphorus nutrition. Plant Soil. 2001, 234: 73-81. 10.1023/A:1010535132296.
Ermolayev V, Weschke W, Manteuffel R: Comparison of Al-induced gene expression in sensitive and tolerant soybean cultivars. J Exp Bot. 2003, 54: 2745-2756. 10.1093/jxb/erg302.
Ma Z, Miyasaka SC: Oxalate exudation by taro in response to Al. Plant Physiol. 1998, 118: 861-865. 10.1104/pp.118.3.861.
Liszkay A, van der ZE, Schopfer P: Production of reactive oxygen intermediates (O(2)(.-), H(2)O(2), and (.)OH) by maize roots and their role in wall loosening and elongation growth. Plant Physiol. 2004, 136: 3114-3123. 10.1104/pp.104.044784.
Schopfer P: Hydroxyl radical-induced cell-wall loosening in vitro and in vivo: implications for the control of elongation growth. Plant J. 2001, 28: 679-688. 10.1046/j.1365-313x.2001.01187.x.
Carmen Cordoba-Pedregosa M, Cordoba F, Villalba JM, Gonzalez-Reyes JA: Zonal changes in ascorbate and hydrogen peroxide contents, peroxidase, and ascorbate-related enzyme activities in onion roots. Plant Physiol. 2003, 131: 697-706. 10.1104/pp.012682.
Cordoba-Pedregosa MC, Villalba JM, Cordoba F, Gonzalez-Reyes JA: Changes in intracellular and apoplastic peroxidase activity, ascorbate redox status, and root elongation induced by enhanced ascorbate content in Allium cepa L. J Exp Bot. 2005, 56: 685-694. 10.1093/jxb/eri051.
Souza IRP, Alves VMC, Parentoni SN, Oliveira AC, Teixeira FF, MacAdam JW, Purcino AC: Change in root apical protein and peroxidase activity in response to aluminum in tolerant and sensitive maize inbred lines. Braz J Plant Physiol. 2002, 14: 219-224. 10.1590/S1677-04202002000300006.
Devi SR, Yamamoto Y, Matsumoto H: An intracellular mechanism of aluminum tolerance associated with high antioxidant status in cultured tobacco cells. J Inorg Biochem. 2003, 97: 59-68. 10.1016/S0162-0134(03)00182-X.
Guo Z, Tan H, Zhu Z, Lu S, Zhou B: Effect of intermediates on ascorbic acid and oxalate biosynthesis of rice and in relation to its stress resistance. Plant Physiol Biochem. 2005, 43: 955-962. 10.1016/j.plaphy.2005.08.007.
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4: 249-264. 10.1093/biostatistics/4.2.249.
Richards KD, Snowden KC, Gardner RC: Wali6 and wali7. Genes induced by aluminum in wheat (Triticum aestivum L.) roots. Plant Physiol. 1994, 105: 1455-1456. 10.1104/pp.105.4.1455.
Delisle G, Champoux M, Houde M: Characterization of oxalate oxidase and cell death in Al-sensitive and tolerant wheat roots. Plant Cell Physiol. 2001, 42: 324-333. 10.1093/pcp/pce041.
Cruz-Ortega R, Cushman JC, Ownby JD: cDNA clones encoding 1,3-beta-glucanase and a fimbrin-like cytoskeletal protein are induced by Al toxicity in wheat roots. Plant Physiol. 1997, 114: 1453-1460. 10.1104/pp.114.4.1453.
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. 10.1007/s00438-006-0169-x.
The authors thank Dr N. Chevrier (UQAM) for bulking up the Atlas66 seed stock. This work was supported by a Natural Sciences and Engineering Research Council of Canada grant (OGP0138557) to M.H.
MH designed the experiment, performed part of the experiments, analyzed the results and prepared the manuscript. AOD performed part of the experiments and their analysis, participated in preparing the figures and critically read the manuscript. All authors have read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Microarray data of candidate stress-associated genes. The differential expression between Al-treated and control (non-treated) samples was analyzed to identify genes that are significantly over-expressed two-fold or more. A first round of selection retained 70 genes for which there is at least two wheat lines with p values > 0.001 AND p values > 0.01 in the other two lines. A second round of selection for which there was at least three lines with p values > 0.001 retained an additional 13 genes (in yellow) for a total of 83 candidate genes. A = Atlas66; B = Bounty; C-T = Century-T; C-S = Century-S. Numbers associated with the line's abbreviation (0, 5 or 50) represents the Al concentration in μM while the following letter indicates the biological replicate sample number (S1 to S3). (XLS 106 KB)
Additional file 2: Microarray data of candidate tolerance-associated genes. The differential expression between Al-treated tolerant and sensitive lines (A 50 – B 5 and C-T 50 – C-S 5) was analyzed to identify genes that are significantly over-expressed two-fold or more in tolerant compared to sensitive lines (ANOVA p values > 0.001 in one pair AND p values > 0.01 in the other pair of wheat lines). Genes on this list that were not differentially expressed between Al-treated and the controls (non-treated) are classified as constitutively expressed. A = Atlas66; B = Bounty; C-T = Century-T; C-S = Century-S. Numbers associated with the cultivar's abbreviation (0, 5 or 50) represents the Al concentration in μM while the following letter indicates the biological replicate sample number (S1 to S3). A: Constitutively expressed candidate genes associated with Al tolerance. B: Up-regulated candidate genes associated with Al tolerance. (XLS 42 KB)
Additional file 3: Primers used for amplification of different transcripts. The transcripts representing the different probesetIDs boxed in Table 1 and Table 2 were used to design unique primers for qRT-PCR. (DOC 24 KB)
About this article
Cite this article
Houde, M., Diallo, A.O. Identification of genes and pathways associated with aluminum stress and tolerance using transcriptome profiling of wheat near-isogenic lines. BMC Genomics 9, 400 (2008). https://doi.org/10.1186/1471-2164-9-400
- Tolerant Cultivar
- Root Growth Inhibition
- Sensitive Cultivar
- Oxalate Oxidase
- C2H2 Zinc Finger