Glutaredoxin GRXS17 Associates with the Cytosolic Iron-Sulfur Cluster Assembly Pathway   Open Access

Abstract

Cytosolic monothiol glutaredoxins (GRXs) are required in iron-sulfur (Fe-S) cluster delivery and iron sensing in yeast and mammals. In plants, it is unclear whether they have similar functions. Arabidopsis (Arabidopsis thaliana) has a sole class II cytosolic monothiol GRX encoded by GRXS17. Here, we used tandem affinity purification to establish that Arabidopsis GRXS17 associates with most known cytosolic Fe-S assembly (CIA) components. Similar to mutant plants with defective CIA components, grxs17 loss-of-function mutants showed some degree of hypersensitivity to DNA damage and elevated expression of DNA damage marker genes. We also found that several putative Fe-S client proteins directly bind to GRXS17, such as XANTHINE DEHYDROGENASE1 (XDH1), involved in the purine salvage pathway, and CYTOSOLIC THIOURIDYLASE SUBUNIT1 and CYTOSOLIC THIOURIDYLASE SUBUNIT2, both essential for the 2-thiolation step of 5-methoxycarbonylmethyl-2-thiouridine (mcm⁵s²U) modification of tRNAs. Correspondingly, profiling of the grxs17-1 mutant pointed to a perturbed flux through the purine degradation pathway and revealed that it phenocopied mutants in the elongator subunit ELO3, essential for the mcm⁵ tRNA modification step, although we did not find XDH1 activity or tRNA thiolation to be markedly reduced in the grxs17-1 mutant. Taken together, our data suggest that plant cytosolic monothiol GRXs associate with the CIA complex, as in other eukaryotes, and contribute to, but are not essential for, the correct functioning of client Fe-S proteins in unchallenged conditions.

A large fraction of plant intracellular iron is incorporated in iron-sulfur (Fe-S) prosthetic groups, which are well suited for electron transfer reactions and are often essential for the catalytic function of several enzymes (Balk and Schaedler, 2014). There are three pathways for the assembly of Fe-S clusters in plants: a mitochondrial, a plastidial, and a cytosolic pathway. The cytosolic Fe-S assembly (CIA) pathway is involved in the maturation of [2Fe-2S] clusters into [4Fe-4S] clusters and is capable of providing [4Fe-4S] clusters to cytosolic and nuclear proteins. Because of the absence of a cytosolic Cys desulfurase, the CIA pathway depends on the integrity of the mitochondrial Fe-S cluster machinery and a functional ABC transporter of the mitochondria (Kispal et al., 1999; Gerber et al., 2004). The three pathways share a common mechanism in which the Fe-S is preassembled on a scaffold protein and then transferred to apoproteins by carriers or targeting factors (Fig. 1A; Couturier et al., 2013; Balk and Schaedler, 2014).

Glutaredoxins (GRXs) together with thioredoxins (TRXs) are thiol oxidoreductases that are able to control the redox state of proteins and are present in most organisms (Herrero and de la Torre-Ruiz, 2007). The yeast GRX proteins Grx3/4 and the mammalian ortholog GRX3/PKC-interacting cousin of TRX (PICOT) have been associated with the CIA pathway and contain themselves [2Fe-2S] clusters (Picciocchi et al., 2007; Haunhorst et al., 2010). Deletion of Grx3/4 in yeast leads to defects in cytosolic and mitochondrial Fe-S assembly, deregulation of iron homeostasis, and defects in proteins containing di-iron centers (Mühlenhoff et al., 2010). Yeast Grx3/4 and human GRX3 belong to the PICOT protein family and contain one N-terminal TRX and one (Grx3/4) or two (GRX3) C-terminal GRX domains, also known as PICOT homology domains (Haunhorst et al., 2010). Because they contain only a single Cys residue in their GRX active sites, they are classified as monothiol GRXs. They are conserved and present in a broad range of organisms, including bacteria, yeasts, plants, and mammals (Isakov et al., 2000). Whereas there are other monothiol GRXs present in mitochondria, Grx3/4 and GRX3 are the only nucleocytosolic-localized monothiol GRXs (Herrero and de la Torre-Ruiz, 2007).

The sole class II Arabidopsis (Arabidopsis thaliana) nucleocytosolic monothiol GRX is GRXS17, which contains one N-terminal TRX and three C-terminal GRX domains. GRXS17 dimers are capable of associating with three [2Fe-2S] clusters in vitro (Knuesting et al., 2015). Its physiological and molecular role in plants is not well understood (Couturier et al., 2013). GRXS17 function has been associated with protection against oxidative stress in Arabidopsis and thermotolerance in Arabidopsis and tomato (Cheng et al., 2011; Wu et al., 2012; Knuesting et al., 2015). Arabidopsis GRXS17 loss-of-function plants (grxs17) are hypersensitive to heat stress and show alterations in auxin sensitivity and polar transport (Cheng et al., 2011). The molecular function of an association of cytosolic monothiol GRX with Fe-S clusters and the CIA pathway has been a subject of debate, and a role in Fe, Fe-S, or oxidative signaling has been proposed, in addition to a role in delivery or repair of Fe-S clusters (Couturier et al., 2015). Recently, delivery of an Fe-S cluster by human GRX3 to the CIA pathway component DRE2/Anamorsin has been demonstrated (Banci et al., 2015).

Cytosolic tRNAs in eukaryotes carry several chemical modifications, often at the anticodon loop. Uridines at the first position of the anticodon (U₃₄) of tRNAs tK(UUU), tE(UUC), and tQ(UUG) are modified to 5-methoxycarbonylmethyl-2-thiouridine (mcm⁵s²U) in eukaryotes. Furthermore, this evolutionarily conserved modification (Mehlgarten et al., 2010) is essential for unperturbed translation and cellular signaling (Zinshteyn and Gilbert, 2013; Scheidt et al., 2014; Nedialkova and Leidel, 2015). The 2-thiolation (s²) step of mcm⁵s²U is catalyzed by the UBIQUITIN-RELATED MODIFIER (URM1) pathway and requires the CIA complex in yeast (Nakai et al., 2007; Leidel et al., 2009). The mcm⁵ modification is catalyzed by the elongator (ELP) pathway and requires the Elp3/ELO3 catalytic subunit, i.e. a [4Fe-4S] protein (Huang et al., 2005; Paraskevopoulou et al., 2006; Selvadurai et al., 2014), together with the Trm9/Trm112 complex, which is a tRNA methyltransferase necessary for the last step of mcm⁵ formation (Kalhor and Clarke, 2003; Chen et al., 2011; Leihne et al., 2011). The presence of the mcm⁵ chain is necessary for an effective thiolation of tRNAs (Leidel et al., 2009; Chen et al., 2011). Although in vivo data are scarce, tRNA modifications might have a regulatory function, because certain open reading frames (ORFs) are enriched in codons recognized by modified tRNAs. In yeast, genes involved in the DNA-damage response are enriched in GAA, AAA, and CAA codons and elongator mutants defective in mcm⁵s² modification are hypersensitive to DNA stress (Chen et al., 2011).

