- Research article
- Open Access
Early iron-deficiency-induced transcriptional changes in Arabidopsis roots as revealed by microarray analyses
BMC Genomics volume 10, Article number: 147 (2009)
Iron (Fe) is an essential nutrient in plants and animals, and Fe deficiency results in decreased vitality and performance. Due to limited bio-availability of Fe, plants have evolved sophisticated adaptive alterations in development, biochemistry and metabolism that are mainly regulated at the transcriptional level. We have investigated the early transcriptional response to Fe deficiency in roots of the model plant Arabidopsis, using a hydroponic system that permitted removal of Fe from the nutrient solution within seconds and transferring large numbers of plants with little or no mechanical damage to the root systems. We feel that this experimental approach offers significant advantages over previous and recent DNA microarray investigations of the Fe-deficiency response by increasing the resolution of the temporal response and by decreasing non-Fe deficiency-induced transcriptional changes, which are common in microarray analyses.
The expression of sixty genes were changed after 6 h of Fe deficiency and 65% of these were found to overlap with a group of seventy-nine genes that were altered after 24 h. A disproportionally high number of transcripts encoding ion transport proteins were found, which function to increase the Fe concentration and decrease the zinc (Zn) concentration in the cytosol. Analysis of global changes in gene expression revealed that changes in Fe availability were associated with the differential expression of genes that encode transporters with presumed function in uptake and distribution of transition metals other than Fe. It appeared that under conditions of Fe deficiency, the capacity for Zn uptake increased, most probably the result of low specificity of the Fe transporter IRT1 that was induced upon Fe deficiency. The transcriptional regulation of several Zn transports under Fe deficiency led presumably to the homeostatic regulation of the cytosolic concentration of Zn and of other transition metal ions such as Mn to avoid toxicity.
The genomic information obtained from this study gives insights into the rapid transcriptional responses to Fe shortage in plants, and is important for understanding how changes in nutrient availability are translated into responses that help to avoid imbalances in ion distribution. We further identified rapidly induced or repressed genes with potential roles in perception and signaling during Fe deficiency which may aid in the elucidation of these processes.
Iron is essential for numerous biological oxidation-reduction reactions in plants and plays a central catalytic role in all primary metabolic processes, including chlorophyll biosynthesis, photosynthetic and respiratory electron transport, nitrogen assimilation, and a large number of other anabolic and catabolic reactions . In an oxidizing environment, the concentration of soluble Fe3+ is low and the Fe that is present is available to the cell typically in the form of ferric chelates. The mechanism for Fe uptake in Arabidopsis and all other dicotyledonous and non-grass monocotyledonous plants involves the reduction of Fe3+, release of Fe2+ from a chelate, and uptake of the liberated Fe2+ by the cell (strategy I; [2, 3]). In Arabidopsis, the plasma membrane-bound Fe3+-chelate reductase FRO2 catalyzes the reduction of Fe3+ at the cell surface . Homologues of FRO2 have been found in other strategy I plants such as pea , tomato , and cucumber . The Fe2+ is then taken up by an Fe-deficiency-regulated transporter of the ZIP family, IRT1 [8, 9]. In addition to transcriptional regulation, IRT1 and FRO2 are subject to posttranslational regulation, indicating the importance of an exact and rapid response by adjusting the amount of protein appropriate to the Fe demand in order to avoid the accumulation of surplus amounts of Fe [10, 11]. The coordinated functions of IRT1 and FRO2 are principally responsible for Fe uptake in Arabidopsis roots [4, 12, 13].
Research in the regulation of the Fe deficiency response in plants has generated intense interest. The first insights into the nature of this regulation were based on the early discovery of the fer mutant in tomato . The fer mutant (t3238fer) was unable to initiate the typical responses to Fe deficiency, including enhanced extrusion of protons, proliferation of root hairs, and Fe3+-chelate reductase activity. Furthermore, relatively low Fe concentrations (e.g. 0.1 μM) in the growth media were lethal . The FER gene encodes a bHLH transcription factor that is expressed in roots and root tips but not in the leaves or upper hypocotyls . The expression of FER was not greatly altered in tomato plants grown with 0.1 or 10 μM Fe. Thus in tomato, FER is expressed in a root-specific and largely Fe-independent manner.
In Arabidopsis, FIT (AtbHLH029) was discovered as a FER homolog responsible for regulating the Fe deficiency response [16–19]. The genes FRO2 and IRT1 were found to be regulated by FIT, which also controls the expression of a total of 72 genes with known or putative functions in Fe homeostasis . The transcriptional factor FIT was itself regulated by Fe deficiency; however, constitutive expression of FIT under the control of a 35S promoter did not alter the expression patterns of IRT1 and FRO2, nor was the Fe content in the FIT over-expressing plants altered compared to controls . It was concluded that FIT is necessary but not sufficient for regulation of the Fe deficiency response in Arabidopsis. Recently, it was demonstrated that FIT interacts with two bHLH transcription factors, AtbHLH038 and AtbHLH039 . The over-expression of FIT with either AtbHLH038 or AtbHLH039 resulted in the constitutive expression of both IRT1 and FRO2 and the increased accumulation of Fe compared to controls. The function of AtbHLH038 and AtbHLH039 was shown by T-DNA insertion mutations to be redundant . These results have led to the conclusion that an as yet unidentified cis-element is responsible for the Fe-deficiency-associated regulation of IRT1 and FRO2 among other genes.
As in Arabidopsis, the mechanism of perception of Fe deficiency in grasses is still unknown. The cis-acting elements IDE1 and IDE2, which are presumably responsible for Fe-deficiency-inducible expression, have been identified in the promoter region of the barley IDS2 gene . A constitutively expressed transcription factor, IDEF1 belonging to the ABI3/VP1 family, has been shown to specifically bind to IDE1, and over-expression of IDEF1 led to the induction of the bHLH transcription factor OsIRO2 . OsIRO2 expression was also induced by Fe deficiency, and OsIRO2 itself has been shown to bind to the cis-element 5'-CACGTGG-3' . This cis-element was found in the promoters of a number of genes involved in Fe uptake in grasses. In particular, OsIRO2 was shown to be required for the transcriptional regulation of the genes involved in phytosiderophore synthesis and two transcription factor genes, OsNAC4 and a gene encoding an AP2 domain-binding protein. These later two transcription factors appeared to be regulated directly by OsIRO2 .
Recently, a transcription factor belonging to the NAC family, IDEF2, was found to bind specifically to the IDE2 element . NAC transcription factors are a plant-specific family, which have been implicated in developmental responses to biotic and abiotic stress . IDEF2 transcripts were constitutively expressed in shoots and roots and were not greatly altered under Fe deficiency. Reduction of IDEF2 transcript abundance by RNAi under conditions of Fe deficiency resulted in a considerable reduction in transcript abundance of the Fe-phytosiderophore transporter OsYSL2. At present, it appears that in grasses, the Fe-deficiency signal is mediated by at least two constitutively expressed transcription factors that under conditions of Fe deficiency, initiate a transcriptional cascade of events that lead to the adaptive responses to Fe deficiency.
There are similarities in the response pathway to Fe deficiency in Arabidopsis and rice (e.g. IDE1-like sequences and involvement of bHLH transcription factors), and based on our current knowledge, it is reasonable to expect a higher complexity in the response chain in Arabidopsis. The elucidation of regulatory components is rendered more difficult by the lack of knowledge on the temporal pattern of their expression. We present here a detailed time-course analysis of the short-term transcriptional changes induced by Fe deficiency by using DNA microarrays. We demonstrate that the majority of transcriptional changes were established within the first six hours after subjecting the plants to Fe deficiency. The transcriptional changes include most of the previously reported transcriptional changes in response to Fe deficiency but also novel transcripts that are likely to be important in transduction of the Fe deficiency response in Arabidopsis.
The characterization of rapid responses to Fe deficiency at the transcriptome level requires a growth system that permits rapid removal of Fe from the media. Recently, various groups have published DNA microarray experiments using plants grown on agarose plates. Initiation of Fe deficiency was accomplished by transferring to Fe-deficient media, often in the presence of a Fe2+ chelator. Although widely used, the inherent drawbacks of the agar plate system include, inducing physical injuries to the plants during transferring, the roots carry-over of residual amounts of Fe-containing agarose, thus making the time-point of Fe deficiency difficult to determine, and lastly, the use of chelators that may affect the partitioning of metals across the plasma membrane independent of Fe deficiency. In the present study, we have used a well-characterized hydroponic system to analyze the Fe-deficiency response . The system allowed the bulk transfer of a large number of plants with intact roots into Fe-sufficient or Fe-deficient media within seconds and with minimal mechanical damage. Our approach allows for a time-dependent monitoring of Fe deficiency-induced transcriptional changes with a clear-cut onset of the exposure to Fe-free media. We have used this system to generate materials used in the present study.
Time-course of Fe deficiency-induced alterations in transcript abundance
To characterize the early response of transcript abundance to the removal of Fe from the growth media, differential gene expression was determined by DNA microarray analysis following 0, 0.5, 1, 6 and 24 h of Fe-deficient growth. For an initial characterization of the response to Fe deficiency, the genes were extracted from the data base that had a p-value of < 0.05 and a fold change of > 50%. These results are summarized in Venn diagrams for up- and down-regulated transcripts at very early (0 to 1 h) and early (1 to 24 h) time-points following removal of Fe (Figure 1). With only one exception, changes in transcript abundance at 0, 0.5 and 1 h were limited to a single time point, and the differentially expressed genes that were detected immediately after transfer or after 0.5 or 1 h of Fe-deficient growth were not persistent over time. An examination of the over- or under-represented gene ontogeny categories (GO) showed no categories that were indicative of a response to Fe deficiency (Table 1). For these reasons, the genes showing differential expression at the 0, 0.5 and 1 h time points were considered to be Fe-unspecific responses.
Based on the total number of genes analyzed, approximately 0.3% of the transcripts were changed in abundance after 24 h of Fe-deficient growth. Furthermore, approximately 50% of the transcripts that showed a significant change in abundance after 24 h were already differentially expressed by 6 h (Figure 1). The GO categories that contained more than two members and that were over-represented at 6 and 24 h Fe deficiency were associated with cation transport processes and Fe homeostasis (Table 1). The 40 gene transcripts that were unique to the 24 h time-point showed no over-representation in ontogeny categories; whereas the transcripts unique to 6 h showed over-represented gene ontogeny categories for Co/Zn detoxification, nicotianamine metabolism and Cu transport (data not shown). It was, however, clear that the early response to Fe deficiency was initiated between 1 and 6 h of growth under Fe-deficient conditions.
Clustering of Fe-deficiency changes in gene expression
To further characterize the response to Fe deficiency in a time-course manner, all data for transcript abundance at 0, 1, 6 and 24 h Fe deficiency were subjected to a cluster analysis using k-means methods with the algorithm as described previously . The cluster stability was determined using the "benhur" function found in the clusterStab package of the Bioconductor software http://www.bioconductor.org/. Apart from a two-cluster model that was not useful in analyzing the response to Fe deficiency, a six-cluster model generated the smallest number of clusters that gave the most stable results following multiple resorting. Therefore, the six-cluster model was used for further analysis. Clusters one and six showed transcripts with decreasing and increasing abundance during Fe-deficient growth, respectively (Figure 2). The significance of these trends was confirmed by a comparison of the transcript frequency in the top 200 statistically most altered transcripts with the frequency of transcripts in each cluster. The transcripts in cluster one and six were approximately 9- and greater than 60-fold overrepresented, respectively. Cluster four was slightly overrepresented in the top 200 transcripts, while the remaining clusters were slightly underrepresented. Clusters one and six clearly contained the majority of transcripts that were responding to Fe deficiency.
The fine structure of clusters one and six was investigated using heatmap diagrams and hierarchical clustering methods using Euclidian distances (Figure 3). Only transcripts with a p < 0.05 were investigated. The transcript abundance in cluster one decreased during Fe-deficient growth. The decrease was weakly apparent after 1 h Fe-deficient growth and clearly evident at 6 and 24 h. The kinetic behavior of individual transcripts was similar, with a small group of three transcripts showing a more pronounced response (Figure 3, cluster one). Prominent among the down-regulated transcripts in this group were the three ferritin genes, bacterial-type hemoglobin (GLB3) and the ZIP3 transporter. These genes have functions in Fe storage, binding and transport (Table 2). The biological processes of Fe transport and stimulus and the molecular functions of Fe binding and Cu and Zn transport were significantly over-represented in the genes showing a decrease in abundance (Figure 4). Interestingly, an over-representation in transcripts that encode proteins presumably localized in plastids was observed. Not surprisingly however, these transcripts included the ferritin genes but also components of the reductive pentose cycle.
The genes found in cluster six, which contained transcripts that responded with increasing abundance following removal of Fe, could be divided into two subgroups by hierarchical clustering of approximately equal size. The first of these responded to Fe deficiency rapidly and with high increases in abundance already after 6 h of Fe deficiency (Figure 3). A tendency for increased abundance was apparent after 1 h; although, this increase was not significant in all cases. Both of these groups contained a number of transcripts that were known to respond to Fe deficiency, including IRT1, IRT2, FRO3 and some bHLH transcription factors (Table 3). Genes encoding transport proteins and transcription factors were prominently represented in this group.
The second subgroup in cluster six responded more slowly to Fe deficiency and showed a significant increase in transcript abundance at 24 h of Fe-deficient growth with only slight increases after 6 h. The increases after 6 h Fe-deficient growth were in most cases not significant (p > 0.05). The behavior of this subgroup to Fe deficiency was similar to approximately one-third of the transcripts grouped into cluster four (Figure 3). When all clusters with a significantly increased abundance under Fe deficiency were re-clustered, the group belonging to cluster four transcripts was associated with cluster six (Figure 3). The genes responding more slowly to Fe deficiency and belonging to cluster six are listed in Table 4 and those belonging to cluster four in Table 5.
Summarizing the general behavior of transcripts with increased abundance under Fe deficiency, GO categories related to the transport of Fe, Cu, oligopeptides and amino acids were strongly over-represented (Figure 5). Transcripts annotated to response to Zn and to detoxification of Zn and Co were predominant in this group. Furthermore, cellular components associated with the vacuole were over-represented, indicating the important role of the vacuole in the Fe-deficiency response.
Analysis of the differential response to Fe deficiency at the transcript level
The transcripts reported in Table 3 showed the greatest response to Fe deficiency based on fold-induction of any of the transcripts that were identified to be differentially regulated by Fe supply. To confirm the response of a portion of these transcripts to Fe deficiency, the transcript abundance of nine randomly chosen genes from Table 3 was analyzed by semi-quantitative PCR following Fe-deficient growth for 24 h and following Fe re-supply. As a well documented Fe response gene, FRO2 was included in the analysis for reference. In agreement with the published data, FRO2 and IRT1 were strongly up-regulated at 24 h following the removal of Fe from the growth medium (Figure 6). With the re-supply of Fe, transcript abundance for these genes remained elevated over the subsequent 4 h compared to the Fe-sufficient controls. In two cases, At5g67370 and At4g16370, the response to Fe deficiency was generally weak and a clear response was not consistently observed. All other genes tested were strongly up-regulated under conditions of Fe deficiency, thus confirming the results of the microarray experiments. Following the re-supply of Fe, the transcript abundance of these individuals decreased rapidly after 1 h of re-supply, and in most cases no or only weak signals were observed after 4 h of Fe re-supply.
A quantitative validation of the microarray experiments was performed by real-time RT-PCR. The time-course of expression of six genes identified as being strongly induced by Fe deficiency (Table 3) was monitored over the experimental period (Figure 7). The changes in transcript abundance corresponded well with the Gene Chip signals, thus validating the data obtained from the microarray experiments.
Transcript abundance changes rapidly in response to Fe deficiency
The general response of plants to Fe deficiency has been well documented, particularly with regards to Arabidopsis and rice . It is still, however, largely unclear how the Fe status of the plant is perceived and which cellular components might be involved in signal perception. In many of the previous studies of gene expression during the response to Fe deficiency in Arabidopsis, plants were cultivated on agar plates and transferred to plates lacking Fe and often supplemented with a ferrous chelator to induce the response to Fe deficiency. In our opinion, the disadvantages to this approach are three fold. First, the agar plate method of induces more mechanical damage to the roots as a result of transfer, second, the inevitable contamination of a variable amount of Fe-containing agar that was carried over during transfer, and last, the harsh side effects of extraction with ferrous chelates. In order to overcome these adverse effects, we have carried out the experiments using a hydroponics system that has been employed previously . This system allows simultaneous transfer of a large number of plants to Fe-free media with little mechanical damage to the root system. Our current analysis of global gene expression during the early phases of the Fe-deficiency response demonstrated that the changes in transcript abundance relevant to the adaptive Fe stress response are occurring at the 6 h time point and time points following (Table 1 and Figure 1). In preliminary analyses that formed the basis of the investigation presented here, we determined the transcript abundance of representative genes that are known to respond to Fe deficiency (e.g. IRT1, FRO2 and FRO3). The result of that study showed that although changes in transcript abundance of Fe responsive genes could be observed at time-point prior to 6 h, the response at and prior to 4 h was variable and the changes in transcript abundance that were observed did not persist (data not shown). Differential expression was observed immediately following transfer and after 0.5 and 1 h of growth in Fe-depleted hydroponic culture media. However, these changes were restricted to individual time-points and the change were not continuous between the early time-points. In fact, a large overlap in significantly changed gene expression was only observed between 6 and 24 h Fe-deficient growth (Figure 1).
The results from our study correspond well with a previous microarray analysis on the fit mutant by Colangelo and Guerinot . Of the 72 genes that were identified by the previous study as being responsive to Fe deficiency with increased transcript abundance, 22 were also found in the present study. These 22 genes were distributed between early and late responding transcripts.
Rapid transcriptional changes are associated with Fe uptake and distribution
Considering the Fe deficiency-dependent and time-dependent changes in transcript abundance, we have identified 16 genes whose abundance was decreased and 51 whose abundance was increased over the 24 h period. The decreased abundance in transcripts encoding two nodulin-like proteins (At3g25190 and At1g21140) was of particular interest, because these proteins contain a domain belonging to the CCC1-like family. The family of proteins with this domain includes a yeast vacuolar transmembrane protein (CCC1) that has been shown to transport Fe and Mn from the cytosol into vacuole  and VIT1, a CCC1 homolog in Arabidopsis that has been shown to be an influx Fe transporter on the vacuolar membrane . If the function of the two nodulin-like proteins were similar to VIT1, then the decreased abundance of these transcripts would be associated with a decreased capacity of the cell to store Fe in the vacuole. The decreased abundance in the three transcripts encoding ferritins was also easily understood in the context of Fe mobilization in response to deficiency.
Somewhat puzzling was the presence of three transcripts encoding chloroplast proteins in the group of most significantly down-regulated transcripts (Table 2; At3g26650, At2g30570 and At2g05070). At3g26650 encodes one of the two subunits of the chloroplast glyceraldehyde-3-phosphate dehydrogenase (GAPDH). This gene has been shown to be coordinately expressed with the phosphoribulokinase gene (PRK) and CP12-2 genes; however, in our study the three CP12-2 and the PRK genes showed no significant changes in response to Fe deficiency (data not shown). At2g30570 encodes a protein similar to subunit W of the photosystem II reaction center, and At2g05070 encodes Lhcb2.2, a member of the family of light-harvesting chlorophyll a/b-binding (LHC). How these transcripts might be involved in the root response to Fe deficiency is not known.
In contrast to the genes that displayed decreasing transcript abundance in response to Fe deficiency, gene which showed an increased in abundance formed a larger more heterogeneous group. They were subdivided further based on cluster analysis. One of these groups responded more rapidly the Fe deficiency (Table 3) than the other (Tables 4 and 5). As stated above, genes encoding Fe transporters, transcription factors and Fe-containing proteins were prominent among these up-regulated transcripts. With regard to metal transporters, the most rapidly responding transcripts were those encoding the Fe uptake transporters (IRT1 and IRT2 [9, 32]), the Fe3+-chelate reductase FRO3 , and the oligopeptide-Fe transporter OPT3 . The transcript abundance of NRAMP4, a vacuolar Fe efflux carrier , and of ferroportin, a putative Fe efflux carrier of unknown location , was increased at the 24 h time point but not before. Thus, based on the results of the present study, we suggest that the response of Arabidopsis to Fe deficiency was to increase and maintain the capacity for Fe uptake in the cytoplasm.
Recently, an investigation of the transcriptional response of Arabidopsis roots to abiotic stress, including the response to Fe deficiency, has been performed by Dinneny et al. using a cell sorting approach to increase the resolution of the analysis . Upon reprocessing the microarray data published in that study using the methods and criteria concurrent with our own investigation, it was revealed that the number of transcripts which show a significant change was greater than 3100, in contrast to less than 100 transcripts in the present study. While a considerable overlap was also found between the two studies, the individuals identified by Dinneny et al.  generally appeared after a longer exposure to Fe deficiency in the present study. In contrast, the differences in the respective studies can most likely be attributed to the different method of culturing the plants used in the experiments. However, possibly as a result of the high resolution of their systems, the majority of the transcripts identified in the Dinneny et al.  study appears to have no clear connection to the known responses to Fe deficiency and may represent, at least partly, biological signals unrelated to the Fe-deficiency response.
Fe deficiency induces rapid changes in Zn homeostasis
In addition to regulation of Fe transporter abundance under Fe deficiency, transcripts corresponding to Zn transporters were regulated under Fe deficiency. In this case, transcripts for the Zn plasma membrane efflux carrier, MTPc3 , and the vacuolar influx carrier, MTPa2 , were strongly and rapidly increased, while transcripts for the Zn plasma membrane influx carrier, ZIP3, were decreased. At the 24 h time-point, transcripts for the Zn vacuolar influx carrier, ZIF1 , were also increased. The effect of Fe deficiency was not restricted to Fe and Zn. A Cu transporter  was rapidly increased and the Fe-regulated Ni transport was slowly increased . ZIP3 (At2g32270) has been shown to be a Zn transporter that catalyzed Zn uptake when heterologically expressed in yeast . ZIP3 responded to Zn deficiency , but a role for ZIP3 in Fe transport has not been shown.
A summary of the effects of Fe deficiency on ion transport in the cell is presented in Figure 8. Assuming that the regulation of transcript abundance is also reflected in altered ion fluxes in the cell, the response of the cell is to increase cytosolic Fe through increased influx into the cytosol and decreased efflux into the vacuole. On the other hand, rapid regulation of the gene transcripts MTPc3 and MTPa2 indicates changes associated with divalent cation homeostasis and have the effect of decreasing divalent cation influx into and increasing efflux out of the cytosol. The response of the cell to an increase cytosolic Fe concentration under conditions of Fe deficiency could be easily understood, whereas the opposite reaction to a reduction in the cytosolic Zn concentration is more difficult to explain. Perhaps under Fe deficiency the capacity for Zn uptake exceeded the nutritional requirement of the cells, and the cytosolic uptake as well as the efflux capacity was adjusted to establish a new Zn homeostasis. The specificity of the plant metal transporters may be the molecular cause for the deregulation of metal ion homeostasis during Fe deficiency. Recently, it was reported that a correlation between Fe deficiency and increased shoot concentrations of Zn, Mn and Co . Fe-deficient Arabidopsis roots accumulated high concentrations of Co, Mn and Zn compared to nutrient sufficient controls [9, 10, 45]. Thus, a side effect of increasing the Fe uptake capacity is a concomitant uptake other metals via the IRT1 transporter, when the availability of Fe in the media is low.
The re-supply of Fe to Fe-deficient plants did not in all cases lead to a rapid down-regulation of genes which are highly induced by Fe-deficiency. In some cases, e.g. for bHLH101 and bHlH39 re-induction was observed 2 h after Fe resupply after an initial repression by Fe after 1 h. Such an effect has been previously described for IRT1 and FRO2 and attributed to in the regulation of Fe-responsive genes by both local and systemic signals . Our data demonstrate that this regulation is not restricted to IRT1 and FRO2, but also include genes with regulatory function.
Plant Material and Growth Conditions
Arabidopsis thaliana L. cv. Landsberg erecta was grown hydroponically in a growth chamber at a constant relative humidity of 75%, under a 10/14 h light/dark cycle (Osram E40/ES Plantastar, 300 μmol cm-2 s-1) at 21°C (day) or 18°C (night). Two to four seeds were sown in 1.5 ml Eppendorf® tubes that were filled with 1 ml of 0.25% agarose. The tip of the tube was cut off with a hot scalpel and the tubes were placed in black-plastic containers in contact with the nutrient solution that was constantly aerated. Each container contained approximately 45 Eppendorf® tubes. At 20 days following sowing, the agarose plug was removed using a gentle stream of water. The nutrient solution was replaced weekly and at 24 h prior to beginning the experiment. The nutrient solution was composed of KNO3 (3 mM), MgSO4 (0.5 mM), CaCl2 (1.5 mM), K2SO4 (1.5 mM), NaH2PO4 (1.5 mM), H3BO3 (25 μM), MnSO4 (1 μM), ZnSO4 (0.5 μM), (NH4)6Mo7O24 (0.05 μM) CuSO4 (0.3 μM), Fe-EDTA (40 μM) with pH adjusted to 6.0 with KOH .
At 30 days after sowing, the plants were washed briefly in a 0.1 mM EDTA solution and transferred into either Fe-free or Fe-replete (controls) nutrient solutions. Following transfer, the roots were harvested after 0, 0.5, 1, 6 and 24 h into RNALater® (Applied Biosystems). The experiment was initiated 1 h into the light cycle. For each time-point analyzed, approximately 20 plants were harvested. The total processing time was approximately 15 s. To confirm that the plants were responding to Fe deficiency, the Fe3+-chelate reductase activity was determined after 24 h Fe deficiency using the method of Moog et al . Only groups of Fe-deficient plants that showed a greater than 80% increase in the Fe3+-chelate reductase at 24 h were used for RNA isolation. The data reported are averages taken from three independent experiments.
In the re-feeding experiments shown in Figure 6, 30-day-old plants were washed in 0.1 mM EDTA as described above and transferred for 24 h to Fe-deficient media. Subsequently, samples were taken before transfer (0 h control), following 24 h growth in Fe-deficient media (24 h -Fe) and following resupply of Fe at 0, 1, 2 and 4 h. Control plants were treated similarly but were grown in the presence of Fe.
RNA Preparation and DNA-Microarray Hybridization
Total RNA was extracted from roots (RNeasy® Midi Kit, Qiagen). The microarray analysis was performed using the Affymetrix ATH1 Arabidopsis GeneChip. The quality control of RNA, preparation of cRNA and hybridization of oligonucleotide chips was performed by the Affymetrix Gene Expression Service Lab. at the IPMB, Academia Sinica http://ipmb.sinica.edu.tw/affy/.
Preprocessing and Statistical Analysis
Preprocessing and statistical analysis of data was performed using the Bioconductor packages AffylmGUI and Limma (http://www.bioconductor.org/). Arrays were normalized using the Robust P robe-level L inear M odel (PLM) provided in the AffylmGUI package using the default parameters. Over- and under-represented gene ontogeny categories (GO) were determined using the program BiNGO 2.3 with Cytoscape 2.6 (http://www.psb.ugent.be/cbd/papers/BiNGO/) and the website Amigo http://amigo.geneontology.org. The correction for false discovery rates was determined with Bingo software using the method of Benjamini and Hochberg . Cluster analysis was performed using the k-means methods with the algorithm of Hartigan and Wong . The cluster stability was determined using the "benhur" function found in the clusterStab program of Bioconductor http://www.bioconductor.org/.
All microarray data from this study have been deposited in NCBI Gene Expression Omnibus (accession number GSE15189).
cDNA was synthesized from one μg total RNA using an oligo-dT primer. The abundance of cDNA was standardized to the abundance of ubiquitin 10. The amount of cDNA was visualized by ethidium bromide staining following agrarose gel electrophoretic separation. The PCR primers for the genes tested were: AA-Transporter (At5g38820) fwd: ctacgcagctcatctcatgc, rev: cgattgcattggagaaaacc, bHLH039 (At3g56980) fwd: gctttggtgtttctgtttcg, rev: tcacttcgttgtcaccaagc, bHLH101 (At5g04150) fwd: ccgccgtagaaaacttaacg, rev: gacgtagcaatctggacagc, cation flux (At3g58060) fwd: ttcgtcaaacgatcatgagg, rev: tcttgtgcagctctctcagc, CYP71B5 (At3g53280) fwd: agaattggcgagagatgagg, rev: tggaaaaggctttcatcacc, FRO2 (At1g01580) fwd: cctaaatcagctgccgcacatgacg, rev: gttcacaaacattatgctcgtcgggc, IRT1 (At4g19690) fwd: ccccatggtcatggtcatgg, rev: attccaccgcacccgagaag, kelch-repeat (At3g07720) fwd: taatttgcagcagagcaacg, rev: ttgagctcctggtgagtcg, Myb10 (At3g12820) fwd: caaaattggcgatctcttcc, rev: tttccccaaactcctcttcc, Myb72 (At1g56160) fwd: gaaaggaagagcaccatgc, rev: tttcccacatctcaacaatcc, OPT3 (At4g16370) fwd: ggaatgctctttgcttttgc, rev: gggctattagggtggtctcc, ubiqutin (At4g05320) fwd: ggaaaacaattggaggatgg, rev: ttagaaaccaccacgaagacg, unknown (At1g47400) fwd: gtcttttgtcgcaaacttgg, rev: taggaaacaatcacgcagca, unknown (At5g55620) fwd: ccaaatgttcggctactcg, rev: ttgatatgcatgtggagagc, unknown (At5g67370), fwd: ggtttgacccacatttgtcc, rev: tgagtttctcgggtcaaagc, Zn binding (At1g74770) fwd: atgtttgcagggaaaagtgc, rev: gtaggagccacaggttgtgc
Real-time RT PCR
Total RNA was isolated from roots of 20 plants with RNeasy Plant Mini Kit (Qiagen) according to the manufacturer's instructions. Nucleic acid quantity was evaluated by using a NanoDrop ND-1000 UV-Vis Spectrophotometer (NanoDrop Technologies, Wilmington, USA). One μg of total DNase-treated RNA (Turbo DNase, Ambion) was reverse-transcribed using Superscript III Reverse Transcriptase (Invitrogen) with oligo dT primers in a total volume of 20 μL. Real-time quantitative PCR was performed using double-stranded DNA binding dye Syber Green PCR Master mix (Applied Biosystems) in an ABI GeneAmp 7000 Sequence Detection System. Each reaction was run in a triplicate and the melting curves were constructed using Dissociation Curves Software (Applied Biosystems), to ensure that only a single product is amplified. Validation experiments were performed to test the efficiency of the target amplification and the efficiency of the reference amplification. Duplicate CT values were analysed with Microsoft Excel using the comparative CT(ΔCT) method as described by the manufacturer (Applied Biosystems). The amount of target (2-ΔΔCT) was obtained by normalizing to an endogenous reference (α-tubulin, At5g19770) and relative to a calibrator (control tissue). The PCR primers for the genes tested were: FRO3 (At1G23020) fwd: ccatcactcctcaatcacttcca, rev: tccagccttgcttgccata, IRT1 (At4g19690) fwd: caccattcggaatagcgttagg, rev: ccagcggagcatgcattta, IRT2 (At4g19680) fwd: cgtatctccggcgatctca, rev: tgatgcacgggttatcaaagc, MYB72 (At1G56165) fwd: cgagaggtaaccaaatcgcaat, rev: tcttcaaatccgcatgggata, At3G12900 fwd: ccagcctctctctgagcgaat, rev: ctgggtagcctcacacgtaagg, At1G47400 fwd: ggcatatagtgagaatggtggtgat, rev: tcacgcagcaggagcataat
Marschner H: Mineral Nutrition of Higher Plants. 1995, London: Academic Press Ltd, 2
Römheld V, Marschner H: Evidence for a specific uptake system for iron phytosiderophres in roots of grasses. Plant Physiol. 1986, 80 (1): 175-180. 10.1104/pp.80.1.175.
Curie C, Briat JF: Iron transport and signaling in plants. Annu Rev Plant Biol. 2003, 54: 183-206. 10.1146/annurev.arplant.54.031902.135018.
Robinson NJ, Procter CM, Connolly EL, Guerinot ML: A ferric-chelate reductase for iron uptake from soils. Nature. 1999, 397: 694-697. 10.1038/17800.
Waters BM, Blevins DG, Eide DJ: Characterization of FRO1, a pea ferric-chelate reductase involved in root iron acquisition. Plant Physiol. 2002, 129: 85-94. 10.1104/pp.010829.
Li L, Cheng X, Ling HQ: Isolation and characterization of Fe(III)-chelate reductase gene LeFRO1 in tomato. Plant Mol Biol. 2004, 54: 125-136. 10.1023/B:PLAN.0000028790.75090.ab.
Waters BM, Lucena C, Romera FJ, Jester GG, Wynn AN, Rojas CL, Alcántara E, Pérez-Vicente R: Ethylene involvement in the regulation of the H(+)-ATPase CsHA1 gene and of the new isolated ferric reductase CsFRO1 and iron transporter CsIRT1 genes in cucumber plants. Plant Physiol Biochem. 2007, 45: 293-301. 10.1016/j.plaphy.2007.03.011.
Eide D, Broderius M, Fett J, Guerinot ML: A novel iron-regulated metal transporter from plants identified by functional expression in yeast. Proc Natl Acad Sci USA. 1996, 93: 5624-5628. 10.1073/pnas.93.11.5624.
Vert G, Grotz N, Dédaldéchamp F, Gaymard F, Guerinot ML, Briat JF, Curie C: IRT1, an Arabidopsis transporter essential for iron uptake from the soil and for plant growth. Plant Cell. 2002, 14: 1223-1233. 10.1105/tpc.001388.
Connolly EL, Fett JP, Guerinot ML: Expression of the IRT1 metal transporter is controlled by metals at the levels of transcript and protein accumulation. Plant Cell. 2002, 14: 347-1357. 10.1105/tpc.001263.
Connolly EL, Campbell NH, Grotz N, Prichard CL, Guerinot ML: Overexpression of the FRO2 ferric chelate reductase confers tolerance to growth on low iron and uncovers posttranscriptional control. Plant Physiol. 2003, 133: 1102-1110. 10.1104/pp.103.025122.
Henriques R, Jásik J, Klein M, Martinoia E, Feller U, Schell J, Pais MS, Koncz C: Knock-out of Arabidopsis metal transporter gene IRT1 results in iron deficiency accompanied by cell differentiation defects. Plant Mol Biol. 2002, 50: 587-597. 10.1023/A:1019942200164.
Varotto C, Maiwald D, Pesaresi P, Jahns P, Salamini F, Leister D: The metal ion transporter IRT1 is necessary for iron homeostasis and efficient photosynthesis in Arabidopsis thaliana. Plant J. 2002, 31: 589-599. 10.1046/j.1365-313X.2002.01381.x.
Brown JC, Chaney RL, Ambler JE: A new tomato mutant inefficient in the transport of iron. Physiol Plant. 1971, 25: 48-53. 10.1111/j.1399-3054.1971.tb01086.x.
Ling HQ, Bauer P, Bereczky Z, Keller B, Ganal M: The tomato fer gene encoding a bHLH protein controls iron-uptake responses in roots. Proc Natl Acad Sci USA. 2002, 99: 13938-13943. 10.1073/pnas.212448699.
Colangelo EP, Guerinot ML: The essential basic helix-loop-helix protein FIT1 is required for the iron deficiency response. Plant Cell. 2004, 16: 3400-3412. 10.1105/tpc.104.024315.
Jakoby M, Wang HY, Reidt W, Weisshaar B, Bauer P: FRU (BHLH029) is required for induction of iron mobilization genes in Arabidopsis thaliana. FEBS Lett. 2004, 577: 528-534. 10.1016/j.febslet.2004.10.062.
Yuan YX, Zhang J, Wang DW, Ling HQ: AtbHLH29 of Arabidopsis thaliana is a functional ortholog of tomato FER involved in controlling iron acquisition in strategy I plants. Cell Res. 2005, 15: 613-621. 10.1038/sj.cr.7290331.
Bauer P, Ling HQ, Guerinot ML: FIT, the FER-LIKE IRON DEFICIENCY INDUCED TRANSCRIPTION FACTOR in Arabidopsis. Plant Physiol Biochem. 2007, 45: 260-261. 10.1016/j.plaphy.2007.03.006.
Yuan Y, Wu H, Wang N, Li J, Zhao W, Du J, Wang D, Ling HQ: FIT interacts with AtbHLH38 and AtbHLH39 in regulating iron uptake gene expression for iron homeostasis in Arabidopsis. Cell Res. 2008, 18: 385-397. 10.1038/cr.2008.26.
Wang HY, Klatte M, Jakoby M, Bäumlein H, Weisshaar B, Bauer P: Iron deficiency-mediated stress regulation of four subgroup Ib BHLH genes in Arabidopsis thaliana. Planta. 2007, 226: 897-908. 10.1007/s00425-007-0535-x.
Kobayashi T, Nakayama Y, Itai RN, Nakanishi H, Yoshihara T, Mori S, Nishizawa NK: Identification of novel cis-acting elements, IDE1 and IDE2, of the barley IDS2 gene promoter conferring iron-deficiency-inducible, root-specific expression in heterogeneous tobacco plants. Plant J. 2003, 36: 780-793. 10.1046/j.1365-313X.2003.01920.x.
Kobayashi T, Ogo Y, Itai RN, Nakanishi H, Takahashi M, Mori S, Nishizawa NK: The transcription factor IDEF1 regulates the response to and tolerance of iron deficiency in plants. Proc Natl Acad Sci USA. 2007, 104: 19150-19155. 10.1073/pnas.0707010104.
Ogo Y, Itai RN, Nakanishi H, Inoue H, Kobayashi T, Suzuki M, Takahashi M, Mori S, Nishizawa NK: Isolation and characterization of IRO2, a novel iron-regulated bHLH transcription factor in graminaceous plants. J Exp Bot. 2006, 57: 2867-2878. 10.1093/jxb/erl054.
Ogo Y, Itai RN, Nakanishi H, Kobayashi T, Takahashi M, Mori S, Nishizawa NK: The rice bHLH protein OsIRO2 is an essential regulator of the genes involved in Fe uptake under Fe-deficient conditions. Plant J. 2007, 51: 366-377. 10.1111/j.1365-313X.2007.03149.x.
Ogo Y, Kobayashi T, Nakanishi Itai R, Nakanishi H, Kakei Y, Takahashi M, Toki S, Mori S, Nishizawa NK: A novel NAC transcription factor, IDEF2, that recognizes the iron deficiency-responsive element 2 regulates the genes involved in iron homeostasis in plants. J Biol Chem. 2008, 283: 13407-13417. 10.1074/jbc.M708732200.
Uauy C, Distelfeld A, Fahima T, Blechl A, Dubcovsky J: A NAC Gene regulating senescence improves grain protein, zinc, and iron content in wheat. Science. 2006, 314: 12398-12301. 10.1126/science.1133649.
Thimm O, Essigmann B, Kloska S, Altmann T, Buckhout TJ: Response of Arabidopsis to iron deficiency stress as revealed by microarray analysis. Plant Physiol. 2001, 127: 1030-1043. 10.1104/pp.010191.
Hartigan A, Wong MA: A k-means clustering algorithm. Appl Statist. 1979, 28: 100-108. 10.2307/2346830.
Li L, Chen OS, McVey Ward D, Kaplan J: CCC1 is a transporter that mediates vacuolar iron storage in yeast. J Biol Chem. 2001, 276: 29515-29519. 10.1074/jbc.M103944200.
Kim SA, Punshon T, Lanzirotti A, Li L, Alonso JM, Ecker JR, Kaplan J, Guerinot ML: Localization of iron in Arabidopsis seed requires the vacuolar membrane transporter VIT1. Science. 2006, 314: 1295-1298. 10.1126/science.1132563.
Vert G, Briat JF, Curie C: Arabidopsis IRT2 gene encodes a root-periphery iron transporter. Plant J. 2001, 26: 181-189. 10.1046/j.1365-313x.2001.01018.x.
Mukherjee I, Campbell NH, Ash JS, Connolly EL: Expression profiling of the Arabidopsis ferric chelate reductase (FRO) gene family reveals differential regulation by iron and copper. Planta. 2006, 223: 1178-1190. 10.1007/s00425-005-0165-0.
Stacey MG, Patel A, McClain WE, Mathieu M, Remley M, Rogers EE, Gassmann W, Blevins DG, Stacey G: The Arabidopsis AtOPT3 protein functions in metal homeostasis and movement of iron to developing seeds. Plant Physiol. 2008, 146: 589-601. 10.1104/pp.107.108183.
Lanquar V, Lelièvre F, Bolte S, Hamès C, Alcon C, Neumann D, Vansuyt G, Curie C, Schröder A, Krämer U, Barbier-Brygoo H, Thomine S: Mobilization of vacuolar iron by AtNRAMP3 and AtNRAMP4 is essential for seed germination on low iron. EMBO J. 2005, 24: 4041-4051. 10.1038/sj.emboj.7600864.
Wintz H, Fox T, Wu YY, Feng V, Chen W, Chang HS, Zhu T, Vulpe C: Expression profiles of Arabidopsis thaliana in mineral deficiencies reveal novel transporters involved in metal homeostasis. J Biol Chem. 2003, 278: 47644-47653. 10.1074/jbc.M309338200.
Dinneny JR, Long TA, Wang JY, Jung JW, Mace D, Pointer S, Barron C, Brady SM, Schiefelbein J, Benfey PN: Cell identity mediates the response of Arabidopsis roots to abiotic stress. Science. 2008, 320: 942-945. 10.1126/science.1153795.
Mäser P, Thomine S, Schroeder JI, Ward JM, Hirschi K, Sze H, Talke IN, Amtmann A, Maathuis FJ, Sanders D, Harper JF, Tchieu J, Gribskov M, Persans MW, Salt DE, Kim SA, Guerinot ML: Phylogenetic Relationships within Cation Transporter Families of Arabidopsis. Plant Physiol. 2001, 126: 1646-1667. 10.1104/pp.126.4.1646.
Arrivault S, Senger T, Krämer U: The Arabidopsis metal tolerance protein AtMTP3 maintains metal homeostasis by mediating Zn exclusion from the shoot under Fe deficiency and Zn oversupply. Plant J. 2006, 46: 861-879. 10.1111/j.1365-313X.2006.02746.x.
Haydon MJ, Cobbett CS: A novel major facilitator superfamily protein at the tonoplast influences zinc tolerance and accumulation in Arabidopsis. Plant Physiol. 2007, 143: 1705-1719. 10.1104/pp.106.092015.
Sancenón V, Puig S, Mira H, Thiele DJ, Peñarrubia L: Identification of a copper transporter family in Arabidopsis thaliana. Plant Mol Biol. 2003, 51: 577-587. 10.1023/A:1022345507112.
Schaaf G, Honsbein A, Meda AR, Kirchner S, Wipf D, von Wirén N: AtIREG2 encodes a tonoplast transport protein involved in iron-dependent nickel detoxification in Arabidopsis thaliana roots. J Biol Chem. 2006, 281: 25532-25540. 10.1074/jbc.M601062200.
Grotz N, Fox T, Connolly E, Park W, Guerinot ML, Eide D: Identification of a family of zinc transporter genes from Arabidopsis that respond to zinc deficiency. Proc Natl Acad Sci USA. 1998, 95: 7220-7224. 10.1073/pnas.95.12.7220.
Baxter IR, Vitek O, Lahner B, Muthukumar B, Borghi M, Morrissey J, Guerinot ML, Salt DE: The leaf ionome as a multivariable system to detect a plant's physiological status. Proc Natl Acad Sci USA. 2008, 105: 2081-12086. 10.1073/pnas.0804175105.
Korshunova YO, Eide D, Clark WG, Guerinot ML, Pakrasi HB: The IRT1 protein from Arabidopsis thaliana is a metal transporter with a broad substrate range. Plant Mol Biol. 1999, 40: 37-44. 10.1023/A:1026438615520.
Vert GA, Briat JF, Curie C: Dual regulation of the Arabidopsis high-affinity root iron uptake system by local and long-distance signals. Plant Physiol. 2003, 132: 796-804. 10.1104/pp.102.016089.
Schmidt W: Root-mediated ferric reduction: responses to iron deficiency, exogenously induced changes in hormonal balance and inhibition of protein synthesis. J Exp Bot. 1994, 45: 725-731. 10.1093/jxb/45.6.725.
Moog PR, Kooij van der TA, Brüggemann W, Schiefelbein JW, Kuiper PJ: Responses to iron deficiency in Arabidopsis thaliana: The Turbo iron reductase does not depend on the formation of root hairs and transfer cells. Planta. 1995, 195: 505-513. 10.1007/BF00195707.
Wettenhall JM, Simpson KM, Satterley K, Smyth GK: A ylmGUI: a graphical user interface for linear modeling of single channel microarray data. Bioinformatics. 2006, 22: 897-899. 10.1093/bioinformatics/btl025.
Maere S, Heymans K, Kuiper M: BiNGO: a Cytoscape plugin to assess overrepresentation of Gene Ontology categories in biological networks. Bioinformatics. 2005, 21: 3448-3449. 10.1093/bioinformatics/bti551.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B. 1995, 57: 289-300.
TJB led the experimental design of the present study, carried out the data analysis and drafted the manuscript. TJWY contributed to the design of the study, prepared RNA samples, conducted experiments to validate the microarray data, and was involved in editing of the manuscript. WS participated in the data analysis, strategic planning of the work, and writing of the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Buckhout, T.J., Yang, T.J. & Schmidt, W. Early iron-deficiency-induced transcriptional changes in Arabidopsis roots as revealed by microarray analyses. BMC Genomics 10, 147 (2009). https://doi.org/10.1186/1471-2164-10-147
- Transcript Abundance
- Arabidopsis Root
- bHLH Transcription Factor
- Ferritin Gene
- Chelate Reductase