# Increases in activity of proteasome and papain-like cysteine protease in Arabidopsis autophagy mutants: back-up compensatory effect or cell-death promoting effect?

Marien Havé, Thierry Balliau, Betty Cottyn-Boitte, Emeline Dérond, Gwendal Cueff, Fabienne Soulay, Aurélia Lornac, Pavel Reichman, Nico Dissmeyer, Jean-Christophe Avice, et al.

► To cite this version:

Marien Havé, Thierry Balliau, Betty Cottyn-Boitte, Emeline Dérond, Gwendal Cueff, et al.. Increases in activity of proteasome and papain-like cysteine protease in Arabidopsis autophagy mutants: back-up compensatory effect or cell-death promoting effect?. Journal of Experimental Botany, 2018, 69 (6), pp.1369-1385. 10.1093/jxb/erx482. hal-02184349

HAL Id: hal-02184349

https://normandie-univ.hal.science/hal-02184349

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH PAPER

# Increases in activity of proteasome and papain-like cysteine protease in Arabidopsis autophagy mutants: back-up compensatory effect or cell-death promoting effect?

Marien Havé<fnref n="1" />, Thierry Balliau<fnref n="2" />, Betty Cottyn-Boitte<fnref n="1" />, Emeline Dérond<fnref n="1" />, Gwendal Cueff<fnref n="1" />, Fabienne Soulay<fnref n="1" />, Aurélia Lornac<fnref n="3" />, Pavel Reichman<fnref n="4" />, Nico Dissmeyer<fnref n="4" />, Jean-Christophe Avice<fnref n="3" />, Patrick Gallois<fnref n="5" />, Loïc Rajjou<fnref n="1" />, Michel Zivy<fnref n="2" /> and Céline Masclaux-Daubresse<fnref n="1" />,*

Received 11 September 2017; Editorial decision 14 December 2017; Accepted 14 December 2017

Editor: Peter Bozhkov, Swedish University of Agricultural Sciences, Sweden

## Abstract

Autophagy is essential for protein degradation, nutrient recycling, and nitrogen remobilization. Autophagy is induced during leaf ageing and in response to nitrogen starvation, and is known to play a fundamental role in nutrient recycling for remobilization and seed filling. Accordingly, ageing leaves of Arabidopsis autophagy mutants (atg) have been shown to over-accumulate proteins and peptides, possibly because of a reduced protein degradation capacity. Surprisingly, atg leaves also displayed higher protease activities. The work reported here aimed at identifying the nature of the proteases and protease activities that accumulated differentially (higher or lower) in the atg mutants. Protease identification was performed using shotgun LC-MS/MS proteome analyses and activity-based protein profiling (ABPP). The results showed that the chloroplast FTSH (FILAMENTATION TEMPERATURE SENSITIVE H) and DEG (DEGRADATION OF PERIPLASMIC PROTEINS) proteases and several extracellular serine proteases [subtilases (SBTs) and serine carboxypeptidase-like (SCPL) proteases] were less abundant in atg5 mutants. By contrast, proteasome-related proteins and cytosolic or vacuole cysteine proteases were more abundant in atg5 mutants. Rubisco degradation assays and ABPP showed that the activities of proteasome and papain-like cysteine protease were increased in atg5 mutants. Whether these proteases play a back-up role in nutrient recycling and remobilization in atg mutants or act to promote cell death is discussed in relation to their accumulation patterns in the atg5 mutant compared with the salicylic acid-depleted atg5/sid2 double-mutant, and in low nitrate compared with high nitrate conditions. Several of the proteins identified are indeed known as senescence- and stress-related proteases or as spontaneous cell-death triggering factors.

**Keywords:** AALP, RD21, CATHB3, SAG12, metacaspase, nitrogen remobilization, senescence.

<footnote n="1">¹ INRA-AgroParisTech, Institut Jean-Pierre Bourgin, UMR1318, ERL CNRS 3559, Saclay Plant Sciences, Versailles, France</footnote>

<footnote n="2">² UMR GQE- le Moulon, INRA, Université Paris-Sud, CNRS, AgroParisTech, Université Paris-Saclay, 91190 Gif-sur-Yvette, France</footnote>

<footnote n="3">³ UCBN, INRA, UMR INRA-UBCN 950 Ecophysiologie Végétale, Agronomie & Nutrition N.C.S., Université de Caen Normandie, F-14032 Caen, France</footnote>

<footnote n="4">⁴ Independent Junior Research Group on Protein Recognition and Degradation, Leibniz Institute of Plant Biochemistry (IPB), Weinberg 3, D-06120 Halle (Saale), Germany and Science Campus Halle – Plant-based Bioeconomy, Betty-Heimann-Strasse 3, D-06120 Halle (Saale), Germany</footnote>

<footnote n="5">⁵ School of Biological Sciences, Faculty of Biology, Medicine and Health, University of Manchester, Michael Smith Building, Oxford Road, Manchester, M139PT, UK</footnote>

<footnote>* Correspondence: celine.masclaux-daubresse@inra.fr</footnote>
## Introduction

Autophagy is a universal mechanism that facilitates the degradation of unwanted cell constituents in the lytic compartments of eukaryotic cells. Autophagy is essential for the recycling of cellular material and controls nitrogen remobilization and grain-filling in Arabidopsis (Guiboileau et al., 2012; Masclaux-Daubresse et al., 2017). The autophagy machinery consists of the formation of cytosolic double-membrane vesicles, termed autophagosomes, that engulf and sequester unwanted cytoplasmic constituents such as damaged organelles and protein aggregates (Liu and Bassham, 2012). Degradation does not occur directly inside the autophagosome as molecules destined for recycling are transported to the lytic vacuoles where proteases and hydrolases operate. Autophagosome formation involves the products of AUTOPHAGY (ATG) genes, previously discovered in yeast by the pioneering work of Professor Yoshinori Ohsumi (Nobel Prize in medicine or physiology 2016; Tsukada and Ohsumi, 1993). Subsequently, homologous genes in plants and animals have been found for almost all ATGs. Among the 50 ATG genes discovered in yeast, 18 are part of the autophagy core machinery (Yang and Bassham, 2015). These genes are absolutely essential to the formation of autophagosomes. In Arabidopsis, the core machinery genes are either single or belong to gene families. In particular, ATG5, involved in the ATG5-ATG12 conjugation system, is essential for the formation of the ATG8-PE (PE, phosphatidylethanolamine) conjugate that features on the autophagosome membrane and helps the expansion of the vesicle membrane (Masclaux-Daubresse et al., 2017). The atg5 Arabidopsis mutants have been intensively studied (Thompson et al., 2005; Yoshimoto et al., 2009; Lenz et al., 2011; Guiboileau et al., 2012, 2013; Masclaux-Daubresse et al., 2014; Floyd et al., 2015; Wang et al., 2015; Fahy et al., 2017). Like the other atg mutants studied so far (atg9, atg7, and atg2), the atg5 mutants present smaller rosettes, hypersensitivity to N and C starvation, and reduced yield. Using ¹⁵N-tracer experiments, Guiboileau et al. (2012) showed that the atg5 and atg9 mutants and the atg18a-RNAi line are strongly affected in N remobilization, and we have verified further that this is also the case for many other atg single-mutants (C. Masclaux-Daubresse, unpublished results). Several lines of evidence show that N and C leaf metabolisms are strongly affected in atg5, atg9, and atg18a-RNAi (Guiboileau et al., 2013; Masclaux-Daubresse et al., 2014). Mutants have lower starch, sucrose, and hexose contents in their leaves and higher amounts of amino acids (especially glutamate, aspartate, and methionine). They also display higher protein and RNA concentrations relative to wild-type plants. In order to try and explain such phenotypes, Guiboileau et al. (2013) measured the endo-, carboxy-, and amino-peptidase activities, and a series of western blots were also carried out in order to determine whether protein accumulation was selective. The results showed that several proteins and peptides were specifically accumulated in the leaves of atg mutants. The absence of transcriptional changes suggested that these proteins and peptides were specific cargoes. Surprisingly, protease activities were significantly increased

in the leaves of *atg* mutants (Guiboileau *et al.*, 2013), showing that the accumulation of the peptides and proteins over-abundant in the *atg* rosettes was not related to lower proteolysis capacity. In was then hypothesized that in the absence of autophagosome trafficking, vacuole proteases and their cytoplasmic substrates cannot co-localize, resulting in protein over-accumulation.

The aim of this study was to identify which protease activities are increased in autophagy mutants in order to compensate for the defects in autophagy-dependent protein recycling and to provide alternative remobilization pathways. Proteomics shotgun LC-MS/MS analyses were carried out in order to identify the proteases and proteins with increased or decreased accumulation in *atg5*. Because part of the phenotypes of autophagy mutants is related to the increase of salicylic acid (SA) in *atg* leaves, the *atg5/sid2* double-mutant was also analysed in order to discriminate SA-dependent and SA-independent changes. As a low-nitrate regime is known to promote nitrogen remobilization and autophagy activity, mutants were grown under both low- and high-nitrate conditions. Low-nitrate is indeed more favourable for the detection of autophagy-related phenotypes (Guiboileau *et al.*, 2012, 2013; Masclaux-Daubresse *et al.*, 2014). In addition to the identification of the proteases differentially accumulated in *atg* mutants, we monitored whether they were active by using specific probes (Morimoto and van der Hoorn, 2016), and performed western blots in order to evaluate the respective amounts of their pro- and active-forms. The results are discussed with regards to the known relationships that exist between protease activities, leaf senescence, N remobilization, and cell death. This work reveals new candidates that may participate in the management of N remobilization and possibly in the last steps of autophagy in Arabidopsis leaves.

## Material and methods

### Plant material and growth conditions

*Arabidopsis thaliana* (L.) Columbia wild-type, *atg5* (SALK_020601), *atg5/sid2*, and *sid2* mutants have been previously characterized by Yoshimoto *et al.* (2009) and Guiboileau *et al.* (2012, 2013). Plants were cultivated according to Guiboileau *et al.* (2013) under low- (LN; 2 mM) and high- (HN; 10 mM) nitrate conditions. Adequate growth conditions for the testing of nitrogen-limitation effects on Arabidopsis plants have previously been established by Loudet *et al.* (2003a, 2003b) and further improved by Lemaître *et al.* (2008), (see Supplementary Method S1 at JXB online for details). Whole rosettes were harvested at 60 d after sowing (DAS). Three independent plant replicates were sampled. For HN, each biological replicate contained four rosettes; for the LN, each biological replicate contained 24 rosettes. Harvests commenced 2 h into the 8-h photoperiod and were completed within 1 h. Samples were stored at –80 °C for subsequent use. Three batches of plants were cultured, providing samples from three independent experiments. The analyses detailed below were then performed on three biological replicates and
repeated on samples from 2-3 independent batches of plants. For validation of the anti-RD21A antibody, we used the homozygous *rd21a* SALK_090550C T-DNA insertion mutant (Wang *et al.*, 2008). For validation of the anti-CATHB3 antibody, the SALK_19630 *cathb3* homozygous T-DNA insertion mutant was used (Ge *et al.*, 2016).

## Shotgun proteomic analysis

