Genome-wide analysis of overlapping genes regulated by iron deficiency and phosphate starvation reveals new interactions in Arabidopsis roots

Iron (Fe) and phosphorus (P) are essential mineral nutrients in plants. Knowledge regarding global changes in the abundance of Fe-responsive genes under Pi deficiency as well as the processes these genes are involved in remains largely unavailable at the genome level. In the current study, we comparatively analyzed RNA sequencing data sets relative to Fe deficiency (NCBI: SRP044814) and Pi starvation (NCBI: SRA050356.1). Analysis showed a total of 579 overlapping genes that are responsible for both Fe deficiency and Pi starvation in Arabidopsis roots. A subset of 137 genes had greater than twofold changes in transcript abundant as a result of the treatments. Gene ontology (GO) analysis showed that the stress-related processes ‘response to salt stress’, ‘response to oxidative stress’, and ‘response to zinc ion’ were enriched in the 579 genes, while Fe response-related processes, including ‘cellular response to nitric oxide’, ‘cellular response to iron ion’, and ‘cellular iron ion homeostasis’, were also enriched in the subset of 137 genes. Co-expression analysis of the 579 genes using the MACCU toolbox yielded a network consisting of 292 nodes (genes). Further analysis revealed that a subset of 90 genes were up-regulated under Fe shortage, but down-regulated under Pi starvation. GO analysis in this group of genes revealed an increased cellular response to iron ion/nitric oxide/ethylene stimuli. Promoter analysis was performed in 35 of the 90 genes with a 1.5-fold or greater change in abundance, showing that 12 genes contained the PHOSPHATE STARVATION RESPONSE1-binding GNATATNC cis-element within their promoter regions. Quantitative real-time PCR showed that the decreased abundance of Fe acquisition genes under Pi deficiency exclusively relied on Fe concentration in Pi-deficient media. Comprehensive analysis of the overlapping genes derived from Fe deficiency and Pi starvation provides more information to understand the link between Pi and Fe homeostasis. Gene clustering and root-specific co-expression analysis revealed several potentially important genes which likely function as putative novel players in response to Fe and Pi deficiency or in cross-talk between Fe-deficient responses and Pi-deficient signaling.


Background
The evolutionary ability of iron (Fe) to change oxidation states between Fe(III) and Fe(II) renders it irreplaceably important in many essential processes associated with basic redox reactions, such as in photosynthesis, respiration and many vital enzymatic reactions [1][2][3][4][5]. Although Fe is abundant in the earth's crust, it is one of the least available elements for plants in aerobic soils with neutral to basic pH [1][2][3][4][5]. Approximately 30 % of the land worldwide consists of alkaline soils, leading to a demand in bioavailable Fe for plant fitness [1,5,6]. As a consequence, Fe deficiency is a major constraint in crop yield and quality [7]. In contrast, in acidic and anaerobic conditions, accumulation of excess Fe is toxic to plant growth and development due to formation of potentially harmful reactive oxygen species (ROS). Plants therefore must tightly regulate cellular Fe homeostasis to allow for effective acquisition, distribution and utilization of Fe [1,8,9].
Under Fe-deficient conditions, Arabidopsis (Arabidopsis thaliana) and other dicotyledonous and non-graminaceous monocotyledonous plants use a reduction strategy, referred to as strategy I [10], to increase Fe bioavailability. In this strategy, acidification of the rhizosphere mediated by the H + -translocating P-type ATPase AHA2 [6,11] occurs as the first step, which leads to an increase in the concentration of chelated Fe(III). Fe(III) is subsequently reduced to soluble Fe(II) by the root surface-localized ferric chelate reductase FERRIC-REDUCTION OXI-DASE2 (FRO2) [12]. Soluble Fe(II) is then transported into epidermal cells by the Fe-REGULATED TRANS-PORTER1 (IRT1) [13]. Being the major components of the Fe acquisition system, FRO2 and IRT1 are the major players controlling entry of Fe from the soil into cells. At the transcriptional level, expression of both genes is coordinately regulated by the basic helix-loop-helix (bHLH) transcription factor FER-LIKE Fe DEFICIENCY-INDUCED TRANSCRIPTION FACTOR (FIT), but not the transcription factor POPEYE, which is also involved in Fe homeostasis [14][15][16][17]. FIT forms heterodimers with bHLH38 and bHLH39 and positively regulates a subset of Fe-responsive genes, including three key genes required for Fe acquisition [12-14, 18, 19]. Recent studies have shown that the transcription factors bHLH100 and bHLH101, which belong to the Ib sub-group of bHLH proteins, are also involved in Arabidopsis Fe deficiency responses by interacting with FIT [20] or via a FIT-independent manner [21].
Studies have shown that expression of FRO2 and IRT1 is tightly controlled both locally and systemically [22,23]. However, in some cases disrupted Fe signaling in several mutants, such as frd3 [24,25], opt3 [26] and the quadruple nicotianamine synthase mutant nas4x-1 [27] in Arabidopsis, dgl and brz mutants [28][29][30] in pea (Pisum sativum) and the chloronerva mutant chln [31] in tomato (Solanum lycopersicum), constitutively activates expression of Fe acquisition genes even under sufficient Fe conditions. By contrast, FRO2 and IRT1 expression has been documented to be decreased under phosphate (Pi)-deficient conditions [32][33][34][35][36]. Currently, the predominate explanation for decreased expression of Fe acquisition genes under Pi-deficient conditions is that Pi deficiency results in enhanced Fe accessibility to plants in the media, which leads to an over accumulation of Fe in plants, subsequently causing down-regulated expression of Fe-responsive genes. However, if the Pi-deficient media without available Fe or with low concentrations of Fe, does the down-regulated expression of Fe-responsive genes occur? A recent report showed that PHOSPHATE STARVATION RESPONSE1 (PHR1), a major regulator of the Pi deficiency response, could bind the promoter of the Fe storage gene Ferritin1 through the imperfect palindromic sequence motif P1BS (PHR1 binding sequences, GNATATNC), strongly supporting the link between Fe and Pi homeostasis [37]. However, it remains an open question whether this link exists or not in phr1 mutant plants.
Moreover, although down-regulation of Fe deficientinduced Fe acquisition genes under Pi deficient conditions has been documented [33,38,39], knowledge regarding genome-wide transcriptional changes of Feresponsive genes under Pi deficiency remains unavailable, and the processes of the genes involved are largely unknown. To provide systemic information about transcriptional changes in Fe-responsive genes under Pi deficiency and to further extend knowledge of the relationship between Fe and Pi at the transcript level, we mined and re-analyzed previous RNA sequencing (RNA-seq) data sets relative to transcriptome profiling in Fe-deficient [40] and Pi-deficient Arabidopsis roots [36], with an emphasis on 579 overlapping genes that respond to both Fe and Pi deficiency. We revealed that a subset of 137 genes had a twofold or greater change in abundance under either of the treatments. A subset of 90 genes with an increased abundance under Fe deficiency, but a decrease under Pi deficiency, may be critical for Fe responses under Pi-deficient conditions. By gene clustering and root-specific co-expression analysis, we revealed several potentially important genes that likely function as putative novel players in response to Fe and Pi deficiency or in the cross talk between Fe deficient responses and phosphate-deficient signaling, which may be determined in follow-up experiments. Finally, we found that FIT-regulated genes were down-regulated by Pi deficiency, and an extent of Fe in the Pi deficient media is required for this down-regulation, suggesting that, besides FIT, PHR1, Fe itself might be a critical factor involved in the transcriptional regulation under both Pi-and Fe-deficiency.

Genes responsible for Fe and Pi deficiency in Arabidopsis roots
Previously published RNA-seq data sets [36,40] relative to Fe and Pi deficiency in Arabidopsis roots were re-analyzed, and differentially expressed genes (P < 0.05) upon Fe deficiency were compared with those (P < 0.05) exposed to Pi deficiency. Subsequent analyses focused on the 579 overlapping genes (Additional file l) as shown in Fig. 1. Of the 579 genes, 137 showed an increase or decrease in transcript abundance, with fold changes greater than twofold under either of the treatments (Table 1; Additional file 2). Fe acquisition genes FRO2 and IRT1, copper transporter COPT2, Fe(II)-dependent oxygenase gene AT3G12900, cytochrome P450 CYP82C4 (AT4G31940), mannose-binding lectin protein gene AT1G52120, glutathione transferase lambda 1 GSTL1 and amino acid transporter gene AT5G38820 showed the strongest induction under Fe deficiency and were up-regulated by more than 50-fold (Table 1; Fig. 2a). Excluding AT1G52120, these genes were among the most repressed under Pi deficiency and were down-regulated by two to more than tenfold (Table 1; Fig. 2a). Genes encoding transcriptional factor bHLH039, ZIP9, zinc binding protein (AT1G74770), ATROPGEF10, receptor like protein 24 RLP24, phloem protein 2-B6 and other functionally unknown proteins were among the second group of highly induced genes following Fe deficiency and were up-regulated by more than fivefold (Table 1; Fig. 2a, b). The most induced genes following Pi deficiency were AtOCT1, an unknown protein gene AT5G20790 and a major facilitator protein gene AT1G30560, which were induced by more than 50-fold (Table 1; Fig. 2c). Highly induced genes under Pi deficiency were ATPS3 (phosphate starvation-induced gene 3), SQD2 (sulfoquinovosyl diacylglycerol 2) and U-box domain-containing protein kinase gene AT5G65500, with changes more than fivefold (Table 1; Fig. 2d, e). Interestingly, a subset of genes involved in lignin biosynthesis was induced by both Fe and Pi deficiency (Table 1; Fig. 2f ).
Gene ontology (GO) analysis of the 579 overlapping genes revealed that stress-related processes, including 'response to salt stress' , 'response to oxidative stress' and 'response to zinc ion' , were enriched (Additional file 3), while analysis of the subset of 137 genes showed that Fe response-related processes, including 'cellular response to nitric oxide' , 'cellular response to Fe ion' and 'cellular Fe ion homeostasis' , were also enriched (Additional file 4).

Gene expression patterns of overlapping genes
Expression patterns of the 579 overlapping genes were divided into four types according to changes at the transcript level under two stress conditions (Fig. 3a). Type one was composed of 223 genes (of which 24 genes were down-regulated by more than twofold) with decreased transcription under both Pi and Fe deficiency. GO analysis of this group of genes revealed that the processes of 'embryo development ending in seed dormancy' , 'microtubule-based process' and 'chloroplast organization' , were most enriched (Fig. 3b). In contrast, transcript abundance of 169 genes in the type two category were shown to be increased following both Pi and Fe deficiency, with processes of 'glucosinolate and leucine biosynthesis' and 'UV response' being enriched (Fig. 3b). Type three was composed of 97 genes with increased transcript abundance under Pi deficiency, but decreased under Fe deficiency. In contrast, transcript abundance of the 90 genes in type four were decreased under Pi deficiency and increased under Fe deficiency. GO enrichment analysis showed that the zinc-related processes 'response to zinc ion' , 'zinc ion transport' and 'galactose metabolic process' were enriched in type three (Fig. 3b). Iron-related processes 'cellular response to Fe ion' , 'cellular response to nitric oxide' , 'cellular response to ethylene stimulus' , 'cellular Fe ion homeostasis' and 'protein import into nucleus' were enriched in type four (Fig. 3b).

Identification of overlapping gene modules by co-expression analysis
Stress-specific variability in gene expression may occur at the individual gene level, but can also occur in a coordinated manner. To determine functional modules, coexpression networks (i.e., groups of genes that show similar expression patterns under diverse conditions) of the 579 overlapping genes were generated using MACCU software [41]. Pairwise co-expressed genes were selected with a Pearson correlation coefficient cutoff of 0.7 [36,41]. The 300 publicly available microarrays that were mined for co-expression analysis discriminated between root-related experiments. As such, the co-expression relationships reported herein are restricted to roots [42][43][44]. Co-expression relationships between these genes were visualized using Cytoscape (http://www.cytoscape. org). This analysis yielded a network composed of 292 nodes (genes) and 1595 edges (correlations between genes; Additional file 5). The network can be further   The largest module was composed of 210 genes, most of which are associated with stress (Additional file 6). GO enrichment analysis revealed that the biological processes 'glucosinolate biosynthetic process' , 'response to cadmium ion' , 'response to salt stress' and 'leucine biosynthetic process' were most enriched in this module (Additional file 7). Via connection to the zinc binding protein gene AT1G74770, two marker genes strongly induced by Fe deficiency, IRT1 (AT4G19690) and CYP82C4 (AT4G31940), were associated with this module. Co-expression analysis of the subset of 137 overlapping genes with changes greater than twofold yielded a network consisting of 48 nodes and 56 edges. The Fe deficiency-regulated marker genes IRT1 and CYP82C4 (AT4G31940) remained in the network (Fig. 4a). This network can be divided into one large (26 genes) and six small clusters (Fig. 4a). Detailed expression information of these genes upon Fe or Pi deficiency is shown in Fig. 4b. GO enrichment analysis of the genes involved in the co-expression network revealed that the biological processes 'cellular response to nitric oxide' , 'cellular response to Fe ion' and 'cellular response to ethylene stimulus' were enriched (Table 2).
To search for potentially functional novel modules, coexpression analysis was applied to the subset of 90 genes that were induced by Fe deficiency, but down-regulated by Pi deficiency (Additional file l). A network containing 26 nodes and 17 edges was created using the same criteria (Fig. 5). The network can be divided into 10 small clusters (none with more than ten nodes), with the largest one containing several Fe-responsive marker genes and one transcriptional factor WRKY 17 (Fig. 5). The second largest cluster was composed of four genes, including the Pi homeostasis regulator SIZ1 (Fig. 5). For 97 genes induced by Pi deficiency but repressed by Fe deficiency, co-expression analysis resulted in a network containing 26 nodes and 29 edges that were divided into one large and two small clusters (Additional file 8).

Analysis of P1BS motif in 35 genes induced by Fe deficiency but down regulated by Pi deficiency
A subset of 35 genes in the type four category had an increase in transcript abundance under Fe deficiency but a decrease under Pi deficiency, with changes more than 1.5-fold (Additional file 9). A 3000 bp sequence upstream of the translation start (named −3000 bp) of these 35 genes was retrieved from TAIR10 and used to search the PHR1 recognition sequence 5′-GNATATNC-3′ (P1BS motif ). Results showed that 12 of the 35 genes contained at least one P1BS motif, and seven of the 12 genes contained the P1BS motif within −1000 bp of their promoter regions (Table 3). In total, 11 P1BS patterns were hit in the 12 genes, with 5′-GTATATGC-3′ and 5′-GTA-TATTC-3′ being the most frequent (5 and 3 hits out of 18 total hits, respectively).

Down-regulation of Fe-acquisition genes upon Pi deficiency is dependent on Fe concentration in the media
To determine how Fe acquisition genes are down-regulated by Pi deficiency and whether this down-regulation is dependent on PHR1, we investigated changes in genes that were most induced by Fe deficiency (including the two Fe acquisition genes FRO2 and IRT1) at the transcript level in wild type and the phr1 mutant under varied growth conditions as follows: Pi deficiency (−Pi, in which   Fig. 6a and in agreement with previously reported results [45], transcriptional expression of PHR1 was not significantly regulated by Pi deficiency in wild type plants and could not be detected in phr1 mutant plants under both Pi sufficient and deficient conditions. As a control, the expression of SPX1 [46,47], a Pi-responsive marker gene, was significantly induced by Pi deficiency. Consistent with our transcriptomic data, transcriptional expression of the Fe acquisition genes IRT1 and FRO2 as well as the Fe deficiency-induced marker gene CYP82C4 was significantly  down-regulated under Pi deficiency in both Col-0 and phr1 roots (Fig. 6a). Because all these Fe-responsive genes tested were mainly regulated by transcription factor FIT [14], we thus tested whether the expression of FIT itself was affected or not by Pi deficiency. As shown in Fig. 6a, the expression level of FIT was significantly lower in Pideficient roots than in Pi-sufficient roots in wild type plants. In addition to FIT, another transcription factor PYE [34], regulating the expression of another subset of Fe-responsive genes, has been reported to be required for plant Fe homeostasis. However, both PYE and its target AT1G74790 were not affected by Pi deficiency (Fig. 6a).
To determine whether this down-regulation is dependent on Fe concentrations in the media, we compared  (Fig. 6b). Similar to results in wild type plants, these genes were significantly induced under Fe deficiency and Fe and Pi deficiency (−Fe−Pi) in the phr1 mutant roots (Fig. 6b). However, down-regulation of gene expression under Pi deficiency was not all significantly attenuated by 5 µM Fe in Pi-deficient media in the phr1 mutant roots (Fig. 6b).

Discussion
As an essential element for all living organisms, particularly as a major constraint in crop yield and quality, Fe deficiency responses in plants have been extensively studied in the last decade [1,6,9]. With the emergence of high throughput research platforms, many genes and proteins have been revealed to be regulated by Fe deficiency [40,[48][49][50][51][52][53][54][55]. Evidence has shown that transcriptional expression of some Fe-responsive genes can be altered due to deficiencies or excesses of mineral elements, including cross-talk between Fe and other mineral elements. For example, the Fe transporter LeIRT1 is reported to be up-regulated by potassium (K) deficiency, as revealed by microarray analysis [56], and expression of the K transporter gene LeKC1 was induced not only by K starvation but also by Fe deficiency [56]. Via comprehensive analysis of Fe-responsive protein kinase (PK) and protein phosphatase (PP) genes, we found that strong over-representation of PK and PP genes that encode proteins is involved in K homeostasis, which supports the link between potassium uptake and Fe deficiency [44]. The ameliorative effect of K supply on Fe-deficient responses was previously reported [57]. Although several lines of evidence have suggested a link between Pi and Fe homeostasis [32-35, 37-39, 58], little genome-wide information on transcriptional expression changes in Fe-responsive genes under Pi deficiency is available, and the biological processes that these genes are involved in remain elusive in Arabidopsis. By mining previous RNA-seq data sets, we present comprehensive information on transcriptional expression of overlapping genes regulated by Fe and Pi deficiency in Arabidopsis roots. In total, 579 overlapping genes, or less than 20 % of all differentially expressed genes evaluated in each treatment, were responsive to both Fe and Pi deficiency. Only 137 of the 579 genes had greater than twofold changes in transcript abundance (Additional file l; Table 1). Many of the most strongly induced genes under Fe deficiency, such as AT3G12900, IRT1, FRO2, CYP82C4 and AT5G38820 [18,41,54], are among the 579 overlapping genes, while most of the strongest induced Pi deficiency-induced marker genes, such as pyridoxal phosphate phosphataserelated protein gene AT1G17710, transposable element gene AT2G04460, ATISP1, SPX3, APT1 and AT4, are not overlapping [33,35]. GO enrichment analysis of the 137 genes with changes greater than twofold (Additional file 4) showed that Fe response-related processes such as 'cellular response to Fe ion' and 'cellular Fe ion homeostasis' were enriched, but none of the Pi response-related processes were pronounced, suggesting that plant responses to Fe deficiency might be more specific than responses to Pi deficiency under the conditions presented herein. The most strongly induced Fe-responsive genes were clustered together and downregulated under Pi deficiency, except for AT1G52120 in which transcript abundance was increased under both stress conditions (Fig. 2a). In this cluster, IRT1, FRO2, and BHLH039 are known to be involved in Fe acquisition and transcriptional regulation, and COPT2 was confirmed to participate in cross talk between Fe deficiency responses and low phosphate signaling in a recent study [59]. Other genes in the group (Fig. 2a), such as AT3G12900, CYP82C4, AT5G38820 and AT1G52120, do not have defined functions currently, but may be involved in responses to Fe deficiency or Pi deficiency or cross-talk between Fe deficiency responses and phosphate-deficient signaling. Another group of interesting genes are AT5G20790, AtOCT1 and AT1G30560 (Fig. 2c) given that their transcriptional expression was among the most highly induced under Pi deficiency. In particular, both AT1G30560 and AtOCT1 were significantly up-regulated upon Fe deficiency, suggesting that these two genes might play important roles in responses to both stresses. In animals, organic cation/carnitine transporters (OCTs) are associated with homeostasis and distribution of various small endogenous amines (e.g. carnitine, choline) and detoxification of xenobiotics like nicotine. AtOCT1 has been reported to be involved in Arabidopsis root development. Knockout of AtOCT1 expression results in a higher degree of root branching compared to the wild type in vitro. This disordered development may be due to an inability to transport carnitine [60]. It has been well established that the number and length of lateral roots are increased under Pi deficiency in Arabidopsis and other plants. Therefore, whether AtOCT1-mediated transport of carnitine or related chemicals is involved in lateral root development under Pi deficiency remains elusive.
GO enrichment analysis of the 579 overlapping genes revealed that these Fe-and Pi-responsive genes were associated with diverse biological processes (Additional file 3), particularly with the GO categories 'response to salt stress' , 'response to oxidative stress' and 'response to zinc ion' (Additional file 3). These results imply that acclimation of plants to Fe and Pi deficiency and possibly other nutritional stresses is associated with profound changes in the transcriptome, including stress-specific responses such as alteration of ribosome composition [43] and other general responses. Only four (AT1G27760, AT3G04720, AT4G11650 and AT5G24090) of the 21 genes associated with 'response to salt stress' had an increase in transcript abundance greater than 1.5-fold, suggesting that this common response to Fe and Pi deficiency might be less important than Fe response-related processes. GO enrichment analysis of the most responsive genes (i.e., those with greater than 1.5-fold change) revealed that Fe response-related processes, were enriched, but none of the Pi response-related processes were (Additional file 4),suggesting that Pi deficiency has more pronounced effects on Fe homeostasis than Fe deficiency has on Pi homeostasis.
Functional annotation of a given gene is the most important goal in modern molecular biology and is essential for understanding how the cell works. All omics studies are discovery tools and are not capable of defining gene function. The actual functions of differentially expressed genes under certain conditions discovered by high throughput platforms require further experimental evidence. However, current research platforms can discover hundreds to thousands of differentially expressed genes in a single run, and most of them are annotated as function unknown. Functional exploration of every differentially expressed gene without selection would be extremely laborious and impossible. Fortunately, coexpression analysis provides the option to choose genes of interest for further study. The basic idea of co-expression analysis is that genes that show transcriptionally coordinated expression patterns under diverse conditions are often functionally related [61], thus allowing functional predictions regarding genes with unknown functions inferred from their co-expression relationships with genes with known functions [62,63]. Using coexpression analysis, we discovered ten, six and ten potentially critical regulatory modules with diverse nodes from inputs of the 579 (total overlapping genes), 137 (genes with changes greater than twofold) and 90 Fe deficiencyinduced, Pi deficiency-repressed genes (Additional file 5;Figs. 4a,5). Unexpectedly, only 50, 35 and 29 %, respectively, of the input genes were associated with formation of co-expression networks, suggesting that the majority of overlapping genes are functionally diverse and involved in a variety of biological processes. The network obtained from the group of 90 Fe deficiency-induced, Pi deficiency-repressed genes (Fig. 3a) is of particular interest. In this network, several genes may play important roles in responses to Fe and Pi deficiency. For instance, the gene AT1G74770 annotated with zinc ion binding protein showed a strong relationship with the Fe transporter IRT1, implying that this gene may be required for a Fe response. Another putative zinc ion binding protein encoding gene, AT1G18910, was shown to be connected to AT1G74770 and the transcription factor gene WRKY17, suggesting that these genes may also be involved in plant adaptation to Fe deficiency or zinc toxicity elicited by excess zinc under Fe deficiency.
It is generally accepted that a group of genes with similar expression patterns might be positively and/or negatively regulated by the same regulator(s). In Arabidopsis, the PHR1 transcription factor (TF) and its homolog PHL1 consist of the central regulatory system controlling transcriptional expression of a subset of Pi deficiency response genes by binding to the P1BS motif in promoter regions; while FIT and PYE are two major TFs regulating transcriptional expression of two subsets of Fe deficiency response genes. 35 out of the 579 genes, including Fe acquisition genes IRT1 and FRO2 and Fe responsive marker gene CYP82C4, were induced under Fe deficiency but down-regulated under Pi deficiency with changes greater than 1.5-fold. qPCR examination (Fig. 6a, b) confirmed that both IRT1 and FRO2 as well as CYP82C4, mainly regulated by FIT in response to Fe deficiency, were down-regulated by Pi deficiency in wild type plants, probably due to the decreased abundance of FIT (Fig. 6a). By contrast, the transcriptional expression of both PYE and its target AT1G74790 was not altered in response to Pi deficiency in both wild type and phr1 mutant plants. Taken together, these results suggest that FIT-regulated but not PYE-regulated Fe-response genes are affected by Pi deficiency and the down-regulation by Pi deficiency might be partially due to the down-regulation of FIT (since although the transcript abundance of FIT was not significantly different between Pi-sufficient and -deficient conditions in phr1 mutant plants, the transcriptional expression of IRT1 and FRO2 as well as CYP82C4 was still significantly down-regulated by Pi deficiency). In addition, under Pi sufficiency, although the transcript level of FIT was not different between wild type and phr1 mutant plants, the transcript abundance of IRT1 and FRO2 as well as CYP82C4 was still significantly down-regulated in the mutant plants (Fig. 6a, b). These results indicate that, besides FIT and P1BS motif (CYP82C4 doesn't contain a P1BS motif in the −3000 bp sequence of its promoter region), some other factors may be involved in the down-regulation of gene expression under Pi deficiency.
Indeed, only 34 % of the genes (12 of 35) contained at least one P1BS motif in promoter regions (−3000 bp sequence upstream of the translation start) and only 20 % (7 of 35) had a P1BS motif within −1000 bp of their promoter regions, further suggesting that other positive or negative regulators might be involved in down-regulation of these Pi-responsive genes. One of these regulators may be Fe itself. It has been reported that Pi deficiency results in enhanced Fe accessibility to plants in the media, which leads to an over accumulation of Fe in plants, subsequently causing down-regulated expression of Fe-responsive genes. This point of view was confirmed by supply of different Fe concentrations in the Pi deficient media (Fig. 6b). If no additional Fe was supplied to the Pi deficient media (−Fe−Pi), the transcriptional expression of all tested genes was induced both in wild type and phr1 mutant plants, an expression pattern similar to the one of Fe deficiency (Fig. 6b). This result suggests that an extent of Fe in the Pi deficient media is required for the downregulation of Fe-responsive genes under Pi deficiency. Indeed, Pi-deficiency caused down-regulation was much enhanced by supply of 5 µM Fe in the Pi deficient media both in wild type and phr1 mutant plants (Fig. 6b). In the future, the dose effects of Fe in the Pi-deficient media on the transcriptional expression of Fe-responsive genes need further validation.

Conclusions
In summary, we provide genome-wide information on the transcriptional expression of 579 overlapping genes that responded to both Fe and Pi deficiency in Arabidopsis roots and the biological processes that the genes are involved in. Gene clustering and root-specific coexpression analysis revealed several potentially important genes, including CYP82C4 and AT5G38820, which likely function as putative novel players in response to Fe and Pi deficiency or in cross-talk between Fe-deficient responses and Pi-deficient signaling. These results imply that Pi deficiency has more pronounced effects on Fe homeostasis than Fe deficiency has on Pi homeostasis.

Quantitative reverse transcription-PCR
Total RNA was isolated using the RNeasy Plant Mini Kit (Qiagen) and treated with DNase using the TURBO DNA-free Kit (Ambion) as suggested by the manufacturer. cDNA was synthesized and qPCR was performed according to a previous report [40] using the SYBR Green PCR Master Mix (Applied Biosystems) with programs recommended by the manufacturer in the ABI Prism 7500 Sequence Detection System (Applied Biosystems). The melting temperature of the primers ranged from 58 to 62 °C. Primer pairs were selected using Primer3 (http://primer3.sourceforge.net/). Elongation factor1-β2 (At5g19510) and Tubulin3 (At5g19770) were used as internal controls (transcript abundance of both genes did not change under Fe and Pi deficiency) for transcript normalization. The primers used in this study are listed in Additional file 10. Three independent replicates were performed for each sample. The delta threshold cycle (∆ct = the ct of a gene−the ct of the reference) was used to determine the relative amount of gene expression. Student's t test (P < 0.05) was used to compare differences between samples grown under treatment and control conditions.

Data collection and processing
Transcriptomic data sets of roots from 13-day-old Arabidopsis seedlings grown in the presence or absence of Fe or Pi by RNA-seq were downloaded from a public database (NCBI: SRP044814, SRA050356.1). The 3106 and 3296 differentially expressed genes (P < 0.05) upon Pi and Fe deficiency were compared, and the resulting 579 overlapping genes were subsequently analyzed as shown in Fig. 1. Microarray data of 2671 ATH1 arrays from the NASCarray database (http://affymetrix.arabidopsis.info/) were downloaded and normalized using the RMA function in the Bioconductor Affy package software. Three hundred root-related arrays were manually identified as