Identification of Autophagosome-associated Proteins and Regulators by Quantitative Proteomic Analysis and Genetic Screens*

Autophagy is one of the major intracellular catabolic pathways, but little is known about the composition of autophagosomes. To study the associated proteins, we isolated autophagosomes from human breast cancer cells using two different biochemical methods and three stimulus types: amino acid deprivation or rapamycin or concanamycin A treatment. The autophagosome-associated proteins were dependent on stimulus, but a core set of proteins was stimulus-independent. Remarkably, proteasomal proteins were abundant among the stimulus-independent common autophagosome-associated proteins, and the activation of autophagy significantly decreased the cellular proteasome level and activity supporting interplay between the two degradation pathways. A screen of yeast strains defective in the orthologs of the human genes encoding for a common set of autophagosome-associated proteins revealed several regulators of autophagy, including subunits of the retromer complex. The combined spatiotemporal proteomic and genetic data sets presented here provide a basis for further characterization of autophagosome biogenesis and cargo selection.

Macroautophagy (hereafter referred to as autophagy) is an evolutionarily conserved lysosomal pathway involved in the turn-over of long-lived proteins, cytoplasm, and whole organelles (1)(2)(3)(4). Its deregulation in mice leads to embryonic or perinatal lethality (5,6), neurodegeneration (7), or cancer (8), emphasizing the physiological importance of this catabolic process. Although autophagy occurs constitutively at a low basal level, starvation, growth factor deprivations, and protein aggregation, as well as other cellular stresses rapidly increase its activity. Under these conditions, autophagy is crucial for generating nutrients or removing damaged cytoplasmic components, thereby serving mainly as a protective cellular response (2).
The autophagic process begins with the nucleation of a flat membrane cistern that enwraps cytoplasmic organelles and/or a portion of the cytosol. The membrane elongates until the edges of the membrane fuse, thereby forming a double membrane structure called an autophagosome, which fuses with endosomes forming an amphisome (9) and subsequently matures to an autolysosome by fusing with lysosomal vesicles. The final degradation of the cargo takes place within autolysosomes, where lysosomal hydrolases digest the luminal content, allowing the recycling of amino acids, nucleotides, and fatty acids (10). The process is controlled by a set of evolutionarily conserved autophagy-related proteins (Atg proteins) initially identified in yeast (11,12). The protein complex consisting of phosphatidylinositol 3-phosphate kinase class 3, p150 myristylated protein kinase, and beclin 1 (Atg6) is essential for the initial membrane assembly, whereas the following membrane elongation depends on two ubiquitinlike conjugation systems. One of them converts microtubuleassociated protein 1 light chain 3 (LC3/Atg8) 1 from its free form (LC3-I) to a phosphatidylethanolamine-conjugated form (LC3-II), which associates with both membranes of the autophagosome (13). This process is frequently used as an autophagy marker because the change in the LC3 staining pattern from diffuse to dotted can be readily visualized.
Autophagy is generally considered an unselective bulk degradation pathway. However, under certain conditions, autophagosomes have been suggested to selectively remove, for example, damaged mitochondria (14), endoplasmic reticulum (ER) (15), peroxisomes (16), ribosomes (17), and the midbody ring at the end of cytokinesis (18). Moreover, our recent proteomic analysis of starved cells shows that cellular proteins decrease in an ordered fashion depending on their subcellular localization (19). Thus, autophagy may also serve as a specific degradation system, similar to the proteasome, which recognizes ubiquitin-coupled proteins for degradation (20). Even though the ubiquitin-proteasome system and autophagy have long been viewed as complementary degradation systems with no point of intersection, it was shown recently that autophagy can act compensatorily when the ubiquitin-proteasome system is impaired in Drosophila melanogaster (21). These data suggest that there might be a link between the two major cellular proteolysis pathways (22)(23)(24).
The aim of this study is to identify proteins associated with the mature autophagosome and to compare the protein composition of autophagosomes induced by different stimuli. For this purpose, we analyzed autophagosomes isolated from MCF7 breast cancer cells following amino acid starvation or treatment with either rapamycin (an inhibitor of the mammalian target of rapamycin complex 1 (mTORC1)) or concanamycin A (an inhibitor of the lysosomal H ϩ -ATPase) by quantitative MS-based proteomics (25) relying on protein correlation profiling (PCP) (26) and stable isotope labeling by amino acids in cell culture (SILAC) (27,28). We identified 728 putative autophagosome-associated proteins from a background of co-purifying proteins. Of these, only a total of 94 proteins were common to all stimuli. To validate that the co-migrating proteins were derived from autophagosomes, we performed GFP pulldowns directed against a GFP-tagged version of the autophagosomal marker LC3 and tested subcellular localization by fluorescent microscopy for selected candidate proteins. To test whether the common autophagosome-associated proteins functioned as autophagy regulators, we turned to yeast genetics and screened the available yeast strains with mutations in the corresponding genes for changes in autophagy by the alkaline phosphatase (ALP) assay (29) identifying inter alia the retromer complex. Human cells depleted for the novel autophagy regulators identified in the yeast screen were subsequently analyzed for changes in autophagic flux by cell assays detecting the turnover of autophagosome-associated proteins (30,31). The data presented provide a quantitative proteomic characterization of autophagosomes identifying a set of autophagosome-associated proteins and regulators.
Transfections-Fusion constructs of DsRed and RHEB or FKBP1A were generated using vectors pDsRed-C1 and -N1, respectively, (Clontech). RHEB and FKBP1A cDNAs were purchased from im-aGenes (Berlin, Germany), and fusion constructs were sequence checked. Stable cell lines were generated by liposomes transfection using Metafectene TM (Biontex, Martinsried, Germany). During transient transfections FuGENE 6 (Roche Applied Science) was used as a transfection agent according to the manufacturer's instructions. siRNA sequences corresponding to the human cDNA sequence are listed in the supplemental material. The cells were transfected with 50 nM siRNA employing Oligofectamine TM (Invitrogen).
Immunodetection-The primary antibodies used for immunoblot analysis are listed in the supplemental "Experimental Procedures." The appropriate peroxidase-conjugated secondary antibodies were from Dako A/S (Glostrup, Denmark) and Vector Laboratories (Burlingame, CA), and the ECL immunoblotting reagents were from Amersham Biosciences.
The cells were fixed and stained applying standard procedures (49). For detection of DsRed fusion proteins, the cells were fixed and permeabilized in 3.7% paraformaldehyde for 7 min followed by methanol overnight. The primary antibodies are listed in supplemental "Experimental Procedures." The respective fluorescent secondary antibodies were from Molecular Probes or Jackson Immunolabs. The images were obtained using either a Zeiss LSM 510 META confocal laser scanning microscope or Zeiss Axiovert 100M confocal laser scanning microscope (Carl Zeiss, Jena, Germany).
Detection of Autophagosomes-Percentages of cells with eGFP-LC3 translocation into dots (a minimum of 100 cells/sample) were counted blindly in cells fixed in 3.7% formaldehyde and 0.19% picric acid (v/v) applying a Zeiss Axiovert 100M confocal laser scanning microscope (Carl Zeiss).
Autophagic Flux-The autophagic flux was measured applying a luciferase-based assay as described previously (30).
Proteasome Activity Assay-The proteasome activity was measured using the 20 S proteasome activity assay kit from Chemicon International (catalog number APT280; Temecula, CA), and the protein levels were measured using the Bio-Rad Dc protein assay kit (Bio-Rad) according to the manufacturer's instructions (supplemental "Experimental Procedures").
Autophagy Measurements in Yeast-The experiments were carried out with BY4741 (MATa his3⌬1 leu2⌬0 met15⌬0 ura3⌬0) and respective null mutants, obtained from Euroscarf (growth conditions described in the supplemental "Experimental Procedures"). For induction of autophagy by nitrogen starvation, the cells were inoculated from fresh overnight cultures to 0.4 A 600 (ϳ4 ϫ 10 6 cells/ml) in Synthetic Complete Dextrose, grown for 4 -5 h to mid-log phase reaching ϳ1.5 A 600 , washed twice in double distilled H 2 O, and incubated in SD-N medium at 1 A 600 for 3 h. For induction of autophagy by rapamycin treatment, the cells were grown to mid-log phase (ϳ1 A 600 ) as above, and rapamycin (AG Scientific) was added to the culture to a final concentration of 0.5 g/ml preceding 3 h of incubation.
Autophagy was determined in cellular extracts of 1-ml culture aliquots by ALP activity according to Ref. 50 using respective strains transformed with and selected for stable insertion of pTN9 HindIII fragment containing the cytosolic form of Pho8p (Pho8p⌬N60). For details see the supplemental "Experimental Procedures".
Autophagosome Purification and Mass Spectrometry-Autophagosome purification was modified from the protocol described in Ref. 51 (for details see the supplemental "Experimental Procedures"). The PCP-SILAC background fractions (Lys 0 /Arg 0 ) were combined and spiked in a 1:1 ratio to the respective (Lys 4 /Arg 6 ) PCP SILAC fractions (52). The resulting fractions were dilluted 1:1 in Homogenization Medium and pelleted at 31,000 ϫ g. The pellets were resuspended in SDS loading buffer and reduced and alkylated. Protein mixtures were separated by SDS-PAGE (4 -12% Bis-Tris gradient gel, NuPAGE; Invitrogen) and in-gel digested with trypsin. In the GFP tag enrichment experiment, vesicles were enriched by centrifugation followed by GFP antibody pulldown (supplemental Fig. S3). The proteins were digested in solution with trypsin and the peptides separated by ion exchange chromatography. The resulting peptide mixtures were analyzed by online C18 reversed phase nanoscale liquid chromatography tandem mass spectrometry essentially as described (53) (for details see the supplemental "Experimental Procedures").
Assigning Peptide Sequences Using MASCOT and Quantitation Using MSQuant-Raw MS2 spectra were centroided and merged into a single peak list file using the in-house-developed software DTASupercharge (default values; version 1.19, msquant.sourceforge.net) and searched against the human MSIPI database (v. 3.34 with 68404 entries) (54) using MASCOT 2.0 (Matrix Science, London, UK) basically as described (19) and the following parameters: carbamidomethyl-cysteine was set as fixed modification, methionine oxidation, deamidation of asparagine and glutamine, and protein aminoterminal acetylation were set as variable modifications. Double or triple SILAC were chosen as quantification mode. Two miss cleavages were allowed, enzyme specificity was trypsin, precursor mass accuracy had to be within 30 ppm for FT-data and 7 ppm for Orbitrap data, and the fragment spectra mass accuracy was set at 0.6 Da. The identified peptides were recalibrated using MSQuant (55), and the results were combined using MGFcombiner (version 1.05) and searched again using MASCOT 2.0 with the above mentioned parameters, except that the precursor mass tolerance was set to 5 ppm. To determine the number of false positive peptide hits, the data were searched against a full-length human MSIPI decoy database, and the MASCOT peptide score was adjusted to yield a number of false positive peptide identifications of less than 1% (calculated as follows: false positive rate (%) ϭ reverse hits ϫ 2 ϫ 100/forward hits). For a protein to be counted as identified, a minimum of two unique peptides (bold red hits, minimum length of 7 amino acids) had to be sequenced and to fulfill the determined criteria. To yield the maximum number of peptides per identified protein, the data were researched against decoy databases consisting of the identified proteins only considering proteins that were identified using a full-length decoy database and that did pass the stringent identification criteria stated above. Thus, the peptide MASCOT identification score yielding less than 1% false positive could be dropped to 15, resulting in considerably more peptide IDs/protein. PCP SILAC profiles were calculated using all bold red and parenthesized peptides. The proteins were only included if respective peptides were sequenced in each of the six fractions yielding six data point profiles. In a few exceptions, the peptides were inserted based in on elution time predicted from the analyses of neighbor fractions. The peptides were quantified by MSQuant (55) using extracted ion chromatograms of the monoisotopic signals. For cluster analyses, the protein ratios were normalized using the base 2 logarithm.
The GFP pulldown data was processed and analyzed using Max-Quant (version 1.0.13.13) (56) and Mascot (version 2.3). The Quant module of MaxQuant was used with parameters set as follows: Arg 6 / Lys 4 and Arg 10 /Lys 8 as triple SILAC labels, a maximum of two missed cleavages, filtering of MS/MS spectra to retain only the six most intense peaks/100 Da, and MS/MS mass tolerance of 0.5 Da. Variable modifications were methionine oxidation and protein amino-terminal acetylation. Cysteine carbamidomethylation was a fixed modification. Protein and peptide FDR of 0.01, PEP based on Mascot score, minimum peptide length of 6, minimum score of 7, minimum unique sequence of 1, minimum peptides of 1, use of both unmodified and modified peptides for protein quantification, and use of both razor and unique peptides for quantification with a minimum ratio count of 2.
Clustering of Profiles and Comparative Analysis-The protein profiles were clustered using GProX (57) to identify proteins with profiles similar to known autophagosomal proteins. Clustering was done using the fuzzy c-means algorithm (38), which is a soft clustering algorithm being noise robust and indicating how well protein profiles are represented by respective clusters. The Fuzzy c-means parameter m varied between 1.25 (Rapa), 2.5 (HBSS), and 3.0 (ConA) for the three stimuli to compensate for differences in resolution of the data and to ensure minimal size of autophagosomal clusters. Proteins were considered to being autophagosomal if they clustered with MAP1LC3, p62, and GABA A receptor-associated protein-like 2 (GABARAPL2). To test the validity of the clusters, it was required that major complexes like the ribosome, the exosome, and the proteasome located to a single cluster (see the supplemental "Experimental Procedures").

