# Identification and interest of molecular markers to monitor plant Pi status

Laura Cuyas, Pascale David, Damien de Craieye, Sophia Ng, Mustapha Arkoun, Claude Plassard, Mohamadi Faharidine, Delphine Hourcade, Francesca Degan, Sylvain Pluchon, et al.

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

To cite this version:

Laura Cuyas, Pascale David, Damien de Craieye, Sophia Ng, Mustapha Arkoun, et al. Identification and interest of molecular markers to monitor plant Pi status. BMC Plant Biology, 2023, 23, pp.401. 10.1186/s12870-023-04411-8. hal-04346516

HAL Id: hal-04346516

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

Submitted on 15 Dec 2023

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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH

Open Access

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

# Identification and interest of molecular markers to monitor plant Pi status

Laura Cuyas<fnref n="1" />†, Pascale David<fnref n="2" />†, Damien de Craieye<fnref n="2" />, Sophia Ng<fnref n="2" />,<fnref n="3" />, Mustapha Arkoun<fnref n="1" />, Claude Plassard<fnref n="4" />, Mohamadi Faharidine<fnref n="5" />, Delphine Hourcade<fnref n="5" />, Francesca Degan<fnref n="5" />, Sylvain Pluchon<fnref n="1" /> and Laurent Nussaume<fnref n="2" />*

**Abstract**

**Background** Inorganic phosphate (Pi) is the sole source of phosphorus for plants. It is a limiting factor for plant yield in most soils worldwide. Due to economic and environmental constraints, the use of Pi fertilizer is and will be more and more limited. Unfortunately, evaluation of Pi bioavailability or Pi starvation traits remains a tedious task, which often does not inform us about the real Pi plant status.

**Results** Here, we identified by transcriptomic studies carried out in the plant model *Arabidopsis thaliana*, early roots- or leaves-conserved molecular markers for Pi starvation, exhibiting fast response to modifications of phosphate nutritional status. We identified their homologues in three crops (wheat, rapeseed, and maize) and demonstrated that they offer a reliable opportunity to monitor the actual plant internal Pi status. They turn out to be very sensitive in the concentration range of 0-50 µM which is the most common case in the vast majority of soils and situations where Pi hardly accumulates in plants. Besides in vitro conditions, they could also be validated for plants growing in the greenhouse or in open field conditions.

**Conclusion** These markers provide valuable physiological tools for plant physiologists and breeders to assess phosphate bio-availability impact on plant growth in their studies. This also offers the opportunity to cope with the rising economical (shortage) and societal problems (pollution) resulting from the management of this critical natural resource.

**Keywords** Arabidopsis, Crops, Maize, Molecular markers, Pi deficiency, Rapeseed, Soil, Wheat

**Background**

Phosphorus (P) is an essential macronutrient for plants. It is absorbed as inorganic orthophosphate (Pi). It is a key component for cellular compartmentalization (membrane lipids), heredity (nucleic acids) and energy metabolism (ATP); it is also a crucial actor for signaling and cellular reactions through phosphorylation based mechanisms [1–6].

Despite its relative abundance on earth (11st most abundant element), Pi is unevenly distributed (Tiessen, 2008). It is estimated that one third of total cultivated soils are lacking phosphorus for optimum plant growth [7, 8]. Indeed, many characteristics of Pi explain the problems of Pi nutrition. First of all, Pi has a very poor mobility in soils, leading to a major part of the Pi fertilizers

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

BMC

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

<footnote>†Laura Cuyas and Pascale David contributed equally to this work.</footnote>

<footnote>*Correspondence:
Laurent Nussaume
lnussaume@cea.fr</footnote>

<footnote n="1">¹TIMAC AGRO, Laboratoire de Nutrition Végétale, Agrolnovation International, 18 Avenue Franklin Roosevelt, 35400 Saint-Malo, France</footnote>

<footnote n="2">²Aix Marseille Univ, CEA, CNRS, BIAM, UMR7265, EBMP, 13115 Saint-Paul Lez Durance, France</footnote>

<footnote n="3">³Centre for AgriBioscience, La Trobe University, 5 Ring Road Bundoora, Victoria 3086, Australia</footnote>

<footnote n="4">⁴INRAE, CIRAD, IRD, Univ Montpellier, Eco&Sols, Institut Agro, 34060 Montpellier, France</footnote>

<footnote n="5">⁵Arvalis, Institut du Végétal, Station Expérimentale, Boigneville, France</footnote>
spread to be recovered by microorganisms at the expense
of the crops. Secondly, it forms insoluble complexes with
many soil cations or chelates with clays reducing its bio-
availability [9–13]. These phenomena explain why it is
estimated that only 20% of the Pi applied as a fertilizer is
actually used by plants [14]. It is interesting to notice that
many microorganisms solubilizing insoluble phosphate
form are reported to have positive effects on crop pro-
ductivity (by increasing Pi nutrition). Nevertheless, such
studies are mostly performed in controlled conditions
and could not be reproduced in fields where solubilized
Pi benefits mainly to microbial biomass [15].

In the 1960s, there was a big concern about how the
world would be able to feed itself and face the upcom-
ing increasing population, which doubled between the
1960s and the 1990s. This has led to significant changes
in agricultural politics known as the Green Revolution. It
was based on the combination of high-yielding varieties,
irrigation, mechanization and use of chemical fertilizers.
This promoted a constant increase of agricultural pro-
ductivity worldwide [16–18]. Nevertheless, the challenge
is still going on: the world population is expected to grow
approximately 30% in the next future years reaching up to
10 billion in 2050 [19–21]. As a result, the United Nations
Food and Agriculture Organization (FAO) estimates that
world food production should increase by 70% in order
to achieve global food security [21, 22]. However, there is
almost no more available cropland for future expansions
[22, 23], transferring the pressure on the increase of agri-
cultural productivity [24, 25].

As a major macronutrient, Pi fertilizers are crucial
actors of this policy to ensure grass and plant crops pro-
ductivity [26]. Nonetheless, improper use of Pi fertiliz-
ers causes severe damages to the environment such as
eutrophication of rivers and lakes due to algae blooms:
the leak of Pi excess to rivers promoted development of
toxic cyanobacteria [27-30], or metal pollution. The vast
majority of Pi rock, having sedimentary origin, contains
high levels of toxic metals such as cadmium, chromium,
mercury, leads, uranium or thorium [31-33]. These are
therefore often found in Pi fertilizers if they are not suf-
ficiently purified, leading to novel regulation. Thus,
recently, the Council of the European Union adopted a
regulation on Pi fertilizers, limiting the cadmium con-
tent at 60 mg/kg P<sub>2</sub>O<sub>5</sub> (Commission Decision (EU)
2020/1205). A value that will be reduced by three over
the next 12 years [34]. Such measure is therefore likely
to limit the use of the Pi fertilizers in the future. This is
all more important since Pi is not a renewable source and
the available resources must be managed in a reasoned
manner [35, 36]. It is therefore essential and mandatory
to rationalize fertilizer uses. In the future, to reach preci-
sion farming, we need a tight control of the amount of Pi

fertilizer supplied to optimize Pi use efficiency and guar-
antee a correct balance between the amount of Pi pro-
vided and the part absorbed by the plant [37, 38].

Pi deficiency is diagnosed by multiple ways on the
field, but no one is fully satisfactory. Regular analyses
of the soil combined with specific fertilization accord-
ing to the crop and soil Pi status remains the most com-
mon approach. This methodology is presenting two clear
disadvantages: soil plots are not homogenous due to low
Pi mobility, so results would depend on the place/depth
soil sample has been collected; and although several
methods to determine current Pi soil stocks and criti-
cal Pi values exist, they are presenting high variabilities
among them and they are not always efficient enough to
diagnose Pi deficiency nor to reflect real bioavailable Pi
for plants [39].

At the plant level, the difficulties resulted from the
absence of clear traits. The main Pi deficiency visual
symptoms resulted from anthocyanins accumulation,
promoting dark green/reddish purple colors in the plant
and reduced plant growth. Unfortunately, the apparition
of visual symptoms implies already an irreversible altered
plant development. Besides, they lack specificity, as
anthocyanins accumulation is a common plant response
to many stresses [40]. Different physical or biochemi-
cal measures to assay Pi content are also available [41].
Unfortunately, measuring anion Pi directly in the plant,
is also not very informative as it accumulates only when
present in excess [42], a situation rarely observed in field
conditions. In addition, Pi values vary depending on the
crop, age of the plant; organs studied or even its crosstalk
with other nutrients, complicating its correlation with
real Pi status of the plant.

Therefore, we decided to investigate if molecular mark-
ers can provide reliable indication for rationalizing the
phosphate application to plants. Indeed, plant responses
to Pi starvation is a tightly regulated process, where tran-
scriptomic regulation plays a major role [2, 43-48]. Using
transcriptomic analysis on model plant *Arabidopsis*, we
looked for markers conserved in crops exhibiting fast and
important sensitivity response to Pi presence. Then we
tested their capacities to monitor plant Pi starvation sta-
tus in *Arabidopsis* but also in three different crops (rape-
seed, maize and wheat) to confirm capacities to extend
their use for plants of agronomic interest.

