Skip to main content
  • Research article
  • Open access
  • Published:

The genetic basis of natural variation for iron homeostasis in the maize IBM population



Iron (Fe) deficiency symptoms in maize (Zea mays subsp. mays) express as leaf chlorosis, growth retardation, as well as yield reduction and are typically observed when plants grow in calcareous soils at alkaline pH. To improve our understanding of genotypical variability in the tolerance to Fe deficiency-induced chlorosis, the objectives of this study were to (i) determine the natural genetic variation of traits related to Fe homeostasis in the maize intermated B73 × Mo17 (IBM) population, (ii) to identify quantitative trait loci (QTLs) for these traits, and (iii) to analyze expression levels of genes known to be involved in Fe homeostasis as well as of candidate genes obtained from the QTL analysis.


In hydroponically-grown maize, a total of 47 and 39 QTLs were detected for the traits recorded under limited and adequate supply of Fe, respectively.


From the QTL results, we were able to identify new putative candidate genes involved in Fe homeostasis under a deficient or adequate Fe nutritional status, like Ferredoxin class gene, putative ferredoxin PETF, metal tolerance protein MTP4, and MTP8. Furthermore, our expression analysis of candidate genes suggested the importance of trans-acting regulation for 2’-deoxymugineic acid synthase 1 (DMAS1), nicotianamine synthase (NAS3, NAS1), formate dehydrogenase 1 (FDH1), methylthioribose-1-phosphate isomerase (IDI2), aspartate/tyrosine/aromatic aminotransferase (IDI4), and methylthioribose kinase (MTK).


Iron (Fe) deficiency in maize (Zea mays subsp. mays) mostly occurs during growth on calcareous or alkaline soils, where Fe becomes sparingly soluble due to its precipitation in form of hydroxides, oxides, or phosphates [1]. Approximately 30% of the world’s arable soils are of high pH and include preferential maize cultivation areas like the river valley of Nebraska with about 0.4 million hectares (ha) [2] and the arid and semi-arid regions of the Great Plains [35]. Therefore, yield reduction in Fe-deficient maize is of agronomic importance [6].

Graminaceous plant species like maize acquire Fe by so-called strategy II mechanisms, which include the release of phytosiderophores, acting as high-affinity hexadentate chelators for ferric Fe, and an elevated expression of transport systems for Fe(III)-phytosiderophores at the root plasma membrane [7]. Comparative studies among different grass species have suggested that in particular the release of phytosiderophores is limiting most an efficient acquisition of Fe from the soil [8]. Relative to other graminaceous plant species, however, maize is generally considered as a weak strategy II plant, i.e. releasing approximately fivefold lower amounts of phytosiderophores than for instance barley [9]. This may explain at least in part the high susceptibility of maize to Fe deficiency-induced chlorosis. Although genotypical variation for chlorosis tolerance and the rate of phytosiderophore release have already been reported [10], attempts to characterize the intraspecific variation in chlorosis tolerance or other Fe efficiency traits across larger populations of maize genotypes have not yet been reported.

Previous studies on other graminaceous plant species have characterized various essential mechanisms and genes involved in Fe efficiency, i.e. the ability of plants to produce less chlorotic leaves, higher biomass, or grain yield under Fe-deficient growth conditions. In graminaceous species Fe deficiency upregulates the transcription factors IDEF1 and IRO2, which leads to an increase of phytosiderophore biosynthesis [11]. The first step in phytosiderohore biosynthesis is the conjugation of three S-adenosyl-methionine (SAM) residues by the enzyme nicotianamine synthase (NAS) to nicotianamine (NA) [12]. The linked methionine salvage pathway restores methionine levels and includes the genes methylthioadenosine/S-adenosyl homocysteine nucleosidase (MTN), methylthioribose kinase (MTK), methylthioribose-1-phosphate isomerase (IDI2), dehydratase-enolase-phosphatase (DEP), aspartate/tyrosine/aromatic aminotransferase (IDI4), and Formate dehydrogenase 1 FDH restores the Methionine for SAM synthesis [13]. NA is then subject to subsequent amino transfer by nicotianamine aminotransferase (NAAT) [14] and a reduction step by deoxymugineic acid synthase (DMAS) to yield deoxy-mugineic acid (DMA), which is the only phytosiderophore species being released by maize plants [15]. Subsequently, DMA is released by the transporter of mugineic acid 1 (TOM1) which is localized at the root plasma membrane [16]. Most of these genes being involved in phytosiderophore biosynthesis and release are subject to upregulation when the Fe nutritional status of the shoot is low [17]. Following metal chelation in the rhizosphere, the uptake of Fe(III)-phytosiderophores into root cells is mediated by membrane proteins of the yellow stripe 1/yellow stripe like (YS1/YSL) family that possess a particularly high affinity for phytosiderophore-chelated ferric Fe [18, 19]. The fairly robust upregulation of YS1 gene expression under Fe-deficiency goes along with an upregulation of ST-type sulfate transporters, most likely due to an enhanced sulphur demand for synthesis of the nicotianamine precursor SAM [20]. Inside root cells, ferric Fe may be reduced and exchange chelated to NA [21] and further transported radially for xylem loading and translocated to the shoot, where the majority of Fe is stored in ferritin (FER). In seeds and young seedlings vacuolar loading and unloading are critical for Fe efficiency too, since Fe loading of the vacuole by the vacuolar iron transporter VIT1 [22] and remobilization therefrom by natural resistance associated macrophage proteins (NRAMP3 and 4) are at least in Arabidopsis essential processes for early seedling development under Fe-limiting growth conditions [23]. While in Arabidopsis and other plants VIT1 is poorly responsive to the Fe nutritional status of the plant, NRAMP3 and NRAMP4 are upregulated under Fe deficiency [22, 23].

Despite this comprehensive knowledge on the functional aspects of Fe acquisition and homeostasis, studies examining the natural variation of Fe efficiency traits in maize and the role of natural allelic variation in determining bottle-necks of Fe efficiency in maize have remained poor. However, such information will be instrumental for the selection and development of chlorosis-tolerant maize cultivars by classical plant breeding methods. Furthermore, when such analyses are linked to molecular marker information, they have the potential to identify new genes mechanistically involved in the trait of interest, that have not been identified using classical functional genetics. This is due to the fact that in contrast to mutant screens, which consider one gene in one genetic background [24], analyses on the natural variation of traits allow discovering multiple gene actions in complex genetic backgrounds [25].