Here, we found the Arabidopsis glutaredoxin GRXS17 to associate with members of the CIA pathway machinery and with several putative Fe-S client proteins, such as XANTHINE DEHYDROGENASE1 (XDH1) and the URM1 pathway components CYTOSOLIC THIOURIDYLASE SUBUNIT1 (CTU1) and CTU2. Accordingly, mutants lacking GRXS17 showed similar phenotypes, at molecular, cellular, and/or physiological levels, as mutants in genes encoding CIA components or proteins involved in mcm⁵s² tRNA modification. Our results endorse the association of cytosolic monothiol GRXs with the CIA complex, Fe-S cluster metabolism, and tRNA modification in plants.

RESULTS

GRXS17 Associates with the CIA Complex and Fe-S Proteins

To characterize the role of Arabidopsis GRXS17 and uncover potential links with Fe-S cluster assembly or delivery to proteins, we used tandem affinity purification (TAP; Van Leene et al., 2015). GRXS17 was fused both to an N-terminal and a C-terminal TAP-tag and expressed both in Arabidopsis cell cultures and seedlings under the control of the cauliflower mosaic virus (CaMV) 35S promoter. Copurified proteins included nearly all known core components of the CIA pathway, two known Fe-S proteins, i.e. XDH1 and BolA2, and a number of proteins implicated in tRNA metabolism (Table I; Supplemental Dataset S1). Notably, we identified more GRXS17-interacting proteins when using the C-terminal TAP-tag than when using the N-terminal TAP-tag, suggesting that the position of the tag affects protein folding or that different domains of GRXS17 interact with proteins related to different functions. It has recently been suggested that the human GRXS17 ortholog GRX3 establishes protein-protein interactions with its N-terminal TRX domain and that this interaction favors the correct orientation of the Fe-S cluster binding domains (GRX donors) to transfer the cluster to an acceptor protein (Banci et al., 2015).

GRXS17 associates with the CIA complex and Fe-S proteins in Arabidopsis cells and seedlings

The CIA pathway is responsible for providing Fe-S clusters to respective apoproteins in the cytosol and the nucleus (Bernard et al., 2013). In the current model, electrons provided by NADPH oxidation are transferred by Top1T722A MUTANT HYPERSENSITIVE and DEREPRESSED FOR RIBOSOMAL PROTEIN S14 EXPRESSION (DRE2) to the scaffold protein NUCLEOTIDE BINDING PROTEIN35. Once the [4Fe-4S] clusters have been assembled on the NUCLEOTIDE BINDING PROTEIN35 scaffold, they are transferred to target apoproteins by dedicated proteins forming the CIA targeting complex (Fig. 1A). In Arabidopsis, this complex is composed of NUCLEAR ARCHITECTURE RELATED1 (NAR1), CIA1, AS1/2 ENHANCER7 (AE7/CIA2), and METHIONINE REQUIRING18/METHYL METHANESULFONATE19 (MET18/MMS19) and locates to both the cytosol and the nucleus (Luo et al., 2012). Grx3, the yeast GRXS17 ortholog, has been demonstrated to interact with DRE2 in yeast (Zhang et al., 2011), and recently, it has been shown that the human GRX3 is able to transfer [2Fe-2S] clusters to anamorsin, the human DRE2 ortholog (Banci et al., 2015). We thus hypothesized that GRXS17 might similarly directly interact with DRE2 in Arabidopsis, which was confirmed by yeast two-hybrid analysis (Y2H; Fig. 1B). We confirmed this interaction in planta by bimolecular fluorescent complementation (BiFC). Proteins were fused with N-terminal or C-terminal fragments of GFP (designated nGFP and cGFP, respectively) and coexpressed transiently in leaves of Nicotiana benthamiana (Supplemental Fig. S1A). A cytosolic signal was obtained with the combination of nGFP-GRXS17 and DRE2-cGFP, which corresponded to the localization of GFP-DRE2 (Supplemental Fig. S1C). Thus, we concluded that Arabidopsis GRXS17 associates with the CIA complex as in yeast and humans.

Loss-of-Function grxs17 Plants Show Some Degree of Hypersensitivity to DNA Damage