**Materials and methods**

**Plant materials and growth conditions**

Wild-type *A. thaliana* Col-0 seeds were sterilized
and grown vertically in Petri dishes [49]. The modi-
fied Murashige and Skoog medium contained 0.47 mM
MgSO₄, 2.1 mM NH₄NO₃, 1.89 mM KNO₃, 0.67 mM
CaCl₂, 0.5 µM KI, 0.79 mM H₃BO₃, 10 µM MnSO₄, 5 µM
ZnSO₄, 1 µM Na₂MoO₄, 0.1 µM CuSO₄, 0.1 µM CoCl₂, 5 g L⁻¹ sucrose and 3.4 mM 2-(N-morpholino) ethanesulfonic acid (MES) buffered at pH 5.8 with KOH. The agar (8 g L⁻¹) for plates was from Sigma-Aldrich (A7921 Lot BCBZ7284, see Table S1 for elemental composition). Plates were supplemented with 2 µM of FeCl₂ and a range of KH₂PO₄ (0 to 1500 µM). The growth chamber conditions were 16-h-light (25 °C)/8-h-dark (22 °C). Seeds of Brassica napus (cv. Adriana), Zea mays (cv. Ronaldinho) and Triticum aesticum (cv. Rubisko, Fluor and Johnson) were used. For greenhouse experiments, seeds were surface-sterilized for 10 min in 70% (v/v) bleach solution supplemented with 0.01% (v/v) Tween-20 (Sigma P1379) and washed five times with bi-distilled water. Seeds were germinated on vermiculite (equilibrated with bi-distilled water) for 3 days in dark and then in white light for 4 days (maize and wheat) or 11 days (rapeseed). Seedlings were then transferred to 2L pots containing a 50:50 mixture of silica sand and vermiculite (one seedling per pot for maize and wheat, and two seedlings per pot for rapeseed) or to a soil growing substrate (1L soil:1L silica sand). We used 3 soils whose main physical and chemical characteristics are given in Table S2. Seeds were removed once transferring to avoid Pi remobilization from seeds. Plants were watered twice per week at 80% water holding capacity with different nutrient solutions according crop species (Table 1).

Solutions were supplemented with different concentrations of KH₂PO₄ (0, 50, 100, 250, 500, 1000 and 1500 µM) and KCl (1500, 1450, 1400, 1250, 1000, 500 and 0 µM) respectively. Plants were grown in the greenhouse under a 16-h photoperiod, at 25 °C light/22 °C dark and at 60% of humidity, for 30 days. Field experiments with wheat were performed on the French commune of Giroussens

located in the South of France. The experimental site is located in one of low ancient fluvial terraces, presenting silty and acid soils (pH: 6.1) locally called 'boulènes', a sub-group of Planosols characterized by the thickness of the silty layer (or the depth of the clay layers). Soil sampling plan followed a regular grid design, with a total of sixty soil sampling points. For each point, the P₂O₅ Olsen mg/kg was measured. The quantification limit is 10 of P₂O₅ Olsen mg/kg and the uncertainties for 30 is ±4,7 mg/kg P₂O₅ Olsen. In the P fertilized plot, the average value of P₂O₅ Olsen is 38 mg/kg (before fertilization) and in the non-fertilizer part of the trial (P-, South-Est) the value is 20 mg/kg P₂O₅ Olsen. Phosphate fertilization of the field was performed by a surface application of triple super phosphate fertilizers (120 kg P₂O₅/ha).

## Physiological measurements

### Fresh weight measurements

Shoots and roots were harvested separately and weighed directly. All seedlings were collected individually except for *Arabidopsis* (harvested in pools of 5 to 20 plantlets) or rapeseed seedlings (harvested in pools of two plantlets). Samples were immediately frozen in liquid nitrogen and stored at -80 °C until other assays were performed.

### Chlorophyll and flavonol estimations

Spectral analyses were performed by using the DUALEX photometer (DUALEX v4.5, Force A, France, https://www.force-a.com/fr/produits/dualex). Measurements were taken from the middle region of the second fully expanded leaf (either from rapeseed, maize, or wheat) the day before harvesting.

### Measurement of total P and free Pi cellular content

Cellular free Pi measurements were performed as previously described using the malachite green assay [41]. For total phosphorus, lyophilized leaf extracts were digested in concentrated 14.5 N nitric acid to convert organic P into mineral Pi. The samples were then diluted with water to reduce the nitric acid concentration below 0.1 N, and the Pi was measured as described above.

### Measures of plant biomass for field experiments

They were performed by a Phantom 4 Pro V2 model, with a mounted camera Micasense RedEdge. The camera calibration was realized thanks to the calibration target provided by Micasense and specific integration on the drone, by means of a 3D printing. There are four steps for image acquisition. First, the control of the ground points before the flight for the positioning of the image is acquired. Then, the control of the constant brightness (clear sky or an overcast sky without clearings) is carried out before launching the flight. During the flight, the images are

<table><caption>**Table 1** Composition and concentrations (mM) of nutrient solutions used to water the pots</caption><thead><tr><th rowspan="2">Salt</th><th colspan="3">Plant species</th></tr><tr><th>Rapeseed</th><th>Maize</th><th>Wheat</th></tr></thead><tbody><tr><td>KNO<sub>3</sub></td><td>5</td><td>2.5</td><td>1.25</td></tr><tr><td>Ca(NO<sub>3</sub>)<sub>2</sub></td><td>3.125</td><td>2.5</td><td>2.5</td></tr><tr><td>MgSO<sub>4</sub></td><td>1.25</td><td>0.5</td><td>0.25</td></tr><tr><td>MgCl<sub>2</sub></td><td>3.75</td><td>-</td><td>-</td></tr><tr><td>CaCl<sub>2</sub></td><td>-</td><td>-</td><td>0.25</td></tr><tr><td>H<sub>3</sub>BO<sub>3</sub></td><td>0.035</td><td>0.0575</td><td>0.0575</td></tr><tr><td>MnSO<sub>4</sub></td><td>0.0125</td><td>0.025</td><td>0.0125</td></tr><tr><td>ZnSO<sub>4</sub></td><td>0.0075</td><td>0.01</td><td>0.005</td></tr><tr><td>CuSO<sub>4</sub></td><td>0.00175</td><td>0.0025</td><td>0.00225</td></tr><tr><td>(NH<sub>4</sub>)<sub>6</sub>Mo<sub>7</sub>O<sub>24</sub></td><td>0.00175</td><td>-</td><td>0.00075</td></tr><tr><td>CoCl<sub>2</sub></td><td>0.00025</td><td>0.00025</td><td>0.00025</td></tr><tr><td>NaFe(III)EDTA</td><td>0.5</td><td>0.75</td><td>0.25</td></tr></tbody></table>
individually calibrated during the acquisition using an on-board Downwelling Light Sensor, called "DLS2". Before and after the flight, an image of the ground target calibration is taken immediately before and after the acquisition. The raw data were processed with the equation indicated below and by averaging the reflectance values by microplot. The reflectance values extracted from the spectroradiometers (band width of ±5 nm around the targeted central wavelength filtered with a gated method) were used to calculate the Normalized Difference Vegetation Index (NDVI = (p<sub>NIR</sub> - p<sub>RED</sub>)/(p<sub>NIR</sub> + p<sub>RED</sub>); Rouse et al., 1974). Where: λRED = 675 nm; λNIR = 785 nm that already demonstrated its usefulness for crop phenotyping (Comar et al., 2012).

## Gene expression analyses

For *Arabidopsis*, the extraction of total RNA from roots or shoots and the RT-qPCR experiments were performed as previously described [50]. RNAseq experiments were performed as previously described [43].

For plant crops, Nucleospin 8 RNA kit (Macherey-Nagel) was used to isolate the RNA from the different plant species tested (50 mg root powder/sample). The quality and concentration of all samples were checked using 4200 Tapestation (Agilent Technologies), followed by DNase treatment and cDNA synthesis from 1 µg of RNA (iScriptTM gDNA). RT-qPCR reactions were performed in a Real-Time PCR Detection System (BIO-RAD) using a total of 10 µl reaction containing 5 µL of Universal SYBR Green Supermix (BIO-RAD) and primers at 0.5 µM. All reactions were performed in technical triplicates. Primers were designed with Primer3. The list of all primers used is provided in Table S3. Relative expression changes were calculated by the ΔCCq method. The exponential expression is calculated as $2^{-ΔCq}$, in which ΔCq is the difference between the Cq of the phosphate starvation induced gene analyzed and the average of Cq obtained from all housekeeping genes used. Values were then normalized to corresponding control.

## Measurement of available Pi in soils

Available Pi was measured either following Olsen extraction of soil performed as described in [41] or after binding to anionic exchange membranes (AEM). In this last case, we used 2×2 cm membranes (Selemion AMV anionic exchange membranes, AGC Engineering) previously gently washed 8 times for 10 min with constant shaking in 20 mL of 0.5 M KCl. They were incubated with 4 g of soil at 80% of water holding capacity for 24 h. AEM were then rinsed with bi-distilled water to remove soil and then immersed in 3.75 mL of 0.5 M KCl and gently shaken for 10 min. This procedure was repeated 6 times and the solution was collected each time. Measurement

of orthophosphate released (either in Pi Olsen or in AEM assays) was carried out using the malachite green method, as previously described by [41].

## In silico analyses

We used the Genevestigator database, using only datasets corresponding to wild-type (Col-0) tissues and selected the conditions in the database: AT_mRNASeq_ARABI/GL-9.

## Statistical analyses

Statistical analyses were conducted by performing one-way ANOVA and significant differences were analyzed by SNK Test. Means are marked by different letters for values that were significantly different (*p* < 0.05). All statistical analyses were performed with RStudio software.

## Results

### Identification of Pi deficiency molecular markers

To date, different molecular markers responding differently to Pi status have been described (Misson et al., 2005; Bustos et al., 2010). Many belong to multigenic families where a majority of their members respond to Pi deficiency such as phosphate transporters, purple acid phosphatases, the SPX transcriptional inhibitors [2, 51–58]. Besides, many genes responding to Pi deficiency are also triggered by additional biotic or abiotic stresses [2, 59]. In addition, many genes have close homologs that may be regulated differently (such as *Arabidopsis PHT1*;1, *PAP3*, *PAP4* or *SPX4*). Therefore, it is necessary to ensure the specificity of the selected markers to avoid false positives due to these other parameters. In order to select optimal molecular markers, we have: 1) performed data mining to identify fast responding genes, acting systemically and specifically to Pi deficiency; 2) selected those which are conserved among different plant species and 3) chosen those which are exhibiting important variation of transcripts between -Pi/ + Pi to provide a good range of sensitivity for the assay.