One first step to reach this goal is quantitative trait loci (QTL) mapping, which provides information on the chromosomal locations contributing to the quantitative variation of complex traits [26, 27]. Besides high resolution mapping of such QTLs, their combination with expression studies of positional candidate genes have the potential to improve our understanding of the QTL of interest.

The objectives of this study were to (i) determine the natural genetic variation of traits related to Fe homeostasis in the maize intermated B73 × Mo17 (IBM) population when these plants were grown under adequate or limiting Fe supply, (ii) identify QTLs for these traits, (iii) analyze Fe-dependent expression levels of genes known to be involved in Fe homeostasis as well as positional candidate genes from QTL analysis.



The heritability represents the genotypic contribution to the phenotypic variation. The variance analysis of the phenotypic data evaluated for each of the 13 traits using 85 intermated recombinant inbred lines (IRILs) provided the genetic and error variance values for the heritability calculation. The broad sense heritabilities for the traits evaluated under Fe deficiency ranged from 0.35 (shoot length (SL)) to 0.80 (SPAD value at leaf 3 (SP3)) (Table 1). A similar trend was observed for the broad sense heritabilities for the Fe-sufficient growth regime, which was lowest for SL (0.28) and highest for SP5 (0.80).

Table 1 Traits recorded in the current study for two Fe conditions (Fe-deficient and Fe-sufficient), where H 2 is the broad sense heritability on an entry means basis

Trait variation in the IBM population

The adjusted entry means (AEM) for the traits SP3 to SP6, root length (RL), root weight (RW), SL, shoot dry weight (SDW), SDW/SL, and lateral root formation (LAT) was on average across all IRILs lower under Fe deficiency than under Fe sufficiency (Figure 1). The trait water (H2O) showed higher AEM values for the IRILs under Fe deficiency, whereas leaf necrosis (NEC) showed on average across all IRILs no difference between both Fe regimes. The parental inbred Mo17 showed for all traits in both Fe regimes a lower AEM compared to the second parental inbred B73. The progenies showed transgressive segregation for all traits.

Figure 1
figure 1

Boxplot of the adjusted entry means of all traits for the 85 maize intermated recombinant inbred lines of the maize IBM population evaluated at Fe-deficient (white) and Fe-sufficient (grey) regimes. The adjusted entry means of the parental inbreds B73 and Mo17 are represented by a square and a triangle, respectively. The line in the box represents the median of the trait. T-test application to examine the difference for a trait between both Fe conditions. *, **, ***: P = 0.05, 0.01, and 0.001, respectively; ns, not significant.

The network analysis based on the partial correlation coefficients between all pairs of traits that focus on the correlation between the residuals revealed eight groups of traits (Figure 2). The highest positive and negative correlation was observed between SDW and SDW/SL (between 1.00 and 0.83) and between SDW/SL and SL (between -1.00 and -0.83), respectively, for both Fe regimes. Furthermore, the Pearson pairwise correlation coefficients between the SPAD value and the shoot Fe concentration increased from SP3 (0.41), SP4 (0.51), SP5 (0.61) to SP6 (0.62) under Fe deficiency. In contrast, the opposite trend was observed under Fe sufficiency where the correlation coefficients ranged from 0.08 (SP3), 0.01 (SP4), -0.28 (SP5) to -0.41 (SP6).

Figure 2
figure 2

Network representation of the partial correlations between all pairs of traits evaluated at Fe-deficient (orange) and Fe-sufficient (green) conditions for the 85 maize intermated recombinant inbred lines. Thickness and color intensity of the lines is proportional to the strength and direction of the partial correlation.

QTL analysis

The QTL analyses for the traits recorded under Fe deficiency regime using the 85 IRILs and their corresponding AEMs for each trait revealed a total of 47 QTLs (Table 2). The highest number of QTLs was detected for SP3 (8) and the lowest for RL, SL, SDW/SL, and LAT (1). The proportion of phenotypic variance explained by the QTL was highest for SP5 QTL2 (34.0%). The maximum of the proportion of phenotypic variance explained in a simultaneous fit by all QTLs for one trait was 59.4% (SP3), where the minimum was 9.6% (SL). The additive effect of the QTLs revealed that at 15 QTLs the allele increasing the trait value was contributed by Mo17.

Table 2 Summary of the quantitative trait loci (QTL) detected using the maize IBM population evaluated in a Fe-deficient nutrient solution, where Chr. is the chromosome, Pos. the position in centi Morgan on the genetic map, Add. the additive effect, %r 2 the percentage of the explained phenotypic variance, and genetic map interval in centi Morgan of the flanking markers with corresponding physical map interval including the number of genes in the corresponding QTL confidence interval according to the filtered gene set B73 RefGen_v2

Under Fe sufficiency, 39 QTLs were detected (Table 3). In dependence of the individual trait, the number of QTLs ranged from 10 (SP5) to 1 (SP4 and SL). The proportion of phenotypic variance explained by the QTL showed for QTL2 of SDW the highest (21.8%) value. The proportion of phenotypic variance explained in a simultaneous fit by all QTLs was maximal for SP5 (65.5%) and minimal for SL (9.8%). The additive effect of the QTLs indicated for 12 QTLs that the trait increasing allele was contributed by Mo17.

Table 3 Summary of the quantitative trait loci (QTL) detected using the maize IBM population evaluated in a Fe-sufficient nutrient solution, where Chr. is the chromosome, Pos. the position in centi Morgan on the genetic map, Add. the additive effect, %r 2 the percentage of the explained phenotypic variance, and genetic map interval in centi Morgan of the flanking markers with corresponding physical map interval including the number of genes in the corresponding QTL confidence interval according to the filtered gene set B73 RefGen_v2

The largest QTL confidence interval was detected for QTL5 of RW (16.5 cM) and the lowest for QTL3 of SP4 and H2O (0.3 cM) under Fe deficiency. For Fe sufficiency, the size of the confidence intervals ranged from 22.3 cM for QTL5 of SP6 to 0.1 cM for QTL3 of H2O. The number of genes detected under Fe deficiency and Fe sufficiency within these confidence intervals using the physical map information ranged from 367 (SP4 QTL4 and RW QTL7) to 0 (H2O QTL5 and QTL6) and between 273 (QTL2 of H2O) and 0 (QTL3 of SP3), respectively.