To obtain insights into the function of GRXS17, we profiled the transcriptome of grxs17 loss-of-function seedlings. A transfer DNA (T-DNA) insertion line (grxs17-1; SALK_021301) described previously (Supplemental Fig. S1; Cheng et al., 2011; Knuesting et al., 2015) was characterized in detail and found to lack both GRXS17 transcripts (Supplemental Fig. S2B) and GRXS17 protein (Supplemental Fig. S2C). Fourteen-day-old seedlings were grown in vitro and used for RNA sequencing (RNA-Seq). One hundred eighty-eight genes were found to be significantly up-regulated, and 209 were down-regulated (Supplemental Table S1). Using ATTED-II (http://atted.jp/), we identified a conserved coexpression network comprising 25 genes up-regulated in grxs17-1 mutants (Fig. 1C; Supplemental Table S1). Many of the genes in this network are known to be induced upon genotoxic stress (Ascencio-Ibáñez et al., 2008). Because Arabidopsis mutant plants compromised in components of the purified CIA complex (Table I), namely AE7/CIA2 (Luo et al., 2012) and MET18/MMS19 (Han et al., 2015), have been reported to have increased expression of DNA-damage response genes (Luo et al., 2012; Han et al., 2015), we analyzed the expression of POLY(ADP-RIBOSE) POLYMERASE2 (PARP2), a frequently used marker gene, in the grxs17-1 mutant, as well as in an antisense line, denominated grxs17-2, that was previously described (Cheng et al., 2011) and which we confirmed to lack GRXS17 protein accumulation (Supplemental Fig. S2C). PARP2 transcript levels were indeed induced in both mutant alleles (Fig. 1D). Therefore, we assessed the sensitivity of the grxs17-1 and grsx17-2 mutant lines to DNA-damage stress. In agreement with the transcriptome data, the grxs17 mutants were found to display some degree of hypersensitivity to the DNA-damage agent MMS (Fig. 1E), though not as marked as the met18 mutants (Fig. 1E; Luo et al., 2012; Han et al., 2015). Together, these data suggest that GRXS17 not only forms part of the CIA complex in Arabidopsis but also may contribute to its function in the DNA-damage response.

Metabolic Profiling of Loss-of-Function grxs17 Plants Points to a Perturbed Flux through the Purine Salvage Pathway

Although human and yeast GRXS17 orthologs are essential for the maturation of cytosolic and nuclear Fe-S proteins (Mühlenhoff et al., 2010; Haunhorst et al., 2013), Arabidopsis GRXS17 has recently been shown to play only a minor role in maintaining the activity of two classes of cytosolic Fe-S enzymes (aconitases and aldehyde oxidases; Knuesting et al., 2015). Nonetheless, we identified another Fe-S enzyme, XDH1, to be associated with GRXS17 (Table I). XDH1 belongs to the family of xanthine oxidoreductases (XORs) that catalyze the oxidation of hypoxanthine and xanthine to uric acid during purine degradation. XDH1 is active as a homodimer in which each monomer has four cofactors: two [2Fe-2S] prosthetic groups, one FAD molecule, and one molybdenum cofactor (Moco), bound to the N-terminal, central, and C-terminal part of the protein, respectively (Zarepour et al., 2010).

Because XDH1 contains two [2Fe-2S] cofactors that are required for electron transfer to reduce the substrates xanthine and hypoxanthine (Zarepour et al., 2010), we investigated the GRXS17-XDH1 interaction in more detail. Interaction between XDH1 and GRXS17 was confirmed in a Y2H assay (Fig. 2A) and by BiFC, showing a cytosolic localization (Supplemental Fig. S1, B and C). To examine the physiological relevance of this interaction, we measured XDH activity in grx17-1 and grxs17-2 seedlings (Fig. 2B) but found no statistically significant difference in XDH activity compared with wild-type seedlings (Fig. 2B). Hence, similar to what was previously found for the enzymes cytosolic aconitase and aldehyde oxidase (Knuesting et al., 2015), our findings indicate that GRXS17 is not essential for XDH1 enzyme activity in plants.

In addition to XDH1, three additional proteins required for purine degradation were identified in complex with GRXS17: URIDINE RIBOHYDROLASE1 and 2 (URH1/2) and the urease accessory protein UreG (Table I). URH1 and URH2 are nucleoside ribohydrolases with inosine and xanthosine hydrolytic activity, leading to hypoxanthine and xanthine generation (Riegler et al., 2011). Urease catalyzes the hydrolysis of urea to ammonia and carbon dioxide, which are the final products of purine degradation. In Arabidopsis, urease activation is mediated by the accessory proteins UreD, UreF, and UreG (Witte et al., 2005). To evaluate the impact of loss-of-function of GRXS17 on purine metabolism, we measured the accumulation of uric acid precursors (hypoxanthine, xanthine, and xanthosine) and of ureides (uric acid, allantoin, allantoic acid, and urea) in wild-type and grxs17 seedlings (Fig. 2C; Supplemental Fig. S3). Although the effects were modest, our results are in agreement with the function of XDH1 in the purine salvage pathway and point to a perturbed flux through this pathway in the GRXS17 loss-of-function plants with enhanced accumulation of hypoxanthine, xanthine, and xanthosine and reduced accumulation of allantoic acid and urea. Together, these results suggest that association of GRXS17 with purine metabolism enzymes may be functionally relevant and may contribute to the purine salvage pathway.

GRXS17 Interacts with tRNA Thiolation Proteins

The TAP of GRXS17 protein complexes revealed interactions with several proteins that are linked to tRNA metabolism (Table I). For instance, we found GRXS17 to associate with the Arabidopsis orthologs of CTU1 and CTU2 required for 2-thiolation of U₃₄ of the three cytosolic tRNA species tK(UUU), tE(UUC), and tQ(UUG). These tRNAs carry an mcm⁵s²U₃₄ modification, which is conserved throughout archae and eukaryotes (Fig. 3A; Leidel et al., 2009; Philipp et al., 2014). Although only the protein encoded by AT2G44270 was described as the Arabidopsis CTU1 ortholog (Leiber et al., 2010), several peptides copurified with GRXS17 could not be distinguished from the product of locus AT1G76170 with 93% sequence similarity to CTU1 (Table I) and was named CTU1-like (CTU1L). To verify the interaction between GRXS17 and CTU1 and CTU2, we performed Y2H and BiFC assays. Y2H indicated only a very weak interaction for both CTU1 and CTU2 (Fig. 3B). Through BiFC, we could confirm a cytosolic-localized GRXS17-CTU2 interaction (Fig. 3C; Supplemental Fig. S4), thus supporting the TAP results and confirming that GRXS17 interacts directly with CTU2. Interaction between GRXS17 and CTU1 could not be detected in BiFC, which may either be a false negative result or an indication of indirect interaction, for instance via heterodimerization with CTU2.

To further investigate the link between GRXS17 and CTUs, we tested if GRXS17 was required for tRNA thiolation in Arabidopsis. Thiolated tRNAs were visualized using their retardation in PAGE in the presence of [(N-acryloylamino)phenyl]mercuric chloride (APM), a compound that interacts with 2-thiouridine (Igloi, 1988). As reported previously, ctu1 and ctu2 T-DNA insertion lines did not contain thiolated tRNA (Fig. 3C; Leiber et al., 2010; Philipp et al., 2014). In contrast, in grxs17-1 plants, levels of thiolated total tRNA were similar to wild type (Fig. 3C). Thus, GRXS17 interacts with CTU subunits but is not essential for their tRNA anticodon thiolation function in Arabidopsis.

GRXS17 and Elongator Functions Are Related in Arabidopsis

Whereas CTU1 and CTU2 are essential for 2-thiolation of mcm⁵-modified tRNA anticodons, the ncm⁵U₃₄ modification itself is dependent on elongator (Fig. 3A; Huang et al., 2005; Selvadurai et al., 2014), a complex that is structurally and functionally conserved between Saccharomyces cerevisiae and Arabidopsis (Mehlgarten et al., 2010). ELO3, the Arabidopsis ortholog of yeast elongator subunit Elp3, contains histone acetyltransferase and radical sterile alpha motif domains, which were recently shown to be catalytically critical for the tRNA modification function of archaeal Elp3 (Selvadurai et al., 2014).

grxs17-1 and grxs17-2 mutant plants exhibit an elongated leaf phenotype with larger leaf length/width ratio as compared to wild type (Supplemental Fig. S5, A and B), which is typical of elongator mutants (Nelissen et al., 2005). Therefore, we investigated the possibility of a link between GRXS17 and elongator. First, we studied the root and leaf phenotypes of grxs17-1, grxs17-2, and elo3-6 in the Columbia-0 (Col-0) background. The elo3-6 and the two grxs17 mutant alleles had a reduced primary root growth (Fig. 4A), which was correlated with a reduced number of cortex cells in the root apical meristem (Fig. 4B), suggesting a faster transition to differentiation. Second, we measured the total leaf area in elo3-6 and the two grxs17 mutant alleles relative to wild-type plants. Notably, all three mutants had a similar growth profile: (1) larger juvenile leaves 1, 2, and 3, which are fully expanded at 24 d after germination (DAG); (2) a larger leaf 4, which is the transition to adult stage; and (3) smaller adult leaves 6, 7, and 8 because of a delay in growth in the mutants (Supplemental Fig. S5A). This is also reflected in the absence of leaf 8 and cauline leaves at the 24 DAG time point in some of the mutant genotypes (Fig. 4C). We examined the cellular basis of the changes in the fully developed leaf 3, of which the total cell number and final cell area are representative of cell proliferation and growth activities during its development. We observed that the cell area in all mutant lines was similar to that of the Col-0 control. However, the calculated number of cells per leaf was significantly higher in the two grxs17 mutant genotypes and in the elo3-6 line as compared to wild type (Supplemental Fig. S5, C and D), indicating that the larger leaf 3 area is the result of more cell proliferation in grxs17 and elo3-6 mutant lines.

Next, we compared our RNA-Seq dataset (Supplemental Table S1) with previously published microarray datasets of the elongator mutants elo2-1 and elo3-1 (Nelissen et al., 2005). Although the latter dataset was performed on mutants in the Landsberg erecta ecotype and with a different experimental setup, we observed a large overlap between the top 100 genes induced in grxs17-1, elo2-1, and elo3-1 (Fig. 4D). The overlap comprised the DNA-damage network indicated above (Fig. 1C). Indeed, loss of ELO3 has been reported to lead to a DNA-damage response (Xu et al., 2012). We assessed this observation further by analyzing gene expression in the grxs17 and elo3-6 mutants and in the T-DNA insertion lines of CTU1 and CTU2, called ctu1 and ctu2, all in the Col-0 background (Fig. 4E; Supplemental Fig. S6). Indeed, all genes tested that were up-regulated in grxs17 mutants were also up-regulated in elo3-6. Similarly, also in ctu1 and ctu2 we observed an upregulation of many grxs17-up-regulated genes, although more modest. In agreement with these results, ctu1 and ctu2 also showed some degree of sensitivity to the DNA-damage agent MMS (Fig. 4F).

Finally, previous reports indicated that the elongator subunits ELP2 and ELP3/ELO3 are involved in the salicylic acid signaling pathway (DeFraia et al., 2010, 2013) and that the elongator complex is required for Arabidopsis resistance to necrotrophic fungal pathogens such as Botrytis cinerea and Alternaria brassicicola (Wang et al., 2015). Therefore, we tested the susceptibility of the grxs17-1 mutant to B. cinerea. Compared with wild-type plants, grxs17-1 plants were more susceptible to B. cinerea infection (Fig. 4G), suggesting that GRXS17 plays a role in the defense response against B. cinerea.

Taken together, our data indicate that grxs17 and elo3 mutants have similar phenotypes at the molecular, cellular, and physiological level, both in growth and in defense, and thus support a joined, but not necessarily identical, role in the functioning of particular processes such as tRNA modification.

GRX and Elongator Functions Are Related, Also in Yeast

Given the conservation of the GRX proteins in eukaryotes, we wanted to test whether the yeast orthologs of Arabidopsis GRXS17, i.e. Grx3/4 (Herrero and de la Torre-Ruiz, 2007), also have similar related functions in tRNA modification by elongator in S. cerevisiae. Because grx3 grx4 double mutants are barely viable in S288c (Wu et al., 2012) and lethal in the background of yeast strain W303-1A, we used a previously described GRX3 deletion strain of W303 that maintains the GRX4 gene under the control of a GAL1 promoter, which is repressible by Glc and inducible by Gal (Mühlenhoff et al., 2010). W303 contains an ochre stop codon (UAA) mutation in the ADE2 locus (ade2-1ochre), which leads to premature translation termination and causes adenine auxotrophy (Fig. 5A). This phenotype can be suppressed by read-through of SUP4, an ochre suppressor tRNA^Tyr carrying a U₃₄ mutation in its anticodon (Goodman et al., 1977; Huang et al., 2005; Jablonowski et al., 2006; Fig. 5A). Because the SUP4 tRNA^Tyr needs an mcm⁵-modified U₃₄ to decode the UAA triplet as Tyr (Huang et al., 2005; Jablonowski et al., 2006), it cannot suppress the ade2-1ochre reporter without proper mcm⁵ formation at U₃₄. SUP4 failed to promote stop codon read-through of ade2-1  ^ochre and hence, conferred adenine auxotrophy in the yeast reporter strain lacking GRX3 function, even if GRX4 was not repressed (Fig. 5A). This indicates that GRX3 gene function in yeast is critical for SUP4 function and thus supports a related role with elongator, as in Arabidopsis.

When we analyzed tRNA thiolation in elo3-6 mutant, we observed, in contrast to grxs17-1, almost no thiolated tRNA (Fig. 3D). These results suggest that grxs17-1 presents no strong defect in mcm⁵ modification in plants based on the absence of a strong thiolation defect. However, in yeast and plants, the GRXs Grx3/4 and GRXS17 may have a potential influence on elongator tRNA modification function, which we have previously shown to be conserved in and functionally exchangeable between S. cerevisiae and Arabidopsis (Mehlgarten et al., 2010). In the light of these findings, it is interesting to note that the GRX3 gene was identified previously as a high-copy suppressor in yeast of growth inhibition by zymocin, a fungal tRNase ribotoxin complex (Jablonowski et al., 2001; Jablonowski and Schaffrath, 2007). Zymocin activity targets elongator-dependent mcm⁵s²U34 modifications in tRNA anticodons (Lu et al., 2005; Jablonowski et al., 2006) and eventually kills S. cerevisiae cells. Therefore, loss of tRNA modification in elongator (elp) or tRNA methyltransferase (trm9) mutants protects against zymocin, making the tRNase a useful tool for diagnosing elongator function and, hence, tRNA modification in yeast (Nandakumar et al., 2008) and potentially in Arabidopsis (Mehlgarten et al., 2010; Leitner et al., 2015). We tested zymocin resistance in a yeast ELP3 deletion strain and in a GRX3 deletion strain having the GRX4 gene under control of the GAL1 promoter. As expected, deletion of ELP3 resulted in resistance to zymocin due to loss of mcm⁵s²U₃₄ modification (Fig. 5B). However, when GRX3  /4 were deleted, yeast cells could not grow in the presence of the toxin. Although these results suggest appropriate tRNA modifications_, our SUP4 assay above (Fig. 5A) together with previous data showing that high-copy GRX3 confers zymocin resistance (Jablonowski et al., 2001), may indicate that the absence of GRX3  /4 function differentially affects mcm⁵ or mcm⁵s² modifications at U₃₄, an option that has to await further elucidation.

DISCUSSION

Cytosolic monothiol GRXs are required in Fe-S cluster delivery in yeast and mammals. In plants, it remained unclear whether they have similar functions. Here, we used interactome analysis and in-depth functional characterization of a loss-of-function mutant to establish that that the sole class II cytosolic monothiol GRX from Arabidopsis associates with CIA components and Fe-S client proteins and contributes to their functioning.

Profiling of Loss-of-Function grxs17 Plants Points to an Elevated DNA-Damage Response

The function of many CIA components in yeast, mammals, and Arabidopsis has been associated with genomic stability, DNA repair, and metabolism, because many proteins necessary for DNA replication and repair are known to contain Fe-S clusters. In yeast and humans, Met-18/MMS19 has been shown to associate not only with several other CIA components, but also with Fe-S target proteins involved in DNA metabolism and to mediate the maturation of certain Fe-S proteins involved in DNA repair and replication (DNA helicases, polymerases, nucleases or glycosylases; Gari et al., 2012; Stehling et al., 2012; van Wietmarschen et al., 2012). The human MMS19 has been shown to be necessary for the maturation of only certain Fe-S proteins, mostly involved in DNA metabolism, but not for the activity of cytosolic aconitase iron regulated protein1 (IRP1) or Gln phosphoribosylpyrophosphate amidotransferase (Stehling et al., 2012). Human GRX3 is able to bind, in addition to [2Fe-2S] clusters, [4Fe-4S] clusters in vitro, which is necessary for the maturation of apo-IRP1 into aconitase (Xia et al., 2015), thus suggesting that GRX3/GRXS17, and not MMS19/MET18, is involved in the transfer of [4Fe-4S] clusters necessary for IRP1 maturation. This hypothesis is in accordance with the decrease in cytosolic aconitase activity observed in the Arabidopsis grxs17-1 mutant (Knuesting et al., 2015).

In Arabidopsis, MET18 also interacts with DNA polymerases and facilitates the maturation of Fe-S clusters on DNA polymerases in the cytosol (Han et al., 2015). Another MET18-interacting protein is the DNA demethylase ROS1, an Fe-S client protein that requires the Fe-S cluster for its activity, linking the CIA pathway with epigenetic modifications (Duan et al., 2015). Accordingly, an elevated DNA-damage response was found in Arabidopsis mutants deficient in CIA complex components, such as ae-7 and met18 (Luo et al., 2012; Han et al., 2015), or in Fe-S cluster containing proteins, such as elo3 (Xu et al., 2012).

In our study, we discovered a network (or regulon) that comprises genes involved in the genotoxic (DNA-damage) stress response to be up-regulated in the grxs17-1 mutant. Accordingly, we found that grxs17 plants present some degree of hypersensitivity to the DNA-alkylating agent MMS. A similar response has previously been reported for the ae-7, met18, and elo3-6 mutants (Luo et al., 2012; Xu et al., 2012; Han et al., 2015). Given that most enzymes implicated in DNA replication and DNA repair contain Fe-S clusters, our data suggest that GRXS17 and its association with the CIA complex may be necessary for their correct functioning and thus in genome integrity maintenance in plants.

GRXS17 Associates with Fe-S Client Proteins

Among the GRXS17-associated proteins that we identified by TAP were CTU1 and CTU2, two proteins essential for the thiolation of uridine at the wobble position of cytosolic tRNA in eukaryotes (Björk et al., 2007; Schlieker et al., 2008; Leidel et al., 2009). CTU1 presents homology with Escherichia coli TtcA, the protein responsible for the thioltransferase activity necessary for s²C₃₂ tRNA thiolation, a tRNA modification not present in eukaryotes (Jäger et al., 2004). E. coli TtcA was shown to bind, through three conserved Cys residues, Fe-S clusters that are essential for its activity (Bouvier et al., 2014). The conserved motifs Cys-X₁-X₂-Cys present in E. coli TtcA and Arabidopsis CTU1 suggest that Arabidopsis CTU1(L) proteins could be Fe-S client proteins. Accordingly, proteins from the CIA machinery are known to be required for 2-thio modification of cytosolic tRNAs, suggesting that at least one cytosolic or nuclear protein containing Fe-S clusters is necessary for thiolation of cytosolic tRNAs (Nakai et al., 2007). Accordingly, GRXS17 could be involved in the transfer of putative Fe-S clusters to a CTU1(L)/CTU2 complex.

Several proteins already known to bind Fe-S clusters were also found in our GRXS17 TAP interactome, including BolA2 and XDH1. Interactions between GRX and BolA proteins are conserved in yeast, humans, and plants. In all of these eukaryotes, it has been demonstrated that the GRX and BolA domains are bridged by the binding of a [2Fe-2S] cluster (Li et al., 2009, 2012; Couturier et al., 2014). However, GRXS17 can also bind Fe-S clusters independently of BolA2 interaction, through the formation of Fe-S bridged homodimers, and it can contribute to the activity of cytosolic Fe-S enzymes (Knuesting et al., 2015).

In yeast, Grx3/4 are localized in the nucleus, where they regulate the nuclear export of Aft1, a transcription factor that regulates Fe-responsive gene expression under low Fe conditions (Pujol-Carrion et al., 2006). Heterodimerization of Grx3/4 with the BolA-like protein Fra2 leads to nuclear export of Aft1 (Pujol-Carrion et al., 2006; Li et al., 2009). In both yeast and humans, GRXs play an important role in Fe utilization. Grx3/4-deficient yeast cells have defects in de novo synthesis of Fe-S clusters and heme, two of the major Fe-consuming processes in these organisms. The Fe content is not decreased in these yeast cells, but the Fe is not bioavailable due to deficient Fe delivery to mitochondria. A similar defect in Fe bioavailability was observed in zebrafish (Danio rerio) and HeLa cells when HsGRX3 (or the zebrafish ortholog) was depleted (Mühlenhoff et al., 2010; Haunhorst et al., 2013). In plants, BolA2 does not seem to play a role in Fe homeostasis, and the biological output of the GRXS17-BolA2 interaction remains an unanswered question (Couturier et al., 2014; Roret et al., 2014).

XDH1, in turn, belongs to the family of XORs and is a central player in purine catabolism. In Arabidopsis, two XOR-encoding genes are present with a strict XDH activity, i.e. XDH1 and XDH2. Although both genes are located close to each other on the same chromosome, XDH2 is expressed constitutively at basal levels, whereas XDH1 is differentially expressed upon several developmental and environmental stimuli (Hesberg et al., 2004). Plants deficient in XDH1 are affected in growth and development. Seedlings with RNAi-mediated silencing of XDH1 are smaller than wild-type plants with shorter flowering stems, smaller fruit size, and higher sterility rate (Nakagawa et al., 2007). In addition, XDH1-deficient plants also show earlier onset of age-dependent, dark-induced leaf senescence (Nakagawa et al., 2007; Brychkova et al., 2008) and decreased tolerance to drought stress (Watanabe et al., 2010). In accordance with the function of XDH1 in purine catabolism, precursors of uric acid (hypoxanthine and xanthine) are significantly more abundant in XDH1-deficient plants, whereas downstream products (allantoic acid and urea) are less abundant (Nakagawa et al., 2007; Brychkova et al., 2008). Quantification of these metabolites in grxs17-1 plants indicated that these purine catabolism intermediates also accumulate differentially in the absence of GRXS17, reflecting a perturbed flux through the purine salvage pathway. However, it is interesting to note that despite the apparent clarity of these results, certain aspects of the changes in the levels of these metabolites were unexpected, namely the 20% increase in uric acid in the mutant in spite of the reduction of the activity of XDH1. The most likely explanation for this is the upregulation of another, as-yet-undefined, route of uric acid production in an unsuccessful (when viewed from the standpoint of the other downstream metabolites) attempt to compensate for the deficit of XDH1. This aspect suggests additional complexity in the regulation of this pathway and thus warrants further experimental study.

Remarkably, no overlap was found between the interacting partners identified in our study and those in a recent study in which affinity chromatography with HIS-tagged GRXS17 was used to identify GRXS17-interacting proteins (Knuesting et al., 2015). Interactors reported in the latter study were more related to signaling and included the transcriptional regulator Nuclear Factor Y Subunit C11/Negative Cofactor 2α (Knuesting et al., 2015). The discrepancy between these two interactomes could be due to differences in technical procedures and suggests that GRXS17 is able to interact with even more proteins than reported here and that it may be involved in different cellular processes.

GRXS17 Function Is Evolutionarily Conserved

The evolutionary conservation of the processes described in this paper is strong: CIA complex, GRX, XOR, and tRNA-modifying proteins are all conserved in plants and humans. Fe-S proteins are thought to be reminiscent of the origin of life, when Fe and S were readily available under a reducing, anaerobic environment. Metabolic pathways that evolved at these early stages of life became essential to all living organisms, and many of them require Fe-S proteins (Sheftel et al., 2010). Because Arabidopsis grxs17 seedlings actually only show mild developmental defects (Cheng et al., 2011; Knuesting et al., 2015) in contrast to the strong developmental defects associated with GRX depletion in yeast or animal cells, plants such as Arabidopsis may have evolved alternative protein(s) or pathway(s) that can take over these evolutionary conserved tasks.

The role of GRX proteins in the CIA pathway is still not entirely understood. In yeast, incorporation of Fe-S clusters into both cytosolic and mitochondrial proteins is dependent on Grx3/4 (Mühlenhoff et al., 2010). We establish here that, in plants, GRXS17 associates with several components of the CIA complex, including DRE2, NAR1, CIA1, AE7/CIA2, and MET18. GRXS17-DRE2 interaction is conserved in yeast and humans, in which the orthologs also interact directly and in which this close interaction allows the transfer of the [2Fe-2S] clusters from the GRX to DRE2 (Banci et al., 2015; Fig. 1A). In humans, disturbance of the DRE2/Anamorsin-GRX interaction has been proposed as a strategy to reduce cell proliferation in solid tumors in which Anamorsin expression is enhanced (Saito et al., 2011). Possibly, because maturation of DRE2 in humans is dependent on GRX3, GRX3 might, analogously to GRXS17, form the bridge between the mitochondrial and cytosolic Fe-S assembly machineries by providing the initial Fe-S cluster to the essential CIA proteins. Alternatively, GRXS17 might specifically contribute to the maturation of [2Fe-2S] cluster-containing proteins in the cytosol and the nucleus, because the clusters assembled by the CIA are of the [4Fe-4S] type. GRXS17 might thus associate with the CIA-targeting complex and function as a [2Fe-2S] specific adaptor for this complex.

MATERIALS AND METHODS

Plant Material and Growth Conditions

All mutant lines used in this study were in the Col-0 ecotype background. grxs17-1 (SALK_021301), grxs17-2 (antisense line), rol5-2/ctu1-2 (GK-709D04), ctu2-2 (GK-686B10), elo3-6 (GK-555H06), met18-1 (SALK_121963), and met18-2 (SALK_147068) mutants were described previously (Leiber et al., 2010; Nelissen et al., 2010; Chen et al., 2011; Luo et al., 2012; Philipp et al., 2014). Arabidopsis (Arabidopsis thaliana) seeds were sterilized by the chlorine gas method and sown on sterile plates containing Murashige and Skoog medium supplemented with 1% (w/v) Suc, 0.8% (w/v) agarose, pH 5.7. Plates were kept 2 d in the dark for stratification at 4°C before being transferred to a growth room at 21°C with a 16-h light/8-h dark regime, with a light intensity of 80 μmol m⁻² s⁻¹, unless mentioned otherwise.

Cloning and Site-Directed Mutagenesis

The coding sequence of GRXS17, DRE2, CTU1, and CTU2 were amplified from cDNA obtained from wild-type Arabidopsis Col-0 seedlings, using specific primers (Supplemental Table S3) containing Gateway recombination attB sites. The amplified ORFs were then recombined with pDONR207 and transferred to the corresponding pDEST Gateway vectors.

Yeast Strains, Media, and SUP4 tRNA Suppressor Assay

Y2H analysis was performed as described previously (Cuéllar Pérez et al., 2013). In brief, the Saccharomyces cerevisiae PJ69-4A strain was cotransformed with pGADT7 and pGBKT7 vectors containing bait and prey, respectively. Transformants were selected on synthetic defined (SD) medium lacking Leu and Trp. Three individual transformants were grown overnight in liquid medium lacking Leu and Trp, and a 10-fold dilution of these cultures was dropped on control and selective solid media additionally lacking His. Empty vectors were used as negative controls and yeast cells were allowed to grow for 2 d at 30°C before interaction was scored.

For the conventional SUP4 tRNA suppressor assay based on ade2-1 ochre read-through (Huang et al., 2005), we used W303-1A strain (MATa ura3-1, ade2-1, trp1-1, his3-11, 15, leu2-3,112; Mortimer and Johnston, 1986) and its pGAL-GRX4/grx3Ɗ derivative (W303-1A but pGRX4::GAL-l-natNT2; grx3Ɗ::LEU2; Mühlenhoff et al., 2010). Routine yeast growth was performed in rich yeast extract, peptone dextrose (YPD), or Gal (YPG) media or on SD medium containing 2% (w/v) dextrose (SD Glu) or 2% (w/v) Gal (SD Gal; Sherman, 2002). Both lacked adenine to check for ade2-1 ochre read-through by the tRNA suppressor SUP4 and formation of adenine prototrophic cells (Jablonowski et al., 2006). The latter involved growth for 3 d at 30°C of 10-fold serial dilutions of yeast strains W303-1A and pGAL-GRX4/grx3Ɗ transformed with SUP4 on a single-copy plasmid pTC3 (Shaw and Olson, 1984; Jablonowski et al., 2009).

TAP

N or C terminally tagged TAP constructs (GS or GSrh tag) were generated as described previously (Van Leene et al., 2015). They were used for the transformation of Arabidopsis PSB-d cell suspension cultures without callus selection and further grown and subcultured as described by Van Leene et al. (2011). Stably transformed cultures were scaled up and harvested 6 d after subculturing. Transgenic Arabidopsis seeds were generated by floral dip (Clough and Bent, 1998), using Col-0 as the background ecotype and the same constructs as for cell culture transformation. Transformants were selected, and purifications were performed as described by Van Leene et al. (2015), with the exception that no benzonase treatment was performed on the cell culture extracts. Expression of TAP-tagged constructs was verified on an aliquot of total protein extract before purification.

BiFC

pCaMV35S:ORF-tag constructs using the N- or C-terminal part of GFP (nGFP and cGFP, respectively) were constructed by triple Gateway reactions using pK7m34GW or pH7m34GW (Karimi et al., 2005) as described previously (Boruc et al., 2010). pCaMV35S::tag-ORF constructs were generated by double Gateway recombination using pH7m24GW2 or pK7m24GW2 (Boruc et al., 2010). The constructs were coexpressed together with a P19-expressing Agrobacterium tumefaciens strain (Voinnet et al., 2003) in Nicotiana benthamiana using Agrobacterium-mediated transient transformation with a modified infiltration buffer (10 mm MgCl₂, 10 mm MES, pH 5.7, 100 μm acetosyringone). Interactions were scored by screening the lower epidermal cells for fluorescence using LEICA SP2 confocal microscopy, 3 d after transformation.

RNA-Seq

Seedlings were grown in vertical plates in three biological repeats for 14 d in Murashige and Skoog medium supplemented with 1% (w/v) Suc. Seedlings were frozen in liquid nitrogen, and total RNA was extracted using the RNeasy plant mini kit (Qiagen) and DNase I (Promega) treatment. A Trueseq RNA-Seq library (Illumina) was compiled and sequenced as 50-bp single read using the Illumina HiSEquation 2000 technology at GATC Biotech. Read quality control, filtering, mapping to The Arabidopsis Information Resource 10 version of the Arabidopsis genome, and read counting were carried out using the Galaxy portal running on an internal server (http://galaxyproject.org/). Sequences were filtered and trimmed with the Filter FASTQ v1 and FASTQ Quality Trimmer v1 tools, respectively, with default settings (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Reads were subsequently mapped to The Arabidopsis Information Resource 10 version of the Arabidopsis genome using GSNAPv2 (Wu and Nacu, 2010), allowing a maximum of five mismatches. The concordantly paired reads that uniquely map to the genome were used for quantification on the gene level with HTSeq-count from the HTSeq python package (Anders et al., 2015). Data were normalized using TMM, implemented in edgeR (Robinson et al., 2010), and common dispersion was then estimated using the conditional maximum likelihood method (Robinson and Smyth, 2008). Differentially expressed genes were defined by a 2-fold difference between mutant lines and the wild-type control with P-value < 0.05. The false discovery rate was limited to 5% according to Benjamini and Hochberg (1995).

Gene Expression Analysis

Seedlings were grown in the same conditions described for RNA-Seq, and total RNA was isolated as mentioned above. One microgram of RNA was used for cDNA synthesis using the iScript kit (Bio-Rad). Quantitative reverse transcriptase-PCR was performed on a LightCycler 480 system (Roche) using the Fast Start SYBR Green I PCR mix (Roche). At least three biological repeats and two technical repeats were used for each analysis. Data were analyzed using the second derivative maximum method, and relative expression levels were determined using the comparative cycle threshold method. Primer sequences are provided in Supplemental Table S3.

tRNA Extraction and Analysis

Seedlings were grown for 14 d in half-strength liquid Murashige and Skoog medium supplemented with 1% (w/v) Suc and 100 mg/L myoinositol. Seedlings were frozen in liquid nitrogen, and total RNA was extracted according to Björk et al. (2001). One-half microgram of total RNA was resolved on 8% acrylamide gels containing 0.5× TBE, 7 m urea, and 50 µg/mL APM. Northern-blot analysis was performed essentially as described in Leidel et al. (2009), with the following probes: 5′-TGGCGCCGTCTGTGGGGATCG-3′ to detect tK(UUU), 5′-TGGAGGTTCTACCGGGAGTCGAACC-3′ to detect tQ(UUG), 5′-TGGCTCCATTGCCGGGAATCGAACC-3′ to detect tE(UUC) and as negative control for nonthiolated tRNAs, 5′-TGGCACACCCGGTGGGACTCG-3′ to detect tR(UCU), and 5′-TGGTGCGTCTGCCGGGAGTCG-3′ to detect tG(UCC).

Immunoblot Analysis

After quantification of the protein content as described by Bradford (1976), the indicated protein samples were loaded on a 4% to 15% TGX gel (Bio-Rad) and ran for 20 min at 300 V. Next, proteins were transferred to 0.2-µm polyvinylidene difluoride membranes (Bio-Rad) with the Transblot Turbo (Bio-Rad). Chemiluminescent detection was performed with western Bright ECL (Isogen).

Anti-GRXS17 Antibody

GRXS17 was cloned into the pDEST17 vector (Invitrogen) and expressed in the Escherichia coli BL21(DE3) strain. The 6xHIS-GRX17 recombinant protein was purified by IMAC using Ni-NTA (Qiagen) resin followed by size exclusion chromatography using an Äkta purifier equipped with a Superdex200 column (GE Healthcare). A rabbit was immunized with the purified recombinant 6xHIS-GRX17 in the service facilities of Agro-Bio (France). Total serum was collected 63 d postimmunization. A 1/5,000 dilution of this antibody and a 1/10,000 dilution of anti-rabbit IgG-HRP (GE) were used for immunoblotting.

DNA-Damage Agent

Sterilized seeds were geminated in Murashige and Skoog medium, transferred after 4 d to Murashige and Skoog plates supplemented or not with 0.01% v/v MMS (Sigma) and scored after 17 d. The experiment was performed in triplicate.

Root Phenotype Analysis

Seedlings grown vertically at 21°C under 24-h light conditions (75 µmol m⁻² s⁻¹) were used for root analysis. The root meristem size was determined 5 DAG as the number of cells in the cortex cell file from the Quiescent center (QC) to the first elongated cell (Casamitjana-Martínez et al., 2003). The samples were mounted with clearing solution (80 g chloral hydrate, 30 mL glycerol, and 10 mL dH₂O) and observed immediately. The main root length was determined 11 DAG using ImageJ software (http://imagej.nih.gov/ij/). At least 23 seedlings of each line were analyzed.

Leaf Phenotype Analysis

Plants, 24 DAG in soil at 21°C under a 16-h light/8-h dark regime with a light intensity of 75 μmol m⁻² s⁻¹were used for leaf series on 1% agar plates, picture taking, and image analysis of leaf lamina area, length, and width as described (Cnops et al., 2004). Leaf 3 was chosen for epidermal cell imaging because of its full expansion at 24 DAG (Fig. 4C; Pyke et al., 1991; Medford et al., 1992). Leaves were fixed overnight in 100% ethanol and mounted with 90% lactic acid. The leaf area was measured with the ImageJ software. The epidermal cells on the abaxial side were drawn with a Leica DMLB microscope equipped with a drawing tube and differential interference contrast objectives. The total number of cells per leaf was calculated as described previously (De Veylder et al., 2001). We estimated the total number of cells per leaf by dividing the leaf area by the mean cell area (averaged between the tip and basal positions). Means between samples were compared by a two-tailed Student’s t test.

Enzymatic Activity Measurement (XDH)

XDH1 activity measurements in plant crude extracts were performed as described in Hesberg et al. (2004). Briefly, total protein extract was obtained from 10-d-old Arabidopsis seedlings. Plant tissue was grown in liquid nitrogen, resuspended in two volumes of extraction buffer (100 mm Tris-HCl pH 7.5, 2.5 mm EDTA, 5 mm DTT), and centrifuged. Supernatants were concentrated using Nanosep centrifugal devices (30K Omega, Pall Life Sciences) and 100 μg of total protein quantified by the method of Bradford (1976) were used for activity assays. Four to sixteen percent native polyacrylamide gels under nonreducing conditions were run at 4°C, followed by in-gel staining with 1 mm hypoxanthine, 1 mm 3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyltetrazolium bromide, and 0.1 mm phenazine methosulfate in 250 mm Tris-HCl, pH 8.5.

Metabolite Profiling by Gas Chromatography/Time of Flight-Mass Spectrometry

Metabolite analysis by gas chromatography-mass spectrometry was performed essentially as described by Lisec et al. (2006) by injecting metabolite extracts from 1 mg fresh weight of plant material into the gas chromatograph/time of flight-mass spectrometer. Chromatograms and mass spectra were evaluated using Xcalibur 2.1 software (Thermo Fisher Scientific). Specific fragment and expected elution time of target metabolites (urea, allantoic acid, allantoin, hypoxanthine, uric acid, xanthine, and xanthosine) were identified by coelution with authentic standards. Peak areas of the mass (m/z) fragments were normalized to the internal standard (ribitol) and fresh weight of the samples. Identification and annotation of detected peaks are shown in Supplemental Table S2 following recent recommendations for reporting metabolite data (Fernie et al., 2011).

Pathogen Infection

Cultivation and spore harvesting of Botrytis cinerea strain B05.10 (provided by Rudi Aerts, Katholieke Hogeschool Kempen, Belgium) was performed as described previously (Broekaert et al., 1990). Arabidopsis wild-type (Col-0) and mutants plants were grown for 4 weeks in soil (“DCM potgrond voor Zaaien en Stekken”; DCM, Sint-Katelijne-Waver, Belgium) in a growth chamber at 21°C, 75% humidity, and a 12-h day/light cycle with a light intensity of approximately 120 μmol m⁻² s⁻¹. A 5-μL drop of a B. cinerea spore suspension (5 × 10⁵/mL in half-strength potato dextrose broth) was inoculated on three leaves per plant. Plants were kept in transparent, sealed boxes to retain almost 100% humidity after inoculation. Disease symptoms were scored by measuring the diameter of the necrotic lesions at 2 and 3 d postinoculation. Thirty-two plants per line were analyzed. Two independent assays with similar results were performed.

Accession Numbers

Accession numbers of the genes used in this study are as follows: GRXS17, AT4G04950; MET18, AT5G48120; DRE2, AT5G18400; NAR1, AT4G16440; CIA1, AT2G26060; AE7/CIA2, AT1G68310; XDH1, AT4G34890; URH1, AT2G36310; URH2, AT1G05620; UREG, AT2G34470; CTU1, AT2G44270; CTU1L, AT1G76170; CTU2, AT4G35910; ELO3, AT5G50320; Leu-tRNA ligase, AT1G09620; PARP2, AT4G02390; ACC2, AT1G36180; GMI, AT5G24280; F4P13.14, AT3G01600; AAA-ATPase, AT2G18193; and DTX3, AT2G04050 and AT5G60250.

Supplemental Data

The following supplemental materials are available.

• Supplemental Figure S1. GRXS17 interacts with DRE2 and XDH1 in planta.

• Supplemental Figure S2. Characterization of the grxs17 mutants used in this study.

• Supplemental Figure S3. Analysis of the functioning of the purine salvage pathway in grxs17-1 seedlings.

• Supplemental Figure S4. BIFC controls.

• Supplemental Figure S5. grxs17 and elo3-6 mutants have similar developmental defects.

• Supplemental Figure S6. grxs17 and elo3-6 mutants show similar molecular defects.

• Supplemental Table S1. RNA-Seq analysis of the grxs17-1 mutant.

• Supplemental Table S2. Identification and annotation of detected purine metabolite peaks.

• Supplemental Table S3. Primers used in this study.

• Supplemental Dataset S1. Mass spectrometry data of TAP experiments.

ACKNOWLEDGMENTS

We thank Florian Bittner for providing the XDH1 clone, Frederik Coppens for support with the RNA-Seq analysis, Pia Neyt for technical support, and Annick Bleys for help in preparing the manuscript.

Glossary

 
• GRX

  glutaredoxin

 
• CIA

  cytosolic Fe-S assembly

 
• TRX

  thioredoxin

 
• TAP

  tandem affinity purification

LITERATURE CITED

Author notes