We performed RNAseq analysis with *Arabidopsis* to identify molecular markers exhibiting fast response to modification of Pi homeostasis [43, 60]. To reach this goal, we starved plants and performed Pi refeeding experiments. After transferring seedlings to plates containing Pi, we harvested roots and leaves after 1 and 3 h for RNA-seq analysis. We selected 6 genes exhibiting a significant two-fold reduction of their transcript level at 3 h (Table 2). In addition, we selected *AtPHR1*, a major regulator for Pi homeostasis that is transcriptionally very stable (Table 2; Rubio et al., 2001) as an additional control for the transcriptomic analyses of this work (besides housekeeping genes used for normalization). We also included *AtIPS1*, a non-protein coding gene that plays
<table><caption>**Table 2** List of genes selected from RNAseq experiment responding fast to Pi replenish. Arabidopsis seedlings were grown under -Pi for 7 days and supplemented with 0.5 mM Pi for 1 and 3 h. Results from seedlings grown for 7 days under +Pi and transferred to +Pi for 3 h are used as control +Pi (ctrl). Values from shoots and roots are shown. Results are expressed in log2FC (foldchange), comparing each value to -Pi. Values correspond to three biologic independent replicates</caption><thead><tr><th rowspan="2">Accession n°</th><th rowspan="2">Name</th><th colspan="3">Shoots log2FC</th><th colspan="3">Roots log2FC</th></tr><tr><th>re 1h</th><th>re 3h</th><th>ctrl + Pi</th><th>re 1h</th><th>re 3h</th><th>ctrl + Pi</th></tr></thead><tbody><tr><td>AT3G09922</td><td>IPS1</td><td>-0.11</td><td>-0.51</td><td>-9.51</td><td>-0.15</td><td>-0.82</td><td>-9.00</td></tr><tr><td>AT4G28610</td><td>PHR1</td><td>-0.06</td><td>0.16</td><td>0.34</td><td>0.15</td><td>0.12</td><td>0.20</td></tr><tr><td>AT2G38940</td><td>PHT1;4</td><td>-0.90</td><td>-2.79</td><td>-4.28</td><td>-1.95</td><td>-1.85</td><td>-4.66</td></tr><tr><td>AT5G20150</td><td>SPX1</td><td>-1.49</td><td>-3.20</td><td>-5.36</td><td>-2.01</td><td>-2.13</td><td>-6.23</td></tr><tr><td>AT2G45130</td><td>SPX3</td><td>-0.95</td><td>-4.61</td><td>-11.05</td><td>-2.61</td><td>-2.37</td><td>-9.56</td></tr><tr><td>AT5G01220</td><td>SQD2</td><td>0.52</td><td>-2.74</td><td>-3.85</td><td>-0.81</td><td>-1.93</td><td>-4.10</td></tr><tr><td>AT5G20790</td><td>UNICORN1</td><td>-1.94</td><td>-3.34</td><td>-6.74</td><td>-3.02</td><td>-2.43</td><td>-6.36</td></tr></tbody></table>

an important role to tune finely Pi homeostasis [61] yet exhibiting a much slower regulation (Table 2). As previously described [43], the regulation was first detected in the root (1 h) but after 3 h the signal amplitude became more important in the aerial part. All of them (except *AtPHR1*) were very tightly repressed after the addition of Pi (around 7 Log2 fold change on average).

As described by [47], it is important to distinguish the expression of systemic genes (determined by internal levels of Pi) from locally regulated ones (dependent on changes in root growth responding to soil composition [62]). For those last ones, low Pi concentration in the medium increases the bioavailability of many cations including metals such as iron or aluminum [10, 50]. Therefore, many genes induced locally by Pi deficiency are in reality not responding directly to Pi concentration present in the soil but to the presence of high bioavailable metal(s) [50]. Therefore, it is essential to choose genes systemically regulated by plant Pi status. This is the case of the markers selected here, which all present P1BS regulatory box(es) in their promoter [47] P1BS regulatory box(es) is/are requested for the binding of PHR1 transcription factor [44, 63], the master gene controlling the main regulations associated with Pi homeostasis.

In order to verify that the selected genes are specific for Pi deficiency, we carried out in silico investigation (Fig. S1). Pi starvation promoted an important induction of all the markers (by 4 to 11-Log2 fold), whereas all other biotic and abiotic stresses investigated, promoted none or very limited modifications (far below the twofold threshold classically used).

### Phosphate starvation response in Arabidopsis

To examine the impact of Pi conditions in *Arabidopsis* and to prove the effectiveness of the selected molecular markers, we grew *Arabidopsis* seedlings in vitro

with 0, 5, 10, 20, 50, 100, 250, 500 and 1500 µM of KH₂PO₄ for 14 days. We analyzed fresh weight as well as Pi content in addition to the expression of selected genes. Results indicated that above 50 µM, plants started to accumulate Pi and reached optimal growth at 250 µM (Fig. 1a, b). Shoot-to-root ratio is a good indicator to identify optimal Pi-condition, where plants prioritize shoot to root development. Below 20 µM Pi present in the medium, plants favor root development. Concomitantly with Pi accumulation in leaves, the increase of the aerial part with respect to the root system was observed (Fig. 1b). We also analyzed the Pi accumulation in the shoots and the roots. This accumulation was only observed above 50 µM in both root and shoot illustrating that plant favors development prior to Pi storage (Fig. 1c). In agreement with leaf development preferred over root development, Pi accumulation occurred faster in the leaves in comparison to the roots.

The expression of selected molecular markers indicates that a supply of 5 to 10 µM KH₂PO₄ was already significantly reducing their expression (Fig. 1d). They all decrease in inverse proportion to the increase of Pi concentration present in the medium: the lowest levels being reached starting from 250 µM. As observed for Pi accumulation or growth, the markers responded faster in the leaves and reached their minimal level as soon as optimal growth was observed. It should be noticed that they are sensitive enough to detect differences between 0 to 20/50 µM Pi whereas neither growth nor Pi accumulation could do the same. *AtPHR1* was used as control in this experiment and showed no modulation of its expression in response to various Pi supplies as previously described [63].

Then, we investigated if we could extend the use of these markers to different crops for monitoring their Pi status.
<figure><img src="image_5.png" /><figcaption>b</figcaption></figure>

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

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

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

**Fig. 1** Growth of Arabidopsis seedlings exposed to different concentrations of Pi for 14 days. a Photos of different modalities. Scale bar = 1 cm. b Shoot and root fresh weight (mg) and shoot-to-root ratio per seedling. c Shoot and root total free-Pi content (ng Pi/mg fresh weight; malachite green method). d Relative gene expression of selected molecular markers. All values are relative to control (0 µM Pi roots or shoots), which is normalized to 1. *AtTUBULIN* was used as housekeeping gene. Bars indicate means ± SD. Different letters indicate significant different means (one-way ANOVA followed by SNK test, *p* < 0.05, *n* = 2–4 biological replicates, each replicate corresponds to a pool of different seedlings, weight has been extrapolated per seedling). \*More independent experiments have been performed that displayed the same tendency, data not shown

## Application of Pi starvation markers in rapeseed, a dicotyledonous crop

To control Pi application during crop development, we grew the plants on inert substrate (silica sand and

vermiculite) in a greenhouse for one month, watered at 80% of water holding capacity with nutrient solution supplemented with different doses of Pi (0, 50, 100, 250, 500, 1000 and 15,000 µM KH₂PO₄ corresponding to 0,
1.16, 2.32, 5.81, 11.62, 23.23, 34.85 mg P/kg substrate) for 30 days. The choice of this inert substrate is related to its low Pi content and the facility it offers to access the root compartment.

