{'fraction_non_alphanumeric': 0.042628966112669975, 'fraction_numerical': 0.03615368012087201, 'mean_word_length': 5.17852975495916, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 1, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
A facile forward-genetic screen for Arabidopsis autophagy mutants reveals twenty-one loss-of-function mutations disrupting six ATG genes

ABSTRACT Macroautophagy is a process through which eukaryotic cells degrade large substrates including organelles, protein aggregates, and invading pathogens. Over 40 autophagy-related (ATG) genes have been identified through forward-genetic screens in yeast. Although homology-based analyses have identified conserved ATG genes in plants, only a few atg mutants have emerged from forward-genetic screens in Arabidopsis thaliana. We developed a screen that consistently recovers Arabidopsis atg mutations by exploiting mutants with defective LON2/At5g47040, a protease implicated in peroxisomal quality control. Arabidopsis lon2 mutants exhibit reduced responsiveness to the peroxisomally-metabolized auxin precursor indole-3-butyric acid (IBA), heightened degradation of several peroxisomal matrix proteins, and impaired processing of proteins harboring N-terminal peroxisomal targeting signals; these defects are ameliorated by preventing autophagy. We optimized a lon2 suppressor screen to expedite recovery of additional atg mutants. After screening mutagenized lon2-2 seedlings for restored IBA responsiveness, we evaluated stabilization and processing of peroxisomal proteins, levels of several ATG proteins, and levels of the selective autophagy receptor NBR1/At4g24690, which accumulates when autophagy is impaired. We recovered 21 alleles disrupting 6 ATG genes: ATG2/At3g19190, ATG3/At5g61500, ATG5/At5g17290, ATG7/At5g45900, ATG16/At5g50230, and ATG18a/At3g62770. Twenty alleles were novel, and 3 of the mutated genes lack T-DNA insertional alleles in publicly available repositories. We also demonstrate that an insertional atg11/At4g30790 allele incompletely suppresses lon2 defects. Finally, we show that NBR1 is not necessary for autophagy of lon2 peroxisomes and that NBR1 overexpression is not sufficient to trigger autophagy of seedling peroxisomes, indicating that Arabidopsis can use an NBR1-independent mechanism to target peroxisomes for autophagic degradation. Abbreviations: ATG: autophagy-related; ATI: ATG8-interacting protein; Col-0: Columbia-0; DSK2: dominant suppressor of KAR2; EMS: ethyl methanesulfonate; GFP: green fluorescent protein; IAA: indole-3-acetic acid; IBA: indole-3-butyric acid; ICL: isocitrate lyase; MLS: malate synthase; NBR1: Next to BRCA1 gene 1; PEX: peroxin; PMDH: peroxisomal malate dehydrogenase; PTS: peroxisomal targeting signal; thiolase: 3-ketoacyl-CoA thiolase; UBA: ubiquitin-associated; WT: wild type


Introduction
Eukaryotes use macroautophagy, hereafter referred to as autophagy, to maintain cellular homeostasis during development and in response to environmental changes. This degradation process targets protein aggregates, pathogens, damaged or superfluous organelles, and other cellular components [1][2][3][4][5][6]. During autophagy, a double-membrane structure known as a phagophore surrounds and engulfs the fated cargo. Once enclosed, the outer membrane of the resulting double-membrane-bound autophagosome fuses with a lysosome (in metazoans) or vacuole (in plants and fungi) to deliver the cargo into the lysosome or vacuole lumen. The cargo is broken down, and constituent amino acids, lipids, and carbohydrates are returned to the cytosol, where these molecules can be recycled to synthesize new organelles and proteins.
In plants, autophagy promotes survival during starvation and harsh environmental conditions. Autophagy is upregulated in plants during nutrient deprivation [7], and plants lacking autophagy are hypersensitive to carbon and nitrogen starvation [8], fungal infection [9], and various abiotic stresses, including heat stress [10], drought and salt stress [11], and oxidative stress [12]. Autophagy is also upregulated during seed maturation [13], suggesting a role during embryo development. Although autophagy is not required for plant survival in optimal growth conditions [8], plants lacking autophagy accumulate the selective autophagy receptor NBR1 [14] and certain peroxisomal proteins [15,16], suggesting that basal autophagy also contributes to general intracellular processes.
A variety of receptors with differing cargo specificity mediate selective autophagy by controlling which substrates are engulfed by autophagosomes. Selective autophagy receptors link organelles, protein aggregates, or other cargo to the autophagy machinery by binding to both the fated cargo and ATG8 [33]. These receptors generally bind ATG8 via an ATG8-interacting motif with the consensus core sequence [W/Y/F]xx[L/I/V] with neighboring acidic (D or E) residues [33,34]. The ability to recognize not only specific types of organelles but also damaged or unnecessary organelles is paramount to maintaining cellular homeostasis, and the molecular mechanisms that enact various selective autophagy pathways continue to be unraveled.
Several selective autophagy receptors have been characterized [4,33]. For example, mammalian NBR1 and SQSTM1/p62 (sequestosome 1) bind both LC3/ATG8 and ubiquitin [33]. Plants encode NBR1 orthologs (known as Joka2 in tobacco), but plant orthologs of SQSTM1/p62 have not been identified [14,35]. Arabidopsis nbr1 mutants display sensitivity to only a subset of the stressors to which atg mutants display sensitivity [10,36], suggesting that NBR1 acts as a selective autophagy receptor for a subset of autophagic cargos. For example, Arabidopsis NBR1 is implicated in clearing ubiquitinated protein aggregates following heat stress [10] and in limiting cauliflower mosaic virus infection by targeting virus particleforming capsid proteins for autophagic degradation [37] but appears not to participate in autophagy of proteasomes [38]. Several plant-specific ATI (ATG8-interacting) proteins [39,40] are implicated in autophagy. For example, ATI1/At2g45980 plays a role in autophagy of plastids [41], and ATI3A/ At1g17780 is implicated in autophagy of the endoplasmic reticulum (ER) [40]. Besides ATI proteins, Arabidopsis uses a ubiquitin receptor to target a transcription factor for autophagic degradation [42] and uses a proteasome subunit as a receptor for autophagy of proteasomes [38].
One target of autophagy is the peroxisome, an organelle that sequesters oxidative reactions. Most notably, peroxisomes house fatty acid β-oxidation, which converts fatty acids into acetyl-CoA and produces hydrogen peroxide (H 2 O 2 ) as a byproduct [43]. Peroxisomes use various enzymes to detoxify H 2 O 2 , including catalase [44], which is itself susceptible to oxidative damage by H 2 O 2 [45]. Pexophagy, the selective autophagy of peroxisomes, contributes to peroxisome homeostasis in plants by degrading damaged or obsolete peroxisomes [46]. Pexophagy occurs at a higher basal rate than autophagy of other organelles in Arabidopsis seedlings, as evidenced by the higher relative accumulation of peroxisomal versus other organellar proteins in atg5 mutants [15,47]. Moreover, Arabidopsis atg mutants accumulate peroxisomes [16] and nonfunctional, oxidized catalase [15,47], suggesting a role for pexophagy in maintaining peroxisome function. However, the molecular mechanisms by which plant cells recognize peroxisomes in need of turnover remain unclear.

Forward-genetic screening strategy for Arabidopsis autophagy mutants
We used a 4-step strategy ( Figure 1) to recover 21 mutants carrying mutations in 6 different Arabidopsis ATG genes (Table 1). First, to identify potential atg mutants, we screened approximately 500,000 8-day-old lon2-2 M 2 seedlings from 104 mutagenized pools for the presence of lateral roots on medium containing 8 µM IBA (Figure 1(a)). Unlike wild-type seedlings, lon2 seedlings fail to form abundant lateral roots in response to IBA [49], and preventing autophagy in lon2 mutants suppresses this phenotype [48]. . Four-step strategy for isolating Arabidopsis autophagy-defective mutants. (a) Initial screening for lon2 suppressors. lon2 forms few lateral roots in response to IBA; putative suppressors that formed lateral roots in the presence of IBA similar to wild-type and lon2-2 atg7-4 were moved to soil for propagation. Top panel: wild-type (WT), lon2-2, and lon2-2 atg7-4 seedlings were grown on media containing 8 µM IBA and imaged at 8 days; scale bar: 5 mm. Bottom row: magnified images of roots outlined in the top panel showing lateral roots, which are absent in lon2-2; scale bar: 1 mm. (b) Secondary screening of putative suppressors. Leaf extracts from approximately 30-day-old adult controls (left of dashed line) or M 2 putative suppressor plants (right of dashed line) were processed for immunoblotting with antibodies to the indicated proteins. NBR1 is a selective autophagy receptor that accumulates in atg mutants [14]. PMDH is synthesized as a precursor (p) that is processed to a mature form (m) in the peroxisome. HSC70 is a loading control. The asterisk indicates a protein cross-reacting with the ATG7 antibody. (c) Retesting progeny of putative suppressors that displayed restored PTS2 processing in the M 2 generation. Top: lateral root density of 8-day-old controls (left of dashed line) or M 3 or M 4 suppressor seedlings (right of dashed line) grown without or with IBA. Error bars show standard deviations (n = 8). Statistically significant (P < 0.0001) differences determined by one-way ANOVA are depicted by different letters above the bars. Bottom: extracts from 6-day-old controls (left of dashed line) or M 3 or M 4 suppressor seedlings (right of dashed line) were processed for immunoblotting. Membranes from duplicate gels were serially probed with antibodies to the indicated proteins to obtain the top 3 and bottom 4 panels. MLS, thiolase, and ICL are peroxisomal proteins that are stabilized when both LON2 and autophagy are defective [48]. Asterisks indicate proteins cross-reacting with the ATG7 or ICL antibodies. (d) Identifying mutations via whole-genome sequencing. The L40 suppressor was backcrossed to the original lon2-2 line, IBA-sensitive F 2 seedlings were selected, and genomic DNA from pooled F 3 seedlings was sequenced. Homozygous singlenucleotide polymorphisms consistent with EMS mutagenesis (G/C to A/T transitions) and causing nonsynonymous mutations in coding regions, altering splice sites, or occurring in introns or untranslated regions are indicated by locus identifiers and displayed to the right of the 5 Arabidopsis chromosomes using The Arabidopsis Information Resource Chromosome Map Tool.
Second, to eliminate false positives and mutants with elevated auxin levels or responsiveness, we moved seedlings with lateral roots to soil and collected leaf tissue from approximately 30-day-old plants for immunoblot analysis. We used antibodies recognizing NBR1, PMDH (peroxisomal malate dehydrogenase/At5g09660), and ATG7 to identify suppressors with impaired autophagy and improved peroxisome function ( Figure 1(b)). The selective autophagy receptor NBR1 is degraded at a basal rate via autophagy and accumulates when autophagy is impaired [14]. PMDH is synthesized as precursor containing an N-terminal peroxisomal-targeting signal 2 (PTS2) that is removed following peroxisomal import. The incomplete PTS2 processing in lon2 plants [49] presumably stems from elevated pexophagy, which reduces the number of peroxisomes; preventing autophagy restores PMDH processing in lon2-2 by increasing the number of importcompetent peroxisomes [48]. M 2 plants with elevated NBR1 levels and restored PMDH processing were prioritized as likely atg mutants, and plants with decreased ATG7 levels were candidate atg7 mutants (Figure 1(b)).
Third, to confirm lon2-2 suppression, we collected M 3 or M 4 progeny of putative suppressors; retested IBA responsiveness; and examined thiolase, MLS or ICL, and ATG3 levels in 6-day-old seedlings ( Figure 1(c)). Thiolase is a peroxisomal enzyme that is destabilized in lon2-2 seedlings [49] and stabilized in lon2-2 atg mutants [48]. Similarly, the peroxisomal enzymes ICL and MLS are stabilized in lon2-2 atg double mutants [48]. Mutants with elevated thiolase and ICL or MLS were prioritized as likely atg mutants, and plants with decreased ATG3 levels were candidate atg3 mutants.
Fourth, we used DNA sequencing to identify causal lesions. We sequenced ATG7 or ATG3 in suppressors with decreased levels of the respective proteins and prepared DNA for wholegenome sequencing for the remaining suppressors ( Figure 1(d)).
New atg3, atg5, and atg16 alleles recovered as lon2 suppressors In addition to ATG7, we isolated novel alleles of 3 other genes involved in ATG8 lipidation. We identified the atg3-2 splicesite mutation by sequencing ATG3, which encodes an E2-like enzyme that conjugates ATG8 to phosphatidylethanolamine [32], from a mutant that lacked detectable ATG3 protein ( Figure 3(c)). We recovered 2 mutations in ATG5, which encodes a protein with 2 ubiquitin-like domains that is conjugated by the ubiquitin-like protein ATG12 to act as an E3-like enzyme in ATG8 lipidation [29,31,62]. We used wholegenome sequencing to identify the atg5-5 splice-site mutation (Figure 1(d), 3(a)) and the atg5-6 missense mutation (Figure 3 (a), Fig. S1A), which altered a conserved acidic residue (Glu323) in the second ubiquitin-like domain to a Lys residue (Fig. S3). We did not detect ATG5 protein in either of these mutants (Figure 3(c)).
We recovered 3 mutations in ATG16, which in yeast encodes a protein that tethers the ATG12-ATG5 conjugate to the membrane to allow ATG8 lipidation and autophagosome formation [31]. Whole-genome sequencing revealed that atg16-1 and atg16-3 harbored nonsense mutations (W47Stop and Q97Stop, respectively) and that atg16-2 harbored both a splice-site mutation and a nearby missense mutation (R80K) (Figure 3(a), Fig. S1C to E). We surmise that the atg16-2 splice site mutation is causal because it is expected to prevent full-length ATG16 accumulation and is 5ʹ of the missense mutation.
Because the roles of ATG3 and ATG16 have not been confirmed in plants, we began this analysis by comparing ATG8 levels in an allele of each atg mutant in a wild-type LON2 background. As previously reported [29], ATG8 levels were elevated in atg5 and atg7 mutants (Figure 3 (d)). We found similarly elevated ATG8 levels in our atg3 and atg16 alleles (Figure 3(d)), which is consistent with the possibility that autophagy is impaired in these mutants.
New atg2 alleles recovered as lon2 suppressors We used whole-genome sequencing to identify 2 mutations in ATG2 (Figure 4(a, b), Fig. S4), which encodes a protein involved in phagophore expansion. Both atg2 mutants fully restored IBA responsiveness to lon2-2 ( Figure 4(d)) and displayed the characteristic immunoblot phenotypes (Figure 4 (e)) of other lon2-2 atg seedlings (Figure 1(c), 2(c), 3(c)). Moreover, atg2-6 displayed elevated ATG8 levels ( Figure 3 (d)). The atg2-7 early nonsense mutation (W90Stop, Figure 4 (b)) is expected to prevent ATG2 accumulation. The atg2-6 mutation appeared in the middle of a cluster of mutations on chromosome 3 (Figure 4(a)) identified by sequencing backcrossed seedlings. However, the mutation was in intron 6 according to the primary prediction of Araport 11, the latest A. thaliana genome annotation. Araport 11 also predicts 2 ATG2 splice variants that would result in atg2-6 interrupting a splice-acceptor site. Indeed, the region of the intron following the mutated nucleotide can be translated in frame with exon 7 and lacks stop codons (Figure 4(b)). Moreover, aligning predicted ATG2 proteins of other Brassicaceae family members revealed that this region of the predicted Arabidopsis intron 6 is conserved and exonic in other Brassicaceae ATG2 orthologs (Figure 4(b)), supporting the conclusion that the mutation in atg2-6 disrupts a splice site in A. thaliana and suggesting that the 'alternative' splice variant encodes an ATG2 isoform that is essential for autophagy. Although the predicted splicing of A. lyrata ATG2 matched that of the A. thaliana ATG2 predicted by Araport 11, the A. lyrata genome was assembled based in part on alignment with A. thaliana [63], meaning that annotation errors in the A. thaliana genome may have been propagated to A. lyrata. Based on the whole-genome sequencing and Brassicaceae ATG2 alignment, we considered atg2-6 as a splice-site mutation and the causal suppressor lesion.
To further probe the importance of ATG18a in pexophagy, we backcrossed lon2-2 atg18a-3 to remove lon2-2 and other unlinked mutations and examined atg18a-3 seedlings constitutively expressing GFP carrying a peroxisomal targeting signal (GFP-PTS1) [66] using confocal microscopy. atg7 and atg5 mutants accumulate peroxisomes in hypocotyls of developing seedlings [16], and atg18a mutants display clusters of oxidized peroxisomes in mesophyll cells [15]. We observed increased peroxisome abundance relative to wild type in atg18a-3 cotyledon epidermal cells, cotyledon mesophyll cells, and hypocotyl cells in 4-to 7-day-old seedlings ( Figure 5), confirming that ATG18a promotes pexophagy during seedling development. , and atg16 mutations identified as lon2 suppressors are shown in pink, red, and orange, respectively; the previously described atg3-1 EMS-derived lon2 suppressor [48] is in gray, and T-DNA insertion alleles [29,72,73] are indicated by triangles. aa, amino acids. (b) Lateral root density of 8-day-old wild type (WT), lon2-2, atg3-1, lon2-2 atg3-2, lon2-2 atg5, and lon2-2 atg16 seedlings grown without or with IBA. Error bars show standard deviations (n = 8). Statistically significant (P < 0.0001) differences determined by one-way ANOVA are depicted by different letters above the bars. (c) Extracts from 6-day-old seedlings were processed for immunoblotting. Membranes from duplicate gels were serially probed with antibodies to the indicated proteins to obtain the top 5 and bottom 3 panels. The asterisk indicates a protein cross-reacting with the ATG5 antibody. (d) Extracts from 6-day-old seedlings were processed for immunoblotting and serially probed with antibodies to NBR1, ATG8A, and HSC70. ATG8 is lipidated and encoded by multiple genes in Arabidopsis.  [48] are in gray, and a T-DNA insertion allele [64] is indicated by a triangle. The partial alignment shows predicted Brassicaceae ATG2 proteins, including 2 A. thaliana ATG2 splice variants predicted by the latest genome annotation (Araport 11); the alternative ATG2 splice variant is interrupted by the atg2-6 mutation. aa, amino acids. (c) ATG18a gene diagram showing atg18a mutations identified as lon2 suppressors (teal) and a previously described T-DNA insertion allele [65] (triangle). (d) Lateral root density of 8-day-old wild type (WT), lon2-2, atg2-4, lon2-2 atg2, and lon2-2 atg18a seedlings grown without or with IBA. Error bars show standard deviations (n = 8). Statistically significant (P < 0.0001) differences determined by one-way ANOVA are depicted by different letters above the bars. (e) Extracts from 6-day-old seedlings were processed for immunoblotting with antibodies to the indicated proteins.
NBR1 is not necessary for pexophagy of lon2 peroxisomes NBR1 is a selective autophagy receptor in Arabidopsis [10,14,36] and is necessary and sufficient for pexophagy in mammals [67], but no nbr1 mutants emerged from our screen for lon2-2 suppressors. To test the necessity of NBR1 for Arabidopsis pexophagy, we obtained T-DNA insertional alleles of nbr1 (Figure 7(a)). These T-DNA insertions dramatically decreased full-length NBR1 protein accumulation (Figure 7(c)); however, the anti-NBR1 antibody was generated to the C-terminal UBA domain [14] and would not detect any N-terminal fragments of the protein that might accumulate in the mutants. Because intronic T-DNAs are sometimes removed by splicing [68,69], we characterized both intronic (nbr1-1) and exonic (nbr1-4) insertional mutants of NBR1 in combination with lon2-2. Both lon2 nbr1 double mutants resembled lon2-2 and did not form lateral roots in response to IBA (Figure 7(b)). Moreover, neither nbr1 allele stabilized MLS or thiolase or ameliorated the PTS2-processing defects of lon2 (Figure 7(c)). This lack of suppression suggested that NBR1 is dispensable for pexophagy in lon2. Because atg18a-3 accumulated peroxisomes ( Figure 5), we also visualized nbr1-4 expressing the peroxisomal marker GFP-PTS1 to evaluate peroxisome abundance. Unlike atg18a-3, peroxisome abundance in hypocotyls of 6-day-old nbr1-4 seedlings resembled wild type (Figure 7(d)). The failure of nbr1-4 mutants to accumulate excess peroxisomes suggested that NBR1 is not necessary for peroxisome turnover under our standard growth conditions.
If excess NBR1 was sufficient to induce pexophagy, then we would expect that overexpressing NBR1 would stimulate pexophagy, as is observed in mammalian cells [67], and phenocopy lon2. We drove untagged, HA-tagged, and YFP-tagged NBR1 from the constitutive cauliflower mosaic virus 35S promoter. Despite elevated NBR1 levels (Figure 7(e)), these lines resembled wild type in lateral root formation (Figure 7(b)), PTS2 processing (Figure 7(e)), and thiolase stability (Figure 7(e)). This failure of excess NBR1 to potentiate seedling pexophagy indicates that NBR1 is not limiting for pexophagy in Arabidopsis, which is consistent with the possibility that NBR1 is not a pexophagy receptor in Arabidopsis seedlings.

A facile screen for Arabidopsis atg mutants
Screening for lon2-2 suppressors revealed 21 alleles (20 novel alleles) in 6 ATG genes (Table 1), providing new mutants for investigating autophagy in Arabidopsis. A pilot lon2-2 suppressor screen has uncovered atg2, atg3, and atg7 mutants [48], and the expanded screen described here uncovered additional alleles of these genes plus atg5, atg16, and atg18a mutants. In contrast to these lon2 suppressor screens, a single atg2 allele emerged from a forward-genetic screen for plants with enhanced cell death in response to fungal infection [22], a single atg2 allele has emerged from a hydroxyurea-resistance screen [70] that primarily recovers cat2 (catalase/At4g35090) mutants [71], and 2 atg2 alleles along with single atg7 and atg18a mutants have been recovered from a microscopy-based screen for aggregated peroxisomes [15]. The lon2-2 suppressor screen combines high throughput with specificity and is the only forward-genetic screen to report atg3, atg5, or atg16 mutants. Although atg5 T-DNA insertional mutants have been characterized [29,72,73], no Arabidopsis atg3 or atg16 mutants have been reported aside from those arising from lon2-2 suppressor screens described here and previously [48], and these alleles will be useful for future studies to determine whether ATG3 and ATG16 function similarly in Arabidopsis as in other organisms. Moreover, our atg7 missense alleles (Figure 2(a), Fig. S2) may inform ATG7 structure-function studies.
Roughly half (11 of 21) of the lon2-2 suppressors harbored a mutation in ATG7 (Table 1); the reason for this abundance of atg7 mutants is mysterious. When adding the atg7 alleles identified in our pilot screen [48], 16 of 29 lon2-2 suppressors were atg7 mutants, and atg7-6 was independently isolated in both iterations of the screen (Table 1) [48]. Although ATG7 encodes a large protein (697 amino acids), ATG2 and ATG11 encode larger proteins (1861 and 1148 amino acids, respectively) and thus would seem to be more likely mutagenic targets. However, our combined screening recovered only 4 atg2 and no atg11 mutants. Moreover, ATG2 has 12 introns compared to 10 introns in ATG7, suggesting that ATG2 would present more opportunities than ATG7 for mutations in splice sites in addition to coding regions. The frequency of atg7 mutants recovered in this screen might be explained by the structure and function of the ATG7 protein. As an E1-like enzyme, ATG7 activates both ATG8 and ATG12 and coordinates transfer to their respective E2-like enzymes, ATG3 and ATG10 [32,[74][75][76]. ATG7 is the interaction hub among these proteins, requiring precise spatial arrangement [32,[74][75][76][77], so it is tempting to speculate that ATG7 might be more susceptible to functional perturbation via missense mutations than other ATG proteins. Indeed, 5 of the 16 atg7 mutants identified as lon2-2 suppressors (atg7-7, atg7-8, atg7-11, atg7-17, and atg7-19; Figure 2(a), Fig. S2) harbor missense mutations whereas our only non-atg7 lon2-2 suppressor to harbor an atg missense mutation is atg5-6 ( Figure 3(a), Fig. S3). However, fewer mutations in other ATG genes were recovered, hindering comparison of the relative proportions of null versus missense alleles; further screening might recover missense alleles of other ATG genes. Moreover, even discounting missense mutations, 11 of our atg7 mutants harbor nonsense or splice-site mutations, which still exceeds the number of alleles recovered for any other ATG gene (Table 1), indicating that the explanation for the apparent bias towards atg7 mutants may exist beyond the structure or function of ATG7. Perhaps the ATG7 locus is more accessible to EMS mutagenesis than other ATG genes.

function in autophagy remains an open question.
Because no atg11 mutants emerged from our screen for lon2-2 suppressors, we crossed the atg11-1 T-DNA insertional allele to lon2-2 and observed only partial suppression of lon2-2 phenotypes (Figure 6), suggesting that pexophagy is incompletely blocked in atg11-1. This conclusion is supported by the observation that atg11-1 accumulates less catalase than other atg mutants [18]. Moreover, ATG8 and NBR1 are only partially stabilized in atg11-1 relative to atg7-4 and other null alleles in core ATG genes (Figure 3(d), 6(c)), consistent with previous reports that autophagy is incompletely blocked in atg11-1 [17,18]. As ATG11 is thought to act as a scaffold for ATG1a and ATG13 within the phagophore assembly site [17], these findings imply that autophagy can still occur, albeit less efficiently, in atg11-1 seedlings. Despite the large size of ATG11 (1148 amino acids), its role as a scaffold may allow some tolerance for missense mutations, which, combined with the observation that atg11-1 only partially restored lon2-2 IBA responsiveness ( Figure 6(b)), may explain why atg11 mutants did not emerge from lon2-2 suppressor screens. In Pichia pastoris, Atg11 binds both Atg8 and the pexophagy receptor Atg30 [78], so it will be interesting to determine whether ATG11 plays a similar role in Arabidopsis pexophagy. Arabidopsis ATG11 also interacts with ATG8 [17], and additional studies are needed to decipher the role of this protein in processes upstream and downstream of ATG8 lipidation.

Arabidopsis pexophagy receptors remain to be identified
In spite of the conservation of autophagy across kingdoms, identified pexophagy receptors differ among species, and a pexophagy receptor has not been reported in plants. Pichia pastoris Atg30 [86] and Saccharomyces cerevisiae Atg36 [87] act as pexophagy receptors, but plants lack Atg30 or Atg36 orthologs. We did not recover mutants of candidate pexophagy receptors in our screen for lon2-2 suppressors. Because NBR1 is a pexophagy receptor in mammalian cells [67] and a selective autophagy receptor in Arabidopsis [10,14,36,37], we used reverse genetics to investigate the role of Arabidopsis NBR1 in pexophagy. Two T-DNA insertional alleles of nbr1 failed to suppress lon2-2 defects (Figure 7(b, c)), and NBR1 overexpression failed to phenocopy lon2 seedling defects (Figure 7(b, e)), which presumably result from heightened pexophagy. Moreover, peroxisome abundance in nbr1-4 seedlings resembled wild type rather than the increased abundance of peroxisomes observed in atg18a-3 (Figure 7(d)). These data indicate that NBR1 is not the sole pexophagy receptor in lon2 seedlings but do not eliminate the possibility that NBR1 acts redundantly in pexophagy or in other contexts. Intriguingly, peroxisomal catalase accumulates in nbr1 mutants during heat stress [36], suggesting that NBR1 might be important for LON2-independent pexophagy in Arabidopsis. Further studies are needed to fully elucidate the roles of NBR1 in Arabidopsis selective autophagy.
Perhaps plants employ a pexophagy receptor that is distinct from those identified in yeast and mammalian cells. The plant-specific ATI (ATG8-interacting) proteins [39,40] are pexophagy receptor candidates but did not emerge from our screen for lon2-2 suppressors. The ubiquitin-binding protein DSK2 [88] is another pexophagy receptor candidate in Arabidopsis. Arabidopsis DSK2 binds to and targets the transcription factor BES1 for autophagic degradation [42]. DSK2 also interacts with 2 peroxisomal membrane proteins [89]. Because Arabidopsis has 2 DSK2 paralogs (At2g17190 and At2g17200) [90], recovering dsk2 mutants as lon2 suppressors is unlikely, and directed approaches will be needed to explore the potential roles of DSK2 and ATI proteins in pexophagy.
Potential roles for peroxins in pexophagy PEX (peroxin) proteins are necessary for peroxisome biogenesis and matrix protein import [91]; many peroxins are found on the peroxisome membrane where they would be accessible to the autophagy machinery. Several peroxins have been implicated in pexophagy in Arabidopsis as well as yeast and mammalian cells. In general, peroxins or other peroxisomal membrane proteins might mediate pexophagy in 3 ways.
Third, peroxins could interact directly with ATG8, bypassing the need for a bridging receptor. Although the relative simplicity of ATG8-interacting motifs has confounded bioinformatic approaches to identifying ATG8-interacting proteins, a recent refinement considering acidic residue positioning reveals that Arabidopsis PEX10/At2g26350 and PEX6/ At1g03000 interact with ATG8 [34]. Interestingly, pexophagy contributes to defects observed in pex1 and pex6 mutants in Arabidopsis [98,99], yeast [100], and humans [101], suggesting that PEX6 prevents rather than promotes pexophagy, perhaps by competing with the pexophagy signal for interaction with ATG8 or by removing ubiquitinated substrates, such as PEX5, from the peroxisomal surface.
Despite the potential importance of peroxins in pexophagy, pex mutants are unlikely to emerge as lon2-2 suppressors because pex mutants are generally IBA resistant [91] and thus not expected to restore IBA-responsive lateral rooting to lon2 (Figure 1(a)). Alternative approaches will be needed to dissect the roles of peroxins in Arabidopsis pexophagy.

Conclusions
We used forward genetics to recover 21 alleles of 6 ATG genes. All 6 genes are core ATG genes, underscoring the functional conservation of these components. These mutant alleles will be useful in future studies interrogating the role of these genes in Arabidopsis autophagy, especially the atg3 and atg16 alleles, which disrupt genes without available T-DNA insertional alleles. Using the lon2 mutant as a system for studying pexophagy in Arabidopsis, we also investigated the role of ATG11 and NBR1 in this process and found that ATG11 plays a limited role in pexophagy whereas NBR1 does not appear to be involved in lon2related pexophagy. lon2 mutants will undoubtedly serve as a useful platform for future studies elucidating the regulation of pexophagy in Arabidopsis.
nbr1-4 (WiscDsLoxHs007_07A) was obtained from the ABRC. The position of the T-DNA insertion was confirmed by PCR of genomic DNA using a gene-specific primer and a T-DNA left-border specific primer (Table S1), and the exact position of the insertion was determined by sequencing the PCR product (Lone Star Labs, Houston, TX, USA). nbr1-4 35S:GFP-PTS1 was obtained by crossing.
Seeds were surface-sterilized using 3% (w:v) sodium hypochlorite solution containing 0.01% (v:v) Triton X-100 (Sigma-Aldrich, T9284), resuspended in 0.1% agar (Fisher Scientific, BD-214030), and stratified at 4°C overnight. Seedlings were grown at 22°C on plant nutrient (PN) medium [106] supplemented with 0.5% (w:v) sucrose (VWR, JT4097-6) (PNS) and solidified with 0.6% (w:v) agar. To quantify IBA sensitivity, seedlings were grown on PNS plates under continuous white light for 4 days, transferred to new PNS plates with either ethanol (mock) or 10 µM IBA (Sigma-Aldrich, I5386; made from a 100-mM stock in ethanol), and grown under light filtered with yellow long-pass filters (to slow photochemical breakdown of indolic compounds) [107] for an additional 4 days. Lateral roots emerged from the primary root were counted using a stereomicroscope (S6E, Leica Microsystems, Buffalo Grove, IL, USA), and primary root lengths were measured using a ruler. Measurements were completed at least twice with similar results.

Mutant isolation
lon2-2 35S:GFP-PTS1 seeds (1.3 g,~65,000 seeds) were mutagenized for 16 h at room temperature with 0.24% to 0.25% (v:v) ethylmethane sulfonate (Sigma-Aldrich, M0880), rinsed extensively with water, and grown in 104 M 1 pools. For screening, 0.1 g of M 2 seeds (~5000 seeds) per pool were surface-sterilized, stratified, and plated on PNS medium supplemented with 8 µM IBA and solidified with 1.0% (w:v) agar in square Petri dishes (89 x 89 mm). Plates were positioned vertically, and seedlings were grown under yellow-filtered light for 8 days, when M 2 seedlings with 3 or more lateral roots were moved to soil. Leaf tissue was collected from approximately 30-day-old M 2 plants for immunoblot analysis using anti-PMDH, anti-NBR1, anti-ATG7, and anti-HSC70 antibodies (as described below). M 3 seeds from M 2 plants with multiple lateral roots, complete PTS2 processing, and elevated NBR1 protein levels were collected and re-tested for IBA resistance (as described above in Plant Materials and Growth Conditions), and 6-day-old tissue was tested by immunoblot analysis using anti-thiolase, anti-ATG3, anti-HSC70, and anti-ICL or anti-MLS (as described below). Lines exhibiting IBA sensitivity and thiolase and ICL or MLS stabilization were retained as lon2 suppressors. All suppressors were genotyped to confirm lon2-2 homozygosity (Table S1). After suppressor mutation identification, PCR-based genotyping markers designed using dCAPS Finder 2.0 [108] were used to confirm the presence of the identified lesions (Table S2) and follow mutations in the progeny of crosses.
Genomic DNA for whole-genome sequencing was prepared as previously described [109]. DNA was sequenced with HiSequation 2000 sequencers (Illumina, San Diego, CA, USA) at the Genome Technology Access Center at Washington University in St. Louis and aligned with the TAIR 10 build of the A. thaliana Col-0 genome using Novoalign (Novocraft; http://novocraft.com). SNPs were identified using SAMtools [110] and annotated with snpEFF [111]. Mutations were filtered using a script that prioritized homozygous canonical EMS-derived mutations (G-to-A and C-to-T) causing non-synonymous amino acid changes or altered splice sites but also retained mutations in introns and heterozygous EMS-consistent mutations. We disregarded mutations that were present in our lab stock of Col-0 or in multiple suppressors from different pools, indicating an origin in the starting line. Chromosomal nucleotide positions and gene identifiers of genes with homozygous EMS-consistent mutations are listed in Data S1.
Sequencing genomic DNA from the L1 isolate revealed mutations in both ATG18a (c463t causing Q155Stop) and ATG18f (g708a causing A198T). To determine which of these lesions was linked to lon2 suppression, we crossed L1 to lon2-2 35S:GFP-PTS 1 and used genotyping markers (Table S2) to isolate F 3 plants that were homozygous for the atg18a lesion but homozygous for wild-type ATG18f and vice versa. The resultant F 4 progeny were tested for IBA responsive lateral rooting and 6-day-old immunoblot phenotypes to reveal that the atg18a mutation suppressed lon2-2 whereas the atg18f mutation did not alter lon2-2 phenotypes. We used the backcrossed lon2-2 atg18a-3 mutant in subsequent analyses.

Individual gene sequencing
ATG7 or ATG3 was sequenced directly in suppressors with reduced ATG7 or ATG3 protein levels, respectively. ATG7 was PCR amplified using the primers indicated in Table S3, and ATG3 was PCR amplified using previously described primers [48]. Amplicons were purified using the DNA Clean & Concentrator kit (Zymo Research, D4004) and sequenced directly (Lone Star Labs, Houston, TX, USA or Genewiz, Houston, TX, USA) with the primers used for amplification. (Figure 7(d)) oil immersion objective. GFP and chlorophyll were excited with a 488-nm argon laser. GFP emission was collected between 493 nm and 572 nm, and chlorophyll autofluorescence was collected between 620 nm and 719 nm. Each image in Figure  5 is an average of 2 exposures using a 47-µm pinhole, corresponding to a 0.8-µm optical section. Each image in Figure 7(d) is an average of 4 exposures using a 47-µm pinhole, corresponding to a 1.2-µm optical section.

Statistical analysis
One-way ANOVA with the Duncan post-hoc test was used to assess statistical significance (SPSS Statistics software, Version 24.0.0.0).