# Dissection of iron signaling and iron accumulation by overexpression of subgroup Ib bHLH039 protein

Maria Augusta Naranjo-Arcos, Felix Maurer, Johannes Meiser, Stephanie

Pateyron, Claudia Fink-Straube, Petra Bauer

► To cite this version:

Maria Augusta Naranjo-Arcos, Felix Maurer, Johannes Meiser, Stephanie Pateyron, Claudia Fink-Straube, et al.. Dissection of iron signaling and iron accumulation by overexpression of subgroup Ib bHLH039 protein. Scientific Reports, 2017, 7 (1), pp.10911. 10.1038/s41598-017-11171-7. hal-05113323

HAL Id: hal-05113323

https://hal.inrae.fr/hal-05113323v1

Submitted on 24 Sep 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Dissection of iron signaling and iron accumulation by overexpression of subgroup Ib bHLH039 protein

Received: 24 January 2017

Accepted: 25 May 2017

Published online: 07 September 2017

Maria Augusta Naranjo-Arcos<fnref n="1" />,², Felix Maurer², Johannes Meiser², Stephanie Pateyron³, Claudia Fink-Straube⁴ & Petra Bauer<fnref n="1" />,²,⁵

Iron is an essential growth determinant for plants, and plants acquire this micronutrient in amounts they need in their environment. Plants can increase iron uptake in response to a regulatory transcription factor cascade. *Arabidopsis thaliana* serves as model plant to identify and characterize iron regulation genes. Here, we show that overexpression of subgroup Ib bHLH transcription factor bHLH039 (39Ox) caused constitutive iron acquisition responses, which resulted in enhanced iron contents in leaves and seeds. Transcriptome analysis demonstrated that 39Ox plants displayed simultaneously gene expression patterns characteristic of iron deficiency and iron stress signaling. Thereby, we could dissect iron deficiency response regulation. The transcription factor FIT, which is required to regulate iron uptake, was essential for the 39Ox phenotype. We provide evidence that subgroup Ib transcription factors are involved in FIT transcriptional regulation. Our findings pose interesting questions to the feedback control of iron homeostasis.

Adequate plant nutrition with iron is critical for plant development and harvest. The presence of this metal cofactor is essential for multiple biochemical electron transfer reactions in all cells. At the same time, this metal may cause deleterious effects and oxidative stress due to radical formation through the Fenton reaction. Since the bioavailability of iron varies with soil types and is mostly rather low, plants are challenged to adapt and develop an efficient iron acquisition system. Two different strategies are distinguished for iron acquisition in plants. These strategies allow plants to mobilize and take up the amounts of iron they need also under conditions of low iron supply<fnref n="1" />. Strategy I, found in plants except in grasses, is based on the uptake of ferrous iron upon soil acidification and iron reduction, aided by the exudation of small organic compounds like phenolics, riboflavin or acids². Strategy II, on the other side, present in grasses, relies on the uptake of ferric iron chelated by mugineic acid-based phytosiderophores³. Both strategies are regulated at the transcriptional level by networks of iron-regulated transcription factors²,³. Attempts are underway to characterize the major players of iron regulation by using model plant systems. The knowledge about iron uptake regulation will enable us to design novel strategies for targeted breeding of iron-efficient crops.

Three types or subfamilies of bHLH transcription factors control iron uptake in dicotyledonous plants. The bHLH proteins FER and FER-LIKE IRON DEFICIENCY-INDUCED TRANSCRIPTION FACTOR (FIT) are encoded by single genes and are essential for iron acquisition in tomato and Arabidopsis⁴⁻⁷. While FER/FIT are a specific Strategy I component, the other bHLH subgroups are present in Strategy I and II plants. The bHLH subgroup Ib comprises bHLH038, bHLH039, bHLH100 and bHLH101 in Arabidopsis and the rice homolog IRO2⁸⁻¹¹. The bHLH subgroup IVc is represented by POPEYE (PYE), bHLH104, IAA-LEUCINE RESISTANT3 (ILR3/bHLH105) and bHLH034 in Arabidopsis and the rice homolog IRO3¹²,¹³.

FIT acts in the root and up-regulates the genes encoding IRON-REGULATED TRANSPORTER1 (IRT1) and FERRIC REDUCTASE OXIDASE2 (FRO2), two major players for ferrous iron uptake in the Strategy I⁴, ⁵, ⁷, ¹⁴, ¹⁵. Coumarins aid to mobilize iron in the soil and their production and export are also under control of FIT as well as some other root cell metal homeostasis functions¹⁶⁻¹⁹. FIT is induced by iron deficiency. The FIT gene is regulated

<footnote n="1">¹Institute of Botany, Heinrich-Heine University, D-40225, Düsseldorf, Germany. ²Department of Biosciences-Plant Biology, Saarland University, D-66123, Saarbrücken, Germany. ³Transcriptomic Platform, Institute of Plant Sciences - Paris-Saclay, Plateau du Moulon, 91190, Gif-sur-Yvette, France. ⁴Leibniz Institute for New Materials gGmbH, Campus D2.2, 66123, Saarbrücken, Germany. ⁵Cluster of Excellence on Plant Sciences (CEPLAS), Heinrich-Heine University, D-40225, Düsseldorf, Germany. Correspondence and requests for materials should be addressed to P.B. (email: petra.bauer@hhu.de)</footnote>
at transcriptional level, which involves FIT itself 11. *FIT* overexpression in plants is not sufficient to promote iron uptake<sup>5, 20</sup>. FIT requires protein activation upon iron deficiency. At least in part this is mediated by FIT protein-protein interaction and stability control, that serve to fine-tune iron acquisition in response to hormone signals and abiotic stress factors<sup>21, 22</sup>. FIT was found in an active state, when *FIT* was overexpressed together with any one of the four *BHLH* subgroup Ib genes<sup>23, 24</sup>.

The bHLH subgroups Ib and IVc are part of a regulatory cascade leading to the onset of iron deficiency responses. *BHLH* genes from the subgroup Ib and *PYE* from the subgroup IVc are co-expressed in Arabidopsis across different plant growth conditions. These *BHLH* genes are up-regulated by iron deficiency along with genes encoding regulators, transporters and enzymes with iron homeostasis functions in shoots and roots, for example, BRUTUS (BTS), FERRIC REDUCTASE OXIDASE3 (FRO3), NATURAL RESISTANCE-ASSOCIATED MACROPHAGE PROTEIN4 (NRAMP4), OLIGOPEPTIDE TRANSPORTER3 (OPT3) and NICOTIANAMINE SYNTHASE4 (NAS4) 16. Subgroup IVc proteins can interact with each other and with the E3 ligase BRUTUS (BTS)<sup>12, 25-27</sup>. BTS is discussed to be an iron/metal sensor since its stability is affected by the binding of iron/metal cofactors<sup>26, 28</sup>, bHLH104, ILR3 or bHLH034 can bind directly to promoters of *BHLH* subgroup Ib genes and activate them<sup>25, 27</sup>. Loss of function of bHLH IVc factors comes with an increased sensitivity to alkalinity-mediated or general iron deficiency, while overexpression causes iron uptake<sup>25, 27</sup>. *BHLH101* overexpression was able to partially rescue the iron deficiency phenotype of *bhlh34* *bhlh104* mutants by promoting *IRT1*, *FRO2* and *FIT* expression<sup>25</sup>. Thus, these findings suggest that bHLH IVc factors act upstream of bHLH Ib and FIT transcription factors. Single knockout mutants of the four subgroup Ib *BHLH* genes did not have any impact on plant growth 11. Double and triple mutants showed different degrees of leaf chlorosis phenotypes upon iron deficiency. On the other hand, no phenotypic differences were apparent under sufficient iron supply<sup>8, 23, 29</sup>. Generally stronger leaf chlorosis and growth phenotypes were found in those double and triple mutants, in which a *bhlh039* knock-out was present<sup>23</sup>. Transcriptome analysis has shown that a subset of FIT target genes was not up-regulated in response to iron deficiency in *bhlh039* *bhlh100* *bhlh101*, but among them were not *FRO2* and *IRT1* 8. These results support that bHLH038, bHLH039, bHLH100 and bHLH101 are not fully functionally redundant.