Similar to *Arabidopsis*, monitoring of plant Pi status was performed through a set of different parameters related to plant growth (shoot and root fresh weight, as well as shoot-to-root ratio), nutrient content (total Pi content either in shoots or in roots), and the expression of corresponding homologs of previously selected molecular markers (Table 2). In addition, we also used analysis provided by DUALEX®, a leaf clip sensor measuring light absorption spectra, to get access to epidermal chlorophylls and flavonols (anthocyanin precursors) content.

The data indicate that this species is very sensitive to Pi deficiency. A strong differential growth was observed between 0 to 500 µM Pi where an optimal growth was observed (Fig. 2a and b). Interestingly, both root and shoot reacted to Pi supply. Unlike *Arabidopsis*, the shoot-to-root ratio could not be discriminated between the ranges of applied Pi (Fig. 2b).

Chlorosis promoted by Pi deficiency was very mild (Fig. S2a) as illustrated by the 15% reduction of chlorophyll content observed in the range of 0/100 µM Pi. The accumulation of anthocyanin precursor was also observed but remained very limited (Fig. S2a).

We then studied the Pi accumulation. Seedlings started to accumulate Pi either in shoots or in roots from 250 µM Pi, reaching the highest accumulation at 1500 µM in both organs. In rapeseed, unlike *Arabidopsis*, plants did not favor Pi accumulation in the leaves versus root and the concentration was fairly similar in both organs except when Pi is very limiting (range 0 to 100 µM Pi here; Fig. 2c).

The use of selected molecular markers provides a clear opportunity to identify plants exhibiting Pi starvation traits (Fig. 2d). It should be noticed that *BnSQD2* (like few other markers involved in phospholipids replacement) was not tested but could be expected to be also a good marker for Pi presence. As it could be expected for Pi transporters (mainly expressed in root), *BnPHT1* turns out to be very effective only in such organ and could significantly distinguish absence of Pi (0 mM) from limiting Pi (50 to 500 µM) and excess of Pi (above 500 µM). All other selected markers distinguished mostly Pi limiting from non-limiting conditions. In leaves, all markers could also identify these two conditions and interestingly, they could also distinguish the plants grown with 250 µM Pi, which for the aerial part discriminated the border between these two categories of conditions. The dynamic of the response of the markers between these two categories was very important (reaching a 1/100 ratio) for the majority of probes with the exception of *BnPHT1*,

which is mainly expressed in roots. It is also interesting to notice that the marker's expression was inversely proportionally correlated with the plant development measured by fresh weight.

## Application of molecular markers to detect Pi starvation in monocotyledonous such as maize or wheat

In maize, the significant changes for fresh weight were very limited (due to early stage of the analysis) and restricted to increased shoot to root growth for plants receiving 250 (or more) µM Pi solution (Fig. 3a and b). Clear accumulation of anthocyanin at the bottom of the stem could be observed. Those pigments turn out to be not present in all tissue in DUALEX analysis, which is restricted to leaves and could not identify any significant changes between samples (Fig. S2b). There was also an absence of chlorosis as chlorophyll pigment did not exhibit noticeable discrepancies between treatments. Nevertheless, analysis of Pi content correlated the 20/30% increase of leaves versus root growth observed with the situation where Pi starts to be accumulated in shoots (250 µM Pi supplied; Fig. 3c). In roots, such accumulation was also observed but delayed (starting when 1000 µM Pi is supplied; Fig. 3c).

All molecular markers could also detect plants accumulating Pi but with clear different sensitivities (Fig. 3d). Whereas *ZmSPX3* exhibited in both roots and leaves a two-Log10 dynamic of response among the range of Pi concentration tested, it was reduced by 35 times in leaves and by 8 times in roots for *ZmSPX1* and even less for *ZmPHT2* or *ZmSQD2*.

Then, we studied wheat plants, which presented a very different response from maize. Indeed, this species turned out to be very sensitive to Pi supply and reached a maximum of growth starting from 50 µM for root and 250 µM for leaves, illustrating here a strategy to favor soil exploration over leaf development (Fig. 4a and b). Like for rapeseed, chlorosis was very mild and only observed when no Pi was provided and flavonol presence turned out to be similar in all samples (Fig. S2c).

As previously observed for most species studied here, the accumulation of Pi was favored in the aerial part, where it was observed as soon as 250 µM Pi is provided, whereas it required the double concentration for the roots (Fig. 4c). Molecular markers reflected well those features with a strong decrease starting once Pi is added to the medium and reaching a minimum at 250 or 500 µM Pi according to the marker used (Fig. 4d). The dynamic was pretty good for some markers such as *TaSPX1* or *TaIPS1*, they became undetectable beyond addition of 500 µM Pi in the nutrient solution.
<caption>a</caption>

µM Pi (KH₂PO₄)

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

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

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

<figure><img src="image_12.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>d</figcaption></figure>

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

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

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

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

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

**Fig. 2** Growth of rapeseed plants exposed to different concentrations of Pi for 30 days. a Photos (two plants/pot). Scale bar = 5 cm. b Shoot and root fresh weight (grams) and shoot-to-root ratio. c Shoot and root total free-Pi content (ng Pi/mg fresh weight; malachite green method). d Relative gene expression of selected molecular markers. All values are relative to control (0 µM Pi roots or shoots), which is normalized to 1. *BnACTIN-7* and *BnEF-1alpha* were used as housekeeping genes. Bars indicate means ± SD. Different letters indicate significant different means (one-way ANOVA followed by SNK test, *p* < 0.05, *n* = 7–14 for b, *n* = 5–7 for c-d)

## Importance of the variations within the plants

It is well known that all organs do not have similar Pi status during the growth of the plants. To investigate this point, we checked different aerial parts on wheat plants grown at 250 µM Pi. We distinguished the leaf sheath and young, intermediary, or old leaves (Fig. 5a). As expected, the old leaves presented senescence traits (Fig. 5b), and they did not store Pi. They exhibited a similar Pi content comparable to plants growing without Pi.

As a consequence, the induction of molecular markers (Fig. 5c) was very important (mostly identical to the ones observed for plants growing on 0 µM Pi). In contrast, the young developing leaves, which exhibited a sink status, stored a lot of Pi and presented a strong repression of the molecular markers of Pi starvation (Fig. 5c). Intermediary leaves exhibited an average situation in terms of both Pi content and level of expression of the molecular markers. The leaf sheath also presented an important Pi
<figure><img src="image_19.png" /></figure>

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

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

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

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

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

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

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

<figure><img src="image_27.png" /><figcaption>**Fig. 3** Growth of maize plants exposed to different concentrations of Pi for 30 days. a Photos (one plant/pot). Scale bar = 5 cm. b Shoot and root fresh weight (grams) and shoot-to-root ratio. c Shoot and root total free-Pi content (ng Pi/mg fresh weight; malachite green). d Relative gene expression of selected molecular markers. ZmEIF4A and ZmTUBULIN were used as housekeeping genes. All values are relative to control (0 µM Pi roots or shoots), which is normalized to 1. Bars indicate means ± SD. Different letters indicate significant different means (one-way ANOVA followed by SNK test, *p* < 0.05, *n* = 5–7)</figcaption></figure>

content probably connected with its essential role for Pi
translocation and distribution from roots to leaves. As a
result, molecular markers were also strongly repressed
in this organ (Fig. 5c). In conclusion, molecular markers
turn out to be, on average, connected with the Pi content
of these different organs. It prompts us to circumvent

the important variability observed between the differ-
ent organs by sampling the entire aerial parts. We nev-
ertheless discarded the senescent leaves which were easy
to recognize to reduce the bias introduced in the meas-
ures by this important source of molecular Pi marker
induction.
<figure><img src="image_28.png" /><figcaption>**Fig. 4** Growth of wheat plants exposed to different concentrations of Pi for 30 days. a Photos (one plant/pot). Scale bar = 10 cm. b Shoot and root fresh weight (grams) and shoot-to-root ratio. c Shoot and root total free-Pi content (ng Pi/mg fresh weight; malachite green). d Relative gene expression of selected molecular markers. TaACT and TaGAPDH were used as housekeeping genes. All values are relative to control (0 µM Pi roots or shoots), which is normalized to 1. Bars indicate means ± SD. Different letters indicate significant different means (one-way ANOVA followed by SNK test, *p* < 0.05, *n* = 4–7)</figcaption></figure>

Validation of molecular Pi markers to monitor Pi status
of different soils

To assay robustness of the markers once using agronomic
soils, we grew wheat plants for 30 days in three different
soils exhibiting similar total phosphorus content (Table
S2), but distinct levels of Pi due to their different prop-
erties (soil 1: clay soil, deficient on Pi, presenting high

binding capacities to Pi; soil 2: chalky soil, with low/
sufficient Pi content; soil 3: sandy soil, rich on Pi). To
ensure that other nutrients were not limiting, the differ-
ent soils were watered twice a week with a Pi-depleted
wheat nutrient solution -Pi (described in Materials and
Methods). Although several methods exist to determine
soil available Pi, they all over- or underestimate the plant
a

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

b

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

c