QTL confidence interval projection

The presentation of the QTL confidence intervals on the genetic map revealed for regions on chromosome 1, 4, 7, and 8 a clustering of QTLs for multiple traits (Figure 3). We observed that the genes involved in Fe homeostasis NAS3, MTN, Aconitase 1 ACO1, DEP, IDI4, FDH1, and VIT1 mapped to QTL confidence intervals.

Figure 3
figure 3

Projection of 18 genes involved in iron homeostasis on the IBM2 genetic map. The genetic position of confidence intervals of the quantitative trait loci detected for the 13 traits of our study at Fe-deficient and Fe-sufficient regimes are represented by orange and green bars, respectively. The trait value increasing alleles determined according to QTL analyses are indicated as blue and red bars for Mo17 and B73, respectively.

Expression analysis

The expression levels of genes relative to ACTIN1 ranged between 0.19 (Mo17 ST1) and 16.20 (Mo17 IDI2) under Fe deficiency and between 0.45 (B73 ST1) and 20.30 (Mo17 NAS1) under Fe sufficiency (Figure 4). No significant differences (α = 0.05) between both Fe regimes were observed for IDEF1, IRO2, mitochondrial iron transporter 1 MIT1 and NRAMP3, whereas expression of the remaining genes was significantly different (α = 0.05). The strongest expression differences between B73 and Mo17 under Fe deficiency were detected for DMAS1, IDI2, IDI4 and MTK. Furthermore, under Fe sufficiency striking differences were observed for FDH1, FER1, IDI2, IDI4, MTK, NAS1, and NAS3.

Figure 4
figure 4

Quantitative transcript levels of 17 genes important for Fe homeostasis under Fe-deficient (-Fe) and Fe-sufficient (+Fe) regimes ± standard error relative to the transcript level of ACTIN1 . Mean values of four technical replications of the parental inbred lines B73 and Mo17 marked with different letters for each gene are significantly (α = 0.05) different expressed.


Maize is often grown on soils with low Fe availability although it is highly sensitive to Fe deficiency-induced chlorosis. Therefore, understanding the genetic architecture of Fe-efficiency in maize is instrumental for the selection and development of chlorosis-tolerant maize cultivars by classical plant breeding methods. Furthermore, when such analyses are linked to molecular marker information they have the potential to identify new genes mechanistically involved in the trait of interest. Therefore, we took a two-step approach to explore the intraspecific variation in Fe-responsive traits in a segregating population of maize. First, traits related to Fe homeostasis were determined that promise to be relevant for Fe-efficiency. In a second step, these traits were used for a quantitative genetic approach and the subsequent determination of new candidate genes involved in Fe homeostasis.

Genotypic variation of traits related to Fe homeostasis

The means of all traits, except NEC, showed a highly significant (α = 0.01) difference between the two examined Fe regimes (Figure 1). This finding is in accordance with results of [28] who observed a considerable reduction in biomass and chlorophyll concentration under Fe-deficient growth conditions. Our observation illustrates the significant influence of Fe supply on the extent of phenotypical changes in the IBM population.

The broad sense heritabilities observed for the traits under consideration were moderate to high at both Fe regimes (Table 1). This is in accordance with the results of [29] who detected high heritabilities for Fe concentration in maize kernels of the IBM population. This observation indicates that the data of our study provide a reliable basis for detecting QTLs for morphological and physiological traits contributing to superior plant performance under different Fe regimes [30].

To verify whether the SPAD value could be used as indirect measure for the Fe nutritional status, we correlated the leaf chlorophyll index (SPAD value) with measured Fe concentrations. Plants grown under adequate Fe supply showed no correlation between leaf greenness measured by SPAD and the Fe concentration. This is not surprising considering the facts that i) Fe partially precipitates in the apoplast contributing to the so-called chlorosis paradox [31], and ii) the chlorophyll concentration decreases under Fe overload in the chloroplast [32]. When Fe provision to shoots is sufficient to saturate chlorophyll biosynthesis, excess Fe will be stored in ferritin [33] to prevent chlorophyll degradation.

In contrast to Fe-sufficient growth conditions, a correlation coefficient of 0.41 to 0.62 was observed for the correlation between Fe concentrations in leaves and the SPAD values obtained from Fe-deficient plants. This finding was in accordance with results of [28] who observed a tight relationship between the chlorophyll concentration and the extractable Fe in maize leaves under Fe-limiting conditions. This relationship indicates that the SPAD values measured in our study were appropriate measures of the Fe nutritional status under Fe-deficient conditions.

With respect to the SPAD measurements, the 3rd leaf of Fe-deficient plants showed higher trait values in comparison to the younger leaves 4 to 6 (Figure 1). Furthermore, the 3rd leaf exhibited the smallest differences between the Fe-deficient and Fe-sufficient regimes. Under Fe-sufficiency, no obvious difference between the means of SP3 to SP6 was observed. Only for leaf no. 3 there was a significant (α = 0.05) correlation between the SPAD measurements under both Fe regimes (Figure 2). These observations are most likely related to the fact that leaf no. 3 was formed during the early vegetative growth phase when Fe was supplied in the preculture. This also allowed plants to build up an Fe reservoir in the root apoplast which can be an important Fe source during subsequent growth [34]. The detection of loci being important for an efficient depletion of the apoplastic Fe reservoir still need to be determined. Furthermore, the loci causing the differences between both Fe regimes in the relative chlorophyll contents of SP4 to SP6 may also merit further investigations respecting as they are likely to contribute to improved Fe efficiency and thus might be valuable in breeding programs. For this reason QTL mapping is the method of choice to identify those loci that allow improving Fe efficiency in maize.

QTL and gene expression analyses

Besides the general analyses of the physiological and morphological traits evaluated under two Fe regimes for 85 IRILs we combined this phenotypic information with 1652 genetic markers for QTL mapping. Under Fe deficiency the traits SP3, RW, and H2O showed the highest number of QTLs (7 - 8) explaining between < 0.1% to 15% of the phenotypic variance (%r2) (Table 3). In contrast, SP4 to SP6 showed with 3 to 4 QTLs the highest %r2 of 34.0%. This observation suggests that under Fe deficiency SP4 to SP6 are of lower genetic complexity which increases the probability to identify in these QTL confidence intervals genes contributing largely to the natural variation in Fe homeostasis-related traits.

