# The involvement of OsPHO1;1 in the regulation of iron transport through integration of phosphate and zinc deficiency signalling

Chorpet saenchai¹, Nadia Bouain², Mushtak Kisko², Chanakan Prom-u-thai¹, Patrick Doumas², Hatem Rouached², ¹*

¹Agronomy Division, Department of Plant and Soil Sciences, Faculty of Agriculture, Chiang Mai University, Thailand, ²INRA, Biochimie et Physiologie Moléculaire des Plantes Research Unit, INRA/CNRS/Montpellier University/Montpellier SupAgro, 2, Place Pierre Viala 34060 Cedex 2, France

**Submitted to Journal:**

Frontiers in Plant Science

**Specialty Section:**

Plant Biotechnology

**ISSN:**

1664-462X

**Article type:**

Original Research Article

**Received on:**

20 Feb 2016

**Accepted on:**

14 Mar 2016

**Provisional PDF published on:**

14 Mar 2016

**Frontiers website link:**

www.frontiersin.org

**Citation:**

Saenchai C, Bouain N, Kisko M, Prom-u-thai C, Doumas P and Rouached H (2016) The involvement of OsPHO1;1 in the regulation of iron transport through integration of phosphate and zinc deficiency signalling. *Front. Plant Sci.* 7:396. doi:10.3389/fpls.2016.00396

**Copyright statement:**