<figure><img src="image_31.png" /><figcaption>**Fig. 5** Compartmental response to Pi status of shoots. a Schema of different compartments analyzed. b Free-Pi content in different aerial parts of the plant (ng Pi/mg fresh weight). c Relative gene expression of selected molecular markers in shoots. TaACT and TaGAPDH were used as housekeeping genes. All values are relative to control (0 µM Pi), which is normalized to 1.0 and 250 design plants grown in silica:vermiculite (1:1) with Pi supply of 0 and 250 µM, respectively. Letters correspond to S (leaf Sheath), Y (Young leaves), I (Intermediate leaves) and O (Old leaves). Bars indicate means ± SD. Different letters indicate significant different means (one-way ANOVA followed by SNK test, *p* < 0.05, *n* = 6–7)</figcaption></figure>

available Pi content. In this work, we have used two different approaches: the Pi Olsen extraction, and the anionic exchange membrane; which respectively are known to putatively overestimate or underestimate the bioavailable Pi present in the soil (Fig. 6a). Both approaches indicated that soil 1 was deficient on Pi, whereas that soil 3 could be considered as rich and soil 2 with intermediary Pi content. The difficulty to study Pi content by using these methods is illustrated by the differences obtained between soils. Both techniques suggested that soil 3 presents 4 to 4.5 more Pi than soil 1. The precision was reduced for soil 2 with ratio soil 2:soil 1 of 2.5 (Pi Olsen) and 1.5 (AEM), making 66% differences between both techniques.

The measure of the shoot fresh weight illustrated the correlation between plant growth and the Pi content: plants growing in soil 2 and 3, exhibited 3 and 5 times more biomass than the plants growing in soil 1 (Fig. 6b). Due to the difficulty to assess the roots in the field (importance of the soil adherence, Hinsinger et al., 2011); we concentrated our analyses on the aerial part. The shoot Pi content (Fig. 6c) failed to provide useful indication on the Pi homeostasis status in these plants. In contrast, all molecular markers perfectly played the expected role. The four markers tested (*TaSPX1*, *TaIPS1*, *TaPT4* and *TaGDPD*) were all expressed at the highest level in soil 1 with closer values compared to our control raised in 0 µM Pi artificial substrate, confirming the
<table><thead><tr><th rowspan="2"></th><th colspan="3">Pi Olsen</th><th colspan="3">Anionic Exchange Membrane (24h)</th></tr><tr><th>mg Pi/kg soil</th><th>µmol Pi/Kg soil</th><th>mg P/Kg soil</th><th>mg Pi/kg soil</th><th>µmol Pi/Kg soil</th><th>mg P/Kg soil</th></tr></thead><tbody><tr><td>soil 1</td><td>11.97 ± 0.98</td><td>126.09 ± 10.35</td><td>3.90 ± 0.32</td><td>0.04 ± 0.01</td><td>0.45 ± 0.12</td><td>0.013 ± 0.003</td></tr><tr><td>soil 2</td><td>47.02 ± 8.90</td><td>495.09 ± 93.72</td><td>15.33 ± 2.90</td><td>0.06 ± 0.01</td><td>0.67 ± 0.13</td><td>0.020 ± 0.003</td></tr><tr><td>soil 3</td><td>120.34 ± 26.77</td><td>1267.11 ± 281.85</td><td>39.23 ± 8.73</td><td>0.17 ± 0.08</td><td>1.77 ± 0.80</td><td>0.055 ± 0.026</td></tr></tbody></table>

<figure><img src="image_32.png" /><figcaption>**Fig. 6** Growth of wheat plants grown during 30 days in three different soils. **a** Pi content in soils after Olsen extraction (Pi Olsen) or after Pi diffusion in anion exchange membranes (RAE) n = 2–3. **b** Shoot fresh weight (g). **c** Free-Pi content in shoots (ng Pi/mg fresh weight). **d** Relative gene expression analyses of selected molecular markers. TaACT and TaGAPDH were used as housekeeping genes. 0 design plants grown in silica:vermiculite (1:1) with Pi supply of 0 µM. Bars indicate means ± SD. Different letters indicate significantly different means (one-way ANOVA followed by SNK test, *p* < 0.05, *n* = 6–7)</figcaption></figure>

extreme poverty of Pi present in this soil sample. All the molecular markers also exhibited the highest reduction of their expression in soil 3 and an intermediary score for samples resulting from soil 2 was obtained. The dynamic of response was found in a range of 27-fold (TaIPS1) to 6.5-fold (TaSPX1) except TaPT4 that was not expressed in samples deriving from soil 3 (and vary 8 times between soil 1 and soil 2). This clearly illustrated the capacities of the selected molecular markers to report accurately the Pi status of the samples.

To confirm the validity of these markers in open field conditions, a trial was launched in wheat (in the south of France) by using different fertilization strategies

(-/+fertilization, 120 mg/kg ha P₂O₅ Olsen). This trial was performed on a silty acid soil presenting a 20–38 mg/Kg P₂O₅ Olsen (corresponding to the intermediate situation between soil 1 (12 mg/kg P₂O₅ Olsen) and soil 2 (47 mg/kg P₂O₅ Olsen) used in Fig. 6). Two commercial wheat varieties presenting different sensitivities to Pi deficiency (cv. Johnson, more sensitive; and cv. Fluor, more resistant) were used. Drone imaging was used to measure the impact of Pi fertilization on live green biomass (Vegetation Index, VI) and determination of shoot Pi content and gene expression was carried out as previously described. As already observed for soils from Fig. 6, no significant differences were obtained for total P and
free-Pi content (Fig. S3A and B). Vegetative index indicated a positive impact of fertilization, being the difference higher for *cv. Johnson* (+70%)— more sensitive to Pi deficiency, than for *cv. Fluor* (+10%) (Fig. S3C). Molecular probes tested (GDPD, SPX1, IPS1) were able to discriminate both regimes in both varieties in a significant way (Fig. S3D). Interestingly whereas *Johnson* and *Fluor* produced similar yield in Pi depleted conditions (36,9 qx/ha), the results differ significantly in Pi fertilized field with 65,4 and 56,3 qx/ha respectively. This illustrated the complexity of Pi use efficiency traits. In *Fluor*, Pi fertilization promoting milder effect (+14% vegetative growth and +52% yield) and yield production (+52%) than for *Johnson* where vegetative growth and yield increased +75% and +76% respectively.

## Discussion

### Pi deficiency significantly impacts biomass production

Pi, being a main macronutrient, strongly impacts growth when it is present in limiting amount in the substrate. The impact of discrepancies between treatments on plant growth increase with time. Therefore, a plant with a short life cycle such as *Arabidopsis* is expected to be more impacted after a short period. Indeed, after only two weeks, a 50-fold difference between plants growing in Pi-depleted versus rich media can be observed. All the chosen crops have longer life cycle (range from 6 to 8 months on average) than *Arabidopsis* and were therefore studied after a slightly longer period of growth (one month). Interestingly, they all present noticeable difference of growth ranging from 7- to tenfold respectively for wheat and rapeseed, but only twofold for maize. For maize, the limited impact is probably a consequence of the important capacity of maize to stock phytate, a storage form of phosphorus. Indeed, a maize seed can contain 0.7 to 2.2% of phytate, whereas it usually do not exceed half of this concentration in wheat [64]. Besides, for most small-grained cereals 90% of phytate is in the aleurone whereas in maize 90% of the phytate is located in the scutellum and therefore this could fully benefit to young plantlets development [65]. This probably explains why most significant impact of Pi deficiency are observed at latter stages in the maize (between the synthesis of leaves 7 to 17 (Plenet et al., 2000) whereas we reported the observation at stage of 5 to 6 leaves).

### The analysis of internal Pi (or P) content imperfectly reflects Pi homeostasis status

Using measurement of Pi content presents a major drawback. Pi starts to accumulate when it is present in large excess, and this is often not observed in nature. The three soils tested for wheat growth illustrated well this point with a six-fold difference of shoot biomass between the

poorest and the richest soils (soil 1 and soil 3 respectively, Fig. 6b) whereas Pi content was statistically not different between samples (Fig. 6c). A similar situation was observed with the artificial substrate where Pi accumulated in wheat starting from 250 µM (Fig. 4c) when maximum biomass was already present (Fig. 4b). Identical conclusions can be raised with the other species, where significant growth modifications often occur compared to Pi-starved plants without modification of Pi content (as illustrated, for example, in *Arabidopsis* plants growing with 50 µM Pi (Fig. 1) or rapeseed growing with 250 µM Pi (Fig. 2)). Therefore, plant growth cannot be directly correlated with internal free-Pi content. Besides, even if the range of response for Pi accumulation (around eight-fold) could appear satisfactory, we should keep in mind that this reflects only an extreme situation in our analysis, rarely encountered in nature. This is illustrated by the experiment performed with soil 3 (Pi-rich soil). If we use parameters such as growth or expression of molecular markers, we can see that this soil provides conditions mimicked by the addition of 250/500 µM Pi in the artificial substrate. These conditions, which remained far from extreme point (1500 µM Pi), were the ones where maximum accumulation of Pi was observed. Such observation is fully confirmed by our field trial. It shows that Pi fertilization clearly improves wheat growth (Vegetation Index, VI) and triggers clear response of the molecular markers but neither Pi nor P content turn out to be discriminant in such conditions (Fig. S3).

