# Sulphur availability modulates Arabidopsis thaliana responses to iron deficiency

Kevin Robe, Fei Gao, Pauline Bonillo, Nicolas Tissot, Frédéric Gaymard, Pierre Fourcroy, Esther Izquierdo, Christian Dubos

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

## To cite this version:

Kevin Robe, Fei Gao, Pauline Bonillo, Nicolas Tissot, Frédéric Gaymard, et al., Sulphur availability modulates Arabidopsis thaliana responses to iron deficiency. PLoS ONE, 2020, 15 (8), pp.e0237998. 10.1371/journal.pone.0237998. hal-02943146

HAL Id: hal-02943146

https://hal.inrae.fr/hal-02943146

Submitted on 25 May 2021

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Sulphur availability modulates *Arabidopsis thaliana* responses to iron deficiency

Kevin Robe*, Fei Gao*, Pauline Bonillo, Nicolas Tissot, Frédéric Gaymard, Pierre Fourcroy†, Esther Izquierdo, Christian Dubos‡*

BPMP, Univ Montpellier, CNRS, INRAE, Institut Agro, Montpellier, France

才是作者的研究贡献。

† Deceased.

* christian.dubos@inrae.fr

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

## Abstract

Among the mineral nutrients that are required for plant metabolism, iron (Fe) and sulphur (S) play a central role as both elements are needed for the activity of several proteins involved in essential cellular processes. A combination of physiological, biochemical and molecular approaches was employed to investigate how S availability influences plant response to Fe deficiency, using the model plant *Arabidopsis thaliana*. We first observed that chlorosis symptom induced by Fe deficiency was less pronounced when S availability was scarce. We thus found that S deficiency inhibited the Fe deficiency induced expression of several genes associated with the maintenance of Fe homeostasis. This includes structural genes involved in Fe uptake (i.e. *IRT1*, *FRO2*, *PDR9*, *NRAMP1*) and transport (i.e. *FRD3*, *NAS4*) as well as a subset of their upstream regulators, namely *BTS*, *PYE* and the four clade Ib bHLH. Last, we found that the over accumulation of manganese (Mn) in response to Fe shortage was reduced under combined Fe and S deficiencies. These data suggest that S deficiency inhibits the Fe deficiency dependent induction of the Fe uptake machinery. This in turn limits the transport into the root and the plant body of potentially toxic divalent cations such as Mn and Zn, thus limiting the deleterious effect of Fe deprivation.

## Introduction

Plant growth and development depends on both the availability of mineral nutrients present in the soil and the capacity of the root system to mine the rhizosphere in search for these nutrients. Therefore, the ability of plants to acquire and assimilate mineral nutrients impinges on crop productivity and the quality of their derived products [1].

Among the mineral nutrients required for plant metabolism, iron (Fe) and sulphur (S) play a central role as both are needed for the activity of several proteins involved in essential cellular processes (e.g. photosynthesis, respiration, biosynthesis of primary and secondary metabolites). Within proteins, Fe and S are found as individual atoms or tightly connected to each other in the form of prosthetic groups, the Fe-S clusters. In plants, Fe-S proteins are mainly involved in electron transfer reactions (e.g. within the photosynthetic and respiratory electron

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

OPEN ACCESS

**Citation:** Robe K, Gao F, Bonillo P, Tissot N, Gaymard F, Fourcroy P, et al. (2020) Sulphur availability modulates *Arabidopsis thaliana* responses to iron deficiency. PLoS ONE 15(8): e0237998. https://doi.org/10.1371/journal.pone.0237998

**Editor:** Abidur Rahman, Iwate University, JAPAN

**Received:** April 26, 2020

**Accepted:** August 6, 2020

**Published:** August 20, 2020

**Copyright:** © 2020 Robe et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** KR was supported by fellowships from the French National Research Agency (ANR-17-CE20-0008) and the INRAE (department of "Biologie et Amélioration des Plantes"). FG was supported by a fellowship from the China Scholarship Council (CSC). NT was supported by a fellowship from the French "Ministère de l'enseignement supérieur et de la recherche".
**Competing interests:** The authors have declared that no competing interests exist.

transport chains in chloroplasts and mitochondria, respectively) or in the catalysis of oxidation-reduction reactions [2–5]. The biosynthesis of cysteine also connects Fe and S homeostasis, since the desulphuration of cysteine to alanine is a reaction in the assembly process of Fe-S clusters [2]. The biosynthesis of methionine, a sulphur-containing amino acid, is another hub that links Fe and S homeostasis. For instance, methionine is the precursor of molecules that play key roles during the plant response to Fe shortage, such as ethylene, nicotianamine (NA) or phytosiderophores (PS). In non-graminaceous species (strategy I plants for Fe acquisition), ethylene and NA are respectively involved in the activation of the transcriptional regulatory cascade controlling the response to Fe deficiency and the transport of Fe throughout the plant [6–8]. In graminaceous species (strategy II plants for Fe acquisition), PS, which derives from NA, are secreted into the rhizosphere in order to chelate and solubilize Fe and thus improve plant Fe nutrition [8]. Considering that Fe shortage arises in one third of the cultivated land at the surface of the planet and that the frequency of soils displaying S deficiency is increasing (mostly due to a decrease of anthropogenic S release; [9]), it became obvious that there is a need for understanding how Fe and S homeostasis regulation are interconnected in plants. This knowledge being necessary if one aims at sustaining plant growth and productivity in such limited environment with limited use of exogenous fertilizers.

This tight connection between Fe and S, at the molecular and physiological levels, as exemplified earlier through the Fe-S cluster biogenesis and methionine biosynthesis, suggested that coordination in the control of both nutrients homeostasis occurred. Indeed, several examples in the literature confirmed this assumption, whether the plants belong to eudicot or monocot clades. For instance, in tomato (*Solanum lycopersicum* L.) or rapeseed (*Brassica napus* L.), the Fe deficiency induced expression and activity of key genes associated with Fe acquisition was inhibited when plants were submitted to S deprivation [10, 11]. Similarly, the capacity of barley (*Hordeum vulgare* L.) seedlings grown under S deficiency to cope with Fe shortage was restored when S was provided [12]. Conversely, it was also demonstrated, in tomato and durum wheat (*Triticum turgidum* L. subsp. *durum*), that Fe shortage triggers responses associated with S deficiency [13–15]. In *Arabidopsis thaliana*, such cross talk between Fe and S homeostasis was also reported, notably through the study of the expression of the main Fe and S transporter present at the root epidermis, namely *IRT1* (*IRON-REGULATED TRANSPORTER 1*) and *SULTR1;1* (*SULPHATE TRANSPORTER 1;1*), respectively [16]. As expected, these studies also highlighted that Fe and S deficiency impact the expression of a common set of genes that are not directly related to Fe and S acquisition and assimilation, which was in adequacy with the central role that play both nutrients in several key metabolic processes [17–20].

The intricate connection between Fe and S homeostasis suggests that complex regulatory mechanism are at play. For instance, if the regulatory mechanisms controlling Fe and S homeostasis are quite well described [21–23], it remains to determine the means by which the information on the availability of one nutrient is integrated into the regulatory pathway controlling the homeostasis of the other one. The nature of the signals controlling this mechanism is also matter of debates. Because Fe and S are required in large amounts in chloroplasts and mitochondria, in particular in the form of Fe-S clusters, it has been proposed that one or more retrograde signal(s) deriving from these organelles might feedback to the nucleus to regulate the homeostasis of both nutrient [18, 24–27]. Among these potential signals, the SAL1/FRY1 (FIERY)-PAP chloroplast retrograde pathway emerged as an interesting candidate [24, 28]. For instance, PAP (3'-phosphoadenosine 5'- phosphate), an intermediate compound into the assimilation of S, plays a central role in plant response to several environmental stresses (e.g. osmotic, cold, drought, high light or response to pathogens) including nutrient availability
[29–32]. Interestingly, it was also reported that SAL1-PAP retrograde signalling pathway was also involved in the maintenance of Fe homeostasis [33].

In this study we have investigated, using the model plant *Arabidopsis thaliana*, how S deficiency affects at the physiological and molecular levels the plant response to Fe shortage and how this signal is integrated into the transcriptional regulatory cascade controlling Fe deficiency responses. We first observed that chlorosis symptom induced by Fe deficiency was less pronounced when S availability was scarce. Because Fe deficiency is mostly regulated at the transcriptional level, we have investigated the expression of key genes involved in Fe homeostasis, including transcription factors (TFs). Expression analyses revealed that the Fe deficiency induced expression of most of the assayed genes was reduced when both Fe and S were limiting. Nevertheless, only a subset of TFs followed this expression pattern whereas the expression of the other set was unaffected by the double deficiency. Using loss-of-function mutants we found that the SAL1-PAP retrograde signalling pathway is most probably not involved in this process. The analysis of rosette leaves and roots micronutrients content highlighted that the increased accumulation of manganese (Mn) in response of Fe shortage was reduced under the combined Fe and S deficiencies. Altogether these data suggest that S deficiency, by inhibiting the Fe acquisition machinery that is induced in response to Fe deficiency, limits the unspecific transport into the root and the plant body of potentially toxic divalent cations, particularly Mn, thus limiting the deleterious effect of Fe deprivation.

