International Journal of Plant Genomics

International Journal of Plant Genomics / 2010 / Article

Research Article | Open Access

Volume 2010 |Article ID 164862 |

Dechassa Duressa, Khairy Soliman, Dongquan Chen, "Identification of Aluminum Responsive Genes in Al-Tolerant Soybean Line PI 416937", International Journal of Plant Genomics, vol. 2010, Article ID 164862, 13 pages, 2010.

Identification of Aluminum Responsive Genes in Al-Tolerant Soybean Line PI 416937

Academic Editor: Akhilesh Kumar Tyagi
Received20 Apr 2010
Revised16 Jul 2010
Accepted07 Aug 2010
Published03 Oct 2010


Soybean is one of the most aluminum (Al) sensitive plants. The complex inheritance of Al tolerance trait has so far undermined breeding efforts to develop Al-tolerant soybeans. Discovering the genetic factors underlying the Al tolerance mechanisms would undoubtedly accelerate the pace of such endeavor. As a first step toward this goal, we analyzed the transcriptome profile in roots of Al-tolerant soybean line PI 416937 comparing Al-treated and untreated control plants using DNA microarrays. Many genes involved in transcription activation, stress response, cell metabolism and signaling were differentially expressed. Patterns of gene expression and mechanisms of Al toxicity and tolerance suggest that Cys2His2 and ADR6 transcription activators, cell wall modifying enzymes, and phytosulfokines growth factor play role in soybean Al tolerance. Our data provide insights into the molecular mechanisms of soybean Al tolerance and will have practical value in genetic improvement of Al tolerance trait.

1. Introduction

Aluminum (Al) toxicity is a major constraint of crop production on acid soils. In view of the fact that 40% of world’s arable land is acidic [1, 2], Al toxicity remains a major hurdle for increasing world food, fiber, and fuel production particularly via expansion of cultivation into acid soils.

Aluminum inflicts a wide range of cellular injuries in plants that ultimately result in reduced root growth, nutrient and water uptake, and productivity [1, 2]. Plants possess some degree of tolerance to Al toxicity that varies among species and genotypes [1, 36]. Al tolerance mechanisms include exclusion and internal detoxification. Al exclusion via rhizosphere Al-organic acid anion complex formation is the most widely documented physiological mechanism of Al tolerance in cultivated and wild plants alike [1, 7]. Root-exuded citrate, malate, and oxalate are the key organic acid anions involved in such mechanism. Genes involved in Al-induced root exudation of malate and citrate have been cloned in wheat [8] and sorghum [5], and their variants are being discovered in several plant species. Internal detoxification mechanisms involve the formation of Al complexes with organic acids, acidic polypeptides, and/or proteins and subsequent sequestration of Al in organelles away from sensitive sites in the cell [9, 10]. The genetic components of the internal detoxification pathways are yet to be elucidated.

In soybean, Al tolerance is a complex trait perhaps involving several genes and pathways [11, 12]. Quantitative trait loci (QTL) mapping in a population derived from Al tolerant PI 416937 and Al sensitive Young has revealed five DNA markers associated with Al tolerance [11]. Most of the alleles were derived from Al-tolerant PI 416937. Other reported soybean Al tolerance genes include phosphoenolpyruvate carboxylase (PEPC), homolog of translationally controlled tumor proteins (TCTPs), inosine 5′-monophosphate dehydrogenases (IMPDHs) [13], aluminum-induced 3-2 (Sali3-2), and aluminum-induced 5-4a (Sali 4-5a) [14]. Ermolayev et al. [13] and Ragland and Soliman [14] used gene expression as a tool to identify the above genes but the techniques used in these experiments were not sensitive enough to detect large number of genes that might be expected from the quantitative nature of soybean Al tolerance trait. The objective of this study was to discover putative Al tolerance genes in Al-tolerant soybean line PI 416937 using DNA microarrays—a robust genome—wide transcript profiling technology. Such an approach was recently employed in wheat [15, 16], maize [17], Arabidopsis [18], and Medicago truncatula [19, 20] to discern the molecular basis of Al tolerance in the respective species.

2. Materials and Methods

2.1. Plant Genotype and Growth Conditions

An Al-tolerant soybean plant introduction (PI 416937) highly characterized for Al response [12, 21] was used in this experiment. Seeds were surface sterilized with 20% household bleach (Clorox) in water for 12 min, rinsed with distilled-deionized water several times, and were germinated in deionized water moistened standard germination paper at 25°C in an incubator for 72 h. Seedlings uniform in tap root length were transferred to black-painted pots filled with approximately 4 L of 800 μM CaCl2 background solution with 10 μM Al added (treated) or no Al added (control) in a Conviron growth chamber (16/8 h light/ dark cycle with respective temp. of 28°C/20°C, photosynthetic photon density of 100 μmol m−2 s−1). The pH of the culture solution was adjusted to 4.3 and maintained at that level for the entire duration of the experiment. After 2, 12, 48, or 72 h of Al treatment 1 cm sections of the primary root tips of approximately 15 plants/pot were harvested, immediately flash frozen in liquid nitrogen, and stored at −70°C for RNA extraction. Three independent replicates were used per treatment.

2.2. RNA Extraction, Microarray Procedure, and Data Analysis

Total RNA was extracted from 100 mg root tissue samples using Qiagen RNeasy plant RNA isolation kit following the manufacturer’s protocol (Qiagen, Inc.). The Affymetrix GeneChip Soybean Genome Array with over 68 000 probe sets, Glycine max L. and wild soybean combined, was used for microarray analysis of the soybean genome for Al tolerance. Three chips were used per treatment. Detailed procedures for RNA labeling and array analysis are described in the Manufacturer’s GeneChip Expression Technical Manual (Affymetrix). Briefly, the quality of total RNA was determined using the RNA 6000 Nano Chip on Agilent BioAnalyzer 2100 prior to double-stranded cDNA synthesis. Total RNA in the amount of 2 μg was used for double-stranded cDNA generation by linear amplification using oligo dT-T7 primer and reverse transcriptase (RT). Subsequently, biotin-labeled cRNA was synthesized by in vitro transcription (IVT) using the ENZO High Yield IVT kit (ENZO). Quality and quantity of cRNA were assessed using the RNA 6000 Nano chip on Agilent BioAnalyzer 2100. Fifteen-microgram cRNA was used for hybridization. Arrays were hybridized overnight at 45°C for 16 h in GeneArray Hybridization Oven 640 (Affymetrix). The next day, arrays were washed and stained in the Fluidics Station 450 (Affymetrix) and scanned by the High Resolution GeneChip Scanner 3000 (Affymetrix).