### Limitation of the detection of anthocyanin presence

Anthocyanins are important molecules protecting the plant against UV radiation and therefore preventing damages resulting from high light. They are also metal chelating agents, limiting toxic effects of these compounds on photosynthesis [66]. Metals are well known to induce these components and their accumulation is promoted by Pi starvation. Indeed, Pi is a strong cation chelator in soil which reduce bioavailability of metals and therefore Pi starvation conditions are well known to favor metal accumulation in plants [10]. Nevertheless, if this may contribute to anthocyanin presence in Pi-deficient plant, a direct link also exists between Pi deficiency and anthocyanin accumulation as *phr1* deficient *Arabidopsis* plant exhibited a strong reduction of anthocyanin accumulation when grown in Pi-deficient medium [63]. They provide a clear indicator for *Arabidopsis* [2, 67] as symptoms of Pi deficiency take place at relatively early stages (between first and second week of growth, [2]). Nevertheless, with the crops used here, it turns out to be much less obvious to use this criterion as anthocyanin accumulation was not exceeding two-fold in the best cases (for rapeseed leaves) or was located in very specific tissues
(bottom of maize stem). Besides, for such last case, it should be noticed that anthocyanins can be observed at all concentrations, and this is only the extend of the area which is modified (Fig. 2a). It is important to notice that anthocyanins were specifically located and did not spray all over the plants as measurement on leaves did not provide any useful clue (Fig. S2b).

More importantly, the major drawback of using the anthocyanins as a Pi-deficient indicator is their lack of specificity, as many stresses promoted their accumulation such as sugar accumulation, nutrient deficiency, metals, cold treatment and high light [66].

## Interest of the molecular markers to detect Pi deficiency

The absorption of Pi through the root is a very rapid phenomenon [68]. It triggers transcriptional modifications for the genes regulated by PHR1 starting from 3 to 5 min following the addition of Pi as recently shown in *Arabidopsis* root [43]. The modifications of transcription in the leaves appear slightly delayed by 30 min in *Arabidopsis* [43]. Pi translocation measured by detection of radiotracers ($^{32}$Pi or $^{33}$Pi) revealed a process taking place within a few minutes (even for plants exhibiting more important development than *Arabidopsis* such as soybean [68]). It should be noticed that the Pi distribution favors the young leaves due to their sink status before reaching the older ones. Such phenomenon is also observed during Pi translocation from old to young active tissues [69, 70] promoting an important heterogeneity of Pi distribution in the aerial part, which is well correlated with the modifications of transcription for selected Pi homeostasis markers. This is a very important parameter, because if Pi is a systemic component for the plant, it indeed triggers a response which is far to be homogenous. Therefore, to limit the bias, it is crucial (i) to avoid collecting small samples and (ii) to discard senescent leaves (exhibiting Pi starved status). From a practical point of view, we advise to favor the harvest of the young leaves (Fig. 5).

## Conclusions

Although numerous reporter genes based on transcriptional or translational fusion [71] have been developed to detect Pi deficiency [60, 63, 72–75], they all required the creation of transgenic plants. This constraint severely limits their use outside the laboratory unlike the markers identified here, which can be used in a wide range of plants as soon as we have access to their genome. In this work, all selected molecular markers turn out to provide excellent and reliable tools to monitor Pi deficiency status. These genes are selected as very early markers due to their short half-life and they are all direct targets of PHR1 [44]. The investigation of other important biotic

and abiotic stresses (Fig. S1) fails to identify conditions triggering their induction, highlighting their specificity of response associated to Pi deficiency. Even if we cannot rule out the possibility that some of them responds to other signals, the use of 4 to 5 independent markers should easily help to discard false positive. Such problem is illustrated on Fig. 2 with BnSPX1 control (0 mM Pi) which exhibited significant but limited induction compared to 50 or 100 mM Pi conditions. This may be due in such conditions (strong Pi deficiency) to the existence of crosstalks with other metabolism such as nitrogen, which is known to affect the members of SPX multigenic family [76, 77].

All markers selected are conserved among plant species. Excepted *AtUNICORN1* (At5g20790, [43]), which turned out to be specific from dicotyledonous species and could not be identified in maize and wheat. As most markers belong to multigenic families we advise to test the different homologues to identify those exhibiting the best specificity and dynamic of response for Pi.

These markers can be used in both leaves and roots, but it should be kept in mind that Pi transporters are expressed at much lower level in the aerial tissues, as previously reported [73, 78]. Therefore, for members of this last gene family, assay may be necessary to identify those exhibiting the highest level of expression in the leaves.

Main signal triggering Pi homeostasis regulation has been recently identified as a Pi metabolite named Inositol 8 pyrophospate [79, 80]. Unfortunately, such labile compound is very difficult to quantify. Robust markers such as the ones identified here provide a very interesting solution to access plant Pi status, as they are extremely sensitive in the low Pi range, which represent almost 30% of the cropland area on earth [8]. Therefore, they represent a valuable tool for physiological studies or breeders to take into account of Pi bioavailability impact on plant growth to cope with a rising economical and societal problem due to shortage and phyto-pollution trouble encountered with this critical resource.

## Abbreviations

<table><tr><td>AEM</td><td>Anionic Exchange Membrane</td></tr><tr><td>At</td><td><i>Arabidopsis thaliana</i></td></tr><tr><td>Bn</td><td><i>Brassica napus</i></td></tr><tr><td>ATP</td><td>Adenosine tri phosphate</td></tr><tr><td>cDNA</td><td>Complementary deoxyribonucleic acid</td></tr><tr><td>EU</td><td>European union</td></tr><tr><td>FAO</td><td>Food and Agriculture Organization of the United Nations</td></tr><tr><td>IPS1</td><td>Induced phosphate starvation 1</td></tr><tr><td>NDVI</td><td>Normalized Difference Vegetative Index</td></tr><tr><td>gDNA</td><td>Genomic deoxyribonucleic acid</td></tr><tr><td>P</td><td>Total phosphorus</td></tr><tr><td>P-</td><td>Field trial with no application of Pi fertilizer</td></tr><tr><td>PAP</td><td>Purple acid phosphatase</td></tr></table>
Pi Inorganic phosphate

PHT1 High affinity phosphate transporter 1

PHR1 Phosphate response 1

P1BS Binding site of transcription factor PHR1

RNA Ribonucleic acid

RNAseq Ribonucleic acid sequencing

RT-qPCR Reverse transcriptase polymerase chain reaction

SPX SPX domain protein

Ta Triticum aestivum

Zm Zea mays

## Supplementary Information

The online version contains supplementary material available at https://doi.org/10.1186/s12870-023-04411-8.

### Additional file 1.

### Acknowledgements

We specially thank Maëlle CABANNE and Marine FOYER for their excellent technical support. We sincerely thank Dr. Nusrat ALI for providing the house-keeping genes for the three plants crops used in this study. We would also like to thank Dr. Frank JAMOIS as Director of the Research Analytical Platform at Agro Innovation International-Timac Agro. We would like to thank our colleagues Eva DESCHAMPS and Yann BRANDT in charge of the coordination, piloting and realization of the field trial and the trainees and interims who participated in the work at Arvalis. We would like to dedicate this article to the memory of Damien De Craieye. He actively participated in this work during his master and despite a very courageous struggle was taken by a devastating illness this summer.

### Authors' contributions

LC, PD, CP, SP and LN: experimental design; LC, PD, DC and SN: performing the laboratory experiments; LC, PD, MA, SP and LN: data analyses; MF, DH and FD: field experiments; LC, PD, SN and LN: writing.

### Funding

This project (Phyto-P) was supported by DRF impulsion program from the CEA. Part of this work was conducted under the ANR Forage and Ultim, supported by the French National Research Agency (ANR 18-CE43 et ANR 21-CE06).

### Availability of data and materials

The data that supports the findings of this study are available in the supplementary material of this article.

### Declarations

### Ethics approval and consent to participate

Field and laboratory studies were conducted in accordance with local legislation. This article does not contain any studies with human participants or animals and did not involve any endangered or protected species. The Arabidopsis WT (Col-0) plants were obtained from ABRC. The deposition number is CS70000. The permission of using these materials in our experiments were obtained. The plant materials sampled and experiments performed in this research complied with institutional, national, and international guidelines and legislation.

### Consent for publication

Not applicable.

### Competing interests

The authors declare no competing interests.

Received: 23 May 2023 Accepted: 14 August 2023

Published online: 23 August 2023

### References

1. Jouhet J, Marechal E, Block MA. Glycerolipid transfer for the building of membranes in plant cells. Prog Lipid Res. 2007;46(1):37-55.

2. Misson J, Raghothama KG, Jain A, Jouhet J, Block MA, Bligny R, Ortet P, Creff A, Somerville S, Rolland N, et al. A genome-wide transcriptional analysis using Arabidopsis thaliana Affymetrix gene chips determined plant responses to phosphate deprivation. Proc Natl Acad Sci U S A. 2005;102(33):11934-9.