## Results

### Iron deficiency induced chlorosis in *Arabidopsis thaliana* is less pronounced when sulphur availability is scarce

One of the most obvious symptoms that plants display when experiencing iron (Fe) deficiency is the yellowing of the leaf tissues. This phenomenon, also called chlorosis, is due to an altered accumulation of chlorophylls present in the interveinal tissues. In order to investigate how sulphur (S) availability may affect the plant response to Fe deficiency, we assessed how contrasting Fe and S nutrition regimes impact chlorophylls accumulation. For this purpose, three-week-old *Arabidopsis thaliana* plants were subjected for 10 days to Fe deficiency (-Fe +S), S deficiency (+Fe -S) and Fe and S deficiencies (-Fe -S), and compared to plants grown in control condition (+Fe +S).

From the visual screening of the plants (Fig 1A), no striking differences were observed between plants grown in control and S deficiency conditions whereas, as expected, plants grown in Fe deficiency condition displayed a clear chlorosis phenotype. Strikingly, we observed that the chlorosis phenotype observed in plants grown under Fe deprivation was less pronounced when S availability was scarce (-Fe +S vs -Fe -S), in particular for the youngest leaves that developed during the course of the treatments. This observation was confirmed by measuring the chlorophylls content present in the young leaves of the plants grown in these four different conditions (Fig 1B). Anthocyanins content, which increases in tissues of plants facing S deficiency [34], was evaluated in order to assess the effect of S deficiency in our experiments. As expected, we found that anthocyanins accumulation was higher in plants submitted to S deficiency (+Fe -S and -Fe -S) when compared to plants grown in control or Fe deficiency conditions (Fig 1C).

At the molecular level, *IRT1* (*IRON-REGULATED TRANSPORTER 1*) expression, which encodes the high affinity Fe transporter present at the root epidermis, was induced in response to Fe deficiency and this induction was decreased in plants facing both Fe and S deficiencies (Fig 1D), in agreement with [16]. Steady state mRNA level of *FER1* (*FERRITIN 1*), which encodes a protein involved in the transient storage of Fe and whose expression is induced
<figure><img src="image_5.png" /><figcaption>**Fig 1.** *Arabidopsis thaliana* iron deficiency induced chlorosis is less pronounced when sulphur availability is scarce. (A) Rosette phenotype of Arabidopsis plants grown for three weeks in presence of 25 µM Fe(III)-EDTA and then transferred for 10 days in four different media: control (+Fe +S), S deficiency (+Fe -S), Fe deficiency (-Fe +S) and Fe and S deficiencies (-Fe -S). Magenta arrows indicate the youngest leaves that formed during the 10 days of treatments. Bar = 1 cm. (B) Chlorophylls and (C) anthocyanins content of the youngest leaves of the rosette presented panel A (Magenta arrows). (D) Quantitative RT-PCR analysis of *IRT1* (-Fe marker) and *FER1* (+Fe marker) mRNA</figcaption></figure>
levels in plant roots following 5 days of treatment. (E) Quantitative RT-PCR analysis of *SULTR1;1* (-S marker), *APR1*, *APK1* and *APK3* mRNA levels in plant roots following 5 days of treatment. (B-E) Means with the same letter are not significantly different according to one-way ANOVA followed by post-hoc Tukey test (*P*< 0.05). n = 3 biological repeats from one representative experiment. Each experiment was repeated three times.

https://doi.org/10.1371/journal.pone.0237998.g001

