Differential and convergent utilization of autophagy components by positive-strand RNA viruses

Many viruses interface with the autophagy pathway, a highly conserved process for recycling cellular components. For three viral infections in which autophagy constituents are proviral (poliovirus, dengue, and Zika), we developed a panel of knockouts (KOs) of autophagy-related genes to test which components of the canonical pathway are utilized. We discovered that each virus uses a distinct set of initiation components; however, all three viruses utilize autophagy-related gene 9 (ATG9), a lipid scavenging protein, and LC3 (light-chain 3), which is involved in membrane curvature. These results show that viruses use noncanonical routes for membrane sculpting and LC3 recruitment. By measuring viral RNA abundance, we also found that poliovirus utilizes these autophagy components for intracellular growth, while dengue and Zika virus only use autophagy components for post-RNA replication processes. Comparing how RNA viruses manipulate the autophagy pathway reveals new noncanonical autophagy routes, explains the exacerbation of disease by starvation, and uncovers common targets for antiviral drugs.


Introduction
Many viruses interact with the autophagy pathway, which can play both antiviral and proviral roles, sometimes even during the same viral infections. Antiviral activities of autophagy include xenophagy, in which intracellular viruses are targeted for degradation as well as the triggering and facilitation of immune responses [27]. As a testament to the power of autophagy as an arm of innate immunity, many successful viruses such as herpesviruses [28], Sindbis virus [29], and vesicular stomatitis virus [30] have evolved strategies to block the induction of the autophagy pathway. On the other hand, utilization of components of the autophagy pathway by many viruses serves proviral functions. The unique ability of the autophagy pathway to generate membranous structures de novo and to allow or disallow their acidification may be desirable features for any virus that survives by manipulating intracellular membranes. All positive-strand RNA viruses associate with cytoplasmic membranes to replicate their genomes. Virally associated membranes come from a variety of cellular compartments, including ER, Golgi, mitochondria, and lysosomes [31][32][33][34]. Many viral proteins are membrane-bound, and RNA amplification occurs on the topologically cytoplasmic surfaces of vesicles. Viruses such as dengue, for instance, induce ER-membrane invaginations to form isolated pockets in which viral proteins congregate for RNA replication outside the reach of cellular antiviral factors.
To determine which components of the autophagy pathway are used by these RNA viruses and whether there are any shared principles, we generated a panel of CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)-Cas9 knockout (KO) human epithelial-derived cell lines (HeLa). Previously, most research investigating viruses and cellular autophagy has involved targeting single genes genetically or pharmaceutically. This does not take into account the potential utilization of noncanonical autophagy pathways or gain-of-function effects of drugs [50]. We found that PV, DENV, and ZIKV all utilize multiple components of the autophagy pathway while bypassing others and that each virus uses a unique set of initiation components. A common feature is that all of the tested viruses require the LC3 protein but bypass its canonical cellular lipidation process, using other means to recruit LC3 to virally induced membranes. This study highlights the importance of assessing the full autophagy pathway when seeking to understand how pathogens manipulate this pathway for purposes of genome replication and spread, as well as potential common drug targets.