Gene expression values were determined using theGeneChip Operating Software (GCOS 1.1, Affymetrix). The expression levels were subjected to data query and data mining in Data Mining Tool (DMT). Statistical Analysis of the data was conducted using the software packages ArrayAssist Enterprise together with Pathway Assist (Stratagene/Agilent, Santa Clara, CA). The raw GeneChip files from GeneChip Operating Software (GCOS, Affymetrix, CA) were uploaded, background-subtracted, variance stabilized, and normalized with GC-RMA method [22]. The control group was used as a baseline to calculate the intensity ratio/fold changes of the treatment versus control. The ratio was log2-transformed before further statistical analysis. The -values were obtained by an unpaired -test assuming unequal variance. Significantly upregulated and downregulated genes were annotated using protein databases accessed by blastx at National Center for Biotechnology Information (NCBI).

2.3. Quantitative Real-Time PCR

Quantitative real-time PCR quantification of transcript levels for representative genes [Gma. 20326: F-5′-tcactccccaccttatcgag-3′, R-5′-tcatgtggtggagtgtggtt-3′; Gma. 6948: F-5′-ttatctccggcgaaaacctc-3′, R-5′-tcgtggtgcagcagtttaag-3′; Gma.12326: F-5′-agccactcaaatggttcagc-3′, R-5′-tctccttgtccttctccttcc-3′; Gma. 24062: F-5′-tgccgaaggatcatctcaac-3′, R-5′-cgagggataatggttgatgg-3′; Gma.26937: F-5′-tacccaaaaggcaggcatac-3′, R-5′-ggccgaggtacaaacacatc-3′; Gma.4156: F-5′-tccaatgctgacaagtgctc-3′, R-5′-tagggacactccgtccaatc-3′; Gma.2577: F-5′-acgcctatgaacgtgaaacc-3′, R-5′-aacatcagcggagagcattc-3′] from microarray experiments was conducted using the Roche Diagnostics light Cycler 480 System with SYBR green detection (Roche Diagnostic, Corp) using beta-tubulin gene (beta-tubulin: R-5′-CCATCAAACCTCAAGGAAGC-3′, F-5′-TGCTGTCCTCTTGGACAATG-3′) as internal control. mRNA was isolated from plants grown under similar experimental conditions as in the microarray experiments. mRNA extraction and quality test was as described above. RNA samples were treated with Applied Biosystems Turbo DNA-free DNase (Ambion, Inc.) to remove DNA contamination. Briefly, 2 μl 10x DNase I buffer and 1 μl rDNase I were added to 20 μl RNA sample, and the mix was incubated at 37°C for 30 minutes in water bath. Subsequently, 2 μl resuspended DNase inactivation reagent was added and the samples mixed well and incubated at room temperature for 3 minutes. Samples were then centrifuged at 10 000 g for 1.5 min (Eppendorf centrifuge 5415 D) in 1.6 ml centrifuge tubes and supernatants transferred to fresh tubes.

cDNA was synthesized from 1 μg DNase-treated RNA samples using the Roche Diagnostics Transcriptor First Strand cDNA Synthesis Kit (Roche Diagnostics, Corp) according to manufacturer’s protocol. cDNA concentration and quality was determined using NanoDrop Spectrophotometer brand ND-1000 (NanoDrop Technologies, Inc.). cDNA samples were diluted with nuclease-free water in varying ratios ranging from 1 : 4 to 1 : 10 depending on sample concentration. A total reaction volume of 11 μl comprising 2 μl cDNA sample, 2 μl each of the reverse and forward primers at 0.2 μM concentration, and 5 μl SYBR mix was prepared in 96-well plates (Roche Diagnostics) in two biological and three technical replicates for each gene. A real-time PCR profile of preincubation at 95°C for 5 min, a 45-cycle amplification at 95°C for 10 second, 55°C for 20 second, and 72°C for 20 second, melting at 95°C for 1 min, 65°C for 1 min, and 95°C continuous, and cooling at 40°C for 30 seconds was used to amplify the samples. Negative controls in which cDNA sample was replaced with PCR grade water for each primer pair were included in each run. Sample wells were individually assessed for data quality by evaluating amplification curves and PCR product specificity was verified by melting curve analysis. The expression level of target genes was normalized using in-run beta-tubulin gene as internal control, and transcript concentration ratios were calculated using the ΔΔCT-Method [23]. The change in gene expression levels (fold change) was calculated as treatment to control ratio and compared with results from microarray.

3. Results and Discussion

3.1. Gene Expression in Response to 2-Hour Al Treatment

A total of 38 genes were identified as differentially expressed in the 10 μM Al-treated experimental plants compared to no Al added controls at 2 h post Al treatment (Figure 1). Thirty-four of them were upregulated and 4 were downregulated with a fold change ranging from 3.08 to 32.55 (Table 1).

Unigene IDFold changeAverage ± SDFunctional categoryAnnotation -value

Gma.1866432.55 (down) Stress responseAnionic peroxidase/ oxidative stress
Gma.415229.89 (up) Stress responseTrypsin and protease inhibitor
Gma.1796129.13 (up) Stress responseSoybean oleosin isoform B
Gma.2698427.751 (up) Stress responsePutative protease inhibitor
Gma1500724.75 (up) MetabolismFerredoxin 2 protein
Gma.2985511.98 (up) MetabolismRibulose-1, 5-bisphosphate carboxylase
Gma.1811011.36 (up) UnclassifiedCp12-2 Protein/peptide cross-linking
Gma.269378.29 (down) UnclassifiedUnknown
Gma.109878.03 (up) MetabolismRibulose bisphosphate carboxylase
BE0240057.58 (up) Stress responseGlutaredoxin family protein- /glutathione-dependent reductase
Gma.25777.3 (down) MetabolismHydrolase family protein
Gma.213546.51 (up) Transcription factorNAc1 domain protein Plant development protein
Gma.41566.41 (up) Transcription factorZinc finger protein/transcription factor (CCCH-type family)
Gma.326585.78 (up) UnclassifiedHypothetical protein
Gma.121215.36 (up) UnclassifiedHypothetical protein
Gma.64874.62 (up) UnclassifiedUnknown
Gma.155384.47 (up) Stress responseGlutaredoxin family protein (arsenate reductase)
Gma.283764.43 (down) Stress responseSyringolide-induced protein B13-1-9 hypersensitive response
Gma.124814.39 (up) UnclassifiedHypothetical protein
Gma.42264.36 (up) Stress responseATPP2-A13 protein/wound response
Gma.12484.09 (up) Signalingnod33 protein (putative phosphatase)
BQ6298213.97 (up) Transcription factorMy family transcription factor
Gma.10433.86 (up) UnclassifiedHypothetical protein
Gma.313823.79 (up) Transcription factorBzip transcription factor (bzip 105)
BQ7857793.75 (up) UnclassifiedUnknown
Gma.47103.54 (up) UnclassifiedHypothetical protein
Gma.270153.53 (up) UnclassifiedOcticosapeptide PB1 domain protein
Gma.238493.51 (up) UnclassifiedUnknown
Gma.42163.39 (up) MetabolismEndo-xyloglucan transferase/hydrolase
AW7334633.38 (up) UnclassifiedUnknown
BF5955653.35 (up) MetabolismSDP1 (sugar dependent 1)/ triacylglycerol lipase
CF8073423.29 (up) UnclassifiedHypothetical protein
Gma.323763.27 (up) Transcription factorBLH1 (embryo sac develop arrest 29)
Gma.278373.16 (up) UnclassifiedHypothetical protein
Gma.41493.13 (down) UnclassifiedUnknown
Gma.199173.12 (up) MetabolismCTP synthase/biosynthesis
Gma.345513.08 (up) SignalingMARD1 (mediator of ABA-regulated Dormancy1)
Gma.34293.02 (up) Metabolism2-oxoisovalerate dehydrogenase
Gma.32595††5.25 (up) Stress responseGlutathione s-transferase
Gma.20326†††3.2 (down) Unknown