From these observations the idea can be put forward that FIT is not only activated at protein level upon iron deficiency, but that *FIT* gene expression is induced downstream of the iron regulatory bHLH transcription factor cascade. Here, we present evidence for this idea by studying iron accumulation caused by bHLH039 overexpression. FIT was found essential for this phenotype, and transcriptomic data analysis gave hints on the hierarchy of regulatory responses upstream and downstream of this subgroup Ib bHLH protein.

## Results

**Physiological phenotypes of bHLH39 overexpression (39Ox) plants in response to iron.**
Previous reports indicated that subgroup Ib bHLH proteins are active in plant and yeast cells when overexpressed together with FIT<sup>23, 24</sup>. However, we show here that overexpression of the subgroup Ib protein bHLH039 in Arabidopsis even without overexpression of FIT can lead to significant differences in iron (Fe) acquisition responses. We have generated stable transgenic Arabidopsis lines expressing *BHLH039* behind the double cauliflower mosaic virus 35 S promoter and tagged with a sequence encoding the hemagglutinine (HA) antigen for protein immunodetection. Root growth was diminished upon sufficient iron supply (+Fe) in four independent lines, when compared to the wild type (Suppl. Figure 1A). In the wild type, *BHLH039* was induced upon iron deficiency in both roots and shoots compared to the Fe sufficiency situation (Suppl. Figure 1B,C), as expected<sup>11</sup>. In the four transgenic lines, *BHLH039* was expressed at higher level than in the respective wild type samples, showing that the transgene caused an over-expression of *BHLH039* in leaves and in roots upon + and -Fe (Suppl. Figure 1B,C). Line JM78-5 with high *BHLH039* expression was selected for detailed physiological and molecular analysis. This line was designated as 39Ox. In this line, HA-bHLH039 protein was present at +Fe and -Fe in roots and leaves of two week-old plants (Suppl. Figure 1D). Hence, we can deduce that bHLH39 protein was both present and functional at + and -Fe.

At first, we investigated whether the short root growth phenotype caused by bHLH039 overexpression was related to iron supply. 39Ox plants were grown in a six day and in a two week growth system along with wild type plants as the control, exposed to + and -Fe and examined for iron response phenotypes. As described above, six day-old 39Ox seedlings exhibited shorter roots than the wild type upon + and -Fe (Fig. 1A,B). 39Ox seedlings had dark cotyledons and red hypocotyls when exposed to +Fe, indicating stress symptoms and anthocyanin production (Fig. 1C). Red pigmentation was not found in in wild type (Fig. 1C). Anthocyanin measurements confirmed that more anthocyanins were present at + than at -Fe in 39Ox and more in 39Ox than in the wild type at +Fe (Fig. 1D).

At -Fe, Arabidopsis plants typically display an enhanced iron reductase activity in the root to mobilize external iron. Root Fe reductase activity was increased in +Fe-grown 39Ox plants compared to the wild type, as shown by qualitative and quantitative measurements (Fig. 1E,F). The two lines had similar levels of Fe reductase activity upon -Fe (Fig. 1E,F). Fe reductase activity was higher at + than at -Fe in the 39Ox plants, perhaps caused by the lack of a negative Fe regulation. Due to the high standard deviations, the differences between + and -Fe in the wild type were not found to be significant. When we analyzed the six-day seedlings at the molecular level, we found that the Fe deficiency marker genes *FRO2* and *IRT1* were expressed at higher level in 39Ox than in wild type at +Fe, conform with the iron reductase activity data (Fig. 1G). On the other hand, at -Fe, *FRO2* and *IRT1* were rather similarly expressed in 39Ox and wild type. In the wild type, *FRO2* and *IRT1* expression were induced at - versus +Fe, which was not the case in 39Ox (Fig. 1G). *FIT* was also up-regulated in 39Ox compared to wild type at +Fe (Fig. 1G). *BHLH039* expression was found about 100 times increased at +Fe and about tenfold at -Fe in 39Ox compared to the wild type, demonstrating again the overexpression effect (Fig. 1G). Interestingly, *FIT*, *FRO2* and *IRT1* were hardly expressed in 39Ox leaves in either Fe supply condition. On the other hand, *BHLH039*, was induced in wild type leaves at -Fe versus +Fe and was overexpressed in 39Ox leaves (Suppl. Figure 2A,B).
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

**Figure 1.** Iron response phenotypes of 39Ox plants. 39Ox and wild type (WT) seedling plants were grown under Fe sufficiency (+Fe) or deficiency (–Fe) in the six-day (A–F,I) or two-week agar plate assay (H left side). (A) Whole seedling growth (bar 1 cm); (B) Root lengths (n = 20); (C) Anthocyanin pigmentation; in the hypocotyl indicated by arrowheads (bar 1 mm); (D) Anthocyanin contents of seedlings (n = 3); (E) Visualization of root iron reductase activity (bar 1 cm); (F) Root Fe reductase activity (n = 4); (G) Seedling gene expression of iron deficiency genes *FRO2*, *IRT1*, *FIT* and *BHLH039* (n = 3); (H) Iron contents of leaves (left) and seeds (right) (n = 3); (I) Hydrogen peroxide contents of seedlings (n = 3). Error bars represent standard deviations. Different letters indicate significant differences between samples (p < 0.05).
FIT, FRO2 and *IRT1* were highly up-regulated at +Fe in 39Ox roots, where *BHLH039* was expressed ten to 100fold higher than in the wild type (Suppl. Figure 2C,D). *BHLH038* and *BHLH100* were not expressed in the 39Ox situation neither at + nor −Fe in leaves and roots (Suppl. Figure 2B,D). Lack of *BHLH038* and *BHLH100* expression could indicate a reduced sensing of Fe deficiency in 39Ox compared to the wild type. On the other hand, *BHLH101* showed a distinct expression compared to *BHLH038* and *BHLH100*. *BHLH101* was expressed in 39Ox similarly as in the wild type at+Fe, but the expression was reduced at −Fe (Suppl. Figure 2B,D).

Because of the generally increased iron mobilization phenotypes, we suspected that 39Ox plants might have accumulated iron. Iron measurements in leaves and seeds of 39Ox plants grown under Fe sufficiency showed a two- and fourfold increased Fe content, respectively, compared to the wild type (Fig. 1H). It was then interesting to test whether 39Ox seedlings might have experienced oxidative stress caused by the elevated Fe contents. Hydrogen peroxide measurements showed that indeed, H₂O₂ levels were higher in 39Ox exposed to +Fe compared to all other samples (Fig. 1I).

