Abstract

The large-scale turnover of intracellular material including organelles is achieved by autophagy-mediated degradation in lysosomes. Initiation of autophagy is controlled by a protein kinase complex consisting of an Atg1-family kinase, Atg13, FIP200/Atg17, and the metazoan-specific subunit Atg101. Here we show that loss of Atg101 impairs both starvation-induced and basal autophagy in Drosophila. This leads to accumulation of protein aggregates containing the selective autophagy cargo ref(2)P/p62. Mapping experiments suggest that Atg101 binds to the N-terminal HORMA domain of Atg13 and may also interact with two unstructured regions of Atg1. Another HORMA domain-containing protein, Mad2, forms a conformational homodimer. We show that Drosophila Atg101 also dimerizes, and it is predicted to fold into a HORMA domain. Atg101 interacts with ref(2)P as well, similar to Atg13, Atg8a, Atg16, Atg18, Keap1, and RagC, a known regulator of Tor kinase which coordinates cell growth and autophagy. These results raise the possibility that the interactions and dimerization of the putative HORMA domain protein Atg101 play critical roles in starvation-induced autophagy and proteostasis, by promoting the formation of protein aggregate-containing autophagosomes.

1. Introduction

Autophagy ensures the lysosome-mediated degradation and recycling of cytoplasmic components including organelles. During the main pathway, a phagophore cistern (also called an isolation membrane) forms and captures cargo destined for breakdown into a double-membrane autophagosome. This vesicle then fuses with a late endosome or lysosome containing acidic hydrolases. Assembly of the phagophore is achieved by the action of Atg proteins. Atg genes were originally discovered in yeast, and the majority of them have clear orthologs in higher eukaryotes including animals [1, 2].

Initiation of autophagy usually begins with the activation of an Atg1 protein kinase complex in animal cells, which contains the serine/threonine kinase Atg1 (its orthologs are called UNC51 in worms and ULK1 and 2 in mammals), Atg13, FIP200/Atg17, and the metazoan-specific subunit Atg101 [3]. This complex directly binds to Tor (target of rapamycin) kinase, which is active when bound to digesting lysosomes and promotes cell growth and inhibits autophagy by phosphorylating Atg1 [4, 5]. Autophagy-inducing stimuli such as starvation rapidly lead to inactivation of Tor and induction of Atg1-dependent autophagy [3, 6]. This results in the removal of inhibitory phosphogroups from Atg1 by poorly characterized phosphatases, which may potentially include PP2A [7]. Atg1 then undergoes autophosphorylation on residues separate from those phosphorylated by Tor and also phosphorylates downstream targets including Atg13 [3, 8, 9]. Overexpression of Atg1 strongly promotes autophagy in Drosopohila, while expression of a kinase dead form suppresses starvation-induced autophagy [10].

It is largely unknown how the Atg1 kinase complex transmits its autophagy-inducing signal to downstream components, which include an autophagy-specific lipid kinase complex, phospholipid effectors such as Atg18, the transmembrane protein Atg9, and two protein conjugation systems that include the ubiquitin-like Atg8 family proteins [1, 3]. Atg8 is conjugated to phosphatidylethanolamine, and thus it is bound to phagophore and autophagosome membranes [11, 12]. Atg8, and its mammalian homologs such as LC3, binds to cargo receptors including p62 to mediate the selective autophagic breakdown of ubiquitinated protein aggregates [13, 14]. Ref(2)P, the fly homolog of p62, is required for the formation of these protein aggregates in Drosopohila [15].

FIP200/Atg17 is thought to act as a scaffold protein in the Atg1 kinase complex, whereas the role of the metazoan-specific subunit Atg101 is poorly characterized [3]. Atg101 (also known as C12orf44) is a subunit of the Atg1/ULK complex in human cells, and shRNA depletion of Atg101 impairs autophagy in mammalian cells [16, 17]. Moreover, the Atg101 ortholog epg-9 has recently been shown to be required for degradation of P granule aggregates in worm embryos as well [18]. As the precise role of Atg101 in autophagy is unclear, we decided to analyze its function in Drosopohila.

