An oomycete effector subverts host vesicle trafficking to channel starvation-induced autophagy to the pathogen interface

Eukaryotic cells deploy autophagy to eliminate invading microbes. In turn, pathogens have evolved effector proteins to counteract antimicrobial autophagy. How adapted pathogens co-opt autophagy for their own benefit is poorly understood. The Irish famine pathogen Phytophthora infestans secretes the effector protein PexRD54 that selectively activates an unknown plant autophagy pathway that antagonizes antimicrobial autophagy at the pathogen interface. Here, we show that PexRD54 induces autophagosome formation by bridging vesicles decorated by the small GTPase Rab8a with autophagic compartments labeled by the core autophagy protein ATG8CL. Rab8a is required for pathogen-triggered and starvation-induced but not antimicrobial autophagy, revealing specific trafficking pathways underpin selective autophagy. By subverting Rab8a-mediated vesicle trafficking, PexRD54 utilizes lipid droplets to facilitate biogenesis of autophagosomes diverted to pathogen feeding sites. Altogether, we show that PexRD54 mimics starvation-induced autophagy to subvert endomembrane trafficking at the host-pathogen interface, revealing how effectors bridge distinct host compartments to expedite colonization.


Introduction
Autophagy is a conserved eukaryotic cellular process that mediates the lysosomal degradation and relocation of cellular cargoes within double-membraned vesicles called autophagosomes (He and Klionsky, 2009;Leidal et al., 2020). Although previously considered to be a bulk catabolic pathway tasked with maintaining cellular homeostasis under normal or stress conditions, it is now clear that autophagy can be highly selective (Zaffagnini and Martens, 2016). Autophagic cargoes are typically captured during autophagosome formation, a complex process regulated by a set of conserved autophagy-related proteins (ATG) as well as specialized autophagy adaptors and cargo receptors (Mizushima et al., 2011). These captured cargoes are sorted within the autophagosome during maturation of the isolation membrane (also known as the phagophore) via the specific interactions between cargo receptors and ATG8, which decorates the phagophore to serve as a docking platform for cargo receptors (Slobodkin and Elazar, 2013;Ryabovol and Minibayeva, 2016).
The source of the phagophore is still under debate, but its primary source is considered to be the endoplasmic reticulum (ER) (Bernard and Klionsky, 2013;Hamasaki et al., 2013). As the cargo is captured, the phagophore undergoes massive expansion and is finally sealed to form a mature autophagosome. Therefore, formation of the autophagosome requires additional lipid supplies that are needed for elongation and final sealing of the phagophore. Supporting this view, the essential autophagy protein ATG2 was recently discovered to have lipid transfer activity (Valverde et al., 2019;Maeda et al., 2019;Osawa et al., 2019). To cope with cellular starvation, cells can rapidly generate hundreds of autophagosomes, conceivably requiring an available supply of lipids. Remarkably, in yeast, lipids mobilized from lipid droplets (LDs) were found to fuel autophagosome biogenesis during starvation-induced autophagy, which employs relatively larger autophagosomes (Shpilka et al., 2015). In contrast, smaller autophagosomes of the cytosol-to-vacuole targeting (Cvt) pathway do not rely on LDs, suggesting that LDs are specifically recruited for starvation-induced autophagy in order to meet the increased demand of lipids required for the biogenesis of larger sized autophagosomes (Shpilka et al., 2015;Dupont et al., 2014).
Although poorly characterized, there is accumulating evidence for autophagosome maturation relying on vesicle transport and membrane expansion events which are regulated by secretory pathways involving Rab GTPases (Ao et al., 2014). For instance, the mammalian Rab8a that regulates polarized secretion and lipid droplet fusion events is also implicated in autophagy (Bansal et al., 2018;Vaibhava et al., 2012). The model plant Arabidopsis genome contains five closely related Rab8a members organized into the RabE group (RabE1a-d), but whether any of the RabE1 members contribute to autophagy is unknown (Rutherford and Moore, 2002). Molecular mechanisms governing autophagosome biogenesis, including the sources of membrane precursors required for autophagosome elongation and the transport routes to position these lipid supplies at autophagosome assembly sites remain to be comprehensively described in plants.
eLife digest With its long filaments reaching deep inside its prey, the tiny fungi-like organism known as Phytophthora infestans has had a disproportionate impact on human history. Latching onto plants and feeding on their cells, it has caused large-scale starvation events such as the Irish or Highland potato famines.
Many specialized proteins allow the parasite to accomplish its feat. For instance, PexRD54 helps P. infestans hijack a cellular process known as autophagy. Healthy cells use this 'self-eating' mechanism to break down invaders or to recycle their components, for example when they require specific nutrients. The process is set in motion by various pathways of molecular events that result in specific sac-like 'vesicles' filled with cargo being transported to specialized compartments for recycling. PexRD54 can take over this mechanism by activating one of the plant autophagy pathways, directing cells to form autophagic vesicles that Phytophthora could then possibly use to feed on or to destroy antimicrobial components. How or why this is the case remains poorly understood.
To examine these questions, Pandey, Leary et al. used a combination of genetic and microscopy techniques and tracked how PexRD54 alters autophagy as P. infestans infects a tobacco-related plant. The results show that PexRD54 works by bridging two proteins: one is present on cellular vesicles filled with cargo, and the other on autophagic structures surrounding the parasite. This allows PexRD54 to direct the vesicles to the feeding sites of P. infestans so the parasite can potentially divert nutrients. Pandey, Leary et al. then went on to develop a molecule called the AIM peptide, which could block autophagy by mimicking part of PexRD54.
These results help to better grasp how a key disease affects crops, potentially leading to new ways to protect plants without the use of pesticides. They also shed light on autophagy: ultimately, a deeper understanding of this fundamental biological process could allow the development of plants which can adapt to changing environments.
The discovery of various autophagy cargo receptors uncovered a multitude of selective autophagy pathways implicated in crucial cellular functions ranging from development to immunity in both plants and animals (Zaffagnini and Martens, 2016). For instance, the plant selective autophagy cargo receptor Joka2/NBR1 mediates antiviral immunity by eliminating viral components (Hafrén et al., 2017;Hafrén et al., 2018;Jung et al., 2020). Joka2/NBR1 is also required for immunity against bacteria and oomycete pathogens, however, the extent to which defense-related autophagy acts against these pathogens is unknown (Dagdas et al., 2018;Ü stün et al., 2018). Consistent with the important role of autophagy in plant immunity, adapted pathogens have evolved strategies to manipulate the host autophagy machinery Ü stün et al., 2018;Leary et al., 2019).
Plants detect pathogens via immune sensors consisting of surface localized pattern recognition receptors (PRRs) and intracellular nucleotide-binding leucine-rich repeat-containing proteins (NLRs). In response, pathogens secrete an arsenal of effector proteins to modulate host immunity and processes to support their virulence. Interestingly, effectors function not only to evade and suppress host immunity but also to mediate nutrient uptake (Win et al., 2012). Some filamentous plant pathogens, including the Irish potato famine pathogen Phytophthora infestans, produce hyphal extensions called haustoria that grow into the host cells to facilitate effector delivery and gain access to host nutrients (Panstruga and Dodds, 2009). A haustorium is a specialized infection structure that remains enveloped by an enigmatic host-derived membrane known as the extra-haustorial membrane (EHM), whose functions and biogenesis are poorly understood (Whisson et al., 2016). Notably, we previously showed that Joka2-mediated defense-related autophagy is diverted to the EHM during P. infestans infection (Dagdas et al., 2018). The pathogen counteracts this by deploying PexRD54, a host-translocated RXLR class of effector with five consecutive WY motifs, that targets plant autophagy . PexRD54 carries a canonical C-terminal ATG8 interacting motif (AIM) that is typically found on autophagy cargo receptors to bind ATG8 . Among the diverse set of potato ATG8 members (Kellner et al., 2017;Zess et al., 2019), PexRD54 preferentially binds the ATG8CL isoform and outcompetes Joka2/NBR1 from ATG8CL complexes, thereby disarming defense-related autophagy at the pathogen interface Dagdas et al., 2018). Intriguingly, PexRD54 does not fully shutdown autophagy as has been shown for animal pathogens that suppress autophagy (Choy et al., 2012;Kimmey and Stallings, 2016;Real et al., 2018;Xu et al., 2019). Instead, it stimulates formation of autophagosomes that accumulate around the pathogen interface (Dagdas et al., 2018). How PexRD54 stimulates autophagy and in what way the pathogen benefits from this remains unknown.
Here, we show that PexRD54 mimics carbon starvation-induced autophagy by coupling the host vesicle transport regulator Rab8a (Speth et al., 2009;Zheng et al., 2005;Pfeffer, 2017) to autophagosome biogenesis at the pathogen interface. Unlike PexRD54 which activates autophagy by recruiting Rab8a to ATG8CL compartments, PexRD54's AIM peptide fails to associate with Rab8a, and instead functions as an autophagy inhibitor. Thus, by using an effector protein and its peptide derivative as molecular tools to perturb host-autophagy, we provide insights into not only how vesicle transport processes selectively support autophagosome formation, but also how the pathogen exploits these pathways to undermine plant immunity.