Taken together, overexpression of bHLH039 protein resulted in high Fe acquisition and high root to shoot and seed translocation of iron in seedlings. Due to this iron accumulation, 39Ox plants developed iron toxicity-based oxidative stress symptoms.

**Comparative transcriptome profiling of 39Ox plants.** As transcription factor, bHLH039 acts in a regulatory iron sensing and iron regulation cascade. We conducted a transcriptomic microarray analysis to investigate the molecular phenotypes of 39Ox in more detail and to determine how 39Ox plants responded internally to the increased bHLH039 and Fe levels. We compared gene expression between wild type and 39Ox six day-old seedlings exposed to + and −Fe. We chose the six-day growth system because we had previously collected transcriptome data sets for this condition using wild type, *fit* mutant and a triple *bhlh* knockout mutant⁸,¹⁸. The transcriptomic study was based on the CATMA two-color microarray and included three biological replicates of the samples. Similar CATMA microarray data were recently collected with *fit* mutant and wild type¹⁸. Genes with at least 1.5fold statistically significant differential regulation were retained for further analysis (Suppl. Table 1).

We found that 3745 genes were differentially regulated at +Fe comparing 39Ox with wild type, comprising 1786 down-regulated and 1959 up-regulated genes. *BHLH039*, *IRT1*, *FRO2* and *FIT* were up-regulated, which confirmed the overexpression effect and the results described in the previous paragraph (Suppl. Table 1A). 1348 genes were differentially expressed between 39Ox and wild type at −Fe, comprising 813 up-regulated and 535 down-regulated genes (Suppl. Table 1B). Several Fe-regulated genes were chosen and successfully used to validate the microarray data by quantitative RT-qPCR (Suppl. Figure 3).

To determine which physiological and cellular pathways were affected in 39Ox based on the transcriptome profiles, we conducted a gene ontology (GO) term enrichment analysis using all genes differentially regulated in 39Ox at +Fe as well as −Fe (Fig. 2, Suppl. Table 2). Genes associated with GO categories related to stress, defense and hormone responses, regulation and metal homeostasis, secondary and cell wall metabolism and hormone biosynthesis were predominantly up-regulated in 39Ox versus wild type at + and −Fe (Fig. 2, Suppl. Table 2A,C). GO categories indicative of stress, light and hormone responses, photosynthesis, rhythm and chlorophyll biosynthesis and lipid biology were predominantly down-regulated in 39Ox versus wild type at +Fe, and to a lesser extent at −Fe (Fig. 2, Suppl. Table 2B,D). At −Fe, we also observed a down-regulation of cell wall biosynthesis and glucosinolate metabolism in 39Ox compared to the wild type (Suppl. Table 2D). The overall transcriptome data therefore indicated that 39Ox suffered from stress symptoms related to metal homeostasis.

We suspected that +Fe-grown 39Ox plants should display gene expression patterns reflecting increased iron mobilization and Fe accumulation. To check this, we compared the list of 3745 differentially expressed genes between 39Ox versus wild type at +Fe with that of 1128 Fe-regulated genes¹⁸ (Suppl. Table 3). We could identify an intersection of 410 regulated genes (Fig. 3, Suppl. Table 3A). 184 out of the 410 genes were down-regulated in the wild type in response to −Fe, a condition in which *BHLH039* was highly expressed in the wild type. Among them, 81 genes were also down-regulated in 39Ox, suggestive of a repression of these genes in the presence of bHLH039 protein. These 81 genes included ribosomal genes, chloroplast function-related genes and *IRT3* (Fig. 3, Suppl. Table 3A,B). However, 103 of the 184 genes were down-regulated at - versus +Fe in the wild type (down-regulated in the presence of bHLH039 at −Fe) and up-regulated in 39Ox versus wild type at +Fe (up-regulated in the presence of bHLH039 at +Fe). These 103 genes could therefore not have been regulated by bHLH039 itself but by the level of Fe. These 103 genes comprised the three ferritin genes *FER1*, *FER3* and *FER4*, typical high iron and oxidative stress marker genes³⁰, and NAS3, *YSL1* and *YSL3*, involved in long-distance phloem-based iron translocation³¹, ³² (Fig. 3, Suppl. Table 3A). Among the 103 genes regulated in this group, a large coexpression cluster was found, related to the GO categories stress and oxidation-reduction regulation (Suppl. Table 3B). 226 out of the 410 genes were up-regulated in the wild type upon −Fe versus +Fe, when *BHLH039* was highly expressed. Among them, 95 genes were also up-regulated in the +Fe 39Ox situation, reflecting indeed an expression pattern conferred by the presence of bHLH039. These genes included −Fe-responsive genes *FIT* and *FIT* targets, e.g. *FRO2*, *IRT1*, *AHA7*, *MTPA2*, *At3g07720*, *MYB10*, *At3g12900*, *F6'H1*, *PDR9*, *COPT2*, *CYP82C3*, as well as Fe homeostasis genes *NAS1*, *IREG2* and *IRT2* 18 (Fig. 3, Suppl. Table 3A). 131 among the 226 −Fe-up-regulated genes were down-regulated in 39Ox, indicative of a response to Fe but not to bHLH039. These genes, being also part of several co-expression clusters, comprised several photo-oxidative stress and circadian rhythm regulators like *LHY1*, *CCA1*, *COL2*, *ELIP*. Very interestingly this list of genes also contained *NAS4* and *OPT3*, two components for phloem-based iron transport and transmission of a long-distance iron signal³², ³³, as well as *At1g47400* and *At3g56360*, that are part of the Fe deficiency-regulated co-expression cluster of *BHLH039* 16 (Fig. 3, Suppl. Table 3A,B). Taken together, nearly one half of the iron-responsive genes differentially expressed in 39Ox versus wild type at +Fe were iron deficiency response markers, whereas more than a half of these iron-responsive genes were suggestive of a response to high internal Fe levels. Hence, 39Ox plants displayed at the same time characteristic symptoms of iron deficiency and symptoms of iron sufficiency and even iron
<figure><img src="image_12.png" /><figcaption>**Figure 2.** Overview of the GO term enrichment analysis of 39Ox versus wild type plant transcriptome data. The heatmap was generated based on the p-values of the up (u)- and down (d)-regulation of the GO categories in four datasets, representing pair-wise comparisons of 39Ox versus wild type under Fe sufficiency (+) and Fe deficiency (−). GO categories were enriched in at least one of the four datasets. On the left, division of GO terms according to GO level 1 descriptions. Boxes and letters on the right, accurate depiction of selected biological functions, namely a, pathogen and defense b, hormone c, abiotic and oxidative stress d, metal e, rhythm, light and photosynthesis f, organic acids. Detailed results can be found in Suppl. Table 2.</figcaption></figure>

stress. In particular, the FIT target coexpression cluster was generally induced in 39Ox at +Fe, indicative of the Fe deficiency response. On the other side, the iron homeostasis co-expression cluster containing *BHLH* Ib and IVc genes, *BTS*, *NAS4* and *OPT3*, was either not regulated or down-regulated, which reflected a high iron response of
<figure><img src="image_13.png" /><figcaption>**Figure 3.** Overview of the comparison of 39Ox and Fe-regulated genes. A subset of 410 genes is contained in the intersection of 39Ox (3745 genes, this work) and Fe-regulated genes (1128 genes, 18). According to up- or down-regulation patterns, these 410 genes can be further subdivided into four groups. Represented are the numbers of regulated genes found in each group and examples of genes identified in the groups. Detailed results can be found in Suppl. Table 3.</figcaption></figure>