Leaf total proteins were extracted, digested, and analysed by LC-MS/MS using the PAPSSO platform (INRA, Le Moulon, Gif sur Yvette; see Supplementary Method S2 for details). Protein identification was performed using X! Tandem Piledriver (version 2015.04.01; http://www.thegpm.org/TANDEM/index.html) by querying the MS/MS data against the TAIR10 protein library together with a custom contaminant database (trypsin, keratins). Identified proteins were filtered and grouped using X! Tandem Pipeline (3.4.1; pappso.inra.fr/bioinfo/xtandempipeline/) (Langella *et al.*, 2017) according to two criteria: (1) a minimum of two different peptides required with an *E*-value smaller than 0.01, and (2) a protein *E*-value (calculated as the product of unique peptide *E*-values) smaller than $10^{-5}$. The false discovery rates (FDRs) at the peptide and protein levels were 0.03% and 0.0%, respectively. Relative peptide quantification by peak-area integration on eXtracted ion chromatogram (XICs) was performed using the MassChroQ software (pappso.inra.fr/bioinfo/masschroq/) (Valot *et al.*, 2011). Relative protein abundance was calculated and defined as the sum of peptide intensities considering only (1) reproducible peptides, (2) specific peptides, and (3) correlated peptides belonging to the same protein (see Supplementary Method S2 for details). When the peptides of a protein were not present or not reproducibly observed in one or more conditions, spectral counting (SC) was used in place of XICs analysis.

## Activity-based protein profiling of papain-like cysteine proteases, proteasome and serine proteases

Papain-like cysteine proteases (PLCPs), serine proteases, and proteasome activity profiling were assayed using DCG-04 (van der Hoorn *et al.*, 2004), desthiobiotin-FP (DFP; Thermo Scientific) and MVB072 (Kolodziejek *et al.*, 2011), respectively. DCG-04 was synthesized as described by Greenbaum *et al.* (2000). The protein extract obtained by homogenizing 200 mg of frozen material with 450 µl of water (1.5% polyvinylpolypyrrolidone, w/v) was centrifuged (20 000 g, 15 min, 4 °C). Protein concentration was determined using the 2D Quant kit (SIGMA-ALDRICH). All labelling reactions were conducted in 100 µl volume containing 100 µg of proteins in the dark at room temperature. For PLCPs, proteins were labelled with 50 mM sodium acetate (NaAc), pH 5.5, 2 mM DTT, and 5 µM DCG-04 for 150 min. For serine proteases, proteins were labelled with 50 mM NaAc, pH 5.5, or 50 mM Hepes-KOH, pH 7.5, and 2 µM desthiobiotin-FP for 1 h. For the proteasome, proteins were labelled with 50 mM Hepes-KOH, pH 7.5, and 1 µM MVB072 for 1 h. For each genotype, an additional sample corresponding to a mixture of an equal

amount of soluble proteins from high-N and low-N was used as a control to check for probe specificity. Control samples were labelled with each probe as already described and were also subjected to competition assays consisting of a pre-treatment with 100 µM E-64, 50 µM DFP, or 50 µM epoxomicin for 30 min before adding DCG-04, desthiobiotin-FP, or MVB072, respectively. As a no-probe-control (NPC), control samples were incubated with an equal volume of probe solvent (DMSO). Labelling reactions were stopped by addition of 25 µl of 4× SDS-PAGE loading buffer and β-mercaptoethanol (355 mM) and heated at 90 °C for 10 min. For each sample, 10 µg of labelled proteins were separated by 12% or 15% SDS-PAGE. Resolved MVB072-labeled proteins were visualized by in-gel fluorescence scanning using a Typhoon 9400 scanner (GE Healthcare Life Science) with excitation and emission wavelengths at 532 and 580 nm, respectively. Resolved DCG-04- and desthiobiotin-FP-labelled proteins were transferred to a polyvinylidene fluoride membrane. Membranes were incubated overnight in 1% BSA (w/v) in PBS-T 0.1% and incubated with streptavidin-HRP (Ultrasensitive) for 1 h. The chemiluminescence signal was visualized using an ImageQuant LAS 4000 (GE Healthcare Life Science). Colloidal Blue-stained gels were used to monitor protein loading. Fluorescence, chemiluminescence, and Coomassie signals were quantified using the ImageQuant software (GE Healthcare Life Science).

## Affinity purification and identification of PLCPs

Only plants grown under low-N were analysed. For each genotype, soluble proteins from four independent rosette bulk samples were extracted as described above. An equal amount of proteins (750 µg) from first and second, and from third and fourth bulk samples were mixed to obtain two mixed protein solutions per genotype. Then, 1.5 mg of proteins were labelled with 50 mM NaAc, pH 5.5, 2 mM DTT, and 5 µM DCG-04 in 1 ml of reaction volume and incubated for 150 min at room temperature in the dark. Labelling reactions were stopped and the biotinylated proteins were purified as described by Chandrasekar *et al.* (2014), except that avidin beads were used (SIGMA-ALDRICH). The elution step was repeated once and the two eluates were mixed and separated by 12% SDS-PAGE. Gels were stained with silver nitrate and bands with a molecular mass between 25 and 40 kDa, corresponding to bands observed in the activity-based protein profiling (ABPP) experiment, were excised with a OneTouch GridCutter (2.0 × 7.0 mm; Gel Company). Each slice was placed into a 0.5-ml tube, digested, and analysed by LC-MS/MS as detailed in Supplementary Method S3. Proteomics data analyses were performed as explained in Supplementary Method S3. Identified proteins were filtered and grouped using X! Tandem Pipeline (3.4.2) according to the following criteria: (1) a minimum of two different peptides required with an *E*-value smaller than 0.05, and (2) a protein *E*-value smaller than $10^{-4}$.

## Western blots

Total proteins were extracted by homogenizing 100 mg of frozen material in 500 µl of buffer [7 M urea, 2 M thiourea,
30 mM Tris-HCl pH 8.8, 4% CHAPS (w/v), 0.2% triton X100 (w/v), 20 mM DTT, and 1× cOmplete protease inhibitor cocktail (Roche)]. Extracts were centrifuged and protein concentration was determined as described above. For each sample, 5 µg of proteins were separated on 12% SDS-PAGE. Resolved proteins were electroblotted (Trans-Blot Turbo transfer system; Bio-Rad). Polypeptide detections were performed using antibodies raised against the RD21A-specific peptide DELPESIDWRKKG, against SAG12 (Agrisera, Vännäs, Sweden), and against the N-terminal end of the CATHB3 mature protein (synthetic peptide LPKAFDARTAWPQC). Primary antibodies were diluted 1:1000 in 5% milk PBS-T 0.1% buffer for membrane incubation. Secondary HRP peroxidase antibodies (1:10000) were used for chemiluminescence detection. Signals were detected using an ImageQuant LAS 4000 (GE Healthcare Life Science). Coomassie Brilliant Blue-stained gels were used to verify that an equal amount of protein was loaded in each lane. RD21A and CATHB3 antibody specificities were tested using *rd21A* and *cathb3* mutants as negative controls (see Supplementary Fig. S1).

## Bioinformatics and statistical analyses

All statistical analyses were performed using the R software (https://www.r-project.org/). Data were analysed using two-way ANOVA with genotype and nutrition conditions as variable factors, followed by a Tukey's *post-hoc* test. For the relative quantification of proteins by the XICs method, statistical analyses were carried out on log₁₀-transformed protein abundance. Proteins corresponding to proteases with a Tukey's *P*-value ≤0.05 were considered to be significantly differentially accumulated. Fold-change (FC) ratios were determined as abundance ratios of *atg5* versus Col and of *atg5/sid2* versus *sid2*. No FC threshold was applied as even a small change in protease abundance can have a strong impact on the proteome due to post-translational regulation and the large substrate spectrum of proteases. In order to identify proteases among the proteins detected by LC-MS/MS, Gene Ontology and functional information for each protein were analysed using the following databases: MEROPS (https://www.ebi.ac.uk/merops/), UniProt (www.uniprot.org/), agriGO (bioinfo.cau.edu.cn/agriGO), and Mapman (mapman.gabipd.org/).

## Accession numbers

AT4G38220; AT2G27020; AT2G24200; AT4G20850; AT5G35590; AT4G31300; AT3G22110; AT2G05840; AT4G14800; AT3G60820; AT1G53750; AT3G05530; AT4G17510; AT5G05780; AT5G58290; AT5G10540; AT1G21720; AT1G56450; AT5G42790; AT4G01610; AT1G47128; AT1G53850; AT5G45890; AT5G51070; AT3G13235; AT1G79340; AT1G50380; AT3G51260; AT5G23540; AT4G38630; AT1G51710; AT5G66140; AT4G30910; AT4G39090; AT5G43060; AT5G10760; AT1G16470; AT1G79210; AT2G14260; AT5G23140; AT4G16190; AT5G26860; AT3G14290; AT3G20630;

AT2G41790; AT4G36760; AT5G36210; AT1G50250; AT5G42270; AT4G36195; AT1G09750; AT3G54400; AT3G02110; AT2G39850; AT1G52510; AT2G47940; AT2G35780; AT5G39830; AT3G27925; AT3G14067; AT3G19170; AT5G42390; AT3G24590; AT1G13270; AT5G05740; AT2G30950; AT4G18370; AT3G52500; AT5G42240; AT1G06430; AT5G23210; AT3G05350; AT5G65760; AT1G67700; AT3G61820; AT1G01300; AT4G34980; AT2G05920; AT3G18490; AT2G33530; AT4G21650; AT5G08260; AT1G09130; AT1G21750; AT1G47710; AT4G16500; AT517290; AT5G60360; AT5G50260; AT1G02305; AT4G35350.

## Results

### LC-MS/MS identification of protease proteins in autophagy mutants

Shotgun LC-MS/MS proteome analyses were carried out to identify the nature of the proteases specifically accumulated in *atg5* mutants. Total proteins were extracted from the autophagy-defective lines (*atg5* and *atg5/sid2*) and their respective controls (Col and *sid2*) grown for 60 d under high-N (HN) and low-N (LN) as previously described by Guiboileau *et al.* (2013). The *atg5/sid2* double-mutant was analysed in addition to *atg5* in order to discriminate against the side effects of salicylic acid (SA) (Guiboileau *et al.*, 2012). With regards to previous studies (Yoshimoto *et al.*, 2009; Guiboileau *et al.*, 2012, 2013; Masclaux-Daubresse *et al.*, 2014), we preferred here to consider *atg5/sid2* versus *sid2* rather than *atg5*/NahG versus NahG in order to discriminate against SA effects because of the lower SA content found in the *sid2* background relative to the NahG background (see Supplementary Fig. S2). The rosette phenotypes at harvest are shown in Fig. 1a, b. Under HN conditions, no leaf senescence phenotype was observed in any of the four genotypes (Fig. 1a). Under LN conditions, a slight senescence phenotype was observed on the oldest leaves of *atg5*; this phenotype was mitigated in *atg5/sid2*. In accordance with Guiboileau *et al.* (2013), we confirmed that the total protein concentration was higher in *atg5* relative to the Col wild-type under LN, and significantly higher in *atg5/sid2* compared to *sid2* under both LN and HN (Fig. 1c).

LC-MS/MS shotgun proteome analyses were then carried out on the total protein extracted from three independent biological repeats. XICs analysis permitted the relative quantification of proteins when peptides were detected in all the genotypes. When a peptide was undetectable in one of the two genotypes under comparison, individual peptide quantification by spectral count (SC) was performed. XICs and SC analyses, respectively, identified increased accumulation of 47 and decreased accumulation of 36 protease-related proteins in the autophagy-defective lines relative to control lines (Tables 1 and 2).

The predicted localization at subcellular level of the proteases and protease-related proteins that either increased or decreased in accumulation in the autophagy-deficient lines
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Autophagy-defective mutants accumulate proteins in their rosette leaves. Plants were cultivated under high- (a) or low- (b) nitrate conditions for 60 d. Total proteins (c) were extracted from rosette leaves using TCA/acetone and suspended in denaturing buffer. Values are means ±SD of three biological replicates. Significant differences from controls (i.e. *atg5* versus *Col*, and *atg5/sid2* versus *sid2*) are indicated: \*\*\*P<0.001 (ANOVA and Tukey's *post-hoc* test for multiple comparisons). Three independent plant cultures showed similar results: only one representative sample is presented here.</figcaption></figure>

are shown in Fig. 2a. Among the 47 protease-related proteins that increased in autophagy mutants, we found mainly proteasome subunits (22 proteins; Fig. 2b) and cytosol-localized proteases (28% of 26 proteases; Fig. 2a). Proteasome subunits were both core proteases (α and β sub-units: 50% and 23%, respectively; Fig. 2b) and regulatory particles (14% for both base and lid complexes). Vacuole-, plastid-, and extracellular-predicted proteases were also identified (17%, 16%, and 14%, respectively; Fig. 2a). Accumulated proteases were mostly cysteine (38%), serine (27%), and metal proteases (23%) (Fig. 2c).

The 36 proteases that had decreased accumulation were mainly plastid-predicted (44%) and extracellular-predicted proteases (42%), and they belong to the serine (55%), metallo- (28%) and, to a lesser extent, aspartate (17%) proteases. There was no cysteine proteases or proteasome sub-units among the decreased proteins.

Venn diagrams of the proteases and proteasome sub-units that were significantly differentially accumulated in autophagy-deficient lines are presented in Fig. 3. Among the proteases that were accumulated compared to the controls (Table 1; Fig. 3a), seven were more abundant in both *atg5* and *atg5/sid2*, and under both nitrate conditions (Class1-up), four proteases were more abundant in *atg5/sid2* under both LN and HN, and also in *atg5* under LN (Class2-up), and 12 were specifically over-accumulated under LN in both *atg5* and *atg5/sid2*, i.e. in a SA-independent manner (Class3-up). The majority of the Class1-up, 2-up, and 3-up proteases

were part of the proteasome machinery (Table 1), either threonine proteases or regulatory sub-units of the proteasome. Several proteasome-associated proteases such as LAP1, TOP2, and TPP2, or protein deubiquitinases such as UCH3 were also identified (Polge et al., 2009). In Class1-up, only AQI (AQUAPORIN INTERACTOR) was not related to the proteasome, and in Class3-up, three well-known cysteine proteases CATHB3 (CATHEPSIN B3), RD21A (RESPONSIVE TO DEHYDRATION 21A), and SAG12 (SENESCENCE ASSOCIATED GENE 12) were found. Class4-up included proteases significantly accumulated compared to the controls in *atg5* only (i.e. SA-dependent) and under both LN and HN. It included the chloroplast CLPD subunit (also referred to as SAG15 or ERD1), the Arabidopsis METACASPASE 4 (AMC4), the ubiquitin-family aspartate protease DDI1 (DNA-DAMAGE INDUCIBLE 1), which is induced under cadmium stress, and an unknown serine protease (At1g50380). The proteins of Class5-up were specifically over-accumulated in *atg5* under LN. Class5-up included two proteasome non-catalytic sub-units (PAD1 and PAD2), two proteasome regulatory proteins (RPN10 and RPN11), the deubiquitinase UBP6, the leucine aminopeptidase LAP3, the two drought resistance-related cysteine proteases RD19A and RD21B, and the AED1 (APOPLASTIC, EDS1-DEPENDENT 1) aspartate protease. Class6-up, which represents proteases over-accumulated in *atg5/sid2* under high N, included PIP (PROLINE IMINOPEPTIDASE) and the two proteasome subunits PAB1 and PAB2. Class7-up
<table><caption>**Table 1.** List of the proteases with significantly increased accumulation in *atg5* versus *Col*, and in *atg5/sid2* versus *sid2*</caption><thead><tr><th rowspan="3"></th><th rowspan="3">M<sup>(1)</sup></th><th rowspan="3">Accession (TAIR10)</th><th rowspan="3">Name<sup>(2)</sup></th><th rowspan="3">Catalytic class<sup>(3)</sup></th><th rowspan="3">Predicted subcellular localization<sup>(4)</sup></th><th colspan="4">Fold-change<sup>(5)</sup></th></tr><tr><th colspan="2">atg5</th><th colspan="2">atg5/sid2</th></tr><tr><th>HN</th><th>LN</th><th>HN</th><th>LN</th></tr></thead><tbody><tr><td rowspan="11">Class 1</td><td>XICs</td><td>AT4G38220</td><td>AQI</td><td>Met</td><td>ER</td><td><strong>1.63</strong></td><td><strong>1.77</strong></td><td><strong>1.84</strong></td><td><strong>2.21</strong></td></tr><tr><td>XICs</td><td>AT2G27020</td><td>PAG1</td><td>Thr</td><td>20S (CP α)</td><td><strong>1.41</strong></td><td><strong>1.79</strong></td><td><strong>1.64</strong></td><td><strong>1.84</strong></td></tr><tr><td>XICs</td><td>AT2G24200</td><td>LAP1</td><td>Met</td><td>PM,C</td><td><strong>1.44</strong></td><td><strong>1.55</strong></td><td><strong>1.42</strong></td><td><strong>1.41</strong></td></tr><tr><td>XICs</td><td>AT4G20850</td><td>TPP2</td><td>Ser</td><td>PI</td><td><strong>1.38</strong></td><td><strong>1.48</strong></td><td><strong>1.38</strong></td><td><strong>1.55</strong></td></tr><tr><td>XICs</td><td>AT5G35590</td><td>PAA1</td><td>Thr</td><td>20S (CP α)</td><td><strong>1.24</strong></td><td><strong>1.58</strong></td><td><strong>1.28</strong></td><td><strong>1.63</strong></td></tr><tr><td>XICs</td><td>AT4G31300</td><td>PBA1</td><td>Thr</td><td>20S (CP β)</td><td><strong>1.33</strong></td><td><strong>1.37</strong></td><td><strong>1.40</strong></td><td><strong>1.41</strong></td></tr><tr><td>XICs</td><td>AT3G22110</td><td>PAC1</td><td>Thr</td><td>20S (CP α)</td><td><strong>1.28</strong></td><td><strong>1.51</strong></td><td><strong>1.20</strong></td><td><strong>1.46</strong></td></tr><tr><td>XICs</td><td>AT2G05840</td><td>PAA2</td><td>Thr</td><td>20S (CP α)</td><td>1.28</td><td><strong>1.73</strong></td><td><strong>1.43</strong></td><td><strong>1.72</strong></td></tr><tr><td>XICs</td><td>AT4G14800</td><td>PBD2</td><td>Thr</td><td>20S (CP β)</td><td>1.19</td><td><strong>1.80</strong></td><td><strong>1.43</strong></td><td><strong>1.65</strong></td></tr><tr><td>XICs</td><td>AT3G60820</td><td>PBF1</td><td>Thr</td><td>20S (CP β)</td><td>1.15</td><td><strong>1.49</strong></td><td><strong>1.23</strong></td><td><strong>1.52</strong></td></tr><tr><td rowspan="11">Class 2</td><td>XICs</td><td>AT1G53750</td><td>RPT1A</td><td>non-catalytic</td><td>19S (RP base)</td><td>1.25</td><td><strong>1.42</strong></td><td><strong>1.29</strong></td><td><strong>1.47</strong></td></tr><tr><td>XICs</td><td>AT3G05530</td><td>RPT5A</td><td>non-catalytic</td><td>19S (RP base)</td><td>1.31</td><td><strong>2.29</strong></td><td>1.39</td><td><strong>2.33</strong></td></tr><tr><td>XICs</td><td>AT4G17510</td><td>UCH3</td><td>Cys</td><td>C</td><td>1.55</td><td><strong>2.27</strong></td><td>1.82</td><td><strong>2.09</strong></td></tr><tr><td>XICs</td><td>AT5G05780</td><td>RPN8A</td><td>non-catalytic</td><td>19S (RP lid)</td><td>1.18</td><td><strong>2.44</strong></td><td>1.31</td><td><strong>1.88</strong></td></tr><tr><td>XICs</td><td>AT5G58290</td><td>RPT3</td><td>non-catalytic</td><td>19S (RP base)</td><td>1.08</td><td><strong>2.17</strong></td><td>1.38</td><td><strong>2.13</strong></td></tr><tr><td>XICs</td><td>AT5G10540</td><td>TOP2</td><td>Met</td><td>C,PI</td><td>1.05</td><td><strong>2.26</strong></td><td>1.44</td><td><strong>2.03</strong></td></tr><tr><td>XICs</td><td>AT1G21720</td><td>PBC1</td><td>Thr</td><td>20S (CP β)</td><td>1.33</td><td><strong>1.96</strong></td><td>1.41</td><td><strong>2.01</strong></td></tr><tr><td>XICs</td><td>AT1G56450</td><td>PBG1</td><td>Thr</td><td>20S (CP β)</td><td>1.18</td><td><strong>1.86</strong></td><td>1.21</td><td><strong>1.78</strong></td></tr><tr><td>XICs</td><td>AT5G42790</td><td>PAF1</td><td>Thr</td><td>20S (CP α)</td><td>1.41</td><td><strong>1.80</strong></td><td>1.14</td><td><strong>1.83</strong></td></tr><tr><td>XICs</td><td>AT4G01610</td><td>CATHB3</td><td>Cys</td><td>E,V</td><td>1.13</td><td><strong>1.84</strong></td><td>0.85</td><td><strong>1.79</strong></td></tr><tr><td>XICs</td><td>AT1G47128</td><td>RD21A</td><td>Cys</td><td>E,V</td><td>0.91</td><td><strong>1.46</strong></td><td>0.74</td><td><strong>1.57</strong></td></tr><tr><td>SC</td><td>AT1G53850</td><td>PAE1</td><td>Thr</td><td>20S (CP α)</td><td>1.75</td><td><strong>1.69</strong></td><td>1.50</td><td><strong>1.67</strong></td></tr><tr><td>SC</td><td>AT5G45890</td><td>SAG12</td><td>Cys</td><td>V</td><td>1.50 (+)</td><td><strong>22.50 (+)</strong></td><td>1.00 (+)</td><td><strong>24.50 (+)</strong></td></tr><tr><td rowspan="11">Class 4</td><td>XICs</td><td>AT5G51070</td><td>CLPD</td><td>Ser</td><td>PI</td><td><strong>1.49 (+)</strong></td><td><strong>2.37 (+)</strong></td><td>1.13</td><td>1.40 (+)</td></tr><tr><td>XICs</td><td>AT3G13235</td><td>DDI1</td><td>Asp</td><td>C</td><td><strong>1.53</strong></td><td><strong>1.61</strong></td><td>1.21</td><td>1.21</td></tr><tr><td>XICs</td><td>AT1G79340</td><td>AMC4</td><td>Cys</td><td>PM,C</td><td><strong>1.41</strong></td><td><strong>1.51</strong></td><td>1.28</td><td>1.36</td></tr><tr><td>XICs</td><td>AT1G50380</td><td>—</td><td>Ser</td><td>C</td><td><strong>1.30</strong></td><td><strong>1.34</strong></td><td>1.13</td><td>1.08</td></tr><tr><td rowspan="11">Class 5</td><td>XICs</td><td>AT3G51260</td><td>PAD1</td><td>Thr</td><td>20S (CP α)</td><td>1.58</td><td><strong>2.53</strong></td><td>1.72</td><td>1.79</td></tr><tr><td>XICs</td><td>AT5G23540</td><td>RPN11</td><td>non-catalytic</td><td>19S (RP lid)</td><td>1.10</td><td><strong>1.34</strong></td><td>1.10</td><td>1.26</td></tr><tr><td>XICs</td><td>AT4G38630</td><td>RPN10</td><td>non-catalytic</td><td>19S (RP lid)</td><td>1.19</td><td><strong>1.31</strong></td><td>1.28</td><td>1.28</td></tr><tr><td>XICs</td><td>AT1G51710</td><td>UBP6</td><td>Cys</td><td>C</td><td>1.13</td><td><strong>1.22</strong></td><td>1.06</td><td>1.05</td></tr><tr><td>SC</td><td>AT5G66140</td><td>PAD2</td><td>Thr</td><td>20S (CP α)</td><td>1.25</td><td><strong>5.00</strong></td><td>2.33</td><td>1.67</td></tr><tr><td>SC</td><td>AT4G30910</td><td>LAP3</td><td>Met</td><td>PI</td><td>0.89</td><td><strong>1.35</strong></td><td>1.20</td><td>1.45</td></tr><tr><td>XICs</td><td>AT4G39090</td><td>RD19A</td><td>Cys</td><td>V,N</td><td>1.15</td><td><strong>1.74</strong></td><td>0.96</td><td>1.56</td></tr><tr><td>XICs</td><td>AT5G43060</td><td>RD21B</td><td>Cys</td><td>E,V</td><td>0.99 (+)</td><td><strong>1.49 (+)</strong></td><td>0.59</td><td>1.36</td></tr><tr><td>SC</td><td>AT5G10760</td><td>AED1</td><td>Asp</td><td>E</td><td>1.80 (+)</td><td><strong>3.33 (+)</strong></td><td>0.67</td><td>1.50 (+)</td></tr><tr><td rowspan="11">Class 6</td><td>SC</td><td>AT1G16470</td><td>PAB1</td><td>Thr</td><td>20S (CP α)</td><td>1.50</td><td>1.40</td><td><strong>2.22</strong></td><td>1.35</td></tr><tr><td>SC</td><td>AT1G79210</td><td>PAB2</td><td>Thr</td><td>20S (CP α)</td><td>1.55</td><td>1.11</td><td><strong>2.00</strong></td><td>1.21</td></tr><tr><td>XICs</td><td>AT2G14260</td><td>PIP</td><td>Ser</td><td>C,PI</td><td>1.18</td><td>1.23</td><td><strong>1.48</strong></td><td>1.24</td></tr><tr><td>XICs</td><td>AT5G23140</td><td>CLPP2</td><td>Ser</td><td>M</td><td>0.93</td><td>2.42</td><td>0.86</td><td><strong>2.77</strong></td></tr><tr><td>XICs</td><td>AT4G16190</td><td>RD19C</td><td>Cys</td><td>E,V</td><td>1.05</td><td>1.46</td><td>1.18</td><td><strong>2.08</strong></td></tr><tr><td>XICs</td><td>AT5G26860</td><td>LON1</td><td>Ser</td><td>M</td><td>1.21</td><td>1.31</td><td>1.02</td><td><strong>1.40</strong></td></tr><tr><td>SC</td><td>AT3G14290</td><td>PAE2</td><td>Thr</td><td>20S (CP α)</td><td>1.40</td><td>1.33</td><td>1.42</td><td><strong>1.63</strong></td></tr><tr><td>SC</td><td>AT3G20630</td><td>UBP14</td><td>Cys</td><td>C</td><td><strong>1.43</strong></td><td>1.40</td><td><strong>1.42</strong></td><td>1.31</td></tr><tr><td>SC</td><td>AT2G41790</td><td>PXM16</td><td>Met</td><td>Pe</td><td><strong>1.54</strong></td><td>1.29</td><td>1.43</td><td>1.11</td></tr><tr><td>SC</td><td>AT4G36760</td><td>APP1</td><td>Met</td><td>PM,C</td><td><strong>1.51</strong></td><td><strong>1.72</strong></td><td>1.27</td><td><strong>1.62</strong></td></tr><tr><td>SC</td><td>AT5G36210</td><td>—</td><td>Ser</td><td>PI</td><td><strong>1.36</strong></td><td>1.24</td><td><strong>1.57</strong></td><td><strong>1.31</strong></td></tr></tbody></table>

Plants were grown under high- (HN) or low- (LN) nitrate conditions. (1) XICs or Spectral counting (SC) methods. (2) Protein names according to UniProt and TAIR. (3) Catalytic classes according to MEROPS: asp, aspartic proteases; cys, cysteine proteases; met, metallo-proteases; ser, serine proteases; thr, threonine proteases; non-catalytic, non-catalytic proteasome regulatory sub-units. (4) Predicted subcellular localizations according to SUBA3 (SUBAcon), Bio-Analytic Ressource for Plant Biology (Cell eFP Viewer) and Marshall et al. (2015). 20S (CP α), proteasome 20S core protease α; 20S (CP β), proteasome 20S core protease β; 19S (RP base), proteasome 19S regulatory particle base; 19S (RP lid), proteasome 19S regulatory particle lid; C, cytosol; ER, endoplasmic reticulum; E, extracellular; M, mitochondrion; Pe, peroxisome; PM, plasma membrane; PI, plastid; V, vacuole. (5) Fold-change ratios were calculated by dividing protein abundance in *atg5* and *atg5/sid2* by protein abundance in *Col* and in *sid2*, respectively. Entries in bold represent significantly increased accumulations, all other entries are non-significant changes (n = 3; *P*<0.05; ANOVA and Tukey's *post-hoc* test for multiple comparisons). (+) significant increase in gene expression (Log₂ fold-change of mRNA levels >1.2 with *P*-value <10⁻⁵, rank product and FDR estimation).
<table><caption>**Table 2.** List of the proteases with significantly decreased accumulation in *atg5* versus *Col*, and *atg5/sid2* versus *sid2*</caption><thead><tr><th rowspan="3"></th><th rowspan="3">M<sup>(1)</sup></th><th rowspan="3">Accession (TAIR10)</th><th rowspan="3">Name<sup>(2)</sup></th><th rowspan="3">Catalytic class<sup>(3)</sup></th><th rowspan="3">Predicted subcellular localization<sup>(4)</sup></th><th colspan="4">Fold-change<sup>(5)</sup></th></tr><tr><th colspan="2">atg5</th><th colspan="2">atg5/sid2</th></tr><tr><th>HN</th><th>LN</th><th>HN</th><th>LN</th></tr></thead><tbody><tr><td rowspan="3">Class 1</td><td>XICs</td><td>AT1G50250</td><td>FTSH1</td><td>met</td><td>PI</td><td><strong>0.76</strong></td><td><strong>0.49</strong></td><td><strong>0.67</strong></td><td><strong>0.42</strong></td></tr><tr><td>XICs</td><td>AT5G42270</td><td>FTSH5</td><td>met</td><td>PI</td><td><strong>0.71</strong></td><td><strong>0.51</strong></td><td><strong>0.67</strong></td><td><strong>0.48</strong></td></tr><tr><td>XICs</td><td>AT4G36195</td><td>-</td><td>ser</td><td>V,PM</td><td><strong>0.66</strong></td><td><strong>0.61</strong></td><td><strong>0.77</strong></td><td><strong>0.52</strong></td></tr><tr><td rowspan="2">Class 2</td><td>XICs</td><td>AT1G09750</td><td>AED3</td><td>asp</td><td>E</td><td>0.78</td><td><strong>0.63</strong></td><td><strong>0.65</strong></td><td><strong>0.61</strong></td></tr><tr><td>XICs</td><td>AT3G54400</td><td>-</td><td>asp</td><td>E,PI</td><td>0.77</td><td><strong>0.62</strong></td><td><strong>0.72</strong></td><td><strong>0.59</strong></td></tr><tr><td rowspan="8">Class 3</td><td>XICs</td><td>AT3G02110</td><td>SCPL25</td><td>ser</td><td>E</td><td><strong>0.60</strong></td><td><strong>0.46</strong></td><td>0.70</td><td><strong>0.44</strong></td></tr><tr><td>XICs</td><td>AT2G39850</td><td>SBT4.1</td><td>ser</td><td>E</td><td><strong>0.60</strong></td><td><strong>0.48</strong></td><td>0.78</td><td><strong>0.59</strong></td></tr><tr><td>XICs</td><td>AT1G52510</td><td>-</td><td>Ser</td><td>PI</td><td><strong>0.69</strong></td><td><strong>0.52</strong></td><td>0.81</td><td><strong>0.54</strong></td></tr><tr><td>XICs</td><td>AT2G47940</td><td>DEGP2</td><td>Ser</td><td>PI</td><td><strong>0.71</strong></td><td><strong>0.56</strong></td><td>0.81</td><td><strong>0.58</strong></td></tr><tr><td>XICs</td><td>AT2G35780</td><td>SCPL26</td><td>Ser</td><td>V,E</td><td><strong>0.81</strong></td><td><strong>0.58</strong></td><td>0.74</td><td><strong>0.46</strong></td></tr><tr><td>XICs</td><td>AT5G39830</td><td>DEGP8</td><td>Ser</td><td>PI</td><td><strong>0.66</strong></td><td><strong>0.57</strong></td><td>0.75</td><td><strong>0.67</strong></td></tr><tr><td>XICs</td><td>AT3G27925</td><td>DEGP1</td><td>Ser</td><td>PI</td><td><strong>0.76</strong></td><td><strong>0.68</strong></td><td>0.77</td><td><strong>0.68</strong></td></tr><tr><td>XICs</td><td>AT3G14067</td><td>SASP</td><td>Ser</td><td>E</td><td><strong>0.78</strong></td><td><strong>0.73</strong></td><td>0.88</td><td><strong>0.74</strong></td></tr><tr><td rowspan="10">Class 4</td><td>XICs</td><td>AT3G19170</td><td>PREP1</td><td>met</td><td>PI</td><td><strong>0.81</strong></td><td><strong>0.74</strong></td><td>0.84</td><td><strong>0.74</strong></td></tr><tr><td>XICs</td><td>AT5G42390</td><td>SPP</td><td>met</td><td>PI</td><td>0.81</td><td><strong>0.38</strong></td><td>0.75</td><td><strong>0.43</strong></td></tr><tr><td>XICs</td><td>AT3G24590</td><td>PLSP1</td><td>Ser</td><td>PI</td><td>0.69</td><td><strong>0.36</strong></td><td>0.75</td><td><strong>0.50</strong></td></tr><tr><td>XICs</td><td>AT1G13270</td><td>MAP1C</td><td>met</td><td>PI</td><td>0.70</td><td><strong>0.56</strong></td><td>0.71</td><td><strong>0.48</strong></td></tr><tr><td>XICs</td><td>AT5G05740</td><td>EGY2</td><td>met</td><td>PI</td><td>0.78</td><td><strong>0.55</strong></td><td>0.85</td><td><strong>0.56</strong></td></tr><tr><td>XICs</td><td>AT2G30950</td><td>FTSH2</td><td>met</td><td>PI</td><td>0.75</td><td><strong>0.57</strong></td><td>0.68</td><td><strong>0.54</strong></td></tr><tr><td>XICs</td><td>AT4G18370</td><td>DEGP5</td><td>Ser</td><td>PI</td><td>0.66</td><td><strong>0.54</strong></td><td>0.83</td><td><strong>0.65</strong></td></tr><tr><td>XICs</td><td>AT3G52500</td><td>-</td><td>Asp</td><td>E</td><td>0.82</td><td><strong>0.65</strong></td><td>0.77</td><td><strong>0.61</strong></td></tr><tr><td>XICs</td><td>AT5G42240</td><td>SCPL42</td><td>Ser</td><td>E</td><td>0.86</td><td><strong>0.66</strong></td><td>0.87</td><td><strong>0.62</strong></td></tr><tr><td rowspan="5">Class 5</td><td>XICs</td><td>AT1G06430</td><td>FTSH8</td><td>met</td><td>PI</td><td>0.77</td><td><strong>0.71</strong></td><td>0.74</td><td><strong>0.58</strong></td></tr><tr><td>XICs</td><td>AT5G23210</td><td>SCPL34</td><td>Ser</td><td>E</td><td>0.86</td><td><strong>0.65</strong></td><td>0.87</td><td><strong>0.68</strong></td></tr><tr><td>XICs</td><td>AT3G05350</td><td>-</td><td>met</td><td>PI</td><td>0.83</td><td><strong>0.66</strong></td><td>0.88</td><td><strong>0.69</strong></td></tr><tr><td>XICs</td><td>AT5G65760</td><td>-</td><td>Ser</td><td>PI,V</td><td>0.86</td><td><strong>0.70</strong></td><td>0.93</td><td><strong>0.68</strong></td></tr><tr><td>XICs</td><td>AT1G67700</td><td>HHL1</td><td>met</td><td>PI</td><td>0.76</td><td><strong>0.39</strong></td><td>0.74</td><td>0.53</td></tr><tr><td rowspan="5">Class 6</td><td>XICs</td><td>AT3G61820</td><td>-</td><td>Asp</td><td>E</td><td>0.84</td><td><strong>0.61</strong></td><td>0.72</td><td>0.65</td></tr><tr><td>XICs</td><td>AT1G01300</td><td>APF2</td><td>Asp</td><td>E,PM</td><td>0.86</td><td><strong>0.64</strong></td><td>0.73</td><td>0.73</td></tr><tr><td>XICs</td><td>AT4G34980</td><td>SLP2</td><td>Ser</td><td>E</td><td>0.81</td><td><strong>0.71</strong></td><td>0.87</td><td>1.16</td></tr><tr><td>XICs</td><td>AT2G05920</td><td>SBT1.8</td><td>Ser</td><td>E,PM</td><td>0.82</td><td><strong>0.74</strong></td><td>0.71</td><td>0.77</td></tr><tr><td>XICs</td><td>AT3G18490</td><td>ASPG1</td><td>Asp</td><td>E</td><td>0.93</td><td>0.81</td><td>0.85</td><td><strong>0.60</strong></td></tr><tr><td rowspan="3">Class 7</td><td>XICs</td><td>AT2G33530</td><td>SCPI46</td><td>Ser</td><td>E</td><td>1.22</td><td>0.97</td><td>0.89</td><td><strong>0.63</strong></td></tr><tr><td>XICs</td><td>AT4G21650</td><td>SBT3.13</td><td>Ser</td><td>E</td><td>0.71</td><td>0.55 (-)</td><td>0.69</td><td><strong>0.65</strong></td></tr><tr><td>XICs</td><td>AT5G08260</td><td>SCPL35</td><td>Ser</td><td>E</td><td>0.90</td><td>0.83</td><td><strong>0.66</strong></td><td><strong>0.65</strong></td></tr><tr><td>Class 8</td><td>XICs</td><td>AT1G09130</td><td>CLPR3</td><td>Ser</td><td>PI</td><td>0.78</td><td>1.00</td><td><strong>0.41</strong></td><td>0.91</td></tr></tbody></table>

Plants were grown under high- (HN) or low- (LN) nitrate conditions. (1) Method for proteome analysis was XICs. (2) Protein names according to UniProt and TAIR. (3) Catalytic classes according to MEROPS: asp, aspartic proteases; met, metallo-proteases; ser, serine proteases.

(4) Predicted subcellular localizations according to SUBA3 (SUBAcon) and Bio-Analytic Ressource for Plant Biology (Cell eFP Viewer). E, extracellular; PM, plasma membrane; PI, plastid; V, vacuole. (5) Fold-change ratios were calculated by dividing protein abundance in *atg5* and in *atg5/sid2* by protein abundance in *Col* and in *sid2*, respectively. Entries in bold represent significantly decreased accumulations, all other entries are non-significant changes (*n* = 3; *P*<0.05; ANOVA and Tukey's post-hoc test for multiple comparisons). (-) significant decrease in gene expression (Log₂ fold-change of mRNA levels < -1.2 with *P*-value <10⁻⁵, rank product and FDR estimation).

grouped together four proteases over-accumulated in *atg5/sid2* under LN. Two of them are essential for mitochondria maintenance (LON1 and CLPP2), whilst the others are a proteasome subunit (PAE2) and a cysteine protease RD19C (RESPONSIVE TO DEHYDRATION 19C). Classes 8-up, 9-up, 10-up, and 11-up included UBP14 (UBIQUITIN-SPECIFIC PROTEASE 14; accumulated in *atg5* and *atg5/sid2* under HN only), PXM16 (metallo-protease of the M16 family; only accumulated in *atg5* under HN), APP1 (accumulated in all conditions and *atg* genetic backgrounds except in *atg5/sid2* under HN), and the unknown serine protease (At5g36210; accumulated compared to the controls in *atg5*

and *atg5/sid2* under HN, and in *atg5/sid2* under LN). The ALEURAIN-LIKE PROTEASE, AALP (AT5G60360), detected by the shotgun LC-MS/MS proteome analysis was found to be slightly, but not significantly, increased in the *atg* lines under LN (see Supplemenatry Fig. S3).

The proteases that had decreased accumulation were also classified. The two FTSH1 and FTSH5 (FILAMENTATION TEMPERATURE SENSITIVE H 1 and 5) chloroplast proteases and an unknown vacuole serine protease (At4g36195) were significantly less abundant in both *atg5* and *atg5/sid2*, and under both LN and HN (Class1-down; Table 2). Class2-down contained two extracellular-predicted proteases [AED3
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Predicted localization and catalytic classes of the proteases and proteasome sub-units in the autophagy-defective lines with significantly increased (up) or decreased (down) accumulation. Plants were grown under low- (LN) or high- (HN) nitrate conditions. Proteases and proteasome sub-units with significantly increased or decreased accumulation in *atg5* versus Col and/or *atg5/sid2* versus *sid2* under LN or HN, or both LN and HN are presented. The numbers of proteins considered for the percentage calculation are indicated in brackets. Significant differences were determined by ANOVA and Tukey's *post-hoc* test for multiple comparisons (*P*<0.05, *n*=3). (a) Predicted subcellular localization of the increased and decreased proteases in the autophagy-defective lines using SUBAcon (SUBA3 database) and Cell eFP Viewer (Bio-Analytic Ressource for Plant Biology). Only catalytic proteasome sub-units (CPβ) are taken into account here. (b) A total of 22 20S core-protease (CP) and 19S regulatory-particle (RP) sub-units were increased in the autophagy-defective lines: The diagram represents their proportions. (c) Catalytic classes of the increased and decreased proteases in the autophagy-defective lines identified using the MEROPS database. Only the CP β sub-units of proteasome were considered here as active proteases.</figcaption></figure>

(APOPLASTIC, EDS1-DEPENDENT 3) and the product of *AT3G54400*. Class3-down (down in *atg5* under both LN and HN, and in *atg5/sid2* under LN) principally included several serine proteases, almost all predicted in the extracellular space (such as SASP, SENESCENCE-ASSOCIATED SUBTILISIN PROTEASE), or in the plastids (such as the DEGRADATION OF PERIPLASMIC PROTEINS DEGP1, DEGP2, and DEGP8 proteases, and the PREP1 zinc metallo-protease). Class4-down clustered together proteases that were decreased in *atg5* and *atg5/sid2* only under LN. They were almost all plastid- (FTSH2, FTSH8, DEGP5) or extracellular-predicted, and they mostly belong to the serine and metallo-protease families. Class5-down grouped together proteases repressed in *atg5* under LN, and mainly included extracellular proteases, but also contained the HHL1 (HYPERSENSITIVE TO HIGH LIGHT 1) metallo-protease that is involved in the protection of photosystem II (Jin et al., 2014). Class6-down contained proteases differentially

repressed in the *atg5/sid2* mutant under LN. They are almost all extracellular-predicted and serine proteases except for ASPG1 (ASPARTIC PROTEASE IN GUARD CELL 1), which is extracellular-predicted, but also localizes in the endoplasmic reticulum in the guard cells (Yao *et al.*, 2012). Class7-down and Class8-down both contained only one member: the extracellular SCPL35 (decreased accumulation in *atg5/sid2* under LN and HN) and the chloroplast CLPR3 (decreased accumulation in *atg5/sid2* under HN), respectively.

From these results, it clearly appears that the protease-related proteins accumulated in *atg5* and/or *atg5/sid2* mostly belonged to the proteasome machinery and to the cytosolic or vacuole cysteine protease families. The proteases with decreased accumulation were mainly ones with chloroplast- or extracellular-predicted localizations and belonged mainly to the serine- and metallo-protease families. It should be noted that the decrease of chloroplast-predicted proteases in *atg5* was independent of SA (observed in both the Col and
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Venn diagrams of the proteases and proteasome sub-units with significantly increased (up) and decreased (down) accumulation in the autophagy-defective lines. Comparisons of *atg5* versus Col (black lines) and *atg5/sid2* versus *sid2* (grey lines) were performed on plants grown under low- (LN) or high- (HN) nitrate conditions. The number of proteins is indicated in brackets. Significant differences were determined by ANOVA and Tukey's *post-hoc* test for multiple comparisons (*P*<0.05, *n*=3).</figcaption></figure>

*sid2* backgrounds) and independent of N (observed under HN as well as under LN). It is thus likely that the decrease in these chloroplast-predicted proteases was not correlated to the slight senescence phenotype observed in *atg5* under low-nitrate conditions (Fig. 1b), but constitutes a more robust phenotype that was shared by the two autophagy-defective lines studied, independently of the nitrate conditions.

In order to determine whether changes in protease contents were linked to a modification of the transcription of their genes or to modifications in their protein turnovers, gene expression levels were evaluated using the transcriptional data obtained by Masclaux-Daubresse et al. (2014), (Supplementary Dataset S1). Only the CLPD, RD21B, AED1, and SAG12 genes were up-regulated in *atg5* and/or *atg5/sid2* (Table 1), and only SBT3.13 was repressed in *atg5* under low nitrate (Table 2). Thus, the changes in the protease protein contents reported in Tables 1 and 2 were more likely due to modifications in their turnover rates.

*Proteasome and papain-like cysteine protease activities are higher in the atg5 and atg5/sid2 lines, but not serine protease activity*

It is well known that protease activities can be regulated at several post-transcriptional levels. In addition, it has been shown in the case of the proteasome, for example, that autophagy could play a role in proteaphagy (Marshall *et al.*, 2015). Thus, in order to determine whether the accumulation or depletion of proteases detected in the *atg5* mutant lines caused modifications in protease activities, we explored changes in specific protease activities.

The degradation of the Rubisco large sub-unit (RBCL) at acidic pH was monitored using specific protease inhibitors according to the procedure described by Poret *et al.* (2016). We observed that RBCL degradation was higher in the two

*atg5* lines relative to controls (see Supplementary Fig. S4a) and that degradation was inhibited when the E-64 cysteine-protease inhibitor or the aprotinine serine protease inhibitor were added before incubation of the extract (Fig. S4b).

ABPP was performed in order to detect changes in the contents of active proteasome sub-units, active cysteine proteases, and active serine proteases (Morimoto and van der Hoorn, 2016). For the proteasome, the fluorescent probe MVB072 was added to the protein samples in order to derivatize the active proteasome protease sub-units. Fluorescent derivatized sub-units were then detected after SDS-PAGE separation (Fig. 4). MVB072 competition assays were carried out to estimate the probe's specificity toward the proteasome sub-units (see Supplementary Fig. S5). Significant increases in the proteasome-active β1, β2, and β5 sub-units (encoded by PBA1, PBB1/PBB2 and PBE1/PBE2 genes respectively) were detected in *atg5* under HN, and in both *atg5* and *atg5/sid2* under LN. By contrast with PBA1 (proteasome-active β1subunit), no PBB1 or PBE2 subunits were identified in the shotgun analyses. These increases suggested that the proteasome activity was increased in autophagy-defective lines and especially under LN.

In order to detect active papain-like cysteine proteases (PLCPs), DCG-04 was added to the protein samples before SDS-PAGE separation. The biotin-tag that linked to the active PLCP was detected on membranes by western blot analysis using streptavidin-HRP chemiluminescence (Fig. 5). Five bands were detected on the membranes and their chemiluminescence signal was clearly higher in *atg5* and *atg5/sid2* relative to Col and *sid2*, respectively, and especially under LN (Fig. 5a). The controls performed to demonstrate the specificity of the synthesized DCG-04 toward cysteine proteases are presented in Supplementary Fig. S6.

The active serine proteases were detected using desthio-biotin-FP, and at both acidic and neutral pH as serine
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Proteasome activity is higher in autophagy-defective lines. Autophagy-defective lines (*atg5* and *atg5/sid2*) and their respective controls (Col and *sid2*) were cultivated under high- (HN) or low-(LN) nitrate conditions for 60 d. Soluble proteins were extracted from rosette leaves using water and labelled with MVB072 at pH 7.5 to detect proteasome activity. Proteasome activity (a) was determined by detecting the MVB072-specific fluorescence after separating 10 µg of proteins by SDS-PAGE. The β1, β2, and β5 catalytic sub-units are indicated by arrows. (b) The same SDS-PAGE gel was stained with Coomassie Brilliant Blue (CBB) to check for equal protein input. The Rubisco large sub-unit (RBCL) is indicated by an arrow. The fluorescence intensity (c) representative of the proteasome activity was quantified using Image Quant. Proteasome activity was measured on six biological repeats from two independent cultures. Means ±SD are shown (n=6) and significant differences between the autophagy-defective lines and the controls (i.e. *atg5* versus Col, and *atg5/sid2* versus *sid2*) are indicated: **P**<0.01, ***P**<0.001 (ANOVA and Tukey's *post-hoc* test for multiple comparisons). The gels presented in this figure are representative of six biological replicates.</figcaption></figure>

proteases are localized in both vacuoles and chloroplasts (see Supplementary Figs. S7 and S8). Surprisingly, no differences between the *atg5* lines and control lines were detected. This suggests that the global serine protease activity was not different in the autophagy-defective lines compared to the control lines (Supplementary Fig. S8). Inconsistent with the data obtained using the aprotinine serine protease inhibitor on RBCL degradation at pH 5.5, this result suggests that the serine protease activities monitored using desthiobiotin-FP were not the same as those that degraded RBCL in our previous assay.

*Identification of active papain-like cysteine proteases in autophagy mutants.*

Given that LC-MS/MS analyses showed that cysteine pro-
tease proteins accumulated in the autophagy mutants, and

<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Activities of papain-like cysteine proteases (PLCPs) are higher in autophagy-defective lines grown under low-nitrate conditions. Autophagy-defective lines (*atg5* and *atg5/sid2*) and their respective controls (Col and *sid2*) were cultivated under high- (HN) or low- (LN) nitrate conditions for 60 d. (a) Active PLCPs labelled with DCG-04 at pH 5.5 (arrows) were detected using streptavidin-HRP chemiluminescence. (b) The same total amount of input proteins was loaded on each lane after incubation, as shown by the Coomassie Brilliant Blue (CBB)-stained gel. The Rubisco large sub-unit (RBCL) is indicated by an arrow. (c) The chemiluminescence signal intensity representative of the total PLCP activities was quantified using Image Quant. Activities of PLCPs were measured on three biological repeats. Means and ±SD are shown (n=3) and significant differences between the autophagy-defective lines and the controls (i.e. *atg5* versus Col, and *sid2/atg5* versus *sid2*) are shown: ***P**<0.001 (ANOVA and Tukey's *post-hoc* test for multiple comparisons). The gels presented in this figure are representative of three biological replicates. Determination of PLCP activities was repeated on samples from two independent cultures, showing similar results.</figcaption></figure>

that DCG-04 revealed several active PLCP bands on mem-
branes, a pull-down of active PLCPs was carried out using
DCG-04. Co-precipitated PLCPs were separated on SDS-
PAGE gels and streptavidin-HRP reactive bands were ana-
lysed using shotgun LC-MS/MS. This assay was only carried
out on extracts from plants grown under LN conditions, as
no differences had been detected on gels between genotypes
when grown under HN (Fig. 5). In this experiment, spec-
tral counting was used to determine the relative amounts
of active PLCPs in the different genetic backgrounds (see
Supplementary Table S1). The results showed that all the over-
accumulated PLCPs previously identified in the autophagy-
defective lines were active. The largest differences detected
between the *atg5* and control lines after DCG-04 pull-down
were for RD21A and SAG12 (Fig. 6b). Active CATHB3,
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Identification of the active papain-like cysteine proteases (PLCPs) accumulated in autophagy-defective lines. Autophagy defective lines (*atg5* and *atg5/sid2*) and their respective controls (Col and *sid2*) were cultivated under low- (LN) nitrate conditions for 60 d. Soluble proteins from two biological replicates were extracted and labelled with DCG-04 at pH 5.5. Purified DCG-04-labelled proteins after pull-down for each genotype were separated on 12% SDS-PAGE gels and stained with silver nitrate (a); two replicates, rep1 and rep2, are shown. Proteins of molecular size between 40 and 25 kDa were excised from the gel (five bands) and analysed by LC-MS/MS. Ten active PLCPs were identified in all the genotypes and the two biological replicates. The total spectral counts obtained for each PLCPs are presented in (b) and (c). Similar results were obtained with the two replicates. Spectral counts in all the different excised bands presented in (a) are given in Supplementary Table S1.</figcaption></figure>

AALP, RD21B, RD19A, and RD19C were also more abundant in *atg5* and *atg5/sid2* (Fig. 6c) Additional cysteine proteases were found that had not been identified by the shotgun LC-MS/MS whole-proteome analysis. They were CEP1 (CYSTEINE ENDOPEPTIDASE 1; At5g50260), which is involved in tapetum cell death and pathogen defence (Zhang *et al.*, 2014), CATHB2 (CATHEPSIN-LIKE B2 CASPASE-LIKE; At1g02305), which is senescence-induced (Ge *et al.*, 2016), and XCP1 (XYLEM CYSTEINE PEPTIDASE 1; At4g35350), which is a xylem cysteine peptidase (Avci *et al.*, 2008). It should be noted that *CEP1* was over-expressed at the transcriptional level in *atg5* and *atg5/sid2*, and under both HN and LN.

### Western blot identification of mature and immature forms of papain-like cysteine proteases

Antibodies raised against RD21A, SAG12, and CATHB3 were used in order to verify the accumulation of their processed forms and to monitor their immature forms in the *atg5* lines under LN and also under HN (Fig. 7). The RD21A intermediate form (i) was detected in both *atg5* and control lines, and under both LN and HN. Its level was very high in all the lines, although it was higher in the *atg5* lines than in controls. The mature forms (m1 and m2; Fig. 7) were much

less abundant than the intermediate forms, but still more abundant in the *atg5* lines than in controls. The m2 form was notably higher in the *atg5* lines under LN than under HN, and in *atg5/sid2* than in *atg5* under LN. For SAG12, the immature signal was weak, but higher in the *atg5* lines grown under LN. The mature form (m) was only observed in the *atg5* lines under low N and was much stronger than the immature form (i), in contrast with RD21A. All the CATHB3 mature and immature forms were more abundant in the *atg5* lines relative to the control lines, Col and *sid2*. The intermediate form with a processed C-terminal domain (ΔC) and the mature form (M) of CATHB3 were especially higher in the *atg5* lines than in the controls under low-nitrate conditions. Therefore, western blots confirmed that the active forms of RD21A, SAG12, and CATHB3 were more abundant in the autophagy-defective lines.

### Inhibitor proteins of cysteine proteases are differentially accumulated in autophagy-defective lines

Several protease inhibitors are small proteins that can form stable complexes with target proteases and block or prevent access to the active site. Among the long list of protease inhibitors available in the MEROPS database, three were found to be differentially accumulated in autophagy-defective
<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Identification of the mature and immature forms of papain-like cysteine proteases (PLCPs) using western blots. Autophagy-defective lines (*atg5* and *atg5/sid2*) and their respective controls (Col and *sid2*) were cultivated under high- (HN) or low- (LN) nitrate conditions for 60 d. Total proteins were extracted from rosette leaves using denaturing buffer. Proteins (5 µg) were separated on 12% SDS-PAGE gels and detected on protein blots using RD21A-, CATHB3-, and SAG12-specific antibodies. Arrows indicate the different PLCP forms: intermediate (i) and mature (m, m1, m2) for SAG12 and RD21A; immature pro-enzyme (i), intermediate form with C-terminal domain processed (ΔC), mature form (M), and short mature form (m) of CATHB3. Equal amounts of proteins were loaded in each lane, as shown by the Coomassie Brilliant Blue (CBB)-stained blot, on which the Rubisco large sub-unit (RBCL) is indicated by an arrow.</figcaption></figure>

lines compared to the controls (Table 3). PDI5 and Serpin-1, that can inhibit RD21A, had increased accumulation in *atg5* under both LN and HN conditions. PDI5 also accumulated in *atg5/sid2*. The extracellular CYS4 cysteine protease inhibitor was less abundant in *atg5/sid2* under LN.

Discussion

The aim of this study was to identify protease activities up-regulated in autophagy mutants in order to further

determine their roles in nitrogen remobilization for seed filling (Guiboileau *et al.*, 2012). Shotgun LC-MS/MS and activity-based protein profiling (ABPP; Morimoto and van der Hoorn, 2016) analyses showed that, while protease genes were not transcriptionally induced or repressed in the autophagy mutants (*atg5* and *atg5/sid2*) relative to their controls (Col and *sid2*, respectively), many proteases were differentially abundant in the *atg* mutants (Fig. 8).

The proteases depleted in the *atg* mutants were mainly chloroplast- and extracellular-predicted proteases. Several chloroplast FTSH and DEG proteases were less abundant in the *atg5* lines (Kato and Sakamoto, 2010; van Wijk, 2015). FTSH-1, -2, -5, and -8 are the four major isomers of the chloroplast FTSH complexes involved in thylakoid maintenance (Yu *et al.*, 2004). Their decrease in the *atg* lines, especially under LN, suggested that thylakoid maintenance was affected. The DEGs ATP-independent serine-type proteases may be involved in the degradation of thylakoid lumen proteins and of PSII (Chassin *et al.*, 2002) and DEGP2 may be responsible for the endoproteolytic cleavage of the protein D1 (Haussühl *et al.*, 2001). Hence, together with the decrease of FTSH proteins, the decrease of DEGP-1, -2, -5, and -8 protein abundance suggests that thylakoid and photosystem maintenance are affected in *atg* lines. The number of chloroplast proteases with decreased accumulation in the autophagy-defective lines was slightly higher under nitrate starvation, but was not dependent on the leaf-senescence phenotype or on SA synthesis. The lower abundance of chloroplast proteases in the autophagy-defective lines suggests that chloroplasts could be less abundant or less active in the *atg5* lines than in the wild-type, thus explaining the lower carbon, starch, and sugar contents measured in *atg5*, *atg9*, and *atg18a*-RNAi by Guiboileau *et al.* (2013) and Masclaux-Daubresse *et al.* (2014).

In addition to chloroplast proteases, many extracellular-predicted subtilisin-like proteases (subtilases; SBTs) and extracellular-predicted serine carboxypeptidase-like proteases (SCPLs) were decreased in the autophagy-defective lines. The SBTs have mainly been described in plant-pathogen or pest-defence responses (Figueiredo *et al.*, 2014). The SCPLs are more likely involved in secondary metabolism and could also have a role in plant defence (Liu *et al.*, 2008; Mugford and Osbourn, 2010). The current understanding of the function of SBTs remains very limited, although a role for some of them in plant immune priming has been proposed (Figueiredo *et al.*, 2014). It is difficult to explain why the *atg5* and *atg5/sid2* lines had reduced amounts of extracellular SBTs and SCPLs. Dysfunction in exocytosis in the *atg* mutants could be one explanation, although a relationship between exocytosis and autophagy remains to be clearly demonstrated in plants (Wang *et al.*, 2010; Lin *et al.*, 2015; Münz, 2017). Nevertheless, the lower extracellular content of SBTs in *atg5* might explain the higher susceptibility of the autophagy mutants to necrotrophic pathogens, as reported by several studies (Lai *et al.*, 2011; Lenz *et al.*, 2011). Martinez *et al.* (2015) recently reported that the extracellular senescence-induced SASP subtilisin protease was involved in branching; the authors hypothesized that SASP could be
<table><caption>**Table 3.** List of protease inhibitors differentially accumulated in *atg5* versus *Col*, and *atg5/sid2* versus *sid2*</caption><thead><tr><th rowspan="3">Accession (TAIR10)</th><th rowspan="3">Name<sup>(1)</sup></th><th rowspan="3">Description<sup>(2)</sup></th><th rowspan="3">Predicted subcellular localization<sup>(3)</sup></th><th colspan="4">Fold-change<sup>(4)</sup></th></tr><tr><th colspan="2">atg5</th><th colspan="2">atg5/sid2</th></tr><tr><th>HN</th><th>LN</th><th>HN</th><th>LN</th></tr></thead><tbody><tr><td>AT1G21750</td><td>PDIL1-1/PDI5</td><td>protein disulfide isomerase</td><td>ER</td><td><strong>1.77</strong></td><td><strong>2.45</strong></td><td><strong>1.39</strong></td><td><strong>1.95</strong></td></tr><tr><td>AT1G47710</td><td>SERPIN-1</td><td>ser-protease inhibitor</td><td>C</td><td><strong>1.23</strong></td><td><strong>1.23</strong></td><td>1.11</td><td>1.13</td></tr><tr><td>AT4G16500</td><td>CYS4</td><td>cys-protease inhibitor</td><td>E</td><td>0.92</td><td>0.77</td><td>0.77</td><td><strong>0.55</strong></td></tr></tbody></table>

Plants were grown under high- (HN) or low- (LN) nitrate conditions. (1) Protein names according to UniProt and TAIR. (2) Description (TAIR10). (3) Predicted subcellular localizations according to SUBA3 (SUBAcon) and Bio-Analytic Ressource for Plant Biology (Cell eFP Viewer). C, cytosol; ER, endoplasmic reticulum; E, extracellular. (4) Fold-change ratios were calculated by dividing protein abundance in *atg5* and *atg5/sid2* by protein abundance in *Col* and *sid2*, respectively. Entries in bold cells represent significantly increased or decreased accumulations, all other entries are non-significant changes (n = 3; *P*<0.05; ANOVA and Tukey's *post-hoc* test for multiple comparisons).

<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Schematic representation of the different proteases and protease-related proteins differentially accumulated in autophagy mutants and of their predicted cellular localization. The proteases with increased accumulation in *atg* are shown in red while those with decreased accumulation are shown in blue.</figcaption></figure>

involved in the processing of extracellular apoplastic signals repressing inflorescence branching and delaying leaf senescence. Although we found decreased accumulation of SASP in the *atg5* lines, no branching phenotype could be observed, but senescence was possibly enhanced.

Among the proteases with increased accumulation in the autophagy-defective lines, it was noticeable that many were proteasome sub-units and that several proteins were related

to the ubiquitin proteasome system (UPS), such as LAP1, TOP2, TPP2, UCH3, and UBP14. The UPS and autophagy are two major pathways that degrade most of the cellular proteins of eukaryotic cells. It is considered that while the proteasome is responsible for the turnover of short-lived proteins, autophagy is more specifically dedicated to the degradation of long-lived or aggregated proteins (Lilienbaum, 2013, and references therein). It appears that the UPS and autophagy are tightly controlled and co-ordinated. Several studies have suggested that autophagy may act as a back-up system, thereby assisting degradation in case the UPS is overloaded (Nedelsky et al., 2008). It seems that most of the mis-folded soluble proteins that are preferentially degraded by the UPS in cases of endoplasmic reticulum (ER) stress can be degraded through the autophagy pathway if the proteasome capacity is exceeded (Johnston et al., 1998). Several studies have also suggested that the proteasome itself could be degraded by autophagy (Waite et al., 2016). In plants, the accumulation of several proteasome sub-units in autophagy-defective mutants has been reported. Marshall et al., (2015) showed that RPN10 could act as a selective autophagy receptor that targets ubiquitinated and inactive 26S proteasomes and binds to ATG8 for autophagy-mediated degradation. We also found that several proteasome sub-units were more abundant in *atg5*, especially under LN, including RPN10 (Table 1). This suggests that proteophagy was increased under LN and was compromised in *atg5*. It thus seems that low-N conditions that trigger N remobilization (Masclaux-Daubresse and Chardon, 2011) can induce selective proteophagy (Waite et al., 2016; Zientara-Rytter and Sirko, 2016).

With regards to the proteasome, our ABPP assay showed that active β1, β2, and β5 catalytic subunits of the proteasome were more abundant in *atg5* under LN and HN and also in the *atg5/sid2* defective line under LN. This suggested that the activity of the proteasome was globally higher in the autophagy-defective lines. It has been shown that ER stress is increased in plant autophagy mutants (Munch et al., 2014; Yang et al., 2016). A role for the proteasome in the endoplasmic reticulum-associated degradation of mis-folded proteins under ER stress has been demonstrated in mammals (Lee et al., 2001). We can therefore suspect that an increase in proteasome activity in *atg5* and *atg5/sid2* could be linked to ER stress.
In addition to the accumulation of UPS proteins, increases of several cysteine, serine, and metallo-proteases were observed in the *atg5* and *atg5/sid2* lines. Activity measurements, performed using protease inhibitors, indicated that cysteine protease activities were significantly different between the *atg5* and *atg5/sid2* lines and the Col and *sid2* controls. As they are suspected to play a fundamental role in N remobilization during senescence and in response to nitrate shortage, we then focused on the papain-like cysteine proteases (PLCPs; Avice and Etienne, 2014; Diaz-Mendoza et al., 2016; Pružinská et al., 2017). Pull-down using DCG-04 showed that active PLCPs were more abundant in the *atg5* and *atg5/sid2* lines (Table 1). Active SAG12, RD21A, CATHB3, and AALP were strongly enriched in *atg5* and *atg5/sid2*, as well as RD19A, RD19C, CEP1, CATHB2, and XCP1, although to a lower extent. With the exception of SAG12 and CEP1, PLCP genes were not differentially expressed in *atg5* versus Col or *atg5/sid2* versus *sid2*. This absence of a significant change at the transcriptional level suggests that the accumulation of PLCPs in *atg5* and/or *atg5/sid2* was mainly due to higher translation or post-translational maintenance. Most of the PLCPs identified were predicted in the vacuole, except for RD21s that have been also found in ER-bodies (Hayashi et al., 2001), SAG12 that was detected in senescence-associated vacuoles (Otegui et al., 2005), and CEP1 that seems only to be located in the ER (Zhang et al., 2014). Many of the PLCPs, such as RD21s, RD19s, and CATHB3, transit through the ER before being delivered into the vacuole lumen where they play their protease functions. Whether autophagy could play a role in the trafficking of these proteases from the ER to the vacuole, as was suggested for seed-storage proteins (Michaeli et al., 2014), is unknown. The western blots performed using the RD21A, CATHB3, and SAG12 antibodies revealed that both mature and immature forms were accumulated the *atg5* and *atg5/sid2* lines, thus suggesting that there was no defect in protease maturation or trafficking in the autophagy mutants.

All the PLCPs identified in this study are known proteases. They have been described in several reports but their regulation, maturation, and physiological roles remain largely enigmatic (van der Hoorn, 2013; Rogers, 2013; Iglesias-Fernández et al., 2014; Zhang et al., 2014; Pathirana et al., 2017). Several are known to be transcriptionally induced during leaf senescence, such as SAG12, CATHB3, CATHB2, AALP, CEP1, RD21A, RD21B, RD19A, and RD19C (Supplementary Table S2; Noh and Amasino, 1999; Havé et al., 2017; Pružinská et al., 2017). However, as none of them except SAG12 and CEP1 were transcriptionally induced in the autophagy-defective lines relative to the controls, the possibility that these PLCPs were more abundant in *atg5* and *atg5/sid2* due to early senescence is weak. All these PLCPs were found to be more abundant in both *atg5* and *atg5/sid2* under LN but not under HN, and hence it is highly possible that they are involved in N remobilisation processes (Masclaux-Daubresse and Chardon, 2011), and are induced under low N and especially when autophagy is impaired. These proteases are thus good candidates as enzymes involved in the N-recycling pathways that lead to the remobilization of organic nitrogen forms from the source leaves to the sink leaves and to the seeds (Desclos et al., 2009).

Their potential roles as back-up proteolytic systems to compensate for weaknesses in autophagy, or as the last step of an autophagy process and cargo degradation inside the vacuole, need verification. Pružinská et al. (2017) have recently shown that only knockout lines lacking AALP present delayed leaf senescence under dark stress, in contrast with the *rd21a*, *rd21b*, *sag12*, and *cathb3* mutants. It therefore seems that AALP may be the best candidate for playing a role in N remobilization, together with, or in parallel to, autophagy.

It is well known that autophagy mutants exhibit spontaneous regions of cell death on their leaves with ageing (Masclaux-Daubresse et al., 2014). The role of the PLCPs identified here in this phenotype of the autophagy mutants can be questioned (Lampl et al., 2013; Iglesias-Fernández et al., 2014; Cai et al., 2018; Zhang et al., 2014). Indeed, several reports suggest that AMC4, CATHB3, PBA1, and RD21s might play a role in the control of cell death (Watanabe and Lam, 2011; Cai et al., 2018), either as cell death-promoting or longevity factors. Interestingly, we found that AMC4 accumulated in *atg5* but not in *atg5/sid2* whatever the nitrate conditions, while RD21A, CATHB3, and SAG12 accumulated in both *atg5* and *atg5/sid2* but only under LN. These discrepancies in accumulation patterns possibly distinguish the different functions of these proteases. Proteases promoting death should be independent of nitrate but SA-dependent, such as AMC4, while the N-remobilization proteases should be SA-independent but nitrate-dependent, such as RD21A, CATHB3, and SAG12. Hypotheses such as these need further investigations with mutants in order to be confirmed.

Changes in protein inhibitors of proteases were also identified. Increases in PDIL1/PDI5 and in Serpin-1 proteins suggested additional controls exerted on the RD21 proteases in the *atg5* and *atg5/sid2* lines to moderate their activity (Lampl et al., 2013; Rustgi et al., 2017). The decrease of CYS4 suggests that its unknown cysteine protease partners are repressed in the autophagy-defective lines (Je et al., 2014).

In conclusion, we identified several proteases that are differentially accumulated in autophagy-defective lines compared to their controls, and that could explain the phenotype of autophagy mutants such as the spontaneous regions of cell death on the leaves and the lower carbon and sugar contents due to reduced chloroplast maintenance. Many of the PLCPs identified could play a compensatory role for protein degradation and nutrient recycling, especially under low nitrate nutrition. Further investigations and biochemical studies are required using protease and autophagy mutants to clarify their contributions.

## Supplementary data

Supplementary data are available at JXB online.

Methods S1–S3. Phenotyping display and growth conditions; details of the shotgun proteomics analysis; and details of the shotgun LC-MS/MS analysis of pulled-down PLCPs.

Fig. S1. α-mRD21A and cathepsin B3-N Term-1610N antibodies are specific for Arabidopsis RD21A and cathepsin B3, respectively.
Fig. S2. Salicylic acid concentrations in *atg5-1/sid2* and *atg5-2/NahG* grown under low-nitrate conditions.

Fig. S3. ALEU protease amounts tend to increase in autophagy-defective mutants.

Fig. S4.: Rubisco degradation in *atg* defective lines and controls.

Fig. S5. MVB072 competition assay to estimate the probe's specificity toward proteasome sub-units.

Fig. S6. DCG-04 competition assay to estimate the probe's specificity toward PLCPs.

Fig. S7. Desthiobiotin-FP competition assay to estimate the probe's specificity toward serine proteases.

Fig. S8. Serine protease activities are stable in autophagy mutants.

Table S1. Details of spectral counts for PLCP peptides in Col, *sid2*, *atg5*, and *atg5/sid2*.

Table S2.: Senescence up- and down-effects on the transcription of protease and protease-inhibitor genes.

Dataset S1. Transcriptomic data obtained by Masclaux-Daubresse et al. (2014).

## Acknowledgements

This work was mainly supported by the financial support of ANR-12-ADAPT-0010-0 to MH and AL. The work of PR was supported by a Ph.D. scholarship from the German Academic Exchange Service (DAAD) and by the grant DI 1794/3-1 of the German Research Foundation (DFG) to ND. The authors are grateful to Renier van der Hoorn (Oxford University, Oxford, UK) for providing MVB072, to Marine Poret (University of Caen, France) for help, and to David Lawrence for English editing.

## Author contributions

MH performed all the experiments with technical help from ED. GC, FS, and TB ran the LC-MS/MS analyses. MH and TB collected and analysed the LC-MS/MS data. BC-B synthetized the DCG-04 probe. PG provided the CATHB3 antibody and performed controls on *cathb3* mutants to certify antibody specificities. ND and PR provided the anti-RD21A antibody and performed controls on the *rd21a* mutant and recombinant mRD21A protein. MZ and LR participated in the design of the biochemical experiments and in discussions. AL and J-CA designed and performed the Rubisco degradation assays. CM-D designed and supervised the research. MH and CM-D interpreted the results and wrote the manuscript.

## References

**Avci U, Petzold HE, Ismail IO, Beers EP, Haigler CH.** 2008. Cysteine proteases XCP1 and XCP2 aid micro-autolysis within the intact central vacuole during xylogenesis in Arabidopsis roots. The Plant Journal **56**, 303–315.

**Avice JC, Etienne P.** 2014. Leaf senescence and nitrogen remobilization efficiency in oilseed rape (*Brassica napus* L.). Journal of Experimental Botany **65**, 3813–3824.

**Cai YM, Yu J, Ge Y, Mironov A, Gallois P.** 2018. Two proteases with caspase-3-like activity, cathepsin B and proteasome, antagonistically control ER-stress-induced programmed cell death in Arabidopsis. New Phytologist. In press. doi:10.1111/nph.14676.

**Chandrasekar B, Colby T, Emon AEK, Jiang JB, Hong TN, Villamor JG, Harzen A, Overkleeft HS, van der Hoorn RAL.** 2014. Broad-range glycosidase activity profiling. Molecular & Cellular Proteomics **13**, 2787–2800.

**Chassin Y, Kapri-Pardes E, Sinvany G, Arad T, Adam Z.** 2002. Expression and characterization of the thylakoid lumen protease DegP1 from Arabidopsis. Plant Physiology **13**, 857–864.

**Desclos M, Etienne P, Coquet L, Jouenne T, Bonnefoy J, Segura R, Reze S, Ourry A, Avice JC.** 2009. A combined 15N tracing/proteomics study in *Brassica napus* reveals the chronology of proteomics events associated with N remobilisation during leaf senescence induced by nitrate limitation or starvation. Proteomics **9**, 3580–3608.

**Diaz-Mendoza M, Velasco-Arroyo B, Santamaria ME, González-Melendi P, Martinez M, Diaz I.** 2016. Plant senescence and proteolysis: two processes with one destiny. Genetics and Molecular Biology **39**, 329–338.

**Fahy D, Sanad MN, Duscha K, et al.** 2017. Impact of salt stress, cell death, and autophagy on peroxisomes: quantitative and morphological analyses using small fluorescent probe N-BODIPY. Scientific Reports **7**, 39069.

**Figueiredo A, Monteiro F, Sebastiana M.** 2014. Subtilisin-like proteases in plant-pathogen recognition and immune priming: a perspective. Frontiers in Plant Science **5**, 739.

**Floyd BE, Morriss SC, MacIntosh GC, Bassham DC.** 2015. Evidence for autophagy-dependent pathways of rRNA turnover in Arabidopsis. Autophagy **11**, 2199–2212.

**Ge Y, Cai YM, Bonneau L, Rotari V, Danon A, McKenzie EA, McLellan H, Mach L, Gallois P.** 2016. Inhibition of cathepsin B by caspase-3 inhibitors blocks programmed cell death in Arabidopsis. Cell Death and Differentiation **23**, 1493–1501.

**Greenbaum D, Medzihradszky KF, Burlingame A, Bogyo M.** 2000. Epoxide electrophiles as activity-dependent cysteine protease profiling and discovery tools. Chemistry & Biology **7**, 569–581.

**Guiboileau A, Avila-Ospina L, Yoshimoto K, Soulay F, Azzopardi M, Marmagne A, Lothier J, Masclaux-Daubresse C.** 2013. Physiological and metabolic consequences of autophagy deficiency for the management of nitrogen and protein resources in Arabidopsis leaves depending on nitrate availability. New Phytologist **199**, 683–694.

**Guiboileau A, Yoshimoto K, Soulay F, Bataillé MP, Avice JC, Masclaux-Daubresse C.** 2012. Autophagy machinery controls nitrogen remobilization at the whole-plant level under both limiting and ample nitrate conditions in Arabidopsis. New Phytologist **194**, 732–740.

**Haussühl K, Andersson B, Adamska I.** 2001. A chloroplast DegP2 protease performs the primary cleavage of the photodamaged D1 protein in plant photosystem II. Embo Journal **20**, 713–722.

**Havé M, Marmagne A, Chardon F, Masclaux-Daubresse C.** 2017. Nitrogen remobilization during leaf senescence: lessons from Arabidopsis to crops. Journal of Experimental Botany **68**, 2513–2529.

**Hayashi Y, Yamada K, Shimada T, Matsushima R, Nishizawa NK, Nishimura M, Hara-Nishimura I.** 2001. A proteinase-storing body that prepares for cell death or stresses in the epidermal cells of Arabidopsis. Plant & Cell Physiology **42**, 894–899.

**Iglesias-Fernández R, Wozny D, Iriondo-de Hond M, Oñate-Sánchez L, Carbonero P, Barrero-Sicilia C.** 2014. The *AtCathB3* gene, encoding a cathepsin B-like protease, is expressed during germination of *Arabidopsis thaliana* and transcriptionally repressed by the basic leucine zipper protein GBF1. Journal of Experimental Botany **65**, 2009–2021.

**Je J, Song C, Hwang JE, Chung WS, Lim CO.** 2014. DREB2C acts as a transcriptional activator of the thermo tolerance-related *phytocystatin 4* (*AtCYS4*) gene. Transgenic Research **23**, 109–123.

**Jin H, Liu B, Luo L, et al.** 2014. HYPERSENSITIVE TO HIGH LIGHT1 interacts with LOW QUANTUM YIELD OF PHOTOSYSTEM I1 and functions in protection of photosystem II from photodamage in Arabidopsis. The Plant Cell **26**, 1213–1229.

**Johnston JA, Ward CL, Kopito RR.** 1998. Aggresomes: a cellular response to misfolded proteins. The Journal of Cell Biology **143**, 1883–1898.

**Kato Y, Sakamoto W.** 2010. New insights into the types and function of proteases in plastids. International Review of Cell and Molecular Biology, **280**, 185–218.

**Kolodziejek I, Misas-Villamil JC, Kaschani F, et al.** 2011. Proteasome activity imaging and profiling characterizes bacterial effector syringolin A. Plant Physiology **155**, 477–489.

**Lai Z, Wang F, Zheng Z, Fan B, Chen Z.** 2011. A critical role of autophagy in plant resistance to necrotrophic fungal pathogens. The Plant Journal **66**, 953–968.
Langella O, Valot B, Balliau T, Blein-Nicolas M, Bonhomme L, Zivy M. 2017. xltandempeline: a tool to manage sequence redundancy for protein inference and phosphosite identification. Journal of Proteome Research **16**, 494–503.

Lampl N, Alkan N, Davydov O, Fluhr R. 2013. Set-point control of RD21 protease activity by AtSerpin1 controls cell death in Arabidopsis. The Plant Journal **74**, 498–510.

Lee RJ, Hart CL, Liu CW, McCracken AA, DeMartino G, Thomas PJ, Romisch K, Brodsky JL. 2001. The 26S proteasome is sufficient to retro-translocate and degrade a soluble ERAD substrate. Molecular Biology of the Cell **12**, 132A–132A.

Lemaître T, Gaufichon L, Boutet-Mercey S, Christ A, Masclaux-Daubresse C. 2008. Enzymatic and metabolic diagnostic of nitrogen deficiency in *Arabidopsis thaliana* Wassileskija accession. Plant & Cell Physiology **49**, 1056–1065.

Lenz HD, Haller E, Melzer E, et al. 2011. Autophagy differentially controls plant basal immunity to biotrophic and necrotrophic pathogens. The Plant Journal **66**, 818–830.

Lilienbaum A. 2013. Relationship between the proteasomal system and autophagy. International Journal of Biochemistry and Molecular Biology **4**, 1–26.

Lin Y, Ding Y, Wang J, et al. 2015. Exocyst-positive organelles and autophagosomes are distinct organelles in plants. Plant Physiology **169**, 1917–1932.

Liu H, Wang X, Zhang H, Yang Y, Ge X, Song F. 2008. A rice serine carboxypeptidase-like gene OsBISCPL1 is involved in regulation of defense responses against biotic and oxidative stress. Gene **420**, 57–65.

Liu Y, Bassham DC. 2012. Autophagy: pathways for self-eating in plant cells. Annual Review of Plant Biology **63**, 215–237.

Loudet O, Chaillou S, Krapp A, Daniel-Vedele F. 2003a. Quantitative trait loci analysis of water and anion contents in interaction with nitrogen availability in *Arabidopsis thaliana*. Genetics **163**, 711–722.

Loudet O, Chaillou S, Merigout P, Talbotec J, Daniel-Vedele F. 2003b. Quantitative trait loci analysis of nitrogen use efficiency in Arabidopsis. Plant Physiology **131**, 345–358.

Marshall RS, Li F, Gemperline DC, Book AJ, Vierstra RD. 2015. Autophagic degradation of the 26S proteasome is mediated by the dual ATG8/ubiquitin receptor RPN10 in Arabidopsis. Molecular Cell **58**, 1053–1066.

Martinez DE, Borniego ML, Battchikova N, Aro EM, Tyystjärvi E, Guiamét JJ. 2015. SASP, a senescence-associated subtilisin protease, is involved in reproductive development and determination of silique number in Arabidopsis. Journal of Experimental Botany **66**, 161–174.

Masclaux-Daubresse C, Chardon F. 2011. Exploring nitrogen remobilization for seed filling using natural variation in *Arabidopsis thaliana*. Journal of Experimental Botany **62**, 2131–2142.

Masclaux-Daubresse C, Chen Q, Havé M. 2017. Regulation of nutrient recycling via autophagy. Current Opinion in Plant Biology **39**, 8–17.

Masclaux-Daubresse C, Clément G, Anne P, Routaboul JM, Guiboileau A, Soulay F, Shirasu K, Yoshimoto K. 2014. Stitching together the multiple dimensions of autophagy using metabolomics and transcriptomics reveals impacts on metabolism, development, and plant responses to the environment in Arabidopsis. The Plant Cell **26**, 1857–1877.

Michaeli S, Avin-Wittenberg T, Galili G. 2014. Involvement of autophagy in the direct ER to vacuole protein trafficking route in plants. Frontiers in Plant Science **5**, 134.

Morimoto K, van der Hoorn RAL. 2016. The increasing impact of activity-based protein profiling in plant science. Plant & Cell Physiology **57**, 446–461.

Mugford ST, Osbourn A. 2010. Evolution of serine carboxypeptidase-like acyltransferases in the monocots. Plant Signaling & Behavior **5**, 193–195.

Munch D, Rodriguez E, Bressendorff S, Park OK, Hofius D, Petersen M. 2014. Autophagy deficiency leads to accumulation of ubiquitinated proteins, ER stress, and cell death in Arabidopsis. Autophagy **10**, 1579–1587.

Münz C. 2017. The autophagic machinery in viral exocytosis. Frontiers in Microbiology **8**, 269.

Nedelsky NB, Todd PK, Taylor JP. 2008. Autophagy and the ubiquitin-proteasome system: Collaborators in neuroprotection. BBA Molecular Basis of Disease **1782**, 691–699.

Noh YS, Amasino RM. 1999. Identification of a promoter region responsible for the senescence-specific expression of SAG12. Plant Molecular Biology **41**, 181–194.

Otegui MS, Noh YS, Martínez DE, Vila Petroff MG, Staehelin LA, Amasino RM, Guiamet JJ. 2005. Senescence-associated vacuoles with intense proteolytic activity develop in leaves of Arabidopsis and soybean. The Plant Journal **41**, 831–844.

Pathirana R, West P, Hedderley D, Eason J. 2017. Cell death patterns in Arabidopsis cells subjected to four physiological stressors indicate multiple signalling pathways and cell cycle phase specificity. Protoplasma **254**, 635–647.

Polge C, Jaquinod M, Holzer F, Bourguignon J, Walling L, Brouquisse R. 2009. Evidence for the existence in *Arabidopsis thaliana* of the proteasome proteolytic pathway: activation in response to cadmium. The Journal of Biological Chemistry **284**, 35412–35424.

Poret M, Chandrasekar B, van der Hoorn RAL, Avice J-C. 2016. Characterization of senescence-associated protease activities involved in the efficient protein remobilization during leaf senescence of winter oilseed rape. Plant Science **246**, 139–153.

Pružinská A, Shindo T, Niessen S, Kaschani F, Tóth R, Millar AH, van der Hoorn RA. 2017. Major Cys protease activities are not essential for senescence in individually darkened Arabidopsis leaves. BMC Plant Biology **17**, 4.

Rogers JC. 2013. Aleurain. In: Rawlings ND, Salvesen GS, eds. Handbook of proteolytic enzymes, vol. 2; 3rd edition. Academic Press, 1888–1891.

Rustgi S, Boex-Fontvieille E, Reinbothe C, von Wettstein D, Reinbothe S. 2017. Serpin1 and WSCP differentially regulate the activity of the cysteine protease RD21 during plant development in *Arabidopsis thaliana*. Proceedings of the National Academy of Sciences, USA **114**, 2212–2217.

Thompson AR, Doelling JH, Suttangkakul A, Vierstra RD. 2005. Autophagic nutrient recycling in Arabidopsis directed by the ATG8 and ATG12 conjugation pathways. Plant Physiology **138**, 2097–2110.

Tsukada M, Ohsumi Y. 1993. Isolation and characterization of autophagy-defective mutants of *Saccharomyces cerevisiae*. FEBS Letters **333**, 169–174.

Valot B, Langella O, Nano E, Zivy M. 2011. MassChroQ: a versatile tool for mass spectrometry quantification. Proteomics **11**, 3572–3577.

van der Hoorn RAL. 2013. RD21 Peptidase. In: Rawlings ND, Salvesen GS, eds. Handbook of proteolytic enzymes, vol. 2; 3rd edition. Academic Press, 1892–1896.

van der Hoorn RAL, Leeuwenburgh MA, Bogyo M, Joosten MH, Peck SC. 2004. Activity profiling of papain-like cysteine proteases in plants. Plant Physiology **135**, 1170–1178.

van Wijk KJ. 2015. Protein maturation and proteolysis in plant plastids, mitochondria, and peroxisomes. Annual Review of Plant Biology **66**, 75–111.

Waite KA, De-La Mota-Peynado A, Vontz G, Roelofs J. 2016. Starvation induces proteasome autophagy with different pathways for core and regulatory particles. The Journal of Biological Chemistry **291**, 3239–3253.

Wang J, Ding Y, Wang J, Hillmer S, Miao Y, Lo SW, Wang X, Robinson DG, Jiang L. 2010. EXPO, an exocyst-positive organelle distinct from multivesicular endosomes and autophagosomes, mediates cytosol to cell wall exocytosis in Arabidopsis and tobacco cells. The Plant Cell **22**, 4009–4030.

Wang JC, Ma C, Zhang M, Yang L, Chen WL. 2015. ATG5 is required to limit cell death induced by *Pseudomonas syringae* in Arabidopsis and may be mediated by the salicylic acid pathway. Acta Physiologiae Plantarum **37**, 1731.

Wang Z, Gu C, Colby T, Shindo T, Balamurugan R, Waldmann H, Kaiser M, van der Hoorn RA. 2008. β-Lactone probes identify a papain-like peptide ligase in *Arabidopsis thaliana*. Nature Chemical Biology **4**, 557–563.

Watanabe N, Lam E. 2011. Arabidopsis metacaspase 2d is a positive mediator of cell death induced during biotic and abiotic stresses. The Plant Journal **66**, 969–982.
**Yang X, Bassham DC.** 2015. New insight into the mechanism and function of autophagy in plant cells. In: Jeon KW, ed. International review of cell and molecular biology, vol. 320. Academic Press, 1-40.

**Yang XC, Srivastava R, Howell SH, Bassham DC.** 2016. Activation of autophagy by unfolded proteins during endoplasmic reticulum stress. The Plant Journal **85**, 83-95.

**Yao X, Xiong W, Ye T, Wu Y.** 2012. Overexpression of the aspartic protease ASPG1 gene confers drought avoidance in Arabidopsis. Journal of Experimental Botany **63**, 2579-2593.

**Yoshimoto K, Jikumaru Y, Kamiya Y, Kusano M, Consonni C, Panstruga R, Ohsumi Y, Shirasu K.** 2009. Autophagy negatively regulates cell death by controlling NPR1-dependent salicylic acid signaling

during senescence and the innate immune response in Arabidopsis. The Plant Cell **21**, 2914-2927.

**Yu F, Park S, Rodermel SR.** 2004. The Arabidopsis FtsH metalloprotease gene family: interchangeability of subunits in chloroplast oligomeric complexes. Plant Journal **37**, 864-876.

**Zhang D, Liu D, Lv X, Wang Y, Xun Z, Liu Z, Li F, Lu H.** 2014. The cysteine protease CEP1, a key executor involved in tapetal programmed cell death, regulates pollen development in Arabidopsis. The Plant Cell **26**, 2939-2961.

**Zientara-Rytter K, Sirko A.** 2016. To deliver or to degrade – an interplay of the ubiquitin-proteasome system, autophagy and vesicular transport in plants. The FEBS Journal **283**, 3534-3555.