3. Poirier Y, Bucher M. Phosphate transport and homeostasis in Arabidopsis. In: Somerville CR, E.M. M, editor. The Arabidopsis Book. Rockville, MD: American Society of Plant Biologists; 2002. p. 1-35.

4. Vance CP, Uhde-Stone C, Allan D. Phosphorus acquisition and use: critical adaptations by plants for securing a nonrenewable resource. New Phytol. 2003;157(3):423-47.

5. Peret B, Clement M, Nussaume L, Desnos T. Root developmental adaptation to phosphate starvation: better safe than sorry. Trends Plant Sci. 2011;16(8):442-50.

6. Chen J, Wang Y, Wang F, Yang J, Gao M, Li C, Liu Y, Liu Y, Yamaji N, Ma JF, et al. The rice CK2 kinase regulates trafficking of phosphate transporters in response to phosphate levels. Plant Cell. 2015;27(3):711-23.

7. Alewell C, Ringeval B, Ballabio C, Robinson DA, Panagos P, Borrelli P. Global phosphorus shortage will be aggravated by soil erosion. Nat Commun. 2020;11(1):4546.

8. MacDonald GK, Bennett EM, Potter PA, Ramankutty N. Agronomic phosphorus imbalances across the world's croplands. Proc Natl Acad Sci U S A. 2011;108(7):3086-91.

9. Holford ICR. Soil phosphorus: Its measurement, and its uptake by plants. Aust J Soil Res. 1997;35(2):227-39.

10. Hirsch J, Marin E, Floriani M, Chiarenza S, Richaud P, Nussaume L, Thibaud MC. Phosphate deficiency promotes modification of iron distribution in Arabidopsis plants. Biochimie. 2006;88(11):1767-71.

11. Schachtman DP, Reid RJ, Ayling SM. Phosphorus uptake by plants: from soil to cell. Plant Physiol. 1998;116(2):447-53.

12. Shen J, Yuan L, Zhang J, Li H, Bai Z, Chen X, Zhang W, Zhang F. Phosphorus dynamics: from soil to plant. Plant Physiol. 2011;156(3):997-1005.

13. Tiessen H. Phosphorus in the global environment. In: White PJH JP, editor. The Ecophysiology of Plant-Phosphorus Interactions. Dordrecht: Springer Netherlands; 2008. p. 1-7.

14. Simpson RJ, Oberson A, Culvenor RA, Ryan MH, Veneklaas EJ, Lambers H, Lynch JP, Ryan PR, Delhaize E, Smith FA, et al. Strategies and agronomic interventions to improve the phosphorus-use efficiency of farming systems. Plant Soil. 2011;349(1-2):89-120.

15. Raymond NS, Gomez-Munoz B, van der Bom FJ, Nybroe O, Jensen LS, Muller-Stover DS, Oberson A, Richardson AE. Phosphate-solubilising microorganisms for improved crop productivity: a critical assessment. New Phytol. 2021;229(3):1268-77.

16. Evenson RE, Gollin D. Assessing the impact of the green revolution, 1960 to 2000. Science. 2003;300(5620):758-62.

17. Khush GS. Green revolution: preparing for the 21st century. Genome. 1999;42(4):646-55.

18. Pingali PL. Green revolution: impacts, limits, and the path ahead. Proc Natl Acad Sci U S A. 2012;109(31):12302-8.

19. Hunter MC, Smith RG, Schipanski ME, Atwood LW, Mortensen DA. Agriculture in 2050: recalibrating targets for sustainable intensification. Bioscience. 2017;67(4):385-90.

20. Kopittke PM, Menzies NW, Wang P, McKenna BA, Lombi E. Soil and the intensification of agriculture for global food security. Environ Int. 2019;132:105078.

21. FAO. Global agriculture towards 2050. In: How to feed the world 2050: High level expert forum: 12-13 Oct 2009. Rome: Office of the Director, Agricultural Development Economics Division Economic and Social Development Department; 2009.

22. Molotoks A, Stehfest E, Doelman J, Albanito F, Fitton N, Dawson TP, Smith P. Global projections of future cropland expansion to 2050 and direct impacts on biodiversity and carbon storage. Glob Chang Biol. 2018;24(12):5895-908.

23. Eitelberg DA, van Vliet J, Verburg PH. A review of global potentially available cropland estimates and their consequences for model-based assessments. Glob Chang Biol. 2015;21(3):1236-48.

24. Etesami H, Jeong BR. Silicon (Si): Review and future prospects on the action mechanisms in alleviating biotic and abiotic stresses in plants. Ecotoxicol Environ Saf. 2018;147:881-96.
25. Glick BR. Bacteria with ACC deaminase can promote plant growth and help to feed the world. Microbiol Res. 2014;169(1):30-9.

26. Hopkins BG, Hansen NC. Phosphorus management in high-yield systems. J Environ Qual. 2019;48(5):1265-80.

27. Bennett EM, Carpenter SR, Caraco NF. Human impact on erodable phosphorus and eutrophication: a global perspective: increasing accumulation of phosphorus in soil threatens rivers, lakes, and coastal oceans with eutrophication. Biosci Biotechnol Biochem. 2001;51(3):227-34.

28. Good AG, Beatty PH. Fertilizing nature: a tragedy of excess in the commons. PLoS Biol. 2011;9(8):e1001124.

29. Pinckney JL, Paeri HW, Tester P, Richardson TL. The role of nutrient loading and eutrophication in estuarine ecology. Environ Health Perspect. 2001;109(Suppl 5):699-706.

30. Food and Agriculture Organization of the United Nations. Soil Resources DaCS: Effects of intensive fertilizer use on the human environment. In: FAO soils bulletin. vol. 16. FAO soils bulletin. The Hague: FAO; 1972.

31. Casacuberta N, Masque P, Garcia-Orellana J. Fluxes of (238)U decay series radionuclides in a dicalcium phosphate industrial plant. J Hazard Mater. 2011;190(1-3):245-52.

32. da Conceicao FT, Antunes ML, Durrant SF. Radionuclide concentrations in raw and purified phosphoric acids from Brazil and their processing wastes: implications for radiation exposures. Environ Geochem Health. 2012;34(1):103-11.

33. Othman I, Al-Masri MS. Impact of phosphate industry on the environment: a case study. Appl Radiat Isot. 2007;65(1):131-41.

34. Ulrich AE. Cadmium governance in Europe's phosphate fertilizers: not so fast? Sci Total Environ. 2019;650:541-5.

35. Gilbert N. Environment: the disappearing nutrient. Nature. 2009;461(7265):716-8.

36. de Ridder M, de Jong S, Polchar J, Lingemann S: Risks and opportunities in the global phosphate rock market: robust strategies in times of uncertainty. In: Strategic Studies. vol. 17/12/12: The Hague Centre for Strategic Studies (HCSS); 2012.

37. Schoumans OF, Bouraoui F, Kabbe C, Oenema O, van Dijk KC. Phosphorus management in Europe in a changing world. Ambio. 2015;44:S180-92.

38. Cordell D, Jackson M, White S. Phosphorus flows through the Australian food system: Identifying intervention points as a roadmap to phosphorus security. Environ Sci Policy. 2013;29:87-102.

39. Nawara S, van Dael T, De Cooman E, Elsen A, Merckx R, Smolders E, Amery F. Testing soil phosphorus in a depleting P scenario: an accelerated soil mining experiment. Eur J Soil Sci. 2018;69(5):804-15.

40. Kalaji HM, Baba W, Gediga K, Goltsev V, Samborska IA, Cetner MD, Dimitrova S, Piszcz U, Bielecki K, Karmowska K, et al. Chlorophyll fluorescence as a tool for nutrient status identification in rapeseed plants. Photosynth Res. 2018;136(3):329-43.

41. Kanno S, Cuyas L, Javot H, Bligny R, Gout E, Dartevelle T, Hanchi M, Nakanishi TM, Thibaud MC, Nussaume L. Performance and limitations of phosphate quantification: guidelines for plant biologists. Plant Cell Physiol. 2016;57(4):690-706.

42. Grennan AK. Phosphate accumulation in plants: signaling. Plant Physiol. 2008;148(1):3-5.

43. Hani S, Cuyas L, David P, Secco D, Whelan J, Thibaud MC, Merret R, Mueller F, Pochon N, Javot H, et al. Live single-cell transcriptional dynamics via RNA labelling during the phosphate response in plants. Nat Plants. 2021;7(8):1050-64.

44. Bustos R, Castrillo G, Linhares F, Puga MI, Rubio V, Perez-Perez J, Solano R, Leyva A, Paz-Ares J. A central regulatory system largely controls transcriptional activation and repression responses to phosphate starvation in Arabidopsis. PLoS Genet. 2010;6(9):e1001102.

45. Calderon-Vazquez C, Ibarra-Laclette E, Caballero-Perez J, Herrera-Estrella L. Transcript profiling of Zea mays roots reveals gene responses to phosphate deficiency at the plant- and species-specific levels. J Exp Bot. 2008;59(9):2479-97.