RESULTS
Quantitative Mass Spectrometric Analysis of Autophagosomes by PCP-SILAC-Extensive efforts have focused on the description of protein networks involved in autophagy signaling (32). However, the protein composition of autophagosomes is still largely unknown, and only a few autophagosome-associated proteins have been identified (33,34). One of them is the frequently used autophagosomal marker LC3-II. In this study, we took advantage of MCF7 breast carcinoma cells stably expressing LC3 fused to enhanced green fluorescent protein (eGFP-LC3) to follow the accumulation of autophagosomes after 7-h treatments by three well characterized stimuli, i.e. amino acid starvation, rapamycin, and concanamycin A (31) (Fig. 1, A and B). This experimental system was then used to characterize the stimuli-dependent protein composition of autophagosomes and the function of the identified candidate proteins using a combination of quantitative proteomics and genetic analyses (Fig. 1A).
To distinguish autophagosomal candidate proteins from nonspecific co-purifying proteins, we applied the PCP-SILAC method (35,36). This method is based on the assumption that proteins residing in the same organelle will have similar distributions through a density gradient. As outlined in Fig. 1C, two preparation of autophagosomes were isolated in parallel from MCF7-eGFP-LC3 cells labeled with different combinations of isotope-encoded lysine and arginine (37). Fractions collected after the final iodixanol density gradient centrifugation (1.05-1.25 g/ml iodixanol) were tested for the presence of organellar markers and known autophagosomal proteins by Notably, the cells expressing a mutated form of eGFP-LC3-G120A unable to become lipid-conjugated did not form dots upon stimulation (data not shown), strongly suggesting that the dots observed are autophagosomes rather than unspecific aggregates. C, to identify autophagosome-associated proteins using the PCP-SILAC methods, the cells were isotope-labeled and subsequently treated for 7 h with 2 nM ConA. Autophagosomes were purified by gradient centrifugation, six fractions were collected, and the Lys 0 /Arg 0 fractions were combined, yielding an internal standard mixture of proteins over the gradient. This internal standard was distributed in a 1:1 ratio to the original Lys 4 /Arg 6 labeled fractions, and the combined samples were separated by SDS-PAGE, in-gel digested by trypsin, and analyzed by MS. The experiment was repeated using 100 nM rapamycin and starvation in HBSS as stimuli. Applying this setup, we were able to identify 7935 proteins Western blot analysis ( Fig. 1E and supplemental Fig. S1). For peptide isotope ratio determination, six fractions containing autophagosomal markers from the unlabeled cells were combined to generate a common internal standard that was equally distributed into the corresponding six fractions from the labeled cells. Proteins in these six samples were then analyzed by LC-MS, and the resulting data were used to determine the relative enrichment of proteins in each fraction (Fig. 1D). The protein enrichment profile of LC3 and Western blot analyses of eGFP-LC3 indicated that the majority of autophagosomes appeared in fractions 2 and 3 regardless of the stimulus (Fig. 1E). In contrast, marker proteins for other organelles showed clearly distinct profiles (supplemental Fig. S1B). In total, 7935 proteins were identified in single representative PCP-SILAC experiments for each of the three stimuli (from a total set of one concanamycin A experiment, two rapamycin experiments, and two starvation experiments; supplemental Fig. S2). Complete protein enrichment profiles were obtained for 4516 proteins detected in all six fractions (supplemental Fig. S2 and Tables S1-S3). These 4516 proteins were grouped into three clusters per experiment based on their profiles using the noise robust fuzzy c-means algorithm (38). Cluster A proteins peaked in fractions 2 and 3 and included known autophagosomal proteins such as LC3, p62 (SQSTM1), and GABARAPL2 (Figs. 2A and supplemental Fig. S2). Cluster B contained proteins that also peaked in fractions 2 and 3 but showed an additional rise in fraction 6. Cluster C included proteins that peaked in fraction 5 and 6 and were regarded as nonspecific copurifying proteins. Autophagosomal clusters (cluster A) from concanamycin A-treated, rapamycin-treated, and starved cells contained 238, 359, and 482 protein profiles, respectively (supplemental Tables S1-S3). These profiles represent 728 different proteins, of which 94 were found in the autophagosomal cluster for all three stimuli (Table I and Fig.  2B). In support of the PCP-SILAC method, we identified 42 of the 94 proteins in cluster A that were known autophagosomal proteins or recently reported interaction partners of autophagy-related proteins (Table I and supplemental Table S4). Identification of autophagy regulators such as FK506-binding protein 1 (FKBP1A) and Ras homolog enriched in brain (RHEB) suggests that they operate at the autophagosome scaffold.
Using Gene Ontology (39), we analyzed the subcellular localization of proteins detected in the same clusters independent of the autophagosome-stimulating treatment. Because proteins may carry more than one Gene Ontology term, the sum of all annotated proteins may be higher than 100%.
The majority of proteins in the autophagosomal cluster A were annotated as cytoplasmic/cytosolic (Fig. 2C). In addition to the cytoplasmic/cytosolic proteins, cluster B was enriched in proteins from the endosomal compartment, ER, Golgi, cellular membranes, and ribosomes, whereas cluster C contained mainly membrane, extracellular, and ER proteins. Hence, we were able to distinguish different vesicular compartments by our PCP-SILAC approach and identified cluster A as the autophagosomal cluster.
Quantitative Mass Spectrometric Analysis of Immunopurified Autophagosomes-To further support the enrichment of autophagosome-associated proteins in cluster A and to distinguish them from co-migrating vesicular structures, we purified autophagosomes from MCF7-eGFP-LC3 cells by immunoprecipitation using an anti-GFP antibody (Fig. 3A) (33). Cell extracts were depleted for the free pool of LC3 by centrifugation to minimize binding of the nonautophagosomal pool of proteins. Specific autophagosome-associated proteins were determined by SILAC-based mass spectrometry comparing the relative abundance of proteins isolated from untreated and starved cells (Fig. 3A). Starvation induces a high autophagic flux (see Fig. 5C) challenging the identification of protein enrichment in affinity-purified autophagosomes. Therefore, a third condition was analyzed: starved cells treated additionally with concanamycin A, which leads to a pronounced accumulation of autophagosomes by blocking lysosomal degradation ( Fig. 3A and data not shown). Contaminating proteins should be equally abundant in all three conditions and hence have 1:1:1 ratios, whereas autophagosomal proteins such as LC3 and p62 display higher ratios in the ConA sample where autophagosomes accumulate (supplemental Table S5). Indeed, common cluster A proteins show a nonrandom distribution and are significantly enriched compared with cluster B and C proteins supporting autophagosome association for cluster A proteins (p Ͻ 0.001; Fig. 3, B and C, and supplemental Fig. S3 and Table S5). A second eGFP-IP was performed on cells either starved or treated with rapamycin and ConA and confirmed the enrichment of cluster A proteins (supplemental Table S5). In total, 65 common cluster A proteins were identified in the eGFP pulldowns (Table I). Furthermore, immunofluorescence microscopy validated the co-localization of the three cluster A proteins p62, RHEB, and FKBP1A with LC3 following all three stimuli ( changes in the composition and abundance of proteins associated with autophagosomes, we compared proteins identified in cluster A between the different PCP-SILAC experiments and performed additional comparative and temporal quantitative experiments. Interestingly, PCP-SILAC revealed stimuli-dependent groups of proteins in the autophagosomal clusters A (Fig. 5A). For instance, starvation triggered autophagosomes associated with plasma membrane proteins and, similar to rapamycin, vesicular proteins. These classes of proteins were largely absent in autophagosomes purified from concanamycin A-treated cells. We compared protein abundances directly between autophagosomes isolated from starved, rapamycin-treated, and concanamycin A-treated cells. Proteins known to be involved in autophagy, such as LC3, p62, and GABARAPL2, were significantly less abundant in autophagosomes from starved cells as compared with rapamycin-and concanamycin A-induced autophagosomes ( Fig. 5B and supplemental Table S6). This suggests that starvation induces the highest level of autophagic flux as confirmed by the rapid depletion of p62 in starved cells (Fig. 5C). Conversely, a few proteins comprising several members of the annexin family were most abundant in autophagosomes from starved cells (Fig. 5D). Proteasomal proteins on the other hand were present in similar amounts regardless of the stimulus (Fig. 5B).
To investigate the recruitment dynamics of proteins to the autophagosome, we starved cells for 3, 6, and 12 h, isolated autophagosomes, and analyzed the changes in the abundance of autophagosome-associated proteins (supplemental Table S7). p62 decreased most rapidly within the time frame of analysis, whereas proteasomal proteins (PSMB1) and, for example, eukaryotic translation elongation factor 1␥ (EEF1G) showed a slight increase (Fig. 5E). Taken together, this suggests that autophagosomes are highly dynamic struc-  Fig. S2. B, Venn diagram of the cluster A proteins identified from the HBSS-, Rapa-, and ConA-treated cells. Whereas 94 proteins were in common for the three stimuli, the majority of proteins were detected by only one or two stimuli. Shown are three representative data sets of five biological replicates. C, subcellular localization of common proteins in clusters A-C based on Gene Ontology.  tures whose composition is influenced by the nature and timing of the stress stimulus. Autophagy Modulates Proteasomal Activity-Our findings suggest that the autophagic machinery and the proteasome are interconnected (Figs. 2C and 5A). Interplay between the two degradation systems has been proposed lately (24). However, because a detailed description of the connection of these two major cellular degradation pathways is still missing, we investigated their interplay in more detail. First, we confirmed the stimulus-independent partial co-localization of proteasomal proteins and LC3 by fluorescence microscopy and Western blot analyses (Fig. 6, A-J). This evidence of proteasomal 20 S subunits associated with autophagosomes led us to speculate whether autophagy might decrease the proteasome level in cells. Indeed, the abundance of proteasomal proteins in whole cell lysates decreased upon amino acid starvation, and importantly this decrease could be blocked by the addition of 3-methyladenine, an inhibitor of autophagy (40) (Fig. 6K). Accordingly, induction of autophagic flux by rapa-mycin or starvation led to a significant decrease in proteasomal activity (Fig. 6L). Notably, concanamycin A, which inhibits the degradation of autophagosomal cargo (Fig. 5C), did not influence the proteasomal activity (Fig. 6L). Furthermore autophagosomes purified by immunoprecipitation of eGFP-LC3 with an anti-GFP antibody contained proteasomal subunits, supporting their association with autophagosomes. Interestingly, incubation of these purified autophagosomes with detergent abolished the pulldown of the proteasomal subunits. This suggests that the proteasomal subunits do not directly interact with LC3 but rather associate with the autophagosome. Contrary, p62 was pulled down with eGFP-LC3 independent of detergent treatment as expected (Fig.  6M). Taken together, these data suggest that the proteasome associates with autophagosomes independently of LC3 and that functional autophagy leads to a decrease in proteasome amount and activity.
Functional Analysis of Common Autophagosome-associated Proteins-Our experimental setup did not allow us to directly distinguish between cargo and core proteins. p62 is required for starvation-induced LC3 translocation (41), and we speculated whether other candidates involved in autophagy regulation are among the 94 common autophagosome-associated proteins. To investigate this further, we performed an ortholog screen in the yeast Saccharomyces cerevisiae whose genome contains 59 orthologs of the human genes encoding for the 94 common autophagosome-associated proteins, giving rise to 32 viable knock-out strains (Table I). In addition, we included two proteins involved in vesicle trafficking, CAP1 and VPS35, which were present in one or two of the three autophagosomal clusters, enriched in the autophagosome pulldowns and also giving rise to viable yeast knock-out strains. All of the strains were starved and tested for their autophagic response by the ALP assay (29) indicative of autophagic flux (Fig. 7A). Strains deficient in the yeast orthologs of mammalian CAP1 (SRV2), EEF1G (TEF4), LC3 (ATG8), NP (PNP1), RHEB (RHB1), and VPS35 (VPS35) had significantly lower autophagic activity, whereas the yeast strain deficient in the ortholog of mammalian GNB2L1 (ASC1) had significantly higher autophagic activity upon starvation as compared with the wild-type strain (Fig. 7A). Except for the yeast strain deficient for TEF4, similar changes in the autophagic flux were observed when autophagy was triggered by rapamycin treatment (Fig. 7B).
The strong phenotype observed for the vacuolar protein sorting protein VPS35 suggested a so far unknown role of the retromer complex in autophagy. This complex mediates retrograde transport of cargo from endosomes to the trans Golgi network and is composed of the cargo selective subunits VPS26, VPS29, and VPS35 and the sorting nexin subunits VPS5 and VPS17 (42). To corroborate the phenotype data for VPS35, we analyzed additional subunits of the retromer complex. Interestingly, strains deficient for the subunits VPS29, VPS5, and VPS17 revealed equally low autophagic activity in cells starved or treated with rapamycin. These data further support that the retromer complex is required for TOR-dependent autophagy in yeast (Fig. 7C).
To test whether the autophagy regulating functions of the yeast proteins were conserved in mammalian cells, we depleted MCF7 cells for the ortholog proteins using RNA interference and measured the level of p62 in control conditions and following autophagy induction by amino acid starvation or etoposide treatment, which we have used recently to establish a luciferase-based autophagy screen (30) (Fig. 7D). No major differences in p62 levels in control conditions were observed upon siRNA-mediated depletion of CAP1, GNB2L1, or VPS35, whereas depletion of RHEB reduced the cellular p62 level, suggesting that RHEB siRNA induced spontaneous autophagy ( Fig. 7D and data not shown). Accordingly, two nonoverlapping RHEB siRNAs increased the autophagic flux as measured by a decrease in the ratio of RLuc-LC3-WT/ RLuc-LC3 G120A (30) (Fig. 7E) and increased the number of cells with eGFP-LC3 translocation (Fig. 7F). As expected for the control experiments, the induction of autophagy by amino acid starvation and etoposide decreased the level of p62, and etoposide decreased the ratio of RLuc-LC3-WT/RLuc-LC3 G120A and increased the number of cells with eGFP-LC3 translocation in cells treated with control siRNA (Fig. 7, D-F). All of these changes indicative of autophagy were inhibited in cells depleted for EEF1G by two nonoverlapping siRNAs (Fig.  7, D-F, and data not shown). Interestingly, depletion of EEF1G also attenuated etoposide-and starvation-induced inhibition of mTORC1 as measured by levels of p70S6K phosphorylation (supplemental Fig. S6). Autophagy induced by starvation but not that induced by rapamycin is defective in yeast deficient for the EEF1G ortholog TEF4 (Fig. 7, A and B). This suggests that EEF1G is a positive regulator of autophagy required for autophagy signaling upstream of mTORC1.

DISCUSSION
The data presented above provide a framework of autophagosome-associated proteins and will help in characterizing underlying cellular processes regulating the biogenesis and cargo selection of autophagosomes. Employing PCP-SILAC, we grouped the identified proteins in three clusters per organelle preparation for cells treated with concanamycin A, rapamycin, and HBSS. The clusters consisted of autophagosomal candidate proteins (cluster A); proteins found both in autophagosomal fractions and in other organelles such as endolysosomes, ER, and Golgi apparatus (cluster B); and nonautophagosomal proteins such as ER and nuclear pro-teins (cluster C). For further investigations, we concentrated on proteins observed in cluster A.
The most comprehensive proteomic analysis of autophagosomes prior to this study identified 39 and 101 proteins localizing to autophagosomal membranes in starved rat hepatocytes and human cell lines, respectively (33,34). We chose a different approach characterizing not only autophagosomal membrane proteins but also autophagosomal cargo proteins to learn more about cellular proteome dynamics during autophagy. With this approach, we identified 728 autophagosomal candidate proteins, of which 94 were found in the autophagosome cluster regardless of the stimulus. Our data suggest that 31 proteins previously linked to autophagy are also associated with the autophagosome, among them the proteasome.
Comparison of the 94 common proteins with the other two reported data sets revealed seven human ortholog proteins shared with the autophagosomal rat proteins identified by Øverbye et al. (34) and a single protein shared with the 101 identified by Gao et al. (33). The rather small overlap might reflect differences in stimuli, isolation procedures, cell types, FIG. 5. Autophagosome-associated protein dynamics. A, comparison of cluster A proteins identified by PCP-SILAC indicates stimulidependent autophagosome-protein association for Rapa-and ConA-treated and starved cells based on Gene Ontology analysis. B, direct comparison of autophagosome-associated proteins by SILAC-based quantitative analysis of autophagosomes isolated from MCF7 cells treated with 100 nM Rapa or 2 nM ConA or starved in HBSS for 7 h and mixed in a ratio of 1:1:1. The relative abundance of proteins per autophagosome is normalized to the ConA isotope signal. Shown are the combined results from two biological replicates. C, p62 levels in MCF7 cells treated with 1 M Rapa or 2 nM ConA or starved in HBSS for the indicated time periods were analyzed by immunoblotting using GAPDH as loading control. D, relative abundances of indicated annexins in autophagosomes from cells treated and analyzed as in B. Shown are the combined results from two biological replicates. E, autophagosome-associated protein dynamics. Autophagosomes were isolated from SILAC labeled cells and combined in a ratio of 1:1:1 after starvation for 3, 6, and 12 h. Protein-dependent targeting dynamics was observed. Shown are the combined results from two biological replicates. and MS methods. Whereas Øverbye et al. did not identify any of the known ATG proteins, we determined complete protein profiles of the two human ATG8 homologs, MAP1LC3B and GABARAPL2. Gao et al. identified MAP1LC3B, ATG7 and ATG9. We identified ATG3 and ATG4B and detected single peptides from ATG9A and ATG5. The abundance of the latter proteins was too low for the extraction of complete protein profiles, and therefore they were omitted from further analysis.
However, the partial profiles of these proteins closely resemble the profiles of proteins in cluster A (data not shown), indicating that these ATGs may also localize to autophagosomes as suggested previously (43). Taken together, ATG proteins appear to be low abundant and are not commonly identified by large scale proteomics studies. For their analyses, dedicated affinity purification protocols appear to be more promising (32). Shown are bands for the 20 S core subunits, which follow the MS profiles in all three stimuli. K, the relative abundance of proteasomal subunits were determined by SILAC-based mass spectrometry of cells left untreated or starved for 12 h with or without 10 mM 3-methyladenine combined in a ratio of 1:1. Shown are the relative changes compared with control cells (average ratio of detected PSMA, PSMB, PSMC, and PSMD proteins; the error bars indicate standard deviations). *, p Ͻ 0.01 as analyzed by a one-sample t test. L, changes in proteasome activity in response to autophagy were analyzed in lysates of MCF7 cells left untreated (control) or treated for 24 h with 2 nM ConA or 1 M Rapa or starved for amino acids in HBSS. The values are percentages of proteasome activity/protein concentration as compared with untreated control samples and represent the averages Ϯ S.D. from four independent experiments. *, p Ͻ 0.01 as analyzed by a one-sample t test. M, proteasome association with LC3 affinity-purified autophagosomes was analyzed by SILAC-based mass spectrometry using MCF7-eGFP-LC3 cells left untreated (control) or stimulated with 2 nM ConA for 7 h. Anti-GFP immunoprecipitations were performed in lysis buffer with or without 1% Nonidet P-40. Without detergent, intact autophagosomes were purified. Under these conditions, enrichment of proteasomal proteins (average ratio of detected PSMA, PSMB, PSMC, and PSMD proteins) was observed similar to p62/SQSTM1. In the presence of detergent, autophagosomes were destroyed, and the proteasomal proteins were no longer enriched in contrast to proteins binding directly to LC3 such as SQSTM1. The values represent the averages from two independent experiments Ϯ S.D. Ctrl, control.

FIG. 7. Genetic screens for novel autophagy regulators in yeast and mammalian cells.
A and B, ALP screen in S. cerevisiae. The ALP activity was measured in indicated S. cerevisiae knock-out strains starved for nitrogen (A) or treated with 500 ng/ml rapamycin (B) for 3 h. A histidine-dependent strain and an ATG7 knock-out strain were used as negative and positive controls, respectively. The values are percentages of ALP activity/g protein as compared with the wild-type control strain and represent the averages Ϯ S.D. from three independent experiments. *, p Ͻ 0.05; **, p Ͻ 0.005 as analyzed by one-way analysis of variance followed by a least significant difference post hoc test. C, functional analysis of yeast retromer subunits. The indicated retromer yeast knock-out strains were treated as in A and B. The values are percentages of ALP activity/g protein as compared with the wild-type control strain and represent the averages Ϯ S.D. from three independent experiments. An ATG7 knock-out strain was used as a positive control. D, immunoblot analysis of p62 and GAPDH (loading control) from lysates of MCF7-eGFP-LC3 cells 64 h after transfection with indicated siRNAs, which had been validated by QRT-PCR (data not shown). The cells were left untreated or treated with 50 M etoposide or starved for amino acids for 16 or 8 h before harvesting, respectively. The experiment was repeated a minimum of three times with essentially same results. E, functional analysis of human candidate proteins. MCF7 cells stably expressing RLuc-LC3wt or RLuc-LC3 G120A were plated in separate wells of a 96-well dish and transfected with indicated siRNAs. In the upper graph, 50 nM EnduRen TM was added 17 h after the transfection, and the luciferase activities were measured at the indicated time points. For the lower graph, 50 nM EnduRen TM was added 54 h after transfection. Two hours later the luciferase activities were By analyzing autophagosomes from cells undergoing stress-induced autophagy by amino acid starvation, rapamycin-induced autophagy, and basal autophagy, which was blocked by concanamycin A, we compared autophagosomes from cells under fundamental different conditions. These analyses revealed proteins that always appear to associate with autophagosomes and others that appear to associate in a cell condition-dependent manner. Hence, stress-induced macroautophagy appears to differ from basal macroautophagy, and further studies are needed to understand stimulus-dependent cargo recruitment to autophagosomes.
An interesting finding in our study was the association between the proteasome and the autophagosome. In the past, the proteasomal and the autophagosomal/lysosomal pathways have been regarded as discrete degradation routes. Lately, a closer connection became more likely (4,24,44), and it has been described that the proteasome can be degraded in lysosomes (45). However, no detailed analyses of the interplay between these two degradation pathways have been reported. Here, we show that the proteasome localizes to autophagosomes regardless of the three stimuli and that the amount and activity of the proteasome are significantly decreased by functional autophagy. It is surprising that the proteasome associates with autophagosomes as the overall protein degradation increases during autophagy. However, several possible reasons might exist. First, autophagy induction is often associated with cellular stress and decreased cell growth. Thus, the removal of the proteasome and the subsequent accumulation of proteasomal substrates such as cyclins may contribute to cell cycle arrest. We tested this hypothesis but were not able to detect increased levels of cyclin D in cell lysate after autophagy induction (data not shown). This might be due to several compensatory mechanisms: 1) inhibition of protein synthesis because of inhibition of mTORC1 during autophagy activation or 2) autophagy might degrade proteins normally removed by the proteasome, e.g. via removal of p62-ubiquitinated protein aggregates (22). Second, one of the major functions of autophagy is to preserve cellular energy and thus degradation of the excess proteasome, which is not needed in circumstances where autophagy is activated, and de novo protein synthesis minimized via the inhibition of mTORC1 may serve as an important source of energy and amino acids. Furthermore, the degradation of proteasome substrates in autolysosomes may be less energy demanding than proteolysis by the proteasome. Last but not least, the proteasome could "help" to degrade ubiquitinated/unfolded cargo in the autophagosome. In favor of this theory, several subunits of PSMA (subunits 2, 3, 5, 6, and 7) and PSMB (subunits 1, 2, 3, 4, 5, 6 and 7) proteins were detected in the autophagosome cluster. Because these subunits form the "active proteasome," the functional proteasome complex appears to be autophagosomal cargo. This is further supported by the detection of protein-conjugated ubiquitin in autophagosomes by PCP-SILAC, immunofluoresence, and Western blot (data not shown) and by the fact that the proteasome activity is not decreased in concanamycin A-treated cells even though autophagosomes containing proteasomes accumulate in these cells.
Our data suggest that there is a balance between the two degradation systems that can be shifted in favor of one or the other, in our case active autophagy leading to reduced proteasome activity. In support of this notion, we observed increased levels of rapamycin-induced autophagy in yeast when the proteasome was inactivated by temperature-sensitive pre1-1 or pre2-1 mutants. The pre1-1/2-1 double mutant displayed synergistic effects (supplemental Fig. S7). Enhanced autophagy upon proteasome impairment suggests a compensatory mechanism as reported previously (21).
Among the new autophagosomal candidate proteins, we identified six autophagy regulators by a functional yeast screen that were not identified in the original screens characterizing Atg genes in yeast (46). This difference might be due to different screening methods. The original paper identified mutants defective in nitrogen starvation-induced autophagy as clones that lost viability faster than wild-type cells and were defective in accumulation of autophagic bodies in the vacuole. In contrast, we used the ALP assay. The orthologs of mammalian CAP1, EEF1G, RHEB, NP, and VPS35 were required for starvation-induced autophagy, whereas the yeast strain lacking the ortholog of mammalian GNB2L1 showed an increased autophagic response compared with the wild-type strain. Interestingly, testing additional subunits of the retromer complex revealed the same phenotype initially observed for the VPS35 subunit. In addition, the subunit VPS26 associated with autophagosomes in starved cells as identified by PCP-SILAC and anti-eGFP pulldowns. These data suggest that the retromer complex and retrograde transport play a role in autophagy, possibly by delivering vesicles or recycling autophagosomal membrane proteins.
Applying siRNA knockdown in human cells, we observed a change in the autophagic response after depleting cells for EEF1G or RHEB. Although we did not observe any major changes in the level of autophagy in cells depleted for CAP1, VPS35, and GNB2L1, we cannot exclude that these might have an effect if a more efficient down-regulation was measured (T ϭ 0), the cells were left untreated (control) or treated with 50 M etoposide, and the luciferase activities were measured at the indicated time points. The experiment was repeated five times with similar results. F, functional analysis of human candidate proteins. MCF7-eGFP-LC3 cells were left untreated or treated with 50 M etoposide for 6 h. Histograms with percentages of green cellular crosssections with over five LC3-positive dots are shown. The values represent the means Ϯ S.D. from three to six independent experiments. *, p value Ͻ 0.05 analyzed by a two-tailed unpaired t test. Ctrl, control; WT, wild type; MM, mock treated control. achieved. As expected, the mTORC1 activator RHEB constitutively blocked autophagy in human cells. The fact that we identify it as an autophagosomal protein fits with our recent data showing that starvation induces its degradation (19,47). Therefore, removal of RHEB by autophagy might serve as a positive feedback mechanism to ensure a massive autophagy response. Surprisingly, the ortholog of RHEB in yeast was found to be required for both starvation-and rapamycininduced autophagy. This controversy could be explained by the fact that contrary to mammals, which have two RHEB isoforms, yeast only has one RHEB isoform (RHB), which is only distantly related to mammalian RHEBs and whose knock-out does not affect the activity of TOR (48). Interestingly, EEF1G was found to be a universal autophagy regulator. Similar to LC3 it locates in/at the autophagosomes and is required for the process. However, contrary to LC3, it functions upstream of mTORC1. The exact molecular mechanism whereby EEF1G regulates mTORC1 activity remains to be revealed.
In summary, the analysis of organellar protein composition in our study highlights the complex character of autophagy shedding light on an elemental metabolic pathway that we are only starting to comprehend. The catalog of autophagosomerelated proteins assembled from PCP-SILAC data and crossexamined by autophagosome immunoprecipitations revealed interplay between the autophagosome and the proteasome systems. In addition, we identified novel autophagosomal components of which EEF1G was found to be an autophagy regulator in both yeast and human cells. ‡ ‡ To whom correspondence should be addressed. E-mail: joern.dengjel@frias.uni-freiburg.de; mj@cancer.dk, jens.andersen@ bmb.sdu.dk.