Significance threshold ( , Fold change >= 3); †† from the 12 h post treatment, ††† from 72 h post treatment; up: upregulated; down: downregulated; -value: the probability that the match between the gene and its annotation has no biological basis. Fold change: absolute value of the ratio of gene expression under Al to gene expression of untreated control. SD: standard deviation.
3.2. Gene Expression in Response to 12- and 72-Hour Al Treatment

At 12 and 72 h post treatment only one gene each showed significant change in expression in response to Al treatment (Figure 1 and Table 1).

3.3. Gene Expression in Response to 48-Hour Al Treatment

The highest number of differentially expressed genes was detected at 48 h post Al treatment (Figures 1 and 2). A total of 542 genes (97.2% upregulated and 2.8% downregulated) were detected. Those exceeding 13-fold changes are presented in Table 2. The marked fold differences observed in the current research are substantially higher in comparison with results obtained by most authors but are comparable to results of [18, 24]. There were two genes in common between the set of genes detected at 2 h and 48 h post treatment (Gma.2577, 7-fold downregulated at 2 h and 8-fold upregulated at 48 h and Gma.26937, 8-fold downregulated at 2 h and 115 upregulated at 48 h). Similar patterns of gene expression were observed in Arabidopsis roots under Al stress with few overlaps between sets of genes detected at 6 h and 48 h post Al treatment [18].

Unigene IDFold changeAverage SDFunctional categoryAnnotation -value

Gma.6089226.57 (up) UnclassifiedUnknown
BM139770176.22 (up) UnclassifiedUnknown
Gma.2586154.32 (up) UnclassifiedUnknown
Gma.1654130.42 (up) TransportCoatomer protein complex subunit 2 protein transporter
Gma.26937115.43 (up) UnclassifiedUnknown
Gma.35222113.94 (up) Stress responseSyringolide-induced protein B13-1-9 defense protein
Gma.2406289.44 (up) UnclassifiedUnknown
Gma.804879.07 (up) UnclassifiedUnknown
Gma.2746672.87 (up) UnclassifiedUnknown
Gma.1232671.14 (up) UnclassifiedUnknown
Gma.694869.75 (up) UnclassifiedUnknown
BU55139765.63 (up) UnclassifiedHypothetical protein
Gma.252361.59 (up) Stress responseSecretory protein (R14 protein soybean-defense protein)
Gma.3560159.51 (up) TransportHeavy-metal transport/detoxification
Gma.694858.09 (up) UnclassifiedUnknown
Gma.1624655.90 (up) UnclassifiedBAP2 (BON associated protein 2)
Gma.3073154.43 (up) UnclassifiedUnknown
BI96787453.43 (up) UnclassifiedUnknown
Gma.452650.59 (up) Transcription factorZinc finger (C2H2 family protein)
Gma.2519147.41 (up) UnclassifiedUnknown
Gma.939747.24 (up) Stress responseSyringolide-induced protein B13-1-9/defense protein
Gma.2546246.34 (up) Transcription factorWRKY19 DNA-binding protein 19
BU57905845.35 (up) MetabolismN-acetyltransferase activity
Gma.2885243.12 (up) MetabolismCytochrome P450
Gma.2751441.91 (up) Stress responseBasic secretory protein/defense protein
Gma.2334741.67 (up) UnclassifiedUnknown
Gma.2207941.46 (up) Stress responseGlutathione s-transferase
Gma.3299440.77 (up) Stress responseAcidic endochitinase (chitinase III-A)
Gma.2774339.93 (up) UnclassifiedUnknown
Gma.162239.56 (up) UnclassifiedHypothetical protein/ABC transporter like
Gma.3675636.68 (up) Transcription factorWRKY17 protein/transcription factor
Gma.562232.77 (up) UnclassifiedUnknown
Gma.2620432.06 (up) MetabolismTransferase/transferase activity
Gma.2723931.70 (up) UnclassifiedUnknown
Gma.3628731.08 (up) MetabolismCarboxylesterase/lipase activity
Gma.2824630.45 (up) UnclassifiedUnknown
Gma.3675330.41 (up) Transcription factorWRKY 30(DNA binding protein
Gma.856529.01 (up) MetabolismHydrolase /xyloglucan endotransglycosylase
Gma.786128.48 (up) UnclassifiedUnknown
Gma.3279028.00 (up) Stress responseBand 7 family protein/hypersensitive inducible reaction protein 1
Gma.769727.48 (down) UnclassifiedHypothetical protein/31 kDa glycoprotein
Gma.3182727.44 (up) Transcription factorWRKY70/DNA-binding protein 70
Gma.2102227.36 (up) UnclassifiedUnknown
Gma.2827327.29 (up) Transcription factorNAC6 /NAC domain protein/apical elongation plant development protein
Gma.3583027.20 (up) SignalingRegulation of gene silencing/calcium-sensor
Gma.826226.92 (up) Stress responseAGc 2-1(oxidative signal-inducible kinase
Gma.1408026.50 (up) TransportSimilar to ATMRP3/multidrug resistance glutathione s-conjugate-exporting ATPase
Gma.2737125.99 (up) UnclassifiedUnknown
BM52373625.34 (up) MetabolismTransferase family protein
Gma.2833025.32 (up) UnclassifiedCalcium-binding protein
Gma.3471724.87 (up) UnclassifiedUnknown
Gma.2668224.58 (up) UnclassifiedUnknown
Gma.1718424.42 (up) TransportGlycolipid-binding/transport protein
Gma.1188824.39 (up) TransportATPP2-B10 (pheloem protien2) carbohydrate binding
Gma.422223.57 (up) UnclassifiedHypothetical protein
Gma.2671223.39 (up) UnclassifiedUnknown
Gma.3332723.16 (up) Transcription factorTranscription factor
Gma.1701923.08 (up) UnclassifiedUnknown
Gma.3317823.07 (up) UnclassifiedPlastocyanin-like domain-containing copper ion binding
Gma 3536422.81 (up) Stress responseFAD-linked oxidoreductase 1/carbohydrate-oxidase
Gma.1583922.63 (down) MetabolismGDSL-motif lipase/hydrolase
Gma.694822.54 (up) UnclassifiedUnknown
Gma.282122.36 (up) Stress responsePR-5 protein (pathogenesis related)
Gma.862821.20 (up) UnclassifiedUnknown protein
BQ47360420.51 (up) UnclassifiedHypothetical protein
Gma.991320.34 (up) UnclassifiedUnknown protein
Gma.14420.01 (up) TransportNodulin protein/transport function
Gma.3409919.67 (up) UnclassifiedHypothetical protein
Gma.430519.63 (up) Stress responseGlutathione s-transferase (GST 15)
Gma.433619.56 (up) UnclassifiedUnknown
Gma.2462519.40 (up) TransportHeavy metal transport/detoxification
Gma.772619.18 (up) SignallingCalcium-binding EF hand family protein
Gma.2653118.95 (up) Transcription factorZinc finger (C3HC4-type ring familyn)
Gma.3471718.94 (up) UnclassifiedUnknown
BE82228218.64 (up) UnclassifiedUnknown
Gma.2706218.22 (up) Transcription factorNAc domain containing protein 2 plant development/apical elongation
Gma.447818.03 (up) UnclassifiedHypothetical protein
Gma.2480718.01 (up) UnclassifiedUnknown
BK000119.117.98 (up) Cell cyclePhytosulfokines 4 precursor/ growth factor cell differentiation, cell proliferation
Gma.2947917.92 (up) UnclassifiedUnknown
Gma.1095617.76 (up) Stress responseSimilar to pathogenesis-related protein (STH-2)
BI96758917.25 (up) UnclassifiedUnknown
Gma.1718417.16 (up) TransportHeavy-metal-associated domain containing protein metal ion transport
Gma.2456117.08 (up) UnclassifiedUnknown
Gma.2173917.05 (up) MetabolismAAA-type ATpase protein/ATpase activity
Gma.1718416.89 (up) TransportGlycolipid-binding protein/glycolipid transport
Gma.436616.62 (up) MetabolismVTc2 (Vitamin C defective 2)/ L-ascorobic-acid biosynthesis
Gma.2640515.37 (up) UnclassifiedUnknown
DQ22298215.09 (up) TransportLipocalin/ fatty acid transport
Gma.1792915.05 (up) MetabolismTransferase family protein
Gma.2151214.99 (up) UnclassifiedUnknown
BE44073214.85 (up) UnclassifiedUnknown
Gma.2966314.64 (up) UnclassifiedUnknown
Gma.3186114.60 (up) UnclassifiedUnknown
Gma.2965514.52 (up) MetabolismCytochromeP50 subfamily B polypeptide 1
Gma.1125714.51 (up) UnclassifiedHypothetical protein exo-1, 3-beta-glucanase precursor
Gma.2664014.48 (up) UnclassifiedUnknown
Gma.2824314.18 (up) UnclassifiedUnknown
CD39441814.12 (up) MetabolismRibulose-1, 5-bisphosphate carboxylase
Gma.152714.09 (down) MetabolismDihydroflavonol reductase (anthocyanin biosynthesis)
Gma.2523414.08 (up) Transcription factorWRKY43 protein
Gma.2805713.62 (up) Transcription factorSali5-4a protein (ADR6)
Gma.848013.58 (up) Stress responseResistance protein LM12
Gma.2875613.53 (up) UnclassifiedUnknown
BM17721813.45 (up) UnclassifiedUnknown
BG55107813.37 (up) UnclassifiedConserved hypothetical protein
Gma.72813.28 (up) UnclassifiedUnknown
Gma.3533213.23 (up) UnclassifiedUnknown
Gma.2671213.17 (up) UnclassifiedUnknown

Significance threshold ( , Fold change >= 3); up: upregulated, down: downregulated; -value: the probability that the match between the gene and its annotation has no biological basis. Fold change: absolute value of the ratio of gene expression under Al to gene expression of untreated control. SD: standard deviation.

The temporal pattern of Al-induced gene expression changes observed in this study diverges from results of other authors. At 12 and 72 h, almost no genes were differentially expressed or detected. The virtually no detection of Al-regulated genes at 12 and 72 h post treatment seems a little odd but it is what is expressed in this soybean genotype at detection thresholds of -value and 3-fold change in an experiment with 3 replications. Gene expression is species and genotype specific [1520] making comparison of results across different studies difficult. The most likely explanation for the 72-hour result is that Al toxicity could have already been neutralized by the 72 h, making differential gene expression unnecessary. The lack of transcriptional response at 12 h, however, is a biological puzzle, and it could represent a very unusual temporal transcriptome response of this soybean genotype to Al stress. Among the few reported Al microarray studies, the results of Kumari et al. [18] in Arabidopsis is the closest to ours with regard to the number of genes detected at early and late time points. They detected 127 genes at 6 h post treatment and 733 genes at 48 h post treatment using a threshold of a 2-fold change whereas we detected 38 genes at 2 h and 542 at 48 h using a 3-fold change.

All of the differentially expressed genes that were functionally annotated by the Genbank nonredundant protein database were grouped into five functional categories based on their putative cellular function. The functional classification showed that stress- and metabolism-related genes constitute the major fractions of Al-regulated genes (Figure 3).

3.4. Quantitative Real-Time PCR Validation of Microarray Expression Levels

The microarray gene expression levels were validated with quantitative real-time PCR for representative genes (Figure 4). In general, the microarray results were in agreement with qRT-PCR but in a few cases quantitative RT-PCR gave higher levels of expression compared to microarray. Such results are obtained by a number of investigators [16, 20, 25]. Detail discussion of factors contributing to the discrepancy between microarray and RT- PCR gene expression levels is covered in [26]. Many authors attribute the phenomenon to the high dynamic range and greater sensitivity of PCR detection. It is worth noting that the gene expression kinetics depicted in Figure 1 shows the efficacy of our experimental design in capturing the full dynamic range of gene expression profiles in the soybean genotype studied. Gene expression peaks at 2 and 48 h suggesting that major savings in microarray experimental expenditure could be realized by limiting sampling to these time points in future experiments.

3.5. Differentially Regulated Genes by Functional Category
3.5.1. Genes Related to Transcription Factors

A number of transcription factors including bZIP, WRKY, MYB, ADR6, and NAc were highly upregulated in the present study (Tables 1 and 2). Members of these families of transcription factors were previously detected under Al stress in several plant species [16, 1820, 27]. Cys2His2-type zinc finger (bZIP) and auxin downregulated (ADR6) factors are particularly interesting from Al tolerance perspective. Cys2His2-type zinc finger (bZIP) protein coregulates molecular response to proton and Al toxicities [28]. It controls the expression of AlMT1—a malate transporter protein that acts in Al exclusion mechanism. In this study, Cys2His2 (Gma.4526, Table 2) was upregulated 51-fold at 48 h post treatment suggesting that malate plays a major role in Al tolerance mechanism of PI 416937 soybean. Earlier physiological study by Silva et al. [29] showed that Al stress increases exudation of both malate and citrate during the first 6 h of exposure to Al in both tolerant and sensitive soybean types. But they concluded that the sustained accumulation and exudation of citrate is mainly responsible for the genotypic differences in Al tolerance. In the present work, 48 h after Al exposure the malate transporter regulator protein was highly expressed in contrast with the observation of Silva et al. [29]. We postulate that Cys2His2 might regulate the expression of other Al tolerance genes in addition to malate transporter. It is also possible that malate biosynthesis becomes a limiting step or malate might indeed play a major role in soybean Al-tolerance contrary to earlier conclusions. ADR6 transcription factors were previously reported as Al tolerance genes [14, 18]. In the present study, ADR6 was highly upregulated (14-fold, Table 2). The plant hormone auxin and ADR6 exhibit opposite behavior in plant roots under Al stress. Al has been shown to inhibit auxin biosynthesis and transport genes as one possible mechanism of its toxicity [18]. On the contrary, ADR6—an auxin downregulated transcription factor is induced under Al stress perhaps mimicking auxin’s role of promoting root growth. These observations suggest that Cys2His2 and ADR6 transcription factors are important modulators of soybean molecular response to Al stress.

3.5.2. Genes Related to Transporters

Transporters, specifically malate (ALMTs) and citrate (MATE) transporters are the first Al tolerance genes cloned in plants and represent the well-characterized Al tolerance mechanism in a wide range of plant species [5, 8]. None of the family members of these two genes were detected in the present study which could be due to constitutive expression. In contrast, an ABC transporter, a multidrug resistance glutathione-S-transferase-exporting ATPase (Gma.14080, Table 2), was upregulated 27-fold at 48 h post treatment in the present study, which could detoxify xenobiotics by transporting glutathione-S-transferase conjugated toxin to the vacuole from sensitive sites in symplast. The involvement of ABC transporters in Al tolerance mechanism is widely documented [15, 18, 30, 31]. Other Al-induced transporters included heavy metal ion transport proteins (Gma.17184 and Gma.24625), lipid transport proteins (DQ222982 and Gma.17184), carbohydrate transport protein (Gma.11888), and coatomer protein complex subunit 2-protien—a polypeptide complex for membrane trafficking (Gma.1654) (Table 2). Heavy metal transport proteins are either located in plasma membrane or subcellular membranes and detoxify heavy metals by exporting metal-ligand complexes out of the cell or by sequestration or compartmentalization of the complex in the vacuole. The internal detoxification mechanism of Al involves formation of Al-organic acid complexes and subsequent transport of the complex by transport proteins to leaf vacuoles in Al hyperaccumulating plants that are adapted to acid soils [1, 9, 10, 32]. Similar mechanism might operate in cultivated plants, and the heavy metal binding proteins upregulated here might function in such pathway.

Lipid and sugar transport proteins are among other transporters detected. Lipid transport proteins transport lipids to cell wall for biosynthesis of cutin layers and surface waxes as a defense mechanism against pathogen attack [33]. They are also induced by abiotic stresses including aluminum [15, 33]. Lipid transport proteins loosen cell wall in a nonhydrolytic mode and enhance cell elongation, a role traditionally attributed to expansins [34]. Aluminum stress inhibits root growth by restricting cell wall extension [1]; hence, there should be a significance to the upregulation of lipid transport proteins under Al stress. Plant sugar transporters have been reported to be induced by pathogen attack and Al stress [18, 35], as is the case in the present study (Gma.11888, Table 2).

3.5.3. Genes Related to Stress Response

Aluminum toxicity has been shown to elicit a wide range of stress-related proteins [19, 20, 36, 37]. In this study, genes known to be responsive to pathogens, oxidative stress, toxins, or Al were classified under this category. Several pathogenesis-related proteins including syringolide-induced protein, acidic endochitinase, PR-5, basic secretory protein, pathogenesis related protein STH-2, and proteinase inhibitors were upregulated at 48 h post Al treatment (Table 2). The confluence between plant molecular response to aluminum toxicity and pathogen infection likely arises from the fact that both cause oxidative stress. However, the role of pathogenesis-related proteins in Al tolerance is equivocal. Overexpression of peroxidase and proteinase inhibitor genes in Arabidopsis did not improve Al tolerance for the transformed plants relative to controls [38]. On the other hand, overexpressing pepper basic pathogenesis-related protein 1 gene in tobacco resulted in enhanced tolerance to heavy metal cadmium and pathogen infection [39].

Other Al-upregulated stress-related genes included carbohydrate oxidase, glutathione-S-transferase, and glutathione-based reductase (Tables 1 and 2). Carbohydrate oxidase and cell wall peroxidases have been reported to provide protection against pathogens by generating hydrogen peroxide from carbohydrate substrates in the apoplast [14]. Hydrogen peroxide has antimicrobial property and also acts as signal molecule for defense genes expression. In the case of aluminum, the activity of these enzymes is correlated with plant Al sensitivity [40, 41]. Glutathione-S-transferase and glutathione-based reductase are the key enzymes of cellular detoxification and antioxidation system [42]. Glutathione reductase catalyses the conversion of oxidized glutathione to reduced form. Glutathione-S-transferase conjugates toxins and electrophilic compounds to reduced glutathione. The glutathione-conjugated toxin is then exported out of the cell or into the vacuole by the ABC transporter proteins discussed above. The concurrent upregulation of glutathione-based reductase, glutathione-S-transferase, and ABC transporter protein suggests that PI 416937 soybean may guard itself against Al by extruding Al out of the cell or by compartmentalization of Al to the vacuole. Yet there are conflicting evidences with respect to the role of the glutathione defense system in plant Al tolerance. Overexpression of glutathione-S-transferase in Arabidopsis thaliana has been shown to enhance plant Al tolerance [38]. On the other hand, Maron et al. [17] found more oxidative stress genes upregulation in Al sensitive cultivar of maize than in Al tolerant cultivar and argue that oxidative stress genes upregulation is a symptom of Al toxicity rather than a tolerance mechanism, an assertion that is supported by findings of [43]. In addition, these genes are responsive to several biotic and abiotic stress factors and, therefore, should not be regarded as major Al tolerance genes while partial role is certainly possible.

3.5.4. Genes Related to Cellular Metabolism

Genes involved in catabolic or biosynthesis of various metabolites were differentially expressed. The most interesting ones from Al tolerance perspective are genes for biosynthesis of ascorbic acid and genes encoding cytochrome P450 and endo-xyloglucan transferases/hydrolases. All were upregulated in the present study, and the last two were previously reported to be upregulated in Arabidopsis [19, 20] and wheat [15, 16] roots under Al stress. Ascorbic acid is an important component of cellular antioxidation system. Oxidative stress is one aspect of Al toxicity, and maintenance of cellular ascrobate homeostasis has been reported to be an essential component of plant Al tolerance [16]. Cytochrome P450 may serve as monooxygenase in the biosynthetic pathways for lignin, defense compounds, hormones, pigments, fatty acids, and signaling molecules or in the detoxification pathway to catalyze the breakdown of numerous endogenous and exogenous toxic compounds [44]. We detected two genes (Gma.28852 upregulated 43-fold and Gma.29655- upregulated 15-fold) which code for cytochrome P450 (Table 2). Gma.28852 encodes protein involved in pathways of ascorbate metabolism, coumarine and phenylpropanoid biosynthesis, and gamma hexachlorohexane degradation. Endoxyloglucan hydrolases are cell wall metabolism enzymes. Members of this family of enzymes have been implicated in Al tolerance [16, 1820]. There is a causal relationship among endoxyloglucan hydrolases, cell wall composition, and Al tolerance. Al induced increases in cell wall pectin and hemicellulose increases plant Al sensitivity [43]. Pectin and hemicellulose form complexes with Al resulting in increased cell wall rigidity and reduced cell extension and growth [27, 43, 45]. Endoxyloglucan hydrolases appear to relax the Al-rigidified cell wall presumably by hydrolyzing the Al-sugar complexes.

3.5.5. Genes Related to Cell Signaling

Perception of stress signal by the cell is the starting point for cascade of events leading to gene expression and change in cell metabolism in response to a stress factor. Aluminum perception and signaling is currently poorly understood. Cell wall-associated receptor kinase (WAK1) was the first Al signaling gene discovered [46], but there is no evidence that demonstrate, WAK1’s major role in Al tolerance [1]. Microarray analyses have shown kinases, phosphates, and EF hand Ca2+ binding proteins as possible components of Al signaling pathway [16, 18]. In the present work, a Ca2+ sensor protein (Gma.35830), calcium-binding EF hand family protein (Gma.7726), oxidative signal kinase (Gma.8262), and a gene for growth factor phytosulfokines precursor (BK0001191) were upregulated 48 h post Al treatment (Table 2). The phytosulfokines growth factor is a novel Al-induced gene, and it is involved in cell proliferation and growth, characteristics that confer Al tolerance.

4. Conclusion

We conducted a transcriptome analysis in Al-tolerant soybean line PI 416937 to identify potential genetic factors underlying Al tolerance trait. Our results uncovered several genes which might potentially have influence on soybean Al tolerance. Among these, two transcription factors, cell wall metabolism enzymes and a cell proliferation gene are particularly interesting from perspective of the physiological and molecular mechanisms of plant Al tolerance. The first transcription factor, Cys2His2 zinc finger protein, coregulates molecular response to proton and aluminum toxicities, the major acid soil stress factors [28]. The second transcription activator, ADR6 is an auxin downregulated gene. Al suppresses auxin biosynthesis and transport in root system which might be one possible mechanism of Al induced root growth inhibition [18]. Conversely, ADR6 is triggered under Al stress probably acting in a parallel pathway to auxin to restore root growth under Al stress. Root cell wall rigidification by Al binding is one principal mechanism of Al toxicity. Cell wall metabolism enzymes and proteins are induced under Al stress and may counteract Al effects on root cell walls. It is increasingly evident that these proteins as well as cell wall pectin and hemicellulose content are important determinants of Al tolerance in cereals [3, 4, 43]. Evidence from this study also implies that cell wall remodeling enzymes and proteins may play role in soybean Al tolerance. Inhibition of cell division and proliferation is another major mechanism of Al toxicity. We identified a novel cell proliferation stimulating gene phytosulfokines growth factor which might reverse this effect of Al. Taken together; our findings provide important insights into the molecular mechanisms of aluminum tolerance in soybean. The genes we identified may guide efforts to improve plant Al tolerance trait.


The authors are grateful to USDA Plant Genetic Resources and Dr. Thomas Carter of North Carolina State University for generous gifts of the seed of soybean genotype used in this study. They thank Dr. Ernst Cebert of Alabama A&M University for his kind assistance in seed multiplication. They would like to express their appreciation to Dr. Michael Crowley, Genomic core facility at University of Alabama at Birmingham for his assistance and training on microarray hybridization in this project. This research was supported in part by NSF/EPSCOR Grant No. 05026. They are also grateful to Alabama Graduate Research Scholars Program for its financial support. Contributed by Agricultural Experiment Station, Alabama A&M University.


  1. L. V. Kochian, O. A. Hoekenga, and M. A. Piñeros, “How do crop plants tolerate acid soils? Mechanisms of aluminum tolerance and phosphorous efficiency,” Annual Review of Plant Biology, vol. 55, pp. 459–493, 2004. View at: Publisher Site | Google Scholar
  2. L. V. Kochian, M. A. Piñeros, and O. A. Hoekenga, “The physiology, genetics and molecular biology of plant aluminum resistance and toxicity,” Plant and Soil, vol. 274, no. 1-2, pp. 175–195, 2005. View at: Publisher Site | Google Scholar
  3. D. Eticha, A. Staß, and W. J. Horst, “Localization of aluminium in the maize root apex: can morin detect cell wall-bound aluminium?” Journal of Experimental Botany, vol. 56, no. 415, pp. 1351–1357, 2005. View at: Publisher Site | Google Scholar
  4. Z. A. K. M. Hossain, H. Koyama, and T. Hara, “Growth and cell wall properties of two wheat cultivars differing in their sensitivity to aluminum stress,” Journal of Plant Physiology, vol. 163, no. 1, pp. 39–47, 2006. View at: Publisher Site | Google Scholar
  5. J. V. Magalhaes, J. Liu, C. T. Guimarães et al., “A gene in the multidrug and toxic compound extrusion (MATE) family confers aluminum tolerance in sorghum,” Nature Genetics, vol. 39, no. 9, pp. 1156–1161, 2007. View at: Publisher Site | Google Scholar
  6. J. L. Yang, Y.Y. Li, Y. J. Zhang et al., “Cell wall polysaccharides are specifically involved in the exclusion of aluminum from the rice root apex,” Plant Physiology, vol. 146, no. 2, pp. 602–611, 2008. View at: Publisher Site | Google Scholar
  7. D. A. Samac and M. Tesfaye, “Plant improvement for tolerance to aluminum in acid soils—a review,” Plant Cell, Tissue and Organ Culture, vol. 75, no. 3, pp. 189–207, 2003. View at: Publisher Site | Google Scholar
  8. T. Sasaki, Y. Yamamoto, B. Ezaki et al., “A wheat gene encoding an aluminum-activated malate transporter,” The Plant Journal, vol. 37, no. 5, pp. 645–653, 2004. View at: Publisher Site | Google Scholar
  9. A. Morita, O. Yanagisawa, S. Takatsu, S. Maeda, and S. Hiradate, “Mechanism for the detoxification of aluminum in roots of tea plant (Camellia sinensis (L.) Kuntze),” Phytochemistry, vol. 69, no. 1, pp. 147–153, 2008. View at: Publisher Site | Google Scholar
  10. T. Watanabe, M. Osaki, H. Yano, and I. Rao, “Internal mechanisms of plant adaptation to aluminum toxicity and phosphorus starvation in three tropical forages,” Journal of Plant Nutrition, vol. 29, no. 7, pp. 1243–1255, 2006. View at: Publisher Site | Google Scholar
  11. C. M. Bianchi-Hall, T. E. Carter Jr., M. A. Bailey et al., “Aluminum tolerance associated with quantitative trait loci derived from soybean PI 416937 in hydroponics,” Crop Science, vol. 40, no. 2, pp. 538–545, 2000. View at: Google Scholar
  12. H. Nian, Z. Yang, H. Huang, X. Yan, and H. Matsumoto, “Citrate secretion induced by aluminum stress may not be a key mechanism responsible for differential aluminum tolerance of some soybean genotypes,” Journal of Plant Nutrition, vol. 27, no. 11, pp. 2047–2066, 2004. View at: Publisher Site | Google Scholar
  13. V. Ermolayev, W. Weschke, and R. Manteuffel, “Comparison of Al-induced gene expression in sensitive and tolerant soybean cultivars,” Journal of Experimental Botany, vol. 54, no. 393, pp. 2745–2756, 2003. View at: Publisher Site | Google Scholar
  14. M. Ragland and K. M. Soliman, “Two genes induced by Al in soybean roots,” Plant Physiology, vol. 114, p. 395, 1997. View at: Google Scholar
  15. P. Guo, G. Bai, B. Carver, R. Li, A. Bernardo, and M. Baum, “Transcriptional analysis between two wheat near-isogenic lines contrasting in aluminum tolerance under aluminum stress,” Molecular Genetics and Genomics, vol. 277, no. 1, pp. 1–12, 2007. View at: Publisher Site | Google Scholar
  16. M. Houde and A. O. Diallo, “Identification of genes and pathways associated with aluminum stress and tolerance using transcriptome profiling of wheat near-isogenic lines,” BMC Genomics, vol. 9, article 400, 2008. View at: Publisher Site | Google Scholar
  17. L. G. Maron, M. Kirst, C. Mao, M. J. Milner, M. Menossi, and L. V. Kochian, “Transcriptional profiling of aluminum toxicity and tolerance responses in maize roots,” New Phytologist, vol. 179, no. 1, pp. 116–128, 2008. View at: Publisher Site | Google Scholar
  18. M. Kumari, G. J. Taylor, and M. K. Deyholos, “Transcriptomic responses to aluminum stress in roots of Arabidopsis thaliana,” Molecular Genetics and Genomics, vol. 279, no. 4, pp. 339–357, 2008. View at: Publisher Site | Google Scholar
  19. D. Chandran, N. Sharopova, S. Ivashuta, J. S. Gantt, K. A. VandenBosch, and D. A. Samac, “Transcriptome profiling identified novel genes associated with aluminum toxicity, resistance and tolerance in Medicago truncatula,” Planta, vol. 228, no. 1, pp. 151–166, 2008. View at: Publisher Site | Google Scholar
  20. D. Chandran, N. Sharopova, K. A. Vandenbosch, D. F. Garvin, and D. A. Samac, “Physiological and molecular characterization of aluminum resistance in Medicago truncatula,” BMC Plant Biology, vol. 8, article 89, 2008. View at: Publisher Site | Google Scholar
  21. I. R. Silva, T. J. Smyth, D. F. Moxley, T. E. Carter, N. S. Allen, and T. W. Rufty, “Aluminum accumulation at nuclei of cells in the root tip. Fluorescence detection using lumogallion and confocal laser scanning microscopy,” Plant Physiology, vol. 123, no. 2, pp. 543–552, 2000. View at: Google Scholar
  22. Z. Wu, R. A. Irizarry, R. Gentleman, F. M. Murillo, and F. Spencer, “A model-based background adjustment for oligonucleotide expression arrays,” Journal of the American Statistical Association, vol. 99, no. 468, pp. 909–917, 2004. View at: Publisher Site | Google Scholar
  23. M. W. Pfaffl, “A new mathematical model for relative quantification in real-time RT-PCR,” Nucleic Acids Research, vol. 29, no. 9, article e45, 2001. View at: Google Scholar
  24. T. L. Maguire, S. Grimmond, A. Forrest, I. Iturbe-Ormaetxe, K. Meksem, and P. Gresshoff, “Tissue-specific gene expression in soybean (Glycine max) detected by cDNA microarray analysis,” Journal of Plant Physiology, vol. 159, no. 12, pp. 1361–1374, 2002. View at: Publisher Site | Google Scholar
  25. Z. Wu, K. M. Soliman, J. J. Bolton, S. Saha, and J. N. Jenkins, “Identification of differentially expressed genes associated with cotton fiber development in a chromosomal substitution line (CS-B22sh),” Functional and Integrative Genomics, vol. 8, no. 2, pp. 165–174, 2008. View at: Publisher Site | Google Scholar
  26. J. S. Morey, J. C. Ryan, and F. M. Van Dolah, “Microarray validation: factors influencing correlation between oligonucleotide microarrays and real-time PCR,” Biological Procedures Online, vol. 8, no. 1, pp. 175–193, 2006. View at: Publisher Site | Google Scholar
  27. J. Zhang, Z. He, H. Tian, G. Zhu, and X. Peng, “Identification of aluminium-responsive genes in rice cultivars with different aluminium sensitivities,” Journal of Experimental Botany, vol. 58, no. 8, pp. 2269–2278, 2007. View at: Publisher Site | Google Scholar
  28. S. Iuchi, H. Koyama, A. Iuchi et al., “Zinc finger protein STOP1 is critical for proton tolerance in Arabidopsis and coregulates a key gene in aluminum tolerance,” Proceedings of the National Academy of Sciences of the United States of America, vol. 104, no. 23, pp. 9900–9905, 2007. View at: Publisher Site | Google Scholar
  29. I. R. Silva, T. J. Smyth, C. D. Raper, T. E. Carter, and T. W. Rufty, “Differential aluminum tolerance in soybean: an evaluation of the role of organic acids,” Physiologia Plantarum, vol. 112, no. 2, pp. 200–210, 2001. View at: Publisher Site | Google Scholar
  30. P. B. Larsen, M. J. B. Geisler, C. A. Jones, K. M. Williams, and J. D. Cancel, “ALS3 encodes a phloem-localized ABC transporter-like protein that is required for aluminum tolerance in Arabidopsis,” The Plant Journal, vol. 41, no. 3, pp. 353–363, 2005. View at: Publisher Site | Google Scholar
  31. T. Sasaki, B. Ezaki, and H. Matsumoto, “A gene encoding multidrug resistance (mdr)-like protein is induced by aluminum and inhibitors of calcium flux in wheat,” Plant and Cell Physiology, vol. 43, no. 2, pp. 177–185, 2002. View at: Google Scholar
  32. J. F. Ma, P. R. Ryan, and E. Delhaize, “Aluminium tolerance in plants and the complexing role of organic acids,” Trends in Plant Science, vol. 6, no. 6, pp. 273–278, 2001. View at: Publisher Site | Google Scholar
  33. J.-C. Kader, “Lipid-transfer proteins: a puzzling family of plant proteins,” Trends in Plant Science, vol. 2, no. 2, pp. 66–70, 1997. View at: Publisher Site | Google Scholar
  34. J. Nieuwland, R. Feron, B. A. H. Huisman et al., “Lipid transfer proteins enhance cell wall extension in tobacco,” The Plant Cell, vol. 17, no. 7, pp. 2009–2019, 2005. View at: Publisher Site | Google Scholar
  35. L. E. Williams, R. Lemoine, and N. Sauer, “Sugar transporters in higher plants—a diversity of roles and complex regulation,” Trends in Plant Science, vol. 5, no. 7, pp. 283–290, 2000. View at: Publisher Site | Google Scholar
  36. M. A. R. Milla, E. D. Butler, A. R. Huete, C. F. Wilson, O. Anderson, and J. P. Gustafson, “Expressed sequence tag-based gene expression analysis under aluminum stress in rye,” Plant Physiology, vol. 130, no. 4, pp. 1706–1716, 2002. View at: Publisher Site | Google Scholar
  37. K. D. Richards, E. J. Schott, Y. K. Sharma, K. R. Davis, and R. C. Gardner, “Aluminum induces oxidative stress genes in Arabidopsis thaliana,” Plant Physiology, vol. 116, no. 1, pp. 409–418, 1998. View at: Google Scholar
  38. B. Ezaki, R. C. Gardner, Y. Ezaki, and H. Matsumoto, “Expression of aluminum-induced genes in transgenic Arabidopsis plants can ameliorate aluminum stress and/or oxidative stress,” Plant Physiology, vol. 122, no. 3, pp. 657–665, 2000. View at: Google Scholar
  39. S. Sarowar, Y. J. Kim, E. N. Kim et al., “Over expression of a pepper basic pathogenesis-related protein 1 gene in tobacco plants enhances resistance to heavy metal and pathogen stresses,” Plant Cell Reporter, vol. 24, pp. 216–224, 2005. View at: Google Scholar
  40. P. R. S. Boscolo, M. Menossi, and R. A. Jorge, “Aluminum-induced oxidative stress in maize,” Phytochemistry, vol. 62, no. 2, pp. 181–189, 2003. View at: Publisher Site | Google Scholar
  41. G. Delisle, M. Champoux, and M. Houde, “Characterization of oxalate oxidase and cell death in Al-sensitive and tolerant wheat roots,” Plant and Cell Physiology, vol. 42, no. 3, pp. 324–333, 2001. View at: Google Scholar
  42. S. Saka, W. Aouacheri, and C. Abdennour, “The capacity of glutathione reductase in cell protection from the toxic effect of heated oils,” Biochimie, vol. 84, no. 7, pp. 661–665, 2002. View at: Publisher Site | Google Scholar
  43. Q. Liu, J. L. Yang, L. S. He, Y. Y. Li, and S. J. Zheng, “Effect of aluminum on cell wall, plasma membrane, antioxidants and root elongation in triticale,” Biologia Plantarum, vol. 52, no. 1, pp. 87–92, 2008. View at: Publisher Site | Google Scholar
  44. M. A. Schuler and D. Werck-Reichhart, “Functional genomics of P450s,” Annual Review of Plant Biology, vol. 54, pp. 629–667, 2003. View at: Publisher Site | Google Scholar
  45. C. Mao, K. Yi, and L. Yang, “Identification of aluminum-regulated genes by cDNA-AFLP in rice (Oryza sativa L.): aluminum-regulated genes for the metabolism of cell wall components,” Journal of Experimental Botany, vol. 55, no. 394, pp. 137–143, 2003. View at: Google Scholar
  46. M. Sivaguru, B. Ezaki, Z.-H. He et al., “Aluminum-induced gene expression and protein localization of a cell wall-associated receptor kinase in Arabidopsis,” Plant Physiology, vol. 132, no. 4, pp. 2256–2266, 2003. View at: Publisher Site | Google Scholar

Copyright © 2010 Dechassa Duressa et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.