39Ox. These patterns of 39Ox-mediated iron-related transcriptome changes were confirmed when the expression patterns of 162 robust Fe-regulated genes were studied in 39Ox (Suppl. Table 3C,D).

Next, we used the transcriptome data for confirmation of downstream target genes that depend for their expression on bHLH039. Previously, 134 genes were identified, that are differentially expressed upon –Fe in six day-old seedlings of the triple knockout mutant line 3xbhlh (= bhlh039 bhlh100 bhlh101) versus wild type 8. This list of 134 genes (denoted groups I and III in ref. 8) contained 52 genes, that we found to be differentially regulated in 39Ox versus wild type. 25 out of the 52 genes were up-regulated in the 3xbhlh and down-regulated in 39Ox versus wild type. They were mostly related to chloroplast, rhythm and photo-oxidative stress responses in leaves (Suppl. Table 3E,F). 18 out of the 52 genes were down-regulated in 3xbhlh and up-regulated in 39Ox. These genes comprised a large coexpression cluster and four FIT and bHLH039 potential downstream target genes, namely At3g07720, At3g12900, CYP82C4 and COBL6 (Suppl. Table 3E,F). Therefore, the 39Ox transcriptome data indicated again that these latter four genes might be targets regulated by bHLH039.

Finally, we investigated the regulation of³⁴ FIT-dependent and root-expressed genes that we had previously described¹⁸, ³³ of them showed differential regulation in 39Ox versus wild type (Suppl. Table 3G). 29 genes were up-regulated in the presence of FIT and in 39Ox and included the FIT-dependent iron uptake genes *IRT1* and *FRO2*. Two genes down-regulated by FIT were also found to be down-regulated in 39Ox, namely *ZIP2* and *SCPL31*. Hence, the regulation of FIT downstream genes in 39Ox plants likely depended on the presence of FIT. Since FIT is an essential protein for iron uptake, the 39Ox effect is well explained by the activation of the FIT response pathway through bHLH039.

**Requirement of FIT for 39Ox phenotype.** If FIT is indeed an important target activated downstream of bHLH039, we predicted that FIT should be required for regulation of the FIT targets and the 39Ox phenotype at +Fe.

To test this, we investigated the 39Ox phenotype in the background of *fit* loss of function mutant plants. Such *fit* mutant plants do not mobilise sufficient iron and suffer from severe leaf chlorosis 5. We found that in contrast to 39Ox, 39Ox/*fit* plants showed similar morphological and physiological phenotypes as *fit* mutants. 39Ox/*fit* plants had leaf chlorosis and did not display the short root phenotype of 39Ox (Fig. 4A,B). Instead, 39Ox/*fit* roots grew similarly as *fit* roots. They were as long as wild type roots at +Fe, and reduced in length at –Fe (Fig. 4B). 39Ox/*fit* plants did also not have the strong anthocyanin pigmentation of 39Ox (Fig. 4C,D). However, plants with *fit* background showed more anthocyanin pigmentation than the wild type at –Fe, probably due to the iron deficiency stress. Iron reductase activity was low in 39Ox/*fit* at + and –Fe like in *fit* plants (Fig. 4C). Furthermore, gene expression analysis confirmed that *BHLH039* was overexpressed. *FRO2* and *IRT1* gene expression were low and not highly induced by –Fe nor by +Fe in 39Ox/*fit*, as it was observed in *fit* mutants (Fig. 4D). Gene expression levels of FIT-dependent genes *AT3g07720*, *At3g12900*, *MTPA2* and *CYP82C* were also comparably low in 39Ox/*fit* and *fit* mutant (Suppl. Figure 4). HA-tagged bHLH039 protein was equally expressed in 39Ox plants with wild type as well as *fit* mutant background (Fig. 4E,F). These results indicate that a functional FIT was required for the action of bHLH039 in 39Ox. FIT did not affect the protein abundance of bHLH039.
<figure><img src="image_14.png" /><figcaption>**Figure 4.** Analysis of the 39Ox effect in the *fit* mutant background. Wild type (WT), *fit* mutant (*fit*), 39Ox and 39Ox with *fit* mutant background (39Ox/*fit*) seedling plants were grown under Fe sufficiency (+Fe) or deficiency (−Fe) in the six-day agar plate assay. (A) Root growth phenotype (bar 1 cm); (B) Root lengths (n = 20); (C) Magnification of the shoots from (A) (bar 1 mm); (D) Anthocyanin contents of seedlings (n = 3); (E) Root iron reductase activity (n = 4); (F) Gene expression of *FIT*, *FRO2*, *IRT1* and *BHLH039* (n = 3); (G) Immunoblots with anti-HA antibodies and actin as a loading control. Arrowheads indicate the positions of the bands corresponding to HA₃-bHLH039 (36 kD) and actin (42 kD). *fit* seedlings grown at +Fe were used as a negative HA-protein detection control. (H) Protein abundance based on immunoblot signals of HA₃-bHLH039 normalized to actin (n = 3). Error bars represent standard deviations. Different letters indicate significant differences between samples (p < 0.05).</figcaption></figure>

Above, it was shown that Fe reductase activity was found along the roots of 39Ox seedlings. Thus, we suspected that *FIT* expression might also have been activated all along the root in response to 39Ox. To investigate in which region of the root the *FIT* gene was activated, we localized *FIT* promoter activity based on *FIT* promoter-driven beta-glucuronidase reporter activity in transgenic pFIT::GUS plants with 39Ox background. In wild type background, *FIT* promoter activity was detected mainly in the epidermis cells of the elongation root zone in Fe-deficient conditions, as expected 5 (Fig. 5). GUS activity was stronger in plants with 39Ox than wild-type background at +Fe (Fig. 5B). GUS activity was detected along the entire 39Ox roots, but not in
<figure><img src="image_15.png" /><figcaption>**Figure 5.** Analysis of *FIT* promoter-driven b-glucuronidase (GUS) activity. pFIT::GUS seedlings with wild type (pFIT::GUS/WT) and 39Ox mutant background (pFIT::GUS/39Ox) were grown under Fe sufficiency (+Fe) or deficiency (–Fe) in the six-day agar plate assay. (A) Qualitative GUS activity assay, on the left, whole plant (bar 1 cm) with squares, representing close-ups on the right (bar 1 mm). The three squares represent the transition from hypocotyl to upper root zone (top), middle zone of the root (middle), and root tips (bottom). (B) Quantitative GUS activity assay (n = 3). Error bars represent standard deviations. Different letters indicate significant differences between samples (p < 0.05).</figcaption></figure>

hypocotyls and cotyledons (Fig. 5A). Thus, overexpression of bHLH039 activated the *FIT* promoter in an ectopic manner in roots.

Discussion

Here we showed, by using overexpression plants and combinations of transgenic and mutant plants, that bHLH039 induces directly or indirectly *FIT* in the Fe deficiency regulatory cascade and that bHLH039 is active even in the presence of sufficient iron.

A remarkable conclusion from this study is that bHLH039 protein was not only abundant in 39Ox plants but also functionally active regardless of Fe supply. On the contrary, FIT and HA-FIT were present at + and -Fe in plant cells in an overexpression situation, but only active at -and not +Fe<sup>5, 20</sup>. It was known that *BHLH039* is manifold induced at -Fe, while at +Fe hardly any transcripts are produced. In contrast, the *FIT* gene is expressed at +Fe and induced about threefold at -Fe 5. Hence, the action of bHLH039 is regulated in the plant mainly via gene expression regulation, while the action of FIT is regulated at transcriptional and protein level. The transcriptional activation of *BHLH039* makes sense in view of the current models about the iron deficiency regulatory cascade. According to these models gene expression of *BHLH039* and of other *BHLH* subgroup Ib genes is induced by subgroup IVc bHLH factors in response to -Fe<sup>25, 27</sup>.

The most striking 39Ox phenotypes were apparent at +Fe but not at -Fe, showing that the utilization of Fe was affected in the 39Ox plants. First, the 39Ox plants displayed high expression of iron acquisition genes, high Fe reductase activity, high Fe contents in leaves and seeds and high ferritin gene expression, in contrast to the wild type. These phenotypes clearly show that 39Ox plants took up more Fe from the growth medium than the wild type. Second, 39Ox plants displayed oxidative stress symptoms at +Fe, like up-regulation of oxidative stress gene expression, anthocyanin pigmentation, short root phenotypes and H₂O₂ production. Both Fe accumulation and oxidative stress phenotypes of 39Ox were Fe-dependent. We interpret the observed 39Ox phenotypes as follows (depicted in a model in Suppl. Figure 5): The transcription factor bHLH039 up-regulated together with low levels of FIT protein the *FIT* gene. FIT and bHLH039 then induced iron acquisition genes like *FRO2* and *IRT1*. The resulting FRO2 activity and IRT1 stimulated iron uptake into the roots. Iron was translocated to shoots and high amounts of Fe accumulated in leaves and seeds. The high Fe levels in roots and shoots then caused toxicity symptoms such as oxidative stress, short root growth and anthocyanin pigmentation. At least half of the Fe-regulated differentially expressed genes in 39Ox versus wild type were indicative of a high iron response that normally serves to alleviate iron stress. For example, ferritin *FER* gene expression was increased, together with *NAS3*, *YSL1* and *YSL3*. These genes were up-regulated in response to the high Fe status in the shoot, indicating that iron was bound and/or circulated in the plant. Ferritin was presumably used to suppress oxidative stress and to store Fe in plastids<sup>30-32</sup>. Even though these protective responses were enhanced, we need to assume that they were not sufficiently effective to alleviate the iron stress in 39Ox. Surprisingly, repressive iron signals, that normally would be expected to delimit excessive iron uptake, were switched off or ineffective in +Fe-grown 39Ox plants. It has been hypothesized that BTS could be a sensor of Fe and negative regulator in the transcription factor cascade<sup>12, 26, 28</sup>. BTS protein function and regulation in plant cells in response to cellular Fe awaits further experimentation 26. BTS was not up-regulated in 39Ox and therefore BTS could not likely have repressed the cascade. Perhaps low BRUTUS (BTS) E3 ligase protein levels or levels of related proteins were active in 39Ox and effectively suppressed
the Fe deficiency cascade at the level of subgroup IVc bHLH proteins. *BHLH* IVc genes are not Fe-regulated, with the exception of *PYE* that we did not detect in our list of 39Ox up-regulated genes either. Moreover, very interestingly, *BHLH038*, *BHLH100* and *BHLH101* were down-regulated in 39Ox at + and -Fe. Since subgroup IVc bHLH proteins induce the promoters of *BHLH* Ib genes, we suspect that the Fe deficiency cascade upstream of *BHLH* subgroup Ib genes was switched off<sup>25, 27</sup>. Furthermore, four genes linked with that cascade and normally induced by -Fe and co-expressed with *BHLH039*, namely *OPT3*, *NAS4*, *At1g47400* and *At3g56360*, were even down-regulated in 39Ox. Overexpression of bHLH039 must have counteracted any suppressor effects by BTS or related proteins. Our results are therefore in agreement with previous published data, such as overexpression of a bHLH subgroup IVc protein leading to iron over-accumulation phenotypes similar as seen in 39Ox or *bts* loss of function plants<sup>25, 27</sup>. A long distance iron suppressor signal transported via the phloem from shoot to root has been hypothesized to communicate the shoot iron status to the root<sup>11, 32, 34, 35</sup>. A long-distance iron signal was suspected to be loaded into the phloem of leaf cells via OPT3<sup>33, 34</sup>. Such a long-distance signal was either inactive in 39Ox or in case it acted at the level of the transcription factor cascade upstream of bHLH039 not effective in 39Ox. The down-regulation of *OPT3* in 39Ox could reflect a potential mechanism for the inhibition of transmission of a long distance iron suppressor signal from leaf to root. However, it could also be that OPT3 may not be involved in the transmission of a signal. Finally, local suppressor mechanisms for iron uptake are described that suppress the activity of individual proteins mediating iron uptake in the root epidermis cell. For example, IRT1 localization is altered in the presence of metals<sup>36</sup>. Under sufficient Fe, IRT1 is more prone to degradation<sup>37</sup>. FRO2 is described to be inactive upon sufficient Fe despite of the presence of this protein in a *FRO2* overexpression situation<sup>38</sup>. It is surprising that these negative regulation mechanisms acting inside root epidermal cells were also not in place in the 39Ox background. This suggests that protein regulation factors for IRT1 and FRO2 might be related to the action of FIT.

bHLH039 action required functional FIT. 39Ox/fit plants lost the Fe-dependent 39Ox phenotypes. Therefore, FIT must have been functionally active in 39Ox plants, and again, potential inhibitory components previously reported to negatively control FIT protein activity²⁰ were not present or in place in 39Ox. This gives support to the idea that the main activator of FIT protein in 39Ox was bHLH039. The same occurs in wild type plants exposed to -Fe (see model in Suppl. Figure 5). Once subgroup Ib bHLH factors are produced in root cells in response to the regulatory Fe deficiency cascade, FIT becomes activated at protein level through direct protein interactions, as shown²³, ²⁴. FIT was found to be involved in its own feed-forward gene expression regulation¹¹. In 39Ox plants, the FIT gene promoter was activated in an ectopic manner in roots but not in leaves. Hence through the interaction with bHLH039 in 39Ox plants, FIT might have further stimulated its own transcription all along the root. This scenario can also explain why molecular Fe acquisition responses like *IRT1* and *FRO2* induction were only switched on in roots but not in leaves of 39Ox. FIT is a root-specific transcription factor 5 and hence a basal level of FIT was required at +Fe for bHLH039 action in 39Ox. Different scenarios of direct or indirect positive bHLH039-mediated FIT promoter induction are conceivable. One possibility is that bHLH039-FIT heterodimers bind to the FIT promoter (direct activation of FIT promoter by FIT-bHLH039 heterodimers, evidence for FIT heterodimerization in²⁴). Alternatively, FIT may bind to its promoter as a homodimer being activated by binding in addition to bHLH039 (direct activation of FIT promoter by multiprotein complexes consisting of FIT homodimers and bHLH039, evidence for FIT homodimerization in²⁴). Finally, an additional yet unknown transcription factor could be induced by bHLH039 and FIT and then serve to activate the FIT promoter and FIT target genes (indirect regulation of FIT promoter by FIT and bHLH039). High Fe stress-related anthocyanin production was not found in 39Ox leaves with *fit* mutant background. Hence, such shoot responses could not have been under direct or indirect control of bHLH039. Therefore, anthocyanin production was more likely the result of secondary effects related to iron accumulation in the shoot rather than the action of bHLH039 in the cells.

In conclusion, bHLH039 is a direct or indirect activator of FIT transcription in the root. Any suppressive Fe signals were inactive in 39Ox and did not inhibit Fe uptake. By manipulation of bHLH039 expression, it was thus possible to enhance iron levels in shoots and seeds. Our findings raise interesting questions with regard to the coupling of iron signals with post-translational negative control mechanisms acting upon iron uptake components. An interesting question to pursue in the future is to tackle the roles of bHLH subgroup Ib factors in the shoot. Since the Fe levels of 39Ox leaves and seeds were enhanced, it is conceivable that bHLH039 was involved in regulating the delivery of Fe to sinks.

## Methods

**Plant material and growth conditions.** The wild type ecotype used was Col-0. *fit-3* (termed *fit* throughout the text) and pFIT::GUS plants were previously described⁵.

bHLH039 overexpression plants were constructed as follows: A genomic DNA fragment containing exons and introns of *BHLH039* was inserted into pALLIGATOR yielding a double cauliflower mosaic virus 35S promoter 2xpaCMV::triple hemagglutinine (HA) tag-gBHLH039 coding cassette, according to a previously described procedure²⁰. Transgenic lines were generated, genotyped and multiplied until homozygous lines were obtained. Four lines were investigated and one line selected for further analysis (here designated 39Ox), as described in the text (see Suppl. Figure 1). 39Ox/fit-3 (termed 39Ox/fit) and pFIT::GUS/39Ox plants were generated after crossing and selection of homozygous lines.

For experiments, seeds were surface sterilized and plants raised in a six-day agar plate assay in the presence of 50 microM Fe (+Fe) or 0 Fe (−Fe) or in a two-week agar plate assay with a three-day exposure to medium containing 50 microM Fe (+Fe) or 50 microM ferrozine and 0 Fe (−Fe) using a modified Hoagland medium, as previously described 5,22 and as indicated in the text and figure legends.
**Morphological and physiological plant analysis.** Seedlings grown on agar plates were pictured with the Axio Zoom.V16 Stereo-Zoom-Microscope (Zeiss, Jena, Germany). Root lengths were measured after photography with the help of the ImageJ software tool.

Anthocyanin extracts were prepared from 15 preweighted seedlings by incubation for 5 minutes at 95° Celsius and then overnight at 25° Celsius with a propanol:HCl:H₂O (18:1:81 v/v/v) mixture, as described³⁹. The extract was photometrically measured at OD535. After background subtraction (OD<sub>650</sub>) the values were normalized to the respective sample weights.

Iron reductase activity was determined using a quantitative liquid and qualitative agar plate assay with ferro-zine as indicator for Fe²⁺, as previously described⁴⁰. Quantitative iron reductase activity was normalized to root length.

For iron content determination, leaves were harvested, dried, powdered and processed for graphite furnace atomic absorption spectroscopy, as previously described²². Seeds were dried, weighed and microwave-digested (Multiwave 3000, Anton Paar GmbH, Graz, Austria) with concentrated HNO₃. Fe concentrations were determined by inductively-coupled plasma-optical emission spectrometry ICP-OES with a Czerny-Turner type monochromator (Ultima 2, Horiba Jobin-Yvon, New Jersey, USA, and Longjumeau, France).

Hydrogen peroxide contents were measured using the Amplex red hydrogen peroxide/peroxidase assay kit (Molecular Probes) as previously described²¹.

**beta-glucuronidase (GUS) assay.** For localization of GUS activity, pFIT::GUS and pFIT::GUS/39Ox plants were incubated in the GUS staining buffer (50 mM sodium phosphate, 2 mM potassium ferrocyanide, 2 mM potassium ferricyanide, 0.2% Triton X-100, and 2 mM GUS substrate 5-bromo-4-chloro-3-indolyl-b-D-glucuronic acid) according to⁴¹. Leaf chlorophyll was removed by incubation in 70% ethanol.

Quantitative GUS activity analysis using total protein extract was conducted as previously described, following a fluorescence measurement method with MUG substrate (4-methylumbeliferyl-beta-D-glucuronide)⁵.

**Gene expression analysis by RT-qPCR.** A detailed protocol of the RT-qPCR procedure is available⁴⁰, ⁴². RNA from about 60 seedlings per biological replicate was extracted. cDNA was synthetized with poly(dT) primer from DNase-treated total RNA (300–500 ng) and then analyzed by real time qPCR using the DyNAmo ColorFlash SYBR Green qPCR Kit (Thermo Scientific, Waltham, USA) and the ICycler (Bio-Rad, Hercules, USA), as described⁴². Primers and primer sequences are listed in⁸,¹⁸. qPCR data accuracy and reliability was checked as described⁴⁰,⁴². The initial transcript amount of two technical replicates was determined by standard curve analysis using mass standards and averaged. The absolute expression data were normalized with the constitutive expression control gene encoding elongation factor 1Bα1 (EF1Bα1). The mean of three biological replicates per sample was calculated and used for statistical treatment of RT-qPCR data.

**Microarray-based transcriptome analysis.** Six-day wild type and 39OX seedlings were raised on +Fe an-Fe in three biological replicates. 60 seedlings were pooled per sample and RNA extracted using the Qiagen plant RNA extraction kit (Qiagen, Hilden, Germany). A two-color CATMAv6 microarray analysis was performed including a dye-swap control, following the published procedure⁴³. After the statistical analysis gene lists were retained with log₂ normalized intensities per samples, log₂ normalized ratios, and p-values (FDR correction). Probes with a Bonferroni p-value of 0.05 or less and a fold change of 1.5 or more were considered differentially expressed. The same RNA preparations were also used for validation by RT-qPCR of selected genes.

Transcriptome data sets were compared using Excel. Gene ontology (GO) term and coexpression analysis was conducted and visualized using VirtualPlant⁴⁴, Revigo⁴⁵ and STRING⁴⁶, as indicated in the text and legends.

**Immunoblot analysis.** Immunoblot analysis was conducted according to²¹. Around 40 six day-old seedlings were harvested, frozen in liquid nitrogen and grinded using the Precellys (PeqLab Biotechnologie GmbH, Erlangen, Germany). Per 1 mg plant powder 1 microl of 2x Laemmli buffer was added. After centrifugation 10 minigro heat-denatured total protein was loaded on a 12% SDS-polyacrylamide gel and proteins separated by electrophoresis. Proteins were electro-transferred to a nitrocellulose membrane (GE Health Care, Amersham) with transfer buffer. HA₃-bHLH039 protein was detected after immunoblotting with anti-HA-horseradish peroxidase (HRP)-coupled high-affinity monoclonal rat antibody (3F10 clone 1:1.000 dilution; Roche Diagnostics GmbH, Basel, Switzerland) by chemiluminescence (ECL chemiluminescence kit, GE Health Care, Little Chalfont, UK) and the FluoChem Q (Biozym, Hamburg, Germany) device. The membrane was stripped and reused for the detection of actin protein with anti-actin antibodies (rabbit anti-ACT, 1:5000 dilution; Agrisera, Vännäs, Sweden) and goat anti-rabbit-HRP-coupled secondary antibody (Pierce, 1:1000 dilution; Thermo Fisher Scientific, Carlsbad, USA).

The intensities of the bands were determined using the software AlphaView (FluorChem Q). The local backgrounds were subtracted. HA signal intensity was normalized to that for actin. Alternatively and as mentioned in the text and figure legend, normalization was conducted to major bands observed and measured after Ponceau S staining of the membrane for total protein.

**Statistical analysis.** Statistical analysis was conducted using the one-way ANOVA test followed by the post-hoc Tukey test for pair-wise comparisons. p-values of 0.05 or lower indicate statistical significance between two mean sample values, represented in the figures by different letters. p-values above 0.05 indicate no statistically significant difference between the two compared sample values, indicated by the same letter above the two sample values in the figures. The diagrams represent the mean values of biological replicates and the error bars reflect standard deviations.
References

1. Marschner, H. & Römheld, V. Strategies of plants for acquisition of iron. *Plant Soil* **165**, 261–274 (1994).

2. Brumbarova, T., Bauer, P. & Ivanov, R. Molecular mechanisms governing Arabidopsis iron uptake. *Trends Plant Sci.* **20**, 124–133, doi:10.1016/j.tplants.2014.11.004 (2015).

3. Kobayashi, T., Nakanishi Itai, R. & Nishizawa, N. K. Iron deficiency responses in rice roots. *Rice* **7**, 27, doi:10.1186/s12284-014-0027-0 (2014).

4. Colangelo, E. & Guerinot, M. The essential basic helix-loop-helix protein FIT1 is required for the iron deficiency response. *Plant Cell* **16**, 3400–3412, doi:10.1105/tpc.104.024315 (2004).

5. Jakoby, M., Wang, H., Reidt, W., Weisshaar, B. & Bauer, P. FRU (bhlh029) is required for induction of iron mobilization genes in *Arabidopsis thaliana*. *FEBS Lett.* **577**, 528–534, doi:10.1016/j.feb.2004.10.062 (2004).

6. Ling, H. Q., Bauer, P., Bereczky, Z., Keller, B. & Ganal, M. The tomato *fer* gene encoding a bhlh protein controls iron uptake responses in roots. *Proc. Natl. Acad. Sci. USA* **99**, 13938–13943, doi:10.1073/pnas.212448699 (2002).

7. Yuan, Y. X., Zhang, J., Wang, D. W. & Ling, H. Q. Atbhlh29 of *Arabidopsis thaliana* is a functional ortholog of tomato FER involved in controlling iron acquisition in strategy in plants. *Cell Res.* **15**, 613–21, doi:10.1038/sj.cr.7290331 (2005).

8. Maurer, F., Naranjo Arcos, M. & Bauer, P. Responses of a triple mutant defective in three iron deficiency-induced basic helix-loop-helix genes of the subgroup Ib(2) to iron deficiency and salicylic acid. *PLoS One* **9**, e99234, doi:10.1371/journal.pone.0099234 (2014).

9. Ogo, Y. *et al.* OsIRO2 is responsible for iron utilization in rice and improves growth and yield in calcareous soil. *Plant Mol. Biol.* **75**, 593–605, doi:10.1007/s11103-011-9752-6 (2011).

10. Ogo, Y. *et al.* Isolation and characterization of IRO2, a novel iron-regulated bhlh transcription factor in graminaceous plants. *J. Exp. Bot.* **57**, 2867–2878, doi:10.1093/jxb/erl054 (2006).

11. Wang, H. Y. *et al.* Iron deficiency-mediated stress regulation of four subgroup Ib *BHLH* genes in *Arabidopsis thaliana*. *Planta* **226**, 897–908, doi:10.1007/s00425-007-0535-x (2007).

12. Long, T. *et al.* The bhlh transcription factor POPEYE regulates response to iron deficiency in Arabidopsis roots. *Plant Cell* **22**, 2219–2236, doi:10.1105/tpc.110.074096 (2010).

13. Zheng, L. *et al.* Identification of a novel iron regulated basic helix-loop-helix protein involved in Fe homeostasis in *Oryza sativa*. *BMC Plant Biol.* **10**, 166, doi:10.1186/1471-2229-10-166 (2010).

14. Robinson, N. J., Procter, C. M., Connolly, E. L. & Guerinot, M. L. A ferric-chelate reductase for iron uptake from soils. *Nature* **397**, 694–697 (1999).

15. Vert, G. *et al.* IRT1, an Arabidopsis transporter essential for iron uptake from the soil and for plant growth. *Plant Cell* **14**, 1223–1233 (2002).

16. Ivanov, R., Brumbarova, T. & Bauer, P. Fitting into the harsh reality: Regulation of iron deficiency responses in dicotyledonous plants. *Mol. Plant* **5**, 27–42, doi:10.1093/mp/ssr065 (2012).

17. Mai, H. J. *et al.* Iron and FER-like iron deficiency-induced transcription factor (FIT)-dependent regulation of proteins and genes in *Arabidopsis thaliana* roots. *Proteomics* **15**, 3030–3047, doi:10.1002/pmic.201400351 (2015).

18. Mai, H. J., Pateyron, S. & Bauer, P. Iron homeostasis in *Arabidopsis thaliana*: transcriptomic analyses reveal novel FIT-regulated genes, iron deficiency marker genes and functional gene networks. *BMC Plant Biol.* **16**, 211, doi:10.1186/s12870-016-0899-9 (2016).

19. Schmid, N. B. *et al.* Feruloyl-coA 6'-hydroxylase1-dependent coumarins mediate iron acquisition from alkaline substrates in Arabidopsis. *Plant Physiol.* **164**, 160–72, doi:10.1104/pp.113.228544 (2014).

20. Meiser, J., Lingam, S. & Bauer, P. Post-translational regulation of the Fe deficiency bhlh transcription factor FIT is affected by iron and nitric oxide. *Plant Physiol.* **157**, 2154–2166, doi:10.1104/pp.111.183285 (2012).

21. Le, C. T. *et al.* ZINC FINGER OF ARABIDOPSIS THALIANA12 (ZAT12) interacts with FER-LIKE IRON DEFICIENCY-INDUCED TRANSCRIPTION FACTOR (FIT) linking iron deficiency and oxidative stress responses. *Plant Physiol.* **170**, 540–57, doi:10.1104/pp.15.01589 (2016).

22. Lingam, S. *et al.* Interaction between the bhlh transcription factor FIT and ethylene INSENSITIVE3/ETHYLENE INSENSITIVE3-LIKE1 reveals molecular linkage between the regulation of iron acquisition and ethylene signaling in Arabidopsis. *Plant Cell* **23**, 1815–1929, doi:10.1105/tpc.111.084715 (2011).

23. Wang, N. *et al.* Requirement and functional redundancy of Ib subgroup bhlh proteins for iron deficiency responses and uptake in *Arabidopsis thaliana*. *Mol. Plant* **6**, 503–13, doi:10.1093/mp/sss089 (2013).

24. Yuan, Y. *et al.* FIT interacts with ATBHLH38 and ATBHLH39 in regulating iron uptake gene expression for iron homeostasis in Arabidopsis. *Cell Res.* **18**, 385–397, doi:10.1038/cr.2008.26 (2008).

25. Li, X., Zhang, H., Ai, Q., Liang, G. & Yu, D. Two bhlh transcription factors, bhlh34 and bhlh104, regulate iron homeostasis in *Arabidopsis thaliana*. *Plant Physiol.* **170**, 2478–93, doi:10.1104/pp.15.01827 (2016).

26. Selote, D., Samira, R., Matthiadis, A., Gillikin, J. W. & Long, T. A. Iron-binding E3 ligase mediates iron response in plants by targeting basic helix-loop-helix transcription factors. *Plant Physiol.* **167**, 273–86, doi:10.1104/pp.114.250837 (2015).

27. Zhang, J. *et al.* The bhlh transcription factor BHLH104 interacts with IAA-LEUCINE RESISTANT3 and modulates iron homeostasis in Arabidopsis. *Plant Cell* **27**, 787–805, doi:10.1105/tpc.114.132704 (2015).

28. Kobayashi, T. *et al.* Iron-binding haemerythrin ring ubiquitin ligases regulate plant iron responses and accumulation. *Nat. Commun.* **4**, 2792, doi:10.1038/ncomms3792 (2013).

29. Sivitz, A. B., Hermand, V., Curie, C. & Vert, G. Arabidopsis bhlh100 and bhlh101 control iron homeostasis via a FIT-independent pathway. *PLoS One* **7**, e44843, doi:10.1371/journal.pone.0044843 (2012).

30. Reyt, G., Boudouf, S., Boucherez, J., Gaymard, F. & Briat, J. F. Iron- and ferritin-dependent reactive oxygen species distribution: impact on Arabidopsis root system architecture. *Mol. Plant* **8**, 439–53, doi:10.1016/j.molp.2014.11.014 (2015).

31. Chu, H. *et al.* Successful reproduction requires the function of Arabidopsis YELLOW STRIPE-LIKE1 and YELLOW STRIPE-LIKE3 metal-nicotianamine transporters in both vegetative and reproductive structures. *Plant Physiol.* **154**, 197–210, doi:10.1104/pp.110 (2010).

32. Schuler, M., Rellán-Álvarez, R., Fink-Straube, C., Abadía, J. & Bauer, P. Nicotianamine functions in the phloem-based transport of iron to sink organs, in pollen development and pollen tube growth in Arabidopsis. *Plant Cell* **24**, 2380–23400, doi:10.1105/tpc.112.099077 (2012).

33. Zhai, Z. *et al.* OPT3 is a phloem-specific iron transporter that is essential for systemic iron signaling and redistribution of iron and cadmium in Arabidopsis. *Plant Cell* **26**, 2249–2264, doi:10.1105/tpc.114.123737 (2014).

34. Garcia, M. J. *et al.* Shoot to root communication is necessary to control the expression of iron-acquisition genes in strategy I plants. *Planta* **237**, 65–75, doi:10.1007/s00425-012-1757-0 (2013).

35. Vert, G., Briat, J. & Curie, C. Dual regulation of the Arabidopsis high-affinity root iron uptake system by local and long-distance signals. *Plant Physiol.* **132**, 796–804, doi:10.1104/pp.102.016089 (2003).

36. Barberon, M. *et al.* Polarization of IRON-REGULATED TRANSPORTER 1 (IRT1) to the plant-soil interface plays crucial role in metal homeostasis. *Proc. Natl. Acad. Sci. USA* **111**, 8293–8298, doi:10.1073/pnas.1402262111 (2014).

37. Shin, L. *et al.* IRT1 DEGRADATION FACTOR1, a ring E3 ubiquitin ligase, regulates the degradation of IRON-REGULATED TRANSPORTER1 in Arabidopsis. *Plant Cell* **25**, 3039–3051, doi:10.1105/tpc.113.115212 (2013).
38. Connolly, E., Campbell, N., Grotz, N., Prichard, C. & Guerinot, M. Overexpression of the FRO2 ferric chelate reductase confers tolerance to growth on low iron and uncovers posttranscriptional control. *Plant Physiol.* **133**, 1102–1110, doi:10.1104/pp.103.025122 (2003).

39. Deikman, J. & Hammer, P. E. Induction of anthocyanin accumulation by cytokinins in *Arabidopsis thaliana*. *Plant Physiol.* **108**, 47–57 (1995).

40. Bauer, P. Regulation of iron acquisition responses in plant roots by a transcription factor. *Biochem. Mol. Biol. Educ.* **44**, 438–49, doi:10.1002/bmb.20967 (2016).

41. Jefferson, R. A., Kavanagh, T. A. & Bevan, M. W. GUS fusions: beta-glucuronidase as a sensitive and versatile gene fusion marker in higher plants. *EMBO J.* **6**, 3901–3907, doi:10.1007/978-1-4939-3115-6_2 (1987).

42. Ben Abdallah, H. & Bauer, P. Quantitative reverse transcription-qPCR-based gene expression analysis in plants. *Methods in Molecular Biology* **1363**, 9–24 (2016).

43. Slep, G. et al. CATMA, a comprehensive genome-scale resource for silencing and transcript profiling of Arabidopsis genes. *BMC Bioinformatics* **8**, 400, doi:10.1186/1471-2105-8-400 (2007).

44. Katari, M. S. et al. Virtualplant: a software platform to support systems biology research. *Plant Physiol.* **152**, 500–15, doi:10.1104/pp.109.147025 (2010).

45. Supek, F., Bosnjak, M., Skunca, N. & Smuc, T. Revigo summarizes and visualizes long lists of gene ontology terms. *PLoS One* **6**, e21800, doi:10.1371/journal.pone.0021800 (2011).

46. Szklarczyk, D. et al. String v10: protein-protein interaction networks, integrated over the tree of life. *Nucleic Acids Res.* **43**, D447–52, doi:10.1093/nar/gku1003 (2015).

## Acknowledgements

The authors express their thanks for excellent technical assistance with plant growth, physiological and gene expression analysis to Elke Wieneke and Angelika Anna, Düsseldorf. The authors are very grateful to Joséé Caius, Saclay, for his technical assistance in transcriptome analysis. This research was funded by Deutsche Forschungsgemeinschaft grants Ba 1610/5-1 and 7-1 to P.B. and a PhD fellowship of the Saarland University to M.A.N.A. The authors are thankful to Tzvetina Brumbarova, Düsseldorf, for critical reading of the manuscript.

## Author Contributions

P.B. conceived the original research plans and supervised the experiments; M.N.A. performed most of the experiments; M.N.A. and P.B. designed the experiments and analyzed the data; J.M. and F.M. generated and initially characterized transgenic lines; S.P. conducted and supervised microarray experiments; C.F.-S. conducted and supervised metal determination; P.B. and M.N.A. wrote the article.

## Additional Information

**Supplementary information** accompanies this paper at doi:10.1038/s41598-017-11171-7

**Competing Interests:** The authors declare that they have no competing interests.

**Accession codes:** CATdb project AU13-07_Ox, GEO accession number GSE93093.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_16.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2017