Results
PexRD54-ATG8 binding is not sufficient for stimulation of autophagosome formation We have previously shown that PexRD54 stimulates formation of ATG8CL-autophagosomes in a macroautophagy dependent manner, and PexRD54 itself is a substrate of host autophagy when expressed in plant cells Maqbool et al., 2016). During infection, these pathogen-induced autophagosomes are diverted to the EHM in a process that relies on the core autophagy machinery (Dagdas et al., 2018). Because AIM-mediated binding of PexRD54 to ATG8CL is essential for the activation of autophagosome formation by the effector protein Maqbool et al., 2016), we reasoned that PexRD54 could stimulate autophagosome formation by either negative regulation of host autophagy suppressors or through positive regulation of host components to the autophagosome biogenesis sites.
To understand how PexRD54 stimulates autophagy, we used N. benthamiana, because: it is widely accepted as a solanaceous model plant in the field of plant-pathogen interactions; enables rapid functional, and cell biology assays using Agrobacterium; and can be infected by P. infestans. To first address whether ATG8 binding is sufficient to trigger autophagosome induction, we generated a PexRD54 truncate comprising only the C terminal AIM peptide (amino acids 350-381, hereafter AIMp), and compared its potency to stimulate autophagosome formation to the full-length protein. Strikingly, instead of stimulating autophagosome formation, AIMp fused to RFP (RFP:AIMp) significantly reduced the number of ATG8CL-autophagosomes in leaf epidermal cells ( Figure 1A-B). Compared to RFP:GUS control, expression of RFP:AIMp reduced the number of GFP:ATG8CL-autophagosomes by~sixfold, whereas cells expressing RFP:PexRD54 had a~fourfold increase in GFP: ATG8CL-autophagosome numbers as has been shown before ( Figure 1A-B; Dagdas et al., 2016). The AIMp interacted with ATG8CL in planta ( Figure 1C-D), as was previously shown through in vitro studies Maqbool et al., 2016). However, this association appears to take place mainly in the cytoplasm as the suppression of autophagosome formation by AIMp is such that we hardly observe any GFP:ATG8CL autophagosomes ( Figure 1A-B,E). In contrast, RFP:PexRD54 and GFP:ATG8CL produced strong overlapping fluorescence signals which peak at mobile, ring-like ATG8CL-autophagosome clusters that are induced by PexRD54 as described previously . However, in cells expressing RFP:GUS control, we did not detect any RFP signal that peaks at ATG8CL-puncta ( Figure 1A-B,E). Taken together these results show that binding of PexRD54 to ATG8CL, although necessary, is not sufficient to activate autophagosome biogenesis. This suggests while the full-length protein stimulates autophagy, PexRD54's AIM peptide functions as an autophagy suppressor.

The AIM peptide of PexRD54 suppresses autophagy
To determine whether AIMp negatively regulates autophagy, we investigated its impact on autophagic flux by monitoring GFP:ATG8CL depletion over time. As there are no specific antibodies available for individual plant ATG8 isoforms, GFP:ATG8 depletion is used to measure autophagic flux in plants rather than quantifying the ratio of lipidated to unlipidated ATG8 isoforms as for other model organisms. Consistent with the AIMp triggered decrease in ATG8CL-autophagosome numbers ( Figure 1A . Western blotting showed that in the presence of RFP:AIMp, GFP:ATG8CL was still able to produce a strong protein signal even after 6 days of transient expression. On the other hand, GFP:ATG8CL protein signal was hardly detectable after just four days in the presence of RFP:PexRD54 or RFP:EV, indicating that the AIMp hampers autophagic flux ( Figure 2A). Consistent with our previous report that PexRD54 acts as a cargo receptor that activates autophagy , we saw rapid decline of RFP:PexRD54 protein levels after 2 days post infiltration (dpi) with only free RFP detectable at six dpi, indicating vacuolar processing of RFP:PexRD54 (Figure 2-figure supplement 1). This was not due lack of PexRD54 expression at later time points as we could still detect mRNA expression of all constructs which are driven by constitutive 35S promoter ( Figure 2-figure supplement 2). In line with this view, confocal microscopy revealed that while RFP:PexRD54 localizes to the cytoplasm at two dpi, RFP signal becomes vacuolar at 4 and 6 dpi ( Figure 2-figure supplement 3), providing further evidence that PexRD54 itself is a substrate of autophagy and is degraded inside the plant vacuole, consistent with the flux assays we reported earlier Maqbool et al., 2016). In contrast, RFP:AIMp remains mostly cytoplasmic at 2-6 dpi consistent with its function as a suppressor of ATG8 autophagy ( Figure 2-figure supplement 3). However, we did not observe any stabilisation of the GFP control by RFP:AIMp, RFP:PexRD54 or RFP:GUS (Figure 2-figure supplement 1), indicating that reduced turnover of GFP:ATG8CL by RFP:AIMp is specific and is not due to altered Agrobacterium-mediated expression efficiency. Moreover, the AIM peptide mutated in the conserved ATG8-interacting region was unable to prevent ATG8CL depletion ( Figure 2-figure supplement  4). Altogether, these data indicate that inhibition of autophagy by AIMp relies on ATG8 binding and AIMp can potentially block ATG8 dependent autophagy pathways.
We next investigated the extent to which AIMp acts on other potato ATG8 isoforms. RFP:AIMp showed a robust stabilization of all six potato ATG8 isoforms ( Figure 2B), demonstrating that AIMp acts as a broad spectrum autophagy suppressor. To support these results, we measured the  Figure 1 continued on next page endogenous NBR1/Joka2 and ATG8 protein levels in N. benthamiana in the presence or absence of the AIMp. Consistently, ectopic expression of RFP:AIMp led to marked increase in both native NbJoka2 and NbATG8(s) levels compared to RFP:GUS expression, whereas RFP:PexRD54 mildly increased levels of native NbJoka2 but not NbATG8s (Figure 2-figure supplement 5A). These results further support the view that PexRD54's AIMp suppresses autophagy non-selectively, whereas PexRD54 activates ATG8CL-autophagy while neutralizing Joka2-mediated autophagy as shown before . Furthermore, AIMp enhanced protein levels of the recombinantly tagged derivatives of Joka2 and ATG8CL that are expressed under constitutive 35S promoter (Figure 2-figure supplement 5B), indicating that AIMp-triggered Joka2 and ATG8CL stabilization is not due to gene expression differences caused by the AIMp.
We then explored the potency of the AIMp in autophagy suppression when applied exogenously. For this we custom synthesized PexRD54's AIM peptide (AIMp syn , 10 amino acids at the C terminus) along with an AIM peptide mutant (mAIMp syn ) that has a two amino acid substitution in the conserved residues of the AIM , fused to cell penetrating peptides. We first tested their activities in roots of a transgenic N. benthamiana line that stably express GFP:ATG8CL to test the uptake and efficacy of the peptide in different plant tissues and as many plant autophagy studies are performed in roots. Although both AIMp syn and mAIMp syn fused to 5-Carboxyfluorescein (CF-AIMp syn and CF-mAIMp syn ) were effectively taken up by the root cells (Figure 2-figure supplement 6), only the wild-type AIMp syn reduced the frequency of GFP:ATG8CL-puncta (by~10-fold) compared to mAIMp syn or buffer control ( Figure 2C-D). We also repeated these assays in leaf epidermal cells successfully, however, peptide translocation efficiency and thus autophagosome reduction by AIMp syn were much lower in leaves compared to root cells ( Figure 2E-F, Figure 2-figure supplement 7). These findings demonstrate that PexRD54's AIM peptide suppresses autophagy in a variety of tissues, likely through binding to plant ATG8 isoforms with a high affinity and limiting their access to the autophagy adaptors that are essential for induction of autophagy. This implies that full-length PexRD54 carries additional features to stimulate autophagosome formation, by for instance, recruiting and/or manipulating other host components.
PexRD54 associates with the host vesicle transport regulator Rab8a independent of ATG8CL binding We next set out to investigate the mechanism of autophagy activation by PexRD54. Although the underlying molecular mechanisms are largely unknown, autophagosome biogenesis relies on vesicle trafficking and fusion events in yeast and animals (Singh et al., 2019;Nair et al., 2011). We therefore reasoned that in addition to binding ATG8CL, PexRD54 could possibly hijack host vesicle transport machinery to stimulate autophagosome biogenesis. Interestingly, our previous proteomics survey identified Rab8a, which shows the most sequence similarity to Arabidopsis RabE1a, a member of the small Ras-related GTPases that mediate vesicle transport and fusion events, as a candidate PexRD54 interactor . We first validated PexRD54-Rab8a association through co-immunoprecipitation assays by co-expressing the potato Rab8a (herein Rab8a) or the N. benthamiana Rab8a (NbRab8a) with PexRD54 in planta ( Figure 3A, Figure 3-figure supplement 1). Notably, the AIM mutant of PexRD54 (PexRD54 AIM ) that cannot bind ATG8CL still interacted with Rab8a to a similar degree as PexRD54 ( Figure 3A), indicating that PexRD54 associates with Rab8a independent of its ATG8CL-binding activity. Consistent with this, the AIMp failed to associate Figure 1 continued data points, color indicates biological repeats. Statistical differences were analyzed by Student t-test in R. Measurements were significant when p<0.05 (*) and highly significant when p<0.0001(****). (C) In planta co-immunoprecipitation between ATG8CL and PexRD54, AIMp or an Empty vector (EV). FLAG:ATG8CL was transiently co-expressed with either RFP:PexRD54, RFP:AIMp, or RFP:EV. IPs were obtained with anti-RFP antiserum. (D) Reverse pull-down between ATG8CL and PexRD54, AIMp, or GFP:EV. IPs were obtained with anti-GFP antiserum. (C-D) Red asterisks indicate expected band sizes. Protein sizes in kDa. (E) Confocal micrographs of RFP:AIMp (top), RFP:PexRD54 (middle), or RFP:GUS (bottom), co-expressed with GFP:ATG8CL. Transects in overlay panel correspond to plot of relative fluorescence over the labeled distance. RFP:PexRD54 co-localizes in discrete puncta with GFP: ATG8CL while RFP:AIMp and RFP:GUS fluorescence signals do not peak at GFP:ATG8CL puncta. Images shown are maximal projections of 12 frames with 1 mm steps. Scale bars represent 10 mm. The online version of this article includes the following source data for figure 1: Protein sizes in kDa. (C) Confocal micrographs of transgenic N. benthamiana root cells stably expressing GFP-ATG8CL infiltrated with cell penetrating peptides or buffer. Images shown are maximal projections of 10 frames with 1 mm steps. (D) Compared to a buffer control (15, N = 16 images quantified), synthesized AIM peptide fused to a cell penetrating peptide (AIMpsyn) suppresses ATG8CL autophagosomes per image in roots (1, N = 18 images quantified), while the AIM peptide mutant mAIMp synAA does not (15, N = 17 images quantified). Statistical differences were analyzed by Student t-test in R. Measurements were significant when p<0.05 (*) and highly significant when p<0.0001(****). (E) Transgenic N. benthamiana leaf epidermal cells stably expressing GFP-ATG8CL infiltrated with cell penetrating peptides or buffer. Images shown are maximal projections of 25 frames with 1 mm steps. Scale bars represent 10 mm. (F) Scatter-boxplot shows exogenous application of cell penetrating AIMp syn in GFP-ATG8CL transgenic N. benthamiana significantly decreases the number of ATG8 puncta per image (8, N = 21 images quantified) compared to cell penetrating mAIM synAA (12, N = 22 images quantified), or Buffer control (12, N = 22 images quantified). Scattered points show individual data points, color indicates biological repeat. Statistical differences were analyzed by Students t-test in R. Measurements were significant when p<0.05 (*) and highly significant when p<0.0001(****). The online version of this article includes the following source data and figure supplement(s) for figure 2: Source data 1. Source data for western blots and autophagosome counts.      with Rab8a in pull down assays, although it still strongly interacted with ATG8CL ( Figure 3B). These results suggest that PexRD54's N-terminal region preceding the C-terminal AIM mediates Rab8a association.
We then investigated the subcellular distribution of PexRD54 and Rab8a through confocal microscopy in leaf epidermal cells. Both stably and transiently expressed Rab8a fused to GFP (GFP:Rab8a) produced fluorescent signals at both the plasma membrane and the vacuolar membrane (tonoplast) (Figure 3-figure supplement 2). In addition, GFP:Rab8a localized to mobile puncta (0.2-0.5 mm in diameter) as well as to larger ring-shaped structures (Figure 3-figure supplement 2, Video 1), indicating that Rab8a could be involved in multiple cellular trafficking events. To determine the subcellular compartment(s) where PexRD54 associates with Rab8a, we performed co-localisation experiments of GFP:Rab8a with RFP:PexRD54, RFP:AIMp or RFP:GUS. In line with the pull-down assays ( Figure 3A-B), a subset of punctate structures labeled by GFP:Rab8a showed a clear overlap with RFP:PexRD54-puncta, whereas we did not detect any RFP signal peaking at GFP:Rab8a puncta in cells expressing RFP:AIMp or RFP:GUS ( Figure 3C). Together, these results demonstrate that PexRD54 associates with Rab8a in an AIM independent manner and raise the possibility that Rab8a could be an important component of PexRD54 driven autophagy.

PexRD54 shows a higher affinity towards Rab8a-S29N mutant
Because Rab GTPases function by converting between GTP and GDP bound states, we decided to generate Rab mutants that mimic the active (GTP) and inactive (GDP bound) conformations, which are helpful for characterization of the Rab GTPase functions. Although earlier work challenged the applicability of these mutations (Langemeyer et al., 2014;Nottingham and Pfeffer, 2014), we reasoned that Rab8a mutants could still be useful to dissect the role of Rab8a in PexRD54 activated autophagy. To determine whether PexRD54 favors a particular form of Rab8a, we produced Rab8a point mutants that we presume to mimic the GTP (Rab8a Q74L ) or GDP (Rab8a S29N ) bound states and investigated their subcellular distribution (Figure 3-figure supplement 3). Unlike GFP:Rab8a, which predominantly labeled the plasma membrane, GFP:Rab8a S29N mutant showed an even distribution at the plasma membrane and the tonoplast (Figure 3-figure supplement 2B-C). In addition, both GFP:Rab8a S29N and GFP:Rab8a marked punctate structures with varying size and shape (Figure 3-figure supplement 3B-C). In contrast, GFP:Rab8a Q74L was mainly trapped in the tonoplast and showed reduced punctate distribution compared to GFP:Rab8a or GFP:Rab8a S29N (Figure 3-figure supplement 3B-D), indicating that the Q74L mutant may not be representing the fully active form of Rab8a as previously reported for other Rab GTPases (Langemeyer et al., 2014;Nottingham and Pfeffer, 2014).
We next examined the extent to which Rab8a mutants colocalize with PexRD54. When coexpressed with BFP:PexRD54, both GFP:Rab8a and GFP:Rab8a S29N consistently produced sharp fluorescence signals that overlap with the typical ring-like autophagosomes marked by PexRD54 ( Figure 3-figure supplement 4A-B). However, GFP:Rab8a Q74L showed a similar localization pattern to the GFP control, and mostly did not produce fluorescence signals that peak at BFP:PexRD54puncta (Figure 3-figure supplement 4C-D). We quantified these observations in multiple independent experiments where GFP:Rab8a and GFP:Rab8a S29N frequently (68%, N = 23) labeled BFP: PexRD54-puncta, whereas GFP:Rab8a Q74L only did so significantly less often (25%, N = 20) (Figure 3-figure supplement 4E). As an additional control, we also checked for colocalization between Rab8a mutants and PexRD54's AIM peptide. However, we did not observe any puncta co-labeled by RFP:AIMp and GFP:Rab8a or any of the Rab8a mutants we tested (Figure 3-figure supplement 5). These observations are consistent with the results that PexRD54's AIM peptide fails to associate with Rab8a ( Figure 3B) and suppresses autophagosome formation (Figures 1-2). We then compared the binding affinity of PexRD54 to Rab8a and its mutants via in planta co-immunoprecipitation. Rab8a S29N pulled-down PexRD54 more than wild type GFP:Rab8a or GFP:Rab8a Q74L in planta  In planta co-immunoprecipitation of PexRD54 with Rab8a, Rab8a Q74L , Rab8a S29N , or GFP. FLAG:PexRD54 was transiently co-expressed with GFP:Rab8a, GFP:Rab8a Q74L , GFP:Rab8a S29N , or GFP:EV. IPs were obtained with anti-GFP antiserum and total protein extracts were immunoblotted with GFP and FLAG antisera. Red asterisks indicate expected band sizes.

PexRD54 recruits Rab8a to autophagosome biogenesis sites
We next explored the potential role of Rab8a in PexRD54-triggered autophagy. To this end, we first investigated the extent to which PexRD54 associates with its two host interactors, Rab8a and ATG8CL. We achieved this through live-cell imaging of Rab8a and ATG8CL co-expressed in combination with either PexRD54, AIMp, or BFP control. This revealed that BFP:PexRD54, but not free BFP or BFP:AIMp, localizes to puncta co-labeled by RFP:ATG8CL and GFP:Rab8a ( Figure 4A-C). Furthermore, GFP:Rab8a localized to ring-shaped RFP:ATG8CL clusters triggered by BFP:PexRD54, whereas no such structures occurred in cells expressing BFP control or BFP:AIMp ( Figure 4A-C). Notably, our quantitative imaging revealed that, even in the absence of PexRD54, more than half of RFP:ATG8CL-puncta (60%, N = 18 images) are positively labeled by GFP:Rab8a ( Figure 4A,D). However, BFP:PexRD54 expression significantly increased the frequency of GFP:Rab8a-positive RFP:ATG8CL-puncta (85%, N = 18 images) ( Figure 4B,D). Conversely, we rarely detected any fluorescent puncta that were co-labeled by GFP:Rab8a and RFP:ATG8CL in the presence of BFP:AIMp (6%, N = 18), which strongly suppresses autophagosome formation ( Figure 4C-D). These results indicate that a subset of Rab8a localizes to autophagy compartments marked by ATG8CL and this is enhanced by PexRD54. Consistently, in plants stably expressing GFP:Rab8a, we observed a similar degree of PexRD54-triggered increase in ATG8CL-Rab8a colocalization in ring-shaped ATG8CL-clusters (Figure 4-figure supplement 1), suggesting that PexRD54 might boost Rab8a recruitment to autophagic compartments. To then gain biochemical evidence for PexRD54-mediated recruitment of Rab8a to ATG8CL compartments, we conducted in planta co-immunoprecipitation assays between Rab8a and ATG8CL in presence of PexRD54, AIMp, or a control. We extracted proteins at an early time point (2 dpi) to minimize differences in ATG8CL levels, especially because we express proteins that alter autophagic   activity (PexRD54 and AIMp). We observed noticeably more RFP:ATG8CL pulled down with GFP: Rab8a in presence of 3xHA:PexRD54 but not 3xHA:AIMp or 3xHA:EV ( Figure 4E). As the 3xHA tag is only 27aa and 3.3 kDa, expression of smaller constructs such as 3xHA:EV and 3xHA:AIMp are not visible on western blots. Therefore, we validated expression of the HA-tagged constructs by RT-PCR using the RNA extracts from the Agroinfiltrated leaf patches (  show variable colocalization between RFP:ATG8CL and GFP:Rab8a. Scale bars represent 10 mm (D) BFP:PexRD54 expression significantly increases punctate colocalization between RFP:ATG8CL and GFP:Rab8 (85%, N = 18 images quantified), while BFP:AIMp significantly reduces colocalization between RFP:ATG8CL and GFP:Rab8 (6%, N = 18 images quantified) compared to the BFP:EV control (60%, N = 18 images quantified). Scattered points show individual data points, colors indicate biological repeats. Statistical differences were analyzed by Welch Two Sample t-test in R. Measurements were significant when p<0.05 (*) and highly significant when p<0.0001(****). (E) In planta co-immunoprecipitation between Rab8A and ATG8CL, and PexRD54 or AIMp. GFP:Rab8A and RFP:ATG8CL were transiently co-expressed with either 3xHA:EV, 3xHA:PexRD54, or 3xHA:AIMp. IPs were obtained with anti-GFP antiserum. Total protein extracts were immunoblotted. Red asterisks indicate expected band sizes. Protein sizes in kDa. The online version of this article includes the following source data and figure supplement(s) for figure 4: Source data 1. Source data for western blots, colocalization data sheets and image analysis plugin.   Altogether, these results show that PexRD54 enhances Rab8a accumulation at ATG8CLautophagosomes.

Rab8a is required for PexRD54-triggered autophagy
We next investigated whether Rab8a is required for PexRD54-mediated autophagy. We measured the impact of NbRab8a silencing on autophagy by quantifying the RFP:ATG8CL-autophagosome numbers. In the N. benthamiana genome, we identified at least four genes encoding full-length Rab8a like proteins (NbRab8a1-4). We first decided to generate a RNA interference (RNAi) construct (RNAi:NbRab8a 1-2 hereafter) that can target the three prime untranslated regions (UTR) of NbRab8a1 and NbRab8a2, the two closest homologs of the potato Rab8a found in N. benthamiana. RNAi:NbRab8a 1-2 showed specific silencing of the NbRab8a1-2 but not NbRab8a3 and NbRab8a4 ( Figure 5-figure supplement 1). In the absence of PexRD54, silencing of NbRab8a1-2 did not alter the number of RFP:ATG8CL puncta ( Figure 5-figure supplement 2). However, following stimulation of autophagy by transient expression of GFP:PexRD54, the number of RFP:ATG8CLpuncta/cell in RNAi:NbRab8a 1-2 background reduced by half compared to a RNAi:GUS control ( Figure 5-figure supplement 2). This suggests that simultaneous knockdown of NbRab8a1 and NbRab8a2 does not affect basal autophagy, but negatively impacts PexRD54triggered autophagy. To validate these results, we set up a complementation assay in which we silenced NbRab8a1-2 in transgenic N. benthamiana lines stably expressing the GFP tagged potato Rab8a, which evades RNA silencing because it lacks the three prime UTR targeted by the RNAi:NbRab8a 1-2 construct. Consistent with the results obtained in Figure 5-figure supplement 2, we detected greater than twofold decrease in the number of HA:PexRD54 triggered RFP:ATG8CL-puncta upon delivery of RNAi:NbRab8a 1-2 construct in wild-type plants compared to RNAi:GUS ( Figure 5A We then set up a genetic complementation assay using stable transgenic N. benthamiana lines expressing the potato GFP:Rab8a protein that is resistant to Video 2. Mobile PexRD54/Rab8a-positive puncta comigrate with ATG9 vesicles. RFP:Rab8a, BFP:PexRD54, and ATG9-GFP are co-expressed three via agroinfiltration in N. benthamiana leaf epidermal cells. Confocal laser scanning microscopy was used to monitor PexRD54/Rab8a-labeled vesicles and ATG9 compartments 3 days post infiltration. The movie represents time-lapse of 28 frames acquired during 4 min 12 s (Frame interval: 9 s). https://elifesciences.org/articles/65285#video2 silencing by RNAi:NbRab8a 1-2 construct. In stable transgenic plants expressing the silencing resistant potato GFP:Rab8a protein, RNAi:NbRab8a 1-2 did not change the number of RFP:ATG8CL puncta with or without HA:PexRD54, compared to cells that express RNAi:GUS control ( Figure 5C We next generated a hairpin-silencing construct (RNAi:NbRab8a1-4) that targets all four Rab8a members in N. benthamiana. The RNAi:NbRab8a1-4 construct substantially silenced NbRab8a1 and NbRab8a3 while silencing NbRab8a2 and NbRab8a4 to a lesser extent (Figure 5-figure supplement 4). RNAi:NbRab8a1-4 did not however silence an unrelated Rab GTPase family member Rab11 ( Figure 5-figure supplement 4). Intriguingly, knockdown of the four Rab8a isoforms using the (11, N = 26 images quantified) compared to GUS silencing control (30, N = 26 images quantified), but in the absence of PexRD54 silencing Rab8a1-2 has no effect on endogenous autophagosome number (10, N = 27 images quantified) compared to silencing control (9, N = 24 images quantified). Statistical differences were analyzed by Welch Two Sample t-test in R. Measurements were significant when p<0.05 (*) and highly significant when p<0.0001(****).          RNAi:NbRab8a 1-4 construct significantly reduced basal ATG8CL ( Figure 5-figure supplement 5) and PexRD54-triggered ATG8CL autophagosome numbers ( Figure 5-figure supplement 6). This suggests a potential redundancy in Rab8a function in basal autophagy. Alternatively, NbRab8a3-4 might be involved in basal autophagy, whereas NbRab8a1-2 are not, which needs to be tested in future.
To gain further genetic evidence for Rab8a's positive role in PexRD54-triggered autophagosome formation, we used the dominant negative Rab8a mutant (N128I) (Essid et al., 2012) and measured its impact on formation of RFP:ATG8CL-autophagosomes in the presence or absence of HA: PexRD54. Consistent with the silencing assays ( Figure 5A-D), GFP:Rab8a N128I led to 50% reduction in PexRD54 triggered ATG8CL-autophagosome numbers compared to wild-type GFP:Rab8a (Figure 5-figure supplement 7). We further validated the dominant negative role of GFP:Rab8a N128I , which significantly reduced PexRD54-triggered ATG8CL puncta compared to GFP control in two independent biological replicates ( Figure 5-figure supplement 8).
Since we found that PexRD54 associates with Rab8a and its mutants with varying affinities ( Figure 3D), we next checked whether ectopic expression of Rab8a and its mutants (S29N and Q74L) have any effect on the formation of PexRD54-autophagosomes. Compared to GFP control, GFP:Rab8a expression led to a slight increase (~1.5 fold) in the number of BFP:PexRD54 puncta (Figure 5-figure supplement 9), suggesting that Rab8a could positively regulate autophagosome formation. Expression of GDP-bound GFP:Rab8a S29N substantially enhanced (~three fold) the frequency of BFP:PexRD54 puncta compared to a GFP control, whereas GTP bound GFP:Rab8a Q74L did not lead to any significant changes in the number of BFP:PexRD54 puncta compared to the GFP control ( Figure 5-figure supplement 9). These results are consistent with the pulldown assays, which revealed stronger interaction between PexRD54 and Rab8a S29N . These findings were initially surprising given the general view that Rab S-to-N mutations in this position lead to less active (or inactive) forms that mimic the GDP-bound state, whereas the Q-to-L mutations in this site are assumed to be locked in GTP-bound state that is more active. However, there are reports which revealed that these mutations cannot be generalized (Langemeyer et al., 2014;Nottingham and Pfeffer, 2014). Consistent with this view, the GDP-bound form of the mammalian Rab8a (T22N mutant) was found to promote lipid droplet (LD) fusions, indicating that S-to-N mutation in this Rab protein functions differently. Therefore, further biochemical evidence is required to determine whether these mutations show perturbed GTPase activities. Nevertheless, together with the data presented in Figure 3D, these findings demonstrate that PexRD54-driven autophagy requires Rab8a.
Rab8a is specifically recruited to PexRD54-autophagosomes and is dispensable for Joka2-mediated autophagy To better characterize the autophagy pathway stimulated by PexRD54, we further investigated the interplay between Rab8a and ATG8CL. The weak interaction of ATG8CL and Rab8a in the absence of PexRD54 ( Figure 4E) suggests for an indirect association potentially mediated through a host autophagy adaptor. Therefore, we explored whether increased ATG8CL and Rab8a association triggered by PexRD54 is a general hallmark of autophagy activation or is a process that is stimulated through plant selective autophagy adaptors. Because the plant autophagy cargo receptor Joka2 that mediates aggrephagy also binds ATG8CL and stimulates autophagosome formation Jung et al., 2020), we tested if Joka2 could interact with Rab8a and enhance Rab8a-ATG8CL association. Unlike PexRD54, Joka2 did not colocalize or interact with Rab8a ( Figure 6A-B). Moreover, our quantitative imaging revealed that Joka2 overexpression leads to a reduction of RFP:ATG8CL puncta positively labeled by GFP:Rab8a ( Figure 6C,D). This sharply contrasts with the positive impact of PexRD54 on ATG8CL-Rab8a association ( Figure 6C,D), indicating that the autophagy pathway mediated by Joka2 is different from the PexRD54 triggered autophagy, and possibly does not require Rab8a function. Supporting this, we did not detect any difference in formation of Joka2-triggered autophagosomes upon NbRab8a1-2 silencing compared to GUS silencing ( Figure 6E-F, Figure 6-figure supplement 1). Collectively, these results indicate that Joka2-mediated autophagy pathway does not involve Rab8a, and the weak association between ATG8CL and Rab8a observed in the absence of PexRD54 is not mediated by Joka2 but potentially through an unknown autophagy adaptor.

PexRD54 triggers autophagy that is reminiscent of carbon-starvationinduced autophagy
Since we found that Joka2-mediated aggrephagy pathway does not necessarily rely on Rab8a (Figure 6), we decided to test whether other plant autophagy pathways employ Rab8a. Because autophagy can be induced through carbon starvation (Huang et al., 2019) and recent studies revealed a link between Rab8a, lipid droplets (LDs) and autophagy induced by carbon starvation in different systems (Shpilka et al., 2015;Fan et al., 2019;Wu et al., 2014), we tested whether Rab8a-ATG8CL association is altered during autophagy activation following light restriction. We detected a slight  yet significant increase in the number of RFP:ATG8CL puncta upon incubation of plants for 24 hr in the dark compared to normal light conditions ( Figure 7A-B). Consistently, we detected slightly lower levels of endogenous ATG8s and transiently expressed RFP:ATG8CL following 24 hr dark treatment (Figure 7-figure supplement 1). However, when GFP:PexRD54 is present, we did not measure any further enhancement of RFP:ATG8CL puncta following 24 hr in the dark, suggesting            that PexRD54-mediated autophagy can override or mask starvation-induced autophagy ( Figure 7A-B). Furthermore, in plants stably expressing GFP:Rab8a that are light restricted, we detected an increased degree of colocalization between RFP:ATG8CL and GFP:Rab8a ( Figure 7C-D, Figure 7figure supplement 2), in a similar fashion to enhanced ATG8CL-Rab8a association mediated by PexRD54 ( Figure 4A-D). Collectively, these data suggest that PexRD54 mimics carbon-starvationinduced autophagy.
Recent studies have revealed that lipid droplets (LDs) contribute to carbon-starvation-induced autophagy during light restriction (Shpilka et al., 2015;Fan et al., 2019). In addition, the GDPbound mutant form of the mammalian Rab8a is enriched at LD contact sites to regulate their fusion . Therefore, we investigated Rab8a association with LDs under normal or starvation conditions. We first checked co-localization between GFP:Rab8a and LDs marked by the orange-red fluorescent fatty acid (FA), BODIPY 558/568 C12 (herein BODIPY-C 12 ). Confocal microscopy revealed that a small fraction of GFP:Rab8a puncta are labeled by the LD marker BODIPY-C 12 under normal light conditions, whereas the frequency of this colocalization increased by~1.4-fold when plants are light restricted ( Figure 7E-F, Figure 7-figure supplement 3). Furthermore, we observed that GFP:Rab8a puncta positive for BODIPY-C 12 are also labeled by BFP: ATG8CL (Figure 7-figure supplement 4). The stronger association of Rab8a and LDs upon light restriction, combined with the finding that LDs are recruited toward autophagosomes during carbon starvation (Fan et al., 2019), suggest that Rab8a-LD association could be a hallmark of starvation-induced autophagy.
Additionally, we observed that ring-shaped PexRD54 clusters labeled with Rab8a also tightly associate with BODIPY-C 12 labeled puncta ( Figure 7G, Figure 7-figure supplements 5-6). However, we did not detect BODIPY-C 12 fluorescence signal filling the lumen of the PexRD54-labeled compartments ( Figure 7G, Figure 7-figure supplements 5-6), indicating that fatty acids are likely not the autophagic cargoes of PexRD54. Rather, we detected a BODIPY-C 12 signal at the periphery of autophagosomes marked by PexRD54, which overlaps with GFP:Rab8a fluorescence signal ( Figure 7G, Figure 7-figure supplements 5-6). Moreover, these PexRD54/Rab8a-clusters are also accompanied by LDs densely labeled with only BODIPY-C 12 as they navigate through the cytoplasm (Video 3). These findings suggest that fatty acids could be one of the potential membrane sources of the autophagosomes stimulated by PexRD54 as observed during carbonstarvation-induced autophagy in other systems (Shpilka et al., 2015). To gain further evidence for this, we investigated the colocalization of PexRD54 and Rab8a with the LD structural membrane protein Oleosin (Siloto et al., 2006;Fan et al., 2019;Singh et al., 2009). Similar to BODIPY-C 12 -labeled puncta (Figure 7-figure supplement 5), Oleosin labeled LDs clustered around PexRD54/Rab8a-positive ring-like autophagosomes (Figure 7-figure supplement 6). Although Oleosin-positive LDs were adjacent to PexRD54 autophagosomes, in contrast to BOD-IPY-C 12 , Oleosin-YFP did not produce fluorescent signal that overlaps with PexRD54/Rab8a ring-like autophagosomes (Figure 7-figure supplement 6), suggesting that FAs but not LD surface proteins are transferred to PexRD54-triggered autophagosomes. Strikingly, stimulation of autophagy by BFP:PexRD54, but not Joka2: BFP, led to enhanced association of BODIPY-C 12 and GFP:Rab8a-labeled puncta, supporting the hypothesis that PexRD54 mimics autophagy induced during carbon starvation ( Figure 7H, Figure 7-figure supplement 7). Together these results show that unlike the aggrephagy receptor Joka2, PexRD54 triggers responses similar to carbon-starvation-induced autophagy including induction of autophagosomes and enhanced association of ATG8CL-autophagosomes with Rab8a and LDs ( Figure 6 and Figure 7). However, it is possible that PexRD54 could activate other host autophagy pathways that favor pathogen virulence.
To investigate Rab8a's role in recruiting lipid droplets to PexRD54 puncta, we knocked down Rab8a with the RNAi:NbRab8a 1-4 construct and quantified the proportion of PexRD54 autophagosomes associated with Oleosin-labeled LDs. We observed that Rab8a knockdown not only reduces the total amount of PexRD54 autophagosomes but also the proportion of those associated with LDs compared to a silencing control (Figure 7-figure supplement 8). Finally to gain further evidence of the positive contribution of LDs to PexRD54-triggered autophagosomes, we designed a RNAi silencing construct to knockdown Seipin that has been shown to contribute to LD biogenesis in yeast, animals, and plants (Cai et al., 2015;Yang et al., 2012;Taurino et al., 2018;Greer et al., 2020). RNAi-mediated knockdown of the NbSeipin-A significantly reduced the amount of Oleosin-labeled LDs compared to a silencing control (Figure 7-figure supplements 9-10) consistent with the established role of Seipin in LD biogenesis. NbSeipin-A knockdown significantly reduced the amount of PexRD54 puncta in the cell compared to the control (Figure 7-figure supplement 11), supporting the notion that PexRD54 relies on host LD resources to stimulate autophagosome formation. Together, these data show that PexRD54 requires Rab8a to recruit lipid droplets and that LDs contribute to its induction of autophagosome formation. These findings are in agreement with the findings in yeast and mammalian cells that LDs supply membrane sources required for autophagosome formation during starvation-induced autophagy autophagosomes (Shpilka et al., 2015;Dupont et al., 2014).

PexRD54 subverts Rab8a to autophagosomes at the pathogen interface
Our recent work revealed that the perihaustorial niche is a hot spot for the formation of ATG8CL autophagosomes stimulated by PexRD54 Dagdas et al., 2018). Therefore, we next examined whether Rab8a-PexRD54 association occurs at perihaustorial ATG8CL-autophagosomes. We first checked GFP:Rab8a localization alone in the haustoriated cells. In infected leaf epidermal cells transiently or stably expressing GFP:Rab8a, we detected varying sizes of GFP:Rab8a puncta around the P. infestans haustoria (Figure 8-figure supplement 1). These structures included ring shaped compartments that are reminiscent of PexRD54autophagosomes as well as smaller densely packed GFP-positive puncta and large vacuole like structures, indicating that Rab8a could regulate diverse trafficking pathways during infection (Figure 8-figure supplement 1, Videos 4 -5). To verify that the perihaustorial Rab8a puncta represent the PexRD54-autophagosomes, we imaged infected plant cells which co-express GFP:Rab8a and the autophagosome marker protein RFP:ATG8CL in combination with BFP: PexRD54, BFP:AIMp, BFP, or Joka2:BFP. Confocal micrographs of haustoriated plant cells showed accumulation of RFP:ATG8CL-autophagosomes around the haustoria which are colabeled with GFP:Rab8a, and are positive for BFP:PexRD54 but not BFP control ( Figure 8A-C). However, formation of perihaustorial puncta co-labeled by RFP:ATG8CL and GFP:Rab8a was suppressed by arresting autophagosome formation through expression of BFP:AIMp ( Figure 8B). Notably, we detected a subset of perihaustorial GFP:Rab8a puncta that are not labeled by RFP:ATG8CL, further supporting an ATG8CL-independent haustorial trafficking role for Rab8a ( Figure 8A-D, Figure 8-figure supplement 2). In agreement with these findings, other Rab8 members including RabE1 family have been reported to localize to distinct subcellar compartments including Golgi and Peroxisomes in Arabidopsis (Zheng et al., 2005;Cui et al., 2013). Therefore, we next investigated the potential colocalization of Rab8a and PexRD54 with markers that label organelles such as Golgi, peroxisomes and mitochondria. Coexpression of Rab8a and PexRD54 with GmMan1 1-49 -mCherry (Golgi marker), ScCOX4 1-29 -mCherry (mitochondria marker), or CFP-GLOX (peroxisome marker) revealed extensive co-localization of Rab8a with the Golgi marker but not with the markers that label mitochondria or peroxisomes ( Figure  8-figure supplement 3A; Nelson et al., 2007). However, PexRD54-Rab8a puncta did not show any labeling by any of the organelle markers tested (Figure 8-figure supplement 3A). Consistent with these findings, Rab8a distinctly localized at either the Golgi or PexRD54 puncta, but not at peroxisomes or mitochondria, at the haustorium interface (Figure 8-figure supplement 3B). Therefore, we conclude that Rab8a localization to PexRD54 autophagosomes cannot be explained by Rab8a's association with these organelles. Intriguingly, we noticed that peroxisomes frequently appear around the PexRD54-Rab8a puncta, suggesting a potential involvement of this organelle in PexRD54-mediated autophagy (Figure 8-figure supplement 3).
We also observed that Rab8a and PexRD54 co-localize with both Bodipy C-12 and oleosin in clusters of vesicles around haustoria, linking Rab8a's emerging role in lipid trafficking with the pathogen effector and haustorial interface (Figure 8-figure supplement 4). On the other hand, in line with our findings that the Joka2 pathway does not employ Rab8a (Figure 6), perihaustorial Joka2:BFP/ RFP:ATG8CL puncta and GFP:Rab8a puncta were exclusive to each other ( Figure 8D). Consistent with our pull down assays ( Figure 3D), we did not detect any sharp GFP:Rab8a Q74L signal at the perihaustorial BFP:PexRD54 puncta (Figure 8-figure supplement 5). In contrast, GFP:Rab8a, and particularly GFP:Rab8a S29N , produced strong fluorescence signals peaking at perihaustorial BFP: PexRD54 puncta (Figure 8-figure supplement 5), indicating that both wild-type Rab8a and Rab8a S29N are enriched at the perihaustorial PexRD54 autophagosomes. Taken together, these results demonstrate that Rab8a localizes to distinct compartments that accumulate around the haustorium and PexRD54 stimulates diversion of Rab8a positive LDs to perihaustorial autophagosomes.
Localization of Rab8a to Golgi around the haustorium is in agreement with the conserved role of Rab8 in yeast and metazoans to mediate polarized secretion of proteins (Nielsen et al., 2008). These findings, combined with our data that PexRD54-mediated subversion of Rab8a to autophagosomes (Figure 8) around the haustoria, imply a potential function of Rab8a in polarized defense-responses. Therefore, we next investigated the possible role of Rab8a in immunity against P. infestans. Simultaneous RNAi-mediated knockdown of all four NbRab8a members led to a consistent increase in disease symptoms and P. infestans hyphal growth (Figure 8-figure  supplement 6). To further determine the role of the Rab8a family in immunity, we conducted a silencing complementation assays using a codon shuffled NbRab8a-1 construct fused to GFP    (GFP:Rab8a syn ) that can evade RNAi (Figure 8-figure supplement 7). The enhanced susceptibility phenotype caused by RNAi:NbRab8a 1-4 silencing was rescued by simultaneous overexpression of GFP:Rab8a syn but not the GFP control, providing further evidence that Rab8a is required for basal resistance against P. infestans (Figure 8-figure supplement 7). Consistent with knockdown assays, overexpression of dominant negative GFP:Rab8a N128I enhanced plant susceptibility to P. infestans compared to a GFP control (Figure 8-figure supplement 8). Collectively, these results indicate that Rab8a contributes to plant immunity, whereas PexRD54 could possibly interfere with the defense-related role of Rab8a by subverting a subpopulation of it to autophagic compartments.

AIM peptide-mediated autophagy arrest leads to reduced pathogen growth
The observed induction of haustoria targeted autophagosomes prompted the hypothesis that P. infestans could benefit by co-opting host autophagy to support its own growth. Therefore, we next explored how activation of autophagy by PexRD54 affects P. infestans host colonization.
Our discovery of the AIMp as an ATG8-specific autophagy inhibitor that can be used to spatiotemporally arrest plant autophagy allowed us to test the impact of autophagy on P. infestans virulence. To this end, we decided to transiently interfere with pathogen-induced autophagy by expressing the AIMp. Our quantitative image analysis revealed that compared to RFP:PexRD54 expression, transient expression of RFP:AIMp led to~threefold decrease in the number of haustoria that are associated with autophagosomes marked by GFP:ATG8CL ( Figure 9A-B). We then measured how autophagy suppression by the AIMp affects P. infestans infection. In multiple independent experiments (six biological replicates), N. benthamiana leaf patches expressing RFP: AIMp showed a consistent reduction of quantitative disease symptoms compared to an RFP vector control ( Figure 9C-D). This indicated that AIMp-mediated arrest of host autophagy negatively impacts P. infestans virulence, supporting the hypothesis that PexRD54-triggered autophagy is beneficial to the pathogen. Collectively, these results suggest that P. infestans relies on host autophagy function to support its virulence. This could explain why the pathogen deploys full-length PexRD54 that can activate specific host autophagy pathways while subverting defense-related autophagy, instead of just the AIM peptide.

Discussion
Dissecting the specialized functions and mechanisms of autophagy in host-microbe interactions has been challenging. This is mainly due to prolonged stress accumulation in autophagy mutants and non-autophagy related roles of the targeted genes (Munch et al., 2014). Pathogen effectors that target specific components of the host autophagy machinery have emerged as alternative tools to unravel the underlying mechanisms of defense-related autophagy Hafrén et al., 2018;Ü stün et al., 2018). Here, by studying the P. infestans effector protein PexRD54, we shed light on the poorly understood mechanism of pathogen-induced autophagy in plants. We unveil a distinct pathogen virulence mechanism in which the effector protein couples host vesicle transport machinery to autophagosome biogenesis. Through its modular domain architecture, PexRD54 employs a diverse set of host proteins such as Rab8a and ATG8CL to stimulate autophagy that is reminiscent of starvation-induced autophagy. We propose a model in which PexRD54 recruits Rab8 and lipid droplets to the pathogen feeding site to stimulate autophagosome biogenesis (Figure 10). Both PexRD54 and carbon starvation drive association of  Rab8a and LDs which are positioned at autophagy compartments marked by ATG8CL. We hypothesize that the pathogen could benefit from this process potentially through subverting the immune function of Rab8a and gaining access to plant resources carried in autophagosomes diverted to haustoria.

Autophagy suppression by the AIM peptide
Effectors are excellent tools to dissect complex biological processes such as autophagy as they often display high target specificity (Schardon et al., 2016). In this study, we not only expanded our knowledge of pathogen-triggered autophagy but also discovered an effector derived peptide (AIMp) that can specifically block autophagy. Here, the AIMp served as an excellent negative control to understand PexRd54 activities and to interpret the impact of autophagy manipulation by the pathogen. We used it as a tool to perturb PexRD54 activities without directly interfering with Rab8a functions. Most ATG mutants show pleiotropic affects and have non-autophagy-related roles, making it difficult to interpret the outcomes of autophagy inhibition. However, the AIMp blocks autophagy specifically, as it directly acts on ATG8. Chemical inhibitors are also often used to measure autophagy flux. However, these inhibitors are mostly inefficient and lack the required specificity. We discovered that PexRD54's AIM peptide is a strong autophagy suppressor effective against all potato ATG8 isoforms ( Figure 2). Our data suggest that AIMp blocks the autophagosome biogenesis step, as both autophagosome formation and autophagic flux are suppressed by the AIMp (Figures 1-2). Conceivably, the AIM peptide competitively inhibits autophagosome biogenesis by occupying the AIM docking site of ATG8 that accommodates host autophagy regulators (Noda et al., 2008). In agreement with our finding, autophagy cargo receptors that bind ATG8 via AIMs were recently shown to stimulate ATG8 (Chang et al., 2021). The AIM peptide is a genetically encodable tool, which can enable spatio-temporal arrest of autophagy when expressed under inducible or tissue specific promoters. Thus, this should be of great interest for autophagy studies in plants and other systems, which can overcome the limitations of chemical autophagy inhibitors and autophagy mutants. We also developed AIM peptide derivatives with cell penetrating features, which allow studying the tissue-specific functions of autophagy. The cell penetrating AIM peptide can be used to study autophagy in plants and other eukaryotic systems that are not amenable to genetic manipulation.

How does PexRD54 activate autophagy?
Autophagosome biogenesis is a complex multi-step process. But how could an effector activate such an intricate process? We uncovered that PexRD54 either directly or indirectly recruits Rab8a to autophagosome biogenesis sites (Figures 3-4). Similar channeling of Rab8a to ATG8CL-autophagosomes occurs during carbon starvation (Figure 7), suggesting that PexRD54 mimics autophagy induction via nutrient deprivation. Interestingly, the mammalian autophagy cargo receptor Optineurin also interacts with both LC3 (mammalian ATG8 isoform) and Rab8a in mice (Bansal et al., 2018;Vaibhava et al., 2012). Although the functional implications of these interactions are unclear, a proposed model suggests that Optineurin mediates pre-autophagosomal membrane elongation through anchoring Rab8a to autophagosome assembly sites. Our data is consistent with this model and suggests that PexRD54 recruits Rab8a to facilitate autophagosome formation. Scattered points indicate individual data points and different colors represent various biological repeats. Statistical differences were analyzed by Welch Two Sample t-test in R. Measurements were significant when p<0.05 (*) and highly significant when p<0.001(***). The online version of this article includes the following source data for figure 9: Source data 1. Source data for infection assays.

The role of Rab8a in autophagy
The membrane elongation step of autophagosome formation relies on direct transport of lipids from various donor compartments. LDs provide a membrane source for autophagosome biogenesis specifically during starvation-induced autophagy (Shpilka et al., 2015). More recently, a conserved acyl-CoA synthetase (ACS) from yeast was shown to be mobilized on nucleated phagophores where it locally mediates transport of fatty acids required for phagophore elongation (Schütter et al., 2020). But how are these lipid sources mobilized to autophagosome biogenesis sites at first? Intriguingly, the GDP bound form of the mammalian Rab8a mediates LD fusion events . In line with these reports, we found that PexRD54 associated more strongly with the Rab8a S29N mutant (presumably the GDP bound form) and recruits Rab8a to autophagosome biogenesis sites to enhance LD localization around the ATG8CL-foci and stimulate autophagosome formation (Figures 4  and 7, Figure 4-figure supplement 3, Figure 7-figure supplements 4-10). However, we did not find any luminal signal with the FA marker, indicating that FAs are not the likely the cargoes of PexRD54. Furthermore, we detected FAs labeled with Bodipy but not the LD membrane marker protein Oleosin localized to periphery of the PexRD54 compartments ( Figure 7G, Figure 7-figure supplements 5-6). This suggests that PexRD54 could employ Rab8a to position LDs at ATG8CL nucleation sites to facilitate lipid transfer for autophagosome biogenesis. Interestingly, PexRD54triggered association of Rab8a and LDs is also enhanced by carbon starvation but not by Joka2 ( Figure 7G-H, Figure 7-figure supplement 6). We propose that upon carbon starvation, the plant uses LDs as an additional membrane source to accommodate for an increase in autophagosome biogenesis, and PexRD54 exploits this process to stimulate autophagy. Nevertheless, Rab8a may not be engaged in all autophagy routes, as it is dispensable for Joka2-mediated autophagy ( Figure 6). This is consistent with the finding that Arabidopsis NBR1 (Joka2) mutants are sensitive to a variety of abiotic stress conditions but not to carbon starvation (Zhou et al., 2013). Our results are consistent with the recent findings in yeast, where LDs specifically contribute to starvation-induced autophagy (Shpilka et al., 2015). Nevertheless, further assays are needed to determine whether Rab8a is implicated in other autophagy pathways. Altogether, combined with previous findings, we conclude that distinct cellular transport pathways feed autophagosome formation during diverse selective autophagy pathways in plants.

The role of Rab8a in immunity
Our data revealed that the Rab8a family contributes to basal resistance against P. infestans (Figure 8-figure supplements 6-8). Rab8a belongs to the Rab8 family of small GTPases that are implicated in polarized secretion events in eukaryotes (Pfeffer, 2017). Another Rab8 member known as RabE1d was found to contribute to bacterial resistance and regulate membrane trafficking in the model plant Arabidopsis (Speth et al., 2009;Zheng et al., 2005). However, the extent to which Rab8 family members function in immunity remains unknown. Our findings revealed that Rab8a is most likely involved in a diverse range of cellular transport pathways including autophagy. Consistent with the evolutionarily conserved role of Rab8 family members in polarized secretion, we detected Rab8a localization to Golgi around the haustorium. These findings suggest that Rab8a has a function in defense-related secretion that possibly contributes to plant focal immune responses. Further research is required to determine the cargoes and the potential defense-related functions of Rab8a during infection.

Why does P. infestans activate autophagy?
Autophagy suppression by the AIM peptide leads to reduced pathogen virulence (Figure 9), indicating that host autophagy inhibition is not favorable to P. infestans. This could explain why P. infestans deploys PexRD54, which can neutralize defense-related autophagy mediated by Joka2, while enabling other autophagy pathways that are somewhat beneficial. Intriguingly, the autophagy pathway primed by PexRD54 resembles autophagy induced by carbon starvation, but not to autophagy induced by the aggrephagy cargo receptor Joka2. In contrast to Joka2, both PexRD54 and light restriction stimulated enhanced Rab8a-ATG8CL and Rab8a-LD associations ( Figure 6 and Figure 7). Furthermore, their effect on autophagosome formation were not additive ( Figure 7A-B). This hints at the possibility that PexRD54 could facilitate nutrient uptake from the host cells by mimicking starvation conditions to stimulate autophagy. This view is further supported by our earlier finding that PexRD54-autophagosomes are diverted to the haustorial interface . However, our data suggests that the fate of PexRD54 autophagosomes in infected and uninfected cells could be different, since in uninfected cells PexRD54 is degraded inside the plant vacuole following 4-6 days of ectopic expression (Figure 2-figure supplements 1-3). Therefore, PexRD54 may or may not be degraded in infected plant cells. PexRD54 is a member of RXLR effectors which typically show high level of expression during 2-5 days of infection. Consequently, even if PexRD54 is depleted via host autophagy machinery during infection, there would be more of it secreted via the haustorium to maintain PexRD54 virulence functions as long as the haustorium is accommodated in the host cells. PexRD54 likely remodels the cargoes engulfed in ATG8CL-autophagosomes targeted to pathogen interface, which are presumably assimilated by the parasite. Whether these cargoes are hydrolyzed in infected plant cells or directly absorbed by the pathogen remains to be determined.
Alternatively, but not mutually exclusively, PexRD54 could also help neutralize defense-related host components by engulfing them in secure membrane-bound autophagy compartments. For instance, PexRD54 could promote diversion of Rab8a to autophagy to intervene with non-autophagy related immune functions of Rab8a ( Figure 8). Therefore, we cannot exclude the possibility that PexRD54 could activate other host autophagy pathways that are beneficial to the pathogen, and further research is needed to determine whether PexRD54 could mimic other autophagy inducing conditions. In summary, our findings demonstrate that to support its virulence, P. infestans manipulates plant cellular degradative and transport systems by deploying an effector protein that imitates carbon starvation conditions. It also demonstrates effectors can act as adaptors to bridge multiple host components to modulate complex cellular processes for the benefit the pathogen. Further research is needed (i) to determine the cargo of these autophagosomes, (ii) whether they are secreted to the pathogen interface and subsequently absorbed by the pathogen, and (iii) the molecular players involved in pathogen subverted autophagy.

Pathogenicity assays
P. infestans 88069 strain was used in this study. Cultures were grown and maintained by routine passing on rye sucrose agar medium at 18˚C in the dark (van West et al., 1998). Zoospores were collected from 10 to 14 days old culture by flooding with cold water and incubation at 4˚C for 90-120 min. Infection of agroinfiltrated leaves was carried out by addition of 10 mL droplets of zoospore solution at 50,000 spores/ml on detached N. benthamiana leaves (Chaparro-Garcia et al., 2011). Infection for microscopic experiments carried out on attached leaves. Inoculated detached leaves or plants were kept in humid conditions. Day light/UV images were taken at 7 days post infection and lesion areas were measured in ImageJ.

Molecular cloning and plasmid constructs
Various constructs used in this study were published previously. GFP:ATG8CL, GFP:PexRD54, GFP: PexRD54aim, RFP:Rem1.3 constructs were previously described in Bozkurt et al., 2015 For YFP:Oleosin, the eYFP fluorophore was split into N-terminal (residue M1-A155) and C-terminal half (residue D156 -K239). The N-terminal split YFP half was used via a linker peptide RPACKIPNDLKQKVMNH and the C-terminal split YFP half via a linker peptide HNMVKQKLDNPIKCAPR. EcoRV restriction site was added at the end of each linker to allow linearization of the vector and provide an insertion site for subsequent cloning. The DNA fragment encoding Oleosin, together with the linker peptides and restriction sites were amplified from N. benthamiana gDNA using primers NYFP-Oleosin F and Oleosin-CYFP R then assembled into pK7WGF2 vector backbone by Gibson assembly. GFP:Rab8a and RFP:Rab8a constructs were generated by PCR amplification from Solanum tuberosum cDNA using primers GW_StRab8a-1_F GW_StRab8a-1_R followed by Gateway cloning into the entry vector pENTR/D/TOPO (Invitrogen) then into the pK7WGF2 (GFP) and pH7WGR2 (RFP) vectors, respectively. RFP:GUS was created from the pENTR-GUS control plasmid provided in the GATEWAY cloning kit and inserted into pH7WGR2 (RFP) via LR reaction. Single residue mutations of Rab8aS29N, Rab8aQ74L and Rab8aN128I were obtained by inverse polymerase chain reaction (PCR) amplification of the StRab8a entry clone with the primer pairs (phosphorylated at five prime ends) carrying desired mutations; (i) Rab8aS29N_F and Rab8aS29N_R; (ii) Rab8aQ74L_F and Rab8aQ74L_R; (iii) Rab8aN128I_F and Rab8aN128I_R. Templates were then eliminated by one-hour Dpn-I (New England Biolabs) restriction digestion at 37˚C and the PCR products of mutants were ligated using standard protocols to obtain circular Gateway entry clones carrying desired mutations. Next, the entry clones of Rab8a mutants were recombined into destination vectors pK7WGF2 or pB7RWG2 by Gateway LR reaction. All remaining constructs were amplified from existing constructs previously described (Bozkurt et al., 2015;Dagdas et al., 2016;Dagdas et al., 2018), using primer pairs GA_RD54_F with GA_RD54_R for PexRD54, GA_RD54_F with GA_LIR2_R for PexRD54aim, GA_AIMp_F with GA_LIR2_R for PexRD54AIMp, GA_ATG8C_F with GA_ATG8C_R for ATG8CL and GA_NbJoka2_1_Fr with GA_NbJoka2_1_Rv. Silencing constructs for Rab8a were amplified using primer combinations NbRab8A_silF1 and NbRa-b8A_silR1, Rab8a1-4 RNAi _F1, Rab8a1-4 RNAi _F2, Rab8a1-4 RNAi _F3, Rab8a1-4 RNAi _R1, Rab8a1-4 RNAi _R2 and Rab8a1-4 RNAi _R3, and cloned into the pRNAiGG vector as described elsewhere (Dagdas et al., 2018). Silencing of Rab8a was verified using RT-PCR. CFP-GLOX (Glycolate oxidase) was generated by PCR amplification from wheat cDNA using primers GLOX-F and GLOX-R followed by Gateway cloning into the entry vector pENTR/D/TOPO (Invitrogen) then into the pGWB445 (CFP) vector.

Synthetic construct generation
The sequence for the silencing resistant synthetic Rab8a construct (Rab8a-1 syn ) was obtained by codon shuffling the sequence of Rab8a-1 through a combination of the Integrated DNA Technologies (IDT) Codon optimization tool and manual codon shuffling. The construct was then gene synthesized and inserted into the vector by Gibson assembly.

RNA isolation, cDNA synthesis, and RT-PCR
For RNA extraction, 100 mg of leaf tissue was excised and frozen in liquid nitrogen. RNA was then extracted using either GeneJET Plant RNA Purification Kit (Thermo Scientific) or TRIzol RNA Isolation Reagent (Invitrogen) according to producers' recommendations. RNA concentration was measured using NanoDrop Lite Spectrophotometer (Thermo Scientific). For Figure 2-figure supplement 2 and mg of extracted RNA underwent DNase treatment using RQ1 RNase-Free DNase (Promega). Two mg of RNA was then used for cDNA synthesis using SuperScript IV Reverse Transcriptase (Invitrogen). cDNA was then amplified using Phusion High-Fidelity DNA Polymerase (New England Biolabs) with the appropriate primer pairs (TKey Resources Table). GAPDH is used to normalize the transcription levels of the genes.

Confocal microscopy and image processing
Microscopy analyses were carried out on live leaf tissue 3-4 days post agroinfiltration. To minimize the damage of live tissue, leaf discs of N. benthamiana were cut using a cork borer and mounted onto Carolina observation gel (Carolina Biological Supply Company). For BODIPY-dodecanoic acid (BODIPY-C12, Invitrogen) staining, 10 mM was infiltrated into the leaf tissue 5 hr prior to observation. For PexRD54 AIM peptide experiments in leaf tissue, a solution of 10 mM of peptide in agroinfiltration buffer or buffer alone was infiltrated in leaves 3 hr prior to observation. For imaging in roots, seedlings were collected at 3 weeks old and the roots placed in 2 mL tubes containing 5 mM peptide solution in agroinfiltration buffer or buffer alone for 3 hr prior to observation. Confocal florescence microscopy was performed using Leica SP5 and SP8 resonant inverted confocal microscope (Leica Microsystems) using 63x and 40x water immersion objective, respectively. In order to excite fluorescent tagged proteins, Diode laser excitation was set to 405 nm, Argon laser to 488 nm and the Helium-Neon laser to 561 nm and their fluorescent emissions detected at 450-480, 495-550, and 570-620 nm to visualize BFP, GFP, and RFP fluorescence, respectively. Sequential scanning between lines was done to avoid spectral mixing from different fluorophores and images acquired using multichannel. Maximum intensity projections of Z-stack images were processed using ImageJ (2.0) to enhance image clarity.

Data analysis and statistics
Images for quantification of autophagosome numbers were obtained from Z stacks consisting of 1.3 mm depth field multi-layered images with similar settings for all samples unless stated otherwise. To detect and quantify punctate structures in one channel (green channel or red channel or blue channel) and to validate colocalization an overlay of two or three channel, where applicable, was acquired (green channel and/or red channel and/or blue channel). Z stacks were separated into individual images with the ImageJ (2.0) program and analyzed. Autopahgosome quantification was done using ATG8CL as a marker protein as described before . Colocalization of puncta was analyzed using ImageJ (2.0) with a modified version of the colocalization macro described elsewhere (Pampliega et al., 2013) with thresholding set as manual to avoid background cytoplasmic signals in each image. Boxplots were generated with mean of punctate numbers generated from stacks obtained in three to six independent biological experiments. Statistical differences were analysed by Welch Two Sample t-test in R. Measurements were significant when p<0.05 (*) and highly significant when p<0.001(***).

Automated puncta counting algorithm through image processing
The image processing algorithm calculates the gradient of the image to identify the boundaries of the puncta. We then algorithmically identify the enclosed regions formed by the boundaries and counted the number of puncta in each figure. For the case of co-localisation, the co-ordinates of the centres of the puncta/clusters from each channel were calculated and compared to see if they lie within a small tolerance for each puncta and channel. The puncta/clusters satisfying the abovementioned conditions were considered to be co-localized and were counted.