According to the genome sequence of the QTL intervals of SP3, SP5, and SP6 measured under Fe deficiency, these intervals include among others a Ferredoxin class gene (GRMZM2G043162) on chromosome 4 (Figure 3, Table 2). Ferredoxin transcript and protein levels strongly decrease under Fe deficiency [35] suggesting that ferredoxins respond sensitively to the Fe nutritional status. The sequence of the Ferredoxin class gene found in the above-mentioned QTL interval was not homologous to the maize ferredoxins FDX1, FDX2, FDX3, or FDX5 characterized by [36] or to FDX6. We assume that the detected Ferredoxin class gene is, like the other homologs, involved in the capture and distribution of reducing equivalents derived from photosynthetic electron transport chain in chloroplasts. However, more than the other homologs, this Ferredoxin class gene may be a candidate gene conferring differential chlorosis tolerance among maize lines. Targeted biochemical and metabolite analyses will be necessary to validate whether Ferredoxin class gene functions are crucial for chlorosis tolerance.

Under adequate Fe supply, we observed in the QTL intervals of SP3 to SP6 another putative ferredoxin gene that was located on chromosome 1 (Figure 3, Table 3). This ferredoxin gene is most likely an ortholog to PETF (GRMZM2G359127) in Chlamydomonas, which is upregulated under Fe sufficiency [37]. The precise function of this ferredoxin homolog in maize is still unclear. However, the corresponding maize protein sequence revealed a homology of 64% with the PETF gene of the cyanobacterium Fischerella (data not shown). As studies in algae have indicated that multiple ferredoxin isoforms allow for the allocation of reduction equivalents to specific metabolic pathways in the chloroplast, the putative maize ferredoxin PETF might be particularly relevant not only for electron transport in the presence or excess of Fe but also for intraspecific variation in this function.

Under Fe-deficient growth conditions, we detected in the QTL intervals for branching at the terminal 5 cm root (BTR), SP5, and SP6 the gene encoding the metal tolerance protein 8 (MTP8) (GRMZM2G116831) and in the QTL interval of NEC the MTP4 (GRMZM2G118497) gene on chromosome 1 (Figure 3, Table 2). Schaaf et al., 2004 [19] showed that the maize phytosiderophore-transporter YS1 is able to transport besides Fe also other phytosiderophore-chelated metals across the membrane, which may contribute to the typical accumulation of metals in Fe-deficient maize plants [38, 39]. Hanikenne et al., 2005 [40] and Talke et al., 2006 [41] showed that MTPs are necessary for the detoxification of excess metals by sequestering them from the cytoplasm to the vacuole. The detection of the MTPs in the QTL confidence intervals supported the notion that maize has to cope not only with Fe deficiency but also with excess accumulation of metals to prevent the cytoplasm from metal toxicity. In this regard, the results of our study suggested that these two MTPs apparently contribute to the genetic variability of Fe deficiency-induced chlorosis in maize leaves. Furthermore, the trait increasing alleles of QTL intervals including MTP8 and MTP4 were provided by B73 as well as Mo17, respectively. Therefore, the combination of trait increasing alleles from both parental inbreds at different loci may cause the transgressive segregation in progenies which can be used to breed maize inbreds with a higher tolerance against excess metal accumulation during Fe starvation (Figure 1; Fe-deficient).

The QTL confidence intervals of SP4, SP5, and SP6 monitored in Fe-sufficient plants included the nicotianamine synthase 3 (NAS3) (GRMZM2G478568) gene on chromosome 1 (Figure 3, Table 3). This observation is in accordance with results of [42] who showed that the NAS3 protein was exclusively present under Fe-sufficient growth conditions. The corresponding protein is important for NA-mediated Fe chelation under adequate or even excess Fe provision to prevent the formation of reactive oxygen species via the Fenton reaction and thus to maintain Fe homeostasis in the cytosol [21, 43]. The gene expression analyses for NAS3 showed a significant (α = 0.05) three times higher transcript level for B73 compared to Mo17 at adequate Fe growth conditions (Figure 4). Furthermore, the transcript level of NAS3 in Mo17 did not vary significantly between adequate and deficient Fe supply. This, in turn, is in accordance with the observation that the trait-increasing allele in the QTL was provided by B73.

It was proposed by [11] that expression patterns of Fe homeostasis-related genes in grasses are regulated according to the Fe status and by sensing mechanisms mediated by the transcription factors IDEF1 and IDEF2. Furthermore, these authors showed that the expression pattern of OsNAS3, which is the ortholog of ZmNAS3 [42], was downregulated during Fe limitation. Taken together, these observations suggest that the transcriptional regulation of the Mo17 allele of NAS3 might carry a disadvantageous mutation, e.g. in the upstream regulators or cis-acting elements of NAS3.

In the QTL interval for the trait NEC under Fe-deficient conditions (chromosome 1) we detected the DMAS1 (GRMZM2G060952) gene, which is essential for the last step of the phytosiderophore biosynthesis in maize [15]. Since the capacity for the synthesis and release of phytosiderophore in graminaceous plants is strongly determining chlorosis tolerance under Fe-limiting conditions [44], DMAS1 is likely to play a key role for efficient Fe acquisition irrespective of which Fe source is provided [45].

With regard to the DMAS1 gene expression, Mo17 showed under Fe deficiency a three times lower transcript level compared to B73 (Figure 4). Bashir et al., 2006 [15] showed that the DMAS genes of rice, barley, wheat, and maize were upregulated during Fe deficiency in roots allowing to enhance the production and secretion of PS. This indicates that Mo17 might have a disadvantageous allele in the promoter region of DMAS1. However, further genes which might be regulated by IDEF1 and IRO2 (cf. [46]) showed differences in expression levels between B73 and Mo17 under both Fe regimes, namely NAS1 and the methionine cycle related genes FDH1, IDI2, IDI4, and MTK (Figure 4). Besides a weaker induction by an upstream regulator, this may also be due to a disequilibirium in Fe sensing leading to a more Fe-inefficient phenotype in Mo17.