46. Secco D, Wang C, Shou H, Schultz MD, Chiarenza S, Nussaume L, Ecker JR, Whelan J, Lister R. Stress induced gene expression drives transient DNA methylation changes at adjacent repetitive elements. eLife. 2015;4:e09343.

47. Thibaud MC, Arrighi JF, Bayle V, Chiarenza S, Creff A, Bustos R, Paz-Ares J, Poirier Y, Nussaume L. Dissection of local and systemic transcriptional responses to phosphate starvation in Arabidopsis. Plant J. 2010;64(5):775-89.

48. Secco D, Jabnoune M, Walker H, Shou H, Wu P, Poirier Y, Whelan J. Spatiotemporal transcript profiling of rice roots and shoots in response to phosphate starvation and recovery. Plant Cell. 2013;25(11):4285-304.

49. Sarrobert C, Thibaud MC, Contard-David P, Gineste S, Bechtold N, Robaglia C, Nussaume L. Identification of an Arabidopsis thaliana mutant accumulating threonine resulting from mutation in a new dihydrodipicolinate synthase gene. Plant J. 2000;24(3):357-67.

50. Godon C, Mercier C, Wang X, David P, Richaud P, Nussaume L, Liu D, Desnos T. Under phosphate starvation conditions, Fe and Al trigger accumulation of the transcription factor STOP1 in the nucleus of Arabidopsis root cells. Plant J. 2019;99(5):937-49.

51. Shin H, Shin HS, Dewbre GR, Harrison MJ. Phosphate transport in Arabidopsis: Pht1;1 and Pht1;4 play a major role in phosphate acquisition from both low- and high-phosphate environments. Plant J. 2004;39(4):629-42.

52. Muchhal US, Pardo JM, Raghothama KG. Phosphate transporters from the higher plant Arabidopsis thaliana. Proc Natl Acad Sci U S A. 1996;93(19):10519-23.

53. Nussaume L, Kanno S, Javot H, Marin E, Pochon N, Ayadi A, Nakanishi TM, Thibaud MC. Phosphate import in plants: focus on the PHT1 transporters. Front Plant Sci. 2011;2:83.

54. Morcuende R, Bari R, Gibon Y, Zheng W, Pant BD, Blasing O, Usadel B, Czechowski T, Udvardi MK, Stitt M, et al. Genome-wide reprogramming of metabolism and regulatory networks of Arabidopsis in response to phosphorus. Plant Cell Environ. 2007;30(1):85-112.

55. Muller R, Morant M, Jarmer H, Nilsson L, Nielsen TH. Genome-wide analysis of the Arabidopsis leaf transcriptome reveals interaction of phosphate and sugar metabolism. Plant Physiol. 2007;143(1):156-71.

56. Osorio MB, Ng S, Berkowitz O, De Clercq I, Mao C, Shou H, Whelan J, Jost R. SPX4 acts on PHR1-dependent and -independent regulation of shoot phosphorus status in Arabidopsis. Plant Physiol. 2019;181(1):332-52.

57. Ayadi A, David P, Arrighi JF, Chiarenza S, Thibaud MC, Nussaume L, Marin E. Reducing the genetic redundancy of Arabidopsis PHOSPHATE TRANSPORTER1 transporters to study phosphate uptake and signaling. Plant Physiol. 2015;167(4):1511-26.

58. Puga MI, Mateos I, Charukesi R, Wang Z, Franco-Zorrilla JM, de Lorenzo L, Irigoyen ML, Masiero S, Bustos R, Rodriguez J, et al. SPX1 is a phosphate-dependent inhibitor of PHOSPHATE STARVATION RESPONSE 1 in Arabidopsis. Proc Natl Acad Sci U S A. 2014;111(41):14947-52.

59. Chevalier F, Cuyas L, Jouhet J, Gros VR, Chiarenza S, Secco D, Whelan J, Seddiki K, Block MA, Nussaume L, et al. Interplay between jasmonic acid, phosphate signaling and the regulation of glycerolipid homeostasis in Arabidopsis. Plant Cell Physiol. 2019;60(6):1260-73.

60. Hanchi M, Thibaud MC, Legeret B, Kuwata K, Pochon N, Beisson F, Cao AQ, Cuyas L, David P, Doerner P, et al. The phosphate fast-responsive genes PEPC1 and PPsPase1 affect phosphocholine and phosphoethanolamine content. Plant Physiol. 2018;176(4):2943-62.

61. Franco-Zorrilla JM, Valli A, Todesco M, Mateos I, Puga MI, Rubio-Somoza I, Leyva A, Weigel D, Garcia JA, Paz-Ares J. Target mimicry provides a new mechanism for regulation of microRNA activity. Nat Genet. 2007;39(8):1033-7.

62. Svistoonoff S, Creff A, Reymond M, Sigoillot-Claude C, Ricaud L, Blanchet A, Nussaume L, Desnos T. Root tip contact with low-phosphate media reprograms plant root architecture. Nat Genet. 2007;39(6):792-6.

63. Rubio V, Linhares F, Solano R, Martin AC, Iglesias J, Leyva A, Paz-Ares J. A conserved MYB transcription factor involved in phosphate starvation signaling both in vascular plants and in unicellular algae. Genes Dev. 2001;15(16):2122-33.

64. Schlemmer U, Frolich W, Prieto RM, Grases F. Phytate in foods and significance for humans: food sources, intake, processing, bioavailability, protective role and analysis. Mol Nutr Food Res. 2009;53(Suppl 2):S330-375.

65. Brinch-Pedersen H, Sorensen LD, Holm PB. Engineering crop plants: getting a handle on phosphate. Trends Plant Sci. 2002;7(3):118-25.

66. Landi M, Tattini M, Gould KS. Multiple functional roles of anthocyanins in plant-environment interactions. Environ Exp Bot. 2015;119:4-17.

67. Arnaud C, Clement M, Thibaud MC, Javot H, Chiarenza S, Delannoy E, Revol J, Soreau P, Balzergue S, Block MA, et al. Identification of phosphatin, a drug alleviating phosphate starvation responses in Arabidopsis. Plant Physiol. 2014;166(3):1479-91.

68. Kanno S, Arrighi JF, Chiarenza S, Bayle V, Berthome R, Peret B, Javot H, Delannoy E, Marin E, Nakanishi TM, et al. A novel role for the root cap in phosphate uptake and homeostasis. eLife. 2016;5:e14577.
69. Maillard A, Diquelou S, Billard V, Laine P, Garnica M, Prudent M, Garcia-Mina JM, Yvin JC, Ourry A. Leaf mineral nutrient remobilization during leaf senescence and modulation by nutrient deficiency. Front Plant Sci. 2015;6:317.

70. Raghothama KG. Phosphate Acquisition. Annu Rev Plant Physiol Plant Mol Biol. 1999;50:665-93.

71. Klimyuk VI, Nussaume L, Harrison K, Jones JD. Novel GUS expression patterns following transposition of an enhancer trap Ds element in Arabidopsis. Mol Gen Genet. 1995;249(4):357-65.

72. Li YT, Gu MA, Zhang X, Zhang J, Fan HM, Li PP, Li ZF, Xu GH. Engineering a sensitive visual-tracking reporter system for real-time monitoring phosphorus deficiency in tobacco. Plant Biotechnol J. 2014;12(6):674-84.

73. Misson J, Thibaud MC, Bechtold N, Raghothama K, Nussaume L. Transcriptional regulation and functional properties of Arabidopsis Pht1;4, a high affinity transporter contributing greatly to phosphate uptake in phosphate deprived plants. Plant Mol Biol. 2004;55(5):727-41.

74. Wang L, Li Z, Qian W, Guo W, Gao X, Huang L, Wang H, Zhu H, Wu JW, Wang D, et al. The Arabidopsis purple acid phosphatase AtPAP10 is predominantly associated with the root surface and plays an important role in plant tolerance to phosphate limitation. Plant Physiol. 2011;157(3):1283-99.

75. Hirsch J, Misson J, Crisp PA, David P, Bayle V, Estavillo GM, Javot H, Chirenza S, Mallory AC, Maizel A, 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(2):e16724.

76. Medici A, Szponarski W, Dangeville P, Safi A, Dissanayake IM, Saenchai C, Emanuel A, Rubio V, Lacombe B, Ruffel S, et al. Identification of molecular integrators shows that nitrogen actively controls the phosphate starvation response in plants. Plant Cell. 2019;31(5):1171-84.

77. Ueda Y, Kiba T, Yanagisawa S. Nitrate-inducible NIGT1 proteins modulate phosphate uptake and starvation signalling via transcriptional regulation of SPX genes. Plant J. 2020;102(3):448-66.

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

79. Dong J, Ma G, Sui L, Wei M, Satheesh V, Zhang R, Ge S, Li J, Zhang TE, Wittwer C, et al. Inositol pyrophosphate InsP8 acts as an intracellular phosphate signal in Arabidopsis. Mol Plant. 2019;12(11):1463-73.

80. Zhu J, Lau K, Puschmann R, Harmel RK, Zhang Y, Pries V, Gaugler P, Broger L, Dutta AK, Jessen HJ, et al. Two bifunctional inositol pyrophosphate kinases/phosphatases control plant phosphate homeostasis. eLife. 2019;8:e43582.

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

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