Diurnal Changes in Transcript and Metabolite Levels during the Iron Deficiency Response of Rice

Background Rice (Oryza sativa L.) is highly susceptible to iron (Fe) deficiency due to low secretion levels of the mugineic acid (MA) family phytosiderophore (PS) 2′-deoxymugineic acid (DMA) into the rhizosphere. The low levels of DMA secreted by rice have proved challenging to measure and, therefore, the pattern of DMA secretion under Fe deficiency has been less extensively studied relative to other graminaceous monocot species that secrete high levels of PS, such as barley (Hordeum vulgare L.). Results Gene expression and metabolite analyses were used to characterise diurnal changes occurring during the Fe deficiency response of rice. Iron deficiency inducible genes involved in root DMA biosynthesis and secretion followed a diurnal pattern with peak induction occurring 3–5 h after the onset of light; a result consistent with that of other Strategy II plant species such as barley and wheat. Furthermore, triple quadrupole mass spectrometry identified 3–5 h after the onset of light as peak time of DMA secretion from Fe-deficient rice roots. Metabolite profiling identified accumulation of amines associated with metal chelation, metal translocation and plant oxidative stress responses occurring with peak induction 10–12 h after the onset of light. Conclusion The results of this study confirmed that rice shares a similar peak time of Fe deficiency associated induction of DMA secretion compared to other Strategy II plant species but has less prominent daily fluctuations of DMA secretion. It also revealed metabolic changes associated with the remediation of Fe deficiency and mitigation of damage from resulting stress in rice roots. This study complements previous studies on the genetic changes in response to Fe deficiency in rice and constitutes an important advance towards our understanding of the molecular mechanisms underlying the rice Fe deficiency response. Electronic supplementary material The online version of this article (doi:10.1186/s12284-017-0152-7) contains supplementary material, which is available to authorized users.


