Abstract
Phosphatidylinositol 3-phosphate (PI3P) is a signaling phospholipid that play a key role in endomembrane trafficking, specifically autophagy and endosomal trafficking. However, the mechanisms underlying the contribution of PI3P downstream effectors to plant autophagy remain unknown. Known PI3P effectors for autophagy in Arabidopsis thaliana include ATG18A (Autophagy-related 18A) and FYVE2 (Fab1p, YOTB, Vac1p, and EEA1 2), which are implicated in autophagosome biogenesis. Here, we report that FYVE3, a paralog of plant-specific FYVE2, plays a role in FYVE2-dependent autophagy. Using yeast two-hybrid and bimolecular fluorescence complementation assays, we determined that the FYVE3 protein was associated with autophagic machinery containing ATG18A and FYVE2, by interacting with ATG8 isoforms. The FYVE3 protein was transported to the vacuole, and the vacuolar delivery of FYVE3 relies on PI3P biosynthesis and the canonical autophagic machinery. Whereas the fyve3 mutation alone barely affects autophagic flux, it suppresses defective autophagy in fyve2 mutants. Based on the molecular genetics and cell biological data, we propose that FYVE3 specifically regulates FYVE2-dependent autophagy.
Keywords: ATG8, autophagosome, autophagy, phosphatidylinositol 3-phosphate, phosphoinositide, vacuolar trafficking
Introduction
Endomembrane trafficking comprises dynamic membrane flow where various membrane lipids and cargo proteins are sorted and transported to their appropriate destinations. Newly synthesized proteins are targeted to various endomembrane compartments via multiple routes, including secretory and vacuolar trafficking. Proteins endocytosed from the cell surface are sorted at the endosomal compartments and either recycled to the plasma membrane or targeted to the vacuole (or lysosome in animals) for degradation. Cytoplasmic proteins may be sequestered by the autophagic membrane and delivered to the vacuole/lysosome for degradation.
The best-characterized type of autophagy is macroautophagy (hereafter referred to as autophagy), in which a membrane sac called the phagophore is expanded and sealed to form the autophagosome sequestering a small portion of the cytoplasm (Marshall and Vierstra, 2018). The autophagosome fuses with the vacuole, releasing its cargo into the vacuolar lumen as an autophagic body. Finally, the autophagic body is rapidly degraded by acid hydrolases in the vacuole. The formation of phagophores and autophagosomes is facilitated by autophagy-related (ATG) proteins, which are conserved in nearly all eukaryotes. ATG8 is a ubiquitin-fold protein conjugated to the membrane lipid phosphatidylethanolamine on the phagophore. As ATG8 resides on the autophagic membrane even after autophagosome fusion, lipidated ATG8 is a useful marker for phagophores, autophagosomes, and autophagic bodies. ATG5 and ATG7 are among the ATG proteins essential for the ATG8 conjugation and general autophagy, whereas ATG2 and ATG18A form a protein complex that transfers membrane lipids to the phagophore and contributes to phagophore growth (Kim et al., 2020).
Phosphatidylinositol 3-phosphate (PI3P) is a signaling phospholipid that plays regulatory roles in endomembrane trafficking (Schink et al., 2013). In eukaryotic cells, PI3P is enriched at the cytosolic leaflet of endosomal, autophagosomal, and vacuolar (or lysosomal) membranes. In budding yeast and plants, PI3P is produced from phosphatidylinositol by the phosphatidylinositol 3-kinase (PI3K) complex containing the catalytic subunit VPS34 (Vacuolar Protein Sorting 34). There are two distinct PI3K complexes which have 3 common and 1 specific subunits. The PI3K complex I, consisting of VPS34, VPS15, VPS30/ATG6, and ATG14, is responsible for PI3P production at the autophagic membrane, while the PI3K complex II, having VPS38 instead of ATG14, generates PI3P on the endosomal membrane (Schink et al., 2013).
Despite its low abundance in biological membranes, PI3P on endosomal and vacuolar membranes is essential for plant development. Null mutations of vps34 and atg6 in Arabidopsis thaliana cause lethality (Fujiki et al., 2007; Lee et al., 2008), and vps38 mutants show pleiotropic phenotypes resulting from defective endosomal sorting and vacuolar trafficking (Lee et al., 2018; Liu et al., 2018; Yuan et al., 2018). Less obvious developmental phenotypes are observed in Arabidopsis atg14a atg14b double mutants, although they are defective in autophagy (Liu et al., 2020). Generally, plants with compromised functions of vacuoles and endosomes display lethality or severe developmental abnormality (González Solís et al., 2022), whereas autophagy-deficient plants do not. Instead, the autophagy mutants are hypersensitive to different types of abiotic and biotic stress (Qi et al., 2021).
The mechanism underlying the contribution of PI3P to autophagy and endosome trafficking remains poorly understood in plant cells (Chung, 2019). The regulatory roles of PI3P in endomembrane trafficking are primarily mediated by PI3P effectors, which bind to PI3P and, directly or indirectly via protein interactions, serve different downstream functions such as cargo selection, membrane bending, vesicle transport, tethering, scaffolding, and signaling (Schink et al., 2013). In Arabidopsis, PI3P-binding proteins have been shown to localize at endosomal and/or autophagic membranes. For example, SNX1 (Pourcher et al., 2010) and SNX2b (Phan et al., 2008) localize to endosomes, and FYVE4 (Liu et al., 2021) is primarily cytosolic but can be recruited to endosomal membranes. FYVE1/FREE1 (Gao et al., 2014; Gao et al., 2015), ATG18A (Zhuang et al., 2017; Kang et al., 2018; Kim et al., 2022) and FYVE2/CFS1 (Cell death-related endosomal FYVE/SYLF protein 1) (Sutipatanasomboon et al., 2017; Kim et al., 2022; Zhao et al., 2022) localize at both endosomes and autophagosomes. Genetic analyses of mutants defective in the PI3P-binding proteins indicated that these proteins are involved in endosomal sorting (e.g., SNX1, SNX2b, FYVE4), autophagic degradation (e.g., ATG18A, FYVE2/CFS1), or both (e.g., FYVE1/FREE1). However, molecular mechanisms for these putative PI3P effectors are not fully defined, and how PI3P binding impacts their functions is unknown.
PI3P effectors typically bind to PI3P via their FYVE, PX, or PH domains (de Jong and Munnik, 2021). In Arabidopsis, the FYVE domain is shared by FYVE1/FREE1, FYVE2/CFS1, FYVE3/CFS2, and FYVE4. Both FYVE2/CFS1 and FYVE3/CFS2 proteins harbor the PI3P-binding FYVE domain and C-terminal SYLF domain, which assists yeast SYLF-containing proteins to bind actin filaments and phosphoinositides (Robertson et al., 2009; Urbanek et al., 2015; Sutipatanasomboon et al., 2017) ( Figure S1A ). FYVE2/CFS1 facilitates the efficient progression of autophagy and interacts with a variety of proteins such as ATG18A, SAR1B, and VPS23A, which are important for autophagosome biogenesis, COPII vesicle trafficking, and membrane scission, respectively (Sutipatanasomboon et al., 2017; Kim et al., 2022; Zhao et al., 2022). By contrast, no apparent defect in autophagic flux was observed in its paralogous fyve3/cfs2 single mutants (Kim et al., 2022; Zhao et al., 2022). Thus, currently no evidence supports the involvement of FYVE3/CFS2 in autophagy and other membrane trafficking routes.
Here, we aimed to assess genetic and cell biology data to ascertain the role of Arabidopsis FYVE3/CFS2 in autophagy. We showed that FYVE3 was incorporated into the autophagic machinery via interaction with ATG8. An autophagy pathway consisting of core ATG proteins and PI3P effectors was proposed, based on our genetic data indicating a specific involvement of FYVE3 in FYVE2-dependent autophagy.
Materials and methods
Plant materials
The Arabidopsis thaliana T-DNA insertional mutants of atg2-1 (Inoue et al., 2006), atg5-1 (Thompson et al., 2005), atg7-2 (Chung et al., 2010), fyve2-2, and fyve3-1 (Kim et al., 2022) and transgenic plants expressing ProUBQ10:GFP-ATG8A (Kim et al., 2013) and ProVHA-a1:VHAa1-GFP (Dettmer et al., 2006) were previously described, and their progenies were used. ProUBQ10:YFP-ARA7, ProUBQ10:YFP-SYP32 (Geldner et al., 2009), and ProUBQ10:Citrine-2xFYVE (Simon et al., 2014) were obtained from the Arabidopsis Biological Resource Center (ABRC).
Clones and DNA constructs
cDNA and expression clones were obtained from the ABRC with the following identifiers: Arabidopsis ATG6 (G19464), FYVE2 (G18828), FYVE3 (GC105091), SAR1B (G19464), SAR1C (G11541), pIX-Halo-ATG8F (HALO_SFI_23-A01), pIX-Halo-ATG8I (HALO_SFI_71-H06), PSAT4-DEST-n(1-174)EYFP-C1 (CD3-1089), PSAT5-DEST-c(175-END)EYFP-C1(B) (CD3-1097). To make entry clones containing FYVE3 cDNA lacking the stop codon and TRAF1A, FYVE3 and TRAF1A cDNA were amplified using primers listed in Supplementary Table 1 . The FYVE3 and TRAF1A PCR products were inserted into the pENTR/TOPO and pDONR221 vectors (Thermo Fisher Scientific), respectively. cDNA entry clones containing ATG8E, ATG8F, ATG8I, and ATG18A was generated as described previously (Kim et al., 2022).
Preparation of transgenic plants
To obtain binary vectors for transgenic plants expressing ProUBQ10:GFP-FYVE3 or ProUBQ10:mCherry-FYVE3, recombination of the entry clone containing FYVE3 cDNA with pMDC99-AtUBQ10p-GFP and pMDC99-AtUBQ10p-mCherry (Suttangkakul et al., 2011) was performed via the LR Clonase II reaction (Thermo Fisher Scientific). The binary vectors were introduced into Agrobacterium tumefaciens strain GV3101. To obtain stable transformants, Agrobacterium transformants were used to infect Arabidopsis using the floral dip method (Clough and Bent, 1998).
Plant growth conditions
Seeds were surface-sterilized with 50%(v/v) bleach and then washed using sterile water. The seeds were germinated on 1× Murashige & Skoog (MS) solid medium [1× MS salt including vitamins (Duchefa, Haarlem, the Netherlands) with 1% (w/v) sucrose, 0.25% (w/v) phytagel (Sigma, St Louis, MO, USA)] or liquid medium with gentle shaking (100 rpm). Plants were incubated under long-day conditions (16 h light and 8 h dark) at 21-23°C.
To induce autophagy, the seedlings incubated in the liquid medium for 7 d were transferred to nitrogen-deficient medium (1 × MS micronutrient solution, 3 mM CaCl2, 1.5 mM MgSO4, 1.25 mM KH2PO4, 5 mM KCl and 1% sucrose, pH 5.7) and further incubated for indicated duration. For drug treatment, seedlings were treated with 0.5 µM AZD 8055 (LC Laboratories), 30 µM wortmannin (Sigma), and 0.5 µM Concanamycin A (Cayman).
Phenotypic analysis of the fixed carbon starvation response was performed as described (Chung et al., 2010). Briefly, after incubating seedlings on 1× MS-suc solid medium (1× MS salt including vitamins, 29 mM mannitol, and 0.25% phytagel) for 14 d, the seedlings were then deprived of light for 10 d. Seedlings resuming growth were counted after 18 d of recovery.
Immunoblot analysis
Seedlings frozen in liquid nitrogen were homogenized using a TNPI lysis buffer [50 mM Tris-Cl, 150 mM NaCl, 0.5% (w/v) sodium dodecyl sulfate (SDS), 0.5% (v/v) Triton-X 100, 5% (v/v) glycerol, 2 mM dithiothreitol, 1 mM phenylmethyl sulfonyl fluoride, and 10 mM iodoacetamide, pH 8.0]. Lysates were clarified by centrifugation at 16,000 × g at 4°C for 10 min, and the supernatant was mixed with 5 × SDS- polyacrylamide gel electrophoresis (PAGE) sample buffer [200 mM Tris-HCl (pH 6.8), 25% (v/v) glycerol, 10% (w/v) SDS, 10% (v/v) 2-mercaptoethanol]. The samples were boiled at 95°C for 10 min to ensure denaturation of protein samples and loaded into SDS-PAGE. The separated proteins were then transferred onto the Immobilon-P polyvinylidene fluoride membranes (Millipore, St Louis, MO, USA).
To prepare membrane proteins, seedlings expressing GFP-FYVE3 were homogenized with ice-chilled pestle and mortar in TNPI buffer [50 mM Tris-Cl, 150 mM NaCl, 1 mM phenylmethyl sulfonyl fluoride, 10 mM iodoacetamide, pH 8.0]. After centrifugation at 500 × g at 4°C for 5 min, 5% of the supernatant was set aside as total extract, and 95% of the supernatant was further separated into soluble and pellet fractions by centrifugation at 20,000 × g at 4°C for 15 min.
Western blot analysis was performed as previously described (Kim et al., 2022). Anti-GFP (Sigma-Aldrich 11814460001), Anti-H3 (Abcam ab1791), Anti-BiP (Agrisera AS09481), and Anti-UDP-glucose pyrophosphorylase (UGPase; Agrisera AS05 086) antibodies were used. The protein band intensity of the immunoblots was quantified by ImageJ (National Institute of Health).
Quantitative real-time PCR analysis
To evaluate RNA abundance, seedlings were grown on 1 x MS medium for 10 d, frozen with liquid nitrogen, and homogenized in TRIsure agent (Bioline, London, UK). The RNA extract was treated with DNase I (New England Biolabs), and cDNA was synthesized using RevertAid H Minus reverse transcriptase (Thermo Fisher Scientific) with oligo(dT)20 primers. Quantitative real-time PCR (qRT-PCR) was performed using Step One Plus Real-Time PCR System (Applied Biosystems) and EvaGreen (SolGent, Daejeon, Korea). The relative transcript abundance of target genes was calculated by the ΔΔCT method using the UBC9 cDNA as an internal control. All primer sequences are provided in Table S1 .
Yeast two-hybrid analysis
The entry clone containing FYVE3 cDNA was recombined with the destination vector pDEST22 (Thermo Fisher Scientific) via the LR Clonase II reaction. The resulting pDEST22-FYVE3 and pDEST32 expression vectors for bait proteins fused to SAR1B, SAR1C, ATG8A, ATG8E, ATG8F, ATG8I, ATG18A and FYVE2 (Kim et al., 2022) were used for transformation of the AH109 strain and Y187 strain through the Yeastmaker Yeast Transformation System 2 (Clontech, Mountain View, CA, USA). To determine protein interaction, yeast colonies isolated from Trp- and Leu- deficient medium were transferred to Trp-, Leu-, and His- deficient medium containing 0 and 1 mM 3-amino-1,2,4-triazole.
Bimolecular fluorescence complementation
The expression vectors PSAT5-c(175-END)EYFP-C1(B)-SAR1B/ATG8A/ATG18A/ATG6 were described previously (Kim et al., 2022). To generate the expression vectors PSAT4-n(1-174)EYFP-C1-FYVE3 and PSAT5-c(175-END)EYFP-C1(B)-FYVE2/TRAF1A, entry clones containing FYVE3, FYVE2, and TRAF1A cDNA were recombined with the destination vectors PSAT4-n(1-174)EYFP-C1 and PSAT5-c(175-END)EYFP-C1(B), respectively, via the LR Clonase II reaction. The BiFC assay using Arabidopsis leaf protoplasts was performed as previously described (Kim et al., 2022).
Confocal microscopy, image processing, and quantification
To observe fluorescence signals from Arabidopsis plants or leaf protoplasts, either an LSM 510 or AxioObserver LSM 800 confocal microscope (Carl Zeiss) was used. For LSM 510, a 488 nm laser and BP500-530IR emission filter were used to detect GFP signals. For LSM 800, a 488 nm laser was used for excitation and fluorescence was detected at the 410-546 nm range to acquire GFP, citrine, and YFP signals. To detect mRFP and mCherry, a 587 nm laser was used for excitation, and fluorescence was detected at the 595-700 nm range.
Quantitative analysis of the confocal microscopy images was conducted ImageJ (NIH) as previously described (Kim et al., 2022). Data were presented as mean ± standard error. Student’s t-test was used to measure significance via VassarStats (http://www.vassarstats.net/).
Accession numbers
ATG2 (At3g19190), ATG5 (At5g17290), ATG6 (At3g61710), ATG7 (At5g45900), ATG8A (At4g21980), ATG8E (At2g45170), ATG8F (At4g16520), ATG8I (At3g15580), ATG18A (At3g62770), RABF2b/ARA7 (At4g19640), FYVE2/CFS1 (At3g43230), FYVE3 (At1g29800), SAR1B (At1g56330), SAR1C (At4g02080), SYP32 (At3g24350), TRAF1A (At5g43560), UBC9 (At4g27960), and VHA-a1 (At2g28520).
Results
FYVE3 interacts with specific ATG8 isoforms
To gain insight into the FYVE3 protein, we studied its physical interaction with other proteins. An Arabidopsis interactome study identified yeast two-hybrid (Y2H) interactions of FYVE3 with ATG8D (Arabidopsis Interactome Mapping Consortium, 2011). To identify additional interactors, we performed Y2H assay using FYVE3 as a prey and select ATG proteins, FYVE2, and SAR1 isoforms as baits. Our Y2H assay demonstrated interaction of FYVE3 with select ATG8 isoforms such as ATG8A and ATG8F but not with ATG8E and ATG8I ( Figure 1A ).
The positive Y2H interactions were further tested using in planta interaction assays using bimolecular fluorescence complementation (BiFC) in Arabidopsis leaf protoplasts. The BiFC interaction of NYFP-FYVE3 with CYFP-ATG8A was verified ( Figure 1B ), whereas no considerable fluorescence was reconstituted from negative controls of NYFP-FYVE3 paired with either CYFP or CYFP-TRAF1A (Qi et al., 2017) ( Figure 1B ). In addition, interactions of FYVE3 with SAR1B, ATG18A, ATG6, and FYVE2 were detected in BiFC experiments ( Figure 1B ), although these interactions were not apparent in Y2H ( Figure 1A ). We previously showed that FYVE2 interacts with SAR1B and ATG18A (Kim et al., 2022). Thus, FYVE3 may be brought into proximity of SAR1B, FYVE2, and ATG18A, for example, by their simultaneous binding to a scaffold protein. These data suggested that FYVE3 associates with ATG8A, FYVE2, SAR1B, ATG18A, and ATG6 in plant cells.
FYVE3 is recruited from the cytosol to autophagic and endosomal membranes
To investigate the subcellular distribution and trafficking of FYVE3 proteins, we prepared Arabidopsis transgenic plants expressing FYVE3 fused to either GFP or mCherry. The membrane fractionation experiment indicated that GFP-FYVE3 was primarily observed in the soluble fraction ( Figure S2 ). We then determined the subcellular distribution of mCherry-FYVE3 by co-localizing it with various organelle markers ( Figure 2 ). Root cells expressing mCherry-FYVE3 had diffuse fluorescence with relatively scarce punctate signals. The mCherry-FYVE3 puncta often overlapped with the puncta of the PI3P biosensor citrine-2xFYVE ( Figures 2A, G ), the late endosome marker YFP-ARA7 ( Figures 2B, G ), and the autophagic marker GFP-ATG8A ( Figures 2C, G ). However, relatively few mCherry-FYVE3 puncta overlapped with markers for the TGN ( Figures 2D, G ) and Golgi stacks ( Figure 2E, G ).
We also located mCherry-FYVE3 puncta relative to GFP-FYVE2 ( Figure 2F ). GFP-FYVE2 were mostly localized at cytoplasmic puncta, some of which were also decorated with mCherry-FYVE3. Although mCherry-FYVE3 puncta were fewer than GFP-FYVE2 puncta, the former mostly overlapped with the latter ( Figures 2F, G ). As FYVE2 was shown to localize on the autophagic membrane and, to a lesser extent, on the endosomal membrane (Kim et al., 2022), these data indicated that mCherry-FYVE3 puncta mainly represent either autophagic or endosomal compartments that are abundant in PI3P, whereas the diffuse signal corresponds to a cytosolic pool of mCherry-FYVE3.
Association of FYVE3 with autophagic membranes requires ATG8 lipidation and PI3P
The subcellular distribution of FYVE3 is consistent with their protein interaction with select ATG8 isoforms, suggesting that FYVE3 binds to ATG8 on phagophores during autophagosome biogenesis. As fluorescent fusions of FYVE3 are largely co-localized with autophagic puncta, we then examined whether these proteins are delivered together to the vacuole for degradation and detected as autophagic bodies. Observation of the autophagic bodies requires the treatment of GFP-ATG8 transgenic plants with concanamycin A (ConA), an inhibitor of vacuolar proton pumps. Transgenic plants expressing both GFP-ATG8A and mCherry-FYVE3 were treated with ConA, and their root cell images were compared with DMSO-treated controls ( Figure 3A ). In ConA-treated roots, 81% of the intravacuolar mCherry-FYVE3 bodies emitted a GFP-ATG8A signal ( Figure 3A ). ConA-treated roots expressing GFP-FYVE2 and mCherry-FYVE3 also showed a moderate level of co-localization in the vacuole ( Figure 3B ). Similarly, ConA-treated root cells expressing mCherry-FYVE3 alone showed intravacuolar puncta resembling autophagic bodies ( Figure S3A ). These data indicated that fluorescent fusions of FYVE3 are degraded in the vacuolar lumen, together with ATG8A and FYVE2.
To define the molecular mechanism underlying the vacuolar trafficking of FYVE3 fluorescent fusions, we tested the effect of ConA and the PI3K inhibitor wortmannin (Wm) treatment on the subcellular localization of GFP-FYVE3 ( Figure 3C ). As expected, autophagic body-like puncta were observed in the vacuole of ConA-treated wild-type (WT) seedlings expressing GFP-FYVE3 ( Figure 3C , WT, + ConA). Accumulation of the GFP-FYVE3 bodies was suppressed by Wm treatment ( Figure 3C , WT, + ConA, Wm), suggesting that GFP-FYVE3 is also targeted to the vacuole via PI3K-dependent trafficking. To test whether core ATG genes are involved in this trafficking, we introgressed the GFP-FYVE3 transgene into atg2-1 and atg5-1 background by a genetic cross. atg5-1 is defective in ATG8 lipidation (Chung et al., 2010), whereas atg2-1 accumulates a high level of lipidated ATG8 (Kang et al., 2018). When we observed their progenies by confocal microscopy, no autophagic body-like GFP-FYVE3 puncta were observed in ConA-treated atg5-1 and atg2-1 vacuoles ( Figure 3C ), demonstrating that canonical autophagy mediates the vacuolar targeting of GFP-FYVE3.
The vacuolar targeting of GFP-FYVE3 was validated by the immunoblot analysis using anti-GFP antibodies, similar to the GFP-ATG8 and GFP-FYVE2 cleavage assays (Kim et al., 2022). In WT, a very faint band corresponding to GFP-FYVE3 was detected, whereas free GFP released from GFP-FYVE3 was readily observed and its intensity increased by ConA ( Figure 3D ) and dark ( Figure S3B ) treatments, indicating that the free GFP band is a partial degradation product in the vacuole. In atg5-1 and atg2-1, the GFP-FYVE3 band intensity was increased, whereas ConA and dark treatments did not affect the intensity of protein bands (indicated by open diamonds in Figures 3D and S3B ) that migrated slightly slower than free GFP bands (solid arrowheads in Figures 3D and S3B ). Combined with confocal microscopy data ( Figure 3C ), these immunoblot results indicated that GFP-FYVE3 is delivered to the vacuole for degradation via canonical, ATG2- and ATG5-dependent autophagy pathway. Consistently, nitrogen starvation increased the ratio of free GFP to GFP-FYVE3 in WT but not in an atg7-2 background ( Figure S3C ). The intensity of the slower-migrating protein band appeared to increase by nitrogen starvation ( Figure S3C ). Although we do not know the nature of this protein band, the same phenomenon was observed in atg7-2 mutants expressing the autophagy marker GFP-ATG8 (Spitzer et al., 2015). These data indicated that the vacuolar trafficking of GFP-FYVE3 showed a typical pattern of the autophagy marker.
To define the molecular requirement for the formation of GFP-FYVE3 puncta, we examined the effect of atg5 and atg2 mutations and Wm treatment. The abundance of GFP-FYVE3 puncta decreases and increases in atg5 and atg2, respectively ( Figure S3D ). Again, Wm effectively diminished the accumulation of GFP-FYVE3 puncta in atg2 ( Figure S3D ). The PI3P-dependent accumulation of GFP-FYVE3 puncta in atg2 cytoplasm is very similar to that of GFP-ATG8A (Kang et al., 2018), and may result from arrested elongation/maturation of phagophores in atg2.
To investigate the effect of atg2 mutation on mCherry-FYVE3 to the autophagic membrane, we co-localized mCherry-FYVE3 with GFP-ATG8A in WT and atg2 background. atg2 root cells over-accumulated cytoplasmic mCherry-FYVE3 puncta, which extensively overlapped with GFP-ATG8A puncta ( Figure 3E ). Thus, ATG2 is essential for the targeting of mCherry-FYVE3 to the vacuole, but is dispensable for mCherry-FYVE3 recruitment to the phagophore.
FYVE3 is required for the over-accumulation of autophagic structures in fyve2 background
To test the potential functions of FYVE3 in autophagy, we employed fyve3-1, a T-DNA insertional allele (Kim et al., 2022). RT-PCR analysis indicated that fyve3-1 is a knock-out allele, and there is no apparent transcriptional feedback between FYVE2 and FYVE3 genes ( Figure S1B ). Because of the FYVE and SYLF domains shared by FYVE2 and FYVE3, we initially postulated a genetic redundancy of FYVE2 and FYVE3. In fyve2 mutant background, phagophores are over-accumulated and autophagic flux is partially compromised (Kim et al., 2022). However, no significant change in autophagic flux was measured in fyve3 single mutants, and fyve2 fyve3 double mutants did not show further reduction in autophagic flux, when compared with fyve2 single mutants (Kim et al., 2022; Figure S4 ). These data were inconsistent with a simple model that FYVE2 and FYVE3 are duplicates. To further investigate the genetic relationship between these two genes, we quantified the abundance of autophagic organelles in N-starved WT, fyve2, fyve3, and fyve2 fyve3 double mutant cells expressing the autophagic marker GFP-ATG8A. Confocal microscopy analysis confirmed previous observation that fyve2 root cells over-accumulated GFP-ATG8A puncta in the cytoplasm, compared with WT ( Figure 4A ; Kim et al., 2022). As previously reported (Kim et al., 2022), Wm effectively blocked the formation of GFP-ATG8A puncta in fyve2 ( Figure 4A ). In contrast, fyve3 single mutants did not accumulate GFP-ATG8A puncta ( Figure 4A ). Notably, less than 10% of GFP-ATG8A puncta were detected in fyve2 fyve3 double mutants, relative to fyve2 mutants ( Figure 4A ). Complementation experiments revealed that the expression of mCherry-FYVE3 transgene in fyve2 fyve3 mutants restored the accumulation of GFP-ATG8A puncta ( Figure S5 ). These results indicate that FYVE3 is responsible for the over-accumulation of immature autophagic structures in fyve2.
Based on autophagy reporter analysis using fyve3 single and fyve2 fyve3 double mutants, we hypothesized that FYVE3 might serve a specific function during autophagosome formation. This hypothesis was tested using a different mutant combination of atg2 and fyve3 ( Figure 4B ). Compared with fyve2, atg2 mutants accumulate even more autophagic puncta, and autophagic flux is compromised almost completely (Kang et al., 2018; Kim et al., 2022). Unlike fyve2 fyve3 double mutants, atg2 fyve3 double mutant cells showed over-accumulation of GFP-ATG8A puncta as atg2 single mutants did (compare Figure 4B with Figure 4A ). As fyve3 did not suppress the over-accumulation of GFP-ATG8A puncta in atg2, FYVE3 appears to play a specific role in the formation of GFP-ATG8A puncta in fyve2. Of note, treatment with Wm only partially suppressed the GFP-ATG8A puncta in atg2 fyve3 double mutants ( Figure 4B ).
Further genetic analysis, including atg2 fyve2 fyve3 triple mutants, also confirmed that FYVE3 is specifically required for the accumulation of GFP-ATG8A puncta in fyve2 ( Figure 5 ). Under N starvation, atg2 fyve2 fyve3 triple mutants accumulated as many GFP-ATG8A puncta as atg2, atg2 fyve2, and atg2 fyve3 mutants ( Figure 5A ), masking the suppressive effects of fyve3 on fyve2 mutation. A similar trend was observed when seedlings were treated with AZD8055 (AZD; Figure 5B ), which induces autophagy by inhibiting Target Of Rapamycin kinase (Soto-Burgos and Bassham, 2017; Dauphinee et al., 2019).
Accumulation of autophagic puncta can result from either induction of autophagy or retarded biogenesis of autophagosomes. Autophagic flux analysis is necessary to distinguish these two possibilities. As fyve3 alone did not significantly affect autophagic flux (Kim et al., 2022; Figure S4 ), we assessed the effect of fyve3 on autophagic flux in atg2 and fyve2 mutant background. To quantify autophagic flux in different mutants, we used the GFP-ATG8A cleavage assay (Chung et al., 2010; Shin et al., 2014), where anti-GFP immunoblot analysis of seedlings expressing GFP-ATG8A is performed, and autophagic flux is determined from the protein band intensity of free GFP moiety, relative to that of GFP-ATG8A. GFP-ATG8A transgenic seedlings in WT, atg2, fyve2, and fyve3, or their double/triple mutant background were exposed to either N starvation ( Figure 5C ) or AZD ( Figure 5D ) and subjected to GFP-ATG8A cleavage assay. In both conditions, autophagic flux was fully inhibited in any genetic background containing atg2 and also greatly interfered in fyve2 ( Figures 5C, D ). In contrast, fyve2 fyve3 double mutant was indistinguishable from WT and fyve3, confirming that fyve3 suppresses fyve2 ( Figures 5C, D ).
Finally, the physiological function of FYVE3 gene expression was tested by phenotypic analysis of fyve3 single and fyve2 fyve3 double mutants. Autophagy-deficient mutants are hypersensitive to carbon starvation (Hanaoka et al., 2002; Chung et al., 2010). Our fixed-carbon starvation experiments showed that fyve2-2 seedlings showed hypersensitivity to light deprivation, whereas fyve3-1 did not ( Figure 5E ). Importantly, fyve2-2 fyve3-1 double mutants showed a WT-like sensitivity, demonstrating that the suppression effect of fyve3 on fyve2 is physiologically significant.
Discussion
FYVE3 is a putative PI3P effector involved in autophagy (Chung, 2019), but its function remains unelucidated. In this study, we investigated the role of FYVE3 in plant autophagy. We demonstrated that FYVE3 interacts with ATG8 isoforms ( Figure 1 ), associates with ATG8-positive autophagic membranes ( Figure 2 ), and is degraded in the vacuole in a manner dependent on core ATG genes and PI3K ( Figure 3 ). As autophagic flux is not significantly altered by fyve3 single mutation under N starvation ( Figures 4 , 5 ), FYVE3 could be considered an autophagic cargo that does not play a critical role in autophagosome biogenesis. However, this scenario is unlikely, because genetic analysis of fyve2 fyve3 double mutants indicated that fyve3 suppressed fyve2 phenotypes, such as over-accumulation of autophagic puncta, reduction in autophagic flux, and starvation response ( Figures 4 , 5 ). This suppression by fyve3 is in contrast with the hyperaccumulation of autophagic markers in atg2 fyve3 double mutants, where no suppression was observed. Based on these genetic interactions, we propose that FYVE3 acts upstream of FYVE2 but downstream of ATG2 during autophagosome formation ( Figure 5F ).
The phenotypic consequence of fyve3 mutation is not evident unless it is combined with fyve2 mutation, which prevented us and others from identifying its involvement in autophagy (Kim et al., 2022; Zhao et al., 2022). In our previous work (Kim et al., 2022), we employed 48-h nitrogen starvation for autophagic flux assay, whereas 24-h starvation was used for the assay in this study. Nevertheless, we observed a similar suppression in the autophagic flux analysis in which 21-d-old leaves of fyve2 fyve3 mutants were used (Kim et al., 2022). Zhao et al. (2022) reported that the steady-state level of the autophagic cargo receptor NBR1 was comparable between fyve2 and fyve2 fyve3, but did not directly assess the autophagic flux of fyve2 fyve3 (i.e., cfs1 cfs2) double mutants. As NBR1 is transcriptionally regulated in an autophagy-independent manner (Jung et al., 2020), further work is needed to investigate how NBR1 is regulated in fyve2 and fyve3 mutants.
The biochemical function of FYVE3 on the autophagic membrane is not clear. As autophagic flux is almost completely inhibited by atg2 but only partially by fyve2 (Kang et al., 2018; Kim et al., 2022), we postulate two separate autophagic routes, FYVE2-dependent and FYVE2-independent routes ( Figure 5F ). FYVE3 may specifically participate in FYVE2-dependent autophagy. Considering the suppression effect indicating the involvement of FYVE3 in FYVE2-dependent autophagy, the absence of fyve3 single mutant phenocopying fyve2 is noteworthy. This can be explained if FYVE3 is a negative regulator of FYVE2-mediated autophagy. For example, FYVE3 may compete with autophagosome biogenesis factors for ATG8 binding. Alternatively, it is possible that FYVE3 both positively mediates the FYVE2-dependent autophagy and negatively regulates the FYVE2-independent route. In the absence of FYVE3, vacuolar delivery of ATG8 through the FYVE2-independent pathway may increase and compensate for a reduced autophagic flux of ATG8 through the FYVE2-dependent autophagy. In the latter model, the interaction of FYVE3 with ATG8 may be responsible for recruiting FYVE2 on expanding phagophores (Kim et al., 2022). Further investigation is needed to reveal the nature of FYVE2-dependent autophagy and a specific role for FYVE3 in this pathway.
In summary, our study uncovers a regulatory module of FYVE2 and FYVE3 and its potential role in autophagosome biogenesis ( Figure 5F ). Our study suggests that FYVE3 participates in FYVE2-dependent autophagy via its interaction with PI3P and ATG8 isoforms. Thus, future investigation should ask what regions of FYVE3 are responsible for these interactions, whether FYVE3 is a positive or negative regulator of FYVE2-dependent autophagy, and how core ATG proteins and PI3P collaborate with FYVE2 and FYVE3.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
JK performed most experiments. JK, HJ, KS, HL and TC developed genetic materials. TC wrote the manuscript with support from JK, HJ and HL. All authors contributed to the article and approved the submitted version.
Funding Statement
This work was supported by the National Research Foundation of Korea (NRF) grants funded by the Korea government (MSIT) (NRF-2019R1A2C1084203 and NRF-2022R1A2C1003658) to TC.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1160162/full#supplementary-material
References
- Arabidopsis Interactome Mapping Consortium . (2011). Evidence for network evolution in an arabidopsis interactome map. Science 333, 601–607. doi: 10.1126/science.1203877 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Chung T. (2019). How phosphoinositides shape autophagy in plant cells. Plant Sci. 281, 146–158. doi: 10.1016/j.plantsci.2019.01.017 [DOI] [PubMed] [Google Scholar]
- Chung T., Phillips A. R., Vierstra R. D. (2010). ATG8 lipidation and ATG8-mediated autophagy in arabidopsis require ATG12 expressed from the differentially controlled ATG12A and ATG12B loci. Plant J. 62, 483–493. doi: 10.1111/j.1365-313X.2010.04166.x [DOI] [PubMed] [Google Scholar]
- Clough S. J., Bent A. F. (1998). Floral dip: A simplified method for agrobacterium-mediated transformation of arabidopsis thaliana. Plant J. 16, 735–743. doi: 10.1046/j.1365-313x.1998.00343.x [DOI] [PubMed] [Google Scholar]
- Dauphinee A. N., Cardoso C., Dalman K., Ohlsson J. A., Fick S. B., Robert S., et al. (2019). Chemical screening pipeline for identification of specific plant autophagy modulators. Plant Physiol. 181, 855–866. doi: 10.1104/pp.19.00647 [DOI] [PMC free article] [PubMed] [Google Scholar]
- de Jong F., Munnik T. (2021). Attracted to membranes: lipid-binding domains in plants. Plant Physiol. 185, 707–723. doi: 10.1093/plphys/kiaa100 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Dettmer J., Hong-Hermesdorf A., Stierhof Y. D., Schumacher K. (2006). Vacuolar h+-ATPase activity is required for endocytic and secretory trafficking in arabidopsis. Plant Cell 18, 715–730. doi: 10.1105/tpc.105.037978 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Fujiki Y., Yoshimoto K., Ohsumi Y. (2007). An arabidopsis homolog of yeast ATG6/VPS30 is essential for pollen germination. Plant Physiol. 143, 1132–1139. doi: 10.1104/pp.106.093864 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Gao C., Luo M., Zhao Q., Yang R., Cui Y., Zeng Y., et al. (2014). A unique plant ESCRT component, FREE1, regulates multivesicular body protein sorting and plant growth. Curr. Biol. 24, 2556–2563. doi: 10.1016/j.cub.2014.09.014 [DOI] [PubMed] [Google Scholar]
- Gao C., Zhuang X., Cui Y., Fu X., He Y., Zhao Q., et al. (2015). Dual roles of an arabidopsis ESCRT component FREE1 in regulating vacuolar protein transport and autophagic degradation. Proc. Natl. Acad. Sci. U. S. A. 112, 1886–1891. doi: 10.1073/pnas.1421271112 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Geldner N., Denervaud-Tendon V., Hyman D. L., Mayer U., Stierhof Y. D., Chory J. (2009). Rapid, combinatorial analysis of membrane compartments in intact plants with a multicolor marker set. Plant J. 59, 169–178. doi: 10.1111/j.1365-313X.2009.03851.x [DOI] [PMC free article] [PubMed] [Google Scholar]
- González Solís A., Berryman E., Otegui M. S. (2022). Plant endosomes as protein sorting hubs. FEBS Lett. 596, 2288–2304. doi: 10.1002/1873-3468.14425 [DOI] [PubMed] [Google Scholar]
- Hanaoka H., Noda T., Shirano Y., Kato T., Hayashi H., Shibata D., et al. (2002). Leaf senescence and starvation-induced chlorosis are accelerated by the disruption of an arabidopsis autophagy gene. Plant Physiol. 129, 1181–1193. doi: 10.1104/pp.011024 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Inoue Y., Suzuki T., Hattori M., Yoshimoto K., Ohsumi Y., Moriyasu Y. (2006). AtATG genes, homologs of yeast autophagy genes, are involved in constitutive autophagy in arabidopsis root tip cells. Plant Cell Physiol. 47, 1641–1652. doi: 10.1093/pcp/pcl031 [DOI] [PubMed] [Google Scholar]
- Jung H., Lee H. N., Marshall R. S., Lomax A. W., Yoon M. J., Kim J., et al. (2020). Arabidopsis cargo receptor NBR1 mediates selective autophagy of defective proteins. J. Exp. Bot. 71, 73–89. doi: 10.1093/jxb/erz404 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Kang S., Shin K. D., Kim J. H., Chung T. (2018). Autophagy-related (ATG) 11, ATG9 and the phosphatidylinositol 3-kinase control ATG2-mediated formation of autophagosomes in arabidopsis. Plant Cell Rep. 37, 653–664. doi: 10.1007/s00299-018-2258-9 [DOI] [PubMed] [Google Scholar]
- Kim J. H., Jung H., Chung T. (2020). Birth, growth, maturation, and demise of plant autophagic vesicles. J. Plant Biol. 63, 155–164. doi: 10.1007/s12374-020-09252-8 [DOI] [Google Scholar]
- Kim J. H., Lee H. N., Huang X., Jung H., Otegui M. S., Li F., et al. (2022). FYVE2, a phosphatidylinositol 3-phosphate effector, interacts with the COPII machinery to control autophagosome formation in arabidopsis. Plant Cell 34, 351–373. doi: 10.1093/plcell/koab263 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Kim J., Lee H., Lee H. N., Kim S. H., Shin K. D., Chung T. (2013). Autophagy-related proteins are required for degradation of peroxisomes in arabidopsis hypocotyls during seedling growth. Plant Cell 25, 4956–4966. doi: 10.1105/tpc.113.117960 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Lee Y., Kim E. S., Choi Y., Hwang I., Staiger C. J., Chung Y. Y., et al. (2008). The arabidopsis phosphatidylinositol 3-kinase is important for pollen development. Plant Physiol. 147, 1886–1897. doi: 10.1104/pp.108.121590 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Lee H. N., Zarza X., Kim J. H., Yoon M. J., Kim S. H., Lee J. H., et al. (2018). Vacuolar trafficking protein VPS38 is dispensable for autophagy. Plant Physiol. 176, 1559–1572. doi: 10.1104/pp.17.01297 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Liu F., Hu W., Li F., Marshall R. S., Zarza X., Munnik T., et al. (2020). AUTOPHAGY-RELATED14 and its associated phosphatidylinositol 3-kinase complex promotes autophagy in arabidopsis. Plant Cell 32, 3939–3960. doi: 10.1105/tpc.20.00285 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Liu F., Hu W., Vierstra R. D. (2018). The vacuolar protein sorting-38 subunit of the arabidopsis phosphatidylinositol-3-Kinase complex plays critical roles in autophagy, endosome sorting, and gravitropism. Front. Plant Sci. 9. doi: 10.3389/fpls.2018.00781 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Liu C., Zeng Y., Li H., Yang C., Shen W., Xu M., et al. (2021). A plant-unique ESCRT component, FYVE4, regulates multivesicular endosome biogenesis and plant growth. New Phytol. 231, 193–209. doi: 10.1111/nph.17358 [DOI] [PubMed] [Google Scholar]
- Marshall R. S., Vierstra R. D. (2018). Autophagy: The master of bulk and selective recycling. Annu. Rev. Plant Biol. 69, 173–208. doi: 10.1146/annurev-arplant-042817-040606 [DOI] [PubMed] [Google Scholar]
- Phan N. Q., Kim S. J., Bassham D. C. (2008). Overexpression of arabidopsis sorting nexin AtSNX2b inhibits endocytic trafficking to the vacuole. Mol. Plant 1, 961–976. doi: 10.1093/mp/ssn057 [DOI] [PubMed] [Google Scholar]
- Pourcher M., Santambrogio M., Thazar N., Thierry A. M., Fobis-Loisy I., Miege C., et al. (2010). Analyses of sorting nexins reveal distinct retromer-subcomplex functions in development and protein sorting in arabidopsis thaliana. Plant Cell 22, 3980–3991. doi: 10.1105/tpc.110.078451 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Qi H., Xia F., Xiao S. (2021). Autophagy in plants: Physiological roles and post-translational regulation. J. Integr. Plant Biol. 63, 161–179. doi: 10.1111/jipb.12941 [DOI] [PubMed] [Google Scholar]
- Qi H., Xia F. N., Xie L. J., Yu L. J., Chen Q. F., Zhuang X. H., et al. (2017). TRAF family proteins regulate autophagy dynamics by modulating AUTOPHAGY PROTEIN6 stability in arabidopsis. Plant Cell 29, 890–911. doi: 10.1105/tpc.17.00056 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Robertson A. S., Allwood E. G., Smith A. P., Gardiner F. C., Costa R., Winder S. J., et al. (2009). The WASP homologue Las17 activates the novel actin-regulatory activity of Ysc84 to promote endocytosis in yeast. Mol. Biol. Cell 20, 1618–1628. doi: 10.1091/mbc.E08-09-0982 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Schink K. O., Raiborg C., Stenmark H. (2013). Phosphatidylinositol 3-phosphate, a lipid that regulates membrane dynamics, protein sorting and cell signalling. Bioessays 35, 900–912. doi: 10.1002/bies.201300064 [DOI] [PubMed] [Google Scholar]
- Shin K. D., Lee H. N., Chung T. (2014). A revised assay for monitoring autophagic flux in arabidopsis thaliana reveals involvement of AUTOPHAGY-RELATED9 in autophagy. Mol. Cells 37, 399–405. doi: 10.14348/molcells.2014.0042 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Simon M. L., Platre M. P., Assil S., van Wijk R., Chen W. Y., Chory J., et al. (2014). A multi-colour/multi-affinity marker set to visualize phosphoinositide dynamics in arabidopsis. Plant J. 77, 322–337. doi: 10.1111/tpj.12358 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Soto-Burgos J., Bassham D. C. (2017). SnRK1 activates autophagy via the TOR signaling pathway in arabidopsis thaliana. PloS One 12, e0182591. doi: 10.1371/journal.pone.0182591 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Spitzer C., Li F., Buono R., Roschzttardtz H., Chung T., Zhang M., et al. (2015). The endosomal protein CHARGED MULTIVESICULAR BODY PROTEIN1 regulates the autophagic turnover of plastids in arabidopsis. Plant Cell 27, 391–402. doi: 10.1105/tpc.114.135939 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Sutipatanasomboon A., Herberth S., Alwood E. G., Haweker H., Muller B., Shahriari M., et al. (2017). Disruption of the plant-specific CFS1 gene impairs autophagosome turnover and triggers EDS1-dependent cell death. Sci. Rep. 7, 8677–017-08577-8. doi: 10.1038/s41598-017-08577-8 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Suttangkakul A., Li F., Chung T., Vierstra R. D. (2011). The ATG1/ATG13 protein kinase complex is both a regulator and a target of autophagic recycling in arabidopsis. Plant Cell 23, 3761–3779. doi: 10.1105/tpc.111.090993 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Thompson A. R., Doelling J. H., Suttangkakul A., Vierstra R. D. (2005). Autophagic nutrient recycling in arabidopsis directed by the ATG8 and ATG12 conjugation pathways. Plant Physiol. 138, 2097–2110. doi: 10.1104/pp.105.060673 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Urbanek A. N., Allwood E. G., Smith A. P., Booth W. I., Ayscough K. R. (2015). Distinct actin and lipid binding sites in Ysc84 are required during early stages of yeast endocytosis. PloS One 10, e0136732. doi: 10.1371/journal.pone.0136732[doi [DOI] [PMC free article] [PubMed] [Google Scholar]
- Yuan R., Lan J., Fang Y., Yu H., Zhang J., Huang J., et al. (2018). The arabidopsis USL1 controls multiple aspects of development by affecting late endosome morphology. New Phytol. 219, 1388–1405. doi: 10.1111/nph.15249 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Zhao J., Bui M. T., Ma J., Künzl F., Picchianti L., de la Concepcion ,.J. C., et al. (2022). Plant autophagosomes mature into amphisomes prior to their delivery to the central vacuole. J. Cell Biol. 221, e202203139. doi: 10.1083/jcb.202203139 [DOI] [PMC free article] [PubMed] [Google Scholar]
- Zhuang X., Chung K. P., Cui Y., Lin W., Gao C., Kang B. H., et al. (2017). ATG9 regulates autophagosome progression from the endoplasmic reticulum in arabidopsis. Proc. Natl. Acad. Sci. USA 114, E426–E435. doi: 10.1073/pnas.1616299114 [DOI] [PMC free article] [PubMed] [Google Scholar]
Associated Data
This section collects any data citations, data availability statements, or supplementary materials included in this article.
Supplementary Materials
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.