The gene coding for the transcription factor IDEF1 was not located in any QTL interval of our study and was not differentially expressed between B73 and Mo17 at either Fe regime (Figure 4). On the one hand, this might be explained by the fact that the resolution power for QTL mapping was too weak for the detection of a QTL including IDEF1 in the confidence interval. However, the present results might also be explained by the presence of another upstream regulator in the IDEF1-dependent or an independent regulatory pathway. This unknown gene could be polymorphic between B73 and Mo17 and could be hidden within the confidence intervals of the corresponding QTL detected in our study in which no other obvious candidate gene was found. The detection of regulatory transcription factors essential for chlorosis-free growth under different Fe regimes would certainly benefit from additional analyses of quantitative trait loci with expression (eQTL) and protein (pQTL) data of genes modulating Fe homeostasis.


With regard to Fe homeostasis, Mo17 contributed some advantageous alleles which caused in combination with the more advantageous genetic background of B73 a transgressive segregation in some IRILs. The morphological and physiological traits determined here indicated a moderate to high dependency on natural genetic variation suggesting a powerful basis for QTL mapping approaches. Based on our QTL mapping results, we were able to identify new putative candidate genes like Ferredoxin 1, putative ferredoxin PETF, MTP4 and MTP8 which have so far not been considered as relevant for efficient Fe homeostasis under both, low or high Fe concentrations. Furthermore, we characterized candidate gene expression and provided an insight into putative trans-acting regulation on candidate genes especially for DMAS1, NAS3, NAS1, FDH1, IDI2, IDI4, and MTK.


Plant material

The intermated recombinant inbred lines (IRILs) of the IBM population was used, which was derived from a cross of the maize parental inbreds B73 and Mo17 [47]. Due to the unavailability of seeds for the IRILs MO040, MO043, MO048, MO057, MO062, MO063, MO076, MO079, and MO344, a total of 85 IRILs were evaluated in our study.

Culture conditions and evaluated traits

Maize seeds were sterilized in a 3% NaClO solution for 3 minutes and then treated with 60°C hot water for 5 minutes. Afterwards, seeds were placed between two filter paper sheets moistened with saturated CaSO4 solution for germination in the dark at room temperature. After 6 days, the germinated seeds were transplanted to a continuously aerated nutrient solution with nutrient concentrations as described by [48]. The plants were supplied with 100 μM Fe(III)-EDTA for 7 days. From day 14 to 28, plants were supplied with 10 (Fe-deficient) or 300 (Fe-sufficient) μM Fe(III)-EDTA. The nutrient solution was exchanged every third day. Plants were cultivated from day 7 to day 28 in a growth chamber at a relative humidity of 60%, a light intensity of 170 μmol m-2 s-1 in the leaf canopy, and a day-night temperature regime of 16 h/24°C and 8 h/22°C, respectively.

Four plants of each genotype were grown in one 5 L pot. All pots were arranged in a split-plot design in the growth chamber, where the two parental genotypes were included as checks. The entire experiment was replicated b = 3 times.

The relative chlorophyll content of leaf 3, 4, 5, and 6 (SP3, SP4, SP5, and SP6) was measured with a SPAD meter (Minolta SPAD 502) at day 25 for each individual plant. Furthermore, stress symptoms like branching at the terminal 5 cm of the root (BTR) and leaf necrosis (NEC) were recorded as a visual score on a scale from 1 (high trait expression) to 9 (low trait expression). Furthermore, the lateral root formation (LAT) was recorded on a scale from 1 (low trait expression) to 9 (high trait expression) at day 26. Additionally, the root length (RL), root weight (RW), and shoot length (SL) were measured for all plants in one pot as one sample and root samples were frozen immediately in liquid nitrogen at harvest on day 28. After drying the shoot material at 70°C, shoot dry weight (SDW), water content (H2O), and the ratio between shoot dry weight and shoot length (SDW/SL) was calculated.

A total of 21 IRILs were selected such that they represented the largest possible variation of trait values for SP5. For each of these IRILs, shoot samples of four plants were pooled so that each IRIL was represented by one sample for each of the three replicates. Afterwards, the samples were ground and Fe concentrations were measured using inductively coupled plasma optical emission spectrometry (iCAP 6000 SERIES, Thermo Fisher) according to [49].

Quantitative RT-PCR analysis

Total RNA was extracted from roots of two replications of the parental inbreds B73 and Mo17 that were collected from both Fe regimes using RNeasy Plant Mini Kit (QIAGEN, Germany). Total RNA was treated with DNase (Ambion DNA-free, Invitrogen). Afterwards, cDNA synthesis was performed (SuperScript VILO, Invitrogen) and primers for candidate genes (Additional file 1: Table S1, Additional file 2: Table S2) were used for quantitative RT-PCR according to the manufacturer’s instructions (DyNAmo ColorFlash SYBR Green qPCR Kit) using ACTIN1 (NM_001155179.1) to normalize relative transcript abundances of candidate genes (Table 2).

Statistical analyses

The data of each Fe treatment were analyzed using the following mixed model:

y ik = μ + g i + r k + e ik ,

where y ik is the mean of four plants of the ith genotype in one pot of the kth replication, μ the general mean, g i the effect of the ith genotype, r k the effect of the kth replication, and e ik the residual error.

To estimate adjusted entry means (AEM) for all genotypes, g i and r k were considered as fixed. Furthermore, g i was considered as random to estimate the genotypic variance ( σ g 2 ) and the error variance ( σ e 2 ). All mixed model calculations were performed with ASReml [50].

The broad sense heritability H2 for each Fe regime was calculated as:

H 2 = σ g 2 σ g 2 + σ e 2 b .

The AEM of all genotypes for all traits and Fe regimes were tested with a Kolmogorov-Smirnov test [51] for their normal distribution. Partial correlation coefficients were assessed between all pairs of traits [52]. Network analyses of the partial correlations were prepared according to [53].

Genetic map