Diverse viruses utilize different components of the autophagy pathway
To provide an in-depth characterization of which autophagy components are utilized by the three different viruses, we targeted either one or two genes from distinct complexes required for canonical autophagy (Fig 1A). Sequencing confirmed the genomic ablation of the targeted genes in ΔFIP200, ΔATG9, ΔBECN1, ΔVPS34, ΔLC3, and ΔATG5 (S1A Fig). For GC-rich ULK1, specific antibodies confirmed a functional loss of protein, as they did for the targeted gene products in all the CRISPR-Cas9 KO cell lines (S1B Fig). To test whether canonical autophagy was inactivated in these KO cell lines, the accumulation of p62/SQSTM (p62), which is degraded in wild-type (WT) cells via constitutive autophagy, was monitored ( Fig 1B). Accumulated p62 was present at low abundance in the parental HeLa cells (WT) and increased in cells treated with chloroquine (CQ) to prevent the formation of autolysosomes and thus the degradation of autophagosomal contents. All KO cell lines showed enhanced constitutive p62 accumulation, whether or not they were treated with CQ, indicating that the canonical basal autophagy pathway was successfully inhibited (Fig 1B). Additionally, to confirm that the canonical, starvation-induced autophagy pathway was obstructed in all of the KO cell lines, we quantified the presence of endogenous LC3 puncta upon starvation. WT cells displayed significantly more puncta upon starvation than any of the KO lines (S1C Fig).
We then took advantage of these generated KO cell lines to determine which components of the autophagy pathway affect viral growth for a picornavirus (PV) and a flavivirus (DENV). PV is a nonenveloped virus that matures intracellularly. The intracellular amplification of PV in WT and KO cell lines was monitored in a single infectious cycle, which comprises cell entry, translation, RNA replication, and packaging but not cell-to-cell spread. As shown in Fig 1C, cells ablated for ULK1, FIP200, ATG9, and LC3B had significant viral growth defects at this time point. However, viral titers in cells ablated for BECN1, VPS34, and ATG5 were indistinguishable from WT. More detailed time courses of a single PV infectious cycle confirmed the reduction in virus yield in ΔULK1, ΔFIP200, ΔATG9, and ΔLC3 cells (Fig 1C). The fact that ATG5 is not required is particularly curious, given that viral growth was drastically reduced in the ΔLC3B line ( Fig 1C). These observations argue that the virus still utilizes the LC3 protein, though perhaps not its lipidated form. We confirmed that adding back autophagy components to our KO cell lines restored viral titers to WT levels for ULK1 and FIP200 (S2A Fig). Additionally, adding back a kinase dead ULK1 mutant (K46I) reduced viral yield in WT cells, consistent with its previously reported dominant negative effects on the autophagy pathway [51]. ULK1-K46I also failed to rescue the viral defect in ULK1 KO cells (S2A Fig).
DENV is an enveloped virus that matures upon release from cells after budding into the ER and traversing the canonical downstream secretion pathway. We monitored a single cycle of DENV growth by quantifying the release of infectious virus into the extracellular medium. All CRISPR-Cas9 KO cells were infected with DENV for a single 24-hour replication cycle at a low multiplicity of infection (MOI) (Fig 1D). We found that DENV does not require the same initiation complex as PV, as viral titers were unaffected in ΔULK1 cells. We additionally tested whether a double knock down of ULK1 and ULK2 affected DENV growth, given that ULK1 and ULK2 can be redundant in some cell types [9,52]. We saw no defect in DENV titers when small interfering RNAs (siRNAs) against ULK2 were used in the ΔULK1 cell line or when ULK1 and ULK2 siRNAs were used in the ΔFIP200 cell line (S2B Fig), suggesting that DENV does not require this complex. Similar to PV, DENV titers were reduced in ΔATG9 cells, perhaps highlighting a common need for lipid scavenging in HeLa cells. DENV utilized VPS34 but not BECN1, both part of the same complex, arguing that DENV might specifically require VPS34 but not the rest of the complex. Additionally, adding back VPS34 to the ΔVPS34 cell line restored DENV titers to WT levels (S2C Fig). Downstream in the autophagy pathway, DENV utilized LC3 but not ATG5, similarly to PV.
DENV and PV infection in cell culture benefit from the induction of the canonical autophagy pathway when treated with autophagy inducers such as rapamycin [46] or loperamide [39]. Recently, we showed that induction of the autophagy pathway by starvation led to an enhancement of PV amplification that was dependent on the canonical pathway, as ΔATG5 and ΔBECN1 cell lines failed to show enhanced viral infection upon starvation [53]. To test whether in vivo PV infection is also enhanced by starvation as well as whether this exacerbation depends on the canonical autophagy pathway, we bred C57BL/6 mice that were homozygous for a transgene expressing the poliovirus receptor (PVR) [54] as well as the Atg5 gene flanked by Lox sites, the targets of Cre recombinase. These mice were bred either to express or not express Cre recombinase under the control of a tamoxifen-inducible promoter [55]. The conditional deletion of floxed Atg5 upon tamoxifen treatment was confirmed by quantitative PCR (qPCR) of a DNA junction expected to be generated during the Cre-mediated excision (S2D Fig). The decrease of Atg5 function in the Cre-expressing mice was confirmed at the protein level by the large increase in abundance of LC3-I, which normally would be converted to LC3-II by Atg5 (S2E Fig). Although a small amount of residual LC3-II remained in the Creexpressing mice, this is typical of the nonabsolute nature of "floxed" gene removal by Cre recombinase [55]. Mice were fasted for 48 hours before intramuscular infection with PV. Fasted mice showed increased PV titers in mice with an intact Atg5 gene ( Fig 1E). However, in the absence of Atg5, the enhanced viral titer seen with fasting was not observed (Fig 1E). Additionally, we treated mice with loperamide to induce the autophagy pathway by a different mechanism [56,57] and found a similar response of enhanced viral titers in WT mice treated with the drug but no difference between untreated and drug-treated mice lacking Atg5 ( Fig  1F). These data argue that the exacerbation of viral infection by starvation or loperamideinduced autophagy requires ATG5 and most likely the canonical pathway, even though poliovirus growth per se uses only portions of this pathway. We speculate that preinduction of the canonical autophagy pathway increases the concentrations of proteins or lipids that facilitate subsequent viral infection.

PV uses components of the autophagy pathway for early replicative events, while DENV and ZIKV use it for post-RNA replication processes
Previous experiments have argued that in the absence of cellular autophagy or its components, intracellular RNA amplification as well as later steps are defective in PV-infected cells [36,39]. On the other hand, DENV-infected cells lacking an intact autophagy pathway synthesize viral RNA, but virion maturation is defective [46]. To distinguish between amplification in the first round of infection and subsequent spread, intracellular RNA accumulation was monitored by reverse transcription quantitative PCR (RT-qPCR) at early and late time points (Fig 2A). Even in the first cycle of PV infection in the KO cell lines, the decline in RNA accumulation was significant in ΔULK1, ΔFIP200, ΔATG9, and ΔLC3 cells, and this effect continued upon subsequent rounds (Fig 2B). To test whether viral entry was affected in the KO cell lines, we performed an entry assay. Virus was allowed to infect cells for 30 minutes, after which virions that had not entered were stripped from the cell surface with an acid wash. RT-qPCR analysis of PV RNA that successfully entered the cell showed no difference between WT and KO cell lines To determine whether autophagy components influence DENV translation and RNA synthesis, viral packaging and spread, or both, viral RNA synthesis was monitored after single (24 hpi [hours post infection]) and multiple (48 hpi) rounds of infection. In the absence of ATG9, VPS34, or LC3, no reduction in DENV RNA was observed in the first cycle (Fig 2B), even though the production of infectious virus was significantly reduced at this time point (Fig 1D). This is consistent with previous observations that in the presence of autophagy inhibitor spautin-1, packaged virions were defective but RNA accumulation was not [46]. Upon additional rounds of infection, the accumulation of DENV viral RNA was greatly reduced (Fig 2B), consistent with the poor infectivity of virions from the first round. In accordance with these data, we observed similar amounts of viral protein accumulation at early time points in WT and ΔATG9 cells but reduced accumulation at later time points (S3E Fig). We conclude that in the absence of a noncanonical autophagy pathway, initial DENV translation and RNA synthesis are not affected, but infection of subsequent rounds is greatly reduced as a result of defective assembly, maturation, or egress.
ZIKV is a flavivirus that is similar to DENV but infects a larger number of cell types in humans [47]. ZIKV RNA was monitored during infection, and like DENV, viral RNA abundance remained similar to WT in the first infectious cycle (24 hpi). However, after additional infectious cycles, both RNA ( Fig 2B) and virus production ( Fig 2C) were greatly reduced. These data argue that for ZIKV, like DENV, autophagy components are only necessary for a postreplication process such as packaging or particle maturation. However, unlike either PV or DENV, the accumulation of infectious ZIKV was reduced in the entire CRISPR-Cas9 KO panel except the ΔATG5 line ( Fig 2C). This argues that ZIKV utilizes both upstream initiation complexes ULK1/FIP200 and BECN1/VPS34, perhaps feeding into many aspects of the canonical autophagy pathway. Nonetheless, the fact that all three viruses bypass the need for ATG5-mediated lipidation of LC3B highlights a shared strategy.
To monitor the accumulation of viral RNAs on a single-cell basis, we measured the number of infected cells during multiple cycles of infection by flow cytometry. Virus-specific probes were designed to hybridize to each positive-strand viral genome. These probes can subsequently be bound by branched DNA structures, with each DNA "tree" containing thousands of fluorophores as a means to amplify the signal, which can allow the detection of single RNA molecules [58,59]. The low background for these probes from uninfected cells was validated in WT HeLa cells by confocal microscopy (Fig 2D). These probes were then used to determine the percentage of cells positive for viral RNA by flow cytometry, using a low MOI (0.1 plaque forming units [PFU]/cell) and multiple cycles of infection. In the cell lines that showed significant defects in growth for each virus, fewer cells were positive for PV, DENV, or ZIKV viral RNA ( Fig 2D).

Differential utilization of upstream autophagy components by RNA viruses
Autophagy can be induced by two different arms of upstream signaling: mTOR inactivation, leading to dephosphorylation of ULK1 and thus its activation, or AMPK activation, leading to a distinct phosphorylation of ULK1 in the absence of mTOR repression [8,60,61]. mTOR typically responds to nutrient signals while AMPK responds to the energy status of the cell. Other viruses have been shown previously to activate the autophagy pathway via AMPK activation [62]. The differential utilization of upstream autophagy components by PV, DENV, and ZIKV could result from the activation of different upstream signaling cascades. To test this possibility, we looked at the phosphorylation status of S6K, which is phosphorylated when mTOR is active and autophagy is repressed [63]. We saw that while rapamycin treatment led to the loss of S6K phosphorylation, PV, DENV, and ZIKV infection did not alter its phosphorylation status ( Fig 3A). This suggests that the activation of autophagy during viral infection is independent of mTOR inactivation and thus, by definition, noncanonical. We next looked at AMPK phosphorylation status and discovered that PV and ZIKV led to AMPK phosphorylation and thus activation but not DENV ( Fig 3B). Since AMPK is known to directly activate ULK1, these data support the idea that PV and ZIKV utilize the ULK1 complex in activating the autophagy pathway, while DENV bypasses both mTOR and AMPK, perhaps avoiding ULK1 activation entirely.
To determine whether ULK1 and FIP200 function in the same way in stimulating PV growth as they do in canonical autophagy, we used a small molecule inhibitor of ULK1 (MBL56), which specifically blocks its kinase activity [64]. In WT HeLa cells, addition of MBL56 decreased viral titers after 6 hours of infection in a dose-dependent manner ( Fig 3C). However, MBL56 had no additional effect on viral titers in infected ΔFIP200, arguing that ULK1 and FIP200 function together just as they do in canonical autophagy.
DENV bypasses the upstream ULK1 complex and instead requires VPS34. To determine if DENV requires a catalytically functional VPS34 protein, we used VPS34-specific inhibitor SAR405 [65]. Addition of SAR405 led to reduced DENV titers at 24 hpi in a dose-dependent manner ( Fig 3D). We see reduced ZIKV titers with addition of SAR405 as well, showing that ZIKV also requires a catalytically functional VPS34 ( Fig 3D). Interestingly, we did not see any reduction in PV titers with addition of SAR405 (Fig 3D), consistent with the previous conclusion that PV does not require VPS34 activity. Therefore, PV and ZIKV induce the autophagy pathway upstream of AMPK, while DENV uses an alternate route to induce autophagy pathway, perhaps through the direct activation of VPS34.

Virally induced membrane rearrangements are altered in cells that lack individual autophagy components
To test whether the extensive membrane rearrangements that occur in cells infected with PV and DENV are altered in cell lines lacking particular autophagy components that affected viral growth, we studied the ultrastructure of infected WT and KO cells. WT, ΔULK1, and ΔATG9 cells were infected with PV at a high MOI (10 PFU/cell) for 6 hours, followed by high-pressure freezing and freeze substitution to preserve membrane structures. Images were visualized by electron microscopy (EM) (Fig 4A). PV-infected WT cells showed the characteristic DMVs shown previously to be induced either by infection or coexpression of PV proteins 2BC and 3A [34,35]. These DMV structures were not present when cells were treated with guanidine, a potent inhibitor of viral replication (Fig 4A and 4B). In ΔULK1 and ΔATG9 cells, qualitatively different membrane rearrangements were observed. Infection of ΔULK1 cells led to an increase in single-membraned vesicles, while infection of ΔATG9 cells led to an increased number of electron dense single-membraned vesicles, both of which diminished when cells were treated with guanidine ( Fig 4A and 4B). Single-membraned vesicles can be the precursors of DMVs in canonical autophagy; this may also be the case in PV-infected cells [66]. Therefore, the single-membraned vesicles observed in ΔULK1 and ΔATG9 cells could represent precursors or off-pathway structures that accumulate during arrest of the functions of the ULK1 complex. The presence of such precursors could also be increased by the lower viral protein abundance in these KO cell lines (S3B Fig). No alterations in membrane structures were seen in any cells in the presence of guanidine, arguing that translation from input RNA is not sufficient to induce these changes.
DENV growth does not use the canonical ULK1 complex but instead requires both VPS34 and ATG9 (Fig 1D). To identify any alteration in membrane structure under these conditions, we visualized DENV-infected WT, ΔVPS34, and ΔATG9 cells by EM. Typically, during a DENV infection, DMVs and large convoluted membrane structures thought to be extensions of the ER are formed by 24 hpi (Fig 4C) [67]. These structures were notably absent in ΔVPS34 and ΔATG9 cells (Fig 4C and 4D), even though comparable amounts of viral RNA and protein were observed in these cell lines at this time point (Fig 2B and S3E Fig), arguing that VPS34 and ATG9 play crucial roles in the architecture of these membrane structures.

LC3 is recruited to membranes independent of lipidation during viral infection
ATG5 is often considered the hallmark of canonical autophagy due to its essential function in the lipidation of LC3. Considering that we observed no defects in the growth of PV, DENV, or ZIKV in ΔATG5 cells but significant dependence on LC3, we were curious about the mechanism of LC3 utilization during these viral infections. To observe the membrane-associated sequestration of LC3 during canonical autophagy and viral infection, we monitored the relocalization of a plasmid expressing a green fluorescent protein (GFP)-LC3 fusion protein [20,68]. Under normal conditions, GFP-LC3 is dispersed throughout the cytoplasm, as shown in the top panel of Fig 5A. Upon induction of the canonical autophagy pathway, GFP-LC3 forms distinct puncta that represent membrane-bound, lipidated LC3. Under starvation conditions, GFP-LC3 puncta form in WT cells but not in ΔATG5 cells. Infection with PV, DENV, or ZIKV also led to the formation of GFP-LC3 puncta in both the presence and absence of ATG5 (Fig 5A and 5C).
To test further whether the induction of LC3 puncta by PV, DENV, and ZIKV was truly independent of LC3 lipidation, we expressed a version of GFP-LC3 that contained a G120A mutation, which lacks the essential glycine residue needed for the attachment of PE [69]. As expected, starvation did not induce the formation of GFP-LC3-G120A puncta in either the presence or absence of ATG5 (Fig 5B, upper panels). The fact that very few GFP-LC3-G120A puncta were observed in WT cells might indicate that the virus preferentially induces puncta formation with lipidated LC3, which is still present as untagged endogenous LC3. However, GFP-LC3-G120A puncta formation was observed in ΔATG5 cells upon infection with PV, DENV, or ZIKV (Fig 5B and 5C). These results suggest that all three viruses can either recruit unlipidated LC3 to membranes or can lead to lipidation of LC3 independently of ATG5.
To determine whether or not LC3 was lipidated during infection, we enriched for membrane-associated proteins and determined the lipidation status of LC3 by immunoblot. Cells were treated with rapamycin, CQ, or infected with PV. As expected, no LC3-II was observed in ΔATG5 cells, and treatment with CQ increased LC3-II abundance in WT but not ΔATG5 cells when whole-cell lysates were tested ( Fig 5D). Proteins most tightly associated with membranes can be enriched by treatment of extracts with saponin [70]. LC3-I was retained in the membrane fraction in PV-infected cells in both WT and ΔATG5 cells (Fig 5D). Thus, PV is capable of recruiting LC3-I to membranes in the absence of lipidation and the canonical lipidation machinery. In agreement with the above data, we observed by EM that PV infection still leads to the formation of DMVs in ΔATG5 cells (Fig 5E and 5F) and that these structures are dependent on active viral replication since they are absent in cells treated with guanidine. However, no DMVs were seen upon PV infection of ΔLC3B cells, but an increase in singlemembraned vesicles was observed. Therefore, PV is capable of forming DMVs without lipidated LC3 but still appears to require the LC3 protein itself. Interestingly, PV infection of ΔATG5 cells resulted in the formation of many large vesicles that resemble multivesicular bodies (MVBs) (Fig 5E and 5F). These might represent an altered curvature capacity of unlipidated LC3 alone [22], as compared to the increased curvature in PV-infected WT cells, which contain both LC3-I and LC3-II ( Fig 5D).

Viral proteins bind LC3
To determine how PV leads to recruitment of unlipidated LC3 to membranes, we used an anti-GFP antibody to capture GFP-LC3 and tightly associated, detergent-resistant proteins during infection. Mass spectrometry (MS) was then used to identify the co-immunoprecipitated (IP) proteins. To confirm that our GFP-LC3 IP was specific, we compared several cellular proteins known to interact with both LC3-I and LC3-II or only with LC3-II (Fig 6A). The LC3 lipidation process starts with a C-terminal cleavage by the proteinase ATG4B and covalent attachment first to ATG7 and then to ATG3 [71][72][73]. The ATG5/12/16L complex then acts as an E3 ligase, transferring the covalent attachment of LC3 from ATG3 to PE [74]. LC3 binds to p62 whether it is lipidated or not [75,76], and, indeed, we observe equal binding to p62 in WT and ΔATG5 cells and with GFP-LC3 or GFP-LC3-G120. However, autophagy-associated protein FYCO1 only binds to LC3-II [77], as there were fewer peptide reads in ΔATG5 cells or when GFP-LC3-G120A was used. We saw little to no binding of ATG4B, which binds the terminal glycine of LC3 [72,73], when GFP-LC3-G120A was used. Interesting, ATG4B-LC3 binding was significantly reduced in ΔATG5 cells, perhaps because LC3-II is also one of its substrates [20,78]. ATG7 and ATG3, which participate in the covalent conjugation of LC3 upstream of the function of ATG5, bound to GFP-LC3 in the presence or absence of ATG5. However, the G120A substitution cannot initiate the cascade of covalent attachments and is not recognized by either ATG7 or ATG3 [79]. Thus, the recovery of peptides from all five of these LC3-interacting proteins exhibits differential dependence on the lipidation status of LC3 in accordance with their biological roles.
To identify viral proteins that directly bind to LC3 as well as those that are found on the same membranous surface, we prepared samples that did or did not contain detergent to disrupt membranes. Plasmids that expressed either GFP-LC3 or GFP-LC3-G120A were transfected into WT or ΔATG5 cells and infected with PV. The viral peptides identified by IP-MS were mapped across the viral genome to visualize the viral proteins pulled down with LC3 ( Fig  6B). In the absence of detergent, peptides from almost every viral protein were recovered (S4A Fig). However, in the presence of detergent, most peptides pulled down from WT cells were derived from capsid proteins and 2C (Fig 6B). When we analyzed PV peptide reads as percentage of total reads, capsid proteins formed a large percentage of the total peptide reads in all conditions, suggesting that these viral proteins bind LC3 or LC3-associated membranes ( Fig  6C). The ratio of reads in the presence of detergent to those in its absence showed that only peptides in 2C in WT cells remained associated with LC3 in the presence of detergent ( Fig  6D). We hypothesize that LC3 binding to viral protein 2C, a known hydrophobic protein, is stabilized by LC3 lipidation.
To determine which form of protein 2C immunoprecipitated with GFP-LC3, we visualized the collected proteins by SDS PAGE and immunoblotted with anti-2C antibody (Fig 6E). Interestingly, precursor 2BC but not 2C was pulled down, arguing that the peptides identified by MS actually derived from 2BC. We additionally found 3A pulled down with LC3, an interesting observation because 2BC and 3A together can induce membrane rearrangements similar to those induced by viral infection [35], and lead to LC3-II accumulation outside of viral infection [80]. We did not find 2B or 3A peptides pulled down in our MS data, perhaps because these are small proteins and thus difficult to detect by MS.
A consensus LC3-interacting region (LIR) has been identified with a common WxxL motif, and contacts between this motif and LC3-binding protein p62 have been visualized by X-ray crystallography [81,82]. Only two potential LIR motifs were identified in the PV genome: WWKL in viral capsid VP2 and WQWL in viroporin 2B, marked in the genome schematic ( Fig 6B). Their high conservation in picornaviruses is shown in S3B  Viruses and autophagy GFP-LC3 IP-MS was also performed on DENV-infected bovine hamster kidney (BHK) cells with WT GFP-LC3 and GFP-LC3-G120A. To confirm that the MS data were specific to LC3 lipidation conditions, we saw both GFP-LC3 and GFP-LC3-G120A bound to p62 but not to ATG3 or ATG7 when LC3-G120A was used (S5A Fig). All DENV proteins were pulled down with WT GFP-LC3 and GFP-LC3-G120A (S5B Fig). This is not surprising, given that DENV proteins are known to directly interact with each other [67,83,84], and suggests that one or more DENV proteins likely bind to LC3. Possible LIR motifs found in NS1, NS2A, and NS5 could potentially be sites of interaction with LC3.

Discussion
By comparing how different RNA viruses manipulate autophagy factors, we show here that there are numerous ways by which a virus can hijack this complex cellular pathway. Each virus, even closely related viruses like DENV and ZIKV, uses slightly different components, although they might ultimately use the subsequent membranes and degradative processes for similar purposes. Each virus requires a different subset of initiation components for efficient viral growth, although all three viruses utilize lipid scavenger ATG9 and recruit LC3 directly to membranes, bypassing the need for ATG5-mediated lipidation (Fig 7). We believe that this comparative study addresses some of the inconsistencies seen with viral induction of the autophagy pathway.
The fact that the ULK1/FIP200 and BECN1/VPS34 complexes can be bypassed by individual viruses is consistent with the idea that noncanonical autophagy pathways can be induced in response to different stress stimuli [85]. ULK1-independent autophagy can be induced by glucose deprivation or excessive ammonia [23,86]. BECN1-independent autophagy can be induced by proapoptotic compounds such as Z18 [24,87]. There is also evidence for a noncanonical autophagy pathway that is VPS34-dependent and BECN1-independent upon induction by arsenic trioxide [88], providing precedent for the DENV-induced noncanonical pathway. These data highlight that functional cellular autophagy is not entirely reliant on every component of the canonical pathway. Perhaps many viruses take advantage of these alternate and possibly redundant pathways to ensure a productive infection [89].
One of the striking findings described here is the dependence of PV, DENV, and ZIKV infection on LC3 but not on its lipidation. It has been shown previously that unlipidated LC3-I is found on viral-associated vesicles during infection with murine hepatitis virus, equine arteritis virus, and Japanese encephalitis virus [90][91][92]. These data, in addition to LC3-dependent and ATG7-independent viral amplification, led these authors to conclude that the virally induced membranes were not derived from the autophagy pathway. Instead, colocalization of ER degradation enhanced by alpha-mannosidase (EDEM1) and LC3 suggested that the virally induced vesicles were ER-associated protein degradation (ERAD) associated [90]. We argue that direct recruitment of LC3 to membranes by viruses can perform autophagosome-like functions. Colocalization with EDEM1 could be a result of the known degradation of EDEM1-containing vesicles by an autophagy-like process [93]. We specifically found unlipidated LC3 associated with membranes during viral infection as well as direct binding of LC3 to PV proteins by IP-MS. A few cases of cellular LC3 lipidation-independent autophagy have been described [25,94]. In these cases, when either ATG5 or ATG7 is ablated, the autophagy pathway can still proceed, leading to the formation of autophagosomes in response to unusual stimuli, such as the drug etoposide [25]. LC3 is not lipidated under these conditions. It remains unknown how these autophagosomes fuse without lipidated LC3. LC3-I has been found associated with membranes during cellular secretory processes, although how LC3-I is recruited to membranes was not tested [93]. It is possible that in these lipidation-independent cases protein-mediated recruitment of LC3-I occurs, as we observed during viral infection.
Recently, the influenza M2 protein was found to contain a functional LIR motif. This LIR motif was shown to be important for viral-mediated relocalization of LC3 to the plasma membrane, presumably to facilitate viral budding [95]. This represents a striking example of autophagosome relocalization by a virus. We have identified potentially functional viral LIR motifs in the capsid protein VP2 and 2BC protein of PV, as well as in the polyproteins of DENV and ZIKV. One hypothesis is that PV proteins bind LC3 as a means to localize PV-replication complexes and capsids specifically to autophagosome-like membranes for assembly purposes. 2BC is localized to replication complexes and is a membrane-associated protein [96], so an interaction with LC3 at the outer surface of the double membrane is possible. The binding of capsid proteins to LC3 was unanticipated but of potential functional interest. Immature capsid precursors are known to be specifically associated with membranes [97] and intact virions can be released nonlytically by a process similar to secretory autophagy [39,40,98]. We propose that LC3 association of 2BC and immature capsids facilitate RNA replication, RNA packaging, and intercellular spread.
Autophagosome-dependent unconventional secretion was first demonstrated in yeast [99][100][101]. In human cells, autophagy-dependent secretion has been demonstrated for interleukin 1β [102] and increasing numbers of other proteins and complexes [26,103]. This suggests that autophagosomes can be rerouted to the plasma membrane and the contents released, providing precedent for virally induced double-membraned structures exiting a cell nonlytically. Infectious picornaviruses such as PV, Coxsackievirus B, and Hepatitis A virus have been found to be released in membrane-bound vesicles, each containing many viral particles [40,41,104]. For both unconventional secretion and the assembly of intracellular membranous components, the cellular autophagy pathway, with its unusual topologies and induction of membrane curvatures, is a rich source of host components for many viral infections.
We observed in mice that fasting could enhance viral growth of PV. We also show that the increase in PV titers is dependent upon the presence of Atg5, linking the starvation phenotype to its role in inducing the autophagy pathway. The data presented here show that PV benefits from the induction of the canonical autophagy pathway but do not require all of the components when virally inducing the pathway. We hypothesize that pre-existing membrane structures from starvation-induced autophagy can be utilized by viruses for growth and maturation, leading to quicker exit from cells and overall more virus production. Malnutrition has been linked to increased death from infectious diseases in children [105], and, although exacerbation of disease by malnutrition is often attributed to confounding effects, another outcome of malnutrition is induction of cellular autophagy. The implications of these data are broad: if viral infection is exacerbated by induction of the autophagy pathway, then starvation and many over-the-counter drugs that induce autophagy could lead to greater disease severity. It also implies that simple supplementation of nutrients to block induction of the autophagy pathway may help abrogate increased viral amplification.

Cells and viruses
Human HeLa cells and BHK21 cells were cultured in Dulbucco's modified eagle's medium (Gibco) supplemented with 10% bovine serum and 1% penicillin/streptomycin. PV type 1 Mahoney was grown from infectious cDNA, as previously described [106]. DENV type 2 (16681) was propagated from an infectious cDNA clone (pD2/IC) in C6/36 mosquito cells. The ZIKV strain used was PRVABC59 (Puerto Rico). The virus was purchased from BEI Resources and propagated on C6/36 mosquito cells.

CRISPR-Cas9 KO cell line generation
CRISPR-Cas9 plasmid px458 (Addgene #48138) was cut with BbsI and ligated to annealed sgRNAs (see S1 Table), as described previously [107]. HeLa cells were transfected with resulting plasmids (Lipofectamine 3000, Invitrogen). Forty-eight hours later, cells were single-cell sorted for GFP+ cells using a BD Aria II sorter (Stanford FACS facility) in a 96-well format. Cells that propagated were tested for gene disruption by harvesting genomic DNA, PCR amplifying the region of interest, and sequencing to look for the presence of frameshift mutations. Clones were also tested for reduced protein expression by western blot. For the ATG9, BECN1, VPS34, and LC3B genes, multiple KO lines were obtained and tested, although subsequent experiments were carried out using one representative line.

Confocal microscopy
Cells were transfected with 50 ng/well of GFP-LC3 in a 24-well plate. After 24 hours, cells were seeded on coverslips and infected or treated the following day. Cells were fixed in 4% paraformaldehyde, stained with DAPI (Invitrogen), and mounted on slides with PermaFluor (Thermo Fisher Scientific). Slides were imaged using a Leica SP8 confocal microscope (Stanford imaging facility). GFP-LC3 puncta were counted per cell.

Flow cytometry and viral-specific probes
Infected cells from a 24-well plate were fixed and stained using the PrimeFlow RNA assay according to the manufacturer's protocol (Invitrogen). Viral RNA-specific probes were used for the positive strand of PV (VF1-10252), DENV2 (VF1-15158), or ZIKV (VF1-20236). Cells were run on a flow cytometer (Scanford, Stanford FACS facility) and analyzed using FlowJo software (v. 10.4).

RT-qPCR
Total RNA was isolated from cells using the RNeasy mini kit (Qiagen). RT-qPCR was performed on an Applied Biosystems 7300 machine using the QuantiTect SYBR Green RT-qPCR Kit (Qiagen). Primers specific to PV, DENV, or ZIKV were used (see S2 Table).

Luciferase assay
Cells were plated in a 96-well plate and transfected with 8 ng of PV replicon RNA (kind gift from Raul Andino). A subset of cells were treated with 2 mM guanidine (Sigma) prior to transfection to inhibit viral replication. Cell lysates were collected at 0, 2, 4, and 6 hours post transfection, and firefly luciferase activity was measured using the Firefly Luciferase Assay System (1000) in a 96-well format with a BioTek Neo2 luciferase reader.

Tranfections and Lentiviral transductions
DNA transfections were done using a Lipofectamine 3000 kit (Invitrogen). The ULK1 plasmid was acquired from Addgene (#31963). The K46I mutation was made by site-directed mutagenesis (QuickChange Lightening kit, Agilent). Lentiviral vectors were cloned using Gibson cloning (NEB), inserting the gene of interest into the pLenti-CMV-Puro-DEST vector, as described previously [108]. 293T cells were transfected with the pLenti vector containing the addback gene, Gag-Pol, VSV-G, and p-Advent. Supernatants were collected 24 hours later and treated with protamine sulphate. Supernatants were added to KO cell lines and allowed to infect for 24-48 hours. Successful transductions were selected by 2 μg/ml puromycin for up to 1 week. Resulting cell lines were tested for rescue of protein levels by western blot. For siRNA knock-down of ULK1 and ULK2, 4 pooled siRNAs targeting human ULK1 or ULK2 (siGEN-OME SMARTpool, Dharmacon #005049 and #005396) were transfected into HeLa cells using Lipofectamine 3000 (Life Technologies). Knockdown was tested by qPCR at 72 hours post transfection.

Electron microscopy
WT and KO cell lines were plated in 10-cm plates and infected with PV (6 hours) or DENV (24 hours) at an MOI of 10 PFU/cell. A subset of samples were infected and treated with 2 mM guanidine (Sigma) simultaneously. Cells were collected and fixed with 4% PFA for 10 minutes.
Cells were washed twice with DMEM and resuspended in 50 ul of 20% BSA in DMEM. Samples were subjected to high-pressure freezing using a Leica EmPACT High Pressure Freezer and freeze substitution to replace water with acetone (Stanford Imaging Facility). Samples were sectioned and imaged on a transmission electron microscope (JEOL JEM1400, Stanford Imaging Facility). Per condition, 10-20 cells were analyzed for cellular structures on blinded images.

Immunoprecipitations and mass spectrometry
GFP immunoprecipitations were performed with GFP-Trap magnetic beads (Chromotek). Cells were harvested with lysis buffer (10 mM Tris pH 7.5, 150 mM NaCl, 0.5 mM EDTA, and 0.5% NP-40) containing EDTA-free protease inhibitor (Roche). Magnetic beads were added and allowed to bind for 2 hours at 4˚C with rotation. Beads were washed 3 times with wash buffer (10mM Tris pH 7.5, 150 mM NaCl, and 0.5 mM EDTA) and eluted with 2 M glycine pH 2.5, followed by neutralization with 1 M Tris pH 10.4, for mass spectrometry. Samples analyzed by immunoblot were eluted by boiling in 2X SDS sample buffer (120 mM Tris pH 6.8, 20% glycerol, 4% SDS, 0.04% bromophenol blue, and 10% BME) and run on an SDS PAGE gel. Mass spectrometry was performed and analyzed by the Stanford University Mass Spectrometry facility (SUMS).

Mouse infections
Embryos of conditional Atg5 flox/+ mice B6.129S (Hara et al., 2006) were obtained from Dr. Noboru Mizushima through the RIKEN BioResource Center (Ibaraki, Japan; stock RBRC02975) and cryorecovered at Stanford University. Atg5 flox/+ mice were further backcrossed to C57BL/6J to N10. B6.Cg-Ndor1Tg(UBC-cre/ERT2)1Ejb/J, a lentitransgenic line harboring a single-copy integrant of a fusion gene consisting of Cre recombinase and a mutant human estrogen receptor (ERT2) under control of the human ubiquitin C promotor (UBC), was obtained from the Jackson Laboratory (Bar Harbor, Maine; stock 008085) at N6 and further backcrossed to C57BL/6J to N10. C57BL/6 mice expressing the human PVR, Atg5 fl/fl , and Cre +/were treated with 75 mg/kg of tamoxifen for 5 days, followed by 2-4 weeks of recovery to induce the expression of Cre and subsequent excision of the floxed Atg5 exon 2. Ten to twelve-week-old mice were fasted for 48 hours prior to infection or given 25 mg/kg loperamide (Sigma) in saline by i.p. injection every 12 hours for the duration of the infection. Mice were infected intramuscularly with PV (1 x 10 7 PFU/50 ul in PBS). Muscle tissue was harvested 4 dpi and processed using the Bullet Blender BBX24 (Next Advance). DNA and protein lysate were also extracted from tissues, using the DNeasy Blood and Tissue kit (Qiagen) and RSB buffer, respectively.

Ethics statement
The mouse experiments in this study were approved by the Administrative Panel on Laboratory Animal Care at Stanford University (approval number 9296). This committee is accredited by the Association for the Assessment and Accreditation of Laboratory Animal Care (AAALAC International). Calf muscle tissue was harvested, and DNA was extracted. qPCR was done for the indicated regions of the Atg5 gene. (E) The same mouse tissue as above was also used to extract protein lysates, run on an SDS PAGE gel and blotted for LC3. All data are represented as mean +/− SEM. � Indicates significant P value of <0.05, �� P value < 0.01, ��� P value < 0.001, ���� P value