Background
Iron (Fe) is utilized in a variety of essential biological processes such as the respiratory (Park et al. 1996) and photosynthetic (Hung et al. 2010) electron-transport chain reactions (Connolly and Guerinot 2002;Reddy 2006), due to its stable redox properties, specifically an ability to donate and accept electrons between both Fe +2 (ferrous) and Fe +3 (ferric) redox states. Although Fe is the fourth most abundant element in the earth's crust, it is not readily available to soil growing plants under aerobic conditions due to its tendency to form insoluble Fe 3+ precipitates (Kim and Guerinot 2007). Insufficient Fe supply to plants leads to Fe depletion of both the mitochondria (Mori et al. 1991) and chloroplast (Stocking 1975;Terry and Low 1982), resulting in impaired efficiency of both cellular respiration and photosynthesis. Additionally, as Fe is required for the formation of the chlorophyll precursor δ-aminolevulinic acid (Marsh et al. 1963;Miller et al. 1984), Fe deficiency results in the characteristic yellowing or chlorosis of the youngest leaves (Terry and Low 1982). Furthermore, Fe deficiency is costly from an agronomic perspective (Hell and Stephan 2003) accounting for more than 25% of the crop yield losses of cereals, legumes and vegetables (Briat et al. 2015).
Iron is also the functional cofactor for the heme containing hemoprotein superfamily of enzymes involved in photosynthesis, respiration, DNA synthesis (Michel and Pistorius 2004) and antioxidants including peroxidases (Zaharieva and Abadía 2003). In cells, antioxidants quench reactive oxygen species (ROS) which are produced as an inevitable consequence of photosynthesis (Parida and Das 2005) and to a lesser extent due to respiration (Gill and Tuteja 2010). Cellular stress induced by Fe deficiency results in impaired antioxidant activity due to absence of the required catalytic Fe, leading to an increased abundance of cellular ROS (Becana et al. 1998;Iturbe-Ormaetxe et al. 1995). Accordingly Fe homeostasis must be tightly regulated in plants to prevent Fe-deficiency-induced oxidative damage to proteins and lipids (Sun et al. 2007).
Plants have evolved mechanisms to increase Fe uptake in response to conditions of limited Fe bioavailability. Dicots and non-graminaceous monocots employ Strategy I uptake, involving proton (H + ) secretion to the rhizosphere and reduction of soil Fe 3+ to Fe 2+ . Secretion of H + reduces soil pH and promotes the solubility of Fe ions and reduction of Fe 3+ to the more bioavailable Fe 2+ promotes increased uptake of Fe 2+ through iron-regulated transporters (IRT) in the plasma membrane of plant root cells Scholz et al. 1992). Graminaceous plant species (grass family Gramineae) utilize Strategy II uptake involving secretion of mugineic acid (MA) family phytosiderophores (PS) into the rhizosphere to chelate and solubilise the soil Fe 3+ , allowing for subsequent uptake of Fe 3+ -PS complexes through yellow stripe 1-like (YSL) transporters in the plasma membrane of plant root cells (Inoue et al. 2009;Nagasaka et al. 2009).
Rice (Oryza sativa L.) is not only the model monocot species but also a staple food that contributes 35-59% of the daily caloric intake for approximately three billion people (Bhullar and Gruissem 2013). Rice is often grown in flooded paddies where microorganisms rapidly deplete dissolved oxygen levels, resulting in anoxic water where Fe has the tendency to be reduced to the more bioavailable Fe 2+ form (Walker and Connolly 2008). Rice has therefore evolved the unusual ability amongst other Strategy II plant species to absorb Fe 2+ via OsIRT transporters in a manner similar to Strategy I plant species (Bughio et al. 2002;Ishimaru et al. 2006). However, whilst genome sequencing demonstrates that the rice genome contains ferric reductase genes (Vasconcelos et al. 2004), rice does not show ferric reductase activity in roots nor the characteristic swelling of roots in response to Fe deficiency as seen in Strategy I plant species (Rong-li et al. 2012). The utilization of aspects from both Strategy I and Strategy II Fe uptake within the rice Fe deficiency response is commonly referred to as the 'Combined Strategy' (Ishimaru et al. 2006;Sperotto et al. 2012). The OsIRT homologous genes HvIRT1 in barley and ZmIRT1 in maize have been reported to be upregulated in response to Fe deficiency. Further, expression of these genes in Fe uptake-defective yeast mutants has been demonstrated to reverse growth defect associated with Fe uptake inability (Li et al. 2013;Pedas et al. 2008). However, in planta uptake of Fe 2+ by Strategy II plants has so far only been confirmed in rice using tracer experiments (Ishimaru et al. 2006;Ricachenevsky and Sperotto 2014).
The amount of PS secreted by Strategy II plants corresponds with plant tolerance to growth on alkaline soils (soils with low Fe bioavailability) as identified by resistance to the onset of Fe-deficiency-associated chlorosis (Marschner et al. 1986). In response to Fe limiting conditions, barley (Hordeum vulgare L.), a plant species able to grow vigorously on high pH soils, secretes large amounts of PS including 2′-deoxymugineic acid (DMA) and the derivatives MA and 3-epihydroxymugineic acid (epiHMA). By contrast, rice, which is highly susceptible to Fe deficiency chlorosis on high pH soils, secretes only DMA in response to Fe deficiency and at much lower levels and shorter durations than barley (Higuchi et al. 1996;Mori et al. 1991). During Fe deficiency, barley (Marschner et al. 1986;Takagi et al. 1984;Walter et al. 1995), wheat (Oburger et al. 2014;Zhang et al. 1991), maize (Ueno et al. 2009), and red fescue (Ma et al. 2003) secrete PS following a diurnal pattern, with higher levels of PS secreted during day (light) periods and lower levels of PS secreted during night (dark) periods.
As Fe is a cofactor in enzymes involved in photosynthesis, this secretion pattern may enable plant Fe uptake to coincide with increased photosynthetic rates during light periods (Takagi et al. 1984). Additionally, as PS are readily degraded by soil microbes as a carbon and/or nitrogen source, the diurnal pattern of PS release may help reduce the degree of PS loss by microbial populations in the rhizosphere (Römheld and Marschner 1990). During Fe deficiency, rice induces expression of DMA biosynthetic genes following a diurnal pattern (Inoue et al. 2009;Nozoye et al. 2004Nozoye et al. , 2011. Additionally, rice secretes higher levels of DMA in the morning than afternoon or night periods (Nozoye et al. 2014). However, the diurnal pattern of transcript and metabolite levels within the rice Fe deficiency response has not been as extensively investigated as the Fe deficiency response of barley.
We studied the diurnal changes of DMA, NA, transcripts and metabolite levels within the rice Fe deficiency response. Rice roots exposed to Fe deficiency were used to investigate diurnal changes in expression levels of genes involved in Strategy II Fe uptake and translocation. A highly sensitive detection method was utilized to quantify in planta NA and DMA levels as well as DMA secretion levels by Fe-deficient rice roots. Additionally, metabolite profiling of Fe-deficient root tissues was performed to understand the molecular mechanisms utilized within the rice Fe deficiency response.

Induction of Fe Deficiency Symptoms
Soil Plant Analysis Development (SPAD) measurements of leaf tissues identified that exposure of 5 week old hydroponically grown rice plants to one week of Fe deficiency induced a significant reduction in the youngest leaf 'greenness' (p ≤ 0.001); however, no significant variation in plant fresh weight was detected (p = 0.972). Accordingly, this treatment was selected as a means of inducing Fe deficiency symptoms in rice without affecting biomass ( Fig. 1).

Biosynthesis and Secretion of DMA in Fe-deficient Rice Roots
Targeted quantification of NA and DMA by liquid chromatography mass spectrometry (LC-MS) from root metabolite extracts collected after one week of Fe deficiency identified that concentrations of DMA (average 1.42-fold higher, p ≤ 0.001), but not NA (average 1.05-fold higher, p = 0.988), were significantly higher in the Fe-deficient roots ( Fig. 3a and b, Additional file 2: Table S2A). Root concentrations of NA and DMA in Fe-sufficient plants showed a similar diurnal pattern, with root NA (p = 0.016) and DMA (p = 0.033) concentrations significantly higher during dark periods than light periods. While NA and DMA concentrations in Fe-deficient roots did not follow a significant diurnal pattern, Fe deficiency induction of root NA and DMA concentrations were both highest during the 11:00-13:00 period (Fig. 3d).
Quantification of DMA by LC-MS from root secretions collected after one week of Fe deficiency, identified that Fe-deficient roots secreted on average 2.4-fold more DMA (p = 0.009) than Fe-sufficient roots across the day (Fig. 3c). Whilst no significant diurnal variation was observed in DMA secretion by Fe-sufficient roots, a diurnal Fig. 1 Induction of Fe deficiency symptoms in rice plants. a Plant fresh weight (n = 9). b SPAD values of the youngest leaf (YL) and 2 nd youngest leaf (2 nd YL) (n = 30). Plant fresh weights and leaf SPAD values of Fe-sufficient (green) and Fe-deficient (red) plants were measured at 0 and 7 days after the Fe deficiency treatment. Values are presented as means ± SE of n biological replicates. Significant differences are indicated by ** for p < 0.01 by Student's t-tests performed between treatments within leaf types and days after Fe deficiency pattern with a peak during 11:00-13:00 (p = 0.001) was observed for DMA secretion by Fe-deficient roots. Accordingly, Fe deficiency promoted increased secretion of DMA following a diurnal pattern with a peak occurring 3-5 h after the onset of light.

Biochemical Pathway Mapping of the Rice Fe Deficiency Response
Concentrations of metabolites in the following pathways were significantly increased in the roots of Fe-deficient rice: α-ketoglutarate family amino acids, 3-phosphoglycerates Fig. 2 Expression of genes associated with DMA biosynthesis in Fe-sufficient and Fe-deficient rice roots. Expression of a OsNAS1, b OsNAS2, c OsNAS3, d OsNAAT1, e OsDMAS1, f OsTOM1 and g OsYSL15 genes in the roots of Fe-sufficient (green) and Fe-deficient (red) rice plants. Data is presented as means of genes of interest (GOI) expression relative to the three gene normalization factor (3GNF) ± SE of 9 biological replicates. Significant differences are indicated by * for p < 0.05, ** for p < 0.01, *** for p ≤ 0.001 by Student's t-tests performed between treatments within collection times (3-P-glycerates), pyruvate derivatives, aromatic amino acids, citrulline-nitric oxide (Cit-NO) cycle, oxaloacetate/aspartate family amino acids and the DMA biosynthetic pathway (Additional file 2: Table S2B). No metabolites were detected at significantly lower concentrations in the roots of the Fe-deficient rice.
Induction factors associated with Fe deficiency were mapped onto custom biochemical pathways (Fig. 5) using the Visualization and Analysis of NeTworks with related Experimental Data program (VANTED) v2.1.0 (Plant Bioinformatics Group, The Leibniz-Institute of Plant Genetics and Crop Plant Research, available at https://immersive-analytics.infotech.monash.edu/vanted/ ). Induced levels of the α-ketoglutarate family amino acids (proline, histidine and glutamine) were observed in Fe-deficient roots following a diurnal pattern with a peak of induction occurring during 18:00-20:00 (Fig. 5a). Induced levels of the 3-P-glycerate derivatives serine and glycine (Fig. 5b) and pyruvate derivatives: valine, leucine, alanine and β-alanine ( Fig. 5c) were also observed in Fe-deficient roots, suggesting an increase in glycolysis during light periods (11:00-13:00 and 18:00-20:00). Increased activity of glycolysis due to Fe deficiency was also supported by a significant increase in concentration of the glucose-to-pyruvate intermediate 3-P-glycerate in Fe-deficient roots (Additional file 3: Figure  S1). Induced levels of the aromatic amino acids: phenylalanine, serotonin, tyrosine, and tryptophan in Fe-deficient roots (Fig. 5d), suggest an increased activity of the shikimate pathway. Induced levels of arginine and citrulline but not ornithine in Fe-deficient roots (Fig. 5e), suggest an increased activity of the citrulline-nitric oxide (Cit-NO) cycle. Additionally, induced levels of the oxaloacetate/aspartate family amino acids: isoleucine, lysine, threonine ( Fig. 5f) and methionine (Fig. 5g) were detected during 6:00-8:00, 11:00-13:00 and 18:00-20:00, but not 22:00-24:00. Activity of the DMA biosynthetic pathway was also increased in Fe-deficient roots based on induced levels of methionine (the dominant precursor of DMA biosynthesis), root DMA and secreted DMA (Fig. 5g). Additionally, the potential methionine precursors serine and cysteine (but not aspartate or homoserine) accumulated in response to Fe deficiency (Additional file 4: Figure S3).

Iron Deficiency Induces Expression of Genes Involved in Strategy II Fe Uptake Following a Diurnal Pattern
Detailed analysis of rice gene expression related to NA and DMA biosynthesis, as well as DMA secretion under Fe deficiency, has been described in several previous investigations (Inoue et al. 2003;Nozoye et al. 2011;Wang et al. 2013). Here we have identified diurnal patterns of expression for the OsNAS1, OsNAS2, OsNAAT1, OsD-MAS1, and OsTOM1 genes in Fe-deficient rice roots with significant induction at the 6:00 and 11:00 sampling times (Fig. 2, Additional file 1: Table S1). These diurnal patterns are consistent with previous investigations (Nozoye et al. 2004(Nozoye et al. , 2011, and may be summarised as a slight induction of expression prior to the onset of light, a peak of induction 3-5 h after the onset of light, and a return to expression levels similar to that of Fe-sufficient roots prior to dark. While the diurnal patterns reported in this study have been previously identified, other studies have reported much higher inductions ratios for Values are presented as means ± SE of n biological replicates. Significant differences are indicated by * for p < 0.05, ** for p < 0.01, *** for p ≤ 0.001 as determined by Student's t-tests performed between treatments within collection times Fe-deficiency-inducible genes including OsNAS1 and OsNAS2 (Inoue et al. 2003;Kobayashi et al. 2009). This may have been due to differences in growing conditions, age of plants and length of exposure to the Fe deficiency treatment. In this study we cultivated rice plants in Fesufficient hydroponic solution for 5 weeks before transferring to Fe-deficient solution for a period of 1 week prior to sampling. Inoue et al. (2003) exposed very young seedlings to Fe-sufficient nutrient solution for 2 weeks, whereas Kobayashi et al. (2009) cultivated seedlings in Fe-sufficient nutrient solution for 3 weeks, before transferring plants to Fe-deficient solution for 1 week. Considering that our plants were more established than these studies prior to Fe deficiency, it is possible that our six week old plants were not fully Fe limited at the time of sampling and this may have reduced the induction ratios that we observed. We also found that expression of OsYSL15 was significantly induced in response to Fe deficiency prior to the onset of  Table S1A). However, the diurnal pattern of OsYSL15 in Fe-deficient roots differed from the diurnal pattern of OsNAS1, OsNAS2, OsNAAT1, OsDMAS1 and OsTOM1 by having a peak time of expression at 18:00-20:00 instead of 11:00-13:00 (Fig. 2). OsYSL15 is an Fe-regulated Fe 3+ -DMA transporter and is thought to be responsible for active uptake of Fe 3+ -DMA complexes during daylight hours (Inoue et al. 2009). Thus, its expression may be induced by the preferential biosynthesis of root DMA at night time (Fig. 3b) and coordinated by diurnal secretion of DMA (Fig. 3c). Whilst increased expression of OsYSL15 after the onset of light has been reported previously (Inoue et al. 2009), the diurnal pattern of OsYSL15 that we observed is not consistent with the daily fluctuations of OsYSL15 described by Inoue et al. (2009). This inconsistency could be explained by differences in rice genotype as well as experimental design. In this study we used rice cultivar Nipponbare whereas Inoue et al. (2009)

Iron Deficiency Promotes DMA Biosynthesis and Accumulation
In this study we observed significant upregulation of OsNAS1, OsNAS2, OsNAAT1 and OsDMAS1 under Fedeficient conditions (Fig. 2), however, only DMA (not NA) concentrations were significantly increased in the Fedeficient rice roots ( Fig. 3a and b); a result that is consistent with previous studies in barley and rice (Higuchi et al. 2001). These results imply that OsNAS upregulation within the Strategy II Fe deficiency response is associated with increased NA production for DMA biosynthesis, rather than NA accumulation. These results also suggest that increased DMA, rather than NA, is utilized for the altered chelation of Fe in the roots of Fe-deficient rice. This is supported by predictive models indicating that at the rice root apoplastic pH of 4.9 (Nishiyama et al. 2008), DMA would outcompete NA as the predominant Fe complex (Hider et al. 2004). Recent speciation studies of metal ligands within chelation complexes in the xylem and phloem sap of rice have identified that DMA is a major chelator of Fe in the phloem sap (Nishiyama et al. 2012;Yoneyama et al. 2015), but neither DMA nor NA are major chelators of Fe in the xylem sap (Ariga et al. 2014). Accordingly DMA-mediated translocation of Fe within the rice Fe deficiency response may occur predominantly via the phloem.
Whilst NA and DMA concentrations in Fe-sufficient roots showed a distinct diurnal pattern with higher concentrations in dark periods than light periods, no significant diurnal pattern of NA or DMA concentrations were observed in Fe-deficient roots. In barley, a diurnal pattern is observed wherein root PS (DMA, MA and epiHMA) concentrations are reduced during periods of heightened PS secretions (Kawai et al. 1993;Nagasaka et al. 2009;Walter et al. 1995). In the current study, Fe deficiency induced a 1.5-fold increase in root DMA concentrations and a 2.4-fold increase in DMA secretions in rice plants (Fig. 3), accounting for root DMA concentrations changing minimally during periods of increased DMA secretion. By contrast in barley, the depletion of root PS concentrations may be explained by Fe deficiency inducing an 8.7-fold increase in root PS concentrations but up to a 250-fold increase in PS secretions (Walter et al. 1995). This difference may relate to the Combined Strategy utilized by rice grown in flooded conditions involving uptake of unbound Fe 2+ via IRT transporters. As this mechanism of Fe uptake is not reliant on DMA secretion, rice may have evolved to retain a larger proportion of the increased DMA produced during Fe deficiency.

Iron Deficiency Promotes Increased Secretion of DMA Following a Diurnal Pattern
Successful modification of the DMA quantification method described in Kakei et al. (2009) from quadrupole time-offlight mass spectrometry to the more sensitive triple quadrupole mass spectrometry provided a 3-fold improvement in sensitivity. Accordingly the developed DMA quantification method was more sensitive than DMA quantification methods described in previous publications including:  (Ma et al. 1995). This improved sensitivity allowed avoidance of the false-negative result wherein DMA within rice root secretions was deemed 'not detected' as the concentration was below the method's limit of detection as in previous publications utilizing nuclear magnetic resonance spectroscopy for analysis (Fan et al. 2001).
Previous studies have identified that rice secretes more DMA during the 6:00-13:00 'morning' period than the 13:00-20:00 'afternoon' or 20:00-6:00 'night' periods (Nozoye et al. 2014). Consistent with this result, our study identified that the peak time of DMA secretions by Fedeficient rice occurred during the 6:00-13:00 'morning' period (Fig. 3). Additionally, the improved sensitivity of the developed method allowed for the analysis of root secretions collected during the 'morning' over two narrower time periods (6:00-8:00 and 11:00-13:00). The ability to analyse root secretion DMA concentrations collected over narrower time periods allowed for further characterisation of rice DMA secretions during this 'morning' period leading to the identification of 11:00-13:00 (3-5 h after the onset of light) as the peak time of DMA secretion from Fe-deficient rice roots (Fig. 3).
Whilst this diurnal pattern of the Fe deficiency response in rice is less distinct than that of other graminaceous plants, the peak time of DMA secretions by rice is consistent with the peak time of PS secretions observed under Fe deficiency in barley (Marschner et al. 1986;Takagi et al. 1984;Walter et al. 1995), wheat (Oburger et al. 2014;Zhang et al. 1991), maize (Ueno et al. 2009), and red fescue (Ma et al. 2003). Consistent with other graminaceous plant species, a major aspect of the rice Fe deficiency response appears to involve re-establishing Fe homeostasis via increased PS mediated Fe uptake and internal translocation with the time of peak secretions occurring 3-5 h after the onset of light.

Iron Deficiency in Rice Induces Changes to Carbon and Energy Metabolism
Previous proteome analyses have identified that the molecular mechanisms induced within plant Fe deficiency responses include changes to plant metabolism associated with: carbon and energy metabolism, S-adenosylmethionine (SAM) metabolism, ion transport, ROS scavenging and cell signalling (Kosová et al. 2011). Proteome analyses of the rice Fe deficiency response specifically identified increased levels of the glycolytic enzyme glyceraldehyde-3-phosphate dehydrogenase in the roots of Fe-deficient rice (Chen et al. 2015). Results of the current study indicated increased glycolysis and respiration within the Fe deficiency response during light periods (11:00-20:00) based on the induced levels of 3-P-glycreate derivatives, branched chain amino acids, pyruvate derivatives (Fig. 5), and 3-P-glycerate (Additional file 3: Figure S1). This is consistent with previous studies on plant stress responses (Di Martino et al. 2003;Wingler et al. 2000;Zoghlami et al. 2011). The rice Fe deficiency response therefore involves changes to carbon and energy metabolism potentially to accommodate the increased energy demand required to drive mechanisms induced within the Fe deficiency response.
Methionine is the dominant precursor of DMA biosynthesis  via the intermediate SAM, which is derived from methionine by SAM synthetase within the methionine cycle ). Transcript analysis of barley suggested that methionine and SAM biosynthesis were induced in response to Fe deficiency (Nagasaka et al. 2009;Negishi et al. 2002). Similarly, transcript analyses of rice suggested that activity of the methionine cycle and DMA biosynthesis are induced 24 h after the onset of Fe deficiency (Itai et al. 2013). In the current study, Fe-deficient rice roots showed induced levels of methionine during 6:00-20:00 (Fig. 5g). However, unlike the diurnal patterns of gene expression detected (Fig. 2), a clear peak of methionine concentrations was not observed during 11:00-13:00. Serine (via cysteine) and/or aspartate (via homoserine) are used by plants to produce homocysteine, the immediate precursor to methionine (Droux et al. 1995). Serine and cysteine levels were induced in Fe-deficient rice roots whilst levels of aspartate and homoserine were unchanged ( Fig. 5b and f; Additional file 4: Figure S4). The co-accumulation of serine, cysteine and methionine suggested that serine may be the dominant precursor to the increased methionine biosynthesis within the rice Fe deficiency response. To confirm this hypothesis, further investigations will be required such as radioisotope tracer experiments previously used to confirm that methionine is the dominant precursor of PS biosynthesis within the Fe deficiency responses of graminaceous plants Nakanishi et al. 1999).

Iron Deficiency in Rice Induces Molecular Mechanisms Associated with Plant Metal-stress Responses
Ion transport, chelation and sequestration are key mechanisms involved in plant metal tolerance and reestablishment of metal homeostasis (Clemens 2001). In response to Fe deficiency, induced levels of histidine, proline and asparagine were detected during 18:00-20:00 ( Fig. 5a and f ). These amines have been identified as chelators of copper, nickel and zinc and facilitate tolerance to these metal stresses via complexation and xylem-mediated translocation (Rai 2002;Sharma and Dietz 2006). Considering that DMA is able to bind a range of metals in the rhizosphere (Scholz et al. 1992), and Fe-deficient rice has previously been shown to uptake increased concentrations of calcium, cadmium and zinc (da Silveira et al. 2007;Zhang et al. 1998), the increased DMA secretions during 11:00-13:00 within the rice Fe deficiency response may promote increased uptake of micronutrients other than Fe. The induced levels of chelators (histidine, proline, asparagine and DMA) during 18:00-20:00 within the rice Fe deficiency response may be associated with increased chelation and translocation of off-target metals (other than Fe) during the 11:00-13:00 period due to increased DMA secretions.

Iron Deficiency in Rice Induces Molecular Mechanisms Involved in Mediation of Oxidative Stress
The increased concentration of ROS from abiotic stresses, including Fe deficiency, results in an "oxidative burst" which in-turn causes damage to cell proteins, lipids, carbohydrates and DNA (Gill and Tuteja 2010;Spinelli et al. 2011). Accordingly, the plant Fe deficiency response includes the induction of molecular mechanisms associated with mediation of the damage due to oxidative stress. Accumulation of amino acids within abiotic stress responses play a protective role as antioxidants which detoxify ROS and prevent damage to the plant cell (Joshi et al. 2010). The antioxidant proline was induced in roots of our Fedeficient rice plants with a peak of induction occurring during 18:00-20:00 (Fig. 5a). This peak time may be associated with off-target metal uptake as the localized increase to metal concentrations may cause oxidative stress via formation of ROS (Ercal et al. 2001). Iron-deficient rice roots also showed induced levels of phenylalanine, tyrosine, serotonin and tryptophan (Fig. 5d) which are antioxidant aromatic amines (Nimalaratne et al. 2011;Ramakrishna and Ravishankar 2011;Stadtman and Levine 2003). Particularly, phenylalanine is the main precursor to a suite of phenolic compounds which have very high antioxidant activity (Ghasemzadeh and Ghasemzadeh 2011;Michalak 2006).
Iron deficiency was associated with a significant decrease in root accumulation of ascorbic acid (Additional file 3: Figure S2), which is typically the primary polar antioxidant in plants (Weber et al. 2008). Previous proteomic studies have identified decreased levels of ascorbate peroxidase in the roots of Fe-deficient rice (Chen et al. 2015). As ascorbate peroxidase is a hemoprotein, Fe deficiency results in a decrease in ascorbate peroxidase activity (Becana et al. 1998;Iturbe-Ormaetxe et al. 1995;Sun et al. 2007). Additionally, previous analyses identified that in rice, but not in other graminaceous plants, Fe deficiency resulted in decreased root activity of glutathione reductase, an enzyme required for the ascorbic acid-glutathione cycle to efficiently mediate oxidative stress (Bashir et al. 2007). Accordingly, the rice Fe deficiency response includes accumulation of non-Fe-dependant antioxidants to scavenge ROS and minimise ROS associated damage.
Nitric oxide (NO) is a highly reactive radical that acts as a chemical signal (Gechev et al. 2006) in multiple plant stress responses including; drought, salt, metal toxicity (Spinelli et al. 2011) and Fe deficiency (Arnaud et al. 2006;Vigani et al. 2013). However, NO may also act as an antioxidant to quench superoxide anion radicals (Arasimowicz and Floryszak-Wieczorek 2007). In this study, Fe deficiency was associated with induced levels of citrulline and arginine (Fig. 5e). These results indicate that the Cit-NO cycle was induced which in turn suggests that NO production was increased (Mori 2007). Treatment of Arabidopsis with putrescine has previously been shown to alleviate Fe deficiency via the accumulation of NO (Zhu et al. 2016). Similarly, treatment with NO directly has been demonstrated to mitigate the chlorosis associated with Fe deficiency in rice ) and maize (Sun et al. 2007) by acting as an antioxidant. Accordingly, increased NO production may occur as part of the rice Fe deficiency response to signal induction of gene expression and/or mitigate damage associated with oxidative stress.

Conclusion
The Strategy II Fe deficiency response of rice has been more difficult to study than that of wheat or barley due to very low levels of DMA secretion. The development of a highly sensitive method for the quantification of DMA and NA in this study allowed for the quantification of DMA in rice root secretions which for many previously published methodologies is below the limit of detection. This study reported a thorough characterization of the rice Fe deficiency response. Gene expression and metabolite profiling identified that Fe deficiency induced DMA biosynthesis and DMA secretions followed a diurnal pattern with a peak of induction occurring during 11:00-13:00. The identification of this peak time period, equating to 3-5 h after the onset of light, verified that the peak time of the rice Fe deficiency response is similar to that of other Strategy II plant species such as barley and wheat. However, the diurnal variation in DMA secretion observed within the rice Fe deficiency response was less prominent than that of other graminaceous species described in previous studies. The results presented herein indicate that the mechanisms utilized by rice in response to Fe deficiency include increased DMA-mediated Fe uptake and accumulation of metabolites associated with plant oxidative stress responses. These results constitute an important progression in our understanding of the rice Fe deficiency response.

Plant Materials and Growth Conditions
Rice (Oryza sativa L. cv. Nipponbare) grains were surface sterilized with 70% ethyl alcohol and 5% sodium hypochlorite and then germinated on filter paper wetted with 7 mL of reverse osmosis (RO) water in a controlled environment chamber with 12-h photoperiod at 24°C. After 10 days, germinated rice seedlings were transferred to 13 L supported hydroponic tanks. The rice seedlings were grown for six weeks at a density of thirty plants per hydroponic tank in a growth room with 12-h photoperiod, and 28°C/24°C day/night temperature (8:00-20:00/20:00-8:00), with nutrient solution as described in Plett et al. (2010). Rice plants were initially acclimated for one week by growth in ¼ strength Fe-replete growth solution followed by ½ strength Fe-replete growth solution for an additional four days (Additional file 5: Table S3). Plants were grown for a further 4 weeks in full strength Fereplete growth solution. Growth solutions were replaced every 3-4 days to maintain pH at 5.5. After the four week growth period, Fe-sufficient control plants were grown for an additional week in full strength Fe-replete growth solution while Fe-deficient treatment plants were grown in growth solution lacking NaFe 3+ EDTA.

Measurement of Plant Size and Chlorophyll Score
Whole plant fresh weight was measured after one week of Fe deficiency treatment. Soil Plant Analysis Development (SPAD) measurements were taken on the youngest leaf (YL) and 2 nd youngest leaf (2 nd YL) at day 0 and 7 of the Fe deficiency treatment using a Minolta SPAD-502Plus (Spectrum Technologies Inc., Plainfield, Illinois, USA) as an indicator of leaf chlorophyll content.

Root Metabolite Extraction
Frozen root tissue was freeze dried using a Labconco Lyph-lock 6 freeze-dry system (Labconco Corp., Kansas City, Missouri, USA), transferred into 2 mL cryo-mill tubes (Lysing Matrix tube, MP Biomedicals, Solon, Ohio, USA), then 200 μL of methanol (MeOH) was added (containing the internal standards: 13 C-Sorbitol at 20 μg/ mL, 13 C 5 -15 N Valine at 20 μg/mL, 2-aminobutyric acid at 10 μg/mL and pentafluorobenzoic acid at 10 μg/mL). The root tissue was then homogenized using a Precellys 24 cryo-mill coupled to a Cryolys (Bertin Technologies, Villeurbanne, France) using the following setting: 3 cycles involving 30 s of milling at 6,800 rpm followed by a 30 s rest period between cycles. The MeOH supernatant was removed and 200 μL of deionized MilliQ water was then added to the cryo-mill tube containing the homogenized root tissue pellet and vortexed for 30 s. The water supernatant was combined with the removed MeOH supernatant and 100 μL of dichloromethane (DCM) was added to the MeOH-water supernatant mixture.

Root Secretion Collection, Purification and Concentration
Root secretions were collected using the method described in Astolfi et al. (2012) with the following modifications: roots from three plants were briefly rinsed with RO water, transferred into 400 mL of autoclaved dH 2 O in 1 L plastic containers, and then wrapped with aluminium foil to prevent exposure to light. After two hours of root secretion collection, 40 μL of Micropur Forte (Katadyn Products Inc.) was immediately added and filtered through Whatman filter discs. Root secretions were purified using cation exchange Amberlite IR120 H resin (Rohm and Haas Company, Philadelphia, Pennsylvania, USA) with a bed volume (BV) of 50 mL. Sample solutions were eluted through the column at two BV per hour, and this eluate was discarded. Recovery of exchangeable cations (including DMA) was performed by eluting 100 mL of 2 M ammonium hydroxide solution through the Amberlite resin at two BV per hour and stored at −80°C. Concentration of root secretions was performed through centrifugal evaporation using a RVC 2-25 (Martin Christ Corporation, Osterode, Germany) at 30°C.

LC-MS Quantification of Amines
Quantification of root metabolite concentrations was performed as described in Boughton et al. (2011) using an LC 1200 series binary pump with autosampler and heated column compartment coupled to a 6410 series Triple Quadrupole MS (Agilent Technologies Inc.). In brief, derivatization of samples involved combining 10 μL of sample, 70 μL of 200 mM borate buffer (containing the internal standard 13 C 5 -15 N-Valine), and 20 μL of fresh 2.85 mg/mL 6-aminoquinolyl-N-hydroxysuccinimidyl carbamate (Aqc) solution. The solution was briefly vortexed, centrifuged at 13,000 rpm for 1 min at 4°C, then incubated in an Eppendorf thermomixer compact (Eppendorf AG) at 1,150 rpm for 10 min at 55°C. Reverse phase (RP) chromatography was performed using a Zorbax Eclipse XDB-C18 Rapid Resolution HS 2.1x100 mm 1.8 Micron column (Agilent Technologies Inc.) with a flow rate of 300 μL/min and a temperature of 30°C. The aqueous mobile phase, solvent ' A' , was prepared as a 0.1% (v/v) formic acid (FA) in dH 2 O, and the organic mobile phase, solvent 'B' , was prepared as 0.1% (v/v) FA in acetonitrile. Chromatography was performed using the following gradient: During MS analyses, all metabolites were quantified based on formation of the m/z = 171 Amq product ion. Metabolites were quantified by multiple reaction monitoring (MRM) according to previously published transitions (Boughton et al. 2011;Callahan et al. 2007).

LC-MS Quantification of DMA and NA
Quantification of 9-fluorenyl methoxycarboxyl chloride (FMOC-Cl) derivatized DMA and NA by RP LC-MS was performed on an LC 1290 series infinity pump with autosampler and heated column compartment coupled to a 6490 series Triple Quadrupole MS (Agilent Technologies Inc.) using a modified version of the method described in Kakei et al. (2009). In short, derivatization of samples involved combining, 5 μL of sample, 10 μL of 1 M sodium borate buffer (pH = 8) containing 42 ng/mL of the internal standard 2-aminoanthracene (2-AA), 10 μL of 50 mM EDTA (pH = 8), and 40 μL of fresh 50 mM FMOC-Cl solution. The solution mixture was incubated at 60°C with constant shaking at 700 rpm, for 15 min then the reaction was quenched via the addition of 8.9 μL of 5% FA solution (pH = 4). A Zorbax Eclipse XDB-C18 Rapid Resolution HS 2.1x100 mm 1.8 Micron column (Agilent Technologies Inc.) was used during chromatography. The aqueous mobile phase, solvent ' A' , was prepared as a 0.1% (v/v) FA in dH 2 O, and the organic mobile phase, solvent 'B' , was prepared as 0.1% (v/v) FA in acetonitrile. Chromatography was performed using the following gradient:

Data Processing and Statistical Analyses
Spectra gathered during LC-MS analyses were initially processed using the MassHunter Workstation software package vB.06 (Agilent Technologies Inc.). For LC-MS analyses, compounds were identified by retention time and MRM match to an authentic standard, integrated peak areas were used for quantification relative to external calibration curves. General linear model and Tukey's test 95% confidence grouping analyses were performed in Minitab 16 statistical software (Minitab Inc., State College, Pennsylvania, USA). Pair-wise statistical analyses were performed using Student's t-tests between sets of biological replicates from different treatments collected during the same time period only. The metabolomics data analysis server Metaboanalyst 3.0 (Wishart Research Group, the University of Alberta, Edmonton, Alberta, Canada) was used to identify clustering of metabolites profiles between samples based on PCA. Autoscaling and cube-root normalization was performed on analyte levels prior to Metaboanalyst analyses. Permutational ANOVA (PERMANOVA) was used to test for significant differences (p < 0.05) in metabolite profiles between Fe treatments and collection times. PERMANOVA tests were run using unrestricted permutation of data and 999 permutations on a resemblance matrix of autoscaled metabolite profiles based on Euclidean distances (Anderson 2005). PERMANOVA tests were performed using PRIMER software (v.6.0) with the PERMANOVA+ add-on (v.1.0.6) (PRIMER-E, Plymouth Marine Laboratory, UK). The Visualization and Analysis of NeTworks with related Experimental Data program, VANTED v2.1.0 was used to map the log base 2 transformations of Fe deficiency associated induction factors of metabolite data on to custom made biochemical pathways maps developed from pathway maps available at the KEGG pathway database (available at http://www.genome.jp/kegg/pathway.html).

Additional files
Additional file 1: Table S1A. Iron deficiency induced gene expression detected by qRT-PCR. Table S1B. Diurnal pattern of Fe deficiency induced gene expression detected by qRT-PCR. (DOCX 21 kb) Additional file 2: Table S2A. Amines concentrations from Fe-sufficient and Fe-deficient rice roots. Table S2B. Iron deficiency associated induction of root amines concentrations. (DOCX 38 kb) Additional file 3: Figure S1. Levels of 3-P-glycerate in Fe-sufficient and Fe-deficient rice roots quantified by LC-qTOF-MS. Figure S2. Levels of ascorbic acid in Fe-sufficient and Fe-deficient rice roots quantified by LC-qTOF-MS. (DOCX 117 kb) Additional file 4: Figure S3. Induction of precursors to methionine and DMA biosynthesis within the rice Fe deficiency response. (DOCX 831 kb) Additional file 5: Table S3. Iron replete hydroponic growth solution.  Tukey's HSD test: Tukey's honest significant difference test; VANTED: Visualization and analysis of networks with related experimental data program; YL: Youngest leaf; YSL: Yellow stripe-like