The publicly available genotypic data ( for the IRILs were used in our study. The genetic map positions of these markers on the IBM2 map are available ( and were the basis of our analyses. 336 markers were excluded that showed a highly significant (P < 0.001) distorted segregation (cf. [54]). The remaining 1652 markers were used for the QTL analyses. Missing genotypic information in our marker set was imputed as described by [55].

QTL analyses

Due to the high number of available markers, cofactors could not be selected using standard stepwise regression. Therefore, the following procedure was applied for each trait. One random marker was selected from each bin. Multiple stepwise regression was used to select cofactors from this set of markers based on the Bayesian information criterion (BIC) [56]. This procedure was repeated 1000 times. The average number of selected cofactors across the 1000 times repetition was used as estimator of the number of bins to study in more detail. Out of these bins 100 markers were chosen randomly and the final set of cofactors based on BIC was selected.

For the QTL mapping adjusted entry means of the 3 experimental replications for each trait, Fe regime, and each of 85 IRILs were used. The QTL analysis was carried out using the multiple QTL mapping (MQM) procedure [57] implemented in the R package ‘qtl’ version 1.21-2 [58]. The QTL detection was performed with a 2 cM (centi Morgan) step size (cf. [59]).

A total of 1000 permutation runs were performed for each trait and Fe regime to determine the α = 0.05 experiment-wise type I error for a QTL [60]. The 95% Bayesian confidence interval was calculated for each QTL location [61]. The confidence interval was expanded to the nearest flanking markers and their physical map localization was derived from B73 RefGen_v2_sequence to be able to extract all putative genes from a defined interval.

If not stated differently, all analyses were performed using the statistical software R [62].


  1. Marschner H: Mineral Nutrition of Higher Plants (Second Edition).UK: Elsevier; 1995.

    Google Scholar 

  2. Hergert GW, Nordquist PT, Petersen JL, Skates BA: Fertilizer and crop management practices for improving maize yields on high pH soils. J Plant Nutr. 1996, 19: 1223-1233. 10.1080/01904169609365193.

    Article  CAS  Google Scholar 

  3. Nordquist PT, Hergert GW, Skates BA, Petersen JL: Genetic improvement of corn for tolerance to high pH soils. J Plant Nutr. 1996, 19: 1309-1317. 10.1080/01904169609365200.

    Article  CAS  Google Scholar 

  4. Hansen NC, Hopkins BG, Ellsworth JW, Jolley VD: Iron nutrition in fieldcrops. In Iron Nutrition in Plants and RhizosphericMicroorganisms. Editedby Barton LL, Abadia J. Netherlands: Springer; 2006:23–59.

    Chapter  Google Scholar 

  5. Bashir K, Ishimaru Y, Nishizawa NK: Iron uptake and loading into rice grains. Rice. 2010, 3: 122-130. 10.1007/s12284-010-9042-y.

    Article  Google Scholar 

  6. Godsey CB, Schmidt JP, Schlegel AJ, Taylor RK, Thompson CR, Gehl RJ: Correcting iron deficiency in corn with seed row-applied iron sulfate. Agron J. 2003, 95: 160-166. 10.2134/agronj2003.0160.

    Article  CAS  Google Scholar 

  7. Kobayashi T, Nishizawa NK: Iron uptake, translocation, and regulation in higher plants. Annu Rev Plant Biol. 2012, 63: 131-152. 10.1146/annurev-arplant-042811-105522.

    Article  CAS  PubMed  Google Scholar 

  8. Römheld V, Marschner H: Mobilization of iron in the rhizosphere ofdifferent plant species. In Advances in Plant Nutrition. Edited by Tinker B,Läuchli A. New York: Praeger Scientific; 1986:155–204.

    Google Scholar 

  9. Römheld V, Marschner H: Genotypical differences among graminaceous species in release of phytosiderophores and uptake of iron phytosiderophores. Plant Soil. 1990, 123: 147-153. 10.1007/BF00011260.

    Article  Google Scholar 

  10. Bernards ML, Von Jolley D, Stevens WB, Hergert GW: Phytosiderophore release from nodal, primary, and complete root systems in maize. Plant Soil. 2002, 241: 105-113. 10.1023/A:1016084023377.

    Article  CAS  Google Scholar 

  11. Kobayashi T, Itai RN, Aung MS, Senoura T, Nakanishi H, Nishizawa NK: The rice transcription factor IDEF1 directly binds to iron and other divalent metals for sensing cellular iron status. Plant J. 2012, 69: 81-91. 10.1111/j.1365-313X.2011.04772.x.

    Article  CAS  PubMed  Google Scholar 

  12. Higuchi K, Nakanishi H, Suzuki K, Nishizawa NK, Mori S: Presence of nicotianamine synthase isozymes and their homologues in the root of graminaceous plants. Soil Sci Plant Nutr. 1999, 45: 681-691. 10.1080/00380768.1999.10415831.

    Article  CAS  Google Scholar 

  13. Nagasaka S, Takahashi M, Nakanishi-Itai R, Bashir K, Nakanishi H, Mori S, Nishizawa N: Time course analysis of gene expression over 24 hours in Fe-deficient barley roots. Plant Mol Biol. 2009, 69: 621-631. 10.1007/s11103-008-9443-0.

    Article  CAS  PubMed  Google Scholar 

  14. Inoue H, Takahashi M, Kobayashi T, Suzuki M, Nakanishi H, Mori S, Nishizawa NK: Identification and localisation of the rice nicotianamine aminotransferase gene OsNAAT1 expression suggests the site of phytosiderophore synthesis in rice. Plant Mol Biol. 2008, 66: 193-203. 10.1007/s11103-007-9262-8.

    Article  CAS  PubMed  Google Scholar 

  15. Bashir K, Inoue H, Nagasaka S, Takahashi M, Nakanishi H, Mori S, Nishizawa NK: Cloning and characterization of deoxymugineic acid synthase genes from graminaceous plants. J Biol Chem. 2006, 281: 32395-32402. 10.1074/jbc.M604133200.

    Article  CAS  PubMed  Google Scholar 

  16. Nozoye T, Nagasaka S, Kobayashi T, Takahashi M, Sato Y, Sato Y, Uozumi N, Nakanishi H, Nishizawa NK: Phytosiderophore efflux transporters are crucial for iron acquisition in graminaceous plants. J Biol Chem. 2011, 286: 5446-5454. 10.1074/jbc.M110.180026.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Giehl RF, Meda AR, von Wirén N: Moving up, down, and everywhere: signaling of micronutrients in plants. Curr Opin Plant Biol. 2009, 12: 320-327. 10.1016/j.pbi.2009.04.006.

    Article  CAS  PubMed  Google Scholar 

  18. Curie C, Panaviene Z, Loulergue C, Dellaporta SL, Briat JF, Walker EL: Maizeyellow stripe1 encodes a membrane protein directly involved in FeIII uptake. Nature. 2001, 409: 346-349. 10.1038/35053080.

    Article  CAS  PubMed  Google Scholar 

  19. Schaaf G, Ludewig U, Erenoglu BE, Mori S, Kitahara T, von Wirén N: ZmYS1 functions as a proton-coupled symporter for phytosiderophore- and nicotianamine-chelated metals. J Biol Chem. 2004, 279: 9091-9096. 10.1074/jbc.M311799200.

    Article  CAS  PubMed  Google Scholar 

  20. Astolfi S, Zuchi S, Hubberten HM, Pinton R, Hoefgen R: Supply of sulphur to S-deficient young barley seedlings restores their capability to cope with iron shortage. J Exp Bot. 2010, 61: 799-806. 10.1093/jxb/erp346.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. von Wirén N, Klair S, Bansal S, Briat JF, Khodr H, Shioiri T, Leigh RA, Hider RC: Nicotianamine chelates both Fe(III) and Fe(II) implications for metal transport in plants. Plant Physiol. 1999, 119: 1107-1114. 10.1104/pp.119.3.1107.

    Article  PubMed Central  PubMed  Google Scholar 

  22. Kim SA, Punshon T, Lanzirotti A, Li A, 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.

    Article  CAS  PubMed  Google Scholar 

  23. 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 Journal. 2005, 24: 4041-4051. 10.1038/sj.emboj.7600864.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Kobayashi Y, Koyama H: QTL analysis of Al tolerance in recombinant inbred lines of Arabidopsis thaliana. Plant Cell Physiol. 2002, 43: 1526-1533. 10.1093/pcp/pcf174.

    Article  CAS  PubMed  Google Scholar 

  25. Salt DE, Baxter I, Lahner B: Ionomics and the study of the plant ionome. Annu Rev Plant Biol. 2008, 59: 709-733. 10.1146/annurev.arplant.59.032607.092942.

    Article  CAS  PubMed  Google Scholar 

  26. Mitchell-Olds T, Pedersen D: The molecular basis of quantitative genetic variation in central and secondary metabolism in Arabidopsis. Genetics. 1998, 149: 739-747.

    PubMed Central  CAS  PubMed  Google Scholar 

  27. Zhang N, Gibon Y, Gur A, Chen C, Lepak N, Höhne M, Zhang Z, Kroon D, Tschoep H, Stitt M, Buckler E: Fine quantitative trait loci mapping of carbon and nitrogen metabolism enzyme activities and seedling biomass in the maize IBM mapping population. Plant Physiol. 2010, 154 (4): 1753-1765. 10.1104/pp.110.165787.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Thoiron S, Pascal N, Briat JF: Impact of iron deficiency and iron re-supply during the early stages of vegetative development in maize (Zea mays L.). Plant Cell Environ. 1997, 20: 1051-1060. 10.1111/j.1365-3040.1997.tb00681.x.

    Article  CAS  Google Scholar 

  29. Lung’aho MG, Mwaniki AM, Szalma SJ, Hart JJ, Rutzke MA, Kochian LV, Glahn RP, Hoekenga OA: Genetic and physiological analysis of iron biofortification in Maize Kernels. PLoS ONE. 2011, 6: e20429-10.1371/journal.pone.0020429.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Xue D, Chen M, Zhang G: Mapping of QTLs associated with cadmium tolerance and accumulation during seedling stage in rice (Oryza sativa L.). Euphytica. 2009, 165: 587-596. 10.1007/s10681-008-9785-3.

    Article  CAS  Google Scholar 

  31. Römheld V: The chlorosis paradox: Fe inactivation as a secondary event in chlorotic leaves of grapevine. J Plant Nutr. 2000, 23: 1629-1643. 10.1080/01904160009382129.

    Article  Google Scholar 

  32. Duy D, Stübe R, Wanner G, Philippar K: The chloroplast permease PIC1 regulates plant growth and development by directing homeostasis and transport of iron. Plant Physiol. 2011, 155: 1709-1722. 10.1104/pp.110.170233.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Briat JF, Lobréaux S, Grignon N, Vansuyt G: Regulation of plant ferritin synthesis: How and why. Cell Mol Life Sci. 1999, 56: 155-166. 10.1007/s000180050014.

    Article  CAS  PubMed  Google Scholar 

  34. Zhang C, Römheld V, Marschner H: Retranslocation of iron from primary leaves of bean plants grown under iron deficiency. J Plant Physiol. 1995, 146: 268-272. 10.1016/S0176-1617(11)82052-3.

    Article  CAS  Google Scholar 

  35. Tognetti VB, Zurbriggen MD, Morandi EN, Fillat MF, Valle EM, Hajirezaei MR, Carrillo N: Enhanced plant tolerance to iron starvation by functional substitution of chloroplast ferredoxin with a bacterial flavodoxin. Proc Natl Acad Sci USA. 2007, 104: 11495-11500. 10.1073/pnas.0704553104.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Hase T, Kimata Y, Yonekura K, Matsumura T, Sakakibara H: Molecular cloning and differential expression of the maize ferredoxin gene family. Plant Physiol. 1991, 96: 77-83. 10.1104/pp.96.1.77.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Terauchi AM, Lu SF, Zaffagnini M, Tappa S, Hirasawa M, Tripathy JN, Knaff DB, Farmer PJ, Lemaire SD, Hase T, Merchant SS: Pattern of expression and substrate specificity of chloroplast ferredoxins from Chlamydomonas reinhardtii. J Biol Chem. 2009, 284: 25867-25878. 10.1074/jbc.M109.023622.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. von Wirén N, Mori S, Marschner H, Römheld V: Iron inefficiency in maize mutant ys1 (Zea mays L. cv Yellow Stripe) is caused by a defect in uptake of iron phytosiderophores. Plant Physiol. 1994, 106: 71-77.

    PubMed Central  PubMed  Google Scholar 

  39. Meda AR, Scheuermann EB, Prechsl UE, Erenoglu B, Schaaf G, Hayen H, Weber G, von Wirén N: Iron acquisition by phytosiderophores contributes to cadmium tolerance. Plant Physiol. 2007, 143: 1761-1773. 10.1104/pp.106.094474.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Hanikenne M, Krämer U, Demoulin V, Baurain D: A comparative inventory of metal transporters in the green alga Chlamydomonas reinhardtii and the red alga Cyanidioschizon merolae. Plant Physiol. 2005, 137: 428-446. 10.1104/pp.104.054189.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. Talke IN, Hanikenne M, Krämer U: Zinc-dependent global transcriptional control, transcriptional deregulation, and higher gene copy number for genes in metal homeostasis of the hyperaccumulator Arabidopsis halleri. Plant Physiol. 2006, 142: 148-167. 10.1104/pp.105.076232.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Mizuno D, Higuchi K, Sakamoto T, Nakanishi H, Mori S, Nishizawa NK: Three nicotianamine synthase genes isolated from maize are differentially regulated by iron nutritional status. Plant Physiol. 2003, 132: 1989-1997. 10.1104/pp.102.019869.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. Ling HQ, Koch G, Bäumlein H, Ganal MW: Map-based cloning of chloronerva, a gene involved in iron uptake of higher plants encoding nicotianamine synthase. Proc Natl Acad Sci USA. 1999, 96: 7098-7103. 10.1073/pnas.96.12.7098.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Marschner H, Römheld V, Cakmak I: Root-induced changes of nutrient availability in the rhizosphere. Plant Nutr. 1987, 10: 1175-1184. 10.1080/01904168709363645.

    Article  CAS  Google Scholar 

  45. Yehuda Z, Shenker M, Römheld V, Marschner H, Hadar Y, Chen Y: The role of ligand exchange in the uptake of iron from microbial siderophores by gramineous plants. Plant Physiol. 1996, 112: 1273-1280.

    PubMed Central  CAS  PubMed  Google Scholar 

  46. Kobayashi T, Itai RN, Ogo Y, Kakei Y, Nakanishi H, Takahashi M, Nishizawa NK: The rice transcription factor IDEF1 is essential for the early response to iron deficiency, and induces vegetative expression of late embryogenesis abundant genes. Plant J. 2009, 60: 948-961. 10.1111/j.1365-313X.2009.04015.x.

    Article  CAS  PubMed  Google Scholar 

  47. Lee M, Sharopova N, Beavis WD, Grant D, Katt M, Blair D, Hallauer A: Expanding the genetic map of maize with the intermated B73 × Mo17 (IBM) population. Plant Mol Biol. 2002, 48: 453-461. 10.1023/A:1014893521186.

    Article  CAS  PubMed  Google Scholar 

  48. von Wirén N, Marschner H, Römheld V: Roots of iron-efficient maize also absorb phytosiderophore-chelated zinc. Plant Physiol. 1996, 111: 1119-1125.

    PubMed Central  PubMed  Google Scholar 

  49. Shi R, Weber G, Köster J, Reza-Hajirezaei M, Zou C, Zhang F, von Wirén N: Senescence-induced iron mobilization in source leaves of barley (Hordeum vulgare) plants. New Phytol. 2012, 195: 372-383. 10.1111/j.1469-8137.2012.04165.x.

    Article  CAS  PubMed  Google Scholar 

  50. Gilmour AR, Gogel BJ, Cullis BR, Thompson R: ASReml User Guide Release2.0. Hemel Hempstead: VSN International Ltd; 2006.

    Google Scholar 

  51. Chakravarti IM, Laha RG, Roy J: Handbook of Methods of Applied Statistics,.Volume I. New York, London, and Sydney: John Wiley and Sons; 1967.

    Google Scholar 

  52. Whittaker J: GraphicalModels in AppliedMultivariate Statistics. Chichester:John Wiley; 1990.

    Google Scholar 

  53. Butts CT: Network: a package for managing relational data in R. J Stat Softw. 2007, 24: 1-36.

    Google Scholar 

  54. Dufour P, Johnsson C, Antoine-Michard S, Cheng R, Murigneux A, Beckert M: Segregation distortion at marker loci: variation during microspore embryogenesis in maize. Theor Appl Genet. 2001, 102: 993-1001. 10.1007/s001220100584.

    Article  CAS  Google Scholar 

  55. Sen S, Churchill GA: A statistical framework for quantitative trait mapping. Genetics. 2001, 159: 371-387.

    PubMed Central  CAS  PubMed  Google Scholar 

  56. Schwarz G: Estimating the dimension of a model. Ann Stat. 1978, 6: 461-464. 10.1214/aos/1176344136.

    Article  Google Scholar 

  57. Jansen RC: A general Monte Carlo method for mapping multiple quantitative trait loci. Genetics. 1996, 142: 305-311.

    PubMed Central  CAS  PubMed  Google Scholar 

  58. Broman KW, Sen S: A Guide to QTL Mapping with R/qtl. Statistics for Biologyand Health. New York: Springer; 2009.

    Book  Google Scholar 

  59. Carson ML, Stuber CW, Senior ML: Quantitative trait loci conditioning resistance to phaeosphaeria leaf spot of maize caused by Phaeosphaeria maydis. Plant Dis. 2005, 89: 571-574. 10.1094/PD-89-0571.

    Article  CAS  Google Scholar 

  60. Churchill GA, Doerge RW: Empirical threshold values for quantitative trait mapping. Genetics. 1994, 138: 963-971.

    PubMed Central  CAS  PubMed  Google Scholar 

  61. Sen S: Confidence intervals for gene location - The effect of modelmisspecification and smoothing. PhD thesis, Department of Statistics,University of Chicago 1998.

    Google Scholar 

  62. R Core Team: R: A Language and Environment for Statistical Computing.Vienna: R Foundation for Statistical Computing; 2012.

    Google Scholar 

Download references


We would like to thank the Maize Genetics Cooperation Stock Center (MGCSC) for providing seeds of the IBM population. We also thank Nicole Kliche-Kamphaus, Andrea Lossow, Nele Kaul, Isabel Scheibert, and Susanne Reiner for the excellent technical support. This work was supported by research grants from the Deutsche Forschungsgemeinschaft (STI596/4-1 and WI1728/16-1) and the Max Planck Society.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Benjamin Stich.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

AB, CU, and JM carried out the hydroponic growth of maize genotypes, tissue collection, and phenotype evaluation. AB and JM performed the quantitative real time PCRs. RS collected the iron concentration data. AB analyzed the data. AB, BS, and NvW drafted the manuscript. All authors read and approved the manuscript.

Electronic supplementary material


Additional file 1: Table S1 Genes involved in Fe homeostasis, which were projected on the IBM2 genetic map of maize (Figure 4) and for quantitative RT-PCR analysis in root tissue.(CSV 2 KB)

Additional file 2: Table S2 Primer list (forward: F, reverse: R) for qRT-PCR.(CSV 1 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Benke, A., Urbany, C., Marsian, J. et al. The genetic basis of natural variation for iron homeostasis in the maize IBM population. BMC Plant Biol 14, 12 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: