Next Article in Journal
Expression Profiling Reveals the Possible Involvement of the Ubiquitin–Proteasome Pathway in Abiotic Stress Regulation in Gracilariopsis lemaneiformis
Next Article in Special Issue
Identification and Potential Participation of Lipases in Autophagic Body Degradation in Embryonic Axes of Lupin (Lupinus spp.) Germinating Seeds
Previous Article in Journal
Synthesis, Properties, and Biomedical Application of Dicationic Gemini Surfactants with Dodecane Spacer and Carbamate Fragments
Previous Article in Special Issue
Sugar Starvation Disrupts Lipid Breakdown by Inducing Autophagy in Embryonic Axes of Lupin (Lupinus spp.) Germinating Seeds
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interactome of Arabidopsis ATG5 Suggests Functions beyond Autophagy

by
Pernilla H. Elander
1,†,
Sanjana Holla
1,†,
Igor Sabljić
1,
Emilio Gutierrez-Beltran
2,3,
Patrick Willems
4,5,
Peter V. Bozhkov
1 and
Elena A. Minina
1,*
1
Department of Molecular Sciences, Uppsala BioCenter, Swedish University of Agricultural Sciences and Linnean Center for Plant Biology, 750-07 Uppsala, Sweden
2
Instituto de Bioquımica Vegetal y Fotosıntesis, Universidad de Sevilla and Consejo Superior de Investigaciones Cientıficas, 41092 Sevilla, Spain
3
Departamento de Bioquimica Vegetal y Biologia Molecular, Facultad de Biologia, Universidad de Sevilla, 41012 Sevilla, Spain
4
Department of Plant Biotechnology and Bioinformatics, Ghent University, 9052 Ghent, Belgium
5
Department of Biomolecular Medicine, Ghent University, 9000 Ghent, Belgium
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2023, 24(15), 12300; https://doi.org/10.3390/ijms241512300
Submission received: 5 July 2023 / Revised: 27 July 2023 / Accepted: 28 July 2023 / Published: 1 August 2023
(This article belongs to the Special Issue New Insight into Signaling and Autophagy in Plants 3.0)

Abstract

:
Autophagy is a catabolic pathway capable of degrading cellular components ranging from individual molecules to organelles. Autophagy helps cells cope with stress by removing superfluous or hazardous material. In a previous work, we demonstrated that transcriptional upregulation of two autophagy-related genes, ATG5 and ATG7, in Arabidopsis thaliana positively affected agronomically important traits: biomass, seed yield, tolerance to pathogens and oxidative stress. Although the occurrence of these traits correlated with enhanced autophagic activity, it is possible that autophagy-independent roles of ATG5 and ATG7 also contributed to the phenotypes. In this study, we employed affinity purification and LC-MS/MS to identify the interactome of wild-type ATG5 and its autophagy-inactive substitution mutant, ATG5K128R Here we present the first interactome of plant ATG5, encompassing not only known autophagy regulators but also stress-response factors, components of the ubiquitin-proteasome system, proteins involved in endomembrane trafficking, and potential partners of the nuclear fraction of ATG5. Furthermore, we discovered post-translational modifications, such as phosphorylation and acetylation present on ATG5 complex components that are likely to play regulatory functions. These results strongly indicate that plant ATG5 complex proteins have roles beyond autophagy itself, opening avenues for further investigations on the complex roles of autophagy in plant growth and stress responses.

1. Introduction

Plants must endure a range of unfavorable environmental conditions to survive and propagate. Heat, drought, water logging, salinity, and pests are but a few examples of environmental stresses, which plants must cope with. By utilizing autophagy, plant cells can dispose of harmful cellular constituents and recycle the material for other purposes. During macroautophagy (hereafter autophagy), cellular constituents are engulfed by a de novo formed double-membrane vesicle called autophagosome and delivered to the lytic vacuole for degradation [1]. The autophagy pathway is regulated by approximately 40 autophagy-related (ATG) proteins, which include about 20 core ATGs generally grouped into five large complexes/systems essential for autophagic activity [2]. Central to the autophagic pathway is the lipidation of ATG8, which comprises a series of reactions enabled by ATG4, ATG7, ATG3, ATG10 and ATG5-ATG12/ATG16 complex culminating in the conjugation of ATG8 with phosphatidylethanolamine (PE). ATG8-PE is anchored on the membrane of the forming autophagosome facilitating the membrane elongation, cargo recognition and autophagosome trafficking [3,4].
In our previous study [5], we upregulated individually ATG5 and ATG7, encoding for two core autophagy proteins required for the lipidation of ATG8. We discovered that upregulation of either of these two genes positively affected autophagic activity without impacting the transcription of the other genes participating in the ATG8-lipidation, indicating that ATG5 and ATG7 are rate-limiting components of the autophagic pathway. Interestingly, overexpression (OE) of these genes promoted plant growth, seed yield and longevity. Furthermore, OE plants developed more inflorescences and exhibited prolonged flowering time, thus producing more seeds. Furthermore, ATG5 and ATG7 OE had a delayed onset of leaf senescence compared to wild-type (WT) plants and improved resistance to necrotrophic pathogens and oxidative stress [5].
Even though the observed phenotypic traits correlated with upregulated autophagic activity in the OE plants, we still could not decidedly claim that these traits were caused solely by upregulated autophagy. The question remained whether ATG5 and ATG7 might also act in other pathways.
Indeed, there are several examples of ATG proteins participating in autophagy-unrelated pathways. For instance, ATG8 orthologs in animal cells partake in intracellular trafficking and Golgi transport [6]. A plant ATG8 ortholog was shown to interact with ABNORMAL SHOOT3 at the late endosome to promote senescence by protein degradation, this interaction did not require ATG8 conjugation with PE [7]. Furthermore, a product of ATG5 proteolytic cleavage was shown to act as an apoptotic effector in animal cells [8].
Arabidopsis knockout (KO) mutants of core ATG genes show no discernible phenotypes at early developmental stages, when grown under standard conditions [5]. However, they undergo early senescence, display increased stress-susceptibility, and compromised immunity to necrotrophs compared to WT plants [5]. Interestingly, KO of different ATG genes display a range of phenotypes under autophagy-inducing stress conditions. For example, while atg5 and atg7 plants have similar phenotypes [5], atg2 plants display more severe senescence and growth stagnation symptoms [9], whereas atg9 plants have less severe symptoms [10]. Since autophagy is abrogated in all these mutants, the difference in the phenotypes might stem from autophagy-unrelated functions of the encoded ATG proteins, as we hypothesised for ATG5 and ATG7 OE plants.
To shed light on autophagy-unrelated functions of the ATG proteins, we began with identifying ATG5 interactome under standard growth conditions that do not boost autophagy. ATG5 is known to form complex with ATG12 and/or ATG16 [11,12]. In Arabidopsis cells most of ATG5 is covalently linked to ATG12 [5] via a C-terminal glycine residue of ATG12 bound to a side chain of a lysine residue of ATG5. An AlphaFold2-generated structure prediction for Arabidopsis ATG5 [13,14], revealed its structural similarities with the better characterized ortholog of ATG5 from yeast [11], enabling a prediction of Lys128 (K128) as the lysine residue binding to ATG12.
ATG16 in turn, binds non-covalently to ATG5 to form an ATG5-ATG12/ATG16 complex, which possesses E3-like ligase activity and is further referred to as “complete complex”. ATG16 recruits the complete complex to the phagophore, where its E3-like ligase activity is implemented to conjugate ATG8 to the resident PE [15]. Interestingly, ATG5 and ATG16 can interact in the absence of ATG12 to form a different complex [12], which in this study is named “minimal complex”.
Here, we identified proteins pulled-down with two types of bait: WT ATG5 capable of forming a complete ATG5-ATG12/ATG16 complex that is required for autophagy and the mutant ATG5K128R that can form only the minimal ATG5-ATG16 complex (Figure 1A). We demonstrated only partial overlap between autophagy-competent ATG5 and autophagy-incompetent ATG5K128R interactomes. This functional network unveils previously unknown interacting partners of plant ATG5 and suggests its autophagy-unrelated functions, thereby providing a new insight into diverse roles ATGs play in plant growth, development, and stress responses.

2. Results

2.1. Generation and Characterization of Arabidopsis ATG5-TAP Lines

To identify interactors of both complete (ATG5-ATG12/ATG16) and minimal (ATG5-ATG16) complexes (Figure 1A), we engineered genetic constructs encoding for Tandem Affinity Purification (TAP) tag [16] fused with the C-terminus of ATG5, ATG5K128R mutant or GFP, placed under the control of 2x35S promoter or ATG5 native promoter. We predicted that ATG5K128R mutation impairs the conjugation of Arabidopsis ATG5 to ATG12 at the early steps of the autophagy pathway (Figure 1A).
The resulting constructs were introduced into Arabidopsis Columbia-0 (Col-0) WT or autophagy deficient atg5-1 background. Two transgenic lines per construct showing expression detectable by immunoblot were selected for further studies. To confirm that ATG5K128R mutation indeed impairs conjugation with ATG12, formation of the complete complex and thus abrogates autophagy, we firstly performed immunodetection of the TAP tag in the transgenic lines expressing WT form of ATG5 or its K128R mutant. Indeed, the band corresponding to ATG5-TAP-ATG12 conjugate was not detectable in the protein extracts from plants expressing ATG5K128R (Figure 1B).
We further verified that the K128R mutant is not able to restore autophagy in ATG5-deficient cells. For this, we isolated protoplasts from atg5-1 plants expressing GFP-ATG8a marker for autophagosomes. Protoplasts were transformed with plasmids encoding TagRFP- or TAP-tagged ATG5 or ATG5K128R and treated with AZD8055 and Concanamycin A (ConA) to induce autophagy and block degradation of autophagic bodies in the vacuole [17]. Upon induction of autophagy, GFP-positive autophagic bodies accumulated in the vacuoles in the presence of ATG5, but not in the presence of ATG5K128R, corroborating the inability of this K128R mutant to form the functional complete complex required for autophagy (Figure 1C).
Finally, we assessed phenotypes of atg5-1 plants complemented with ATG5-TAP or ATG5K128R-TAP driven by the native ATG5 promoter. In contrast to WT and ATG5-TAP complemented line, the ATG5K128R-TAP complemented lines failed to rescue autophagy deficiency phenotype manifested in early senescence, smaller rosette leaves and reduced number of inflorescences (Figure 1D).

2.2. Interactomes of ATG5 and Its Complete and Minimal Complexes

2.2.1. Affinity Purification

The aforementioned Arabidopsis lines expressing ATG5-TAP, ATG5K128R-TAP or GFP-TAP under control of 2x35S promoter in the WT (Col-0) background were used for the pull-down assay. Tandem affinity purification was performed as described earlier [18] on the protein extracts obtained from fully expanded leaves harvested at the bolting stage of plants grown under standard long day conditions. This developmental stage was previously shown to correlate with a low basal autophagy level [5] and was chosen to enrich for autophagy-independent interactors of ATG5. Pull-down for ATG5 and ATG5K128R was performed in three biological replicates, each comprising a pool of more than 15 plants. The pull-down for the negative control (GFP-TAP) was performed in two biological replicates sampled in an identical manner. Purified proteins were analysed using liquid chromatography-tandem mass spectrometry (LC-MS/MS) and resulting data was analyses using MaxQuant and MSqRobSum [19].

2.2.2. Detection of Post-Translational Modifications

ATG proteins are constitutively expressed and undergo various post-translational modifications (PTMs) allowing rapid adjustment of their activity during switches between normal and stress conditions. For instance, inhibitory phosphorylation of mammalian ATG5 and ATG12 by ATG1 kinase governs spatio-temporal control of ATG8 lipidation during phagophore expansion [20], while inhibitory acetylation of mammalian ATG5, ATG7, ATG8 and ATG12 by the p300 acetyltransferase [21] is reversed by deacetylase Sirt1 during autophagy induction [22].
The role of PTMs in plant autophagy regulation is scarce and even less is known about the role of PTMs in the autophagy-unrelated functions of the plant ATG proteins. Therefore, we searched our LC-MS/MS data for the PTMs reliably detectable by mass spectrometry and also relevant for autophagy regulation, i.e., phosphorylation and acetylation, and detected PTMs of ATG5 and ATG12 (Figure 2A, Tables S1–S3). Interestingly, phosphorylation of Ser187 was identified only in ATG5, while acetylation of Lys183 was detectable on both ATG5 and ATG5K128R indicating that the latter modification might also play a role in autophagy-independent function of ATG5. Additionally, phosphorylation of Ser7 was present at detectable levels on the ATG12A protein, indicating potential conservation of inhibitory phosphorylation of plant ATG12 under nutrient-rich conditions.

2.2.3. Criteria for Identifying ATG5 Interactomes

In order to identify specific and common interactors of the minimal and the complete complexes, we selected proteins enriched in ATG5K128R and/or in ATG5 samples compared to the GFP-control. Significant interactors of ATG5/ATG5K128R were identified using MSqRobSum (see Methods), requiring a minimal fold change of 1.5 (compared to GFP control) and p value < 0.05. After filtering, we retained 85 interactors only present in the ATG5 pull-down, 82 only in the ATG5K128R pull-down and 88 common interactors in both ATG5 and ATG5K128R pull-downs (Figure 2B). From here on, we grouped potential interactors into three categories: proteins pulled-down with ATG5 only (complete complex-specific), with ATG5K128R only (potential interactors of individual ATG5) and with both ATG5 and ATG5K128R (minimal complex-specific) (Figure 2C,D).
In accordance with conditions used for the affinity purification experiment, the typical autophagy-related interactors of ATG5, including ATG3, ATG10 and ATG8 [12], were not detectable in our samples, confirming that the identified candidate interactors were likely enriched for autophagy-unrelated pathways. Expectedly, ATG12A and ATG12B were only detected in the interactome of the complete complex, whereas ATG16 was detected in both interactomes (Figure 2C,D).

2.2.4. Putative Interactors Acting as Stress Sensors

Interestingly, we identified PP2A-A3, the scaffolding subunit of the Protein Phosphatase 2A (PP2A), as an interactor of the minimal complex, while a catalytic subunit of PP2A (PP2A-3) was pulled down with autophagy-incompetent ATG5K128R. In mammalian cells PP2A is known to be an important regulator of the above mentioned PTMs, namely it dephosphorylates ATG13 upon inhibition of the TORC1 kinase complex enabling formation of ATG1-ATG13 complex and induction of autophagy [23]. Furthermore, different subunit compositions of the heterotrimeric PP2A have been shown to play alternative roles in animal autophagy [24]. Although plant PP2A is known to be important for abiotic and biotic stress response [25], its role in plant autophagy is still unknown and our observations provide the first insight on its potential participation and indicate a possible implication of autophagy-unrelated crosstalk between ATG5 and PP2A.
In addition to PP2A, we identified hexokinase 1 (HXK1, Table S3) as a potential interactor of the minimal complex. HXK1 was previously shown to suppress plant autophagy under glucose-rich conditions via an unknown mechanism [26]. Our finding suggests that HXK1 might not be acting via TORC1 as suggested previously but rather through the direct interaction with ATG5 and/or ATG16.

2.2.5. Putative Interactors belonging to Endomembrane Trafficking System

Autophagy is an integral part of the endomembrane trafficking system [27,28] and so far, ATG8 has been the best characterized molecular link between plant autophagosomal structures and other components of the endomembrane trafficking system [29]. Remarkably, we discovered a large set of proteins playing a role in endomembrane trafficking that interact with ATG5 and its complete and minimal complexes: proteins involved in COPII-mediated endoplasmic reticulum (ER) to Golgi transport (SAR1B, SAR1C, SEC24C, TRAPPC3), phagophore formation (TRAPPC11, SH3P2), ER-PM contact sites (SYT1, TPLATE subunit TML), retromer complex subunits (SNX1 and SNX2a) and a CLASP protein involved in membrane loading of SNX1 (Figure 2C,D, Table S3). Most interestingly, our results showed that homologous proteins encoded by multi-member gene families, with not yet fully resolved redundancy, show specificity towards either the WT or autophagy incompetent ATG5 bait. For example, SAR1C was detectable in the ATG5K128R pull-down only, unlike SAR1B that was found in both. Similarly, SNX1 was detected only in ATG5 pull-down samples, while SNX2a only in ATG5K128R (Figure 2C,D, Table S3).

2.2.6. Putative Interactors belonging to Ubiquitin-Proteasome System

Autophagy is tightly interlinked with another catabolic pathway, governed by the ubiquitin-proteasome system (UPS) [30,31]. Furthermore, mammalian ATG5 was shown to directly interact with UPS components to aid mitophagy [32]. Our affinity purification assay allowed detection of Ubiquitin Fusion Degradation 1 (UFD), a component of the CDC48 complex that was previously suggested to crosstalk with autophagy to help maintenance of chloroplastic proteins during oxidative stress [33] (Table S3). Furthermore, surprisingly, we discovered selective interaction between 20S proteasome subunit alpha E2 and complete ATG5 complex, while subunits alpha D1 and beta G1 showed preference towards the minimal complex (Table S3).

2.2.7. Putative Interactors of Nuclear-Localized ATG5 and ATG12

Interestingly, we also identified two components of nuclear pore complex as potential interactors of individual ATG5: Nucleoporin 155 (NUP155) and SEC13B (Figure 2D, Table S3) [34]. The mammalian ortholog of ATG5 has been shown to translocate to the nucleus under stress conditions and play a role in arresting cell division [35]. Furthermore, deacetylation-regulated nuclear export of a mammalian ATG8 ortholog was suggested to be implicated in autophagy [36]. In addition, plant ATG8s were previously observed in both cytoplasm and nuclei [5,17,37]. Therefore, we decided to investigate if the plant ATG5 complex components might also be localizing to the nuclei. For this we compared localization of ATG5, ATG8, ATG12 (A and B isoforms), and ATG16 fluorescent fusions transiently expressed in Nicotiana benthamiana epidermal leaf cells (Figure 3). To our surprise, ATG5 and both ATG12 isoforms, but not ATG16, could indeed localize to the nuclei under normal conditions, similarly to the previously observed localization of ATG8 (Figure 3A).
To further elucidate whether ATG5 and ATG12 localize to the nucleus in a form of the stable ATG5-ATG12 conjugate, we checked localization of the ATG5K128R defective in conjugation to ATG12 and of the of ATG12A-GFP, in which GFP impedes conjugation with ATG5. We observed that ATG5K128R and ATG12A-GFP still could be found in the nuclei (Figure 3B), suggesting that conjugation is dispensable for nuclear localization of these proteins. These experiments were performed using GFP, YFP and TagRFP fusions of the proteins to ensure that observed localization was not an artefact caused by the fluorescent tag (only representative data for one of the fusion types is shown for each protein of interest). Finally, we used Western Blot analysis to exclude the possibility that nuclear signal observed by confocal microscopy resulted from passive diffusion of the fluorescent tag cleaved off the expressed fusion proteins. Detection of GFP and TagRFP in the protein extracts from leaves transiently expressing fluorescent fusions of ATG5, ATG12 and ATG16 confirmed the presence of intact fluorescent fusions (Figure 3C). In sum, these results demonstrated that at least three components of the autophagy ubiquitin-like conjugation system (ATG5, ATG12 and ATG8) localize to the nuclei. Furthermore, our affinity purification assay indicates that ATG5 might shuttle between nuclei and cytoplasm through interaction with the nuclear pore complex components NUP155 and SEC13b.
The molecular mechanisms underpinning typical autophagy-deficient plant phenotypes such as premature senescence and early onset and cessation of flowering have been only partially explained [9,38,39]. Our pull-down assay indicated interaction between ATG5 and nuclear proteins involved in stress response, flowering and photomorphogenesis (Table S3). Those include HAM1, a catalytic subunit of NuA4 acetyltransferase complex with a role in chromatin remodelling during environmental response of Arabidopsis [40], lectin EULS3 involved in osmotic stress response and ABA signalling [41] and Nuclear Factor YC protein 3 (NF-YC3) playing a role in photomorphogenesis, drought response, flowering and ABA signalling [42,43], COP9 signalosome subunit 1 involved in photomorphogenesis [44] and CULLIN-ASSOCIATED AND NEDDYLATION DISSOCIATED 1 (CAND1), a known mediator of auxin signaling and flowering [45].

3. Discussion

3.1. Evolutionary Context of Autophagy

Compartmentalization of biosynthesis and catabolism in eukaryotic cells lead to development of an intricate endomembrane trafficking system comprising a set of membrane-bound organelles, each having a specialized function in the multistep delivery of cellular materials to their respective destinations. For example, proteins synthesized on the ER can be delivered to the lytic compartment or the plasma membrane (PM) via stepwise transport through the Golgi apparatus and the endosomal system [46]. Similarly, degradation of the cellular content might also rely on directed trafficking towards a specialized compartment [46,47].
Autophagy and the ubiquitin-proteasome system (UPS) are the two major catabolic pathways in eukaryotic cells that counterbalance biosynthesis. UPS, which relies solely on protein-based molecular machinery for substrate recognition, labelling and degradation [48], is likely a more ancient catabolic process compared to the autophagic pathway. Unlike UPS, autophagy [48] is an integral part of the endomembrane trafficking system and requires formation of a membrane compartment (autophagosomes) to sequester cargo and deliver it for degradation to the lytic compartment, and thus would not be feasible in prokaryotic organisms lacking endomembrane compartments. Indeed, a prototype of the eukaryotic UPS-like system was discovered in Archaea [49], while autophagy is known to be typical for eukaryotic organisms [43].
A significant increase in cell structure complexity of eukaryotes created a necessity for a catabolic pathway able to degrade such cargoes as large as complete organelles. Autophagy was possibly established as such endomembrane trafficking-based high-throughput catabolic mechanism in the last eukaryotic common ancestor (LECA) [43], which already possessed a functioning UPS. Subsequently, although UPS and autophagy are often referred to as separate catabolic pathways, they might actually comprise two branches of a single proteolytic network in which autophagy was established after UPS. Co-evolution of both pathways sculpted intricate mechanisms of mutual regulation, implementing common molecular tools such as ubiquitin-like conjugation systems and overlapping degrons [50].
Since ATG proteins have also evolved as an integral part of the endomembrane trafficking system, it is not unlikely they might participate in the routes of the endomembrane trafficking that are not directly linked to autophagy, in its classical depiction. One such example is the cytoplasm-to-vacuole targeting(cvt) pathway in yeast. This is the only known biosynthetic pathway that utilizes autophagic machinery, where core ATG proteins, with a minor change in complex formation deliver enzymes to the vacuole [51]. Furthermore, the endomembrane systems of plant, animal and fungi have significant structural and functional differences, e.g., the lytic compartments are represented either by vacuoles or by lysosomes, the Golgi apparatus is either a static ribbon like structure or mobile stack or even individual cisternae, and the formation and function of endosomes varies significantly [52,53,54,55]. The evolution of autophagy as the integral part of such endomembrane systems likely diverged to develop kingdom-specific features. There are two pending questions. First, how did the evolution of the endomembrane system impact autophagy-unrelated functions of ATG proteins while preserving their role in autophagy and cross-talk with UPS? Second, did ATG proteins carry over some autophagy-unrelated functions from their distant homologs found in prokaryotes [43]?

3.2. Closely Related Components of Endomembrane Trafficking Pathways Show Selectivity towards Different ATG5 Baits

In this study we initiated a discovery proteomics approach to expand our understanding of plant ATG5 in autophagic pathway and beyond it. By comparing proteins pulled-down using WT or autophagy-incompetent K128R mutant of ATG5 as a bait we identified potential interactors of complete and minimal ATG5 complexes, and of individual ATG5.
Remarkably, we identified examples of close homologs presumably involved in the same pathway interacting with ATG5 but showing different preferences towards autophagy-competent or autophagy-incompetent forms of the bait. These results can be further developed into a new tool to elucidate specific functions of these homologs and bring important insights on the roles of ATG5 in and beyond autophagy.
For example, Secretion-Associated-Ras-related GTPase (SAR1) plays a crucial role in the initiation of COPII vesicles formation and thus enables ER to Golgi trafficking [56]. Arabidopsis genome encodes five paralogs of SAR1A-E with suggested tissue specific expression and at least partial functional diversity [57,58,59]. Only SAR1B and SAR1D were shown to play a role in autophagosome maturation [60] and autophagosome biogenesis [61], respectively. Our pull-down assay identified SAR1B as an interactor of the minimal ATG5 complex and SAR1C as an interactor of specifically ATG5K128R. The latter interaction was most likely to be an example of autophagy-independent role of ATG5, as SAR1C was previously shown to not co-localize with autophagic structures under autophagy-inducing conditions [61].
Another example of such specificity towards the two ATG5 baits are sorting nexins (SNXs), subunits of a conserved among eukaryotes retromer complex, which is essential for retrograde trafficking and autophagy [62,63]. Elucidating the interactome of plant SNXs is especially interesting, as, despite being conserved, they seem to be dispensable for the functionality of the plant retromer complex [64] and might have evolved plant-specific functions instead [65]. In this study, we observed SNX1 being pulled-down only with the WT form of ATG5 (complete complex) and SNX2A only with the ATG5K128R (individual ATG5). Additionally, we also identified CLASP protein required for SNX1 membrane association [66] as a potential interactor of the ATG5K128R. Future studies will help to elucidate the mechanism behind these intriguing specific preferences of SNX1 and SNX2A towards either WT or K128R mutant of ATG5, respectively, and potentially bring better understanding of the individual roles of SNXs.
Additionally, in this study, out of three existing Arabidopsis paralogs of SEC24 [57] which is a COPII coat complex component regulating vesicle formation during ER to Golgi trafficking, we identified only SEC24C as an interactor of complete ATG5 complex. Interestingly, the closest mammalian ortholog of SEC24C was shown to play a role in selective autophagy of ER (ER-phagy) [67]. Future studies involving ER-phagy-inducing conditions will help to assess the potential role of SEC24C in plant ER-phagy and the role of ATG5 in it.

3.3. Cross-Talk of ATG5 and Clathrin-Mediated Trafficking Is Conserved

Lastly, we uncovered a set of ATG5 interactors that are associated with ER-PMCS (endoplasmic reticulum- plasma membrane contact sites) and clathrin-mediated vesicular trafficking (TPLATE COMPLEX MUNISCIN-LIKE (TML), SH3 DOMAIN-CONTAINING PROTEIN2 (SH3P2), and SYNAPTOTAGMIN 1 (SYT). All three proteins are known to partake in autophagy [68,69,70]. However, our study provides novel evidence of their interaction with the ATG5 complexes. Interestingly, a recent publication on the ATG5 interactome in mice [12] also presented several interactors of ATG5 related to clathrin-mediated trafficking, linking this pathway to animal autophagy.

4. Methods

4.1. Plasmids Construction

For generation of 2x35S:ATG5-TAP constructs, the ATG5 gene AT5G17290.1 without UTRs was amplified using PE7 and PE8 primers (Table S4) and cloned into pDONR/Zeo vector using Gateway cloning system (Invitrogen/Thermo Fisher Scientific, Waltham, MA USA). The point mutation was introduced into this entry clone using primers PE10 and PE11 (Table S4) and QuikChange II Site-Directed Mutagenesis Kit (200523, Agilent Technologies, Inc., Santa Clara, CA, USA), following the kit’s instructions. The obtained entry clones with WT and K128R gene versions of ATG5 were then recombined into pCTAP vector (pYL436) under the control of 2x35S promoter [16] using Gateway cloning system (Invitrogen).
For generation of ATG5pr::ATG5-TAP constructs, the ATG5 gene AT5G17290.1 together with the promoter region was amplified using primers PE3 and PE4 (Table S4). The TAP tag was amplified from the pCTAP vector [16] using primers PE5 and PE9 (Table S4). The ATG5 gene and TAP tag were then fused using overlay PCR with primers PE3 and PE9 (Table S4). The obtained amplicon was cloned into pDONR/Zeo vector using Gateway cloning system (Invitrogen). The point mutation was introduced into this entry clone using primers PE10 and PE11 (Table S4) and QuikChange II Site-Directed Mutagenesis Kit (Agilent, 200523), according to the standard protocol. The obtained entry clones with WT and K128R gene versions of ATG5 genes including promoter region were then recombined into pGWB401 vector (Adgene, Plasmid #74795) using Gateway cloning system (Invitrogen).
Constructs for expressing fluorescent fusions of ATG5, ATG8, ATG16, ATG12A and ATG12B were produced using Gateway cloning system (Invitrogen). The corresponding genes or Coding DNA sequences were lifted from genomic DNA or total cDNA of Arabidopsis using primers provided in the Table S4 (PE8, SH139-145, AM 475, 476). Amplicons were recombined into pDONR/Zeo vecor to produce entry clones (see Table S4: SH 210-213, SH 231, 232 and AM 655). The entry clones were later recombined into pGWB-series destination vectors carrying 2x35S promoter and GFP, YFP or TagRFP tags (Adgene) (see Table S4, destination clones for expressing fluorescent fusions in plants).
The destination clones were used to transform Agrobacterium tumefaciens strain GV3101.

4.2. Protoplast and Plant Transformation

4.2.1. Transient Expression of Arabidopsis Protoplasts

Protoplasts were isolated from leaves of four-weeks-old plants expressing GFP-ATG8 in atg5-1 using the “Tape-Arabidopsis Sandwich” method described in [71]. The isolated protoplasts were transformed using 10–20 µg of each plasmid (Table S4) using the method described in [72]. The transfected protoplasts were incubated in 24-well glass bottom plates (VWR CORN4441) for 16 h in light. The protoplasts were further treated with 5 µM AZD-8055 (364424, Santa-Cruz Biotech, Dallas TX, USA) and 0.5 µM Concanamycin A (202111A, Santa-Cruz Biotech) for 24 h, where applicable. The transformed protoplasts were imaged using CLSM800 (Carl Zeiss AG, Oberkochen BW, Germany), objective C-Apochromat 40×/1.2 W, excitation light 488 nm and 561nm and emission ranges of (515–560 nm) and (570–650 nm) for GFP and TagRFP, respectively. Images were analyzed using ZEN blue software (Carl Zeiss).

4.2.2. Transient Expression in Nicotiana Benthamiana

N. benthamiana plants were grown in 8 cm3 pots filled with soil S-Jord (Hasselfors) under controlled growth conditions of 16 h light 8 h dark cycles, 70% relative humidity, light intensity of 150 μE m–2 s–1, and day and night temperature of 22 °C and 20 °C, respectively. Transformed A. tumefaciens strain GV3101 (Rifampicin, 100 μg/mL) carrying the constructs of interest were grown in 5 mL of Luria-Broth high salt medium (L1704, Duchefa Biochemi, Haarlem, The Netherlands), supplemented with appropriate antibiotics (see Table S4, destination clones for expressing fluorescent fusions in plants). Liquid bacterial cultures were shaken at 200 rpm, 28 °C overnight, and then sedimented at 4000 G for five minutes. The resulting pellets were resuspended in MQ water with 150 μM Acetosyringone to the final OD600 = 0.15 and infiltrated in the abaxial side of leaves of five-week-old N. benthamiana plants. The leaves were imaged on the third day post-infiltration using CLSM800 (Carl Zeiss), objective C-Apochromat 40×/1.2 W, excitation light 488 nm and 561 nm and emission ranges of (515–560 nm) and (570–650 nm) for GFP and RFP, respectively. Images were analyzed using ZEN black software (Carl Zeiss).

4.2.3. Arabidopsis Thaliana Growth and Transformation

Arabidopsis plants were grown in 8 cm3 pots filled with soil S-Jord (Hasselfors), under long day conditions: 150 µE m−2 s−1 light for 16 h, 8 h dark, 22 °C, 70% humidity.
Arabidopsis Col-0 wild-type plants and the previously described atg5-1 mutant [73] were transformed using the standard “floral dip”-method [74]. A. tumefaciens strain GV3101 carrying the pYL436 and pGWB401 constructs (Table S4) was used for transformation. Transgenic plants were selected on Murashige and Skoog (MS) medium containing 40 µg mL−1 gentamycin and 50 µg mL−1 kanamycin, respectively. The plants were genotyped to confirm the presence of the transgenes using PE7 and PE8 primers (Table S4). Expression of the transgenes was confirmed using Western blot analyses (as described in the Section 4.3).

4.3. Immunoblotting

Plant material was powdered in liquid nitrogen, mixed with 2 vol. of hot 2× Laemmli buffer, and boiled for 10 min. Debris was pelleted for 5 min at 17,000 G. Proteins were separated on Mini-PROTEAN TGX stain-free Bio-Rad gels, 7.5% (Bio-Rad, Hercules, CA, USA) and transferred onto PVDF membranes. Membranes were blotted anti-actin 1:2000 (AS13 2640, Agrisera AB Vännäs, SWEDEN.), anti-myc 1:1000 (11667203001, Roche, Basel, Switzerland). Reactions were developed using Amersham ECL Prime kit (RPN2232, Cytiva, Marlborough, MA, USA) and detected using BioRad ChemiDoc.

4.4. Tandem Affinity Purification

Plants for affinity purification assay were grown in 8 cm3 pots filled with soil S-Jord (Hasselfors) under long day conditions: 150 µE m−2 s−1 light for 16 h, 8 h dark, 22 °C, 70% humidity.
Plant rosette leaves were harvested at bolting stage and flash frozen in liquid nitrogen. Sample weight was >15 g of frozen tissue. Subsequent procedure was performed as previously described [18] using the following consumables: Protease inhibitor cocktail (P9599-5ML, Sigma-Aldrich, Saint Louis, MO, USA), PreScission protease (GE27-0843-01, Sigma-Aldrich), Ig-G sepharose beads (GE17-0969-01, Sigma-Aldrich),), Ni-sepharose beads (GE17-5318-01, Sigma-Aldrich).
Mass spectrometry data was acquired using a data-dependent acquisition procedure with a cyclic series of a full scan from 350–1500 with resolution of 120,000 control (AGC) target 1E6, maximum injection time 100 ms. The top S (3 s) and dynamic exclusion of 30 s were used for selection of Parent ions for MSMS in the HCD cell with, relative collision energy 30% and scanned in the orbitrap with resolution of 30,000.

4.5. Mass Spectrometry

4.5.1. Liquid Chromatography and Mass Spectrometry Assay

LC-MS/MS was performed at the Mass Spectrometry Facility of Rutgers Center for Advanced Biotechnology and Medicine, US. Samples were loaded on to a fused silica trap column Acclaim PepMap100, 75 µm × 2 cm (Thermo Fisher Scientific, Waltham, MA USA).). After washing for 5 min at 5 μL/min with 0.1% TFA, the trap column was brought in-line with an analytical column (NanoeaseMZ peptideBEH C18, 130A, 1.7 µm, 75 mm × 250 mm, Waters, Milford, MA, USA) for LC-MS/MS. Peptides were fractionated at 300 nL/min using a segmented linear gradient 4–15% B in 30 min (where A: 0.2% formic acid, and B: 0.16% formic acid, 80% acetonitrile), 15–25% B in 40 min, 25–50% B in 44 min, and 50–90% B in 11 min. The column was re-equilibrated with 4% B for 5 min prior to the next run.

4.5.2. LC-MS/MS Data Analysis

Raw data files from the LC-MS/MS analysis were processed using MaxQuant (ver. 1.6.17.0) [19]. Default settings were used except enabling the protein label-free quantification (LFQ) and matching-between-runs options on default settings. Proteomics data was searched against the 26,755 representative Araport11 proteins [75] supplemented with the contaminant proteins list included within MaxQuant.
For the search including PTMs, phosphorylation (STY), acetylation of the side chain of Lys residues and protein N-terminal acetylation were specified in MaxQuant as variable modifications, while cysteine carbamidomethylation was set as fixed. Additionally, multiple events of methionine oxidation were detected, but not taken into consideration for this study, as it was not possible to discern whether they occurred in vivo or during samples processing [76].

4.5.3. Statistical Analysis of Interactors—MSqRobSum

The obtained MaxQuant peptide and protein result tables were used for statistical analysis in R (ver. 4.02) by MSqRobSum (ver. 0.0.0.9000) [77]. Proteins only identified by site and contaminant/decoy proteins were filtered. After default preprocessing, protein intensity summaries were estimated by robust regression, fitting log2 intensities in function of condition (WT ATG5, PM ATG5, and GFP control). Afterwards, user-defined contrasts of interest were tested to identify wild-type ATG5 (WT ATG5/GFP control) and mutated ATG (PM ATG5/GFP control) interactors. Interactors were filtered as having a fold change ≥ 1.5 and p value ≤ 0.05. In addition, proteins not quantified in GFP control conditions but identified in at least 2 out of 3 replicates of a ATG5 pulldown were manually curated as putative ATG5 interactors.

4.5.4. STRING Analysis

STRING analysis was made using https://string-db.org/ (accessed on 22 November 2022) website [78] with following parameters: full STRING network, active interaction sources: all possible; minimum required interaction score: medium confidence (0.400). Such settings were selected to include also possible connections that are not yet properly verified, which we found most suitable in the context of our search for unconventional roles of the ATG proteins. Suggested by STRING analysis connections were manually verified in the existing literature. When clustering was used it was with MCL clustering with an inflation parameter of 3.

5. Conclusions

This study provides novel evidence of PTMs occurring on core ATG proteins, i.e., ATG5 and ATG12 under non-inducing conditions. This discovery opens up exciting opportunities for further investigation and deeper understanding of the regulatory mechanisms governing ATG proteins in plants. Moreover, the comparative affinity purification assay revealed a compelling list of interactors for complete and minimal ATG5 complexes, indicating the possibility of crosstalk between autophagy-related proteins and multiple components of the endomembrane trafficking system and UPS. These findings strongly suggest the existence of shared and coordinated mechanisms that integrate these crucial cellular processes. Significantly, the identification of unique interactors for both ATG5 and ATG5K128R suggests that ATG5 may possess roles and functions that extend beyond its classical involvement in autophagy. This finding underscores the complexity and versatility of ATG proteins and implies their potential contributions to diverse cellular pathways. To deepen our functional understanding of the complete and minimal ATG5 complexes, future studies should focus on verifying the identified interactors and conducting functional analyses of the candidate protein complexes. Additionally, we present novel evidence of nuclear localization for individual ATG5 and ATG12, that does not depend on their conjugation. Furthermore, we identified a set of nuclear-localized proteins, potential interactors of the nuclear fraction of ATG5 and ATG12. This discovery significantly expands the range of suggested functions associated with the components of the ATG5 complex. We hope that this study will pave the way for further explorations and place autophagic pathway into a broader context.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms241512300/s1, Table S1: MSqRob_without PTM; Table S2: MSqRob with PTM; Table S3: List of identified hits; Table S4: primers, plasmids and transgenic lines used in this study.

Author Contributions

Conceptualization, E.A.M. and P.V.B.; methodology, E.A.M., P.H.E. and P.W.; software, P.W.; validation, E.A.M., S.H. and P.H.E.; formal analysis, P.W., E.A.M., S.H. and P.H.E.; investigation S.H., P.H.E. and E.A.M.; resources, E.G.-B.; data curation, P.W., S.H. and E.A.M.; writing—original draft preparation, S.H., P.H.E. and E.A.M.; writing—review and editing, I.S., P.V.B. and P.W.; visualization, S.H., P.H.E., E.A.M. and I.S.; supervision, E.A.M. and P.V.B.; project administration, E.A.M. and P.V.B.; funding acquisition, E.A.M., P.V.B., P.W. and E.G.-B.; All authors have read and agreed to the published version of the manuscript.

Funding

This project was supported by the grants from The Swedish Research Council Formas (2021-01812 to EAM; 2017-00541 to PVB), Carl Tryggers Foundation (CTS 20:287 to EAM), MSCA IF (799433, to EAM), The Swedish Foundation for Strategic Research (“Oil crops for the future”, to PVB), and “Crops for the Future” Research Programme at The Swedish University of Agricultural Science. Contribution of EGB to this work was supported by Ministerio de Ciencia e Innovacion from Spanish Government, Grant PID2020–119737GA–I00 (MCIN/AEI/10.13039/501100011033) and Junta de Andalucia (ProyExcel_00587). Contribution of PW was supported by Research Foundation-Flanders (FWO) Postdoctoral fellowship 12T1722N.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw LC-MS/MS data is available on the Zenodo: 10.5281/zenodo.8199375.

Acknowledgments

We would like to thank George Cheng (Dept of MolSci, SLU, Sweden) for his help with proofreading the manuscript and Haiyan Zheng (MS facility of Rutgers CAMB, US) for help with the LC-MS/MS assay.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Parzych, K.R.; Klionsky, D.J. Vacuolar hydrolysis and efflux: Current knowledge and unanswered questions. Autophagy 2019, 15, 212–227. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Su, T.; Li, X.; Yang, M.; Shao, Q.; Zhao, Y.; Ma, C.; Wang, P. Autophagy: An Intracellular Degradation Pathway Regulating Plant Survival and Stress Response. Front. Plant Sci. 2020, 11, 164. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Ketelaar, T.; Voss, C.; Dimmock, S.A.; Thumm, M.; Hussey, P.J. Arabidopsis homologues of the autophagy protein Atg8 are a novel family of microtubule binding proteins. FEBS Lett. 2004, 567, 302–306. [Google Scholar] [CrossRef] [Green Version]
  4. Nakatogawa, H.; Ichimura, Y.; Ohsumi, Y. Atg8, a Ubiquitin-like Protein Required for Autophagosome Formation, Mediates Membrane Tethering and Hemifusion. Cell 2007, 130, 165–178. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Minina, E.A.; Moschou, P.N.; Vetukuri, R.R.; Sanchez-Vera, V.; Cardoso, C.; Liu, Q.; Elander, P.H.; Dalman, K.; Beganovic, M.; Yilmaz, J.L.; et al. Transcriptional stimulation of rate-limiting components of the autophagic pathway improves plant fitness. J. Exp. Bot. 2018, 69, 1415–1432. [Google Scholar] [CrossRef]
  6. Bu, F.; Yang, M.; Guo, X.; Huang, W.; Chen, L. Multiple Functions of ATG8 Family Proteins in Plant Autophagy. Front. Cell Dev. Biol. 2020, 8, 466. [Google Scholar] [CrossRef]
  7. Jia, M.; Liu, X.; Xue, H.; Wu, Y.; Shi, L.; Wang, R.; Chen, Y.; Xu, N.; Zhao, J.; Shao, J.; et al. Noncanonical ATG8-ABS3 interaction controls senescence in plants. Nat. Plants 2019, 5, 212–224. [Google Scholar] [CrossRef]
  8. Codogno, P.; Meijer, A.J. Atg5: More than an autophagy factor. Nat. Cell Biol. 2006, 8, 1045–1047. [Google Scholar] [CrossRef]
  9. Yoshimoto, K.; Jikumaru, Y.; Kamiya, Y.; Kusano, M.; Consonni, C.; Panstruga, R.; Ohsumi, Y.; Shirasu, K. Autophagy Negatively Regulates Cell Death by Controlling NPR1-Dependent Salicylic Acid Signaling during Senescence and the Innate Immune Response in Arabidopsis. Plant Cell Online 2009, 21, 2914–2927. [Google Scholar] [CrossRef] [Green Version]
  10. Barros, J.A.; Cavalcanti, J.H.F.; Medeiros, D.B.; Nunes-Nesi, A.; Avin-Wittenberg, T.; Fernie, A.R.; Araújo, W.L. Autophagy deficiency compromises alternative pathways of respiration following energy deprivation in Arabidopsis thaliana. Plant Physiol. 2017, 175, 62–76. [Google Scholar] [CrossRef] [Green Version]
  11. Matsushita, M.; Suzuki, N.N.; Obara, K.; Fujioka, Y.; Ohsumi, Y.; Inagaki, F. Structure of Atg5.Atg16, a complex essential for autophagy. J. Biol. Chem. 2007, 282, 6763–6772. [Google Scholar] [CrossRef] [Green Version]
  12. Baines, K.; Yoshioka, K.; Takuwa, Y.; Lane, J.D. The ATG5 interactome links clathrin-mediated vesicular trafficking with the autophagosome assembly machinery. Autophagy Rep. 2022, 1, 88–118. [Google Scholar] [CrossRef] [PubMed]
  13. Jumper, J.; Evans, R.; Pritzel, A.; Green, T.; Figurnov, M.; Ronneberger, O.; Tunyasuvunakool, K.; Bates, R.; Žídek, A.; Potapenko, A.; et al. Highly accurate protein structure prediction with AlphaFold. Nature 2021, 596, 583–589. [Google Scholar] [CrossRef] [PubMed]
  14. Varadi, M.; Anyango, S.; Deshpande, M.; Nair, S.; Natassia, C.; Yordanova, G.; Yuan, D.; Stroe, O.; Wood, G.; Laydon, A.; et al. AlphaFold Protein Structure Database: Massively expanding the structural coverage of protein-sequence space with high-accuracy models. Nucleic Acids Res. 2022, 50, D439–D444. [Google Scholar] [CrossRef] [PubMed]
  15. Romanov, J.; Walczak, M.; Ibiricu, I.; Schüchner, S.; Ogris, E.; Kraft, C.; Martens, S. Mechanism and functions of membrane binding by the Atg5-Atg12/Atg16 complex during autophagosome formation. EMBO J. 2012, 31, 4304–4317. [Google Scholar] [CrossRef]
  16. Rubio, V.; Shen, Y.; Saijo, Y.; Liu, Y.; Gusmaroli, G.; Dinesh-Kumar, S.P.; Deng, X.W. An alternative tandem affinity purification strategy applied to Arabidopsis protein complex isolation. Plant J. 2005, 41, 767–778. [Google Scholar] [CrossRef]
  17. Dauphinee, A.N.; Cardoso, C.; Dalman, K.; Ohlsson, J.A.; Fick, S.B.; Robert, S.; Hicks, G.R.; Bozhkov, P.V.; Minina, E.A. Chemical screening pipeline for identification of specific plant autophagy modulators. Plant Physiol. 2019, 181, 855–866. [Google Scholar] [CrossRef] [Green Version]
  18. Gutierrez-Beltran, E.; Elander, P.H.; Dalman, K.; Dayhoff, G.W.; Moschou, P.N.; Uversky, V.N.; Crespo, J.L.; Bozhkov, P.V. Tudor staphylococcal nuclease is a docking platform for stress granule components and is essential for SnRK1 activation in Arabidopsis. EMBO J. 2021, 40, e105043. [Google Scholar] [CrossRef]
  19. Cox, J.; Mann, M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat. Biotechnol. 2008, 26, 1367–1372. [Google Scholar] [CrossRef]
  20. Schreiber, A.; Collins, B.C.; Davis, C.; Enchev, R.I.; Sedra, A.; D’antuono, R.; Aebersold, R.; Peter, M. Multilayered regulation of autophagy by the Atg1 kinase orchestrates spatial and temporal control of autophagosome formation. Mol. Cell 2021, 81, 5066–5081.e10. [Google Scholar] [CrossRef]
  21. Lee, I.H.; Finkel, T. Regulation of autophagy by the p300 acetyltransferase. J. Biol. Chem. 2009, 284, 6322–6328. [Google Scholar] [CrossRef] [Green Version]
  22. Lee, I.H.; Cao, L.; Mostoslavsky, R.; Lombard, D.B.; Liu, J.; Bruns, N.E.; Tsokos, M.; Alt, F.W.; Finkel, T. A role for the NAD-dependent deacetylase Sirt1 in the regulation of autophagy. Proc. Natl. Acad. Sci. USA 2008, 105, 3374–3379. [Google Scholar] [CrossRef]
  23. Yeasmin, A.M.S.T.; Waliullah, T.M.; Kondo, A.; Kaneko, A.; Koike, N.; Ushimaru, T. Orchestrated action of PP2A antagonizes atg13 phosphorylation and promotes autophagy after the inactivation of TORC1. PLoS ONE 2016, 11, e0166636. [Google Scholar] [CrossRef] [Green Version]
  24. Bánréti, Á.; Lukácsovich, T.; Csikós, G.; Erdélyi, M.; Sass, M. PP2A regulates autophagy in two alternative ways in Drosophila. Autophagy 2012, 8, 623–636. [Google Scholar] [CrossRef] [Green Version]
  25. Khan, Z.H.; Agarwal, S.; Rai, A.; Memaya, M.B.; Mehrotra, S.; Mehrotra, R. Co-expression network analysis of protein phosphatase 2A (PP2A) genes with stress-responsive genes in Arabidopsis thaliana reveals 13 key regulators. Sci. Rep. 2020, 10, 21480. [Google Scholar] [CrossRef] [PubMed]
  26. Huang, L.; Yu, L.-J.; Zhang, X.; Fan, B.; Wang, F.-Z.; Dai, Y.-S.; Qi, H.; Zhou, Y.; Xie, L.-J.; Xiao, S. Autophagy regulates glucose-mediated root meristem activity by modulating ROS production in Arabidopsis. Autophagy 2019, 15, 407–422. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Li, B.; Zeng, Y.; Jiang, L. COPII vesicles in plant autophagy and endomembrane trafficking. FEBS Lett. 2022, 596, 2314–2323. [Google Scholar] [CrossRef]
  28. Itakura, E.; Kishi-Itakura, C.; Mizushima, N. The hairpin-type tail-anchored SNARE syntaxin 17 targets to autophagosomes for fusion with endosomes/lysosomes. Cell 2012, 151, 1256–1269. [Google Scholar] [CrossRef] [Green Version]
  29. Aniento, F.; Hernández, V.S.d.M.; Dagdas, Y.; Rojas-Pierce, M.; Russinova, E. Molecular mechanisms of endomembrane trafficking in plants. Plant Cell 2022, 34, 146–173. [Google Scholar] [CrossRef] [PubMed]
  30. Minina, E.A.; Moschou, P.N.; Bozhkov, P.V. Limited and digestive proteolysis: Crosstalk between evolutionary conserved pathways. New Phytol. 2017, 215, 958–964. [Google Scholar] [CrossRef] [Green Version]
  31. Li, Y.; Li, S.; Wu, H. Ubiquitination-proteasome system (UPS) and autophagy two main protein degradation machineries in response to cell stress. Cells 2022, 11, 851. [Google Scholar] [CrossRef]
  32. Kocaturk, N.M.; Peker, N.; Eberhart, K.; Akkoc, Y.; Deveci, G.; Dengjel, J.; Gozuacik, D. Novel protein complexes containing autophagy and UPS components regulate proteasome-dependent PARK2 recruitment onto mitochondria and PARK2-PARK6 activity during mitophagy. Cell Death Dis. 2022, 13, 947. [Google Scholar] [CrossRef] [PubMed]
  33. Li, J.; Yuan, J.; Li, Y.; Sun, H.; Ma, T.; Huai, J.; Yang, W.; Zhang, W.; Lin, R. The CDC48 complex mediates ubiquitin-dependent degradation of intra-chloroplast proteins in plants. Cell Rep. 2022, 39, 110664. [Google Scholar] [CrossRef] [PubMed]
  34. Tamura, K.; Fukao, Y.; Iwamoto, M.; Haraguchi, T.; Hara-Nishimura, I. Identification and characterization of nuclear pore complex components in Arabidopsis thaliana. Plant Cell 2010, 22, 4084–4097. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Ye, X.; Zhou, X.-J.; Zhang, H. Exploring the role of autophagy-related gene 5 (ATG5) yields important insights into autophagy in autoimmune/autoinflammatory diseases. Front. Immunol. 2018, 9, 2334. [Google Scholar] [CrossRef]
  36. Huang, R.; Xu, Y.; Wan, W.; Shou, X.; Qian, J.; You, Z.; Liu, B.; Chang, C.; Zhou, T.; Lippincott-Schwartz, J.; et al. Deacetylation of nuclear LC3 drives autophagy initiation under starvation. Mol. Cell 2015, 57, 456–467. [Google Scholar] [CrossRef] [Green Version]
  37. Li, F.; Zhang, M.; Zhang, C.; Zhou, X. Nuclear autophagy degrades a geminivirus nuclear protein to restrict viral infection in solanaceous plants. New Phytol. 2020, 225, 1746–1761. [Google Scholar] [CrossRef]
  38. Minina, E.A.; Sanchez-Vera, V.; Moschou, P.N.; Suarez, M.F.; Sundberg, E.; Weih, M.; Bozhkov, P.V. Autophagy mediates caloric restriction-induced lifespan extension in Arabidopsis. Aging Cell 2013, 12, 327–329. [Google Scholar] [CrossRef]
  39. Su, W.; Bao, Y.; Yu, X.; Xia, X.; Liu, C.; Yin, W. Autophagy and its regulators in response to stress in plants. Int. J. Mol. Sci. 2020, 21, 8889. [Google Scholar] [CrossRef] [PubMed]
  40. Bieluszewski, T.; Sura, W.; Dziegielewski, W.; Bieluszewska, A.; Lachance, C.; Kabza, M.; Szymanska-Lejman, M.; Abram, M.; Wlodzimierz, P.; De Winne, N.; et al. NuA4 and H2A.Z control environmental responses and autotrophic growth in Arabidopsis. Nat. Commun. 2022, 13, 277. [Google Scholar] [CrossRef] [PubMed]
  41. Dubiel, M.; Beeckman, T.; Smagghe, G.; Van Damme, E.J.M. Arabidopsis lectin EULS3 is involved in ABA signaling in roots. Front. Plant Sci. 2020, 11, 437. [Google Scholar] [CrossRef]
  42. Hwang, K.; Susila, H.; Nasim, Z.; Jung, J.-Y.; Ahn, J.H. Arabidopsis ABF3 and ABF4 Transcription Factors Act with the NF-YC Complex to Regulate SOC1 Expression and Mediate Drought-Accelerated Flowering. Mol. Plant 2019, 12, 489–505. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Zhang, S.; Hama, Y.; Mizushima, N. The evolution of autophagy proteins—Diversification in eukaryotes and potential ancestors in prokaryotes. J. Cell Sci. 2021, 134, jcs233742. [Google Scholar] [CrossRef] [PubMed]
  44. Nakasone, A.; Fujiwara, M.; Fukao, Y.; Biswas, K.K.; Rahman, A.; Kawai-Yamada, M.; Narumi, I.; Uchimiya, H.; Oono, Y. SMALL ACIDIC PROTEIN1 acts with RUB modification components, the COP9 signalosome, and AXR1 to regulate growth and development of Arabidopsis. Plant Physiol. 2012, 160, 93–105. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Feng, S.H.; Shen, Y.P.; Sullivan, J.A.; Rubio, V.; Xiong, Y.; Sun, T.-P.; Deng, X.W. Arabidopsis CAND1, an unmodified CUL1-interacting protein, is involved in multiple developmental pathways controlled by ubiquitin/proteasome-mediated protein degradation. Plant Cell 2004, 16, 1870–1882. [Google Scholar] [CrossRef] [Green Version]
  46. Huotari, J.; Helenius, A. Endosome maturation. EMBO J. 2011, 30, 3481–3500. [Google Scholar] [CrossRef]
  47. Di Sansebastiano, G.P.; Barozzi, F.; Piro, G.; Denecke, J.; de Marcos Lousa, C. Trafficking routes to the plant vacuole: Connecting alternative and classical pathways. J. Exp. Bot. 2018, 69, 79–90. [Google Scholar] [CrossRef] [Green Version]
  48. Hershko, A.; Heller, H.; Elias, S.; Ciechanover, A. Components of ubiquitin-protein ligase system. Resolution, affinity purification, and role in protein breakdown. J. Biol. Chem. 1983, 258, 8206–8214. [Google Scholar] [CrossRef]
  49. Nunoura, T.; Takaki, Y.; Kakuta, J.; Nishi, S.; Sugahara, J.; Kazama, H.; Chee, G.-J.; Hattori, M.; Kanai, A.; Atomi, H.; et al. Insights into the evolution of Archaea and eukaryotic protein modifier systems revealed by the genome of a novel archaeal group. Nucleic Acids Res. 2011, 39, 3204–3223. [Google Scholar] [CrossRef] [Green Version]
  50. Ji, C.H.; Kwon, Y.T. Crosstalk and interplay between the ubiquitin-proteasome system and autophagy. Mol. Cells 2017, 40, 441. [Google Scholar]
  51. Yamasaki, A.; Noda, N.N. Structural Biology of the Cvt Pathway. J. Mol. Biol. 2017, 429, 531–542. [Google Scholar] [CrossRef] [PubMed]
  52. Nakano, A. The Golgi apparatus and its next-door neighbors. Front. Cell Dev. Biol. 2022, 10, 884360. [Google Scholar] [CrossRef]
  53. Klionsky, D.J.; Eskelinen, E.-L. The vacuole vs. the lysosome. Autophagy 2014, 10, 185–187. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Robinson, D.G. The Golgi Apparatus and the Plant Secretory Pathway; Blackwell: Oxford, UK, 2003; Volume 9. [Google Scholar]
  55. Contento, A.L.; Bassham, D.C. Structure and function of endosomes in plant cells. J. Cell Sci. 2012, 125, 3511–3518. [Google Scholar] [CrossRef] [Green Version]
  56. Yorimitsu, T.; Sato, K.; Takeuchi, M. Molecular mechanisms of Sar/Arf GTPases in vesicular trafficking in yeast and plants. Front. Plant Sci. 2014, 5, 411. [Google Scholar] [CrossRef] [Green Version]
  57. Chung, K.P.; Zeng, Y.; Jiang, L. COPII paralogs in plants: Functional redundancy or diversity? Trends Plant Sci. 2016, 21, 758–769. [Google Scholar] [CrossRef]
  58. Hanton, S.L.; Chatre, L.; Matheson, L.A.; Rossi, M.; Held, M.A.; Brandizzi, F. Plant Sar1 isoforms with near-identical protein sequences exhibit different localisations and effects on secretion. Plant Mol. Biol. 2008, 67, 283–294. [Google Scholar] [CrossRef]
  59. Liang, X.; Li, S.-W.; Gong, L.-M.; Li, S.; Zhang, Y. COPII components Sar1b and Sar1c play distinct yet interchangeable roles in pollen development. Plant Physiol. 2020, 183, 974–985. [Google Scholar] [CrossRef] [Green Version]
  60. Kim, J.H.; Lee, H.N.; Huang, X.; Jung, H.; Otegui, M.S.; Li, F.; Chung, T. FYVE2, a phosphatidylinositol 3-phosphate effector, interacts with the COPII machinery to control autophagosome formation in Arabidopsis. Plant Cell 2022, 34, 351–373. [Google Scholar] [CrossRef] [PubMed]
  61. Li, B.; Zeng, Y.; Cao, W.; Zhang, W.; Cheng, L.; Yin, H.; Wu, Q.; Wang, X.; Huang, Y.; Lau, W.C.Y.; et al. A distinct giant coat protein complex II vesicle population in Arabidopsis thaliana. Nat. Plants 2021, 7, 1335–1346. [Google Scholar] [CrossRef] [PubMed]
  62. Zavodszky, E.; Seaman, M.N.; Moreau, K.; Jimenez-Sanchez, M.; Breusegem, S.Y.; Harbour, M.E.; Rubinsztein, D.C. Mutation in VPS35 associated with Parkinson’s disease impairs WASH complex association and inhibits autophagy. Nat. Commun. 2014, 5, 3828. [Google Scholar] [CrossRef] [Green Version]
  63. Rodriguez-Furlan, C.; Domozych, D.; Qian, W.; Enquist, P.-A.; Li, X.; Zhang, C.; Schenk, R.; Winbigler, H.S.; Jackson, W.; Raikhel, N.V.; et al. Interaction between VPS35 and RABG3f is necessary as a checkpoint to control fusion of late compartments with the vacuole. Proc. Natl. Acad. Sci. USA 2019, 116, 21291–21301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Pourcher, M.; Santambrogio, M.; Thazar, N.; Thierry, A.-M.; Fobis-Loisy, I.; Miège, C.; Jaillais, Y.; Gaude, T. Analyses of sorting nexins reveal distinct retromer-subcomplex functions in development and protein sorting in Arabidopsis thaliana. Plant Cell 2010, 22, 3980–3991. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Heucken, N.; Ivanov, R. The retromer, sorting nexins and the plant endomembrane protein trafficking. J. Cell Sci. 2018, 131, jcs203695. [Google Scholar] [CrossRef] [Green Version]
  66. Ambrose, C.; Ruan, Y.; Gardiner, J.; Tamblyn, L.M.; Catching, A.; Kirik, V.; Marc, J.; Overall, R.; Wasteneys, G.O. CLASP interacts with sorting nexin 1 to link microtubules and auxin transport via PIN2 recycling in Arabidopsis thaliana. Dev. Cell 2013, 24, 649–659. [Google Scholar] [CrossRef] [Green Version]
  67. Cui, Y.; Parashar, S.; Zahoor, M.; Needham, P.G.; Mari, M.; Zhu, M.; Chen, S.; Ho, H.-C.; Reggiori, F.; Farhan, H.; et al. A COPII subunit acts with an autophagy receptor to target endoplasmic reticulum for degradation. Science (1979) 2019, 365, 53–60. [Google Scholar] [CrossRef]
  68. Wang, P.; Pleskot, R.; Zang, J.; Winkler, J.; Wang, J.; Yperman, K.; Zhang, T.; Wang, K.; Gong, J.; Guan, Y.; et al. Plant AtEH/Pan1 proteins drive autophagosome formation at ER-PM contact sites with actin and endocytic machinery. Nat. Commun. 2019, 10, 5132. [Google Scholar] [CrossRef] [Green Version]
  69. Zhuang, X.; Wang, H.; Lam, S.K.; Gao, C.; Wang, X.; Cai, Y.; Jiang, L. A BAR-domain protein SH3P2, which binds to phosphatidylinositol 3-phosphate and ATG8, regulates autophagosome formation in Arabidopsis. Plant Cell 2013, 25, 4596–4615. [Google Scholar] [CrossRef] [Green Version]
  70. Nascimbeni, A.C.; Giordano, F.; Dupont, N.; Grasso, D.; I Vaccaro, M.; Codogno, P.; Morel, E. ER–plasma membrane contact sites contribute to autophagosome biogenesis by regulation of local PI 3P synthesis. EMBO J. 2017, 36, 2018–2033. [Google Scholar] [CrossRef] [PubMed]
  71. Wu, F.-H.; Shen, S.-C.; Lee, L.-Y.; Lee, S.-H.; Chan, M.-T.; Lin, C.-S. Tape-Arabidopsis Sandwich-a simpler Arabidopsis protoplast isolation method. Plant Methods 2009, 5, 16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Yoo, S.-D.; Cho, Y.-H.; Sheen, J. Arabidopsis mesophyll protoplasts: A versatile cell system for transient gene expression analysis. Nat. Protoc. 2007, 2, 1565–1572. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Thompson, A.R.; Doelling, J.H.; Suttangkakul, A.; Vierstra, R.D. Autophagic nutrient recycling in Arabidopsis directed by the ATG8 and ATG12 conjugation pathways. Plant Physiol. 2005, 138, 2097–2110. [Google Scholar] [CrossRef] [Green Version]
  74. Clough, S.J.; Bent, A.F. Floral dip: A simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana. Plant J. 1998, 16, 735–743. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Cheng, C.; Krishnakumar, V.; Chan, A.P.; Thibaud-Nissen, F.; Schobel, S.; Town, C.D. Araport11: A complete reannotation of the Arabidopsis thaliana reference genome. Plant J. 2017, 89, 789–804. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Chen, M.; Cook, K.D. Oxidation artifacts in the electrospray mass spectrometry of Aβ peptide. Anal. Chem. 2007, 79, 2031–2036. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Sticker, A.; Goeminne, L.; Martens, L.; Clement, L. Robust summarization and inference in proteome-wide label-free quantification. Mol. Cell. Proteom. 2020, 19, 1209–1219. [Google Scholar] [CrossRef] [Green Version]
  78. Szklarczyk, D.; Gable, A.L.; Nastou, K.C.; Lyon, D.; Kirsch, R.; Pyysalo, S.; Doncheva, N.T.; Legeay, M.; Fang, T.; Bork, P.; et al. The STRING database in 2021: Customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021, 49, D605–D612. [Google Scholar] [CrossRef]