© 2016 Saenchai, Bouain, Kisko, Prom-u-thai, Doumas and Rouached. This is an open-access article distributed under the terms of the [Creative Commons Attribution License (CC BY)](https://creativecommons.org/licenses/by/4.0/). The use, distribution and reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
# The Involvement of OsPHO1;1 in the Regulation of Iron Transport Through Integration of Phosphate and Zinc Deficiency Signaling.

Chorpet Saenchai, Nadia Bouain, Mushtak Kisko, Chanakan Prom-U-Thai, Patrick Doumas, Hatem Rouached

► To cite this version:

Chorpet Saenchai, Nadia Bouain, Mushtak Kisko, Chanakan Prom-U-Thai, Patrick Doumas, et al. The Involvement of OsPHO1;1 in the Regulation of Iron Transport Through Integration of Phosphate and Zinc Deficiency Signaling. Frontiers in Plant Science, 2016, 7, pp.396. 10.3389/fpls.2016.00396. hal-01594473

HAL Id: hal-01594473

https://hal.science/hal-01594473

Submitted on 26 Sep 2017

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
This Provisional PDF corresponds to the article as it appeared upon acceptance, after peer-review. Fully formatted PDF and full text (HTML) versions will be made available soon.
# The Involvement of OsPHO1;1 in the Regulation of Iron Transport Through Integration of Phosphate and Zinc Deficiency Signaling

Chorpet Saenchai¹,²†, Nadia Bouain¹†, Mushtak Kisko¹, Chanakan Prom-u-thai², Patrick Doumas¹ and Hatem Rouached¹,²*

¹ Biochimie et Physiologie Moléculaire des Plantes Research Unit, Institut National de la Recherche Agronomique - Centre National de la Recherche Scientifique - Montpellier University, Montpellier, France, ² Agronomy Division, Department of Plant and Soil Sciences, Faculty of Agriculture, Chiang Mai University, Chiang Mai, Thailand

Plants survival depends on their ability to cope with multiple nutrient stresses that often occur simultaneously, such as the limited availability of essential elements inorganic phosphate (Pi), zinc (Zn), and iron (Fe). Previous research has provided information on the genes involved in efforts by plants to maintain homeostasis when a single nutrient (Pi, Zn, or Fe) is depleted. Recent findings on nutritional stress suggest that plant growth capacity is influenced by a complex tripartite interaction between Pi, Zn, and Fe homeostasis. However, despite its importance, how plants integrate multiple nutritional stimuli into complex developmental programs, and which genes are involved in this tripartite (Pi ZnFe) interaction is still not clear. The aim of this study was to examine the physiological and molecular responses of rice (Oriza sativa L.) to a combination of Pi, Zn, and/or Fe deficiency stress conditions. Results showed that Fe deficiency had the most drastic single-nutrient effect on biomass, while the Zn deficiency-effect depended on the presence of Pi in the medium. Interestingly, the observed negative effect of Fe starvation was alleviated by concomitant Pi or PiZn depletion. Members of the OsPHO1 family showed a differential transcriptional regulation in response PiZnFe combinatory stress conditions. Particularly, the transcripts of the OsPHO1;1 sense and its natural antisense cis-NatPHO1;1 showed the highest accumulation under PiZn deficiency. In this condition, the Ospho1;1 mutants showed over-accumulation of Fe in roots compared to wild type plants. These data reveal coordination between pathways involved in Fe transport and PiZn signaling in rice which involves the OsPHO1;1, and support the hypothesis of a genetic basis for Pi, Zn, and Fe signaling interactions in plants.

**Keywords:** ions homeostasis, iron, phosphate, rice, signaling crosstalk, zinc

## INTRODUCTION

Phosphorus (P) is an essential nutrient for optimal plant growth, development, and productivity. The effect of inorganic phosphate (Pi) deficiency on crop yield has become a worldwide concern in recent years due to the pressing problem of food availability (Holford, 1997; Abelson, 1999; Neset and Cordell, 2012). Accordingly, this issue has received increasing attention. To date, research has

**OPEN ACCESS**

**Edited by:**
Nokwanda Makunga,
Stellenbosch University, South Africa

**Reviewed by:**
Prateek Tripathi,
University of Southern California, USA
David Secco,
Plant Energy Biology, Australia

*Correspondence:
Hatem Rouached
hatem.rouached@supagro.inra.fr

† These authors have contributed
equally to this work.

**Specialty section:**
This article was submitted to
Plant Biotechnology,
a section of the journal
Frontiers in Plant Science

Received: 20 February 2016

Accepted: 14 March 2016

Published: 06 April 2016

**Citation:**
Saenchai C, Bouain N, Kisko M, Prom-u-thai C, Doumas P and Rouached H (2016) The Involvement of OsPHO1;1 in the Regulation of Iron Transport Through Integration of Phosphate and Zinc Deficiency Signaling. *Front. Plant Sci.* 7:396. doi: 10.3389/fpls.2016.00396
The involvement of OsPHO1;1 in the regulation of iron transport through integration of
phosphate and zinc deficiency signalling

Chorpet Saenchai¹,²,*, Nadia Bouain¹,*, Mushtak Kisko¹, Chanakan Prom-u-thai², Patrick Doumas¹, Hatem Rouached ¹,²*

¹INRA, Biochimie et Physiologie Moléculaire des Plantes Research Unit,
INRA/CNRS/Montpellier University/Montpellier SupAgro, 2, Place Pierre Viala 34060
Montpellier Cedex 2, France

² Agronomy Division, Department of Plant and Soil Sciences, Faculty of Agriculture, Chiang Mai University, Chiang Mai, Thailand.

*Correspondance: Dr. Hatem ROUACHED
UMR Biochimie & Physiologie Moléculaire des Plantes - (B&PMP)
CNRS-INRA-SUPAGRO-UM, Montpellier, 34060 France
hatem.rouached@supagro.inra.fr
Fax : 0033 (0) 4 67 52 57 37

¶: Authors equally contributed to this work.
provided a comprehensive view of the adaptive strategies used by plants to cope with Pi deficiency (Rouached et al., 2010a; Thibaud et al., 2010; Nussaume et al., 2011; Lin et al., 2014). Lists of genes differentially expressed in Pi deficient conditions compared to control conditions have been produced, and signaling pathways involved in maintaining the homeostasis of Pi have been proposed (Misson et al., 2005; Lin et al., 2008, 2014). However, the practical application of such knowledge is hindered by complex crosstalk linking Pi nutrition and nutrition of other essential micronutrients, e.g., elements are not assimilated, but are taken up zinc (Zn) and iron (Fe) (Huang et al., 2000; Bouain et al., 2014a; Khan et al., 2014). Such interconnections may account for the shortcomings of current agronomic models that typically focus on improving the assimilation of individual elements (Rouached et al., 2011a; Ova et al., 2015). Yet despite their fundamental importance, the molecular bases, biological significance, and agronomical repercussions of these interactions remain unknown.

Our current understanding of genes involved in the transport of Pi into and within the plant has considerably progressed in recent decades (Poirier and Bucher, 2002; Bayle et al., 2011; Nussaume et al., 2011; Sun et al., 2012; Poirier and Jung, 2015). The first step of Pi transport in plants is its acquisition at the soil-root interface. Genes belonging to the *PHT1* family, which are involved in Pi uptake and transport, have been cloned, and characterized in both monocotyledons and dicotyledons (Nussaume et al., 2011; Sun et al., 2012). The increased interest in genes involved in the second rate-limiting step in Pi transport within the plant, i.e., its transfer to xylem vessels, led to the identification of *PHO1* gene in *Arabidopsis* as a key player in this process (Hamburger et al., 2002; Secco et al., 2010). *Atpho1* mutants showed defects in Pi homeostasis in the form a marked reduction in Pi transfer from the roots to the shoot, which was accompanied by low-shoot and high-root Pi concentrations (Hamburger et al., 2002; Rouached et al., 2011b; Poirier and Jung, 2015). The PHO1 protein feature is characterized by the presence of two distinct domains, named SPX and EXS, which may play a role in either Pi transport/ sensing and sorting proteins to endomembranes, respectively (Wang et al., 2004; Wege et al., 2016). Thus, the PHO1 protein does not share protein homology with any other characterized Pi transporters (Wang et al., 2004). The *Arabidopsis* genome contains 11 members of the PHO1 family (Wang et al., 2004). So far, only three members, AtPHO1, AtPHO1;H1 (Stefanovic et al., 2007) and AtPHO1;H3 (Khan et al., 2014), have been reported to be involved in Pi root-to-shoot translocation. These data indicate that the role of members of the *Arabidopsis* PHO1 family is not restricted to Pi transport and homeostasis. For example, the *PHO1;H4* homolog is involved in the response of hypocotyls to blue light (Kang and Ni, 2006), as well as in seed size (Zhou et al., 2009) and flowering (Zhou and Ni, 2009). The role of the *PHO1* gene family has been investigated in a few crop species, including rice (*Oryza sativa* L.) (Secco et al., 2010). Three rice *PHO1* genes, named *OsPHO1;1*, *OsPHO1;2*, and *OsPHO1;3* have been identified. So far, only *OsPHO1;2* has been shown to play a key role in the transfer of Pi from roots to shoots and regulated by Pi deficiency (Secco et al., 2010). The role of the other two

rice *PHO1* genes, *OsPHO1;1*, and *OsPHO1;3*, remains to be investigated.

Pi is known to interact with many micronutrients including Zn and Fe (for a review, see Bouain et al., 2014a; Briat et al., 2015). The effects of Pi, Zn, or Fe availability on plant physiology and transcriptomic have typically been studied individually (Misson et al., 2005; Van De Mortel et al., 2006; Zheng et al., 2009). Although this is an important way to address the question, it can only provide a partial view of the situation. The interactions between Pi, Zn, and Fe metabolisms have been demonstrated both in graminaceous and non-graminaceous plants (Zheng et al., 2009; Briat et al., 2015). These interactions have started being documented at a molecular level, showing that Fe or Zn deficiency modifies the expression of genes involved in Pi transport and assimilation, and vice-versa (Zheng et al., 2009; Briat et al., 2015). Nevertheless, the characterization of these interactions is still in its infancy. In addition to the well-established antagonistic relationship between Fe and Zn in plants (Haydon et al., 2012; Shanmugam et al., 2012; Briat et al., 2015), the existence of an interaction between the homeostasis of Pi and Zn (Huang et al., 2000; Bouain et al., 2014b; Khan et al., 2014; Ova et al., 2015) and Pi and Fe (Zheng et al., 2009) was recently reported, pointing to the existence of a complex tripartite PZnFe signaling interaction to regulate plant growth. Therefore, undertaking an integrative study to clarify the exact molecular mechanisms that coordinate the Pi, Fe, and Zn deficiency pathways and thus shape the plants development during multiple PZnFe deficiency stress is of great interest.

Rice (*Oryza sativa* L.) is one of the most important cereal crops for human consumption, and feeds about half the world population<fnref n="1" />. The effect of individual Pi, Zn, or Fe deficiencies on rice shoot and root growth is fairly well documented (e.g., Wissuwa et al., 2006; Zheng et al., 2009; Secco et al., 2010), but the effect of their combined stresses on plant biomass has been the subject of fewer investigations. In this context, the primary aim of the present study was to broaden our understanding of the long-term effect of different combinations of stresses caused by a deficiency in these essential macro- and micronutrients in rice. Wild type plants and selected mutant lines affected in members of *OsPHO1* gene family, namely *Ospho1;1*, were included in this work. Pi, Zn, and Fe ion concentrations and biomass were determined in rice plants grown in seven different combinations of growing conditions. Physiological and gene expression analyses revealed the involvement of *OsPHO1;1* in the interaction between Fe transport and PiZn signaling pathways in rice.

## RESULTS

### Effects of a Combination of Stresses Caused by Pi, Zn, and Fe Deficiency on Biomass in Rice (Nipponbare)

Wild type rice plants (variety Nipponbare) were grown hydroponically in seven different growing conditions including

<footnote n="1">¹http://faostat3.fao.org/browse/Q/QC/E</footnote>
**Abstract.**

Plants survival depends on their ability to cope with multiple nutrient stresses that often occur simultaneously, such as the limited availability of essential elements inorganic phosphate (Pᵢ), zinc (Zn) and iron (Fe). Previous research has provided information on the genes involved in efforts by plants to maintain homeostasis when a single nutrient (Pᵢ, Zn or Fe) is depleted. Recent findings on nutritional stress suggest that plant growth capacity is influenced by a complex tripartite interaction between Pᵢ, Zn and Fe homeostasis. However, despite its importance, how plants integrate multiple nutritional stimuli into complex developmental programs, and which genes are involved in this tripartite (Pᵢ ZnFe) interaction is still not clear. The aim of this study was to examine the physiological and molecular responses of rice (Oriza sativa L.) to a combination of Pᵢ, Zn and/or Fe deficiency stress conditions. Results showed that Fe deficiency had the most drastic single-nutrient effect on biomass, while the Zn deficiency-effect depended on the presence of Pᵢ in the medium. Interestingly, the observed negative effect of Fe starvation was alleviated by concomitant Pᵢ or PᵢZn depletion. Members of the OsPHO1 family showed a differential transcriptional regulation in response PᵢZnFe combinatory stress conditions. Particularly, the transcripts of the OsPHO1;1 sense and its natural antisense cis-NatPHO1;1 showed the highest accumulation under PᵢZn deficiency. In this condition, the Ospho1;1 mutants showed over-accumulation of Fe in roots compared to wild type plants. These data reveal coordination between pathways involved in Fe transport and PᵢZn signalling in rice which involves the OsPHO1;1, and support the hypothesis of a genetic basis for Pᵢ, Zn and Fe signalling interactions in plants.

**Keywords:** Ions homeostasis, iron, phosphate, rice, signaling crosstalk, zinc.

**Running head:** Involvement of *PHO1;H1* in PiZnFe nutrition interaction in rice
single (−Pi, −Zn, and −Fe), double (−Pi−Fe, −Pi−Zn, and
−Zn−Fe), and triple (−Pi−Fe−Zn) ion deficiency stresses.
**Figure 1** shows the phenotypes of 21-days old wild type
rice plants. The most severe effect on shoot growth was
observed for Fe deficiency in comparison to the effects
of single Pi or Zn deficiency stress. Interestingly, the Fe
deficiency stress was alleviated when combined with Pi limitation
(−Pi−Fe or −Pi−Fe−Zn) (**Figure 1**). This is reminiscent
of data showing that the striking aspect of the inhibition
of primary root elongation in *Arabidopsis* was recovered
when the concentrations of Fe and Pi in medium were
simultaneously reduced (Ward et al., 2008). This also suggests
that manipulating Pi availability to a rice plant could be a
valuable strategy for improving a plant's ability to tolerate Fe
deficiency. To gain quantitative insights into these combined
stresses effect, the rice biomass was assessed by measuring
the roots and shoots dry weight (**Figures 2A,B**). Stress
caused by Fe deficiency had the strongest negative influence
on the shoot and root dry biomass, whereas Zn deficiency
had a slight yet significant positive effect (**Figures 2A,B**).

Pi limitation alone had no effect on these parameters.
Change in Pi together with changes in Fe or Zn had
striking effects. Simultaneous Pi and Fe stress (PiFe) alleviated
the severe stress caused by the Fe deficiency alone, while
stress caused by a simultaneous deficiency in Pi and Zn
had a negative effect. Double (FeZn) or triple (PiFeZn)
deficiencies had a negative effect on biomass similar to the
effect of Fe deficiency alone. Shoot to root dry weight
ratio showed that except for -Pi stress, all other nutrient
stress conditions tested had a negative effect on biomass
(Figure 2C).

Taken together, the results indicate that supply of one
mineral element has a sometimes strong effect on the plant
response to supply of other elements. For example, as shown
in **Figure 1**, Fe deficient plants were smaller and chlorotic,
but Fe deficient plants also lacking P and Zn looked healthier
in comparison. These data underline thus the importance
of the coordination between Pi, Zn, and Fe homeostasis
in the control of rice growth capacity and shoot and root
biomass.

<figure><img src="image_2.png" /><figcaption>**FIGURE 1** | Phenotypes of wild type (WT) rice plant and *Ospho1;1-1* and *Ospho1;1-2* mutant lines grown under individual and combinatory ion deficiency stresses. Seedlings were grown for 21 days on Yoshida medium with different ionic concentrations: Complete Yoshida media (Ct) or lacking either individual phosphate (−Pi), zinc (Zn), or iron (Fe); double elements −Pi−Zn, −Pi−Fe, −Zn−Fe, or triple elements −Pi−Fe−Zn. Photographs were taken at 21 days old. The experiment was repeated three times.</figcaption></figure>
INTRODUCTION

Phosphorus (P) is an essential nutrient for optimal plant growth, development, and productivity. The effect of inorganic phosphate (Pᵢ) deficiency on crop yield has become a worldwide concern in recent years due to the pressing problem of food availability (Holford, 1997; Abelson, 1999; Neset and Cordell, 2012). Accordingly, this issue has received increasing attention. To date, research has provided a comprehensive view of the adaptive strategies used by plants to cope with Pi deficiency (Rouached et al., 2010; Nussaume et al., 2011; Lin et al., 2014). Lists of genes differentially expressed in Pᵢ deficient conditions compared to control conditions have been produced, and signalling pathways involved in maintaining the homeostasis of Pᵢ have been proposed (Misson et al., 2005; Lin et al., 2008; Lin et al., 2014). However, the practical application of such knowledge is hindered by complex crosstalk linking Pi nutrition and nutrition of other essential micronutrients, e.g., elements are not assimilated, but are taken up zinc (Zn) and iron (Fe) (Huang et al., 2000; Bouain et al., 2014a; Khan et al., 2014). Such interconnections may account for the shortcomings of current agronomic models that typically focus on improving the assimilation of individual elements (Rouached et al., 2011; Ova et al., 2015). Yet despite their fundamental importance, the molecular bases, biological significance and agronomical repercussions of these interactions remain unknown.

Our current understanding of genes involved in the transport of Pᵢ into and within the plant has considerably progressed in recent decades (Poirier and Bucher, 2002; Nussaume et al., 2011; Bayle et al., 2011; Sun et al., 2012; Poirier and Jung, 2015). The first step of Pᵢ transport in plants is its acquisition at the soil-root interface. Genes belonging to the *PHT1* family, which are involved in Pi uptake and transport, have been cloned and characterized in both monocotyledons and dicotyledons (Nussaume et al., 2011; Sun et al., 2012). The increased interest in genes involved in the second rate-limiting step in Pᵢ transport within the plant, i.e. its transfer to xylem vessels, led to the identification of *PHO1* gene in Arabidopsis as a key player in this process (Hamburger et al., 2002; Secco et al., 2010). *Atpho1* mutants showed defects in Pi homeostasis in the form a marked reduction in Pᵢ transfer from the roots to the shoot, which was accompanied by low-shoot and high-root Pᵢ concentrations (Hamburger et al., 2002; Rouached et al., 2011; Poirier and Jung, 2015). The PHO1 protein feature is characterized by the presence of two distinct domains, named SPX and EXS, which may play a role in either Pᵢ transport/ sensing and sorting proteins to endomembranes respectively (Wang et al., 2004; Wege et al., 2015). Thus, the PHO1 protein does not share
Effects of a Combination of Nutrient
Availabilities on Pi, Zn, and Fe
Concentrations in Rice (Nipponbare)

The effects of a combined Pi, Zn, and/or Fe deficiency stress on the accumulation of each of these elements in rice were assessed in wild type rice plants grown for 21 days on complete solution or in solution deficient in Pi, Fe, or Zn, and combinations of two and all three deficiencies. The concentration of each element was determined both in shoot and roots tissues as indicated in material section. Under Pi deficiency, the Pi concentration was reduced by ≈1/4-fold compared to rice plants grown on complete medium composition, thereby confirming the plants suffered from Pi deficiency. A comparable effect (≈1/4-fold decrease) was observed in response to PiFe, PiZn double stresses

or PiZnFe triple deficiency stresses. Pi concentrations the shoot
increase by 1.5-fold in response Fe deficiency, but only by
≈0.5 by Zn or FeZn deficiency. The Zn concentrations in the
shoots and in the roots in wild type rice plants subjected to
the seven conditions detailed above decreased considerable in
all conditions in which Zn was omitted. Zn concentration in
the shoot was reduced by Pi deficiency stress, whereas combined
Pi and Fe deficiency had no effect. However, the individual Fe
deficiency considerable increased Zn content in both the shoot
(> threefold increase) and to a lesser extent in the roots (> 1.5-fold
increase). Like Zn, the Fe content in the shoot and roots of rice
plants subjected to nutritional deficiency was strongly reduced
when this element was omitted. Fe content increased slightly
in the shoots of Pi, Zn, or Pi/Zn deficient plants (Figure 3).

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

**FIGURE 2** | Biomass of wild typeerice plant and *Ospho1;1-1* and *Ospho1;1-2* mutant lines grown on different mineral element supply. (A) Wild type (WT), (B) *Ospho1;1-1*, and (C) *Ospho1;1-2* mutant seedlings were grown on Yoshida medium with different ionic concentrations: Complete medium (Ct) or lacking either individual element phosphate (–Pi), zinc (Zn), or iron (Fe); double elements –Pi–Zn, –Pi–Fe, –Zn–Fe, or triple elements –Pi–Fe–Zn. FW; Fresh weight is expressed as g/plant. DW, Dry weight. The experiment was repeated three times. The data are given as means ± SE (n = 9). Asterisk indicates statistical significance, *P* < 0.05. Double asterisk indicates statistical significance, *P* < 0.01.
protein homology with any other characterized Pᵢ transporters (Wang *et al.*, 2004). The Arabidopsis genome contains 11 members of the PHO1 family (Wang *et al.*, 2004). So far, only three members, AtPHO1, AtPHO1;H1 (Stefanovic *et al.*, 2007) and AtPHO1;H3 (Khan *et al.*, 2014), have been reported to be involved in Pᵢ root-to-shoot translocation. These data indicate that the role of members of the Arabidopsis PHO1 family is not restricted to Pi transport and homeostasis. For example, the *PHO1;H4* homolog is involved in the response of hypocotyls to blue light (Kang and Ni, 2006), as well as in seed size (Zhou et al., 2009) and flowering (Zhou and Ni, 2009). The role of the *PHO1* gene family has been investigated in a few crop species, including rice (*Oryza sativa* L.) (Secco *et al.*, 2010). Three rice *PHO1* genes, named *OsPHO1;1*, *OsPHO1;2* and *OsPHO1;3* have been identified. So far, only *OsPHO1;2* has been shown to play a key role in the transfer of Pi from roots to shoots and regulated by Pi deficiency (Secco *et al.*, 2010). The role of the other two rice *PHO1* genes, *OsPHO1;1* and *OsPHO1;3*, remains to be investigated.

Pi is known to interact with many micronutrients including Zn and Fe (for a review, see Bouain *et al.*, 2014a; Briat *et al.*, 2015). The effects of Pi, Zn or Fe availability on plant physiology and transcriptomic have typically been studied individually (Misson et al., 2005; Van DeMortel *et al.*, 2010; Zheng *et al.*, 2009). Although this is an important way to address the question, it can only provide a partial view of the situation. The interactions between Pi, Zn and Fe metabolisms have been demonstrated both in graminaceous and non-graminaceous plants (Zheng *et al.*, 2009; Briat *et al.*, 2015). These interactions have started being documented at a molecular level, showing that Fe or Zn deficiency modifies the expression of genes involved in Pi transport and assimilation, and vice-versa (Zheng *et al.*, 2009; Briat *et al.*, 2015). Nevertheless, the characterization of these interactions is still in its infancy. In addition to the well-established antagonistic relationship between Fe and Zn in plants (Haydon *et al.*, 2012; Shanmugam *et al.*, 2012; Briat *et al.*, 2015), the existence of an interaction between the homeostasis of Pi and Zn (Huang *et al.*, 2000; Khan *et al.*, 2014; Bouain *et al.*, 2014b; Ova *et al.*, 2015) and Pi and Fe (Zheng *et al.*, 2009) was recently reported, pointing to the existence of a complex tripartite PZnFe signalling interaction to regulate plant growth. Therefore, undertaking an integrative study to clarify the exact molecular mechanisms that coordinate the Pi, Fe and Zn deficiency pathways and thus shape the plants development during multiple PZnFe deficiency stress is of great interest.
A striking increase in Fe content was observed in the roots of plants under combined -Pi-Zn stress, which was not observed in any case of single stress. This clearly demonstrates that Pi and Zn homeostasis not only interact but also influence Fe nutrition.

## Effects of Pi, Zn, and Fe Availability on the Expression of Members of the OsPHO1 Gene Family

The previous ion concentrations analysis showed that rice plants exposed to either Fe or Zn deficiency or to FeZn simultaneous stress overaccumulated Pi in the shoot (Figure 3). Previous research work by Secco et al. (2010) suggests that members of the OsPHO1 gene family (OsPHO1;1, OsPHO1;2, and OsPHO1;3) are good candidates for the transfer of Pi from the roots to the shoot. Each of these three OsPHO1 genes has its own cis-natural antisense transcript (cis-NATs) (Secco et al., 2010), which could play a role by stabilizing the transfer of its sense gene, as demonstrated by Jabnoune et al. (2013) for OsPHO1;2. In the present study, the regulation of the sense and antisense transcripts of the OsPHO1 gene family member was analyzed in wild type rice roots subjected to the each of the above-mentioned growing ions-deficient conditions. Results in Figure 4 show that in most cases, there was an increase in sense or antisense transcript level for the three OsPHO1 gene family members in response to Pi deficiency. Previous results by Secco et al., (2010) showed the most abundantly expressed gene was OsPHO1;2 in the roots, for both sense and antisense transcripts. Our results revealed an interesting profile was observed for OsPHO1;1 gene in response to the stress caused by combined PiZn deficiency (Figure 4). It is worth noting that the high antisense:sense ratio of OsPHO1;1 observed under this condition was not observed when rice plants were treated with single stress caused by Pi or Zn deficiency. This result suggests that the expression of OsPHO1;1 depends on simultaneous PiZn deficiency signaling, and raise the questions on the involvement of the cis-natural antisense its regulation, in similar manner to OsPHO1;2 (Secco et al., 2010).

## Effects of a Combination of Pi, Zn, and Fe Deficiency Stresses on the Production of Biomass and the Concentration of Ions in the Ospho1;1 Rice Mutant

The role of OsPHO1;1 in rice is not known to date. The over-concentration of its transcripts in the combined Pi/Zn deficiency stress prompted us to assess the effect the mutation of OsPHO1;1 in shoot and root biomass when grown in the above-mentioned combinations of nutrient stresses. Two independent lines, named Ospho1;1-1 and Ospho1;1-2 were used for in this analysis (Secco et al., 2010). Figure 1 shows the phenotypes of 21 days old rice Ospho1;1 mutant lines grown on different medium compositions. The results presented in Figures 1 and 2 for the two Ospho1;1 mutant lines were comparable to the wild type in term of growth and biomass. Nevertheless, analysis of the concentration of ions (Pi, Zn, and Fe) (Figures 3A–F) in rice tissues revealed an over-accumulation of Fe (twofold)

in the roots of Ospho1;1 mutants compared to wild type plant. This difference in Fe concentration in roots was not observed in the shoots of the rice plants (Figures 3E,F). These data show that the regulation of Fe content root was altered in the Ospho1;1 mutant lines in response to the double PiZn deficiency, supporting a role for OsPHO1;1 in the regulation of Fe homeostasis.

## DISCUSSION

Compared with a growing list of documented effects of individual nutritional stress in plants, little is known about the effects of combinations of macro- and micronutrient deficiencies in plants (Mathan and Amberger, 1977; Hirsch et al., 2006; Bournier et al., 2013; Jain et al., 2013; Bouain et al., 2014a; Briat et al., 2015; Rai et al., 2015). To date research on Pi, Zn, and Fe nutrition has focused on the regulation of the homeostasis of each element individually. The effect of one element on the concentration of the other one is just starting to be documented, and point to the existence of a complex tripartite interaction between these three nutrients (Huang et al., 2000; Zheng et al., 2009; Bouain et al., 2014a; Briat et al., 2015; Ova et al., 2015). But despite its importance, our understanding of how the signaling pathways of these nutrients are “wired” into functional networks to control plant growth is still very limited.

The antagonistic relationship between Fe and Zn in plants on one hand, and between Zn or Fe starvation and Pi concentration on the other hand has been documented in many plant species (Huang et al., 2000; Zheng et al., 2009; Bouain et al., 2014b; Khan et al., 2014; Briat et al., 2015). The present study provides additional evidence for the existence and for the importance of two-by-two interactions between these ions in rice. For example, the observed effects of Pi or Zn deficiency on root biomass were only visible when the rice plants were subjected to double PiZn deficiency. Interestingly, in this condition, the reduction in root biomass occurred concomitantly with a marked increase in the concentration of Fe. The consequence of the accumulation of Fe on plant development has been described in Arabidopsis (Ward et al., 2008). These authors proposed that Fe excess in Pi deficient Arabidopsis plants, grown in agar plates, was the cause of a severe reduction in the primary root growth. Recently, it has been proposed that iron-dependent callose deposition adjusts root meristem maintenance to Pi availability (Muller et al., 2015). Consistently, simultaneous Pi and Fe deficiency (PiFe) restored primary root growth. In the present study, while the Fe deficiency had the biggest single-nutrient effects in rice plants grown hydroponically, this effect appeared to depend on the presence of Pi, or Zn, or of both elements simultaneously, in the medium (Figure 1). Indeed, stress caused by the simultaneous deficiency of two nutrients, namely PiFe or ZnFe, alleviates the effect of single Fe deficiency (Fe) on biomass, a phenomenon that was even clearer under triple stress (PiZnFe) (Figure 1). Thus, Fe deficiency cannot be considered as being solely responsible for the observed reduction in growth, but rather the interaction between the
Rice (*Oryza sativa* L.) is one of the most important cereal crops for human consumption, and feeds about half the world population (http://faostat3.fao.org/browse/Q/QC/E). The effect of individual Pi, Zn or Fe deficiencies on rice shoot and root growth is fairly well documented (e;g Wissuwa et al., 2006; Zheng *et al.*, 2009; Secco *et al.*, 2010), but the effect of their combined stresses on plant biomass has been the subject of fewer investigations. In this context, the primary aim of the present study was to broaden our understanding of the long-term effect of different combinations of stresses caused by a deficiency in these essential macro- and micronutrients in rice. Wild type plants and selected mutant lines affected in members of *OsPHO1* gene family, namely *Ospho1;1*, were included in this work. Pi, Zn and Fe ion concentrations and biomass were determined in rice plants grown in seven different combinations of growing conditions. Physiological and gene expression analyses revealed the involvement of *OsPHO1;1* in the interaction between Fe transport and PᵢZn signalling pathways in rice.

## RESULTS

### Effects of a combination of stresses caused by Pi, Zn and Fe deficiency on biomass in rice (Nipponbare).

Wild type rice plants (variety Nipponbare) were grown hydroponically in seven different growing conditions including single (-Pᵢ, -Zn and -Fe), double (-Pᵢ-Fe, -Pᵢ-Zn and -Zn-Fe), and triple (-Pᵢ-Fe-Zn) ion deficiency stresses. Figure 1A shows the phenotypes of 21-days old wild type rice plants. The most severe effect on shoot growth was observed for Fe deficiency in comparison to the effects of single Pᵢ or Zn deficiency stress. Interestingly, the Fe deficiency stress was alleviated when combined with Pi limitation (-Pᵢ-Fe or -Pᵢ-Fe-Zn) (Figure 1). This is reminiscent of data showing that the striking aspect of the inhibition of primary root elongation in Arabidopsis was recovered when the concentrations of Fe and Pi in medium were simultaneously reduced (Ward et al., 2008). This also suggests that manipulating Pᵢ availability to a rice plant could be a valuable strategy for improving a plant's ability to tolerate Fe deficiency. To gain quantitative insights into these combined stresses effect, the rice biomass was assessed by measuring the roots and shoots dry weight (Figure 2 A-B). Stress caused by Fe deficiency had the strongest negative influence on the shoot and root dry biomass, whereas Zn deficiency had a slight yet significant positive effect (Figure 2A-B). Pᵢ limitation alone had no effect on these parameters. Change in Pᵢ together with
<figure><img src="image_4.png" /></figure>

FIGURE 3 | Phosphate, Zinc, and Iron concentrations in wild type (WT) rice and *OsPHO1;1* mutant lines plant grown under individual and combinatory ions deficiency stresses. (A,C,E) Present Pi, Zn, and Fe concentrations in the shoots, respectively. (B,D,F) Present Pi, Zn, and Fe concentrations in the roots, respectively. After a 10 days of germination in 1/4 full-strength of Yoshida media, plants were grown for 3 weeks on full-strength Yoshida medium with different ionic concentrations: Complete medium (Ct) or lacking either individual element Pi, Zn, or Fe; double elements PiZn, PiFe, ZnFe, or triple elements PIFeZn. Shoots and roots were harvested separately and the Pi and Zn and Fe were extracted and their concentrations were determined as indicated in material and methods. Pi concentrations are calculated per fresh weight. Zn and Fe concentrations are calculated per dry weight. The experiment was repeated three times. The data are given as means ± SE (n = 9). Asterisk indicates statistical significance, *P* < 0.05. Double asterisk indicates statistical significance, **P** < 0.01.

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

FIGURE 4 | Transcripts abundance of sense and anti-sense of OsPHO1;1 and OsPHO1;2 in OsPHO1;3 as determined by qRT-PCR. OsActin1 was used as an internal reference. qRT quantitative real-time-PCR was performed to analyze the OsPHO1;1 and OsPHO1;2 in OsPHO1;3 mRNA levels in response to under multiple nutritional stresses: complete Yoshida medium (Ct) or lacking to either individual element phosphate (–Pi), zinc (–Zn), or iron (–Fe); double elements –Pi–Zn, –Pi–Fe, –Zn–Fe, or triple elements –Pi–Fe–Zn. All values are relative to wild-type (WT) roots. The experiment was repeated three times. The data are given as means ± SE (n = 9). Asterisk indicates statistical significance, *P* < 0.05. Double asterisk indicates statistical significance, **P** < 0.01.
changes in Fe or Zn had striking effects. Simultaneous Pᵢ and Fe stress (PᵢFe) alleviated the severe stress caused by the Fe deficiency alone, while stress caused by a simultaneous deficiency in Pᵢ and Zn had a negative effect. Double (FeZn) or triple (PᵢFeZn) deficiencies had a negative effect on biomass similar to the effect of Fe deficiency alone. Shoot to root dry weight ratio showed that except for -Pᵢ stress, all other nutrient stress conditions tested had a negative effect on biomass (Figure 2C).

Taken together, the results indicate that supply of one mineral element has a sometimes strong effect on the plant response to supply of other elements. For example, as shown in Figure 1, Fe deficient plants were smaller and chlorotic, but Fe deficient plants also lacking P and Zn looked healthier in comparison. These data underline thus the importance of the coordination between Pᵢ, Zn and Fe homeostasis in the control of rice growth capacity and shoot and root biomass.

## Effects of a combination of nutrient availabilities on Pᵢ, Zn and Fe concentrations in rice (Nipponbare).

The effects of a combined Pᵢ, Zn and/or Fe deficiency stress on the accumulation of each of these elements in rice were assessed in wild type rice plants grown for 21 days on complete solution or in solution deficient in Pᵢ, Fe, or Zn, and combinations of two and all three deficiencies. The concentration of each element was determined both in shoot and roots tissues as indicated in material section. Under Pᵢ deficiency, the Pᵢ concentration was reduced by ≈ 1/4 fold compared to rice plants grown on complete medium composition, thereby confirming the plants suffered from Pᵢ deficiency. A comparable effect (≈ 1/4 fold decrease) was observed in response to PᵢFe, PᵢZn double stresses or PᵢZnFe triple deficiency stresses. Pᵢ concentrations the shoot increase by 1.5 fold in response Fe deficiency, but only by ≈ 0.5 by Zn or FeZn deficiency. The Zn concentrations in the shoots and in the roots in wild type rice plants subjected to the seven conditions detailed above decreased considerably in all conditions in which Zn was omitted. Zn concentration in the shoot was reduced by Pᵢ deficiency stress, whereas combined Pᵢ and Fe deficiency had no effect. However, the individual Fe deficiency considerably increased Zn content in both the shoot (>3 fold increase) and to a lesser extent in the roots (>1.5 fold increase). Like Zn, the Fe content in the shoot and roots of rice plants subjected to nutritional deficiency was strongly reduced when this element was omitted. Fe content increased slightly in the shoots of Pᵢ, Zn or PᵢZn deficient plants (Figure 3). A striking increase in Fe content was observed in the roots of
homeostasis of the three nutrients PiZnFe. Together, the data presented in this study shed light on the existence of a tripartite interaction between these macro-and micronutrients in plants and provide the basis for investigating how plants integrate multiple nutrient deficiencies signaling in complex plant developmental programs.

Plants must adapt to a fluctuating environment, and changes in gene expression are a key part of this adaptation. The identification of key genes involved in the co-regulation of micro- and macronutrient homeostasis has become the focus of major interest (Zheng et al., 2009; Rouached et al., 2010b; Kellermeier et al., 2014; Khan et al., 2014). Recently, a member of the AtPHO1 gene family, *AtPHO1;H3*, was identified as an important player in the crosstalk between Zn deficiency signaling and the regulation of Pi homeostasis in *Arabidopsis* (Khan et al., 2014). In the present work, the regulation of the expression of a member of *OsPHO1* gene family, *OsPHO1;1*, also appears to depend on simultaneous PiZn deficiency signaling, and mutations of *Ospho1;1* lead to over-accumulation of Fe mainly in roots compared to wild type plants. An attractive hypothesis to explain the results of this work would be that the effect of Ospho1.1 mutation interfere with the transport of Fe in roots under simultaneous -Pi-Zn deficiency. This hypothesis is supported by observations in yeast (*S. cerevisiae*) where the high-affinity Pi transporter PHO84 can regulate or transport metal ions such as “Metal-HPO4” (Luk et al., 2003a,b). Plants like *Arabidopsis* and rice possess various Pi transporters with significant similarity to PHO84, such as members of the *PHO1* gene family (e.g., *OsPHO1;1*). There is currently no evidence for either Fe or Fe/Pi transport by a plant Pi transporter, though this remains an intriguing possibility. The regulation of *OsPHO1;1* by its *cis-NATPHO1;1* can be also proposed for examination in light of the following arguments. First, the natural antisense of *OsPHO1;1* (*cis-NATPHO1;1*) was up-regulated by the combined -Pi-Zn deficiency stress (Figure 4). Second, a recent study by Jabnoune et al. (2013) reported an unexpected role for *cis-NATPHO1;2* in promoting *OsPHO1;2* translation and affecting Pi homeostasis and plant fitness. Considering this recent discovery (Jabnoune et al., 2013) and the results of the present study, one can hypothesize that *OsPHO1;1* could be a target in the -Pi-Zn deficiency signaling pathway and involved in the regulation of the Fe homeostasis.

## CONCLUSION

This study is the first documentation of the effect multiple nutritional stresses on rice growth and biomass, which is illustrated through the interaction between the homeostasis of essential macro(Pi) - and micronutrients (Zn and Fe). This obtained data provides evidence for the existence of a genetic basis for this tripartite nutrient PiZnFe interaction in plant. A biological function for the OsPHO1;1 in this tripartite interaction is proposed, and particularly in the regulation of Fe transport through integration of Pi and Zn deficiency signaling.

# MATERIALS AND METHODS

## Plant Growth Conditions

The rice (*Oryza sativa* cv Nipponbare), wild type (WT) and mutant lines: *Ospho1; 1-1* and *Ospho1;1-2*, was grown hydroponically under non-sterile conditions. Seeds were soaked in deionized water over night in dark then transferred to 1/4 full-strength Yoshida media for 10 days (Yoshida et al., 1976). Seedlings were transferred to the following modified Yoshida nutrient solution: 1.43 mM NH₄NO₃; 1.64 mM MgSO₄; 0.75 mM CaCl₂; 0.51 mM K₂SO₄; 0.33 mM NaH₂PO₄; 0.02 mM H₃BO₃; 0.01 mM MnCl₂; 0.04 mM Fe-NaEDTA; 2.5 µM ZnSO₄; 0.16 µM CuSO₄; 0.08 µM (NH₄)₆Mo₇O₂₄; 2.5 µM MES buffer, adjust to pH 5.5. The nutrient solution was renewed every 5 days. Plants were grown in a growth chamber under the following environmental conditions: light/dark cycle of 14/10 h, temperature of 28/25°C, and RH of 80%. For the control condition, plants were kept in nutrient solution with the above-mentioned composition. The treatments were applied to 10 day-old plants. The concentration of individual nutrients in the complete media were 0.3 mM P, 0.025 mM Zn, and 0.04 mM Fe and zero in the deficient media for -Pi, -Zn, -Fe for single, double or triple ion deficiency combinations. The daytime temperature was 28°C and the night temperature 25°C with a 14/10 h light/dark regime. Plants were harvested at 3 weeks after transferred to treatments. Roots and shoots were separated, and immediately frozen in liquid N, and stored at -80°C.

## Real-Time Quantitative Reverse-Transcription PCR

DNA-free total RNA extraction from frozen root tissues was performed using Plant RNeasy extraction kit and RNAse-free DNaseI (SIGMA-ALDRCH, St Louis, MO, USA). Two micrograms of total RNA were reverse transcribed using ThermoTM script RT (Invitrogen) according to the manufacturer's protocol. Complementary DNA (cDNA) was used for real-time reverse-transcription PCR, which was performed with (LightCycler® 480; Roche Diagnostics). PCR reactions containing 12.5 µL of the LightCycler® 480 SYBR Green I Master mix (Roch, IN, USA), each of the forward and reverse primers, and 5 µL of a 1:50 cDNA dilution in a final volume of 25 µL were considered for gene expression analysis. All PCR reactions were performed in triplicate. In addition to the three *OsPHO1* genes, *OsActin1* gene was considered for the standardization of real-time PCR data. These genes were chosen based on previous analyses (Secco et al., 2010). Quantification of the relative transcripts levels was performed using the comparative CT method (Livak and Schmittgen, 2001; Rouached et al., 2008).

## Zinc and Iron Measurements

Shoots and roots material were dried at 72°C for at least 7 days. The dried tissues were ground and subjected (~20 mg) to acid hydrolysis. Zn and Fe concentrations in the samples were determined using the Agilent 4200 Microwave Plasma-Atomic Emission Spectrometer (MP-AES 4200, Agilent, USA).
plants under combined -Pi-Zn stress, which was not observed in any case of single stress. This clearly demonstrates that Pi and Zn homeostasis not only interact but also influence Fe nutrition.

**Effects of Pi, Zn and Fe availability on the expression of members of the OsPHO1 gene family.**

The previous ion concentrations analysis showed that rice plants exposed to either Fe or Zn deficiency or to FeZn simultaneous stress overaccumulated $P_i$ in the shoot (Figure 3). Previous research work by Secco et al. (2010) suggests that members of the OsPHO1 gene family (*OsPHO1;1*, *OsPHO1;2* and *OsPHO1;3*) are good candidates for the transfer of Pi from the roots to the shoot. Each of these three *OsPHO1* genes has its own cis-natural antisense transcript (cis-NATs) (Secco et al., 2010), which could play a role by stabilizing the transfer of its sense gene, as demonstrated by Jebnoune et al. (2013) for *OsPHO1;2*. In the present study, the regulation of the sense and antisense transcripts of the *OsPHO1* gene family member was analyzed in wild type rice roots subjected to the each of the above-mentioned growing ions-deficient conditions. Results in Figure 4 show that in most cases, there was an increase in sense or antisense transcript level for the three *OsPHO1* gene family members in response to Pi deficiency. Previous results by Secco et al., (2010) showed the most abundantly expressed gene was OsPHO1;2 in the roots, for both sense and antisense transcripts. Our results revealed an interesting profile was observed for OsPHO1;1 gene in response to the stress caused by combined PiZn deficiency (Fig 4A). It is worth noting that the high antisense:sense ratio of *OsPHO1;1* observed under this condition was not observed when rice plants were treated with single stress caused by Pi or Zn deficiency. This result suggests that the expression of *OsPHO1;1* depends on simultaneous PiZn deficiency signalling, and raise the questions on the involvement of the cis-natural antisense its regulation, in similar manner to OsPHO1;2 (Secco et al., 2010).

**Effects of a combination of Pi, Zn and Fe deficiency stresses on the production of biomass and the concentration of ions in the *Ospho1;1* rice mutant.**

The role of *OsPHO1;1* in rice is not known to date. The over-concentration of its transcripts in the combined Pi/Zn deficiency stress prompted us to assess the effect the mutation of *OsPHO1;1* in shoot and root biomass when grown in the above-mentioned combinations of nutrient stresses. Two independent lines, named *Ospho1;1-1* and *Ospho1;1-2* were used for in this analysis (Secco et al., 2010). Figure 1B-C shows the phenotypes of 21 days old rice
# Pi Measurements

For Pi extractions, weighed fresh shoots and roots were collected separately and incubated material in water for 3 h at 70°C. Pi measurements were performed as described by Rouached et al. (2010a). The quantification of Pi was performed by the molybdate assay according to Ames (1966).

# Statistical analysis

We carried out statistical analysis using The Analysis Toolpak add-in program for Excel for Macintosh (Microsoft Incorporated, USA). For all the t test analyses the difference is considered statistically significant with a probability of *p* < 0.05.

# AUTHOR CONTRIBUTIONS

HR conceived and designed the experiment. CS, NB, and MK performing collected the experimental data. CS and NB equally contributed to this work (co-first authors). CP, PD, and HR analyzed data. HR wrote the manuscript. CP and PD assisted with the interpretation of the results and provided editorial support for

the manuscript. All authors have read, edited, and approved the current version of the manuscript.

# FUNDING

This work was funded by the Institut National de la Recherche Agronomique (INRA) to HR, a Franco-Thailand Fellowship to CS.

# ACKNOWLEDGMENTS

The authors are grateful to Pr. Yves Poirier (Lausanne University, Switzerland) for providing seeds of the OsPHO1 family mutants. Authors wish to thank Dr. Emmanuel Guiderdoni (CIRAD, Montpellier, France) for providing plant material for this study and Dr. Sue Rhee (Carneigie Institution for Science, Stanford, CA, USA) for helpful discussion on this research topic. We are grateful to Pr. Pierre Berthomieu for providing the necessary infrastructure to finalize this work in his group.

# REFERENCES

* Abelson, P. H. (1999). A potential phosphate crisis. *Science* 283, 2015. doi: 10.1126/science.283.5410.2015

* Ames, B. N. (1966). Assay of inorganic phosphate, total phosphate and phosphatases. *Methods Enzymol.* 8, 115-118. doi: 10.1016/0076-6879(66)08014-5

* Bayle, V., Arrighi, J. F., Creff, A., Nespoulous, C., Vialaret, J., Rossignol, M., et al. (2011). *Arabidopsis thaliana* high-affinity phosphate transporters exhibit multiple levels of posttranslational regulation. *Plant Cell* 23, 1523-1535. doi: 10.1105/tpc.110.081067

* Bouain, N., Kisko, M., Rouached, A., Dauzat, M., Lacombe, B., Belgaroui, N., et al. (2014a). Phosphate/zinc interaction analysis in two lettuce varieties reveals contrasting effects on biomass, photosynthesis, and dynamics of Pi transport. *Biomed. Res. Int.* 2014:548254. doi: 10.1155/2014/548254

* Bouain, N., Shahzad, Z., Rouached, A., Khan, G. A., Berthomieu, P., Abdelly, C., et al. (2014b). Phosphate and zinc transport and signalling in plants: toward a better understanding of their homeostasis interaction. *J. Exp. Bot.* 65, 5725-5741. doi: 10.1093/jxb/eru314

* Bournier, M., Tissot, N., Mari, S., Boucherez, J., Lacombe, E., Briat, J. F., et al. (2013). *Arabidopsis* ferritin 1 (AtFer1) gene regulation by the phosphate starvation response 1 (AtPHR1) transcription factor reveals a direct molecular link between iron and phosphate homeostasis. *J. Biol. Chem.* 288, 22670-22680. doi: 10.1074/jbc.M113.482281

* Briat, J., Rouached, H., Tissot, N., Gaymard, F., and Dubos, C. (2015). Integration of P, S, Fe and Zn nutrition signals in *Arabidopsis thaliana*: potential involvement of PHOSPHATE STARVATION RESPONSE 1 (PHR1). *Front. Plant Sci.* 6:290. doi: 10.3389/fpls.2015.00290

* Hamburger, D., Rezzonico, E., Macdonald-Comber Petetot, J., Somerville, C., and Poirier, Y. (2002). Identification and characterization of the *Arabidopsis* PHO1 gene involved in phosphate loading to the xylem. *Plant Cell* 14, 889-902. doi: 10.1105/tpc.000745

* Haydon, M. J., Kawachi, M., Wirtz, M., Hillmer, S., Hell, R., and Kramer, U. (2012). Vacuolar nicotianamine has critical and distinct roles under iron deficiency and for zinc sequestration in *Arabidopsis*. *Plant Cell* 24, 724-737. doi: 10.1105/tpc.111.095042

* Hirsch, J., Marin, E., Floriani, M., Chiarenza, S., Richaud, P., Nussaume, L., et al. (2006). Phosphate deficiency promotes modification of iron distribution in *Arabidopsis* plants. *Biochimie* 88, 1767-1771. doi: 10.1016/j.biochi.2006.05.007

* Holford, I. (1997). Soil phosphorus: its measurement, and its uptake by plants. *Aust. J. Soil Res.* 35, 227-240. doi: 10.1071/S96047

* Huang, C., Barker, S. J., Langridge, P., Smith, F. W., and Graham, R. D. (2000). Zinc deficiency up-regulates expression of high-affinity phosphate transporter genes in both phosphate-sufficient and -deficient barley roots. *Plant Physiol.* 124, 415-422. doi: 10.1104/pp.124.1.415

* Jabnoune, M., Secco, D., Lecampion, C., Robaglia, C., Shu, Q., and Poirier, Y. (2013). A rice cis-natural antisense RNA acts as a translational enhancer for its cognate mRNA and contributes to phosphate homeostasis and plant fitness. *Plant Cell* 25, 4166-4182. doi: 10.1105/tpc.113.116251

* Jain, A., Sinilal, B., Dhandapani, G., Meagher, R. B., and Sahi, S. V. (2013). Effects of deficiency and excess of zinc on morphophysiological traits and spatiotemporal regulation of zinc-responsive genes reveal incidence of cross talk between micro- and macronutrients. *Environ. Sci. Technol.* 47, 5327-5335. doi: 10.1021/es400113y

* Kang, X., and Ni, M. (2006). *Arabidopsis* SHORT HYPOCOTYL UNDER BLUE1 contains SPX and EXS domains and acts in cryptochrome signaling. *Plant Cell* 18, 921-934. doi: 10.1105/tpc.105.037879

* Kellermeier, F., Armengaud, P., Seditas, T. J., Danku, J., Salt, D. E., and Amtmann, A. (2014). Analysis of the root system architecture of *Arabidopsis* provides a quantitative readout of crosstalk between nutritional signals. *Plant Cell* 26, 1480-1496. doi: 10.1105/tpc.113.122101

* Khan, G. A., Bouraine, S., Wege, S., Li, Y., De Carbonnel, M., Berthomieu, P., et al. (2014). Coordination between zinc and phosphate homeostasis involves the transcription factor PHR1, the phosphate exporter PHO1, and its homologue PHO1;H3 in *Arabidopsis*. *J. Exp. Bot.* 65, 871-884. doi: 10.1093/jxb/ert444

* Lin, S. I., Chiang, S. F., Lin, W. Y., Chen, J. W., Tseng, C. Y., Wu, P. C., et al. (2008). Regulatory network of microRNA399 and PHO2 by systemic signaling. *Plant Physiol.* 147, 732-746. doi: 10.1104/pp.108.116269

* Lin, W. Y., Huang, T. K., Leong, S. J., and Chiou, T. J. (2014). Long-distance call from phosphate: systemic regulation of phosphate starvation responses. *J. Exp. Bot.* 65, 1817-1827. doi: 10.1093/jxb/ert431

* Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. *Methods* 25, 402-408. doi: 10.1006/meth.2001.1262

* Luk, E., Carroll, M., Baker, M., and Culotta, V. C. (2003a). Manganese activation of superoxide dismutase 2 in *Saccharomyces cerevisiae* requires MTM1, a
*Ospho1;1* mutant lines grown on different medium compositions. The results presented in Figures 1 and 2 for the two *Ospho1;1* mutant lines were comparable to the wild type in term of growth and biomass. Nevertheless, analysis of the concentration of ions (Pi, Zn and Fe) (Figure 3 A-F) in rice tissues revealed an over-accumulation of Fe (2-fold) in the roots of *Ospho1;1* mutants compared to wild type plant. This difference in Fe concentration in roots was not observed in the shoots of the rice plants (Figure 3 E-F). These data show that the regulation of Fe content root was altered in the *Ospho1;1* mutant lines in response to the double PiZn deficiency, supporting a role for *OsPHO1;1* in the regulation of Fe homeostasis.
member of the mitochondrial carrier family. *Proc. Natl. Acad. Sci. U.S.A.* 100, 10353–10357. doi: 10.1073/pnas.1632471100

Luk, E., Jensen, L. T., and Culotta, V. C. (2003b). The many highways for intracellular trafficking of metals. *J. Biol. Inorg. Chem.* 8, 803–809. doi: 10.1007/s00775-003-0482-3

Mathan, K. K., and Amberger, A. (1977). Influence of iron on the uptake of phosphorous by maize. *Plant Soil* 48, 413–422. doi: 10.1007/BF00010097

Misson, J., Raghothama, K. G., Jain, A., Jouhet, J., Block, M. A., Bligny, R., et al. (2005). A genome-wide transcriptional analysis using *Arabidopsis thaliana* Affymetrix gene chips determined plant responses to phosphate deprivation. *Proc. Natl. Acad. Sci. U.S.A.* 102, 11934–11939. doi: 10.1073/pnas.0505266102

Muller, J., Toev, T., Heisters, M., Teller, J., Moore, K. L., Hause, G., et al. (2015). Iron-dependent callose deposition adjusts root meristem maintenance to phosphate availability. *Dev. Cell* 33, 216–230. doi: 10.1016/j.devcel.2015.02.007

Neset, T. S., and Cordell, D. (2012). Global phosphorus scarcity: identifying synergies for a sustainable future. *J. Sci. Food Agric.* 92, 2–6. doi: 10.1002/jsfa.4650

Nussaume, L., Kanno, S., Javot, H., Marin, E., Pochon, N., Ayadi, A., et al. (2011). Phosphate import in plants: focus on the PHT1 transporters. *Front. Plant Sci.* 2:83. doi: 10.3389/fpls.2011.00083

Ova, E. A., Kutman, U. B., Ozturk, L., and Cakmak, I. (2015). High phosphorus supply reduced zinc concentration of wheat in native soil but not in autoclaved soil or nutrient solution. *Plant Soil* 393, 147–162. doi: 10.1007/s11104-015-2483-8

Poirier, Y., and Bucher, M. (2002). Phosphate transport and homeostasis in *Arabidopsis*. *Arabidopsis Book* 1:e0024. doi: 10.1199/tab.0024

Poirier, Y., and Jung, J. Y. (2015). Phosphate transporters. *Annu. Plant Rev.* 48, 125–159.

Rai, V., Sanagala, R., Sinilal, B., Yadav, S., Sarkar, A. K., Dantu, P. K., et al. (2015). Iron availability affects phosphate deficiency-mediated responses, and evidence of cross-talk with auxin and zinc in *Arabidopsis*. *Plant Cell Physiol.* 56, 1107–1123. doi: 10.1093/pcp/pcv035

Rouached, H., Arpat, A. B., and Poirier, Y. (2010a). Regulation of phosphate starvation responses in plants: signaling players and cross-talks. *Mol. Plant* 3, 288–299. doi: 10.1093/mp/ssp120

Rouached, H., Secco, D., and Arpat, B. A. (2010b). Regulation of ion homeostasis in plants: current approaches and future challenges. *Plant Signal. Behav.* 5, 501–502. doi: 10.4161/psb.11027

Rouached, H., Secco, D., Arpat, B., and Poirier, Y. (2011a). The transcription factor PHR1 plays a key role in the regulation of sulfate shoot-to-root flux upon phosphate starvation in *Arabidopsis*. *BMC Plant Biol.* 11:19. doi: 10.1186/1471-2229-11-19

Rouached, H., Stefanovic, A., Secco, D., Bulak Arpat, A., Gout, E., Bligny, R., et al. (2011b). Uncoupling phosphate deficiency from its major effects on growth and transcriptome via PHO1 expression in *Arabidopsis*. *Plant J.* 65, 557–570. doi: 10.1111/j.1365-313X.2010.04442.x

Rouached, H., Wirtz, M., Alary, R., Hell, R., Arpat, A. B., Davidian, J. C., et al. (2008). Differential regulation of the expression of two high-affinity sulfate transporters, SULTR1.1 and SULTR1.2, in *Arabidopsis*. *Plant Physiol.* 147, 897–911. doi: 10.1104/pp.108.118612

Secco, D., Baumann, A., and Poirier, Y. (2010). Characterization of the rice PHO1 gene family reveals a key role for OsPHO1;2 in phosphate homeostasis and the evolution of a distinct clade in dicotyledons. *Plant Physiol.* 152, 1693–1704. doi: 10.1104/pp.109.149872

Shanmugam, V., Tsednee, M., and Yeh, K. C. (2012). ZINC TOLERANCE INDUCED BY IRON 1 reveals the importance of glutathione in the cross-homeostasis between zinc and iron in *Arabidopsis thaliana*. *Plant J.* 69, 1006–1017. doi: 10.1111/j.1365-313X.2011.04850.x

Stefanovic, A., Ribot, C., Rouached, H., Wang, Y., Chong, J., Belbahri, L., et al. (2007). Members of the PHO1 gene family show limited functional redundancy in phosphate transfer to the shoot, and are regulated by phosphate deficiency via distinct pathways. *Plant J.* 50, 982–994. doi: 10.1111/j.1365-313X.2007.03108.x

Sun, S., Gu, M., Cao, Y., Huang, X., Zhang, X., Ai, P., et al. (2012). A constitutive expressed phosphate transporter, OsPht1;1, modulates phosphate uptake and translocation in phosphate-replete rice. *Plant Physiol.* 159, 1571–1581. doi: 10.1104/pp.112.196345

Thibaud, M. C., Arrighi, J. F., Bayle, V., Chiarenza, S., Creff, A., Bustos, R., et al. (2010). Dissection of local and systemic transcriptional responses to phosphate starvation in *Arabidopsis*. *Plant J.* 64, 775–789.

Van De Mortel, J. E., Almar Villanueva, L., Schat, H., Kwekkeboom, J., Coughlan, S., Moerland, P. D., et al. (2006). Large expression differences in genes for iron and zinc homeostasis, stress response, and lignin biosynthesis distinguish roots of *Arabidopsis thaliana* and the related metal hyperaccumulator *Thlaspi caerulescens*. *Plant Physiol.* 142, 1127–1147. doi: 10.1104/pp.106.082073

Wang, Y., Ribot, C., Rezzonico, E., and Poirier, Y. (2004). Structure and expression profile of the *Arabidopsis* PHO1 gene family indicates a broad role in inorganic phosphate homeostasis. *Plant Physiol.* 135, 400–411. doi: 10.1104/pp.103.037945

Ward, J. T., Lahner, B., Yakubova, E., Salt, D. E., and Raghothama, K. G. (2008). The effect of iron on the primary root elongation of *Arabidopsis* during phosphate deficiency. *Plant Physiol.* 147, 1181–1191. doi: 10.1104/pp.108.118562

Wege, S., Khan, G. A., Jung, J. Y., Vogiatzaki, E., Pradervand, S., Aller, I., et al. (2016). The EXS domain of PHO1 participates in the response of shoots to phosphate deficiency via a root-to-shoot signal. *Plant Physiol.* 170, 385–400. doi: 10.1104/pp.15.00975

Wissuwa, M., Ismail, A. M., and Yanagihara, S. (2006). Effects of zinc deficiency on rice growth and genetic factors contributing to tolerance. *Plant Physiol.* 142, 731–741. doi: 10.1104/pp.106.085225

Yoshida, S., Foorno, D. A., Cock, J. H., and Gomez, K. A. (1976). *Laboratory Manual for Physiological Studies of Rice*, 3rd Edn. Los Baños: International Rice Research Institute.

Zheng, L., Huang, F., Narsai, R., Wu, J., Giraud, E., He, F., et al. (2009). Physiological and transcriptome analysis of iron and phosphorus interaction in rice seedlings. *Plant Physiol.* 151, 262–274. doi: 10.1104/pp.109.141051

Zhou, Y., and Ni, M. (2009). SHB1 plays dual roles in photoperiodic and autonomous flowering. *Dev. Biol.* 331, 50–57. doi: 10.1016/j.ydbio.2009.04.023

Zhou, Y., Zhang, X., Kang, X., Zhao, X., Zhang, X., and Ni, M. (2009). SHORT HYPOCOTYL UNDER BLUE1 associates with MINISEED3 and HAIKU2 promoters in vivo to regulate *Arabidopsis* seed development. *Plant Cell* 21, 106–117. doi: 10.1105/tpc.108.064972

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2016 Saenchai, Bouain, Kisko, Prom-u-thai, Doumas and Rouached. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
DISCUSSION

Compared with a growing list of documented effects of individual nutritional stress in plants,
little is known about the effects of combinations of macro- and micronutrient deficiencies in
plants (Mathan and Amberger, 1977; Hirsch et al., 2006; Bournier et al., 2013; Jain et al.,
2013; Bouain et al., 2014a; Briat et al., 2015; Rai et al., 2015). To date research on Pᵢ, Zn and
Fe nutrition has focused on the regulation of the homeostasis of each element individually.
The effect of one element on the concentration of the other one is just starting to be
documented, and point to the existence of a complex tripartite interaction between these three
nutrients (Huang et al., 2000; Zheng et al., 2009; Bouain et al., 2014; Briat et al., 2015; Ova
et al., 2015). But despite its importance, our understanding of how the signalling pathways of
these nutrients are “wired” into functional networks to control plant growth is still very
limited.

The antagonistic relationship between Fe and Zn in plants on one hand, and between Zn or Fe
starvation and Pᵢ concentration on the other hand has been documented in many plant species
(Huang et al., 2000; Zheng et al., 2009; Bouain et al., 2014; Khan et al., 2014; Briat et al.,
2015). The present study provides additional evidence for the existence and for the
importance of two-by-two interactions between these ions in rice. For example, the observed
effects of Pᵢ or Zn deficiency on root biomass were only visible when the rice plants were
subjected to double PᵢZn deficiency. Interestingly, in this condition, the reduction in root
biomass occurred concomitantly with a marked increase in the concentration of Fe. The
consequence of the accumulation of Fe on plant development has been described in
Arabidopsis (Ward et al., 2008). These authors proposed that Fe excess in Pᵢ deficient
Arabidopsis plants, grown in agar plates, was the cause of a severe reduction in the primary
root growth. Recently, it has been proposed that iron-dependent callose deposition adjusts
root meristem maintenance to Pᵢ availability (Muller et al., 2015). Consistently, simultaneous
Pᵢ and Fe deficiency (PᵢFe) restored primary root growth. In the present study, while the Fe
deficiency had the biggest single-nutrient effects in rice plants grown hydroponically, this
effect appeared to depend on the presence of Pᵢ, or Zn, or of both elements simultaneously, in
the medium (Figure 1). Indeed, stress caused by the simultaneous deficiency of two nutrients,
namely PᵢFe or ZnFe, alleviates the effect of single Fe deficiency (Fe) on biomass, a
phenomenon that was even clearer under triple stress (PᵢZnFe) (Figure 1). Thus, Fe deficiency
cannot be considered as being solely responsible for the observed reduction in growth, but
rather the interaction between the homeostasis of the three nutrients PᵢZnFe. Together, the
data presented in this study shed light on the existence of a tripartite interaction between these
macro-and micronutrients in plants and provide the basis for investigating how plants
integrate multiple nutrient deficiencies signalling in complex plant developmental programs.

Plants must adapt to a fluctuating environment, and changes in gene expression are a key part
of this adaptation. The identification of key genes involved in the co-regulation of micro- and
macronutrient homeostasis has become the focus of major interest (Rouached et al., 2009;
Zheng et al., 2009; Khan et al., 2014; Kellermeier et al., 2014). Recently, a member of the
AtPHO1 gene family, AtPHO1;H3, was identified as an important player in the crosstalk
between Zn deficiency signalling and the regulation of Pi homeostasis in Arabidopsis (Khan
et al., 2014). In the present work, the regulation of the expression of a member of OsPHO1
gene family, OsPHO1;1, also appears to depend on simultaneous PᵢZn deficiency signalling,
and mutations of Ospho1;1 lead to over-accumulation of Fe mainly in roots compared to wild
type plants. An attractive hypothesis to explain the results of this work would be that the
effect of Ospho1.1 mutation interfere with the transport of Fe in roots under simultaneous -Pi-
Zn deficiency". This hypothesis is supported by observations in yeast (S. cerevisiae) where
the high-affinity Pi transporter PHO84 can regulate or transport metal ions such as "Metal-
HPO4" (Luk et al., 2003a,b). Plants like Arabidopsis and rice possess various Pi transporters
with significant similarity to PHO84, such as members of the PHO1 gene family (e.g.
OsPHO1;1). There is currently no evidence for either Fe or Fe/Pi transport by a plant Pi
transporter, though this remains an intriguing possibility. The regulation of OsPHO1;1 by its
cis-NATPHO1;1 can be also proposed for examination in light of the following arguments.
First, the natural antisense of OsPHO1;1 (cis-NATPHO1;1) was up-regulated by the
combined -Pi-Zn deficiency stress (Figure 4 A). Second, a recent study by Jabnoune et al.
(2013) reported an unexpected role for cis-NATPHO1;2 in promoting OsPHO1;2 translation
and affecting Pi homeostasis and plant fitness. Considering this recent discovery (Jabnoune et
al., 2013) and the results of the present study, one can hypothesize that OsPHO1;1 could be a
target in the -Pi-Zn deficiency signalling pathway and involved in the regulation of the Fe
homeostasis.

In conclusion, this study is the first documentation of the effect multiple nutritional stresses
on rice growth and biomass, which is illustrated through the interaction between the
homeostasis of essential macro(Pᵢ) - and micronutrients (Zn and Fe). This obtained data
provides evidence for the existence of a genetic basis for this tripartite nutrient PiZnFe
interaction in plant. A biological function for the OsPHO1;1 in this tripartite interaction is proposed, and particularly in the regulation of Fe transport through integration of Pᵢ and Zn deficiency signalling.

## MATERIALS AND METHODS:

### Plant Growth Conditions

The rice (Oryza sativa cv Nipponbare), wild type (WT) and mutant lines: Ospho1; 1-1 and Ospho1;1-2, was grown hydroponically under non-sterile conditions. Seeds were soaked in deionized water over night in dark then transferred to ¼ full-strength Yoshida media for 10 days (Yoshida, 1976). Seedlings were transferred to the following modified Yoshida nutrient solution: 1.43 mM NH₄NO₃; 1.64 mM MgSO₄; 0.75 mM CaCl₂; 0.51 mM K₂SO₄; 0.33 mM NaH₂PO₄; 0.02 mM H₃BO₃; 0.01 mM MnCl₂; 0.04 mM Fe-NaEDTA; 2.5 µM ZnSO₄; 0.16 µM CuSO₄; 0.08 µM (NH₄)₆Mo₇O₂₄; 2.5 µM MES buffer, adjust to pH 5.5 The nutrient solution was renewed every 5 days. Plants were grown in a growth chamber under the following environmental conditions: light/dark cycle of 14/10 h, temperature of 28/25 °C, and RH of 80 %. For the control condition, plants were kept in nutrient solution with the above-mentioned composition. The treatments were applied to 10 day-old plants. The concentration of individual nutrients in the complete media were 0.3 mM P, 0.025 mM Zn, and 0.04 mM Fe and zero in the deficient media for -Pᵢ, -Zn, -Fe for single, double or triple ion deficiency combinations. The daytime temperature was 28 °C and the night temperature 25 °C with a 14/10 h light/dark regime. Plants were harvested at 3 weeks after transferred to treatments. Roots and shoots were separated, and immediately frozen in liquid N, and stored at -80 °C.

### Real-Time Quantitative Reverse-Transcription PCR

DNA-free total RNA extraction from frozen root tissues was performed using Plant RNeasy extraction kit and RNAse-free DNAsel (SIGMA-ALDRCH, St Louis, USA). Two micrograms of total RNA were reverse transcribed using Thermo™ script RT (Invitrogen) according to the manufacturer's protocol. Complementary DNA (cDNA) was used for real-time reverse-transcription PCR, which was performed with (LightCycler® 480; Roche Diagnostics). PCR reactions containing 12.5 µL of the LightCycler® 480 SYBR Green I Master mix (Roch, IN, USA), each of the forward and reverse primers, and 5 µL of a 1:50 cDNA dilution in a final volume of 25 µL were considered for gene expression analysis. All PCR reactions were performed in triplicate. In addition to the three OsPHO1 genes, OsActin1 gene was considered for the standardization of real-time PCR data. These genes were chosen
based on previous analyses (Lee *et al.*, 2009; Secco *et al.*, 2010; Wu *et al.*, 2011).

Quantification of the relative transcripts levels was performed using the comparative CT method (Livak and Schmittgen, 2001; Rouached *et al.*, 2008).

## Zinc and Iron Measurements

Shoots and roots material were dried at 72 °C for at least seven days. The dried tissues were ground and subjected (approximately 20 mg) to acid hydrolysis. Zn and Fe concentrations in the samples were determined using the Agilent 4200 Microwave Plasma-Atomic Emission Spectrometer (MP-AES 4200, Agilent, USA).

## Pi Measurements

For Pi extractions, weighed fresh shoots and roots were collected separately and incubated material in water for 3 h at 70 °C. Pi measurements were performed as described by Rouached *et al.* (2010). The quantification of Pi was performed by the molybdate assay according to Ames (1966).

## Statistical analysis

We carried out statistical analysis using The Analysis Toolpak add-in program for Excel for Macintosh (Microsoft Incorporated, USA). For all the t test analyses the difference is considered statistically significant with a probability of p < 0.05
## ACKNOWLEDGMENT:

The authors are grateful to Pr. Yves Poirier (Lausanne University, Switzerland) for providing seeds of the OsPHO1 family mutants. Authors wish to thank Dr. Emmanuel Guiderdoni (CIRAD, Montpellier, France) for providing plant material for this study and Dr. Sue Rhee (Carneigie Institution for Science, Stanford, USA) for helpful discussion on this research topic. We are grateful to Pr. Pierre Berthomieu for providing the necessary infrastructure to finalize this work in his group.

Provisional
REFERENCE :

Abelson, P.H. (1999). A potential phosphate crisis. Science 283, 201-205.

Ames, B.N. (1966). Ames, B.N. (1966) Assay of inorganic phosphate, total phosphate and phosphatases. Methods Enzymol 8, 115-118.

Bayle, V., Arrighi, J.F., Creff, A., Nespoulous, C., Vialaret, J., Rossignol, M., Gonzalez, E., Paz-Ares, J., and Nussaume, L. (2011). Arabidopsis thaliana high-affinity phosphate transporters exhibit multiple levels of posttranslational regulation. Plant Cell 23, 1523-1535.

Bouain, N., Kisko, M., Rouached, A., Dauzat, M., Lacombe, B., Belgaroui, N., Ghnaya, T., Davidian, J.C., Berthomieu, P., Abdelly, C., and Rouached, H. (2014). Phosphate/zinc interaction analysis in two lettuce varieties reveals contrasting effects on biomass, photosynthesis, and dynamics of Pi transport. Biomed Res Int 2014, 548254.

Bouain, N., Shahzad, Z., Rouached, A., Khan, G.A., Berthomieu, P., Abdelly, C., Poirier, Y., and Rouached, H. (2014). Phosphate and zinc transport and signalling in plants: toward a better understanding of their homeostasis interaction. J Exp Bot 65, 5725-5741.

Bournier, M., Tissot, N., Mari, S., Boucherez, J., Lacombe, E., Briat, J.F., and Gaymard, F. (2013). Arabidopsis ferritin 1 (AtFer1) gene regulation by the phosphate starvation response 1 (AtPHR1) transcription factor reveals a direct molecular link between iron and phosphate homeostasis. J Biol Chem 288, 22670-22680.

Briat, J., Rouached, H., Tissot, N., Gaymard, F., and Dubos, C. (2015). Integration of P, S, Fe and Zn nutrition signals in Arabidopsis thaliana : potential involvement of PHOSPHATE STARVATION RESPONSE 1 (PHR1). Front Plant Sci 6.

Hamburger, D., Rezzonico, E., Macdonald-Comber Petetot, J., Somerville, C., and Poirier, Y. (2002). Identification and characterization of the Arabidopsis PHO1 gene involved in phosphate loading to the xylem. Plant Cell 14, 889-902.

Haydon, M.J., Kawachi, M., Wirtz, M., Hillmer, S., Hell, R., and Kramer, U. (2012). Vacuolar nicotianamine has critical and distinct roles under iron deficiency and for zinc sequestration in Arabidopsis. Plant Cell 24, 724-737.

Hirsch, J., Marin, E., Floriani, M., Chiarenza, S., Richaud, P., Nussaume, L., and Thibaud, M.C. (2006). Phosphate deficiency promotes modification of iron distribution in Arabidopsis plants. Biochimie 88, 1767-1771.

Holford, I. (1997). Soil phosphorus: its measurement, and its uptake by plants. Australian Journal of Soil Research 35, 227-240.

http://Faostat3.Fao.Org/Browse/Q/Qc/E, F.

Huang, C., Barker, S.J., Langridge, P., Smith, F.W., and Graham, R.D. (2000). Zinc deficiency up-regulates expression of high-affinity phosphate transporter genes in both phosphate-sufficient and -deficient barley roots. Plant Physiol 124, 415-422.

Jabnoune, M., Secco, D., Lecampion, C., Robaglia, C., Shu, Q., and Poirier, Y. (2013). A rice cis-natural antisense RNA acts as a translational enhancer for its cognate mRNA and contributes to phosphate homeostasis and plant fitness. Plant Cell 25, 4166-4182.

Jain, A., Sinilal, B., Dhandapani, G., Meagher, R.B., and Sahi, S.V. (2013). Effects of deficiency and excess of zinc on morphophysiological traits and spatiotemporal regulation of zinc-responsive genes reveal incidence of cross talk between micro- and macronutrients. Environ Sci Technol 47, 5327-5335.

Kang, X., and Ni, M. (2006). Arabidopsis SHORT HYPOCOTYL UNDER BLUE1 contains SPX and EXS domains and acts in cryptochrome signaling. Plant Cell 18, 921-934.

Kellermeier, F., Armengaud, P., Seditas, T.J., Danku, J., Salt, D.E., and Amtmann, A. (2014). Analysis of the Root System Architecture of Arabidopsis Provides a Quantitative Readout of Crosstalk between Nutritional Signals. Plant Cell 26, 1480-1496.
Khan, G.A., Bouraine, S., Wege, S., Li, Y., De Carbonnel, M., Berthomieu, P., Poirier, Y., and Rouached, H. (2014). Coordination between zinc and phosphate homeostasis involves the transcription factor PHR1, the phosphate exporter PHO1, and its homologue PHO1:H3 in Arabidopsis. *J Exp Bot* 65, 871-884.

Kj., L., and Td., S. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. *Methods*, 402–408.

Lin, S.I., Chiang, S.F., Lin, W.Y., Chen, J.W., Tseng, C.Y., Wu, P.C., and Chiou, T.J. (2008). Regulatory network of microRNA399 and PHO2 by systemic signaling. *Plant Physiol* 147, 732-746.

Lin, W.Y., Huang, T.K., Leong, S.J., and Chiou, T.J. (2014). Long-distance call from phosphate: systemic regulation of phosphate starvation responses. *J Exp Bot* 65, 1817-1827.

Luk, E., Carroll, M., Baker, M., and Culotta, V.C. (2003). Manganese activation of superoxide dismutase 2 in Saccharomyces cerevisiae requires MTM1, a member of the mitochondrial carrier family. *Proc Natl Acad Sci U S A* 100, 10353-10357.

Luk, E., Jensen, L.T., and Culotta, V.C. (2003). The many highways for intracellular trafficking of metals. *J Biol Inorg Chem* 8, 803-809.

Mathan, K.K., and Amberger, A. (1977). Influence of iron on the uptake of phosphorous by maize. *Plant Soil* 48, 413-422.

Misson, J., Raghothama, K.G., Jain, A., Jouhet, J., Block, M.A., Bligny, R., Ortet, P., Creff, A., Somerville, S., Rolland, N., Doumas, P., Nacry, P., Herrerra-Estrella, L., Nussaume, L., and Thibaud, M.C. (2005). A genome-wide transcriptional analysis using Arabidopsis thaliana Affymetrix gene chips determined plant responses to phosphate deprivation. *Proc Natl Acad Sci U S A* 102, 11934-11939.

Muller, J., Toev, T., Heisters, M., Teller, J., Moore, K.L., Hause, G., Dinesh, D.C., Burstenbinder, K., and Abel, S. (2015). Iron-dependent callose deposition adjusts root meristem maintenance to phosphate availability. *Dev Cell* 33, 216-230.

Neset, T.S., and Cordell, D. (2012). Global phosphorus scarcity: identifying synergies for a sustainable future. *J Sci Food Agric* 92, 2-6.

Nussaume, L., Kanno, S., Javot, H., Marin, E., Pochon, N., Ayadi, A., Nakanishi, T.M., and Thibaud, M.C. (2011). Phosphate Import in Plants: Focus on the PHT1 Transporters. *Front Plant Sci* 2, 83.

Ova, E.A., Kutman, U.B., Ozturk, L., and Cakmak, I. (2015). High phosphorus supply reduced zinc concentration of wheat in native soil but not in autoclaved soil or nutrient solution. *Plant and Soil*, 16.

Poirier, Y., and Bucher, M. (2002). Phosphate transport and homeostasis in Arabidopsis. *Arabidopsis Book* 1, e0024.

Poirier, Y., and Jung, J.Y. (2015). Phosphate Transporters. *Annual Plant Reviews, Phosphorus Metabolism in Plants*, 125.

Rai, V., Sanagala, R., Sinilal, B., Yadav, S., Sarkar, A.K., Dantu, P.K., and Jain, A. (2015). Iron Availability Affects Phosphate Deficiency-Mediated Responses, and Evidence of Cross-Talk with Auxin and Zinc in Arabidopsis. *Plant Cell Physiol* 56, 1107-1123.

Rai, V., Sanagala, R., Sinilal, B., Yadav, S., Sarkar, A.K., Dantu, P.K., and Jain, A. (2015). Iron Availability Affects Phosphate Deficiency-Mediated Responses, and Evidence of Cross-Talk with Auxin and Zinc in Arabidopsis. *Plant Cell Physiol* 56, 1107-1123.

Rouached, H., Arpat, A.B., and Poirier, Y. (2010). Regulation of phosphate starvation responses in plants: signaling players and cross-talks. *Mol Plant* 3, 288-299.

Rouached, H., Secco, D., Arpat, B., and Poirier, Y. (2011). The transcription factor PHR1 plays a key role in the regulation of sulfate shoot-to-root flux upon phosphate starvation in Arabidopsis. *BMC Plant Biol* 11, 19.
Rouached, H., Stefanovic, A., Secco, D., Bulak Arpat, A., Gout, E., Bligny, R., and Poirier, Y. (2011). Uncoupling phosphate deficiency from its major effects on growth and transcriptome via PHO1 expression in Arabidopsis. Plant J 65, 557-570.

Secco, D., Baumann, A., and Poirier, Y. (2010). Characterization of the rice PHO1 gene family reveals a key role for OsPHO1;2 in phosphate homeostasis and the evolution of a distinct clade in dicotyledons. Plant Physiol 152, 1693-1704.

Shanmugam, V., Tsednee, M., and Yeh, K.C. (2012). ZINC TOLERANCE INDUCED BY IRON 1 reveals the importance of glutathione in the cross-homeostasis between zinc and iron in Arabidopsis thaliana. Plant J 69, 1006-1017.

Stefanovic, A., Ribot, C., Rouached, H., Wang, Y., Chong, J., Belbahri, L., Delessert, S., and Poirier, Y. (2007). Members of the PHO1 gene family show limited functional redundancy in phosphate transfer to the shoot, and are regulated by phosphate deficiency via distinct pathways. Plant J 50, 982-994.

Sun, S., Gu, M., Cao, Y., Huang, X., Zhang, X., Ai, P., Zhao, J., Fan, X., and Xu, G. (2012). A constitutive expressed phosphate transporter, OsPht1;1, modulates phosphate uptake and translocation in phosphate-replete rice. Plant Physiol 159, 1571-1581.

Thibaud, M.C., Arrighi, J.F., Bayle, V., Chiarenza, S., Creff, A., Bustos, R., Paz-Ares, J., Poirier, Y., and Nussaume, L. (2010). Dissection of local and systemic transcriptional responses to phosphate starvation in Arabidopsis. Plant J 64, 775-789.

Van De Mortel, J.E., Almar Villanueva, L., Schat, H., Kwekkeboom, J., Coughlan, S., Moerland, P.D., Ver Loren Van Themaat, E., Koornneef, M., and Aarts, M.G. (2006). Large expression differences in genes for iron and zinc homeostasis, stress response, and lignin biosynthesis distinguish roots of Arabidopsis thaliana and the related metal hyperaccumulator Thlaspi caerulescens. Plant Physiol 142, 1127-1147.

Wang, Y., Ribot, C., Rezzonico, E., and Poirier, Y. (2004). Structure and expression profile of the Arabidopsis PHO1 gene family indicates a broad role in inorganic phosphate homeostasis. Plant Physiol 135, 400-411.

Ward, J.T., Lahner, B., Yakubova, E., Salt, D.E., and Raghothama, K.G. (2008). The effect of iron on the primary root elongation of Arabidopsis during phosphate deficiency. Plant Physiol 147, 1181-1191.

Wege, S., Khan, G.A., Jung, J.Y., Vogiatzaki, E., Pradervand, S., Aller, I., Meyer, A.J., and Poirier, Y. (2016). The EXS Domain of PHO1 Participates in the Response of Shoots to Phosphate Deficiency via a Root-to-Shoot Signal. Plant Physiol 170, 385-400.

Zheng, L., Huang, F., Narsai, R., Wu, J., Giraud, E., He, F., Cheng, L., Wang, F., Wu, P., Whelan, J., and Shou, H. (2009). Physiological and transcriptome analysis of iron and phosphorus interaction in rice seedlings. Plant Physiol 151, 262-274.

Zhou, Y., and Ni, M. (2009). SHB1 plays dual roles in photoperiodic and autonomous flowering. Dev Biol 331, 50-57.

Zhou, Y., Zhang, X., Kang, X., Zhao, X., Zhang, X., and Ni, M. (2009). SHORT HYPOCOTYL UNDER BLUE1 associates with MINISEED3 and HAIKU2 promoters in vivo to regulate Arabidopsis seed development. Plant Cell 21, 106-117.

Wissuwa, M., Ismail, A.M., and Yanagihara, S. (2006). Effects of zinc deficiency on rice growth and genetic factors contributing to tolerance. Plant Physiol 142, 731-741.
FIGURES LEGEND

**Figure 1.** Phenotypes of wild type (WT) rice plant and *Ospho1;1-1* and *Ospho1;1-2* mutant lines grown under individual and combinatorial ion deficiency stresses. Seedlings were grown for 21 days on Yoshida medium with different ionic concentrations: Complete Yoshida media (Ct) or lacking either individual phosphate (-Pᵢ), zinc (Zn) or iron (Fe); double elements -Pᵢ. Zn, -Pᵢ-Fe, -Zn-Fe or triple elements -Pᵢ-Fe-Zn. Photographs were taken at 21 days old. The experiment was repeated three times.

**Figure 2.** Biomass of wild type rice plant and *Ospho1;1-1* and *Ospho1;1-2* mutant lines grown on different mineral element supply. (A) wild type (WT), (B) *Ospho1;1-1* and (C) *Ospho1;1-2* mutant seedlings were grown on Yoshida medium with different ionic concentrations: Complete medium (Ct) or lacking either individual element phosphate (-Pᵢ), zinc (Zn) or iron (Fe); double elements -Pᵢ.Zn, -Pᵢ-Fe, -Zn-Fe or triple elements -Pᵢ-Fe-Zn. FW; Fresh weight is expressed as g/plant. DW; Dry weight. The experiment was repeated three times. The data are given as means ± se (n = 9). Asterisk indicates statistical significance, P < 0.05. Double asterisk indicates statistical significance, P < 0.01.

**Figure 3.** Phosphate, Zinc and Iron concentrations in wild type (WT) rice and *Ospho1;1* mutant lines plant grown under individual and combinatorial ions deficiency stresses. A, C and E present Pᵢ, Zn and Fe concentrations in the shoots respectively. B, D and F present Pᵢ, Zn and Fe concentrations in the roots respectively. After a 10 days of germination in ¼ full-strength of Yoshida media, plants were grown for 3 weeks on full-strength Yoshida medium with different ionic concentrations: Complete medium (Ct) or lacking either individual element Pᵢ, Zn or Fe; double elements PᵢZn, PᵢFe, ZnFe or triple elements PᵢFeZn. Shoots and roots were harvested separately and the Pᵢ and Zn and Fe were extracted and their concentrations were determined as indicated in material and methods. Pi concentrations are calculated per fresh weight. Zn and Fe concentrations are calculated per dry weight. The experiment was repeated three times. The data are given as means ± se (n = 9). Asterisk indicates statistical significance, P < 0.05. Double asterisk indicates statistical significance, P < 0.01.
**Figure 4.** Transcripts abundance of sense and anti-sense of *OsPHO1;1* and *OsPHO1;2* in *OsPHO1;3* as determined by qRT-PCR. *OsActin1* was used as an internal reference. qRT quantitative real-time-PCR was performed to analyze the *OsPHO1;1* and *OsPHO1;2* in *OsPHO1;3* mRNA levels in response to under multiple nutritional stresses: complete Yoshida medium (Ct) or lacking to either individual element phosphate (-Pi), zinc (-Zn) or iron (-Fe); double elements -Pi-Zn, -Pi-Fe, -Zn-Fe or triple elements -Pi-Fe-Zn. All values are relative to wild-type (WT) roots. The experiment was repeated three times. The data are given as means ± se (n = 9). Asterisk indicates statistical significance, P < 0.05. Double asterisk indicates statistical significance, P < 0.01.
<figure><img src="image_1.png" /><figcaption>Figure 1</figcaption></figure>
<figure><img src="image_2.png" /><figcaption>Figure 2</figcaption></figure>
<figure><img src="image_3.png" /><figcaption>Figure 3</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>Figure 4</figcaption></figure>