under excess Fe, was also measured. *FER1* expression decreased in plants grown under Fe deficiency, independently of the amount of S present in the medium (Fig 1D). *SULTR1;1* (*SULPHATE TRANSPORTER 1;1*) expression, which encodes the high affinity sulphate transporter present at the root epidermis, was induced in response to S deficiency, still in agreement with [16] (Fig 1E). The expression of *APR1* (*APS REDUCTASE 1*), *APK1* (*APS KINASE 1*) and *APK3* (*APS KINASE 3*), was also assayed as the encoded proteins play key roles in the partitioning of the assimilated S, in the form of APS (adenosine 5'-phosphosulfate), between the biosynthesis of cysteine (APR1), and thus methionine, and the formation of several secondary metabolites such as glucosinolates (APK1 and APK3). The expression of *APR1*, and to a lesser extent *APK1*, was similar to that one of *SULTR1;1*, whereas the expression of *APK3* remained unaffected by the treatments (Fig 1E).

We then investigated how Fe and S deficiencies were affecting parameters associated with plant photosynthetic activity. We first analysed the accumulation of key proteins involved in the photosynthetic chain, namely the PHOTOSYSTEM II SUBUNIT A (PsbA), the CYTOCHROME b6/f COMPLEX SUBUNIT (b6/f) and the PHOTOSYSTEM I SUBUNIT D (PsaD). We confirmed that under Fe deficiency the three proteins accumulated less than in Fe sufficiency, in a manner that was independent of S availability (Fig 2A–2C) [35]. In order to get some insights into how the treatments affected the photosynthetic capacities of the plants, we focused our analysis on the youngest leaves, the ones that developed during the course of the treatments. For this purpose, we have measured the maximum quantum yield of PSII (Fv/Fm), which reflects the potential efficiency of PSII (Fig 2D), and the quantum yield (QY), that reflects the electron flux that goes through the photosynthetic chain (Fig 2E). Fv/Fm and QY were unchanged in S deficiency when compared with control condition whereas both were significantly diminished by Fe deficiency. In contrast to Fv/Fm, QY was significantly higher in Fe and S dual deficiency than that in Fe deficiency alone.

## Sulphur deficiency inhibits the induction of the expression, in response to iron deficiency, of key Arabidopsis genes involved in the maintenance of iron homeostasis

The expression of the high affinity Fe transporter *IRT1* in response to Fe deficiency (-Fe +S) is inhibited by S deficiency (-Fe -S) whereas S deficiency alone (+Fe -S) does not significantly affect its expression when compared to control (+Fe +S) condition (Fig 1D). Whether the expression of two other main genes involved in Fe acquisition, namely *FRO2* (*FERRIC REDUCTION OXIDASE 2*) and *PDR9* (*PLEIOTROPIC DRUG RESISTANCE 9/ABCG37*) [36], follow a similar expression pattern remained to be determined. We found that the expression pattern of both genes was comparable to that of *IRT1* (Fig 3A). This was also the case for the expression of the low affinity Fe transporter *NRAMP1* (*NATURAL RESISTANCE-ASSOCIATED MACROPHAGE PROTEIN 1*) (Fig 3A) [37] as well as *FRD3* (*FERRIC REDUCTASE DEFECTIVE 3*) and *NAS4* (*NICOTIANAMINE SYNTHASE 4*), two key genes involved in the transport of Fe through the xylem and phloem conducting tissues, respectively (Fig 3B) [38, 39]. It is noteworthy that S deficiency alone was sufficient to partly inhibit *NRAMP1* expression (Fig 3A). Altogether, these data suggested that the S deficiency signal might be integrated
<figure><img src="image_6.png" /><figcaption>**Fig 2.** Fe deficiency effect on *Arabidopsis thaliana* rosette leaves photosynthetic activity is less pronounced when S availability is scarce. (A-C) Western blot analyses using total protein extracts from whole rosette leaves of Arabidopsis plants grown for three weeks in presence of 25 µM Fe(III)-EDTA and then transferred for 10 days in four different media: control (+Fe +S), S deficiency (+Fe -S), Fe deficiency (-Fe +S) and Fe and S deficiencies (-Fe -S). Three proteins related to the photosynthetic light reactions were analysed: (A) PsbA (photosystem II), (B) b6f (cytochrome b6f complex) and (C) PsaD (photosystem I). (D-E) Chlorophyll fluorescence parameters: (D) Fv/Fm of dark-adapted plants and (E) quantum yield of photochemical energy conversion. (A-E) Means with the same letter are not significantly different according to one-way ANOVA followed by post-hoc Tukey test (*P*< 0.05). n = 3 biological repeats from one representative experiment. Each experiment was repeated three times.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0237998.g002

early in the transcriptional regulatory cascade governing the response to Fe deficiency in
Arabidopsis.

We thus examined the expression pattern of the main genes involved in the transcriptional
regulatory cascade that govern the Arabidopsis response to Fe deficiency [21, 22, 40–46]. The
main part of this regulatory network involves 10 bHLH transcription factors (TFs) belonging
to four different clades, namely clade Ib (i.e. bHLH38, bHLH39, bHLH100 and bHLH101), IIIa
(FIT, FE-DEFICIENCY INDUCED TRANSCRIPTION FACTOR/bHLH29), IVb (PYE,
POPEYE/bHLH47) and IVc (bHLH34, bHLH104, bHLH105/ILR3 and bHLH115), and two
homologous R2R3-MYBs (MYB10 and MYB72). Upstream this molecular network, BTS
(BRUTUS) encodes a Fe-binding E3 ligase that mediates the degradation, through the 26S pro-
teasome, of the clade IVc bHLH TFs [47]. We found that the induction of the expression of all
the tested genes, but clade IVc bHLHs and FIT, in response to Fe deficiency was inhibited by a
low S availability (Figs 4 and 5). This later observation was in agreement with the pattern of
<figure><img src="image_7.png" /><figcaption>**Fig 3.** Expression analysis of genes involved in Fe uptake and transport in *Arabidopsis thaliana*. cDNA were synthesized from RNA extracted from roots of Arabidopsis plants grown for three weeks in presence of 25 µM Fe(III)-EDTA and transferred to four different media: control (+Fe +S), S deficiency (+Fe -S), Fe deficiency (-Fe +S) and Fe and S deficiencies (-Fe -S). Root samples were harvested five days after plants were transferred to the assayed media. Quantitative RT-PCR analyses were carried out on genes involved in (A) Fe uptake (*FRO2*, *PDR9* and *NRAMP1*) and (B) Fe transport (*FRD3* and *NAS4*). (A-B) Means with the same letter are not significantly different according to one-way ANOVA followed by post-hoc Tukey test (*P*< 0.05) n = 3 biological repeats from one representative experiment. Each experiment was repeated three times.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0237998.g003

expression of all the structural genes involved in the maintenance of Fe homeostasis that were assayed (i.e. *IRT1*, *FRO2*, *PDR9*, *FRD3* and *NAS4*) and *NRAMP1* (Fig 3).

**The impact of sulphur scarcity on the chlorosis symptoms associated with iron deficiency do not involve SAL1-PAP in Arabidopsis**

The way how S deficiency influences plant response to Fe deficiency is still to be discovered. Recent studies have evidenced that PAP (3'-phosphoadenosine 5'- phosphate), an intermediate compound into the assimilation of S, plays a central role in plant response to several environmental stresses (e.g. osmotic, cold, drought, high light, response to pathogens or cadmium excess) [30–32], in the chloroplast retrograde signalling pathway [24] and in the maintenance of Fe homeostasis [33]. PAP is formed from PAPS (3'-phosphoadenosine 5'-phosphosulfate) after the transfer of the sulphate group to target molecules (e.g. glucosinolates). PAP is an inhibitor of nuclear RNase (i.e. XRN2, XRN3 and XRN4, 5'-3' exonucleases) and affects the expression of several stress inducible genes. PAP accumulation is regulated by the phosphatase
<figure><img src="image_8.png" /><figcaption>**Fig 4.** Expression analysis of genes involved in the early steps of the transcriptional regulatory network that controls *Arabidopsis thaliana* response to Fe deficiency. Quantitative RT-PCR analyses were carried out using cDNA synthesized from RNA extracted from roots of Arabidopsis plants grown for three weeks in presence of 25 µM Fe(III)-EDTA and transferred to four different media: control (+Fe +S), S deficiency (+Fe -S), Fe deficiency (-Fe +S) and Fe and S deficiencies (-Fe -S). Root samples were harvested five days after plants were transferred to the assayed media. The assayed genes encode a Fe-binding E3 ligase (*BTS*, *BRUTUS*), six bHLH (*FIT/bHLH29*, *bHLH34*, *bHLH104*, *bHLH105* and *bHLH115*) transcription factors. Means with the same letter are not significantly different according to one-way ANOVA followed by post-hoc Tukey test (*P*< 0.05). n = 3 biological repeats from one representative experiment. Each experiment was repeated three times.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0237998.g004

SAL1/FRY1 (FIERY1) in the chloroplast where it is catabolised into AMP and inorganic phosphate (Pi) [28]. In order to determine whether the inhibition of the Fe deficiency symptom when plants are subjected to both Fe and S deficiencies could be due to PAP, we have grown two loss-of-function mutant alleles of SAL1/FRY1 (i.e. *fry1-6* and *alx8*) as well as the triple *xrn2 xrn3 xrn4* mutant in the presence or absence of Fe and/or S and compared them to wild type plants (S1 Fig). No differences in term of chlorosis (S1A Fig) and chlorophylls accumulation pattern (S1B Fig), in response to the four nutritional regimes that were assayed, were found between wild type plants and the three mutants. This result indicated that the cross talk occurring between Fe and S deficiency responses is most probably not involving PAP in Arabidopsis.

## Sulphur deficiency inhibits the iron deficiency induced accumulation of manganese in Arabidopsis plant tissues

Since leaf chlorophylls content is a proxy of the amount of Fe present in the tissues, we investigated whether Fe content in rosette leaves was increased when plants were facing Fe and S
<figure><img src="image_9.png" /><figcaption>**Fig 5.** Expression analysis of genes involved in late steps of the transcriptional regulatory network that controls *Arabidopsis thaliana* response to Fe deficiency. Quantitative RT-PCR analyses were carried out using cDNA synthesized from RNA extracted from roots of Arabidopsis plants grown for three weeks in presence of 25 µM Fe(III)-EDTA and transferred to four different media: control (+Fe +S), S deficiency (+Fe -S), Fe deficiency (-Fe +S) and Fe and S deficiencies (-Fe -S). Root samples were harvested five days after plants were transferred to the assayed media. The assayed genes encode five bHLH (*PYE/bHLH47: POPEYE*, *bHLH38*, *bHLH39*, *bHLH100* and *bHLH101*) and two R2R3-MYB (*MYB10* and *MYB72*) transcription factors. Means with the same letter are not significantly different according to one-way ANOVA followed by post-hoc Tukey test (*P*< 0.05). n = 3 biological repeats from one representative experiment. Each experiment was repeated three times.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0237998.g005
<figure><img src="image_10.png" /><figcaption>**Fig 6.** The increased accumulation of manganese in *Arabidopsis thaliana* plants in response to Fe deficiency is inhibited when S availability is scarce. (A) Fe, (B) Zn (zinc) and (C) Mn (manganese) content of rosette leaves and roots of Arabidopsis plants grown for three weeks in presence of 25 µM Fe(III)-EDTA and then transferred for 10 days in four different media: control (+Fe +S), S deficiency (+Fe -S), Fe deficiency (-Fe +S) and Fe and S deficiencies (-Fe -S). Means with the same letter are not significantly different according to one-way ANOVA followed by post-hoc Tukey test (*P*< 0.05). n = 3 biological repeats from one representative experiment. Each experiment was repeated three times.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0237998.g006

deficiencies when compared to plants facing Fe deficiency alone. We first observed that, under Fe deficiency, the amount of Fe present in the aerial tissues was lower than that of plants grown in control condition or under S deficiency (Fig 6A). We also found that S availability was not affecting the amount of Fe present in the aerial tissues, whether the plants were grown in presence or absence of Fe. Similar observations were made in roots.

We then investigated whether the observed phenotype could be due to some differential accumulation of zinc (Zn) and/or manganese (Mn) when plants are grown under Fe deficiency or under Fe and S deficiencies. This is because it is well established that in addition to transport Fe, IRT1 has the ability to transport Zn and Mn [48]. This is also because NRAMP1 is, in addition to being a low affinity Fe transporter, the major high affinity Mn transporter in Arabidopsis [49]. This hypothesis was supported by the expression patterns observed for *IRT1* and
NRAMP1 in Fe deficiency compared to Fe and S dual deficiency (Figs 1D and 3). We found that Zn accumulation was increased in both rosette leaves and roots when plants were grown under Fe deficiency and that S deficiency was not reverting this over-accumulation (Fig 6B). In contrast, we found that the increased accumulation of Mn in rosette leaves and roots associated with Fe deficiency was reverted when S availability was scarce (Fig 6C).

These observations suggested that part of the Fe deficiency induced chlorosis symptoms might be due to an increased accumulation of Mn in plant tissues and additional divalent cations such as Zn. In order to test this hypothesis, plants were grown under Fe deficiency alone or under Fe deficiency with decreased concentration (i.e. 1/2, 1/5 and 1/10 of the initial concentration) of Mn, Zn or Mn and Zn. Fe deficiency induced chlorosis symptoms were similar when plants were grown under Fe deficiency or under Fe deficiency with decreased concentration of Mn or Zn (S2 Fig). In contrast, plants grown under Fe deficiency with reduced Mn and Zn concentrations in the growth medium displayed decreased Fe deficiency induced chlorosis symptoms that were positively correlated with the dilution factor (Fig 7A). Mn and Zn contents in rosette leaves were inversely correlated with chlorophylls content.

To investigate the effect of divalent cations accumulation on Fe deficiency induced chlorosis symptoms, Arabidopsis plants were then grown under Fe deficiency at different density (1, 4, 8 and 16 plants / 400 ml growth solution). This experiment showed that the Fe deficiency induced chlorosis was decreasing when the plant density was increasing, phenocopying the effect of S deprivation (Fig 7B). Mn and Zn contents in rosette leaves were also inversely correlated with chlorophylls content in this experiment.

Altogether these data suggest that part of the Fe deficiency induced chlorosis symptoms observed in Arabidopsis might be due to an accumulation of Mn, Zn and most probably other divalent cations. These data also suggest that the inhibition of the Fe deficiency induced chlorosis, when Arabidopsis plants are grown in the absence of S, could be due to an inhibition of, at least, Mn accumulation in plant tissues.

Discussion

Iron (Fe) and sulphur (S) are two essential mineral nutrients for plant growth and development whose availability in soils affects crops productivity and the quality of their derived products [1]. Within the plant cell, Fe and S are tightly connected. The biogenesis of the Fe-S clusters, which are prosthetic groups essential for several proteins involved in key metabolic processes (e.g. photosynthesis, respiration), best exemplifies this interaction [2, 3, 5]. On a wider scale, Fe deficiency limits plant growth on one third of the cultivated land at the surface of the planet whereas, concomitantly, the surface of soils displaying S deficiency is increasing (because of a decrease of S release related to anthropogenic activities; [9]). Such observations indicate that the occurrence of Fe and S dual deficiency is likely to increase in the coming years and that sustaining plant growth and productivity in such limited environment, with lowering the use of exogenous fertilizers, will necessitate decrypting the mechanisms that govern plant Fe and S homeostasis interconnection. Pioneer work, carried out at the physiological and molecular levels, was initiated in order to document Fe and S interaction in different dicots (i.e. Arabidopsis, rapeseed, tomato) and mococots (i.e. barley, durum wheat, wheat) species [10–20]. Altogether, these studies have highlighted that modifications in the availability of one nutrient affects the homeostasis of the other, by notably modulating the expression and activity of several genes and proteins involved in these processes.

In order to get further insight into the interconnection between Fe and S homeostasis, we investigated, using Arabidopsis as model plant, how S deficiency signal was integrated into the transcriptional regulatory cascade controlling Fe deficiency responses [21, 22]. For this
<figure><img src="image_11.png" /><figcaption>**Fig 7.** Manganese and zinc availability participates to the Fe deficiency induced chlorosis. (A) Effect of Mn and Zn availability on rosette leaves Fe deficiency induced chlorosis. Left panel: rosette phenotype of Arabidopsis plants grown for three weeks in presence of 25 µM Fe(III)-EDTA and then transferred for 10 days in Fe deficiency (-Fe +S) condition with various concentration of Mn and Zn. Mn and Zn concentration in Fe deficiency condition was 1/2, 1/5 and 1/10 of the initial -Fe +S medium. Bar = 1cm. Upper right panels: chlorophylls content of the youngest leaves of the rosettes presented upper left panel. Lower panels: Fe, Zn and Mn content of rosette leaves. (B) Effect of plant density on rosette leaves Fe deficiency induced chlorosis. Upper left panel: rosette phenotype of Arabidopsis plants grown for three weeks in presence of 25 µM Fe(III)-EDTA and then transferred for 10 days in Fe deficiency (-Fe +S) condition. Plants were grown at the density of 1, 4, 8 and 16 plants per 400 ml media. Bar = 1cm. Upper right panel: chlorophylls content of the youngest leaves of the rosettes presented upper left</figcaption></figure>
panel. Lower panels: Fe, Zn and Mn content of rosette leaves. (A-B) Means with the same letter are not significantly different according to one-way ANOVA followed by post-hoc Tukey test (P< 0.05). *: significant differences, t-test p < 0.05. n = 3 biological repeats from one representative experiment. Experiments displayed panel A and B were repeated three and two times, respectively.

https://doi.org/10.1371/journal.pone.0237998.g007

purpose, plants were submitted for 10 days to Fe, S or Fe and S deficiencies. We first observed, mostly in the youngest rosette leaves, that Fe deficiency induced chlorosis was less pronounced when plants were also subjected to S deprivation (Fig 1A and 1B). This observation was not in agreement with previous studies reporting that chlorophyll content of plant grown under both Fe and S deficiency was similar or lower to that of plant grown under Fe deprivation [10–13, 19]. However, similar observations to the ones described in this study were recently made in tomato [50]. In contrast, we found that S deficiency was not inducing any significant chlorosis, as described in wheat [13, 19] and tomato [50]. Nevertheless, S deficiency was reported to induce chlorosis in other plant species (i.e. barley, rapeseed, tomato). The use of specific physiological (anthocyanin accumulation) and molecular (i.e. SULTR1;1, APR1, APK1 and APK3 genes) markers confirmed that in our experiments the plants perceived the application of S deficiency (Fig 1C and 1E) and thus suggested that some species- or developmental stage-specific (i.e. seedlings vs mature plants, expanding vs mature leaves) mechanisms could be at play in the plant response to Fe, S or Fe and S deficiencies. Indeed, it could not be excluded that the growth parameters (e.g. day length, composition of the growth medium, stress application) may also influence these responses. However, as expected, we found a strong correlation between the chlorosis state of the plants and their photosynthetic capacities [10, 35] (Figs 1A, 1B and 2).

At the molecular level, the analysis of *IRT1*, *FRO2* and *PDR9* (three key genes involved in the Fe uptake machinery in Arabidopsis) mRNA abundance in roots unravelled common patterns of expression in response to Fe and S availability (Figs 1D and 3A). Their expression was strongly induced in response to Fe deficiency when compared to control condition and this induction was reduced when S availability was scarce (dual Fe and S deficiency). Such pattern of expression was previously reported for *IRT1* and *FRO2* in different species [10, 11, 16] even if this expression profile may fluctuate depending on the experimental conditions [15, 17]. Strikingly, although *IRT1*, *FRO2* and *PDR9* mRNA levels appeared lowered in S deficiency when compared to control condition, none of the observed differences was statistically significant. This is in contrast with the above-mentioned studies where *IRT1* and/or *FRO2* expression was either increased [10, 15] or decreased [11, 16, 17] in response to S deficiency when compared to control condition. Several reasons may explain these differences. In Arabidopsis, *IRT1* and *FRO2* expression was shown to decrease in response to S deficiency [11, 16, 17]. Such discrepancies might be due to the more stringent statistical test that was used in our study. Indeed, it cannot be excluded that it could also be the result of differences in the growth conditions that were used, in particular the composition of the growth medium. Indeed, one may expect variations that could be plant species dependent, even if similar discrepancies were reported in tomato [10, 15]. Our experiment revealed also that the expression pattern of *FRD3* and *NAS4* (two key genes involved in the transport of Fe through the xylem and phloem conducting tissues, respectively) was similar to that of *IRT1* (Fig 3B). This later result indicated that S deficiency modulates in a coordinated manner the expression of genes involved in Fe acquisition (*IRT1*, *FRO2* and *PDR9*) and transport (*FRD3* and *NAS4*), indicating that S deficiency inhibits Arabidopsis response to Fe shortage by modulating the expression of genes encoding regulatory proteins.

This hypothesis relies on the fact that the regulatory network controlling the Arabidopsis response to Fe shortage is well described [21, 22, 44, 45, 58] (Fig 8). Upstream this network,
<figure><img src="image_12.png" /><figcaption>**Fig 8.** Schematic representation of the Fe deficiency response pathway in *Arabidopsis thaliana* and its connection with S shortage. Arabidopsis response to Fe shortage is essentially controlled at the transcriptional level and involves a complex regulatory network among which bHLH transcription factor (TF), from different clades [57] play a key role. The uppermost bHLH TF in this pathway belong to the clade IVc. The expression of these TF is induced in response to Fe deficiency. The stability of clade IVc encoded proteins is regulated by a Fe-binding E3 ligase, BTS (BRUTUS), whose expression is induced by clade IVc in response to Fe shortage. The expression of *PYE* (*POPEYE*, *bHLH47*), a transcriptional repressor, is also induced by clade IVc bHLH. PYE represses its target genes by</figcaption></figure>
interacting with bHLH105/ILR3 (IAA-LEUCINE RESISTANT3) [40, 46, 47]. In addition, clade IVc bHLH act as positive regulators [44, 58], in particular of clade Ib bHLH TF expression. Clade Ib bHLH act as transcriptional activators by forming heterodimers with another key bHLH TF, namely FIT (FE-DEFICIENCY INDUCED TRANSCRIPTION FACTOR, bHLH29). The expression of FIT is also induced in response to Fe shortage (indirect activation by the clade IVc bHLH TF). Downstream this pathway, two R2R3-MYB whose expression is dependent on FIT and clade IVc bHLH, MYB10 and MYB72, also act as transcriptional activators. The combined action of PYE, clade Ib bHLH TF, FIT, MYB10 and MYB72 contribute to maintain Fe homeostasis when Fe availability is scarce by regulating the expression of structural genes involved in Fe uptake, transport and storage. In green are the genes for which the Fe deficiency induced expression is reduced when both Fe and S deficiencies occurs.

https://doi.org/10.1371/journal.pone.0237998.g008

four basic helix-loop-helix (clade IVc: bHLH34, bHLH104, bHLH105/ILR3 and bHLH115) transcription factors (TFs) directly activate the expression, when Fe availability is scarce, of a second set of bHLH TFs, namely *bHLH38*, *bHLH39*, *bHLH100*, *bHLH101* (calde Ib) and *bHLH47/PYE*. In parallel, clade IVc bHLHs also induce (indirect regulation) the expression of *bHLH29/FIT* [51]. FIT-dependent protein complexes are then formed, through heterodimerization with clade Ib bHLHs, leading to the direct activation of the expression of structural genes involved in the maintenance of Fe homeostasis such as *IRT1* and *FRO2*. Among the indirect targets of FIT, MYB10 and MYB72 (two R2R3-MYB TFs) play also an important role in maintaining plant growth under limited Fe availability notably by modulating the PDR9-dependent Fe uptake process [52]. In contrast, PYE acts as a repressor of genes especially involved in Fe transport such as *NAS4* [40, 53]. PYE activity relies on its interaction with ILR3 [40, 46, 47]. BTS, a Fe-binding E3 ligase, plays a key role in regulating this transcriptional regulatory cascade by modulating the stability of clade IVc bHLH TFs [40, 46, 47]. The study of the expression pattern of the different member of this regulatory network in response to Fe and/or S deficiency highlighted two types of mRNA accumulation profiles. First, we found that clade IVc bHLHs and *FIT* expression was unaffected by S deficiency, whether plants were subjected to Fe deficiency or not (Fig 5). The second pattern of expression was similar to that one described earlier for *IRT1* and concerned the downstream targets of clade IVc bHLHs (i.e. clade Ib bHLH and *PYE*) and *BTS* (Fig 4). It is noteworthy that the expression of *FIT*, *PYE*, *bHLH39* and *bHLH100* in response to Fe or Fe and S dual deficiency when compared to control condition was similar to what was reported in a previous study focusing on Arabidopsis [17]; however, in our study, and unlike what was previously reported, the inhibitory effect of S deficiency alone on the expression of these four genes was, like for *IRT1* and *FRO2*, not statistically significant. The inhibition by S deficiency of the Fe deficiency induced expression of the most downstream TFs in this regulatory network is intriguing (Fig 7). One would have expected such inhibition to occur upstream in this regulatory cascade, at the clade IVc bHLH and FIT levels. A hypothesis that could explain such observation would be that these TF play additional roles in the context of S deficiency. For instance, bHLH105/ILR3 was proposed to have a wider role in the plant response to various biotic and abiotic stresses [54]. Another hypothesis would be that clade IVc bHLHs and FIT regulation occurs at the posttranslational level, leading to their degradation via the 26S proteasome. In support of this latter hypothesis, some studies have highlighted the importance of this mechanism in the regulation of FIT [41, 55]. This would also involve a BTS independent mechanism for clade IVc bHLHs, as it is the case for MdbHLH104, an ortologue TF from apple [56]. Additional studies will be required to fully address this question in both Arabidopsis and other dicot crops. Altogether our data demonstrate that S deficiency integrates into the Arabidopsis Fe deficiency transcriptional regulatory network by regulating the expression of specific TFs (Fig 7) and most probably involves either additional transcriptional regulators competing with clade IVc bHLHs and FIT (repressors) or prosttranscriptional regulatory mechanisms affecting the stability of *BTS*, clade Ib
bHLHs, PYE, MYB10 and MYB72. Nevertheless, it is unlikely that this process involves the chloroplast SAL1-PAP retrograde signalling pathway (S1 Fig).

Because we observed that Fe deficiency induced chlorosis was less pronounced when plants were facing Fe and S dual deficiency, one could expect that Fe remobilisation into the aerial part would be increased. Measurements of Fe content in rosette leaves indicated that S deficiency was not influencing the amount of Fe measured in these tissues (Fig 6A). Similarly, we found that the amount of Fe present in roots was not dependent on S availability. Such pattern of Fe accumulation was in agreement with FER1 expression (Fig 1D). Interestingly, the effect of S deficiency on the ability of dicot plants to accumulate Fe seems to be variable. For instance, in rapeseed, it was reported that S deficiency inhibits the accumulation of Fe in leaves and roots [10]. In tomato, the effect of S deficiency on the accumulation of Fe in leaves was reported to vary from a strong reduction (stronger than Fe deficiency itself) to no effect (as reported for Arabidopsis in this study) [11, 15, 18]. These observations indicate that the impact of S deficiency on Fe accumulation might be specific to the plant species, the developmental stage and/or the growth conditions. Another hypothesis to explain that Fe deficiency induced chlorosis was less pronounced when plants were facing Fe and S dual deficiency would be that the amount of other transition metals present in the plant tissues might fluctuate depending on S availability, and thus affect plant Fe deficiency induced chlorosis. In support of this hypothesis, it was recently demonstrated that excess of zinc (Zn), and to a lesser extent manganese (Mn), were phenocopying the Fe deficiency induced responses at the physiological and molecular levels [59]. In addition, it has been reported that Mn excess negatively impacts PSI activity but barely affects PSII [60]. The effect that we observed in the photosynthetic capacities indicates that the dual Fe and S deficiency reduced Mn toxicity since it improved QY (representative of total electron flux) but had no effect on Fv/Fm (that measures PSII photochemical efficiency) [60] (Fig 2D and 2E). Another argument in support of this hypothesis is that the expression of IRT1 (a low affinity Zn and Mn transporter; [48]) and NRAMP1 (the major high affinity Mn transporter; [49]) followed, in our experiments, a similar expression pattern (i.e. lower induction in response to Fe deficiency when S availability is scarce) (Figs 1D and 3A). We first found that Zn and Mn content were significantly increased in rosette leaves of plants grown under Fe deficiency, when compared to plants grown in control condition, and that S deficiency alone was not affecting the amount of both metals, matching previous report observations [17, 59] (Fig 6B and 6C). In addition, we found that the amount of Mn present in rosette leaves of plants grown under Fe and S dual deficiency was lower than that of plants grown under Fe deficiency alone, and comparable to control condition (Fig 6C). Interestingly, we found that reducing the amount of Mn and Zn present in the growth medium of plants subjected to Fe deficiency had a positive effect on the Fe deficiency induced chlorosis (Fig 7A). Similarly, increasing the density of plants subjected to Fe deficiency was also having a positive effect on the Fe deficiency induced chlorosis that was accompanied by a decrease of Mn and Zn in rosette leaves (Fig 7B). Altogether, these results suggested that the effect of S deprivation on the Fe deficiency symptoms could be due to a differential accumulation of Mn and other divalent cations such as Zn in rosette leaves.

Data presented in this study show that S deficiency, by inhibiting part of the Arabidopsis Fe deficiency transcriptional regulatory network, negatively regulates the Fe uptake machinery that is induced in response to Fe deficiency (Fig 8). This process most probably limits the unspecific transport into the root and the plant body of potentially toxic divalent cations such as Mn and Zn by IRT1 and NRAMP1, thus limiting the deleterious effect of Fe deprivation. Indeed, it cannot be excluded that yet undiscovered transport mechanisms might also be at play. Whether S deficiency modulates the Fe deficiency transcriptional regulatory network in a similar manner in other dicots, in particular crops, will have to be investigated. This work
emphasizes also the necessity to study the intricate connection that exists in plants between the homeostasis of macro- and micronutrients if one aims at sustaining plant growth and productivity in nutrients limited environment without the use of exogenous fertilizers. This is also of importance if one aims at improving nutritional and health quality of plant derived products by increasing specific micronutrient content, such as Fe, without accumulating undesired toxic metals.

# Materials and methods

## Plant material

*Arabidopsis thaliana* ecotype Columbia (Col-0) was used as wild type. The following mutant lines were used in this study: *fry1-6* (N520882), *alx8* (N66977) and the *xrn2 xrn3 xrn4* triple mutant [29].

## Growth conditions

Plants were hydroponically grown in a controlled growth chamber (23°C, 70% relative humidity, 200 µmol.m⁻².s⁻¹ light intensity) under short-day condition (8h light/16h dark photoperiod) as described elsewhere [36]. One plant per 100 ml Hoagland solution (S1 Table) density was used unless otherwise stated. Following 3 weeks of growth, in presence of 25 µM Fe(III)-EDTA, nutrient deficiencies were applied for 10 days. Nutrient deficiencies included: absence of iron (-Fe), absence of sulphur (-S), absence of iron and sulphur (-Fe -S). Following 5 days of nutrient deficiencies, roots were sampled and frozen in liquid nitrogen for subsequent genes expression analysis (qRT-PCR). This time point was chosen because it corresponds to the maximum of expression of genes involved in both Fe uptake and partitioning and the regulation of these processes [36]. Samples used for determining chlorophylls, anthocyanins and ions (i.e. iron: Fe; zinc: Zn; manganese: Mn) content as well as for western blot and chlorophyll fluorescence analysis were collected at the end of the experiment, after 10 days of nutrient deficiencies.

## Biochemical analyses

**Chlorophylls content.** Chlorophylls from 5 leaf discs (diameter: 0.35 cm) were extracted in 1 ml 100% acetone in the dark under agitation. The absorbance (A) at 661.8 and 644.8 nm was then measured. Total chlorophylls content was assessed using the following equations: $Chl a + Chl b = 7.05*A_{661.6} + 18.09*A_{644.8}$ and expressed as µg/cm² [61].

**Anthocyanins content.** Anthocyanins from 50 mg of leaves (FW: fresh weight) homogenized in 1 ml acidic methanol (1% HCl, w/v) were extracted in the dark at 4°C (moderate shaking, overnight). Following centrifugation (5 min, 14.000 rpm, room temperature), the absorbance of the supernatant (A) at 530 and 657 nm was measured. Total anthocyanins content was assessed using the following equations: $Anthocyanins = (A_{530} - (0.25*A_{657}))$ and expressed as $A_{Anthocyanins}$/mg FW.

**Iron, manganese and zinc determination.** About 15 mg (DW: dry weight) of ground sample (rosette leaves or roots) were mixed with 750 µl nitric oxide (65%) and 250 µl hydrogen peroxide 30% before homogenization. Samples were left over night at room temperature and then mineralized 7h at 84°C. Once mineralized, the nitric oxide proportion present in the samples was adjusted to 5 to 10% of the final volume by adding ultrapure water. Minerals content present in the samples was then measured by microwave plasma atomic emission spectroscopy (MP-AES, Agilent).
## Chlorophyll fluorescence

A kinetic imaging fluorometer (FluorCam FC 800-O, Photon Systems Instruments, Czech Republic) was used to capture chlorophyll fluorescence images and to estimate the maximal quantum yield [F<sub>v</sub>/F<sub>m</sub> = (F<sub>m</sub>−F<sub>o</sub>)/F<sub>m</sub>] and the effective quantum yield of the photosystem II (φ<sub>PSII</sub>; under 200 µmol·m<sup>-2</sup>·s<sup>-1</sup> actinic light). Plants were first dark-adapted for 30 min and saturating flashes (0.8 s, 5000 µmol·m<sup>-2</sup>·s<sup>-1</sup>) were then shot to measure the maximum fluorescence. Measurements were performed at room temperature.

## Gene expression analysis (qRT-PCR)

Total RNAs was extracted using the Tri-Reagent (Molecular Research Center) method. Briefly, each sample was homogenized in 1 ml Tri-Reagent solution mixed with 160 µl of chloroform: isoamyl alcohol (24:1). Following centrifugation (10 min, 13.000 rpm, 4°C) total RNAs present in the aqueous phase were precipitated by the addition of 400 µl of isopropanol followed by another centrifugation. Pellets were then washed twice with ethanol 70% and dried prior resuspension in RNAse-free water. For each sample, 1 µg of total RNA treated with DNase was reverse transcribed into cDNA using the RevertAid kit (Thermo scientific). Quantitative PCR analyses were carried out using a LightCycler® 480 (Roche) and the LC480-SYBR-Green master I reaction mix (Roche). PP2AA3 (*PROTEIN PHOSPHATASE 2A SUBUNIT A3*) was used as reference gene [62]. Expression levels were calculated using the comparative threshold cycle method. All the primers used are described S2 Table.

## Western blot analysis

Total proteins were extracted from 100 mg of samples (rosette leaves) grinded in liquid nitrogen and homogenized in 250 µl of 1x extraction buffer (100 mM Tris-HCl pH 8, 1% LDS, 25 mM EDTA, 5% β-mercaptoethanol (v/v), 2 mM PMSF). Following two consecutive centrifugations the supernatant was collected and stored on ice prior use. Protein content was determined using the Bradford assay [63]. Proteins were then separated by 12% SDS-PAGE (Laemmli system) and transferred to a PVDF membrane. Membranes were blocked for overnight in TBS buffer (50 mM Tris pH 7.5, 150 mM NaCl, 0.1% Tween-20, 1% BSA) and then incubated 1h with primary antibody and then 1 h with secondary antibody, both diluted in TBS containing 1% BSA. Dilutions of primary antibodies applied were: rabbit anti-PsaD 1:10.000 [64], rabbit anti-b6/f 1:5.000 [64], rabbit anti-PsbA 1:20.000 (Agrisera, AS05084) and goat anti-rabbit HRP conjugated (Promega, W4011) 1:10.000. Immunodetection (LAS-3000 Imaging System, Fuji) was performed using the Luminata Forte (Deutscher) solution.

## Studied Arabidopsis gene IDs

APK1/AKN1, At2g14750; APK3, At3g03900; APR1/PRH19, At4g04610; bHLH29/FIT, At2g28160; bHLH34, At3g23210; bHLH38, At3g56970; bHLH39, At3g56980; bHLH47/PYE, At3g47640; bHLH100, At2g41240; bHLH101, At5g04150; bHLH104, At4g14410; bHLH105/ILR3, At5g54680; bHLH115, At1g51070; BTS, At3g18290; FER1/AtFER1, At5g01600; FRD3, At3g08040; FRO2, At1g01580; FRY1/SAL1/ALX8/HOS2/RON1/RON1/SUPO1, At5g63980; IRT1, At4g19690; MYB10, At3g12820; MYB72, At1g56160; NAS4, At1g56430; NRAMP1, At1g80830; PDR9/ABCG37, At3g53480; PP2AA3, At1g13320; SULTR1;1, At4g08620; XRN2, At5g42540; XRN3, At1g75660; XRN4, At1g54490.
**Supporting information**

S1 Fig. SAL1-PAP retrograde signal is not involved in the attenuation of *Arabidopsis thaliana* Fe deficiency induced chlorosis when S availability is scarce. (A) Rosette phenotype of Arabidopsis plants grown for three weeks in presence of 25 µM Fe(III)-EDTA and then transferred for 10 days in four different media: control (+Fe +S), S deficiency (+Fe -S), Fe deficiency (-Fe +S) and Fe and S deficiencies (-Fe -S). Wild type (WT) plant and three mutants were analysed, namely two loss-of-function alleles of the phosphatase SAL1/FRY1 (*fry1-6* and *alx8*) and a triple mutant of the nuclear 5'-3' exoribonucleases XRN2, XRN3 and XRN4 (*xrn2 xrn3 xrn4*) whose activity is regulated by SAL1/FRY1. Bar = 1 cm. (B) Chlorophylls content of the youngest leaves of the rosette presented panel A. Means with the same letter are not significantly different according to one-way ANOVA followed by post-hoc Tukey test (*P*< 0.05). n = 3 biological repeats from one representative experiment. Each experiment was repeated three times.

(PDF)

S2 Fig. Effect of manganese and zinc deficiency on rosette leaves Fe deficiency induced chlorosis. Arabidopsis plants were grown for three weeks in presence of 25 µM Fe(III)-EDTA and then transferred for 10 days in Fe deficiency (-Fe +S) condition alone or with a concentration of Mn or Zn that was 1/10 of the initial -Fe +S medium. Bar = 1cm.

(PDF)

S1 Table. Composition of the Hoagland solutions used in this study.

(PDF)

S2 Table. Primers used in this study.

(PDF)

**Acknowledgments**

We thank Drs Berger N and Touraine B for help in preparing this manuscript as well as Bou-
cherez J and Maghiaoui A for excellent technical support. *xrn2 xrn3 xrn4* triple mutant seeds
were kindly provided by Dr Vaucheret H (IJPB, INRA Versailles, France). This paper is dedi-
cated to the Memory of Fourcroy P. Fourcroy P passed away before the submission of the final
version of this manuscript. Dubos C accepts responsibility for the integrity and validity of the
data collected and analyzed.

**Author Contributions**

**Conceptualization:** Kevin Robe, Fei Gao, Pauline Bonillo, Nicolas Tissot, Frédéric Gaymard, Pierre Fourcroy, Esther Izquierdo, Christian Dubos.

**Funding acquisition:** Frédéric Gaymard, Christian Dubos.

**Investigation:** Kevin Robe, Fei Gao, Pauline Bonillo, Nicolas Tissot, Pierre Fourcroy, Esther Izquierdo, Christian Dubos.

**Methodology:** Kevin Robe, Fei Gao, Pauline Bonillo, Nicolas Tissot, Pierre Fourcroy, Esther Izquierdo, Christian Dubos.

**Writing – original draft:** Kevin Robe, Fei Gao, Esther Izquierdo, Christian Dubos.

**Writing – review & editing:** Kevin Robe, Fei Gao, Esther Izquierdo, Christian Dubos.
References

1. Briat JF, Dubos C, Gaymard F. Iron nutrition, biomass production, and plant product quality. Trends Plant Sci. 2015; 20: 33–40. https://doi.org/10.1016/j.tplants.2014.07.005 PMID: 25153038

2. Couturier J, Touraine B, Briat JF, Gaymard F, Rouhier N. The iron-sulfur cluster assembly machineries in plants: current knowledge and open questions. Front Plant Sci. 2013; 4:259. https://doi.org/10.3389/fpls.2013.00259 PMID: 23898337

3. Berger N, Vignols F, Przybyla-Toscano J, Roland M, Rofidal V, Touraine B, et al. Identification of client iron-sulfur proteins of the chloroplastic NFU2 transfer protein in Arabidopsis thaliana. J Exp Bot. 2020; eraa166. https://doi.org/10.1093/jxb/eraa166 PMID: 32240305

4. Roland M, Przybyla-Toscano J, Vignols F, Berger N, Azam T, Christ L, et al. The plastidial Arabidopsis thaliana NFU1 protein binds and delivers [4Fe-4S] clusters to specific client proteins. J Biol Chem. 2020; 295: 1727–1742. https://doi.org/10.1074/jbc.RA119.011034 PMID: 31911438

5. Touraine B, Vignols F, Przybyla-Toscano J, Ischebeck T, Dhalleine T, Wu HC, et al. Iron-sulfur protein NFU2 is required for branched-chain amino acid synthesis in Arabidopsis roots. J Exp Bot. 2019; 70: 1875–1889. https://doi.org/10.1093/jxb/erz050 PMID: 30785184

6. Brumbarova T, Bauer P, Ivanov R. Molecular mechanisms governing Arabidopsis iron uptake. Trends Plant Sci. 2015; 20: 124–133. https://doi.org/10.1016/j.tplants.2014.11.004 PMID: 25499025

7. Gayomba SR, Zhai Z, Jung HI, Vatamaniuk OK. Local and systemic signaling of iron status and its interactions with homeostasis of other essential elements. Front Plant Sci. 2015; 6:716. https://doi.org/10.3389/fpls.2015.00716 PMID: 26442030

8. Kobayashi T, Nishizawa NK. Iron uptake, translocation, and regulation in higher plants. Annu Rev Plant Biol. 2012; 63: 131–152. https://doi.org/10.1146/annurev-arplant-042811-105522 PMID: 22404471

9. Scherer HW. Sulphur in crop production. Eur J Agron. 2001; 14: 81–111.

10. Muneer S, Lee BR, Kim KY, Park SH, Zhang Q, Kim TH. Involvement of sulphur nutrition in modulating iron deficiency responses in photosynthetic organelles of oilseed rape (Brassica napus L.). Photosynth Res. 2014; 119: 319–329. https://doi.org/10.1007/s11120-013-9953-8 PMID: 24264737

11. Zuchi S, Cesco S, Varanini Z, Pinton R, Astolfi S. Sulphur deprivation limits Fe-deficiency responses in tomato plants. Planta. 2009; 230: 85–94. https://doi.org/10.1007/s00425-009-0919-1 PMID: 19350269

12. 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. https://doi.org/10.1093/jxb/erp346 PMID: 20018904

13. Ciaffi M, Paolacci AR, Celletti S, Catarcione G, Kopriva S, Astolfi S. Transcriptional and physiological changes in the S assimilation pathway due to single or combined S and Fe deprivation in durum wheat (Triticum durum L.) seedlings. J Exp Bot. 2013; 64: 1663–1675. https://doi.org/10.1093/jxb/ert027 PMID: 23390290

14. Paolacci AR, Celletti S, Catarcione G, Hawkesford MJ, Astolfi S, Ciaffi M. Iron deprivation results in a rapid but not sustained increase of the expression of genes involved in iron metabolism and sulfate uptake in tomato (Solanum lycopersicum L.) seedlings. J Integr Plant Biol. 2014; 56: 88–100. https://doi.org/10.1111/jipb.12110 PMID: 24119307

15. Zuchi S, Watanabe M, Hubberten HM, Bromke M, Osorio S, Fernie AR, et al. The Interplay between Sulfur and Iron Nutrition in Tomato. Plant Physiol. 2015; 169: 2624–2639. https://doi.org/10.1104/pp.15.00995 PMID: 26438787

16. Forieri I, Wirtz M, Hell R. Toward new perspectives on the interaction of iron and sulfur metabolism in plants. Front Plant Sci. 2013; 4:357. https://doi.org/10.3389/fpls.2013.00357 PMID: 24106494

17. Forieri I, Sticht C, Reichelt M, Gretz N, Hawkesford MJ, Malagoli M, et al. System analysis of metabolism and the transcriptome in Arabidopsis thaliana roots reveals differential co-regulation upon iron, sulfur and potassium deficiency. Plant Cell Environ. 2017; 40: 95–107. https://doi.org/10.1111/pce.12842 PMID: 27726154

18. Vigani G, Pii Y, Celletti S, Maver M, Mimmo T, Cesco S, et al. Mitochondria dysfunctions under Fe and S deficiency: is citric acid involved in the regulation of adaptive responses? Plant Physiol Biochem. 2018; 126: 86–96. https://doi.org/10.1016/j.plaphy.2018.02.022 PMID: 29514113

19. Zamboni A, Celletti S, Zenonia S, Astolfi S, Varaninia Z. Root physiological and transcriptional response to single and combined S and Fe deficiency in durum wheat. Environmental and Experimental Botany. 2017; 143: 172–184.

20. Mendoza-Cozatl DG, Gokul A, Carelse MF, Jobe TO, Long TA, Keyster M. Keep talking: crosstalk between iron and sulfur networks fine-tunes growth and development to promote survival under iron limitation. J Exp Bot. 2019; 70: 4197–4210. https://doi.org/10.1093/jxb/erz290 PMID: 31231775
21. 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. 2016; 170: 2478-2493. https://doi.org/10.1104/pp.15.01827 PMID: 26921305

22. Liang G, Zhang H, Li X, Ai Q, Yu D. bHLH transcription factor bHLH115 regulates iron homeostasis in Arabidopsis thaliana. J Exp Bot. 2017; 68:1743-1755. https://doi.org/10.1093/jxb/erx043 PMID: 28369511

23. Takahashi H, Kopriva S, Giordano M, Saito K, Hell R. Sulfur assimilation in photosynthetic organisms: molecular functions and regulations of transporters and assimilatory enzymes. Annu Rev Plant Biol. 2011; 62: 157-184. https://doi.org/10.1146/annurev-arplant-042110-103921 PMID: 21370978

24. Estavillo GM, Crisp PA, Pornsiriwong W, Wirtz M, Collinge D, Carrie C, et al. Evidence for a SAL1-PAP chloroplast retrograde pathway that functions in drought and high light signaling in Arabidopsis. Plant Cell. 2011; 23: 3992-4012. https://doi.org/10.1105/tpc.111.091033 PMID: 22128124

25. Hell R, Wirtz M. Molecular Biology, Biochemistry and Cellular Physiology of Cysteine Metabolism in Arabidopsis thaliana. Arabidopsis Book. 2011; 9:e0154. https://doi.org/10.1199/tab.0154 PMID: 22303278

26. Vigani G, Briat JF. Impairment of Respiratory Chain under Nutrient Deficiency in Plants: Does it Play a Role in the Regulation of Iron and Sulfur Responsive Genes? Front Plant Sci. 2015; 6:1185. https://doi.org/10.3389/fpls.2015.01185 PMID: 26779219

27. Vigani G, Zocchi G, Bashir K, Philippar K, Briat JF. Signals from chloroplasts and mitochondria for iron homeostasis regulation. Trends Plant Sci. 2013; 18: 305-311. https://doi.org/10.1016/j.tplants.2013.01.006 PMID: 23462548

28. Estavillo GM, Chan KX, Phua SY, Pogson BJ. Reconsidering the nature and mode of action of metabolite retrograde signals from the chloroplast. Front Plant Sci. 2012; 3:300. https://doi.org/10.3389/fpls.2012.00300 PMID: 23316207

29. Hirsch J, Misson J, Crisp PA, David P, Bayle V, Estavillo GM, et al. A novel fry1 allele reveals the existence of a mutant phenotype unrelated to 5'->3' exoribonuclease (XRN) activities in Arabidopsis thaliana roots. PLoS One. 2011; 6:e16724. https://doi.org/10.1371/journal.pone.0016724 PMID: 21304819

30. Ishiga Y, Watanabe M, Ishiga T, Tohge T, Matsuura T, Ikeda Y, et al. The SAL-PAP Chloroplast Retrograde Pathway Contributes to Plant Immunity by Regulating Glucosinolate Pathway and Phytohormone Signaling. Mol Plant Microbe Interact. 2017; 30: 829-841. https://doi.org/10.1094/MPMI-03-17-0055-R PMID: 28703028

31. Pornsiriwong W, Estavillo GM, Chan KX, Tee EE, Ganguly D, Crisp PA, et al. A chloroplast retrograde signal, 3'-phosphoadenosine 5'-phosphate, acts as a secondary messenger in abscisic acid signaling in stomatal closure and germination. Elife. 2017; 6:e23361. https://doi.org/10.7554/eLife.23361 PMID: 28323614

32. Xi H, Xu H, Xu W, He Z, Ma M. A SAL1 Loss-of-Function Arabidopsis Mutant Exhibits Enhanced Cadmium Tolerance in Association with Alleviation of Endoplasmic Reticulum Stress. Plant Cell Physiol. 2016; 57: 1210-1219. https://doi.org/10.1093/pcp/pcw069 PMID: 27044671

33. Balparda M, Armas AM, Estavillo GM, Roschzttardtz H, Pagani MA, Gomez-Casati DF. The PAP/SAL1 retrograde signaling pathway is involved in iron homeostasis. Plant Mol Biol. 2020; 102: 323-337. https://doi.org/10.1007/s11103-019-00950-7 PMID: 31900819

34. Nikiforova VJ, Kopka J, Tolstikov V, Fiehn O, Hopkins L, Hawkesford MJ, et al. Systems rebalancing of metabolism in response to sulfur deprivation, as revealed by metabolome analysis of Arabidopsis plants. Plant Physiol. 2005; 138: 304-318. https://doi.org/10.1104/pp.104.053793 PMID: 15834012

35. Hantzis LJ, Kroh GE, Jahn CE, Cantrell M, Peers G, Pilon M, et al. A Program for Iron Economy during Deficiency Targets Specific Fe Proteins. Plant Physiol. 2018; 176: 596-610. https://doi.org/10.1104/pp.17.01497 PMID: 29150559

36. Fourcroy P, Tissot N, Gaymard F, Briat JF, Dubos C. Facilitated Fe Nutrition by Phenolic Compounds Excreted by the Arabidopsis ABCG37/PDR9 Transporter Requires the IRT1/FRO2 High-Affinity Root Fe(2+) Transport System. Mol Plant. 2016; 9: 485-488. https://doi.org/10.1016/j.molp.2015.09.010 PMID: 26415695

37. Curie C, Alonso JM, Le Jean M, Ecker JR, Briat JF. Involvement of NRAMP1 from Arabidopsis thaliana in iron transport. Biochem J. 2000; 347: 749-755. PMID: 10769179

38. Durrett TP, Gassmann W, Rogers EE. The FRD3-mediated efflux of citrate into the root vasculature is necessary for efficient iron translocation. Plant Physiol. 2007; 144: 197-205. https://doi.org/10.1104/pp.107.097162 PMID: 17351051

39. Schuler M, Rellan-Alvarez R, Fink-Straube C, Abadia 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. 2012; 24: 2380-2400. https://doi.org/10.1105/tpc.112.099077 PMID: 22706286
40. Long TA, Tsukagoshi H, Busch W, Lahner B, Salt DE, Benfey PN. The bHLH transcription factor POPEYE regulates response to iron deficiency in Arabidopsis roots. Plant Cell. 2010; 22: 2219–2236.
https://doi.org/10.1105/tpc.110.074096 PMID: 20675571

41. Sivitz AB, Hermand V, Curie C, Vert G. Arabidopsis bHLH100 and bHLH101 control iron homeostasis via a FIT-independent pathway. PLoS One. 2012; 7:e44843. https://doi.org/10.1371/journal.pone.0044843 PMID: 22984573

42. Wang HY, Klatte M, Jakoby M, Baumlein H, Weisshaar B, Bauer P. Iron deficiency-mediated stress regulation of four subgroup Ib BHLH genes in Arabidopsis thaliana. Planta. 2007; 226: 897–908. https://doi.org/10.1007/s00425-007-0535-x PMID: 17516080

43. Zhang J, Liu B, Li M, Feng D, Jin H, Wang P, et al. The bHLH transcription factor bHLH104 interacts with IAA-LEUCINE RESISTANT3 and modulates iron homeostasis in Arabidopsis. Plant Cell. 2015; 27: 787–805. https://doi.org/10.1105/tpc.114.132704 PMID: 25794933

44. Gao F, Robe K, Bettembourg M, Navarro N, Rofidal V, Santoni V, et al. The Transcription Factor bHLH121 Interacts with bHLH105 (ILR3) and Its Closest Homologs to Regulate Iron Homeostasis in Arabidopsis. Plant Cell. 2020; 32: 508–524. https://doi.org/10.1105/tpc.19.00541 PMID: 31776233

45. Gao F, Robe K, Gaymard F, Izquierdo E, Dubos C. The Transcriptional Control of Iron Homeostasis in Plants: A Tale of bHLH Transcription Factors? Front Plant Sci. 2019; 10:6. https://doi.org/10.3389/fpls.2019.00006 PMID: 30713541

46. Tissot N, Robe K, Gao F, Grant-Grant S, Boucherez J, Bellegarde F, et al. Transcriptional integration of the responses to iron availability in Arabidopsis by the bHLH factor ILR3. New Phytol. 2019; 223: 1433–1446. https://doi.org/10.1111/nph.15753 PMID: 30773647

47. Matthiadis A, Long TA. Further insight into BRUTUS domain composition and functionality. Plant Signal Behav. 2016; 11:e1204508. https://doi.org/10.1080/15592324.2016.1204508 PMID: 27359166

48. Vert G, Grotz N, Dedaldechamp F, Gaymard F, Guerinot ML, Briat JF, et al. IRT1, an Arabidopsis transporter essential for iron uptake from the soil and for plant growth. Plant Cell. 2002; 14: 1223–1233. https://doi.org/10.1105/tpc.001388 PMID: 12084823

49. Cailliatte R, Schikora A, Briat JF, Mari S, Curie C. High-affinity manganese uptake by the metal transporter NRAMP1 is essential for Arabidopsis growth in low manganese conditions. Plant Cell. 2010; 22: 904–917. https://doi.org/10.1105/tpc.109.073023 PMID: 20228245

50. Astolfi S, Pii Y, Mimmo T, Lucini L, Miras-Moreno MB, Coppa E, et al. Single and Combined Fe and S Deficiency Differentially Modulate Root Exudate Composition in Tomato: A Double Strategy for Fe Acquisition? Int J Mol Sci. 2020; 21:E4038. https://doi.org/10.3390/ijms2114038 PMID: 32516916

51. Colangelo EP, Guerinot ML. The essential basic helix-loop-helix protein FIT1 is required for the iron deficiency response. Plant Cell. 2004; 16: 3400–3412. https://doi.org/10.1105/tpc.104.024315 PMID: 15539473

52. Zamioudis C, Hanson J, Pieterse CM. beta-Glucosidase BGLU42 is a MYB72-dependent key regulator of rhizobacteria-induced systemic resistance and modulates iron deficiency responses in Arabidopsis roots. New Phytol. 2014; 204: 368–379. https://doi.org/10.1111/nph.12980

53. Klatte M, Schuler M, Wirtz M, Fink-Straube C, Hell R, Bauer P. The analysis of Arabidopsis nicotianamine synthase mutants reveals functions for nicotianamine in seed iron loading and iron deficiency responses. Plant Physiol. 2009; 150: 257–271. https://doi.org/10.1104/pp.109.136374 PMID: 19304929

54. Samira R, Li B, Kliebenstein D, Li C, Davis E, Gillikin JW, et al. The bHLH transcription factor ILR3 modulates multiple stress responses in Arabidopsis. Plant Mol Biol. 2018; 97: 297–309. https://doi.org/10.1007/s11103-018-0735-8 PMID: 29882068

55. Cui Y, Chen CL, Cui M, Zhou WJ, Wu HL, Ling HQ. Four IVa bHLH Transcription Factors Are Novel Interactors of FIT and Mediate JA Inhibition of Iron Uptake in Arabidopsis. Mol Plant. 2018; 11: 1166–1183. https://doi.org/10.1016/j.molp.2018.06.005 PMID: 29960107

56. Zhao Q, Ren YR, Wang QJ, Wang XF, You CX, Hao YJ. Ubiquitination-Related MdBT Scaffold Proteins Target a bHLH Transcription Factor for Iron Homeostasis. Plant Physiol. 2016; 172: 1973–1988. https://doi.org/10.1104/pp.16.01323 PMID: 27660166

57. Heim MA, Jakoby M, Werber M, Martin C, Weisshaar B, Bailey PC. The basic helix-loop-helix transcription factor family in plants: a genome-wide study of protein structure and functional diversity. Mol Biol Evol. 2003; 20: 735–747. https://doi.org/10.1093/molbev/msg088 PMID: 12679534

58. Gao F, Robe K, Dubos C. Further insights into the role of bHLH121 in the regulation of iron homeostasis in Arabidopsis thaliana. Plant Signal Behav. 2020.

59. Leskova A, Giehl RFH, Hartmann A, Fargasova A, von Wiren N. Heavy Metals Induce Iron Deficiency Responses at Different Hierarchic and Regulatory Levels. Plant Physiol. 2017; 174: 1648–1668. https://doi.org/10.1104/pp.16.01916 PMID: 28500270
60. Millaleo R, Reyes-Diaz M, Alberdi M, Ivanov AG, Krol M, Huner NP. Excess manganese differentially inhibits photosystem I versus II in Arabidopsis thaliana. J Exp Bot. 2013; 64: 343–354. https://doi.org/10.1093/jxb/ers339 PMID: 23183256

61. Lichtenthaler FW. Karl Freudenberg, Burckhardt Helferich, Hermann O. L. Fischer: a centennial tribute. Carbohydr Res. 1987; 164: 1–22.

62. Czechowski T, Stitt M, Altmann T, Udvardi MK, Scheible WR. Genome-wide identification and testing of superior reference genes for transcript normalization in Arabidopsis. Plant Physiol. 2005; 139: 5–17. https://doi.org/10.1104/pp.105.063743 PMID: 16166256

63. Bradford MM. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. 1976; 72: 248–254. https://doi.org/10.1006/abio.1976.9999 PMID: 942051

64. Touraine B, Boutin JP, Marion-Poll A, Briat JF, Peltier G, Lobreaux S. Nfu2: a scaffold protein required for [4Fe-4S] and ferredoxin iron-sulphur cluster assembly in Arabidopsis chloroplasts. Plant J. 2004; 40: 101–11. https://doi.org/10.1111/j.1365-313X.2004.02189.x PMID: 15361144