2. Materials and Methods

2.1. Drosopohila Genetics

Flies were reared on standard yeast-cornmeal-agar medium (fed), and mid-L3 stage larvae were transferred to a 20% sucrose solution for 3 h for starvation experiments. Fat body cell clones expressing RNAi constructs Atg101[KK106176] and Atg101[HMS01349] were generated spontaneously, as described in detail previously [1922].

2.2. Histology

Dissected larval carcasses were incubated in a solution of LysoTracker Red (LTR) and DAPI as before [1923]. Fixed samples were processed for indirect immunofluorescence using rat anti-Atg8a [9] and rabbit anti-ref(2)P [22] and imaged as described in detail previously [19, 20, 22]. Dissected fat body lobes containing Atg101 RNAi clones were attached to poly-L-lysine-coated coverslips in PBS, photographed live to record the position of GFP-positive cells, fixed and embedded into Durcupan (Fluka), and then sectioned and analyzed by transmission electron microscopy as described [19, 21]. Statistical analysis was carried out as described [19, 20]. N refers to the number of animals, and multiple cells were evaluated from each animal.

2.3. RT-PCR

Total RNA was isolated from starved L3 stage larvae of the genotypes w (used as wild type), Atg101[KK106176]/+, Act-Gal4/+, and Atg101[HMS01349]/Act-Gal4 using PureLink RNA Mini Kit (Invitrogen), followed by preparation of cDNAs using RevertAid First Strand cDNA Synthesis Kit (Thermo) and RNase free DNase I (Sigma). PCR reactions on cDNA samples were performed with the following primers: TACTCCTCCCGACACAAAGC, CTGGGTCATCTTCTCACGGT for Actin5C (23 amplification cycles), and ATGAACGCGCGTTCGCAG, TCACATTGCGAGCGTTTCCT for Atg101 (25 amplification cycles). Parallel reactions were performed without adding reverse transcriptase to control for DNA contamination. As expected, no PCR products were obtained in these experiments.

2.4. Cell Culture and Immunoprecipitations

Drosopohila D.Mel-2 cells (Invitrogen) were used for transfection and immunoprecipitation experiments as recently described elsewhere [9, 19, 20]. The following constructs were used in coimmunoprecipitation experiments: an N-terminally truncated 3xFLAG-ref(2)P lacking the PB1 domain mediating self-aggregation, 3xHA-Keap1 [20], 3xHA-Atg18 [24], 3xHA-Atg101, full-length 3xHA-Atg13 and its N-terminal, middle, and C-terminal fragments [9], and kinase dead myc-Atg1 [8]. 3xHA-Atg8a, 3xHA-Atg16, 3xHA-RagC, 3xFLAG-Atg101, and Atg101-3xFLAG were generated by PCR amplifying the full-length coding sequences from a cDNA sample and cloning these into appropriate UAS vectors, respectively. GST coding sequence was PCR amplified and cloned into a UAS-3xHA vector downstream of the HA tag. 3xHA-GST-Atg1 constructs were generated by PCR amplifying and cloning the appropriate Atg1 fragments downstream of GST.

2.5. Bioinformatics

Sequence analysis and motif search were performed using the NCBI BLAST service as well as the Pfam (version 27.0) and Prosite (release 20.97) databases and their associated search tools. Alignments of orthologous sequences were extracted from OMA groups as cross-linked to the corresponding SwissProt entries [25]. For structure prediction, the I-TASSER and Phyre2 servers were used [26, 27]. Pairwise and multiple structural alignments were done with DaliLite and MAMMOTH-Mult, respectively [28, 29]. The ANCHOR and IUPred servers were used to predict unstructured regions and potential binding sites within these regions in fly Atg1, Atg13, and Atg101 [30, 31]. Structure visualization was performed using UCSF Chimera [32] and the Prosite HORMA logo was reproduced with WebLogo [33]. The alignment of human, fly, and worm Atg101 was generated using CLUSTALW and colored with TEXSHADE at the Biology Workbench server [34].

3. Results and Discussion

3.1. Drosopohila Atg101 Is Required for Starvation-Induced Autophagy

BLAST searches reveal that Drosopohila Atg101 shares 51% amino acid identity (111/218) and 71% similarity (156/218) with human Atg101 and 33% amino acid identity (81/244) and 46% similarity (113/244) with isoform a of worm Atg101/EPG-9, respectively (Figure 1), suggesting that it is an ortholog of Atg101 proteins. Two separate transgenic RNAi lines are available from public stock centers that allow inducible silencing of Drosopohila Atg101, both integrated into specific, nonrandom landing sites in the genome. Atg101[KK] is based on a long hairpin, while Atg101[HMS] contains a short, microRNA-based duplex RNA that targets a 22-nucleotide sequence in the 3′ untranslated region of the endogenous mRNA. LysoTracker Red (LTR) is widely used for the labeling of acidic autolysosomes in the larval Drosophila fat body, as this vital dye shows little to no punctate staining in fat body cells of well-fed larvae [9, 19, 23, 35]. Expression of either of the Atg101 RNAi constructs in GFP-marked cell clones prevents starvation-induced LysoTracker Red- (LTR-) positive autolysosome formation (Figures 2(a), 2(b), and 2(c)). As expected, systemic expression of either of these transgenic RNAi lines strongly decreases endogenous Atg101 mRNA levels (Figure 2(d)). Depletion of Atg101 also impairs the distribution of endogenous Atg8a, as instead of the numerous small Atg8a-positive autophagosomes seen in control cells, fewer but larger Atg8a structures are observed in GFP-marked Atg101 knockdown cells (Figures 3(a), 3(b), and 3(d)). This may not be simply due to incomplete gene silencing, as these aberrant Atg8a-positive structures are also seen upon simultaneous expression of both Atg101 RNAi constructs (Figures 3(c) and 3(d)). This phenotype is very similar to the accumulation of Atg8/LGG-1 aggregates reported in Atg101/epg-9 mutant worms [18]. The enlarged Atg8a structures observed in Atg101 RNAi cells of starved animals colocalize with protein aggregates containing the specific autophagy cargo ref(2)P (Figure 3(e)). Although it is not possible to conclude that Atg101 is not required for Atg8a recruitment to ref(2)P aggregates solely based on RNAi experiments, this seems to be a possibility, which is also supported by data from mutant worms. Interestingly, of the core Atg genes, Atg2 has also been shown to be dispensable for Atg8a recruitment to the phagophore assembly site and protein aggregates in worm, fly, and mammalian cells [19, 24, 36, 37]. To exclude the possibility that these Atg8a dots might represent autophagosomes, we carried out electron microscopy. Ultrastructural analysis indeed revealed that autophagosome formation is blocked in Atg101 RNAi cells (Figure 4). A potential explanation for the presence of a few enlarged Atg8a-positive dots in Atg101 loss-of-function cells but not in FIP200 null mutants [9] may be that mammalian Atg8 homologs have been shown to directly bind to FIP200, Atg13, and Atg1/ULKs [38]. These Atg1 kinase subunits may facilitate the recruitment of Atg8a to aggregates containing ref(2)P and ubiquitinated proteins, which have been proposed to act as scaffolds for autophagosome biogenesis, both in mammals and Drosophila [9, 24, 39].

3.2. Drosopohila Atg101 Is Required for Selective Basal Autophagy of Ref(2)P-Containing Protein Aggregates

Ref(2)P appears to be selectively degraded by autophagy all the time, independent of the induction conditions such as starvation or developmental contexts, similar to its mammalian homolog p62 [24, 40]. Continuous basal autophagy is usually difficult to be visualized directly, due to its low levels. For this reason, assessing basal autophagic activity by looking at levels of ref(2)P has become a standard test, similar to mammals [14, 19, 22, 41, 42]. Expression of either of the two Atg101 silencing constructs results in large-scale accumulation of ref(2)P aggregates in GFP-marked cells, when compared to surrounding wild-type tissue in well-fed larvae (Figures 5(a), 5(b), and 5(c)), indicating defects in selective basal autophagy.

3.3. Atg101 Interacts with Atg1

The Atg13, FIP200, and ULK1 subunits of the mammalian Atg1/ULK kinase complex have been shown to bind to Atg101 in coimmunoprecipitation experiments on the level of endogenous proteins [16]. Atg101 directly binds to Atg13 in mammals, and the Atg101 ortholog EPG-9 was recently proposed to directly interact with the Atg1 ortholog UNC-51 as well in C. elegans [1618]. In line with these studies, we found that myc-tagged, kinase dead Atg1 coimmunoprecipitates with Atg101 in cultured Drosopohila cells (Figure 6(a)). We used the kinase dead form because wild-type Atg1 is expressed poorly in cultured cells, and it also reduces the expression of other transgenes, likely due to Atg1-mediated feedback inhibition of Tor-dependent translation [9, 10, 43]. To determine which regions of the 835 amino acid long Atg1 protein are involved in its interaction with Atg101, we generated four HA-GST-tagged constructs for different Atg1 fragments. The N-terminal kinase domain was not included in these experiments, as again it is expressed poorly and strongly impairs the expression of coexpressed constructs [9]. We decided to search for potential Atg101 binding sites in the middle region (amino acids 233–725) of Atg1, because the N-terminal kinase domain is the only known domain in this protein, and the C-terminal region is involved in its binding to Atg13 [44]. The rest of the protein is predicted to be mostly unstructured by the IUPRED server, and the ANCHOR method identifies numerous potential binding sites within these unstructured regions. We found that Atg101-FLAG coimmunoprecipitates with two Atg1 fragments containing either amino acids 233–360 or 565–725 but not with fragments containing amino acids 364–459 or 461–570 (Figures 6(b) and 6(c)). These results suggest that Atg101 interacts with multiple regions in Atg1 either directly or indirectly.

3.4. Atg101 Binds to the HORMA Domain of Atg13, and It Can Dimerize

Coimmunoprecipitations showed that full-length HA-Atg13 binds to Atg101-FLAG (Figure 7(b)), in line with data from mammals [16, 17]. Mapping experiments revealed that only those regions of Atg13 coprecipitate with Atg101 that contain amino acids 1–230, whereas the unstructured middle and C-terminal regions do not show binding (Figures 7(a) and 7(b)). The N-terminal part of Atg13 folds into a HORMA (Hop1, Rev1, and Mad2) domain similar to that of Mad2, a spindle checkpoint protein [45]. Mad2 forms a conformational homodimer, as its dimerization requires the binding of two different stable conformations to each other: open, O-Mad2, and closed, C-Mad2 [46]. Since the Atg13 HORMA domain is similar to C-Mad2 [45], we hypothesized that Atg101 may potentially be a HORMA domain protein as well, perhaps capable of forming either an O-Mad2 state or both O-Mad2 and C-Mad2 states. To study this model further, we first determined the interaction of Atg101 with itself. We found that Atg101-FLAG coprecipitates with HA-Atg101, and vice versa HA-Atg101 coprecipitates with Atg101 tagged by FLAG on either its N- or C-terminus (Figure 7(c)), suggesting that Atg101 dimerizes.

3.5. Atg101 Is Potentially a HORMA Domain Protein

Atg101 is classified as a single-domain protein with a domain of unknown function (DUF1649 in Pfam), with no trivially detectable relatives of known structure. Basic structural predictions suggest that it is a structured globular cytoplasmic protein. Results obtained with I-TASSER and Phyre2 indicate that the protein might have a HORMA domain fold. Both servers identified templates with HORMA fold and three of I-TASSER’s 5 predicted structures have a fold similar to each other as well as the Atg13 and Mad2 HORMA domains (Figure 8). Moreover, the model predicted by Phyre2 also has a HORMA-like structural core (not shown). To assess the plausibility of this hypothesis, we created a multiple structure alignment of the best I-TASSER model of Atg101 with the known structures of human Mad2 [46] and yeast Atg13 [45], and manually compared it to the HORMA domain profile composed of 33 proteins in Prosite (entry PS50815). Our analysis reveals a number of corresponding sites showing conservation both in the HORMA domains in Prosite and within Atg101 homologs (OMA group 413828 containing Atg101 proteins from 66 different species). In addition, a characteristic pattern of largely conserved Leu/Ile/Val residues in Atg101 proteins and matching hydrophobic residues in the Prosite profile is apparent (Figure 9).

3.6. Atg101 Interacts with Ref(2)P/p62

Drosopohila Atg101 was suggested to bind to ref(2)P and also to several subunits of the lysosomal proton pump v-ATPase complex (Vha26, Vha36-1, Vha44, Vha55, and Vha68-2) in a large-scale proteomic study [47]. As we have recently found that the Atg1 kinase subunit FIP200 frequently localizes to ref(2)P aggregates near lysosomes [9], we wished to confirm the proposed interaction of Atg101 with ref(2)P. Indeed, FLAG-tagged ref(2)P coprecipitated with HA-Atg101 (Figure 10(a)). We used HA-Atg8a as a positive control in this experiment, as Atg8 homologs are established binding partners of p62 in mammalian cells, and fly Atg8a has been proposed to bind to ref(2)P through a conserved Atg8 interacting region as well [14, 15, 48]. These data raised the possibility that the entire Atg1 kinase complex may bind to aggregates of ref(2)P and ubiquitinated proteins. Indeed, we found that Atg13 also shows an interaction with ref(2)P, similar to Atg101 (Figure 10(b)). Positive controls that we used in these coimmunoprecipitations include Keap1 and Atg18, both of which have recently been shown to bind to ref(2)P [20, 24]. Ref(2)P coprecipitates with additional critical regulators of autophagy as well: the ubiquitin-binding protein Atg16 [49] and the small GTPase RagC, a known Vha complex-associated regulator of Tor kinase [4], which shows a particularly strong interaction with ref(2)P as inferred from a large-scale proteomic study [47] (Figure 10(b)). These data are in line with the potential role of lysosome-associated protein aggregates as scaffolds for autophagosome biogenesis [9, 24, 39].

4. Conclusions

Our data establish that Atg101 is required for autophagy in Drosopohila, and that it is a subunit of the Atg1 kinase complex, similar to mammals and worms [1618]. The potential binding of Atg101 to Atg1 may involve multiple regions in Atg1, whereas the HORMA domain of Atg13 appears to be the only region important for its interaction with Atg101. Atg101 is a small protein of 218 amino acids, and it is predicted to fold into a globular domain, potentially a HORMA domain. Proteins with such structure often have numerous binding partners, the identity of which determines the open or closed conformation in the case of Mad2 [46]. Our bioinformatics analysis predicts a possible HORMA domain structure for Atg101, and we show that it can indeed dimerize. If Atg101 functions similar to Mad2, then it may perhaps be capable of switching between open and closed conformations and form a conformational homodimer as well. This would represent an elegant way of determining which proteins can bind to Atg101 at one time. This model predicts that Atg101 binds to Atg13 through direct interaction of the two HORMA domains, whereas for some of its other partners (potentially including Atg1), certain oligopeptides located in unstructured regions may engage in an interaction with Atg101, as seen in the case of Mad2 as well [46].

It is important to emphasize that the HORMA domain structure of Atg101 is not supported experimentally. Further analysis of Atg101 binding partners, fine mapping of the binding sites mediating the interaction of each partner with Atg101, and analysis of the Atg101 domain by structural biochemistry will be necessary to test our hypothetical model presented in this paper. This will be particularly important, as the absence of an Atg101 ortholog in yeast clearly indicates that a metazoan-specific regulatory mechanism has evolved to control the initiation step of autophagy in animal cells.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank Sarolta Pálfia for her excellent technical assistance, László Nyitray, Viktor Kis, and Rita Pancsa for helpful discussions, and stock centers BDSC and VDRC for providing Atg101 RNAi lines. This work was supported by the Wellcome Trust (087518/Z/08/Z), the Hungarian Scientific Research Fund (OTKA K83509), the European Union (TÁMOP-4.2.1.B-11/2/KMR-2011-0002), and the Hungarian Academy of Sciences (BO/00552/11).