Figure 1. ATG5 K128R mutation abrogates conjugation to ATG12 and formation of the complete complex. (A). A schematic representation of the two complexes used to identify putative autophagy-unrelated interactors of ATG5: the complete complex encompassing ATG5-ATG12/ATG16 and a minimal complex comprising ATG5-ATG16. The K128R point mutation was predicted to impede covalent conjugation of ATG5 to ATG12 leading to the formation of only minimal ATG5-ATG16 complex. The complete complex is essential for autophagy, while the minimal complex is more likely to have autophagy-unrelated functions. (B). Western blot analysis demonstrating the absence of ATG5-ATG12 conjugation in Arabidopsis transgenic line expressing ATG5K128R -TAP. α-actin was used for protein loading control. (C). The ATG5K128R mutant fails to restore autophagy in ATG5-deficient cells. Confocal microscopy demonstrating accumulation of autophagic bodies (white arrowheads) in the vacuole of Arabidopsis protoplasts. Protoplasts were isolated from atg5-1 Arabidopsis plants expressing GFP-ATG8a and transformed with plasmids encoding ATG5-TAP, ATG5K128R-TAP, ATG5-TagRFP, and ATG5K128R-TagRFP. To induce autophagy, protoplasts were treated with 5µM AZD and 1 µM ConA for 24 h prior to imaging. Scale bars, 10 µm. (D). Representative pictures of two-month-old Arabidopsis plants under normal growth conditions (16 h 150 μM light, 22 °C). The ATG5-TAP fusion protein is expressed under the native ATG5 promoter complements the autophagy-deficient phenotype of the atg5-1 mutant, unlike the ATG5K128R-TAPa mutant.
Figure 1. ATG5 K128R mutation abrogates conjugation to ATG12 and formation of the complete complex. (A). A schematic representation of the two complexes used to identify putative autophagy-unrelated interactors of ATG5: the complete complex encompassing ATG5-ATG12/ATG16 and a minimal complex comprising ATG5-ATG16. The K128R point mutation was predicted to impede covalent conjugation of ATG5 to ATG12 leading to the formation of only minimal ATG5-ATG16 complex. The complete complex is essential for autophagy, while the minimal complex is more likely to have autophagy-unrelated functions. (B). Western blot analysis demonstrating the absence of ATG5-ATG12 conjugation in Arabidopsis transgenic line expressing ATG5K128R -TAP. α-actin was used for protein loading control. (C). The ATG5K128R mutant fails to restore autophagy in ATG5-deficient cells. Confocal microscopy demonstrating accumulation of autophagic bodies (white arrowheads) in the vacuole of Arabidopsis protoplasts. Protoplasts were isolated from atg5-1 Arabidopsis plants expressing GFP-ATG8a and transformed with plasmids encoding ATG5-TAP, ATG5K128R-TAP, ATG5-TagRFP, and ATG5K128R-TagRFP. To induce autophagy, protoplasts were treated with 5µM AZD and 1 µM ConA for 24 h prior to imaging. Scale bars, 10 µm. (D). Representative pictures of two-month-old Arabidopsis plants under normal growth conditions (16 h 150 μM light, 22 °C). The ATG5-TAP fusion protein is expressed under the native ATG5 promoter complements the autophagy-deficient phenotype of the atg5-1 mutant, unlike the ATG5K128R-TAPa mutant.
Ijms 24 12300 g001
Figure 2. ATG5 and ATG5K128R interactomes overlap only partially. (A). Post-translational modifications identified using LC-MS/MS. Phosphorylation of Lys183 was detected for both ATG5 and ATG5K128R, while phosphorylation of Ser187 was detected for ATG5 only. Additionally, we observed phosphorylation of Ser7 for ATG12A. (B). Venn diagram displaying partially overlapping interactome networks of ATG5 and ATG5K128R. Eighty-five proteins were found in ATG5 pull-down samples only, indicating that those might be interactors of either complete complex or of ATG12. Eighty-two proteins were pulled down only with ATG5K128R indicating that they might be interacting with the individual ATG5 form, potentially the region masked in the ATG5-ATG12 conjugate. Eighty-eight proteins were shared between ATG5 and ATG5K128R pull-downs, suggesting that those might be interactors of the minimal complex. (C). STRING analysis of the selected interactors pulled down with ATG5 only, ATG5K128R only and with both types of baits. (D). List of selected proteins pulled down with ATG5 only, ATG5K128R only and with both types of baits.
Figure 2. ATG5 and ATG5K128R interactomes overlap only partially. (A). Post-translational modifications identified using LC-MS/MS. Phosphorylation of Lys183 was detected for both ATG5 and ATG5K128R, while phosphorylation of Ser187 was detected for ATG5 only. Additionally, we observed phosphorylation of Ser7 for ATG12A. (B). Venn diagram displaying partially overlapping interactome networks of ATG5 and ATG5K128R. Eighty-five proteins were found in ATG5 pull-down samples only, indicating that those might be interactors of either complete complex or of ATG12. Eighty-two proteins were pulled down only with ATG5K128R indicating that they might be interacting with the individual ATG5 form, potentially the region masked in the ATG5-ATG12 conjugate. Eighty-eight proteins were shared between ATG5 and ATG5K128R pull-downs, suggesting that those might be interactors of the minimal complex. (C). STRING analysis of the selected interactors pulled down with ATG5 only, ATG5K128R only and with both types of baits. (D). List of selected proteins pulled down with ATG5 only, ATG5K128R only and with both types of baits.
Ijms 24 12300 g002
Figure 3. Conjugation is dispensable for nuclear localization of Arabidopsis ATG5 and ATG12. (A). Localization of fluorescently labelled ATG5, ATG8, ATG12A, ATG12B and ATG16 transiently expressed in Nicotiana benthamiana epidermal leaf cells. White arrowheads indicate nuclei position. Scale bar, 20 mm. (B). Localization of fluorescently labelled conjugation-incompetent protein forms of ATG5 and ATG12 transiently expressed in N. benthamiana epidermal leaf cells. Both, ATG5K128R and ATG12-GFP are detectable in the nuclei, indicating that conjugation is dispensable for nuclear localization of these proteins. White arrowheads indicate nuclei position. Scale bar, 20 mm. (C). Western blot analysis of total protein extracts from N. benthamiana leaves expressing fluorescently labelled ATG5, ATG12 and ATG16 confirms formation of the covalent ATG5-ATG12 conjugate if the WT, but not K128R mutant of ATG5 is expressed, and shows no presence of the free fluorescent tag in the samples.
Figure 3. Conjugation is dispensable for nuclear localization of Arabidopsis ATG5 and ATG12. (A). Localization of fluorescently labelled ATG5, ATG8, ATG12A, ATG12B and ATG16 transiently expressed in Nicotiana benthamiana epidermal leaf cells. White arrowheads indicate nuclei position. Scale bar, 20 mm. (B). Localization of fluorescently labelled conjugation-incompetent protein forms of ATG5 and ATG12 transiently expressed in N. benthamiana epidermal leaf cells. Both, ATG5K128R and ATG12-GFP are detectable in the nuclei, indicating that conjugation is dispensable for nuclear localization of these proteins. White arrowheads indicate nuclei position. Scale bar, 20 mm. (C). Western blot analysis of total protein extracts from N. benthamiana leaves expressing fluorescently labelled ATG5, ATG12 and ATG16 confirms formation of the covalent ATG5-ATG12 conjugate if the WT, but not K128R mutant of ATG5 is expressed, and shows no presence of the free fluorescent tag in the samples.
Ijms 24 12300 g003
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Elander, P.H.; Holla, S.; Sabljić, I.; Gutierrez-Beltran, E.; Willems, P.; Bozhkov, P.V.; Minina, E.A. Interactome of Arabidopsis ATG5 Suggests Functions beyond Autophagy. Int. J. Mol. Sci. 2023, 24, 12300. https://doi.org/10.3390/ijms241512300

AMA Style

Elander PH, Holla S, Sabljić I, Gutierrez-Beltran E, Willems P, Bozhkov PV, Minina EA. Interactome of Arabidopsis ATG5 Suggests Functions beyond Autophagy. International Journal of Molecular Sciences. 2023; 24(15):12300. https://doi.org/10.3390/ijms241512300

Chicago/Turabian Style

Elander, Pernilla H., Sanjana Holla, Igor Sabljić, Emilio Gutierrez-Beltran, Patrick Willems, Peter V. Bozhkov, and Elena A. Minina. 2023. "Interactome of Arabidopsis ATG5 Suggests Functions beyond Autophagy" International Journal of Molecular Sciences 24, no. 15: 12300. https://doi.org/10.3390/ijms241512300

APA Style

Elander, P. H., Holla, S., Sabljić, I., Gutierrez-Beltran, E., Willems, P., Bozhkov, P. V., & Minina, E. A. (2023). Interactome of Arabidopsis ATG5 Suggests Functions beyond Autophagy. International Journal of Molecular Sciences, 24(15), 12300. https://doi.org/10.3390/ijms241512300

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop