Comparative expression profiling reveals a role of the root apoplast in local phosphate response
- Wolfgang Hoehenwarter†1,
- Susann Mönchgesang†2,
- Steffen Neumann2,
- Petra Majovsky1,
- Steffen Abel3, 4, 5 and
- Jens Müller†3Email author
© Hoehenwarter et al. 2016
Received: 19 January 2016
Accepted: 18 April 2016
Published: 28 April 2016
Plant adaptation to limited phosphate availability comprises a wide range of responses to conserve and remobilize internal phosphate sources and to enhance phosphate acquisition. Vigorous restructuring of root system architecture provides a developmental strategy for topsoil exploration and phosphate scavenging. Changes in external phosphate availability are locally sensed at root tips and adjust root growth by modulating cell expansion and cell division. The functionally interacting Arabidopsis genes, LOW PHOSPHATE RESPONSE 1 and 2 (LPR1/LPR2) and PHOSPHATE DEFICIENCY RESPONSE 2 (PDR2), are key components of root phosphate sensing. We recently demonstrated that the LOW PHOSPHATE RESPONSE 1 - PHOSPHATE DEFICIENCY RESPONSE 2 (LPR1-PDR2) module mediates apoplastic deposition of ferric iron (Fe3+) in the growing root tip during phosphate limitation. Iron deposition coincides with sites of reactive oxygen species generation and triggers cell wall thickening and callose accumulation, which interfere with cell-to-cell communication and inhibit root growth.
We took advantage of the opposite phosphate-conditional root phenotype of the phosphate deficiency response 2 mutant (hypersensitive) and low phosphate response 1 and 2 double mutant (insensitive) to investigate the phosphate dependent regulation of gene and protein expression in roots using genome-wide transcriptome and proteome analysis. We observed an overrepresentation of genes and proteins that are involved in the regulation of iron homeostasis, cell wall remodeling and reactive oxygen species formation, and we highlight a number of candidate genes with a potential function in root adaptation to limited phosphate availability. Our experiments reveal that FERRIC REDUCTASE DEFECTIVE 3 mediated, apoplastic iron redistribution, but not intracellular iron uptake and iron storage, triggers phosphate-dependent root growth modulation. We further highlight expressional changes of several cell wall-modifying enzymes and provide evidence for adjustment of the pectin network at sites of iron accumulation in the root.
Our study reveals new aspects of the elaborate interplay between phosphate starvation responses and changes in iron homeostasis. The results emphasize the importance of apoplastic iron redistribution to mediate phosphate-dependent root growth adjustment and suggest an important role for citrate in phosphate-dependent apoplastic iron transport. We further demonstrate that root growth modulation correlates with an altered expression of cell wall modifying enzymes and changes in the pectin network of the phosphate-deprived root tip, supporting the hypothesis that pectins are involved in iron binding and/or phosphate mobilization.
KeywordsArabidopsis thaliana Phosphate deficiency Root growth Proteomics Transcriptomics Iron transport Cell wall Pectin
Inorganic phosphate (Pi) is an essential macronutrient for plant growth and development. Despite its high abundance in the rhizosphere, bioavailability of Pi is typically limited because its majority is bound in organic compounds or complexed with metal ions such as Ca (alkaline soils), Fe or Al (acidic soils) . Thus, plants evolved strategies to enhance Pi acquisition and to conserve or remobilize Pi from internal sources to adapt to Pi limiting conditions. Previous efforts elucidated some of these adaptive responses, including the identification of high-affinity Pi transport systems, the characterization of diverse metabolic bypass reactions, the reutilization of Pi from phospholipids, and many more . Most of the Pi starvation response (PSR) genes involved in these systemic adjustments are regulated by the myb transcription factor PHR1 (PHOSPHATE STARVATION RESPONSE1) [3–6].
Dynamic redesign of the root system architecture (RSA) provides another strategy to maintain cellular Pi supply. In Arabidopsis, low external Pi availability is locally sensed by the growing root tip, which causes reduction of cell elongation and meristematic activity at the site of Pi depletion. The resultant inhibition of root growth is accompanied by accelerated formation of root hairs and development of lateral roots to increase the absorptive surface for topsoil exploration [7, 8]. The development of a densely branched and/or shallow root systems increases Pi starvation tolerance in several plant species, including agronomically important crops such as barley, lupin, soybean or common bean . Several Arabidopsis mutants with altered Pi dependent root growth responses have been described [10–18]. However, for most of the underlying genes only little information is available how they affect Pi sensing and root growth modulation. LPR1 (LOW PHOSPHATE ROOT1), its closely related paralog LPR2, and PDR2 (PHOSPHATE DEFICIENCY RESPONSE2) have been identified as central players in local root Pi sensing [11, 13, 19]. PDR2, which codes for the single P5-type ATPase of unknown substrate-specificity (AtP5A), and LPR1, encoding a multicopper oxidase, are expressed in overlapping domains of the root apical meristem (RAM). LPR1 and PDR2 interact genetically and are required for meristem maintenance and cell elongation in Pi-deprived roots. Importantly, the lpr1lpr2 mutation impedes local root growth inhibition under Pi limitation and suppresses the hypersensitive short-root phenotype of pdr2 plants, indicating that they act in the same pathway [11, 13].
Previous work revealed that external Fe availability modifies local Pi sensing [11, 13, 20]. A number of studies observed that Pi-starved Arabidopsis and rice plants accumulate elevated levels of Fe in the root and the shoot [20–23], which has been suggested as a proactive strategy to mobilize Pi from insoluble Fe complexes . Fe participates in the formation of reactive oxygen species (ROS) and it has been proposed that Fe toxicity causes local root growth inhibition . We recently provided evidence for apoplastic LPR1 ferroxidase activity and uncovered a major role of the LPR1-PDR2 module for root tip-specific deposition of Fe3+ in cell walls (CW) of the RAM and elongation zone (EZ) during Pi limitation . We further showed that Fe accumulation in the RAM is massively enhanced in Pi-starved pdr2 roots, but suppressed in the insensitive lpr1lpr2 line. Fe deposition coincides with sites of ROS generation and triggers CW thickening and callose accumulation, which interferes with cell-to-cell communication, RAM maintenance, and cell elongation.
In recent years, a set of transcriptome profiling studies provided significant insights into the transcriptional changes upon Pi deficiency in Arabidopsis [6, 21, 24–28]. In addition, a complementary transcriptome and proteome study highlighted the convergence of mRNA and protein expression profiles on lipid remodeling and glucose metabolism upon Pi-deprivation . In this study, we performed comparative transcriptome and proteome expression profiling on roots of Pi-replete and Pi-starved wild-type (Col-0), pdr2, and lpr1lpr2 plants in combination with a set of physiological and cell biological experiments. Our analysis emphasizes the importance of root Fe uptake and redistribution under Pi limitation. We highlight the potential role of so far unknown players in the regulation of Pi-dependent Fe-redistribution and demonstrate that apoplastic but not intracellular Fe accumulation triggers Pi-dependent root growth modulation. Consistently, we observed regulation of several CW modifying enzymes, which correlates with an increased deposition of pectin at sites of Fe accumulation. The potential role of pectin in Pi-dependent root Fe storage and Pi mobilization is discussed.
Differential gene expression correlates with genotype-specific Pi sensitivity
Identification of genotype-independent Pi-responsive genes
Pi depletion alters expression of cell wall-related genes
GO term analysis also revealed overrepresentation of genes encoding tetrapyrrole- and heme-binding proteins (GO:0046906 and GO:0020037) with oxidoreductase activity (GO:0016491) (Additional file 2: Figure S1D). This group codes for 29 peroxidases and most of those (28) belong to the 73 member-family of class III peroxidases (CIII Prx) (Additional file 7: Table S6), which are extracellular enzymes with partly antagonistic functions in ROS formation and CW dynamics . While Pi-responsive expression of 8 CIII Prx-encoding genes was similar between wild-type and pdr2 roots, 7 genes were regulated independently in each line under low Pi, and only three CIII Prx genes responded significantly to Pi limitation in lpr1lpr2 plants (Additional file 7: Table S6). Again, 19 CIII Prx genes were deregulated in pdr2 under + Pi. Thus, peroxidases may be an important link between ROS formation and CW remodeling upon Pi starvation.
Proteomics supports regulation of Pi-responsive genes in pdr2 and lpr1lpr2 mutants
Genotype-specific changes in the root proteome upon Pi deficiency were monitored in an unlabeled approach using a fast scanning high resolution accurate mass (HRAM) LC-MS system. Three biological and three technical replicates were measured for each genotype under + Pi and –Pi conditions (54 samples) yielding 3,328,368 MS/MS spectra (individual peptide measurements). 726,944 spectra could be annotated to a peptide sequence (peptide spectral match, PSM) with a global false discovery rate (FDR) threshold of 0.01 %. These PSMs were used to identify 5110 protein groups (unique proteins), each with at least one unique peptide and a global FDR threshold of 1 % (Additional file 8: Table S7). Protein abundance was inferred based on peptide abundance determined by peptide ion signal peak integration using the PROGENESIS software. Pairwise comparison of all genotypes under both growth regimes revealed 2439 differentially regulated proteins (p ≤ 0.05). Based on this list, we identified 1304 proteins that were either Pi-responsive in at least one genotype or which were already deregulated in one of the mutant lines grown on Pi-replete conditions (0.769 ≥ FC ≥ 1.3) (Additional file 9: Table S8).
Correlation of proteome and transcriptome analysis
Next, we performed GO term analysis to identify groups of proteins involved in genotype-specific Pi responsiveness. Most proteins have assigned metabolic functions in wild-type and lpr1lpr2, probably reflecting processes related to Pi recycling and mobilization. Strikingly, in + Pi condition and upon transfer to –Pi, the pdr2 line showed a significant regulation of proteins assigned as response to metal ion (GO:0010038) and oxidoreductase activity (GO:0016491). A closer examination revealed repression of 15 peroxidases in pdr2 in + Pi and induction of 9 peroxidases in –Pi condition. Within the group of repressed proteins we identified 14 CIII Prxs of which 3 enzymes were regulated at the transcript level. Only one and six Pi-responsive CIII Prx were detected in wild-type and lpr1lpr2 root extracts, respectively (Additional file 10: Table S9).
We identified the 4 genes, including FER1 and PGIP1, that were co-regulated on mRNA and protein level across all genotypes in response to Pi depletion (Additional file 11: Table S10). In wild-type, we noticed induction of PPa4 (PYROPHOSPHORYLASE 4), a candidate for Pi mobilization, and PCK1 (PHOSPHOENOLPYRUVATE CARBOXYKINASE 1), which is involved in metabolic adjustment to Pi deprivation . We further identified two hemicellulose modifying enzymes, XTH8 and XTH31 (XYLOGLUCAN ENDOTRANSGLUCOSYLASE/HYDROLASE), which were slightly decreased in low Pi. Interestingly, both enzymes were previously shown to be regulated by SIZ1 , a SUMO E3-ligase involved in Pi dependent root growth remodeling [47, 48].
GO term analysis of the 211 mRNA/protein pairs altered in pdr2 revealed an overrepresentation of metabolic processes. The second most significant term (response to metal ion) is consistent with altered metal homeostasis in pdr2 plants . For example, we noticed induction of FER3 and proteins potentially involved in detoxification of metal ion-induced ROS formation, including several GLUTATHIONE-S-TRANSFERASEs (GSTs) (Additional file 11: Table S10). We also identified F6’H1 (feruloyl-CoA 6’-hydroxylase 1), which is involved in coumarin biosynthesis and Fe-mobilization in alkaline soils [49–51]. Our datasets revealed anti-correlation of F6’H1 expression in pdr2, showing elevated protein but decreased transcript levels in –Pi and an inverse relation in + Pi (Additional file 1: Table S1, Additional file 9: Table S8), which indicates stringent regulation of F6’H1 expression in pdr2. In addition, protein level of CCoAOMT1 (caffeoyl coenzyme A O-methyltransferase 1), which converts caffeoyl-CoA to feruloyl-CoA, the substrate of F6’H1 , was also elevated in pdr2 (Additional file 9: Table S8). Thus, coumarin-mediated mobilization of Fe may be involved in Pi dependent Fe accumulation.
Integrative spCCA analysis supports Pi-dependent metal redistribution
We integrated the two –omics approaches to uncover relationships that are supported by both individual datasets. We performed a supervised penalized canonical correlation analysis (spCCA), which searches for correlations between a set of transcripts and proteins . The experimental design was integrated into the analysis to allow for biological interpretation of the derived canonical variables. The experimental factors (i.e., genotype, Pi condition, replicate sample) were provided as a binary matrix of design vectors that uniquely characterize each sample (Additional file 12: Figure S2). The supervised correlation approach seeks a linear combination of a feature subset from each -omics dataset that correlates maximally with a subset of experimental design factors. To maximize stringency, only varying transcripts and proteins were considered for spCCA. For transcriptomics, we choose a list of 1143 ANOVA filtered genes (p ≤ 0.05, var ≥ 0.12) and for proteomics a list of 47 proteins (p ≤ 0.05, var ≥ 0.4). Our analysis revealed distinct canonical variables (CVs), each representing a specific pattern correlating with a subset of proteins and/or transcripts. The first two CVs revealed structured patterns (Fig. 3f, g), while a third CV was disordered and therefore not further examined (Additional file 12: Figure S2B). The first CV mainly represented genes/transcripts (g/t) that were differentially expressed in pdr2 compared to wild-type and lpr1lpr2 independently of Pi status (Fig. 3f). We examined the top 100 g/t in this variable and found several Fe-related candidates (Additional file 13: Table S11), such as Fe chelate reductase 3 (FRO3)  and MYB10, which is required for growth in Fe deficiency . MYB10 and MYB72 mediate Fe-dependent induction of NICOTIANAMINE SYNTHASE 4 (NAS4) , which is also present in this group. NAS proteins synthesize nicotianamine, a Fe-chelator essential for Fe-remobilization in the root . We further identified a member of the ALUMINUM ACTIVATED MALATE TRANSPORTER (ALMT) family. It is of note that ALMT1 is most highly induced in all three genotypes during Pi depletion (Table 1).
The second CV mainly represented g/t that were similarly expressed in Pi-replete pdr2 and lpr1lpr2 roots but slightly differed from the wild-type. In contrast to the first CV, the majority of these g/t were Pi responsive in all genotypes. As expected, we found several known Pi acquisition g/t, including SPX1, CAX3, the phosphate transporter PT2, and the Pi starvation-inducible inorganic pyrophosphatase 1 (Additional file 13: Table S11). Many other g/t are implicated in metal homeostasis, e.g., the Fe/Zn transporters IRT1 and IRT3, the Ni transporter IREG2, the Zn/Cd transporter HMA2 or the NA transporter YSL2, further supporting our observation that metal homeostasis is strictly controlled in all genotypes upon Pi starvation.
Root growth inhibition in low Pi is independent of general Fe uptake and cellular storage
We previously reported that LPR1-dependent Fe accumulation and distribution in root tips controls RAM activity in response to low Pi . Our comparative transcriptomics and proteomics analysis of entire roots revealed Pi-responsive expression of Fe-related genes, notably FER1 and IRT1 (Table 1, Table 4), which correlated with Fe overload in Pi-starved roots of the three genotypes under study  (Additional file 14: Figure S3). To further investigate the role of Fe during the local response of roots to Pi availability, we analyzed the impact of FER1 and IRT1 loss-of-function mutants on Fe-distribution and root growth inhibition upon Pi deprivation.
Similarly, we performed Perls/DAB Fe-staining to examine Fe distribution in wild type and irt1 roots. Compared with Pi-replete wild-type seedlings, the irt1 mutant showed more intense Fe staining on the root surface of the mature root zone (Additional file 15: Figure S4B), which is in accordance with impaired Fe uptake from the rhizosphere. However, both lines displayed similar Fe staining in the RAM and EZ, which is consistent with predominant IRT1 expression in the differentiation zone  and confirms our previous study . Under Pi depletion, Fe staining increased strongly and comparably in all segments of wild-type and irt1 roots, indicating that Fe accumulation and distribution in root tips is independent of IRT1. We generated homozygous pdr2irt1 double and lpr1lpr2irt1 triple mutants and monitored primary root growth on + Pi and –Pi agar. As expected, the irt1 mutation did not affect the Pi-dependent root growth response of pdr2 and lpr1lpr2 plants (Fig. 4b), indicating IRT1-independent Fe accumulation in the root tip in response to low Pi.
Apoplastic Fe redistribution modifies Pi-dependent root growth adaptation
Long distance apoplastic Fe transport and distribution in symplastically disconnected tissues are mediated by the citrate exporter FERRIC REDICTASE DEFECTIVE 3 (FRD3) [57, 58]. Intriguingly, a previous study reported that frd3 plants display a hypersensitive short-root phenotype when grown on –Pi medium . To examine a potential role of FRD3 for mediating Pi-dependent Fe distribution via Fe-citrate complexes, we performed Perls Fe-staining (without DAB intensification to avoid oversaturation) on wild-type and frd3 roots. As previously reported [58–60], Pi-replete frd3 roots overaccumulated Fe in the vascular tissue (Fig. 4c). Within 20 h after transfer to –Pi, wild-type plants accumulated Fe in the outer cell layers, whereas frd3 roots showed enhanced Fe staining in the vasculature, particularly in differentiated root segments. Importantly, only minor differences were noted in the root tip, where Fe accumulation was slightly increased in frd3 (Fig. 4c); However, extended growth on –Pi (up to 6 days) progressively increased this difference, finally causing massive overaccumulation of Fe within the EZ and early differentiation zone of frd3 roots (Fig. 4d).
We previously showed that Pi-dependent Fe accumulation correlates with callose formation at the sites of Fe deposition (<2 days) . After transfer to –Pi (2 days), callose deposition at sites of Fe accumulation and resultant root growth inhibition were similar for wild-type and frd3 plants (Additional file 16: Figure S5A). However, extended exposure (6 days) caused callose overproduction in frd3 roots which correlated with an enhanced growth inhibition (Fig. 4d, e).
Based on our observations, we assumed that mobilization of apoplastic Fe-citrate complexes might be involved in the Pi dependent modulation of root growth. To test this, we transferred wild-type plants from + Pi conditions to + Pi or –Pi medium, supplemented with citrate, which was previously shown to restore Fe mobilization on frd3 mutants  and monitored their growth behavior. Indeed, addition of 100–250 μM citrate promotes root growth within the first two days after transfer to –Pi. However, this effect was transient and external supply of citrate eventually suppressed root growth on low Pi (Additional file 16: Figure S5B, C).
Pi deprivation modifies pectins at Fe accumulation sites
Our comparative expression profiling pointed to a role for pectin-modifying enzymes. Therefore, we studied Pi-dependent changes in the pectineous CW by using Ruthenium Red (RR), an inorganic dye that stains unesterified pectins [61, 62].
We also visualized the distribution of methyl-esterified pectin by using the hydroxylamine ferric chloride (HFC) reagent, which specifically reacts with methyl esters of pectin and results in a yellow to red coloration [62–64]. Only weak staining was evident in roots on + Pi (Fig. 5b) and transfer to –Pi did not significantly change the staining pattern in the differentiating EZ. However, higher magnification images revealed increased staining in the RAM of wild-type, with the highest intensity in the quiescent center (QC) and the cortical cell layer at the transition zone, which demarcates the border between the RAM and EZ. In contrast, pdr2 seedlings showed enhanced staining in the RAM, particularly within the QC region, but no distinct labeling of the cortical cell layer. No differences in pectin staining were detected in lpr1lpr2 roots after transfer to –Pi medium. High magnification images of RR- and HFC-stained roots revealed simultaneous accumulation of acidic and methyl-esterified pectin in the meristem of the two sensitive lines. In particular, after transfer to –Pi, strong HFC and RR staining was evident in the cortex cell layer of wild-type and in the QC region of pdr2 roots (Fig. 5c), which co-localized with major sites of Fe deposition.
Plant adaptation to Pi limitation depends on coordinated transcriptional and translational regulation of gene expression [6, 21, 24–28]. While comparative transcriptome analysis proved to be a viable approach to distinguish between local and systemic regulation in Pi-starved plants [5, 6], only little information is available on the regulation of genes and proteins associated with the Pi-dependent adaption of root system architecture. Previous work revealed that PDR2 and LPR genes act together in the local response to Pi availability by regulating cell type-specific deposition of Fe and callose in the root tip [11, 13, 19]. Here, we took advantage of the contrasting Pi-dependent root phenotype of pdr2 and lpr1lpr2 plants to investigate the associated changes in steady-state transcript and protein levels in a comparative approach. Genotype independent regulation of several PSR genes demonstrated the validity of our experiments and revealed that pdr2 and lpr1lpr2 mutants are likely not affected in the systemic response to Pi limitation (Table 1). Further analysis of our dataset revealed a number of candidate genes that are possibly involved in the Pi-dependent regulation of Fe storage and Fe redistribution as well as in the modulation of CW dynamics and/or ROS formation within the root.
Pi depletion modulates root Fe distribution
Our study revealed genotype-independent repression of numerous Fe-responsive and IRT1-coregulated genes upon transfer to Pi limitation, which likely reflects feedback regulation as a consequence of elevated Fe accumulation in Pi-starved differentiated roots. On the other hand, de-repression of Fe-related genes in Pi-replete pdr2 plants may sensitize Fe overaccumulation in limiting Pi  (Table 2, Additional file 12: Figure S2).
FER1 and related ferritins are plastid-localized Fe storage proteins protecting cells from Fe-mediated oxidative stress . Using Perls/DAB Fe staining, Reyt et al.  recently reported dot-like structures in root cells of wild-type plants that likely display ferritin-bound Fe because they are absent in fer1-3-4 roots . A previous study showed that FER1 is induced by PHR1 in low Pi independent of external Fe , indicating that FER1 may play a role in Pi-dependent Fe distribution.
Our comparative analysis revealed induction of FER1 expression on mRNA and protein level in all three lines (Table 4). Detection of Fe accumulation in dot-like structures supports the notion of intracellular Fe storage under Pi limitation, possibly as ferritin Fe (Figure S4). Importantly, Pi-dependent root growth was not affected by loss of ferritins (fer1--3-4 mutants) or loss of IRT1 in pdr2 (pdr2irt1) and lpr (lpr1lpr2irt1) mutants (Fig. 4), indicating that Pi dependent root growth modulation is independent of intracellular Fe accumulation. Our data are consistent with a recent study reporting indistinguishable primary root growth of fer1-3-4 and wild-type plants on high Fe .
Fe mobilization from the rhizosphere is facilitated by chelators such as carboxylates (e.g., citrate and malate) and coumarins, and apoplastic long distance Fe trafficking is mediated by Fe-citrate complexes [49–51, 67]. FRD3 exports citrate and the frd3 mutant is defective in apoplastic Fe translocation, causing Fe hyperaccumulation in root stele tissues [58–60]. Importantly, frd3 mutants show a hypersensitive short root phenotype in low Pi  and we demonstrated Fe overaccumulation in Pi-deprived frd3 roots (Fig. 4), which indicates that citrate secretion is required for proper Fe-distribution under Pi limitation. Interestingly, citrate application transiently promoted frd3 root growth in low Pi (Additional file 16: Figure S5), indicating that the Pi-dependent short root phenotype of frd3 is likely a consequence of altered Fe redistribution in the growing root.
Transcript analysis and spCCA (Additional file 1: Table S1, Additional file 13: Table S11) revealed regulation of ALMT genes, including a strong Pi-dependent induction of ALMT1 (Table 1), which was previously shown to exude malate into the rhizosphere as a strategy to cope with aluminum toxicity . Earlier studies revealed PHR1-dependent accumulation of malate and citrate in Pi-depleted plants [24, 69]. Interestingly, exudation of both carboxylates into the rhizosphere was shown to facilitate mobilization of Pi and Fe in several plant species that do not form mycorrhiza .
We also noticed deregulation of coumarin biosynthesis-related genes, F6’H1 and CCoAOMT1, in pdr2 roots (Additional file 1: Table S1, Additional file 9: Table S8). Several studies showed that coumarins (scopoletin and esculetin) are exuded into the rhizosphere to mobilize Fe in alkaline soils [49–51]. A recent report showed that esculetin accumulates in roots of Pi-starved wild-type plants but was suppressed in the phr1 mutant, which lacks the induction of PSR genes upon Pi deficiency . Moreover, using a non-targeted approach to identify metabolites from Pi-starved Arabidopsis root exudates, we recently confirmed Pi-dependent regulation of coumarin secretion . Thus, our analysis implicates additional Fe-chelators in the regulation of Pi-dependent Fe accumulation and/or distribution in roots.
Pi depletion modulates root pectins
Inhibition of root cell elongation, formation of root hairs and induction of lateral roots are the most robust local responses to Pi deficiency [7, 8], which all require extensive reorganization of the CW. Our analysis revealed Pi-dependent regulation of CW-modifying enzymes, particularly in the sensitive wild-type and pdr2 plants and to a lesser extent in lpr1lpr2 roots (Table 3). Consistent with a previous transcriptome study , we identified several putative pectin esterases and esterase inhibitors. Pectins are secreted into the apoplast in a highly methylesterified state. In the CW, pectin methylesterases (PME) may remove methyl groups, generating free carboxylate functions on the surface of pectin polymers. Crosslinking of these carboxylate-groups by Ca2+ reduces CW extensibility and regulates cell expansion . Our experiments revealed low Pi-induced accumulation of non-methylated pectin, specifically within the EZ of wild-type and pdr2 roots (Fig. 5), which might contribute to rapid inhibition of cell elongation in these lines. In addition, there is growing evidence that plants exchange Ca2+ ions for other divalent and trivalent metal ions to prevent metal uptake and ROS formation . Gessa et al.  showed in vitro Fe3+ binding to carboxylate groups on polygalacturonic acids (PGA), and two studies in Arabidopsis and rice demonstrated the ability of PGA to mobilize Pi from FePO4 complexes and clay [73, 74]. Interestingly, a decrease in pectins in the Arabidopsis qua1-2 mutant causes a hypersensitive short root phenotype upon Pi depletion . Here, we show that accumulation of pectin in the root meristem coincides with the sites of Fe accumulation (Fig. 5a, c). Local pectin deposition might be a strategy to mobilize Pi from Fe-phosphate complexes. The data support our previous observations of CW thickening and callose deposition at sites of Fe accumulation in the root tip .
A recent study of the Arabidopsis flower transcriptome revealed deregulation of PGIP1 and other CW-modifying enzymes in the ferritin1-3-4 triple mutant . PGIP1 is a member of the leucine-rich repeat (LRR) protein superfamily and inhibits fungal and bacterial polygalacturonases, which cleave non-methylated pectin residues in infected tissues . It further regulates germination by inhibiting the breakdown of seed coat pectins . Intriguingly, our analysis revealed co-regulation of PGIP1 and FER1 on transcript and protein level in all lines upon Pi-depletion (Table 4), further indicating a potential link between the Pi-dependent regulation of Fe distribution and the modification of pectin in the CW.
Peroxidases may modulate ROS formation and cell wall dynamics
We identified 41 CIII Prxs (56 % of the 73-member family) that were regulated on the mRNA and/or protein level, either in response to Pi depletion (23 members) or as a consequence of the pdr2 and lpr mutations (Additional file 7: Table S6, Additional file 10: Table S9). Interestingly, the majority of CIII Prx mRNAs/proteins (30) were deregulated in pdr2 in Pi replete conditions. CIII Prxs are involved in superoxide formation by transferring electrons from NADH to O2 as well as in the Fe catalyzed generation of hydroxyl radicals [76, 77]. ROS formation is likely responsible for the cleavage of CW polysaccharides to promote cell expansion. On the other hand, oxidation of monolignols by CIII Prxs is the predominant mechanism of monolignol polymerization (lignification) which rigidifies the CW and degrades H2O2 . The potential role of CIII Prxs for modulating ROS levels and CW dynamics and their strong deregulation in pdr2 mutants points to a function in local root growth adaptation. A comprehensive analysis of available transcriptome and proteome data revealed that most CIII Prxs are mainly expressed in the root . Two of those, Prx33 and Prx34, bind to Ca2+ polygalacturonates and mediate root growth in Arabidopsis . A more recent study demonstrated that prx33 and prx34 knock-down lines exhibited reduced ROS and callose formation upon treatment with microbe-associated molecular patterns (MAMPs), implicating a direct role of these gene products in ROS formation . Using specific ROS indicators, we recently demonstrated the formation of apoplastic ROS at the site of –Pi induced Fe deposition . The underlying mechanism remains elusive but CIII Prxs may constitute a missing link between Pi dependent ROS formation and CW remodeling.
Comparative transcriptome and proteome analysis allows in-depth dissection of gene expression
Our comparative transcriptome and proteome analysis revealed a highly significant but relatively low positive correlation for the abundance of PSR proteins and their cognate transcripts in all three genotypes tested (Fig. 3a, e). The majority of mRNA/protein pairs in our dataset showed discordant changes, which has been previously observed and discussed in Arabidopsis and other organisms like mice and humans and which is likely explained by (post-) translational regulation and/or a temporal delay between the regulation of transcript and protein abundance. In addition, technical limitations in the efficiency of protein identification (e.g., low abundant proteins and transmembrane proteins) may restrict the detection of proteins relative to their cognate transcripts [25, 81–84].
Correlation values significantly increased when gene activity was subcategorized. For example, we observed a strong positive correlation between protein and mRNA abundance when we focused on proteins that were Pi-responsive in all genotypes (Table 4). Similarly, we found an enhanced positive correlation when we compared only significantly regulated genes with their cognate proteins (Fig. 3b, c, d, e). Moreover, our observations suggest that the integration of transcriptome and proteome datasets can be used as a valuable complementary approach. For example, we identified 28 and 18 CIII Prx, regulated on the transcript and/or protein level, respectively. Only 5 of those showed correlative expression changes in both datasets (Additional file 7; Table S6, Additional file 10: Table S9). However, the integration of both approaches revealed regulation of 41 CIII Prx, suggesting that the majority of CIII Prx are involved in the response to Pi deprivation.
We demonstrate that spCCA is a useful tool to integrate all experimental factors in our investigation, including the proteome and transcriptome data, Pi-status and genotype in order to elucidate unknown correlations in this multidimensional dataset. Interestingly, the first two CVs of our spCCA indicated a prominent role of genes and proteins that were differentially regulated in Pi-replete pdr2 seedlings (Fig. 3f, g). Indeed, detailed analysis of our datasets revealed that the majority of Pi-responsive genes was not significantly deregulated in pdr2, compared to the wild type (Additional file 2: Figure S1C, Additional file 4: Table S3). On the other hand, several Fe-related genes, CIII Prx and pectin modifying enzymes were differentially regulated in Pi-replete pdr2 plants (Table 2, Table 3, Additional file 7: Table S6), indicating that conditional hypersensitivity in pdr2 might be a cause of constitutive de-repression or sensitization of these genes/proteins. P5-type ATPases are orphan, membrane localized ER proteins with unknown substrate specificity . Mutant studies on yeast SPF1 and Arabidopsis MIA/PDR2 strongly suggest a function in ER quality control, protein folding and regulation of secretory processes [13, 86–88]. Hyperaccumulation of pectin and callose in the CW of Pi-depleted pdr2 roots  (and this study) support a function of PDR2 in regulating ER-dependent secretion.
We performed complementary transcriptomics and proteomics approaches to monitor changes in steady-state transcript and protein levels upon Pi deprivation of Arabidopsis wild-type, pdr2 and lpr1lpr2 roots. Our analysis reveals a set of genes and proteins that are involved in the regulation of Fe homeostasis, cell wall remodeling and ROS formation. We observed increased FER1 and decreased IRT1 expression in all genotypes, which are consistent with intracellular Fe accumulation and feed-back inhibited Fe uptake in Pi-depleted roots, respectively. Analysis of fer1-3-4, fer1-2-3-4 and irt1 mutants demonstrates that cellular Fe uptake and Fe storage in ferritin are not involved in Pi-dependent modulation of root growth. We provide evidence for the importance of apoplastic Fe redistribution to maintain root growth upon Pi-depletion and for a role of FRD3 in this process. Our data further reveal Pi-dependent regulation of cell wall-modifying enzyme expression and changes in the deposition of pectins in Pi-deprived roots. The high correlation between sites of Fe deposition and enhanced pectin accumulation suggests that pectins might be involved in Fe binding and/or Pi mobilization from Fe-P complexes.
Plant material and growth conditions
Arabidopsis thaliana accession Columbia (Col-0) and Col lines pdr2-1, lpr1-1lpr2-1, irt1-1, frd3-7, fer1-3-4 and fer1-2-3-4 were previously described [11, 13, 58, 89, 90]. The pdr2-1 mutant was identified and characterized by our group [12, 13, 19]. The irt1-1 (SALK_024525) and frd3-7 (SALK_122235) lines were obtained from the European Arabidopsis Stock Center (NASC). The lpr1-1lpr2-1 double mutant and the ferritin mutants (fer1-(2)-3-4) were kindly provided by T. Desnos  and J.F. Briat , respectively. Seeds were surface-sterilized and germinated on 0.8 % (w/v) Phyto-Agar (Duchefa) containing 50 μM Fe-EDTA and 2.5 mM KH2PO4, pH 5.6 (high or + Pi medium) or no Pi supplement (low or –Pi medium) as reported [13, 19].
Root growth measurement
The position of the root tip was marked on the back of the agar plate directly after seedling transfer from + Pi to + Pi or –Pi medium. Images were taken on a stereomicroscope and total increment of primary root length was calculated at the according time point using ImageJ software. For daily growth rate measurements, the root tip position was marked every 24 h. The distance between two marker-points defines the daily root growth.
Accumulation and distribution of Fe and callose in roots was monitored as previously described . De-methyl esterified pectins were stained for 5–10 min in 0.05 % (w/v) Ruthenium Red solution (Applichem). Hydroxylamine-ferric chloride staining was adapted from Hornatowska and Reeve [63, 64]. Seedlings were initially incubated for 5–10 min in freshly prepared hydroxylamine solution (0.7 % NaOH, 0.7 % hydroxylamine hydrochloride in 60 % EtOH), followed by the addition of an equal (or higher) volume of a solution containing concentrated HCl/EtOH 95 % (1:2 ratio). The solution was removed and ferric chloride was added (10 % FeCl3 in 60 % EtOH containing 0.1 N HCl). Seedlings were cleared using chloral hydrate solution (7:7:1 chloral hydrate:ddH2O:glycerol). Samples were analyzed using a multizoom stereomicroscope (Nikon AZ100) for overview images and a Zeiss AxioImager bright field microscope for detail images.
RNA preparation and microarray hybridization
Seedlings (4-days-old) were transferred from + Pi to either + Pi or –Pi medium and roots were harvested after 20 h. RNA was extracted using the RNeasy Plant Mini Kit from Qiagen followed by an on-column DNA digestion (40 min) using Qiagen RNase-free DNase Set. Quality control and hybridization to ATH1 Arabidopsis GeneChips was done by NASC’s Affymetrix Service (http://affymetrix.arabidopsis.info/).
Statistical analysis of mRNA expression data
Data preprocessing, generation of Venn diagrams and heat maps was performed using Arraystar 4.1 software (DNASTAR). Arrays were normalized with robust multiarray analysis (RMA) and quantile background correction. Pairwise comparisons were performed using a fold-change cutoff value of ≥ 1.5 for increased and of ≤0.66 for decreased transcript levels (p ≤ 0.05; Student’s t-test, no multiple testing correction). Gene ontology analysis was done with the preassigned settings of the Arraystar software using a cutoff value p ≤ 0.05 and FDR (Benjamini Hochberg) correction. Hierarchical clustering was performed with 4870 ANOVA-filtered genes using the hclust package of the R software v.3.0.0. . The ATTED-II database (http://atted.jp) was used to generate a list of IRT1 co-regulated genes based on ATTED’s mutual ranking. All other calculations and graphics were prepared using Microsoft Excel 2010 software.
Preparation of protein samples and LC-MS analysis
Plants were grown as for mRNA analysis. Proteins were extracted from root tissue and digested with trypsin. Peptides were injected into an EASY-nLC II nano liquid chromatography system, equipped with a Nanospray Flex ion source (Thermo Fisher Scientific) and electrosprayed into an Orbitrap Velos Pro mass spectrometer (Thermo Fisher Scientific). Details are described in Additional file 17.
Protein identification and relative quantification
The raw data was imported into Proteome Discoverer v.1.4 (PD). Peak lists generated with a precursor signal to noise ratio of 1.5 with PD were used to search the TAIR10 database amended with common contaminants (35,394 sequences, 14,486,974 residues) with the Mascot algorithm v.2.5 on an in-house Mascot server. The enzyme specificity was set to trypsin and two missed cleavages were tolerated. Carbamidomethylation of cysteine was set as a fixed modification and oxidation of methionine as a variable modification. The precursor tolerance was set to 7 ppm and the product ion mass tolerance was set to 0.8 Da. A decoy database search was performed to determine the peptide false discovery rate (FDR). The search results were imported into the Scaffold Q+ software v.4.1.1 (Proteome Software, Inc.). Peptide and protein FDRs were calculated and the identity thresholds set to 0.01 and 1 % respectively to control the family wise error rate of peptide and protein identifications.
The raw data was imported into Progenesis LC-MS v.4.1 (Nonlinear Dynamics) for relative protein quantification between LC-MS analyses. The peptide ion signal peak landscapes of LC-MS analyses were aligned using the analysis as a reference that gave the highest minimum and maximum number of vectors in the aligned set of analyses when each analysis was used as a reference. Ratiometric normalization in log space to a selected reference analysis over all aligned peptide ion signals was performed. The summed intensities of peptide ion signal peak isotope envelopes over time were used as a measure of peptide abundance. A coefficient of variance (CV) of peptide abundance of less than 50 % for a peptide in all LC-MS analyses of a biological condition (three replicate analyses of each of three biological replicates for a total of 9) was required for a peptide to be quantified. Protein abundance was inferred by the sum of all unique peptides mapping to a given protein (non-conflicting peptides). Protein abundance fold changes and corresponding p-values between the biological conditions were calculated.
Multidimensional scaling (MDS) analysis
Multidimensional scaling was conducted using the isoMDS function from the MASS package version 7.3-29 . Technical replicates of the proteome analysis were averaged, reducing the original dataset to 18 biological replicate samples. Missing values were either imputed by half of the minimum intensity or excluded from further analysis. The resulting matrix of 3849x18 proteins was subjected to ANOVA (p < 0.05) revealing 412 consistent proteins. Intensities were log-transformed.
Supervised penalized canonical correlation analysis (spCCA)
SpCCA analysis was done according to . ANOVA filtered transcriptome and proteome data sets were reduced to signals with a variance of ≥0.12 and ≥0.4 resulting in 1143 transcripts and 47 proteins. The experimental design consisted of a binaric matrix of 18 samples x 8 experimental factors (three genotypes: Col, pdr2, lpr1lpr2; two growth media: +Pi, −Pi agar; and three replicates). SpCCA was conducted with 25 resampling runs (n.r = 25) and 25 random start vectors (max.counter.test = 25) to optimize sparsity parameters in a grid search between (0,0,0) and (0.6,0.5,1) with small step sizes (0.05,0.05,0.1) for transcriptomics, proteomics and design dataset.
Ethics (and consent to participate)
Consent for publication
Availability of data and materials
Microarray data sets with the reference number NASCARRAYS-648 were deposited on the NASCArrays database (http://affymetrix.arabidopsis.info/). The proteomics data have been deposited to the ProteomeXchange Consortium  via the PRIDE partner repository with the dataset identifier PXD003449 and 10.6019/PXD003449 (http://www.ebi.ac.uk/pride/archive/).
We thank J.F. Briat for ferritin mutant seeds. Research at the Leibniz Institute of Plant Biochemistry was supported by institutional core funding provided by the state of Saxony-Anhalt and the Federal Republic of Germany. The publication of this article was funded by the Open Access fund of the Leibniz Association.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Lin WY, Huang TK, Leong SJ, Chiou TJ. Long-distance call from phosphate: systemic regulation of phosphate starvation responses. J Exp Bot. 2014;65(7):1817–27.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Bari R, Datt Pant B, Stitt M, Scheible WR. PHO2, microRNA399, and PHR1 define a phosphate-signaling pathway in plants. Plant Physiol. 2006;141(3):988–99.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Abel S. Phosphate sensing in root development. Curr Opin Plant Biol. 2011;14(3):303–9.View ArticlePubMedGoogle Scholar
- Hinsinger P, Herrmann L, Lesueur D, Robin A, Trap J, Waithaisong K, Plassard C. Impact of roots, microorganisms and microfauna on the fate of soil phosphorus in the rhizosphere. In: Annual Plant Reviews. Volume 48. John Chichester, UK: Wiley & Sons, Inc.; 2015: 375–407.Google Scholar
- Reymond M, Svistoonoff S, Loudet O, Nussaume L, Desnos T. Identification of QTL controlling root growth response to phosphate starvation in Arabidopsis thaliana. Plant Cell Environ. 2006;29(1):115–25.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Ticconi CA, Delatorre CA, Lahner B, Salt DE, Abel S. Arabidopsis pdr2 reveals a phosphate-sensitive checkpoint in root development. Plant J. 2004;37(6):801–14.View ArticlePubMedGoogle Scholar
- Ticconi CA, Lucero RD, Sakhonwasee S, Adamson AW, Creff A, Nussaume L, Desnos T, Abel S. ER-resident proteins PDR2 and LPR1 mediate the developmental response of root meristems to phosphate availability. Proc Natl Acad Sci U S A. 2009;106(33):14174–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Sanchez-Calderon L, Lopez-Bucio J, Chacon-Lopez A, Gutierrez-Ortega A, Hernandez-Abreu E, Herrera-Estrella L. Characterization of low phosphorus insensitive mutants reveals a crosstalk between low phosphorus-induced determinate root development and the activation of genes involved in the adaptation of Arabidopsis to phosphorus deficiency. Plant Physiol. 2006;140(3):879–89.View ArticlePubMedPubMed CentralGoogle Scholar
- Camacho-Cristobal JJ, Rexach J, Conejero G, Al-Ghazi Y, Nacry P, Doumas P. PRD, an Arabidopsis AINTEGUMENTA-like gene, is involved in root architectural changes in response to phosphate starvation. Planta. 2008;228(3):511–22.View ArticlePubMedGoogle Scholar
- Yu H, Luo N, Sun L, Liu D. HPS4/SABRE regulates plant responses to phosphate starvation through antagonistic interaction with ethylene signalling. J Exp Bot. 2012;63(12):4527–38.View ArticlePubMedPubMed CentralGoogle Scholar
- Gonzalez-Mendoza V, Zurita-Silva A, Sanchez-Calderon L, Sanchez-Sandoval ME, Oropeza-Aburto A, Gutierrez-Alanis D, Alatorre-Cobos F, Herrera-Estrella L. APSR1, a novel gene required for meristem maintenance, is negatively regulated by low phosphate availability. Plant Sci. 2013;205–206:2–12.View ArticlePubMedGoogle Scholar
- Karthikeyan AS, Jain A, Nagarajan VK, Sinilal B, Sahi SV, Raghothama KG. Arabidopsis thaliana mutant lpsi reveals impairment in the root responses to local phosphate availability. Plant Physiol Biochem. 2014;77:60–72.View ArticlePubMedGoogle Scholar
- Müller J, Toev T, Heisters M, Teller J, Moore KL, Hause G, Dinesh DC, Bürstenbinder K, Abel S. Iron-dependent callose deposition adjusts root meristem maintenance to phosphate availability. Dev Cell. 2015;33(2):216–30.View ArticlePubMedGoogle Scholar
- Ward JT, Lahner B, Yakubova E, Salt DE, Raghothama KG. The effect of iron on the primary root elongation of Arabidopsis during phosphate deficiency. Plant Physiol. 2008;147(3):1181–91.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Zheng L, Huang F, Narsai R, Wu J, Giraud E, He F, Cheng L, Wang F, Wu P, Whelan J et al. Physiological and transcriptome analysis of iron and phosphorus interaction in rice seedlings. Plant Physiol. 2009;151(1):262–74.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Lan P, Li W, Schmidt W. Complementary proteome and transcriptome profiling in phosphate-deficient arabidopsis roots reveals multiple levels of gene regulation. Mol Cell Proteomics. 2012;11(11):1156–66.View ArticlePubMedPubMed CentralGoogle Scholar
- Woo J, Macpherson CR, Liu J, Wang H, Kiba T, Hannah MA, Wang XJ, Bajic VB, Chua NH. The response and recovery of the Arabidopsis thaliana transcriptome to phosphate starvation. BMC Plant Biol. 2012;12:62.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang J, Lan P, Gao H, Zheng L, Li W, Schmidt W. Expression changes of ribosomal proteins in phosphate- and iron-deficient Arabidopsis roots predict stress-specific alterations in ribosome composition. BMC Genomics. 2013;14:783.View ArticlePubMedPubMed CentralGoogle Scholar
- Lin WD, Liao YY, Yang TJ, Pan CY, Buckhout TJ, Schmidt W. Coexpression-based clustering of Arabidopsis root genes predicts functional modules in early phosphate deficiency signaling. Plant Physiol. 2011;155(3):1383–402.View ArticlePubMedPubMed CentralGoogle Scholar
- del Pozo JC, Allona I, Rubio V, Leyva A, de la Pena A, Aragoncillo C, Paz-Ares. A type 5 acid phosphatase gene from Arabidopsis thaliana is induced by phosphate starvation and by some other types of phosphate mobilising/oxidative stress conditions. Plant J. 1999;19(5):579–89.Google Scholar
- Cheng Y, Zhou W, El Sheery NI, Peters C, Li M, Wang X, Huang J. Characterization of the Arabidopsis glycerophosphodiester phosphodiesterase (GDPD) family reveals a role of the plastid-localized AtGDPD1 in maintaining cellular phosphate homeostasis under phosphate starvation. Plant J. 2011;66(5):781–95.View ArticlePubMedGoogle Scholar
- Liu TY, Aung K, Tseng CY, Chang TY, Chen YS, Chiou TJ. Vacuolar Ca2+/H+ transport activity is required for systemic phosphate homeostasis involving shoot-to-root signaling in arabidopsis. Plant Physiol. 2011;156(3):1176–89.View ArticlePubMedPubMed CentralGoogle Scholar
- Vert G, Grotz N, Dedaldechamp F, Gaymard F, Guerinot ML, Briat JF, et al. IRT1, an Arabidopsis transporter essential for iron uptake from the soil and for plant growth. Plant Cell. 2002;14(6):1223–33.Google Scholar
- Thomine S, Vert G. Iron transport in plants: better be safe than sorry. Curr Opin Plant Biol. 2013;16(3):322–7.View ArticlePubMedGoogle Scholar
- Colangelo EP, Guerinot ML. The essential basic helix-loop-helix protein FIT1 is required for the iron deficiency response. Plant Cell. 2004;16(12):3400–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Buckhout TJ, Yang TJ, Schmidt W. Early iron-deficiency-induced transcriptional changes in Arabidopsis roots as revealed by microarray analyses. BMC Genomics. 2009;10:147.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang TJ, Lin WD, Schmidt W. Transcriptional profiling of the Arabidopsis iron deficiency response reveals conserved transition metal homeostasis networks. Plant Physiol. 2010;152(4):2130–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Palmer CM, Hindt MN, Schmidt H, Clemens S, Guerinot ML. MYB10 and MYB72 are required for growth under iron-limiting conditions. PLoS Genet. 2013;9(11):e1003953.View ArticlePubMedPubMed CentralGoogle Scholar
- Yuan Y, Wu H, Wang N, Li J, Zhao W, Du J, Wang D, Ling HQ. FIT interacts with AtbHLH38 and AtbHLH39 in regulating iron uptake gene expression for iron homeostasis in Arabidopsis. Cell Res. 2008;18(3):385–97.View ArticlePubMedGoogle Scholar
- Wang N, Cui Y, Liu Y, Fan H, Du J, Huang Z, Yuan Y, Wu H, Ling HQ. Requirement and functional redundancy of Ib subgroup bHLH proteins for iron deficiency responses and uptake in Arabidopsis thaliana. Mol Plant. 2013;6(2):503–13.View ArticlePubMedGoogle Scholar
- Cosgrove DJ. Growth of the plant cell wall. Nat Rev Mol Cell Biol. 2005;6(11):850–61.View ArticlePubMedGoogle Scholar
- Wolf S, Greiner S. Growth control by cell wall pectins. Protoplasma. 2012;249(2):169–75.View ArticleGoogle Scholar
- Francoz E, Ranocha P, Nguyen-Kim H, Jamet E, Burlat V, Dunand C. Roles of cell wall peroxidases in plant development. Phytochemistry. 2014;112:15–21.View ArticlePubMedGoogle Scholar
- Ferrari S, Vairo D, Ausubel FM, Cervone F, De Lorenzo G. Tandemly duplicated Arabidopsis genes that encode polygalacturonase-inhibiting proteins are regulated coordinately by different signal transduction pathways in response to fungal infection. Plant Cell. 2003;15(1):93–106.View ArticlePubMedPubMed CentralGoogle Scholar
- Kanai M, Nishimura M, Hayashi M. A peroxisomal ABC transporter promotes seed germination by inducing pectin degradation under the control of ABI5. Plant J. 2010;62(6):936–47.PubMedGoogle Scholar
- Chen ZH, Nimmo GA, Jenkins GI, Nimmo HG. BHLH32 modulates several biochemical and morphological processes that respond to Pi starvation in Arabidopsis. Biochem J. 2007;405(1):191–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Miura K, Lee J, Miura T, Hasegawa PM. SIZ1 controls cell growth and plant development in Arabidopsis through salicylic acid. Plant Cell Physiol. 2010;51(1):103–13.View ArticlePubMedGoogle Scholar
- Miura K, Rus A, Sharkhuu A, Yokoi S, Karthikeyan AS, Raghothama KG, Baek D, Koo YD, Jin JB, Bressan RA et al. The Arabidopsis SUMO E3 ligase SIZ1 controls phosphate deficiency responses. Proc Natl Acad Sci U S A. 2005;102(21):7760–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Miura K, Lee J, Gong Q, Ma S, Jin JB, Yoo CY, Miura T, Sato A, Bohnert HJ, Hasegawa PM. SIZ1 regulation of phosphate starvation-induced root architecture remodeling involves the control of auxin accumulation. Plant Physiol. 2011;155(2):1000–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Schmid NB, Giehl RF, Doll S, Mock HP, Strehmel N, Scheel D, Kong X, Hider RC, von Wiren N. Feruloyl-CoA 6’-Hydroxylase1-dependent coumarins mediate iron acquisition from alkaline substrates in arabidopsis. Plant Physiol. 2014;164(1):160–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Rodriguez-Celma J, Lin WD, Fu GM, Abadia J, Lopez-Millan AF, Schmidt W. Mutually exclusive alterations in secondary metabolism are critical for the uptake of insoluble iron compounds by Arabidopsis and medicago truncatula. Plant Physiol. 2013;162(3):1473–85.View ArticlePubMedPubMed CentralGoogle Scholar
- Fourcroy P, Siso-Terraza P, Sudre D, Saviron M, Reyt G, Gaymard F, Abadia A, Abadia J, Alvarez-Fernandez A, Briat JF. Involvement of the ABCG37 transporter in secretion of scopoletin and derivatives by Arabidopsis roots in response to iron deficiency. New Phytol. 2014;201(1):155–67.View ArticlePubMedGoogle Scholar
- Kai K, Mizutani M, Kawamura N, Yamamoto R, Tamai M, Yamaguchi H, Sakata K, Shimizu B. Scopoletin is biosynthesized via ortho-hydroxylation of feruloyl CoA by a 2-oxoglutarate-dependent dioxygenase in Arabidopsis thaliana. Plant J. 2008;55(6):989–99.View ArticlePubMedGoogle Scholar
- Thum A, Mönchgesang S, Westphal L, Lübken T, Rosahl S, Neumann S, Posch S. Supervised Penalized Canonical Correlation Analysis. arXiv. 2014;preprint arXiv:1405.1534.Google Scholar
- Jain A, Wilson GT, Connolly EL. The diverse roles of FRO family metalloreductases in iron and copper homeostasis. Front Plant Sci. 2014;5:100.View ArticlePubMedPubMed CentralGoogle Scholar
- Curie C, Cassin G, Couch D, Divol F, Higuchi K, Le Jean M, Misson J, Schikora A, Czernic P, Mari S. Metal movement within the plant: contribution of nicotianamine and yellow stripe 1-like transporters. Ann Bot. 2009;103(1):1–11.View ArticlePubMedPubMed CentralGoogle Scholar
- Reyt G, Boudouf S, Boucherez J, Gaymard F, Briat JF. Iron and ferritin dependent ROS distribution impact Arabidopsis root system architecture. Mol Plant. 2014;8(3):439–53.View ArticlePubMedGoogle Scholar
- Durrett TP, Gassmann W, Rogers EE. The FRD3-mediated efflux of citrate into the root vasculature is necessary for efficient iron translocation. Plant Physiol. 2007;144(1):197–205.View ArticlePubMedPubMed CentralGoogle Scholar
- Roschzttardtz H, Seguela-Arnaud M, Briat JF, Vert G, Curie C. The FRD3 citrate effluxer promotes iron nutrition between symplastically disconnected tissues throughout Arabidopsis development. Plant Cell. 2011;23(7):2725–37.View ArticlePubMedPubMed CentralGoogle Scholar
- Green LS, Rogers EE. FRD3 controls iron localization in Arabidopsis. Plant Physiol. 2004;136(1):2523–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Rogers EE, Guerinot ML. FRD3, a member of the multidrug and toxin efflux family, controls iron deficiency responses in Arabidopsis. Plant Cell. 2002;14(8):1787–99.View ArticlePubMedPubMed CentralGoogle Scholar
- Sterling C. Crystal-structure of ruthenium red and stereochemistry of its pectic stain. Am J Bot. 1970;57(2):172–5.View ArticleGoogle Scholar
- Krishnamurthy K. Methods in cell wall cytochemistry. Boca Raton, Florida, USA: CRC press; 1999.Google Scholar
- Reeve RM. A specific hydroxylamine-ferric chloride reaction for histochemical localization of pectin. Stain Technol. 1959;34(4):209–11.View ArticlePubMedGoogle Scholar
- Hornatowska J. Visualisation of pectins and proteins by microscopy. 2005. STFI-Packforsk report.Google Scholar
- Ravet K, Touraine B, Boucherez J, Briat JF, Gaymard F, Cellier F. Ferritins control interaction between iron homeostasis and oxidative stress in Arabidopsis. Plant J. 2009;57(3):400–12.View ArticlePubMedGoogle Scholar
- Bournier M, Tissot N, Mari S, Boucherez J, Lacombe E, Briat JF, Gaymard F. Arabidopsis ferritin 1 (AtFer1) gene regulation by the phosphate starvation response 1 (AtPHR1) transcription factor reveals a direct molecular link between iron and phosphate homeostasis. J Biol Chem. 2013;288(31):22670–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Meyer S, De Angeli A, Fernie AR, Martinoia E. Intra- and extra-cellular excretion of carboxylates. Trends Plant Sci. 2010;15(1):40–7.View ArticlePubMedGoogle Scholar
- Kobayashi Y, Hoekenga OA, Itoh H, Nakashima M, Saito S, Shaff JE, Maron LG, Pineros MA, Kochian LV, Koyama H. Characterization of AtALMT1 expression in aluminum-inducible malate release and its role for rhizotoxic stress tolerance in Arabidopsis. Plant Physiol. 2007;145(3):843–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Pant BD, Pant P, Erban A, Huhman D, Kopka J, Scheible WR. Identification of primary and secondary metabolites with phosphorus status-dependent abundance in Arabidopsis, and of the transcription factor PHR1 as a major regulator of metabolic changes during phosphorus limitation. Plant Cell Environ. 2015;38(1):172–87.View ArticlePubMedGoogle Scholar
- Ziegler J, Schmidt S, Chutia R, Muller J, Bottcher C, Strehmel N, Scheel D, Abel S. Non-targeted profiling of semi-polar metabolites in Arabidopsis root exudates uncovers a role for coumarin secretion and lignification during the local response to phosphate limitation. J Exp Bot. 2016;67(5):1421–32.View ArticlePubMedPubMed CentralGoogle Scholar
- Krzeslowska M. The cell wall in plant cell response to trace metals: polysaccharide remodeling and its role in defense strategy. Acta Physiol Plant. 2011;33(1):35–51.View ArticleGoogle Scholar
- Gessa C, Deiana S, Premoli A, Ciurli A. Redox activity of caffeic acid towards iron(III) complexed in a polygalacturonate network. Plant Soil. 1997;190(2):289–99.View ArticleGoogle Scholar
- Nagarajah S, Posner AM, Quirk JP. Competitive adsorption of phosphate with polygalacturonate and other organic anions on kaolinite and oxide surfaces. Nature. 1970;228(5266):83–5.View ArticlePubMedGoogle Scholar
- Zhu XF, Wang ZW, Wan JX, Sun Y, Wu YR, Li GX, Shen RF, Zheng SJ. Pectin enhances rice (Oryza sativa) root phosphorus remobilization. J Exp Bot. 2015;66(3):1017–24.View ArticlePubMedGoogle Scholar
- Sudre D, Gutierrez-Carbonell E, Lattanzio G, Rellan-Alvarez R, Gaymard F, Wohlgemuth G, Fiehn O, Alvarez-Fernandez A, Zamarreno AM, Bacaicoa E et al. Iron-dependent modifications of the flower transcriptome, proteome, metabolome, and hormonal content in an Arabidopsis ferritin mutant. J Exp Bot. 2013;64(10):2665–88.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen SX, Schopfer P. Hydroxyl-radical production in physiological reactions. A novel function of peroxidase. Eur J Biochem. 1999;260(3):726–35.View ArticlePubMedGoogle Scholar
- Passardi F, Penel C, Dunand C. Performing the paradoxical: how plant peroxidases modify the cell wall. Trends Plant Sci. 2004;9(11):534–40.View ArticlePubMedGoogle Scholar
- Marjamaa K, Kukkola EM, Fagerstedt KV. The role of xylem class III peroxidases in lignification. J Exp Bot. 2009;60(2):367–76.View ArticlePubMedGoogle Scholar
- Passardi F, Tognolli M, De Meyer M, Penel C, Dunand C. Two cell wall associated peroxidases from Arabidopsis influence root elongation. Planta. 2006;223(5):965–74.View ArticlePubMedGoogle Scholar
- Daudi A, Cheng Z, O’Brien JA, Mammarella N, Khan S, Ausubel FM, Bolwell GP. The apoplastic oxidative burst peroxidase in Arabidopsis is a major component of pattern-triggered immunity. Plant Cell. 2012;24(1):275–87.View ArticlePubMedPubMed CentralGoogle Scholar
- Robles MS, Cox J, Mann M. In-vivo quantitative proteomics reveals a key contribution of post-transcriptional mechanisms to the circadian regulation of liver metabolism. PLoS Genet. 2014;10(1):e1004047.View ArticlePubMedPubMed CentralGoogle Scholar
- Baerenfaller K, Massonnet C, Walsh S, Baginsky S, Buhlmann P, Hennig L, Hirsch-Hoffmann M, Howell KA, Kahlau S, Radziejwoski A et al. Systems-based analysis of Arabidopsis leaf growth reveals adaptation to water deficit. Mol Syst Biol. 2012;8:606.View ArticlePubMedPubMed CentralGoogle Scholar
- Wilhelm M, Schlegl J, Hahne H, Moghaddas Gholami A, Lieberenz M, Savitski MM, Ziegler E, Butzmann L, Gessulat S, Marx H et al. Mass-spectrometry-based draft of the human proteome. Nature. 2014;509(7502):582–7.View ArticlePubMedGoogle Scholar
- Schwanhausser B, Busse D, Li N, Dittmar G, Schuchhardt J, Wolf J, Chen W, Selbach M. Global quantification of mammalian gene expression control. Nature. 2011;473(7347):337–42.View ArticlePubMedGoogle Scholar
- Palmgren MG, Nissen P. P-type ATPases. Annu Rev Biophys. 2011;40:243–66.View ArticlePubMedGoogle Scholar
- Cronin SR, Rao R, Hampton RY. Cod1p/Spf1p is a P-type ATPase involved in ER function and Ca2+ homeostasis. J Cell Biol. 2002;157(6):1017–28.View ArticlePubMedPubMed CentralGoogle Scholar
- Vashist S, Frank CG, Jakob CA, Ng DT. Two distinctly localized p-type ATPases collaborate to maintain organelle homeostasis required for glycoprotein processing and quality control. Mol Biol Cell. 2002;13(11):3955–66.View ArticlePubMedPubMed CentralGoogle Scholar
- Jakobsen MK, Poulsen LR, Schulz A, Fleurat-Lessard P, Moller A, Husted S, Schiott M, Amtmann A, Palmgren MG. Pollen development and fertilization in Arabidopsis is dependent on the MALE GAMETOGENESIS IMPAIRED ANTHERS gene encoding a type V P-type ATPase. Genes Dev. 2005;19(22):2757–69.View ArticlePubMedPubMed CentralGoogle Scholar
- Nishida S, Tsuzuki C, Kato A, Aisu A, Yoshida J, Mizuno T. AtIRT1, the primary iron uptake transporter in the root, mediates excess nickel accumulation in Arabidopsis thaliana. Plant Cell Physiol. 2011;52(8):1433–42.View ArticlePubMedGoogle Scholar
- Ravet K, Touraine B, Kim SA, Cellier F, Thomine S, Guerinot ML, Briat JF, Gaymard F. Post-translational regulation of AtFER2 ferritin in response to intracellular iron trafficking during fruit development in Arabidopsis. Mol Plant. 2009;2(5):1095–106.View ArticlePubMedGoogle Scholar
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5(10):R80.View ArticlePubMedPubMed CentralGoogle Scholar
- Venables W, Ripley BD. Modern applied statistics with S. New York, USA: Springer-Verlag; 2002.Google Scholar
- Vizcaíno JA, Deutsch EW, Wang R, Csordas A, Reisinger F, Ríos D, Dianes JA, Sun Z, Farrah T, Bandeira N, Binz PA,Xenarios I, Eisenacher M, Mayer G, Gatto L, Campos A, Chalkley RJ, Kraus HJ, Albar JP, Martinez-Bartolomé S, Apweiler R,Omenn GS, Martens L, Jones AR, Hermjakob H. ProteomeXchange provides globally co-ordinated proteomics data submission and dissemination. Nat Biotechnol. 2014;30(3):223–6.View ArticleGoogle Scholar