Microarray analysis of iron deficiency chlorosis in near-isogenic soybean lines
© O'Rourke et al. 2007
Received: 23 March 2007
Accepted: 21 December 2007
Published: 21 December 2007
Skip to main content
© O'Rourke et al. 2007
Received: 23 March 2007
Accepted: 21 December 2007
Published: 21 December 2007
Iron is one of fourteen mineral elements required for proper plant growth and development of soybean (Glycine max L. Merr.). Soybeans grown on calcareous soils, which are prevalent in the upper Midwest of the United States, often exhibit symptoms indicative of iron deficiency chlorosis (IDC). Yield loss has a positive linear correlation with increasing severity of chlorotic symptoms. As soybean is an important agronomic crop, it is essential to understand the genetics and physiology of traits affecting plant yield. Soybean cultivars vary greatly in their ability to respond successfully to iron deficiency stress. Microarray analyses permit the identification of genes and physiological processes involved in soybean's response to iron stress.
RNA isolated from the roots of two near isogenic lines, which differ in iron efficiency, PI 548533 (Clark; iron efficient) and PI 547430 (IsoClark; iron inefficient), were compared on a spotted microarray slide containing 9,728 cDNAs from root specific EST libraries. A comparison of RNA transcripts isolated from plants grown under iron limiting hydroponic conditions for two weeks revealed 43 genes as differentially expressed. A single linkage clustering analysis of these 43 genes showed 57% of them possessed high sequence similarity to known stress induced genes. A control experiment comparing plants grown under adequate iron hydroponic conditions showed no differences in gene expression between the two near isogenic lines. Expression levels of a subset of the differentially expressed genes were also compared by real time reverse transcriptase PCR (RT-PCR). The RT-PCR experiments confirmed differential expression between the iron efficient and iron inefficient plants for 9 of 10 randomly chosen genes examined. To gain further insight into the iron physiological status of the plants, the root iron reductase activity was measured in both iron efficient and inefficient genotypes for plants grown under iron sufficient and iron limited conditions. Iron inefficient plants failed to respond to decreased iron availability with increased activity of Fe reductase.
These experiments have identified genes involved in the soybean iron deficiency chlorosis response under iron deficient conditions. Single linkage cluster analysis suggests iron limited soybeans mount a general stress response as well as a specialized iron deficiency stress response. Root membrane bound reductase capacity is often correlated with iron efficiency. Under iron-limited conditions, the iron efficient plant had high root bound membrane reductase capacity while the iron inefficient plants reductase levels remained low, further limiting iron uptake through the root. Many of the genes up-regulated in the iron inefficient NIL are involved in known stress induced pathways. The most striking response of the iron inefficient genotype to iron deficiency stress was the induction of a profusion of signaling and regulatory genes, presumably in an attempt to establish and maintain cellular homeostasis. Genes were up-regulated that point toward an increased transport of molecules through membranes. Genes associated with reactive oxidative species and an ROS-defensive enzyme were also induced. The up-regulation of genes involved in DNA repair and RNA stability reflect the inhospitable cellular environment resulting from iron deficiency stress. Other genes were induced that are involved in protein and lipid catabolism; perhaps as an effort to maintain carbon flow and scavenge energy. The under-expression of a key glycolitic gene may result in the iron-inefficient genotype being energetically challenged to maintain a stable cellular environment. These experiments have identified candidate genes and processes for further experimentation to increase our understanding of soybeans' response to iron deficiency stress.
The ability of iron (Fe) to serve as an electron acceptor makes Fe a valuable cofactor in a variety of plant processes including photosynthesis, respiration, and seed development. In the soil matrix, Fe exists in one of two forms, Fe2+or Fe3+. However, many environmental conditions, including the high pH of calcareous soils, can result in little Fe2+ availability [1–4]. To survive in iron limiting environments, plants have evolved two iron uptake strategies, Strategy I and II . Dicot species, including soybean, utilize the Strategy I mechanism to take up the Fe2+ ion. Strategy I plants utilize an ATPase to secrete protons from the roots to acidify the rhizosphere [1, 6–9]. This acidification aids in the release of Fe from chelating agents in the soil. A root membrane reductase reduces the prevalent Fe3+ion to the biologically usable Fe2+ ion. This Fe2+ can then be transported into the roots of the plant where it is available for use in various cellular processes. For strategy I plants, the iron reduction by plant roots has been identified as the rate-limiting step in iron deficiency . Strategy II plants, monocot species, release phytosiderophores from the roots that chelate Fe3+ ions. The entire phytosiderophore iron complex is then transported into the root system of the plant.
Complex genetic and environmental interactions have made soybean IDC an extremely difficult trait to study in field trials [11, 12]. Low Fe availability exacerbates chlorosis levels in many cultivars. This is true in the calcareous soils prevalent in the upper U.S. Midwest farmlands . As plants are subjected to Fe deficiency stress, they respond in a characteristic manner. Developing trifoliates exhibit interveinal chlorosis, growth is stunted, and yield is reduced. Yield reduction has a positive linear correlation with increasing chlorosis levels . To minimize the environmental effect on the plant phenotype, visual phenotypic studies have been conducted with plants grown in a nutrient solution hydroponics system. The hydroponics experiments identified the same QTLs identified in field grown trials  making this a viable system in which to study the effects of IDC on soybean while minimizing environmental effects. The comparison of expression profiles, via utilization of cDNA microarrays, of RNA from Fe efficient and inefficient soybean near isogenic lines (NILs) grown under Fe limited hydroponic conditions will identify differentially expressed transcripts related to iron stress. This will provide clues to the physiological differences between iron efficient and inefficient cultivars
Genes differentially expressed between iron near-isogenic lines that cluster with other stress induced genes.
Associated TIGR TC
TBLASTX UniRef DB Annotation
UniRef Blast E-Value
1Fe, 1PS, 1ST
Serine/Threonine Protein Kinase
3Fe, 118PS, 107ST
Peroxissomal Copper Containing Oxidase
1Fe, 2PS, 1ST
Initiation Factor eIF-4 gamma
3Fe, 118PS, 107ST
Zinc Finger Protein, Cys3His
Zinc Finger Protein, H2
2Fe, 6PS, 11ST
Zinc Finger Protein
2Fe, 6PS, 11ST
Phosphatase type 2C/
1Fe, 5PS, 6ST/1Fe, 9PS, 5ST
Ethylene Responsive Transcription Factor
2Fe, 8PS, 33ST
Ethylene Responsive Transcription Factor/Ubiquitin
2Fe, 8PS, 33ST/2Fe
Response Regulator Protein (ARR)
1Fe, 1PS, 7ST
RuBisCo Activase Protein
1Fe, 1PS, 1ST
2Fe, 31PS, 15ST
1Fe, 22PS, 5ST
Dihydroflavonol 4 Reductase
1Fe, 4PS, 3ST
Heat Shock Protein Hsp70
1Fe, 2ST/1Fe, 2PS, 1ST
RNA Binding Protein
1Fe, 1PS, 1ST
Aquaporin Protein PIP1
2Fe, 31PS, 15ST
No UniRef Hit E < 10E-4
1Fe, 2PS, 2ST
Map Protein Kinase
3Fe, 118PS, 107ST
Genes differentially expressed between near-isogenic lines that did not cluster with other iron or stress induced genes.
Associated TIGR TC
TBLASTX UniProt DB Annotation
UniProt Blast E-Value
Transcription Factor DP1
26S Proteasome non-ATPase Regulatory Subunit
Putative arsA Homolog hASNA-1
Large Polyprotein 2
Plant Specific Membrane Protein
No UniRef Hit E < 10-4
No UniRef Hit E < 10-4
Ubiquitin Conjugating Enzyme
Zinc Ring Finger Protein
TIR-NBS-LRR-TIR Type Disease Resistance Protein
Nonsense Mediated Decay Protein UPF3
No UniRef Hit E < 10E-4
Phagocytosis and Cell Motility Protein
Zinc Ring Finger Protein
As controls, the NILs were also grown in Fe sufficient hydroponics solutions (100 uM Fe(NO3)3) and analyzed on cDNA arrays containing the original 9,728 genes from root specific cDNA libraries examined plus an additional 9,272 genes from seed coat, seedling, cotyledon, flower, and pod cDNA libraries for a more global transcript analysis. An analysis of three biological replicates, with two technical replicates apiece for a total of six hybridizations, showed no genes with consistent differential expression between the NILs under Fe sufficient conditions. Thus, the differential expression seen under Fe deficient conditions is likely a result of the differential response of the NILs to the Fe limited environment rather than inherent genetic differences between the NILs .
Real Time PCR results confirming differential expression identified by microarray analysis
Fold Change Identified by Real Time RT PCR
Fe Efficient Standard Error
Fe Inefficient Standard Error
To determine a probable function of the 43 differentially expressed genes, the GenBank accession of the Expressed Sequence Tag (EST) for the gene was queried against The Institute for Genomic Research (TIGR) database soybean gene index (Version 12.0)  to identify the tentative consensus (TC) sequence containing the respective EST. The TC sequence was compared to the UniProt protein database (February 2006)  using BLASTX  and an E-value cutoff of E < 10-4, to assign a putative function (Tables 1 and 2). Eight of the forty-three sequences examined had no homology to the UniProt protein database. Therefore, the eight individual EST sequences (Gm-c1028-8183, Gm-c1028-8336, Gm-c1009-2578, Gm-c1028-4530, Gm-c1028-1850, Gm-c1028-5360, Gm-c1028-963, and Gm-c1028-4123) were queried against a database of available Soybean Whole Genome Shotgun (WGS) using megaBLAST BLASTN with an E-value cutoff of E < 10-100 to identify genomic sequence that could extend the EST sequence. Identical sequence reads which were at least 500 nucleotides in length and shared 100% nucleotide identity to the EST were assembled into a multiple sequence alignment with the EST. If any of the identified sequences extended the ends of the EST a new consensus was generated for the EST. The new consensus was then compared to the UniProt database by BLASTX with an E-value cutoff of E < 10-4 to assign a putative function.
Genes known to be involved in the Fe deficiency response have been identified and characterized in model organisms such as Arabidopsis thaliana. To determine if homologs of these genes were present on the soybean cDNA array, 33 members of six Arabidopsis gene families known to be involved in Fe uptake and homeostasis (IRT, FRO, FRD, FIT, NRAMP and YSL) were compared to the soybean EST database by BLASTN comparison (E < 10-4). Soybean EST sequences belonging to the Gm-r1021 and Gm-r1083 libraries, and thus putatively represented on the cDNA array, were identified. The soybean sequences were then compared (BLASTN) back to the Arabidopsis genome to determine if they were the reciprocal best match to the original Arabidopsis iron genes and likely functional orthologs. This approach demonstrated that only one soybean ortholog of an Arabidopsis iron uptake gene was represented on the array. Soybean EST Gm-r1083-2131 is the homolog of Arabidopsis Yellow Stripe-Like 6. However, this gene was not differentially expressed in the microarray experiment.
A single linkage cluster analysis  was performed to identify any Fe induced genes with sequence homology (E < 10-4) to other stress induced genes. Twenty-four of the 43 Fe deficiency induced genes clustered with known stress-induced genes (Table 1). Most clusters contain only one Fe induced gene and a number of other stress induced genes. However, one cluster was composed of only two genes (Gm-c1009-2900 and Gm-c1028-6890) which showed homology to each other and were differentially expressed under Fe deficient conditions, but show no significant homology to other stress induced genes. The remaining nineteen Fe deficiency induced genes showed no sequence homology to known stress induced genes, nor to the other Fe deficiency induced genes identified by the microarray experiment (Table 2).
In calcareous soils iron-inefficient soybean genotypes often display symptoms of iron deficiency stress (interveinal chlorosis and reduced yield). These symptoms are exacerbated by the cool wet conditions prevalent in early spring. Under field conditions, if the young soybean plant survives the initial iron stress, the plant continues to grow, albeit slowly, and eventually, as the plant matures and the environmental conditions change, the phenotypic effects of iron stress disappear . Soybean cultivars differ in their ability to respond successfully to iron stress. Results of this study have provided clues to understand some of the physiological differences between iron-efficient cultivars and iron-inefficient cultivars.
In this study, NILs developed especially for their iron deficiency response by the USDA , were used in an established hydroponics system to compare gene expression profiles between iron efficient (Clark) and iron inefficient (IsoClark) NILs. The NILs are phenotypically identical, except in their chlorotic response under iron stress conditions. Clark remains a healthy green under iron deficient conditions while IsoClark exhibits severe interveinal chlorosis.
Growing the NILs in an established hydroponics system allowed for a comparison of gene expression profiles of the roots of iron efficient (Clark) and inefficient (IsoClark) plants to identify differentially expressed genes between the NILs to better understand the physiological responses of soybean to iron stress. Most changes in gene expression identified under iron-limited conditions are negated upon the re-supply of iron to the system , thus confirming these genes as induced by iron deficiency. This comparison allows us to confidently report these forty-three genes as differentially expressed between the NILs in response to iron deficiency.
Induction of the Fe(III) chelate reductase is a key iron stress response in Strategy I plants . Without reductase activity, the available Fe+2 for uptake of iron into the root is extremely limited. Because the gene encoding iron reductase (FRO2) was not found on the cDNA array used in this study, we conducted an iron reductase assay to assess this uniquely-Strategy I response in both the iron efficient and inefficient genotypes. The iron inefficient genotype used in this study failed to respond to iron deficiency stress by induction of increased ferric reductase activity. However, the iron-efficient genotype responded to reduced iron availability by increasing its ability to reduce Fe+3 to the usable Fe+2. Reduction of ferric iron by the roots is considered to be a limiting factor in successful response to reduced iron availability . The lack of induction of the iron reductase in the inefficient isoline is likely a major factor contributing to the severe iron deficiency stress symptoms observed, relative to the iron efficient genotype.
With the advent of microarray technology, researchers can now identify a broad range of genes that work in concert to protect the plant from abiotic and biotic stresses. While some genes may be specific to a particular pathogen, stress, or plant species, others may be part of a general stress response shared across multiple plant species or multiple stresses. We developed an in-house sequence database that contains genes identified from the literature that are significantly differentially regulated in response to abiotic or biotic stresses. Some of the sequences are differentially expressed in response to pathogen attack [24, 25] while the majority, are differentially expressed in response to a variety of abiotic stresses including oxygen deprivation , drought , salt stress [28, 29], Fe deficiency , oxidative stress , phosphate deficiency and others [32–34]. Included in this database are the 43 genes identified in our experiments as differentially expressed in response to limited Fe.
A single linkage cluster analysis  was performed to determine if the genes identified as differentially expressed in our microarray experiment exhibited significant (10E-4) sequence similarity to genes differentially expressed under other abiotic stress conditions. Twenty-four of the 43 identified genes showed significant sequence similarity to other genes whose expression levels are altered by some form of abiotic stress (Table 1). The remaining 19 genes showed no sequence homology to known stress induced genes (Table 2). These 19 genes may be unique to soybeans' iron response. Three sequences show no sequence homology to any of the genes characterized in the UniProt database, or to other genes identified under iron limited conditions. The unique sequence of these three genes suggests they may be unique to legumes. The two groupings of genes (Tables 1 and 2) identified under Fe limited conditions suggest both a universal stress response and an Fe specific stress response are induced upon Fe deficient conditions.
Plants respond to iron stress through an impressive number of metabolic adaptations and adjustments. The iron-inefficient soybean isoline used in this study failed to respond to reduced iron availability by increased activity of Fe(III) chelate reductase. Thus, the reduced availability of the iron in the growth medium created a severe iron stress for the inefficient plants.
The most striking response of the inefficient isoline to iron stress was the dramatic increase in transcripts of genes involved in signaling and hormonal regulation. Increased signaling is likely an attempt on the part of the stressed plant to maintain metabolic homeostasis in a decreasingly sustainable environment. For example, MAP kinase and a SNARE protein are well known signaling proteins that were induced in the inefficient line. In addition, RNA mediating genes for RNA methyltransferase and an RNA binding protein were also induced upon iron stress, as were several DNA-binding zinc finger protein genes and ethylene receptors.
Ethylene is a signaling molecule often associated with root hair development, pathogen infection, wounding and other abiotic stresses  including iron stress [36–38]. In this study, transcripts encoding an ethylene receptor protein and two ethylene responsive transcription factors were up-regulated. These transcription factors are one of the largest families of transcription factors and can be induced by abiotic stresses . A MAP kinase protein was also induced by iron stress in this study and has previously been shown to serve as a relay in the ethylene-signaling pathway . Another up-regulated gene, the response regulator protein, has been shown to be involved in ethylene and other hormone signaling . The identification of so many genes (>16% of all transcripts identified) encoding ethylene response-protein gene transcripts under our experimental conditions strongly indicates the ethylene signaling pathway is involved in the soybean Fe deficiency stress response, probably serving a myriad of duties .
The increase of signaling transcripts in the severely stressed genotype likely accounts for the up-regulation of genes involved in anion transport (endomembrane protein and the putative arsA Homolog hASNA-1) and an aquaporin protein. Aquaporins are known to be induced by iron deficiency and other abiotic stresses . The induction of an aquaporin gene points toward the necessity for the cells to move nutrients and metabolites. Because of the ability of aquaporins to transport small molecules they may also be serving in the movement of additional cellular signals in stress pathways.
UDP-glucosyltransferase calalyzes the transfer of a glucosyl group from UDP-glucose to an acceptor molecule. UDP-glucosyltransferase has recently been shown to be a key enzyme in the production of isoflavones in Glycine max . Because of the role of isoflavones in the soybean stress response the induction of this gene may simply reflect a generalized response to the iron stress. However, the induction of this gene may be indicative of glucosylation of proteins for export through cellular membranes, or synthesis of oligosaccharides from cellular starch or sugars . Glucosylation of protein-linked oligosaccharides may protect them from degradation . Because the endoplasmic reticulum is the main site at which glucosylation of oligosaccharides takes place, induction of an RER1-like gene (functions in returning membrane proteins to the endoplasmic reticulum (ER)) supports this scenario.
A single aldolase gene, Fructose-bisphosphate aldolase, was under-expressed in the inefficient genotype relative to the efficient genotype. The uniqueness of this response under our experimental conditions warrants discussion. The reduced amount of this catalytic gene product may have several outcomes. Fructose-bisphosphate aldolase is an early step in the glycolysis pathway. The products of this pathway are ATP and pyruvic acid (PVA). It is unlikely that suppression of this gene during severe iron stress and chlorosis means the inefficient isoline has an adequate energy source from photosynthesis and therefore does not require the breakdown of glucose. The possible slowdown of glycolysis could result in an energetically challenged cellular environment, thus contributing further to the iron stress. The lack of evidence for increased glycolysis would also suggest that glucose levels are not depleted, leaving that molecule available for other activities (see above).
Although less supported, under-expression of the aldolase gene may result in failure to induce a critical iron homeostasis response in the inefficient genotype. The reduced amount of aldolase transcript in the inefficient genotype suggests this may not have been adequate to respond to reduced iron in the environment. Therefore, in addition to the lack of ferric reductase induction, maintenance of iron homeostasis within the cell may have been further impaired, compounding the iron stress placed on the plant. This has important implications for understanding genetic variation in soybeans' response to iron stress.
These two results (lack of induction of ferric chelate reductase and reduced aldolase transcript) probably play major roles in creating the severe stress seen in the inefficient genotype. Most of the other differentially expressed transcripts can be explained by the soybean physiological responses to the stress.
Under adverse environment conditions, such as iron stress, plants are known to produce reactive oxygen species (ROS). In this study, a peroxisomal copper containing oxidase was up-regulated in the inefficient genotype. Peroxisomal copper containing oxidase catalyzes the oxidation of amines to aldehyde, NH3 and H2O2 . ROS such as H2O2 can cause damage to proteins and lipids . The hydroponic conditions maintaining severe iron deficiency stress invoke the oxidative stress response. In a seemingly defensive reaction, the up-regulation of a peroxidase precursor points to the soybean plant responding to the increased ROS (H2O2) by increasing the amount of ROS-scavenging enzyme(s). This is not unusual. Other Strategy I plants, such as sunflower and sugar beet, also have been shown to respond to iron stress through changes in components of their antioxidative systems [49, 50].
Several of the up-regulated genes in iron stressed roots identified in this study are related to the ubiquitin/proteosome degradation pathway. These include ubiquitin, ubiquitin conjugating enzyme, and a 26S proteasoeme regulatory subunit. The up-regulation of genes in the ubiquitin/proteasome pathway plus the up-regulation of a gene for phagocytosis and a cell motility protein suggests a breakdown of cellular membranes and general deterioration of cellular health of root tissue due to iron deficiency.
Nutrient deprivation in plants has shown to induce both ubiquitin/proteasome and vacuolar degradation of proteins and lipids [43, 51]. Homologs of these genes in other species have been shown to be involved in recycling non-essential proteins and the utilization of the degraded products to maintain vital cellular function . Ubiquitin conjugating enzymes have been shown to be induced under stress conditions  including heavy metal stress. The ubiquitin response has also been associated with the regulation and downstream signaling of resistance genes . The by products of this catabolism are thought to be re-mobilized to sustain growth under stress conditions . The remobilization of the byproducts by the iron inefficient plants may provide carbon and nutrients to rapidly expanding leaves. Thimm et al. suggested a similar physiological response to iron stress, to maintain carbon flow. Garbarino et al. suggested abiotic stress results in improperly folded proteins, which are targeted for degradation by ubiquitinization. Interestingly, one of the over-expressed genes in the inefficient genotype encoded a chaperonin protein and chaperonins are needed for proper folding of nascent proteins.
The transcript for the ubiquitin conjugating enzyme was shown to be up regulated in iron inefficient plants under iron limiting conditions. In other species this enzyme has been shown to require the interaction of zinc ring finger proteins. In this study, five zinc finger protein genes were induced in the iron-stressed genotype. These zinc finger proteins may be acting as transcription factors in the regulation of the ubiquitin pathway in soybean, or they may be involved in the post translational modification of other genes known to be involved in iron homeostasis [56–58].
It is unlikely that protein and lipids are the only cellular components modified by the physiological conditions created from the iron stress. The increased expression of the 2-oxoglutarate (2OG) Fe(II)-dependent oxygenase suggest that the physiological changes brought about by iron deficiency stress has resulted in damage to DNA or modification of RNA and the soybean plant is responding to those challenges by increasing DNA repair and RNA stability. The 2OGFe(II)-dependent oxygenase has been predicted to detoxify methylated bases of ssDNA and reverse methylase modification of RNA, thus creating less toxic base derivatives, and enzymes of this family are also known to catalyze the formation of the plant hormone ethylene .
Many of the changes in transcript level observed in iron-stressed soybean correspond to general stress responses. For example, a TIR-NBS-LRR-TIR gene, a common motif in known resistance genes, was found to be up-regulated in the stressed iron-inefficient genotype. The over-expression of this gene suggests soybean responds to iron stress in a manner akin to the way it would combat pathogenic infection. Similarly, when iron-stressed, other plants such as Arabidopsis and rice, show expression changes of genes involved in wounding, abiotic and biotic stresses , as well as reproduction [61, 62]. These types of genes may all be members of common cascades involved in physiological stress responses.
It is important to note that four of the genes we identified in this experiment had no BLAST homology to the UniProt protein database (Tables 1 and 2). While this makes it difficult to determine their function, the fact that they are induced 2.9 to 3.5 fold, suggest they have very important roles and are worthy of further functional analyses.
The use of cDNA arrays has allowed us to identify transcripts differentially expressed in soybean under Fe stress conditions. Some of the genes identified are similar to general stress response genes while others may be specific to Fe stress response in soybean. It is important to note that the genes found on the cDNA array used in this study represent only a small subset of the total genic component of soybean. As such, the genes identified as differentially expressed in this study represent only a fragmented snapshot of changes occurring in the soybean physiology in response to iron deficiency stress.
However, we have been able to confirm and extend previous knowledge of soybean's iron stress responses and draw important inferences for genetic and physiological differences between soybean iron-efficient and iron-inefficient genotypes. Relative to inefficient soybean genotypes, iron-efficient genotypes may have an increased ability to respond to reduced iron availability in the environment through efficient induction of iron reductase. Root membrane bound reductase capacity is often correlated with iron efficiency. In this study, under iron limited conditions, the iron efficient plant had high root membrane reductase capacity while the iron inefficient plants reductase levels remained low, further limiting iron uptake through the root. Additionally, iron-efficient genotypes may have an efficient induction of catalytic enzymes necessary to release ATP and provide a much needed energy source to help maintain homeostasis.
Many of the genes induced in the iron inefficient NIL are involved in known stress induced pathways. The most striking response of the iron inefficient genotype to iron deficiency stress was the induction of a profusion of signaling and regulatory genes in an attempt to establish and maintain cellular homeostasis. Genes were induced that point toward an increased transport of molecules through membranes. A suppression of a key catalytic gene suggests the iron-inefficient genotype may be energetically challenged to maintain a stable cellular environment.
Many of the induced genes were obviously up-regulated in response to decreasing metabolic integrity and cellular damage. Enzymes were induced that point toward production of protein and lipid-damaging reactive oxidative species and a concomitant induction of an ROS-defensive enzyme. Genes involved in DNA repair and RNA stability were induced. Other genes were induced that are involved in protein and lipid catabolism; perhaps as an effort to maintain carbon flow and scavenge energy. These experiments have identified candidate genes and processes for further experimentation to increase our understanding of soybeans' response to iron deficiency stress.
These transcripts should serve as a starting point for future research to both understand and improve iron uptake and utilization as a step in improving overall plant health. Understanding the role these gene products play in soybean Fe metabolism could help alleviate yield loss for crops grown in calcareous soils. Further manipulation of these genes could lead to higher Fe content or increased Fe bioavailability for soybeans and other Strategy I food crops.
Near isogenic soybean lines (NILs) were developed by the USDA in 1972  specifically for their response to Fe deficiency. Fe efficient PI 548533 (Clark) was crossed with Fe inefficient PI 54619 (T203). Progeny were selfed and resulting F2 plants were screened for Fe inefficiency. The Fe inefficient progeny were backcrossed to the Fe efficient PI 548533 for six generations , resulting in an Fe inefficient plant with the Clark genetic background. The Fe inefficient isoline was released as PI 547430 (IsoClark) . Both the Fe efficient PI 548533 (Clark) and Fe inefficient PI 547430 (IsoClark) lines were grown in the Ames, Iowa USDA greenhouse under 16 hr photoperiods. Plants were germinated in sterile vermiculite with distilled deionized water. After one week they were transplanted into a DTPA nutrient buffered hydroponics system  containing all minerals necessary for normal growth. Experimental 10L systems to induce Fe deficiency stress had 50 uM Fe(NO3)3 Fe levels while systems for the control experiment contained 100 uM Fe(NO3)3. Additionally, each 10 L system contained 2 mM MgSO4*7H2O, 3 mM Mg(NO3)2*6H2O, 2.5 mM KNO3, 1 mM CaCl2*2H2O, 4.0 mM Ca(NO3)2*4H2O, 0.020 mM KH2PO4, 542.5 uM KOH, 217 uM DTPA, 1.52 uM MnCl2*4H2O, 4.6 uM ZnSO4*7H2O, 2 uM CuSO4*5H2O, 0.20 uM NaMoO4*2H2O, 1 uM CoSO4*7H2O, 1 uM NiSO4*6H2O, 10 uM H3BO3, and 20 mM HCO3. A pH of 7.8 was maintained by the aeration of a 3% CO2: air mixture. A supplemental nutrient solution containing 16 mM potassium phosphate, 0.287 mM boric acid and 355 mM ammonium nitrate was added daily to maintain proper plant nutrition. To ensure the chlorosis was due to Fe deficiency stress, A15, an Fe efficient plant, and T203, Fe inefficient plant, were included with each experimental replication. Plants were grown in the hydroponics system for two weeks, until they reached the V3 stage , at which point tissue was harvested for RNA extraction. This experiment was replicated three times, for three independent biological replicates, each with two technical replicates.
Total RNA from Fe deficient plants was extracted from root tissue of three biological replicates, each with two technical replicates, for a total of six slide hybridizations using a modified phenol:chloroform extraction with a lithium chloride precipitation . Total RNA for control samples was extracted from root tissue following the QiagenRNeasy protocol for three biological replicates each with two technical replicates for a total of six slide hybridizations. All samples were composed of root tissue from four individual plants, all grown in the same hydroponic unit. RNA purity was determined by spectroscopic readings at A260 and A280 and by formaldehyde gel visualization. Experimental samples were further purified using the RNeasy kits from Qiagen. Purified RNA was then re-analyzed to determine purity and final concentration. Each sample yielded 180 ug of purified RNA, 90 ug of purified RNA was used for each of the dye swap pairs of cDNA slides. The cDNA array for experimental samples consisted of 9,728 total cDNAs of unigene sets Gm-r1021 and Gm-r1083 spotted onto amine coated glass slides  and entered as platform GPL1013 in NCBIs Gene Expression Omnibus database [64, 65]. The cDNA array for the control samples consisted of the original 9,728 total cDNAs from unigene sets Gm-r1021 and Gm-r1083 plus an additional 9,272 total cDNAs of unigene set Gm-r1070 and entered as platform GPL3015 in GEO.
Purified RNA samples were split into 90 ug aliquots and concentrated to 10 uL in a Savant Speed Vac. The concentrated purified RNA and oligo dT was heated together for 10 min. 20 uL of 1 × Buffer, 10 mM DTT, 500 uM low T dNTPs, 100 uM Cy3 or Cy5 (Amersham Biosciences), and 13 u/uL SuperScriptII (Invitrogen) was added to each RNA/Oligo dT sample then placed at 42°C for 2 hours. Remaining RNA was degraded with an RnaseA/H treatment. The three biological replicates of the Fe deficient samples formed six technical replicates, the raw data has been deposited in GEO [64, 65] and is accessible through GEO series accession number GSE7290. The three biological replicates of the control samples formed six technical replicates, again, the raw data has been deposited in GEO [64, 65] and is accessible through GEO series accession number GSE7325. The labeled Clark (Fe-efficient) and IsoClark (Fe-inefficient) cDNA samples were mixed in a balanced dye swap design. The combined samples were purified with QIAquick PCR purification kits (Qiagen) labeled with PolyA DNA and hybridized for 18 hours at 42°C. After overnight hybridization, slides were washed (wash 1: 1 × SSC, 0.2%SDS, wash 2: 0.2 × SSC, 0.2%SDS, wash 3: 0.1 × SDS) to remove unbound cDNAs. Slides were scanned with ScanArray Express (Stratagene) and resulting images were overlaid and spots identified by the ImaGene program. An analysis program developed at the University of Illinois  was used to identify differentially expressed cDNAs. For our purposes, differential expression is defined as a minimum of two fold over or under expression in the cDNA of IsoClark (Fe-inefficient) relative to Clark (Fe-efficient).
For the RT-Real Time PCR experiments, 200 ng of RNA extracted from root tissue of plants collected over a 48-hour time course was added as initial template for each sample with Time 0 representing the time at which tissue was collected for the microarray experiment. Primers (Table 3) were designed to produce a 250 bp amplicon based on the sequences available from GenBank. Stratagene's Brilliant qRT-PCR kit was used with each 25 uL reaction assembled as described by the Stratagene instruction manual (Catalog #600532) with 2.5 uL of 50 mM MgCl2, and 2 uL of 50 nM Forward and Reverse primers as determined experimentally to optimize the reactions. Cycling protocols consisted of a 45 min. at 42°C for the reverse transcription, 10 min at 95°C to disable any remaining StrataScript, then 40 cycles of 30 sec at 95°C, 1 min at proper annealing temperature for each primer pair, 30 sec at 72°C. The PCR reactions were run in the Stratagene Mx3000P followed by a dissociation curve, taking a fluorescent reading at every degree between 55°C and 95°C to ensure only one PCR product was amplifying. The Stratagene analysis system established a threshold fluorescence level where amplicon fluorescence levels were statistically higher than background fluorescence; this threshold level is referred to as the Ct value, the cycle at which the samples fluorescence is above threshold. To be considered differentially expressed, the Fe efficient and Fe inefficient plants at the same time point had to differ in where they crossed the Ct by more than 1 cycle. One cycle difference in the RT-PCR experiments corresponds to the two-fold difference in gene transcripts between the NILs examined by the microarray experiment. The fold change was calculated from the differences in Ct using the 2ΔCt method [20, 66]. As controls, a passive reference dye was added to each sample, to ensure recorded fluorescence levels were due to SYBR green incorporation. Additionally, each sample was also run in triplicate and each sample was also normalized against tubulin amplification, see primers in Table 3, to ensure the differential expression was not due to differing amounts of initial RNA template added to each sample. As an additional negative control each sample was also run without reverse transcriptase to ensure amplification was due from RNA template.
The single linkage clustering analysis performed in this work used techniques first reported by Graham et al. . In brief, the nucleotide sequences of the 43 cDNAs identified as differentially regulated under Fe chlorosis conditions were added to a data set containing the nucleotide sequences of plant genes known to be involved in general stress responses. These general stress genes were identified based on micro/macro and bioinformatic analyses from the following published works: [52, 24–26, 67, 28, 68, 33, 30, 31, 20, 34]. Of the total 430 sequences used for clustering, 221 were derived from phosphate-starved tissues of Arabidopsis , Medicago truncatula, soybean and Phaseolus vulgaris . The remaining 209 sequences came from a variety of plant stresses [52, 24–28, 68, 33, 30, 31]. Each sequence was given a unique identifier to allow identification of the source treatment. The entire data set was then compared to itself using TBLASTX  with a minimum E-value cutoff of 10E-4. The single linkage clustering perl scripts generated by  were used to assign homologous sequences to a cluster. Note that sequences with no UniProt hit, can cluster to sequences with known annotation. Thus, clustering can be used to imply annotation.
Seeds of iron efficient and inefficient plants were germinated on germination paper for 7 days before being transplanted into the hydroponics system described above, with either 50 or 100 uM Fe(NO3)3. Plants were grown for 2 weeks in the hydroponics system. Cotyledons were removed after 7 days in hydroponics to ensure a uniform chlorotic response. Root reductase activity of the plants was measured with intact roots that were submerged for 30 min in an aerated assay solution containing 1.5 mM KNO3, 1 mM Ca(NO3)2, 3.75 mM NH4H2PO4, 0.25 mM MgSO4, 25 uM CaCl2, 25 uM H3BO3, 2 uM MnSO4, 2 uM ZnSO4, 0.5 uM CuSO4, 0.5 uM H2MoO4, 0.1 uM NiSO4, 100 uM Fe(III)-HEDTA, 100 uM BPDS (bathophenanthroline disulfonic acid), and 1 mM MES, pH 6. Iron reduction was quantified spectrophotometrically, by measuring the formation of the red-colored product, Fe(II)-BPDS3; absorbance was measured at 535 nm. An aliquot of the solution with no roots submerged in it is used as the blank. A molar co-extinction coefficient of 22.14 mM-1cm-1 was used with the measured absorbance reading to calculate the rate of reduction. There were two replicates of the experiment, each with three plants per genotype per iron concentration.
This project was supported by the North Central Soybean Research Program (grant no. 58-3625-2-459).
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.