Post-Transcriptional Coordination of the Arabidopsis Iron Deficiency Response is Partially Dependent on the E3 Ligases RING DOMAIN LIGASE1 (RGLG1) and RING DOMAIN LIGASE2 (RGLG2)*

Acclimation to changing environmental conditions is mediated by proteins, the abundance of which is carefully tuned by an elaborate interplay of DNA-templated and post-transcriptional processes. To dissect the mechanisms that control and mediate cellular iron homeostasis, we conducted quantitative high-resolution iTRAQ proteomics and microarray-based transcriptomic profiling of iron-deficient Arabidopsis thaliana plants. A total of 13,706 and 12,124 proteins was identified with a quadrupole-Orbitrap hybrid mass spectrometer in roots and leaves, respectively. This deep proteomic coverage allowed accurate estimates of post-transcriptional regulation in response to iron deficiency. Similarly regulated transcripts were detected in only 13% (roots) and 11% (leaves) of the 886 proteins that differentially accumulated between iron-sufficient and iron-deficient plants, indicating that the majority of the iron-responsive proteins was post-transcriptionally regulated. Mutants harboring defects in the RING DOMAIN LIGASE1 (RGLG1)1 and RING DOMAIN LIGASE2 (RGLG2) showed a pleiotropic phenotype that resembled iron-deficient plants with reduced trichome density and the formation of branched root hairs. Proteomic and transcriptomic profiling of rglg1 rglg2 double mutants revealed that the functional RGLG protein is required for the regulation of a large set of iron-responsive proteins including the coordinated expression of ribosomal proteins. This integrative analysis provides a detailed catalog of post-transcriptionally regulated proteins and allows the concept of a chiefly transcriptionally regulated iron deficiency response to be revisited. Protein data are available via ProteomeXchange with identifier PXD002126.

Acclimation to changing environmental conditions is mediated by proteins, the abundance of which is carefully tuned by an elaborate interplay of DNA-templated and post-transcriptional processes. To dissect the mechanisms that control and mediate cellular iron homeostasis, we conducted quantitative high-resolution iTRAQ proteomics and microarray-based transcriptomic profiling of iron-deficient Arabidopsis thaliana plants. A total of 13,706 and 12,124 proteins was identified with a quadrupole-Orbitrap hybrid mass spectrometer in roots and leaves, respectively. This deep proteomic coverage allowed accurate estimates of post-transcriptional regulation in response to iron deficiency. Similarly regulated transcripts were detected in only 13% (roots) and 11% (leaves) of the 886 proteins that differentially accumulated between iron-sufficient and iron-deficient plants, indicating that the majority of the iron-responsive proteins was post-transcriptionally regulated. Mutants harboring defects in the RING DOMAIN LIGASE1 (RGLG1) 1  With the exception of some bacteria that substitute manganese for iron or that rely on obligate parasitism (1,2), iron is an essential element for all organisms, catalyzing numerous redox reactions by virtue of its unique electrochemical properties. The cellular iron concentration must be carefully balanced. An inadequate supply compromises essential cellular processes, whereas excess amounts trigger the formation of toxic reactive oxygen species through Fenton chemistry. Iron deficiency causes substantial yield losses and decreases the nutritional quality of crop plants. Low iron levels in edible plant parts are the major cause of iron deficiency-induced anemia, affecting approximately two billion people worldwide particularly in areas where plant-based diets are the major source of iron (3).
Iron shortage leads to pronounced changes in transcriptomic and proteomic profiles and causes reprogramming of metabolic and developmental pathways, which are aimed at improving iron acquisition and distribution (4 -8). In Arabidopsis, massive up-regulation of the phenylpropanoid and scopoletin pathways leads to increased production of iron-binding coumarins that chelate poorly soluble Fe 3ϩ oxyhyrates and thereby improve the acquisition of iron from pools of sparingly soluble iron (4, 9 -11). Induction of the proton-translocating (H ϩ )-ATPASE2 (AHA2) decreases the rhizospheric pH, thereby increasing the solubility of iron by a factor of 1000 for each pH unit (12). Fe 3ϩ chelates are reduced by FERRIC REDUCTION OXIDASE2 (FRO2), and the released Fe 2ϩ is taken up by IRON-REGULATED TRANSPORTER1 (IRT1) (13)(14)(15). In leaves, a rapid coordinated down-regulation of the protoporphyrin biosynthetic pathway and up-regulation of the expression of proteins involved in the detoxification of reac-tive oxygen species avoid the accumulation of potentially toxic tetrapyrrole intermediates and prevent photo-oxidative damage in anticipation of decreased chlorophyll production and diminished photosynthesis rates that occur during prolonged periods of iron deficiency.
Reciprocal grafting experiments with the pea mutant dgl, which displays constitutive expression of iron deficiency responses, indicate that the iron status of the leaves dictates the iron uptake rates of roots (16). In Arabidopsis, OLIGOPEP-TIDE TRANSPORTER3 (OPT3) is critical in shoot-to-root signaling by mediating the loading of iron into the phloem (17)(18)(19). In roots, the expression of the key genes mediating the acquisition of iron from the rhizosphere, namely IRT1, FRO2, AHA2, and the feruloyl CoA ortho-hydroxylase 1 (F6ЈH1), which converts feruloyl CoA into 6Ј-hydroxyferuloyl CoA, a key step in the production of iron-binding compounds, are transcriptionally regulated by the bHLH transcription factor FER-LIKE IRON DEFICIENCY-INDUCED TRANSCRIPTION FACTOR (FIT). FIT forms a heterodimer with the subgroup 1b bHLH proteins bHLH38, bHLH39, and probably also with bHLH100, and bHLH101 (20 -26). Another bHLH transcription factor, POPEYE (PYE), is negatively regulating a separate set of genes by interacting with ILR3 and bHLH115 (27). The ubiquitin ligase BRUTUS (BTS) also interacts with IAA-LEU-CINE RESISTANT3 (ILR3) and bHLH115 and has opposite effects on gene expression, thereby balancing PYE action (27,28).
Limited information is available regarding post-transcriptional processes that control cellular iron homeostasis. A large set of genes is differentially alternatively spliced in response to iron deficiency, probably as a means to tune the amount of functional proteins to changing demand (29). The activity of two key players in iron uptake, FIT and IRT1, has been shown to be controlled both transcriptionally and by post-transcriptional processes. Proteolytic turnover and stabilization of FIT by interaction with the ethylene signaling transcription factors ETHYLENE-INSENSITIVE3 (EIN3) and ETHYLENE-INSENSI-TIVE3-LIKE (EIL1) is critical for FIT activity, which is rapidly degraded after binding to its targets (30,31). Monoubiquitination and endocytic cycling are important for IRT1 function, calibrating the amount of IRT1 protein on the plasma membrane (32). Endocytic sorting of IRT1 is dependent on SORT-ING NEXIN1, a component in a retromer-like protein complex involved in endosome to lysosome protein transport (33). A RING-type ubiquitin ligase, IRT1 DEGRADATION FACTOR1, is involved in the proteolytic turnover of IRT1 (34).
The ubiquitin conjugase UBC13 is required for the formation of branched root hairs, a response of Arabidopsis to iron deficiency that has been interpreted as a strategy to increase the surface area of the roots (35,36). The sequence of UBC13 is conserved among eukaryotes and its function has been related to the error tolerance branch of the DNA repair pathway in yeast, mammals, and plants. In Arabidopsis, UBC13 is encoded by two highly similar paralogs, UBC13A and UBC13B (37). ubc13a mutants fail to induce branched root hairs in response to iron deficiency and show deregulation of several iron-responsive genes (35). ubc13a ubc13b double mutants display severe defects in root hair formation also under iron-sufficient conditions (35,38), indicating that functional UBC13 is critical in this process. UBC13 is the only known protein that can mediate the formation of ubiquitin chains linked to lysine 63 (K63) (39). The RING domain ligase RGLG1 and its close sequelog RGLG2 can interact with UBC13 and, together with UBC13, catalyze the formation of K63-linked polyubiquitin chains (40). RGLG1 and RGLG2 possess ubiquitin ligase activity and can also mediate the formation of canonical, K48-linked polyubiquitin chains that target proteins for degradation (41). Interestingly, rglg1 rglg2 double mutants show a constitutively branched root hair phenotype (35), which invited some speculation as to whether and how the double mutation can be linked to UBC13-mediated processes (42,43). It has been suggested that under iron deficiency UBC13 recruits RGLG to the nucleus, leading to a decrease in RGLG protein in the cytoplasm and ultimately to branching of the root hairs (43).
To dissect post-transcriptional responses of Arabidopsis to iron deficiency and to determine the possible impact of RGLG on the regulation of these responses, we conducted genomewide proteomic and transcriptomic surveys of leaves and roots from iron-sufficient and iron-deficient Col-0 wild-type plants and rglg1 rglg2 double mutants. This analysis showed that post-transcriptional regulation has a stronger influence on the proteomic readout than transcriptional control, affecting proteins that in turn control post-transcriptional processes. In particular, the expression of ribosomal proteins is strongly post-transcriptionally affected by iron deficiency, putatively leading to a bias in translation by prioritizing subsets of mRNAs that are critical to the acclimation to low iron availability. RGLG has a dramatic influence on the proteomic profile of iron-deficient plants, affecting both protein abundance and, most likely as a secondary affect, the transcriptional profile by targeting transcription factors and other proteins involved in transcriptional regulation. The combined analysis permits an integration of several regulatory layers involved in adapting plants to low iron availability and underscores the importance of protein turnover in this process.

MATERIALS AND METHODS
Plant Growth Conditions-Arabidopsis thaliana (L.) Heynh, Columbia (Col-0) ecotype was used as the wild-type control. rglg1 rglg2 mutants were kindly provided by A. Bachmair, University of Vienna. Plants were grown in a growth chamber on medium as described by Estelle and Somerville (44). Seeds were surface sterilized and germinated on a media containing, KNO 3 (5 mM), MgSO 4 (2 mM), Ca(NO 3 ) 2 (2 mM), KH 2 PO 4 (2.5 mM), H 3 (40 M), solidified with 0.4% Gelrite pure (Kelco), 1.5% sucrose and 1 g/L MES. The pH was adjusted to 5.5 with KOH. Seeds were sown on Petri plates and stratified for 1 day in 4°C in the dark before being transferred to a growth chamber and grown at 21°C under continu-ous illumination (50 mol m Ϫ2 s Ϫ1 ). After 10 d of precultivation, plants were transferred to fresh agar medium either with 40 M FeEDTA (ϩFe plants) or without Fe and with 100 M 3-(2-pyridyl)-5,6-diphenyl-1,2,4-triazine sulfonate (-Fe plants). For iron quantification and determination of trichome density, plants were grown for 14 d on media supplemented with either 40 M or 0.5 M FeEDTA.
Protein Extraction-Roots from control and iron-deficient plants (13-day-old) were ground in liquid nitrogen and suspended in 10ϫ volume of precooled acetone (Ϫ20°C) containing 10% (w/v) TCA and 0.07% (v/v) 2-mercaptoethanol. Proteins were then thoroughly mixed and precipitated for 2 h at Ϫ20°C. Proteins were collected by centrifuging at 35,000 ϫ g (JA-20 108 rotor; Beckman J2-HS) at 4°C for 30 min. The supernatant was carefully removed, and the protein pellets were washed twice with cold acetone containing 0.07% (v/v) 2-mercaptoethanol and 1 mM phenylmethanesulfonyl fluoride and a third time with cold acetone without 2-mercaptoethanol. Protein pellets were dried by lyophilization and stored at Ϫ80°C or immediately extracted using protein extraction buffer composed of 8 M urea, 50 mM Tris, pH 8.5, for 1 h at 6°C under constant shaking. Protein extracts were centrifuged at 19,000 ϫ g for 20 min at 10°C. The supernatant was then collected, and the protein concentration was determined using a protein assay kit (Pierce).
In-Solution Trypsin Digestion and iTRAQ Labeling-Total protein (100 g) was reduced by adding dithiothreitol to a final concentration of 10 mM and incubated for 1 h at room temperature. Subsequently, iodoacetamide was added to a final concentration of 50 mM, and the mixture was incubated for 30 min at room temperature in the dark. Then, dithiothreitol (30 mM) was added to the mixture to consume any free iodoacetamide by incubating the mixture for 1 h at room temperature in the dark. Proteins were then diluted by 50 mM Tris, pH 8.5, to reduce the urea concentration to 4 M and digested with 0.5 g of Lys-C (Wako) for 4 h at room temperature. The Lys-C digested protein solution was further digested with 20 g of modified trypsin (Promega) at 37°C overnight after the solution was further diluted with 50 mM Tris, pH 8.0, to reduce the urea concentration to less than 1 M. The resulting peptide solution was acidified with 10% trifluoroacetic acid and desalted on a C18 solid-phase extraction cartridge.
Desalted peptides were then labeled with iTRAQ reagents (Applied Biosystems) according to the manufacturer's instructions. Control samples (proteins extracted from roots of control plants) were labeled with reagent 114; samples from iron-deficient roots were labeled with reagent 117. Three independent biological experiments with two technical repeats each were performed. The reaction was allowed to proceed for 1 h at room temperature. Subsequently, treated and control peptides were combined and further fractionated offline using high-resolution strong cation-exchange chromatography (Poly-Sulfoethyl A, 5 m, 200-Å beads). In total, 50 fractions were collected and combined into 16 final fractions. Each final fraction was lyophilized in a centrifugal speed vacuum concentrator. Samples were stored at Ϫ80°C.
LC-MS/MS Analysis-Liquid chromatography was performed on a Dionex UltiMate 3000 RSLCnano System coupled to a Q Exactive hybrid quadrupole-Orbitrap mass spectrometer (Thermo Scientific) equipped with a Nanospray Flex Ion Source. Peptide mixtures were loaded onto a 75 m ϫ 250 mm Acclaim PepMap RSLC column (Thermo Scientific) and separated using a segmented gradient in 120 min from 5 to 40% solvent B (100% acetonitrile with 0.1% formic acid) at a flow rate of 300 nl/min. Solvent A was 0.1% formic acid in water. The samples were maintained at 8°C in the autosampler. The LTQ Orbitrap was operated in the positive ion mode with the following acquisition cycle: a full scan (m/z 350ϳ1600) recorded in the Orbitrap analyzer at resolution R 70,000 was followed by MS/MS of the 10 most intense peptide ions with HCD acquisition of the same precursor ion. HDC was done with a collision energy of 30%. HCD-generated ions were detected in the Orbitrap at resolution 17,500.
Database Search-Two search algorithms, Mascot (version 2.4, Matrix Science) and SEQUEST, which is integrated in Proteome Discoverer software (version 1.4, Thermo Scientific), were used to simultaneously identify and quantify proteins. Searches were made against the Arabidopsis protein database (TAIR10 20110103, 27416 sequences; ftp://ftp.arabidopsis.org/home/tair/Sequences/blast_datasets/ TAIR10_blastsets/TAIR10 pep 20110103 representative gene model) concatenated with a decoy database containing the reversed sequences of the original database. The protein sequences in the database were searched with trypsin digestion at both ends and two missed cleavages allowed, fixed modifications of carbamidomethylation at Cys, iTRAQ 4plex at N terminus and Lys, variable modifications of oxidation at Met and iTRAQ 4plex at Tyr; peptide tolerance was set at 10 ppm, and MS/MS tolerance was set at 0.05 Da. iTRAQ 4plex was chosen for quantification during the search simultaneously. The search results were passed through additional filters, peptide confidence more than 95% (p Ͻ 0.05), before exporting the data. For protein quantitation, only unique peptides were used to quantify proteins. These filters resulted in a false discovery rate of less than 5% after decoy database searches were performed. For biological repeats, spectra from the two technical repeats were combined into one file and searched.
Statistical Analysis-Proteins identified and quantified in at least two biological repeats were considered to further analyze the abundance change in response to iron deficiency using a method described by Cox and Mann (45). In brief, the log 2 ratios of the 9,110 (roots) and 8,303 (leaves) quantified proteins overlapping in at least two biological repeats were calculated and analyzed for normal distribution. For a given protein in one biological repeat, the ratio was calculated as the inverse log 2 of the median of the log 2 value of all peptide ratios and averaged across the biological replicates. Next, mean and S.D. were calculated and 95% confidence (Z score ϭ 1.96) was used to select those proteins whose distribution was far from the main distribution. For the down-regulated proteins, the confidence interval was Ϫ0.475072001 (0.040469377, mean ratio of the 9,110 proteins in roots) Ϫ1.96 ϫ 0.263031316, corresponding to a protein ratio of 0.83. Similarly, for the up-regulated proteins, the confidence interval was calculated (mean ratio ϩ1.96ϫ S.D.), corresponding to a protein ratio of 1.29. Protein ratios outside this range were defined as being significantly different at p ϭ 0.05.
RNA Extraction-Total RNA was isolated from roots of 13-d-old plants with the RNeasy Plant Mini Kit (Qiagen) according to the manufacturer's instructions. Nucleic acid quantity was analyzed with a NanoDrop ND-1000 UV-Vis spectrophotometer (NanoDrop Technologies, Wilmington).
Microarray Analysis-The Affymetrix gene chip Arabidopsis ATH1 Genome Array was used for microarray analysis. Total RNA samples were prepared as described above. All RNA samples were quality assessed by using the Agilent Bioanalyzer 2100 (Agilent, Santa Clara, CA). Complementary RNA synthesis was performed by use of the GeneChip One-Cycle Target Labeling Kit (Affymetrix, Santa Clara, CA). Hybridization, washing, staining, and scanning procedures were performed as described in the Affymetrix technical manual.
Gene expression data were imported directly into GeneSpring (version 11.5, Agilent). The software was used to normalize the data per chip to the 50 th percentile and per gene to the control samples. Genes that were flagged as absent in two replicates were not considered in the analysis. p values for the Benjamini and Hochberg method (false discovery rates; FDRs) were calculated by GeneSpring. Transcripts were defined as differentially expressed that showed delta signal changes larger than the mean expression value of the whole data set or twofold changes upon iron deficiency with p Ͻ 0.05. Microscopy-For cross-sections, root samples were fixed, dehydrated, and then embedded in Technovit 7100 (Heraeus Kulzer, Wehrheim) resin in gelatin capsules. Transverse sections (30 m) were cut using a RM 2255 Leica microtome (Leica, Nussloch, Germany). Sections were dried and stained with toluidine blue (0.05%) on glass slides and examined using bright-field on an Imager Z1 microscope (Zeiss, Jena, Germany). Trichomes were observed on rosette leaves that were attached to the stage and cooled in liquid nitrogen. Observations were made in low vacuum with a scanning electron microscope (FEI Quanta 200).
Iron Concentration-Plants were dried in an oven at 60°C, mineralized with 225 l 65% nitric acid (HNO) at 96°C for 6 h, and oxidized with 150 l 30% hydrogen peroxide (H 2 O 2 ) at 56°C for 2 h. Iron concentrations were determined against a standard curve made with FeCl 3 that was treated in the same way as the plant materials. Samples were measured spectrophotometrically at A 535 nm in an assay solution that contained 1 mM bathophenanthroline disulfonic, 0.6 M sodium acetate and 0.48 M hydroxylamine hydrochloride.

Protein and mRNA Expression are Largely Uncoupled in
Iron-Deficient Plants-To dissect global transcriptional and post-transcriptional regulation of gene expression in iron-deficient plants, changes in proteomic and transcriptomic profiles were assessed by liquid chromatography combined with tandem mass spectroscopy on a Q Exactive hybrid quadrupole-Orbitrap mass spectrometer and by Affymetrix ATH1 microarray analysis. To gain insights into the post-transcriptional regulation of iron-responsive genes, responses of Col-0 wild-type plants were compared with mutants defective in the RING domain ligases RGLG1 and RGLG2. We previously reported that homozygous rglg1 rglg2 double mutants showed a constitutive branched root hair phenotype, a trait that is inducible in the wild type by subjecting plants to iron deficiency (35). It can thus be assumed that RGLG is involved in the regulation of the activity of iron-responsive genes, presumably acting as a post-translational regulator. It was shown previously that neither rglg1 nor rglg2 mutants produce full length RGLG transcript (40).
Proteins were extracted, digested in solution, iTRAQ-labeled and identified using the Mascot and SEQUEST search algorithms. With this technology, we identified a total of 13,706 and 12,124 proteins in roots and leaves of Col-0 wild-type plants, respectively. Peptides from 9,110 proteins in roots and 8,303 proteins in leaves were determined in at least two of the three biological replicates and were defined as reliably detected (Fig. 1A, 1B, supplemental Data Set S1, supplemental Tables S1 and S2).
Based on a 95% confidence (Z score ϭ 1.96), we defined subsets of 481 (roots) and 445 (leaves) reliably detected proteins as differentially expressed between iron-sufficient and iron-deficient plants (Fig. 1C). Only 40 proteins accumulated differentially in both roots and leaves, among them several proteins that have been previously associated with iron metabolism, including the Feruloyl CoA ortho-hydroxylase 1 (F6ЈH1), NICOTIANAMINE SYNTHASE4 (NAS4), IRON-RESPONSIVE PROTEIN3 (IRP3), a protein with unknown function that is highly expressed in leaves and roots of irondeficient plants (5), and the METAL TOLERANCE PROTEIN A2 (MTPA2) (46). Notably, of the 40 differentially expressed proteins in common between roots and leaves, 10 are structural constituents of ribosomes.
Subsets of 369 and 477 transcripts showed either delta signal changes larger than the mean expression value of the whole data set or twofold changes upon iron deficiency with p Ͻ 0.05 and were defined as differentially expressed. Sixtytwo mRNAs were iron-regulated in both roots and leaves ( Fig.  1D; supplemental Table S2). A comparison with previously published RNA-seq data sets for roots and leaves of plants grown under similar conditions (4, 5) showed a comparable fold-change of the differentially expressed genes despite the use of different technologies, showing the robustness of the current transcriptional analysis (Fig. 2). A comparison of the differentially expressed proteins in leaves and roots revealed that only a relatively small subset of mRNAs had and RGLG2 were considerably smaller than wild-type plants when grown on both iron-replete and low iron media (0.5 M FeEDTA; Fig. 3A). The rglg1 rglg2 double mutants formed leaves that were curled and more pointed than those of the wild type (Fig. 3B), a trait also observed in mutants defective in the expression of genes encoding ribosomal proteins (47,48). Constitutive root hair branching in rglg1 rglg2 plants was confirmed in the present study (Fig. 3C, 3D). In roots of rglg1 rglg2 plants, additional cell divisions were frequently observed (Fig. 3D) (35). In addition to the root hair phenotype, we observed a massive reduction in trichome density in rglg1 rglg2 plants (Fig. 3E-3G). Under control conditions, trichome density in rglg1 rglg2 plants was 84% lower when compared with wild-type plants. Iron deficiency significantly decreased the number of trichomes in the wild type (p Ͻ 0.05%, Fig. 3G), but did not further alter the trichome density in rglg1 rglg2 mutants (Fig. 3E).
The Lack of Functional RGLG Proteins has Far-Reaching Consequences on Global Gene Expression-A total of 443 and 1,351 transcripts were differentially expressed between the wild type and the mutant in roots and leaves, respectively (Fig. 4). Subsets of 285 and 445 proteins accumulated differentially between the two genotypes in roots and leaves. Notably, several transcripts that were iron-responsive in the wild type were differentially expressed in rglg1 rglg2 plants independent of the iron supply. In roots, this set comprises 124 genes and includes genes with well-defined roles in iron metabolism such as AHA2, PAL1, SAM2, and the superoxide dismutases FSD1 and CSD1. Also, the metallothionein MT2B was constitutively (up)-regulated in rglg1 rglg2 plants. In leaves, 289 iron-responsive genes were differentially expressed between Col-0 and rglg1 rglg2 plants under control conditions. Among them are several key genes involved in photosynthesis such as PSI/PSII light harvesting complexes, PSI/PSII subunits, as well as genes encoding proteins involved in plastid-to-nucleus signal transduction (GUN5) and protoporphyrin synthesis (PORB).

FIG. 2. Comparison of gene expression changes between microarray data (this study) and RNA-seq data take from (4) (roots) and (5) (leaves).
Similar to what has been observed at the transcript level, relatively large subsets of proteins that were iron-responsive in the wild type had a constitutive iron response in rglg1 rglg2 mutants grown under iron-sufficient conditions ( Fig. 4; supplemental Data Set S1, supplemental Tables S1 and S2).
Although some constitutively regulated proteins showed similarly increased or decreased abundance in iron-deficient wild type and iron-sufficient rglg1 rglg2 plants, several exceptions to this pattern were observed. In particular, in roots, several proteins accumulated massively in rglg1 rglg2 plants under control conditions whereas, similar to iron-deficient wild-type plants, they were down-regulated when mutant plants were subjected to iron deficiency. All three ferritin genes showed reduced expression in roots and leaves of rglg1 rglg2 plants at both the RNA and protein level, indicative of a higher iron status of rglg1 rglg2 relative to Col-0 plants in both organs and growth types.
Severe repression of two transcription factors in the leaves of rglg1 rglg2 plants, the class IV homeodomain leucine zipper GLABRA2 (GL2) and the R2R3 MYB-family protein MYB23, may account for the reduced trichome phenotype of the mutant. Both proteins are critical for proper trichome development (49,50) and act together in a positive feedback loop in which GL2 is required for transcriptional activation of MYB23 (51). In rglg1 rglg2 plants, expression of both genes was decreased by more than fourfold. In roots, the expression of three genes encoding proteins involved in cell wall loosening, the ␤-galactosidase glycosyl hydrolases BGAL1, BGAL4, and BGAL8 were strongly decreased in rglg1 rglg2 plants when compared with the wild type. BGAL4 and BGAL8 transcripts are enriched in root hairs (29) and showed decreased abundance in iron-deficient wild-type plants that also develop branched root hairs, albeit with a lower frequency (35,36). Another cell wall modifying enzyme, PECTIN METHYL-ESTERASE17 (PME17), was down-regulated in iron-sufficient rglg1 rglg2 and iron-deficient wild-type plants. Similar to BGAL4 and BGAL8, PME17 is predominantly expressed in root hairs (29). In addition, PME17 is part of a core set of root hair genes that have been compiled from genetic, genomic, and computational analyses (52) and are thus putatively involved in root hair morphogenesis. Decreased expression of PME17, possibly in combination with a lack of functional BGAL protein, may contribute to the branched root hair phenotype.
The Majority of Iron-Responsive Proteins is Post-Transcriptionally Regulated-Proteins that are differentially expressed under iron-deficient conditions can be divided into populations of chiefly transcriptionally regulated and mainly or entirely post-transcriptionally regulated proteins. This distinction is based on the detection of a cognate, differentially expressed transcript that is similarly regulated in the protein.
Transcription data are taken from either the current microarray analysis or from previous RNA-seq-based transcriptomic surveys conducted on plants grown under similar conditions (4, 5) that were surveyed to avoid false negatives because of ambiguous or absent probe sets in the ATH1 gene chip. It should be noted that the term "transcriptionally regulated" as it is used here does not exclude the possibility of additional Categorizing the responses in roots and leaves of wild-type plants at the protein level revealed that in both organs transcriptional regulation appears to be more diverse, i.e. more gene ontology (GO) categories are represented when compared with post-transcriptionally regulated proteins (Fig. 5). In roots, GOs related to cellular iron homeostasis such as "response to reactive oxygen species," "response to iron homeostasis," "iron transport," "response to nitric oxide," "ethylene stimulus," and "phenylpropanoid biosynthetic pathway" are represented by several proteins. Post-transcriptionally regulated proteins in roots comprise the GOs "DNA-templated transcription," "N-terminal myristoylation," and "ribosome biogenesis." In leaves, "response to cadmium ion," "ribosome biogenesis," and "protein folding" represent the GOs comprising most transcriptionally regulated proteins. Among the post-transcriptionally regulated proteins, "translation" and "ribosome biogenesis" are most dominant.
Fifty-seven and 51 similarly regulated transcripts and proteins were detected in roots and leaves of wild-type plants, respectively. These subsets account for 13 and 11% of the differentially expressed proteins, indicative of a pronounced post-transcriptional component in the "proteoferrome." To estimate the involvement of transcriptional regulation, we analyzed the concordance of mRNA and protein expression for genes in which both the mRNA and the protein accumulated differentially between iron-sufficient and iron-deficient plants.
Using the microarray-derived transcript data, in roots a Pearson correlation coefficient of r ϭ 0.79 was calculated (Fig. 6). A similar correlation was determined using the RNA-seq data set (r ϭ 0.81). In leaves, almost identical correlation coefficient values were observed with data derived from microarray (r ϭ 0.83) and RNA-seq analysis (r ϭ 0.81), indicating that changes in mRNA abundance account for ϳ80% of the changes in protein abundance in the wild type for transcriptionally regulated proteins. When all protein data were taken into account for which a differentially expressed transcript was detected regardless of the direction of regulation, Pearson correlation values for roots were r ϭ 0.62 and r ϭ 0.65 for ATH1 and RNA-seq data, respectively, those for leaves were r ϭ 0.74 and r ϭ 0.68 (Fig. 6).
The Expression of Key Proteins in Iron Homeostasis is Transcriptionally Regulated-In roots of the wild type, the group of transcriptionally regulated proteins is composed of several enzymes that mediate the uptake of iron (IRT1 and FRO2), proteins involved in the production of iron-binding phenylpropanoids (PAL1, PAL2, F6ЈH1, 4CL1, 4CL2, BGLU42), en- zymes mediating the biosynthesis of the iron chelator nicotianamine (NAS1, NAS4), iron storage proteins (FER1, FER3), and key players in the metabolism of iron-deficient plants (SAM1, MS2, G6PD6, PEPC1). Several of the genes encoding these proteins in roots are regulated by the transcription factor FIT (21). The transcriptionally regulated proteins contain 20 members of the "ferrome," a core set of iron-responsive genes inferred from public microarray data (53). The function of transcriptionally regulated iron-responsive proteins and the mechanisms controlling the expression of the corresponding genes in roots have been reviewed extensively and will not be further discussed here (26, 54 -65).
In leaves, the iron transporter OPT3, all three expressed ferritins (FER1, FER3, FER4), proteins involved in protoporphyrin biosynthesis (HEMA1, PORB), as well as proteins with putative roles in antioxidative stress responses (ENH1, CGLD27) accumulated differentially with similarly directed changes at the mRNA level. Also, NONINTRINSIC ABC PRO-TEIN1 (NAP1), a protein with homology to prokaryotic SufB proteins (66), is a member of this group. Fold-changes of the differentially expressed iron-responsive proteins discussed in the text and changes in the abundance of their cognate mRNAs are shown in Table I.
Post-Transcriptionally Regulated Proteins are Involved in Disparate Cellular Processes-A large subset of post-transcriptionally regulated proteins can be assigned to processes that themselves regulate post-transcriptional processes such as premRNA processing or translation. This positive feedback loop amplifies the extent of nontranscriptional regulation of gene expression and largely disconnects transcription from the proteomic readout. Several of the strongly up-regulated proteins harbor RNA-binding domains and may play putative regulatory roles. PUMILO22 (PUM22), a Puf-family protein that controls gene expression post-transcriptionally by pro-   (67,68), showed increased abundance in iron-deficient plants. Another member of this family, PUM4, was also up-regulated by iron-deficiency. Indicative of a regulatory role in iron homeostasis, the RRM/RBD/RNP motif containing protein At4g35785 was differentially phosphorylated upon iron deficiency in a phosphoproteomic survey (7), and was post-transcriptionally up-regulated in the present study. Increased abundance upon iron deficiency was also observed for the premRNA splicing factor RNA-DIRECTED DNA METHYLATION16 (RDM16). DICER-LIKE4 (DCL4), acting in processing tasiRNAs in a small RNA generating pathway, showed markedly decreased abundance in roots of iron-deficient plants.
Similar to what has been observed in roots, a large proportion of post-transcriptionally regulated, iron-responsive proteins in leaves is involved in post-transcriptional processes. The most obvious change in post-transcriptionally regulated proteins in leaves is the large number of ribosomal proteins (r-proteins) that showed increased abundance in response to iron deficiency. Furthermore, up-regulation of HISTONE DEACETYLATION COMPLEX1, two so-called "histone readers" that bind to histone H3 tails (the PHD finger containing proteins SHORT LIFE and ALFIN-LIKE3), and strong induction of the histone variants HTA6 and HTB2 indicates epigenetic changes induced by iron deficiency. Interestingly, the pool of up-regulated proteins also con- tains RAPTOR2, a key node in stress signal integration and cell growth (69,70). The majority of the post-transcriptionally down-regulated proteins is located in the chloroplasts. Two partners of NAP1, the plastidic SufD homolog NAP6 and the SufC homolog NAP7 that form a complex in chloroplasts involved in the transfer of Fe-S clusters to apoproteins (66, 71), were among the post-transcriptionally regulated proteins. Interestingly, NAP1 was transcriptionally regulated, whereas NAP6 and NAP7 were post-transcriptionally regulated, indicating that the assembly of functional modules could be regulated at different levels. The abundance of another protein involved in chloroplast Fe-S cluster assembly, the scaffold protein CHLOROPLAST-LOCALIZED ISCA-LIKE PROTEIN (CPISCA) (72), was also decreased upon iron deficiency.
Similar to the transcriptionally regulated CGLD27 (73, 74), LUTEIN DEFICIENT2 (LUT2) was predicted to function in carotenoid-xanthophyll metabolism related to the detoxification of reactive oxygen species. The down-regulated proteins also comprise HIGH CHLOROPHYLL FLUORESCENCE153 (HCF153), which is involved in the post-translational steps of cytochrome b6f complex biogenesis (75) and the NEET GROUP PROTEIN (AT-NEET) for which a role in iron metabolism and reactive oxygen species metabolism was proposed (76). Also down-regulated were the ribonuclease III family protein At1g55140, which was associated with RNA metabolism in the chloroplast (77), and the chlorophyllide a oxygenase CHLORINA1 (CH1), which is required for chlorophyll b synthesis. Consistent with the high iron content of PS1, many of the down-regulated proteins are associated with this photosystem. Also, PHOTOSYNTHETIC ELECTRON TRANSFER C, encoding the Rieske Fe-S center of cytochrome b6f complex (78), showed decreased abundance in leaves of irondeficient plants. The strongest post-transcriptional down-regulation was observed for FERREDOXIN C1.
RNA-Binding Proteins Occupy Central Hubs in PPI Networks of Post-Transcriptionally Regulated Proteins-Constructing protein-protein interaction (PPI) networks of posttranscriptionally regulated proteins resulted in large clusters dominated by r-proteins in both roots and shoots (Fig. 7). Although a high connectivity of r-proteins can be anticipated from their cooperativity in ribosomal function, smaller neighboring clusters and differentially expressed proteins associated with these central clusters are of particular interest.
Subcluster R2 contains, among other stress-responsive proteins, three DNAJ-like proteins, J2, J3, and J20, homologs of the cochaperone DNAJ protein from E. coli that were all up-regulated upon iron deficiency. J3 acts as an activator of the plasma membrane H ϩ -ATPase by interacting with and inhibiting the protein kinase PKS5 (81). We showed previously that pks5 mutants have increased proton pumping activity (12), supporting a scenario in which J3 activates AHA2 activity under iron-deficient conditions. Subcluster R3 contains several DNA-binding proteins (TAF6, ATR, RAD51C, CHR40, DTF2), which are involved in the regulation of transcription, chromatin remodeling and DNA repair, and proteins related to pathogen defense (At1g17600, CRK11, GTR1, BCS1). Several mitochondrial proteins associated with respiration (SDH2-2, At5g25450, At1g16700) and MITOCHONDRIAL EDITING FACTOR20 are part of this cluster. ABC TRANSPORTER OF THE MITO-CHONDRION1 (ATM1), a mitochondrial ATP-binding cassette transporter homologous to the yeast mitochondrial membrane protein ScATM1 involved in the mitochondrial export of Fe-S clusters or their precursors (82), was down-regulated upon Fe-deficiency but SDH2-2 and At5g25450 were upregulated, probably compensating for the decrease in mitochondrial respiration under iron-deficient conditions as a consequence of impaired heme and Fe-S cluster synthesis (83,84). The Fe-S cluster-binding protein At1g16700, a component of complex I, showed decreased abundance upon iron deficiency. A smaller cluster, R4, contains mainly plastid genes.
In leaves, the network is more complex. A larger cluster (L1) is dominated by r-proteins and also contains proteins related to photosynthesis, including, among others, a PSI subunit (PSAH2) and the chlorophyllide a oxygenase CH1 (Fig. 7B). Interestingly, a putative positive regulator of transcription (At5g42070) has a central position in this cluster. LUT2, several proteins involved in Fe-S cluster assembly (NAP6, NAP7, NFU2, SUFE2, HCF101, and CPISCA) and the Fe-S binding protein At4g36390, which was strongly down-regulated and highly interconnected with the Fe-S cluster assembly machinery, are members of this cluster.
Post-Transcriptional Regulation of Iron-Responsive Proteins is Partly Dependent on RGLG-Relatively large sets of proteins were iron-regulated solely in wild-type plants but not in rglg1 rglg2 mutants (275 proteins in roots and 233 in leaves) and were thus classified as RGLG-dependent. In addition, iron-responsive proteins that were deregulated in rglg1 rglg2 plants under iron-deficient conditions (i.e. regulated in the opposite direction) were included in this group. In roots of wild-type plants, 164 RGLG-dependent proteins were upregulated by iron deficiency, among them are many proteins with RNA-binding properties including r-proteins and proteins involved in RNA processing or metabolism. Several transporters for mineral nutrients localized on the plasma or vacuolar membrane were among the 114 down-regulated proteins, including NITRATE TRANSPORTER1 (NRT1), PHOSPHATE TRANSPORTER2 (PT2), MTPA2, and IRON-REGULATED PROTEIN2 (IREG2) (Fig. 8). The major plasma membrane transporters mediating the uptake of nitrate and phosphate in roots, NRT1 and PT2, were also down-regulated in iron-deficient wild-type plants. Uptake of anions that are acquired by proton-anion cotransport (nH ϩ /A Ϫ ) such as nitrate and phosphate (85,86) decreases the electrochemical potential of the plasma membrane and counteracts iron solubilization in the rhizosphere via ATPase-mediated proton extrusion (87). A more acidic uptake pattern (i.e. an increased cation/anion ratio), achieved by decreased anion uptake is thus advantageous for the plant under iron-deficient conditions. Remark-ably, the regulation of these transporters at the protein level was not observed in rglg1 rglg2 mutants.
An uncharacterized bicarbonate (HCO 3 Ϫ ) transporter, At3g06450, was highly up-regulated in iron-deficient plants at the protein level, but the transcript abundance remained unchanged (Fig. 8). No regulation of the protein abundance was observed in the mutant. Dark fixation of CO 2 via phosphoenolpyruvate carboxylase (PEPC), catalyzing the assimilation of bicarbonate to phosphoenolpyruvate to form oxaloacetate and inorganic phosphate, has been described as an adaptive response of roots to iron deficiency as a consequence of increased ATPase-mediated proton extrusion and a higher demand for organic acids (88,89). Similar to what was reported previously for Arabidopsis and other strategy I plants, PEPC1 was highly induced upon iron starvation both at the transcript and protein level. An iron-responsive bicarbonate transporter has not yet been described. It is thus tempting to speculate that At3g06450 fulfills such a function and has escaped detection in transcriptional surveys because of a lack of transcriptional regulation of the gene.
FIG. 8. RGLG-dependent regulation of transport proteins. Proton-cotransport of bicarbonate, nitrate and phosphate at the plasma membrane is post-transcriptionally down-regulated to maintain an acidic uptake pattern that facilitates the mobilization of Fe(III) oxides in the rhizosphere by ATPase-mediated proton extrusion. The transcript levels of these transporters remained unchanged. At the tonoplast, IREG2 and MTPA2 protein was decreased under iron-deficient conditions, whereas the genes were highly up-regulated at the transcriptional level. Numbers denote log 2 fold changes of the respective proteins. n.c. ϭ not changed.
In contrast to their cognate proteins that showed decreased abundance upon iron deficiency, IREG2 and MTPA2 transcripts were strongly enriched in roots of iron-deficient plants, indicating post-transcriptional regulation of these transporters. Reduced abundance of MTPA2 protein was not observed in rglg1 rglg2 mutants, indicating that post-transcriptional regulation of MTPA2 is dependent on functional RGLG protein. MTPA2, a vacuolar zinc transporter, is critical for the detoxification of surplus zinc in root cells of iron-deficient plants that accumulates intracellularly as a result of the low substrate specificity of the iron transporter IRT1 (46). We previously showed that the transcriptional regulation of MTPA2 is independent of the zinc supply and solely regulated by the iron status of the plants (81). Together, the present results that iron-dependent regulation of MTPA2 occurs at the transcriptional level, whereas protein abundance is post-transcriptionally controlled, presumably by the zinc status of the cell (Fig. 8).
The tonoplast transport protein IREG2, an ortholog of the mammalian iron efflux transporter ferroportin, was regulated in a similar manner, i.e. massive up-regulation at the transcriptional level and decreased protein abundance in response to iron deficiency. IREG2 sequesters nickel into the vacuole, but increased concentrations of nickel did not alter the expression of IREG2 (90), indicating that transcriptional regulation is independent of the actual concentration of potentially toxic metal ions. Co-expression analysis revealed closely coregulated, iron-dependent expression of IRT1, MTPA2, and IREG2 (81), suggesting that an "iron regulon" dictates transcriptional expression according to the iron status. All three genes are regulated by FIT (21), supporting this supposition. The pronounced up-regulation of MTPA2 and IREG2 mRNA may be interpreted in terms of "priming" the plant to a rapid production of the cognate proteins when the metal concentrations reach a critical level. Tuning of protein abundance and thus transporter activity to the prevailing metal concentration occurs post-transcriptionally and is dependent on RGLG (Fig. 8).
Does Iron Deficiency Induce Ribosome Specialization?-The present data show that iron deficiency post-transcriptionally induces changes in r-protein composition, leading to the formation of iron deficiency-specific ribosomes, a response that could not be deduced from a transcriptional analysis. Changes in r-protein expression were mainly post-transcriptionally controlled and thus masked in transcriptional surveys. Under iron-sufficient conditions, a large set of r-proteins showed significantly higher abundance in rglg1 rglg2 mutant plants compared with the wild type (five in roots and 16 in leaves, one r-protein down-regulated; Fig. 9). With one ex-ception, the higher abundance of r-proteins relative to the wild type was not accompanied by increased transcript abundance, indicating that RGLG is a negative, post-transcriptional regulator of r-protein expression.
In roots of wild-type plants, 15 r-proteins showed increased abundance upon iron deficiency, all of which were not accompanied by changes in mRNA level. In leaves, a subset of 54 r-proteins showed increased expression in iron-deficient relative to iron-sufficient plants. Similar to roots, all r-proteins were post-transcriptionally regulated. Several of the r-proteins that were differentially expressed in the wild type were either not iron-responsive in rglg1 rglg2 mutants (i.e. positively regulated by RGLG), or constitutively up-regulated, a group that we defined as negatively RGLG-regulated. With these criteria, we defined 50 r-proteins (eight in roots and 42 in leaves) as positively RGLG-regulated and nine r-proteins (two FIG. 10. Regulatory control of iron-responsive proteins. Processes directly related to iron homeostasis such as the acquisition, uptake, chelation, and storage of iron, as well as adaptive responses in the plastid and shoot-to-root communication are chiefly transcriptional regulated. Proteins involved in RNA-directed DNA methylation, RNA decay and changes in oxidative phosphorylation are primarily posttranscriptionally regulated. RGLG affects predominantly processes involved in protein translation and degradation, but also post-transcriptional modifications of histones. RGLG is also critical for tuning the abundance of transport proteins. Proteins listed in the figure are discussed in the text. in roots and seven in leaves) as negatively regulated by RGLG (Fig. 9).

DISCUSSION
Transcriptional Regulation of Iron Homeostasis: The Tip of the Iceberg?-In summary, here we provide a high-resolution, nearly complete inventory of the Arabidopsis "proteoferrome": proteins that accumulate differentially in roots and leaves between iron-sufficient and iron-deficient plants. A graphical overview of the mode of regulations of the proteins discussed here is given in Fig. 10. By uncovering an extensive and largely unexplored suite of iron-responsive proteins that are not or only partially transcriptionally regulated, we emphasize that a broad understanding of plant behavior in response to changing iron supply can only be achieved when post-transcriptional processes are taken into consideration. We show further that estimates made from changes in transcript levels may lead to incorrect predictions regarding the abundance of the cognate proteins. Concordance between mRNA and protein expression has often been used to estimate the contribution of transcriptional regulation on changes in protein abundance, and values around 40% have been inferred for several organisms from such calculations reviewed by Vogel and Marcotte (91). However, without a near-complete coverage of changes in the proteomic readout, such estimates rely on a relatively small proportion of transcripts and proteins for which the corresponding partners have been detected. In the present study, this fraction was estimated to account for 13 and 11% of the differentially expressed proteins in roots and leaves, respectively, resulting in a total contribution of transcription of 11 and 9% when the whole subset of differentially accumulated proteins and the concordance of mRNA and protein expression of this subset are taken into consideration. Thus, ϳ90% of the differentially expressed genes are not surveyed in transcriptional profiling studies, that several crit-ical regulatory mechanisms and key players are yet to be uncovered (Fig. 11).
Ribosome Composition -A novel Regulator of Translation?-In Arabidopsis, r-proteins are encoded by gene families comprising two to seven members, producing orthologs that for the most part are sufficiently unique to be distinguished by mass spectrometric analysis and that may have different properties. This stands in juxtaposition to mammals where r-proteins are encoded by a single gene and yeast where the diversity of r-proteins is less pronounced compared with plants. The heterogeneity of r-protein orthologs gives rise to speculations that specific, "made-to-order" ribosomes are generated that can capture mRNAs differentially and adapt the plant to the prevailing conditions by prioritizing a subset of mRNAs for translation (92)(93)(94)(95)(96). This hypothesis is particularly appealing from the perspective that translation in general and ribosome synthesis in particular are highly energy-consuming processes that need to be calibrated under stress conditions to secure plant vigor and productivity. Changes in r-gene expression have been described as being stress-specific (5,97), but the current analysis shows that alterations in r-protein expression are post-transcriptionally regulated, making a proteomic analysis mandatory to identify such changes. At the protein level, it was shown that growth condition can alter the relative incorporation of r-protein paralogs into ribosomes (98). Regulation of protein abundance at the level of translation was also supported by translatome profiling studies (99,100). These studies also support the assumption of a strong effect of "translational fitness" of specific mRNAs, a trait that is not necessarily associated with high steady-state mRNA levels (101,102). Several mechanisms at the translation level such as internal ribosome entry site (IRES)-mediated translation, ribosome pausing, and modification of translation initiation factors (103-105) may contribute to post-transcriptional control of gene activity, although their individual contribution FIG. 11. Estimate of the concordance between mRNA and protein expression in iron-deficient plants. Only 13% (roots) and 11% (leaves) of the proteins that are differentially expressed between iron-deficient and iron-sufficient plants (DE proteins) are associated with similarly regulated transcripts. For proteins in this subset, changes in transcript level account for circa 80% of the observed changes in protein abundance. in plants remains to be shown. The present data strongly support the idea of a regulatory intervention at the level of translation, aiding in translating a subset of mRNAs that are critical in adapting proteomic profiles to the prevailing conditions. Changes in ribosome composition, resulting in altered mRNA efficiency, may represent a key component in a sophisticatedly coordinated response to environmental cues.
RGLG Plays Critical Roles in The Orchestration of The Iron Deficiency Responses-Our results further show that the RING ligases RGLG1 and RGLG2 are required for the orchestration of the post-transcriptional responses, underlining the complexity of acclimation to abiotic stresses. This includes the whole repertoire of post-transcriptional processes comprising mRNA maturation, generation of small RNAs, mRNA export from the nucleus as well as changes at the epigenetic level such as DNA and histone modifications that are mediated by post-transcriptionally regulated proteins. RGLG can affect protein abundance directly by targeting proteins for proteasome-mediated degradation, indirectly via nonproteolytic reversible modifications such as monoubiquitilation and, in conjunction with UBC13, by promoting the formation of K63-linked polyubiquitin chains, which can affect diverse cellular processes such as transcription, histone function, membrane trafficking, DNA repair, and DNA replication (106 -109).
Interestingly, the percentage of proteins that were regulated without changes in corresponding mRNA abundance was lower for RGLG-regulated proteins than for the proteins that were regulated in common in both genotypes. In leaves, cognate transcripts for only 13 of the 240 RGLG-dependent proteins were identified. The concordance between mRNA and proteins was also markedly lower than that observed for the common set of genes (r ϭ 0.03). In roots, a similar picture was observed. Only 6% of the proteins (18 out of 278) were associated with a differentially expressed transcript, with a weak correlation between mRNAs and proteins (r ϭ 0.07). Together, these data indicate that RGLG is more engaged in post-transcriptional regulation of iron-responsive proteins and plays a minor and rather indirect role in the control of transcription.
Although the iceberg analogy would discredit the big steps forward already made by approaches targeting transcriptional regulation, the current investigation invites the exploration of a large portion of genes with altered expression upon iron starvation that may aid in generating iron-efficient germplasm and to combat iron deficiency-induced anemia.
Acknowledgements-We thank Drs. Louis Grillet and Isabel Cristina Vé lez-Bermú dez for critical reading of the manuscript. * Affymetrix GeneChip assays were performed by the Affymetrix Gene Expression Service Laboratory (http://ipmb.sinica.edu.tw/affy/) supported by Academia Sinica. Work in the Schmidt laboratory is supported by grants from the Ministry of Science and Technology (MoST) and Academia Sinica. rglg1 rglg2 mutants were kindly provided by Andreas Bachmair (University of Vienna). Proteomic mass spectrometry analyses and in-gel digestion was performed by the Proteomics Core Laboratory sponsored by the Institute of Plant and Microbial Biology and the Agricultural Biotechnology Research Center, Academia Sinica. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium [110] via the PRIDE partner repository with the dataset identifier PXD002126.
□ S This article contains supplemental Data Set S1 and Tables S1 and S2. § § Present address for I-Chun Pan: Department of Horticulture, National Chung Hsing University, Taichung, Taiwan.