Reliable delivery of presynaptic material, including active zone and synaptic vesicle proteins from neuronal somata to synaptic terminals, is prerequisite for successful synaptogenesis and neurotransmission. However, molecular mechanisms controlling the somatic assembly of presynaptic precursors remain insufficiently understood. We show here that in mutants of the small GTPase Rab2, both active zone and synaptic vesicle proteins accumulated in the neuronal cell body at the trans-Golgi and were, consequently, depleted at synaptic terminals, provoking neurotransmission deficits. Ectopic presynaptic material accumulations consisted of heterogeneous vesicles and short tubules of 40 × 60 nm, segregating in subfractions either positive for active zone or synaptic vesicle proteins and LAMP1, a lysosomal membrane protein. Genetically, Rab2 acts upstream of Arl8, a lysosomal adaptor controlling axonal export of precursors. Collectively, we identified a Golgi-associated assembly sequence of presynaptic precursor biogenesis dependent on a Rab2-regulated protein export and sorting step at the trans-Golgi.
Synaptogenesis is based on the tightly controlled delivery of presynaptic material. Synaptic proteins traffic on presynaptic precursor vesicles from the neuronal soma, where they are produced, to the site of consumption, the synaptic terminal (Ahmari et al., 2000; Chia et al., 2013; Petzoldt et al., 2016; Südhof, 2018). Impairments of presynaptic biogenesis cause severe neurodegenerative or neurodevelopmental diseases, such as Parkinson’s disease, autism spectrum disorder, Alzheimer’s disease, and micro- and macrocephaly (Parrini et al., 2016; Salinas et al., 2008; Waites and Garner, 2011). Different presynaptic cargo proteins are transported along axonal microtubules to the assembling presynapses (Maeder et al., 2014b). This includes active zone (AZ) scaffold proteins, synaptic vesicle (SV) proteins, release factors, and voltage-gated ion channels (Fejtova and Gundelfinger, 2006; Petzoldt et al., 2016; Südhof, 2012). These proteins are thought to traffic in defined stoichiometric ratios as preassembled clusters associated with the precursor membrane to permit effective integration into the nascent presynapse (Ahmari et al., 2000; Goldstein et al., 2008; Siebert et al., 2015; Tao-Cheng, 2007). Despite the fundamental importance of presynapse assembly, the molecular mechanisms underlying presynaptic precursor biogenesis are incompletely understood. The organellar origin and cargo composition of precursors are particularly debated (Chia et al., 2013; Jin and Garner, 2008; Maeder et al., 2014a). Earlier studies suggested that precursors may originate at the Golgi complex (Dresbach et al., 2006; Maas et al., 2012; Shapira et al., 2003). We recently provided evidence that presynaptic precursors possess a lysosomal membrane identity as they acquire lysosomal membrane proteins (LMPs) in a hitherto unexplained maturation process (Vukoja et al., 2018). Transport of these presynaptic lysosome-related vesicles (PLVs) depends on the small Arf-like GTPase Arl8 (Vukoja et al., 2018), a conserved adaptor of kinesin motor proteins implicated in the anterograde transport of lysosomes, lysosome-related organelles (Khatter et al., 2015; Rosa-Ferreira and Munro, 2011), and axonal SV precursors (Klassen et al., 2010).
In the present study, we aim to identify novel regulators of precursor biosynthesis. We hypothesized that precursor biogenesis requires membrane remodeling enzymes, such as Rab (Ras-related in brain) proteins, small GTPases controlling vesicle budding, and fusion via effector protein recruitment (Stenmark, 2009; Zhen and Stenmark, 2015). Rab proteins function as molecular switches by changing between a GTP-bound active state to a GDP-bound inactive state (Stenmark, 2009; Zhen and Stenmark, 2015). We here show that presynaptic precursors originate from the trans-Golgi in an early step of the biosynthetic assembly pathway. This process requires the highly conserved, Golgi-related small GTPase Rab2 for precursor assembly and maturation. Rab2 is reported as a Golgi resident and acts bidirectionally in the ER to Golgi trafficking (Liu and Storrie, 2012; Saraste, 2016; Tisdale and Balch, 1996). In conjunction with the previously shown function of Rab2 and Rab2 effectors in the promotion of dense core vesicle (DCV) maturation in neurons (Ailion et al., 2014; Edwards et al., 2009; Hannemann et al., 2012; Sumakovic et al., 2009) and the association of rab2 mutations with neurodevelopmental defects, e.g., autism spectrum disorders, schizophrenia (Kiral et al., 2018; Takata et al., 2016), and memory and prefrontal morphology defects (Li et al., 2015), in humans, our data identify Rab2 as a crucial factor in early presynaptic precursor biogenesis.
We performed a motoneuron-specific RNAi knockdown of neuronally expressed Rab proteins (Chan et al., 2011) analyzing AZ and SV protein distribution in motoneuron somata and terminals of Drosophila third instar larvae, identifying exclusively Rab2 as a potential regulator of precursor biogenesis.
Presynaptic proteins accumulate in the cell bodies of rab2−/−-deficient neurons
Subsequently, we used a rab2−/− mutant created by ORF excision (Chan et al., 2011) to establish a complete loss-of-function background. Mutant Rab2 protein level was evaluated by Western blot analysis (Fig. S1 A), proving the absence of detectable Rab2 protein in the rab2−/− null mutant. Mutant animals died at the early third larval stage, consistent with rab2−/− null mutants previously described (Lőrincz et al., 2017). We performed immunofluorescence analysis in motoneuron cell bodies located in the cortex of larval ventral nerve cords (VNC) of rab2−/− mutants. We detected a striking ectopic accumulation of the AZ scaffold protein Bruchpilot (BRP; Kittel et al., 2006; Wagh et al., 2006) with a strong increase of BRP protein level (Fig. 1, A and B). BRP did not accumulate diffusely in the cytoplasm, but in distinct, large aggregates not observed in controls (Fig. 1 A, zoom; and Fig. 1 C). We further analyzed the rab2−/− mutants for a set of presynaptic proteins, including another AZ scaffold protein, RIM-binding protein (RIM-BP; Liu et al., 2011). RIM-BP also accumulated in rab2−/− mutant somata (Fig. 1, D and E; and Fig. S1, B and G), and the large, ectopic RIM-BP accumulations overlapped with the BRP signal, while smaller RIM-BP dots also observed in the control were not BRP-positive. The same phenotype was observed for the release factor (m)UNC13A (Böhme et al., 2016; Fig. S2, A and B). Furthermore, SV proteins, here the vesicular glutamate transporter protein (VGlut; Santos et al., 2009), were also strongly accumulating in the rab2−/− motoneuron somata (Fig. 1, F and G; and Fig. S1, C and H). However, in contrast to RIM-BP or UNC13A, VGlut aggregates were 2.5-fold larger than BRP aggregates (Fig. S1 D) and, interestingly, localized rather adjacent to the BRP accumulations. Other SV proteins, Synaptotagmin-1 (Syt-1; Zhou et al., 2017; Fig. S2, C and D) and Dap160/intersectin (Pechstein et al., 2010; Fig. S2, E and F), accumulated with a similar distribution. RNAi mediated knockdown of Rab2 verified by Western blot analysis (Fig. S3 A) consistently provoked significant ectopic accumulations of both AZ and SV proteins (BRP, VGlut, and Syt-1) in motoneuronal cell bodies (Fig. S3, B–H). Unlike SV and AZ proteins, levels of the endogenous mitochondrial ATP synthase, which also traffics from the cell body to the neuronal terminal, were not increased in rab2−/− mutant motoneuronal somata (Fig. S2, G and H), implying that loss of Rab2 could specifically impair presynaptic protein export, while other trafficking pathways are not affected.
We showed previously that axonally trafficking BRP-positive precursors are copositive for the lysosomal-associated membrane protein LAMP1 (Vukoja et al., 2018). We tested for LAMP1 protein levels in rab2−/− mutants by expressing a GFP-LAMP1 construct (Pulipparacharuvil et al., 2005) in motoneurons and observed a strong increase of LAMP1 in rab2−/− mutant cell bodies (Fig. 1, H and I; and Fig. S1, E and I). Similar to SV proteins, LAMP1-positive aggregates were larger and localized rather adjacent to ectopic BRP accumulations. Our findings agree with previous observations showing an increase of LAMP1 signal in larval brains of rab2−/− mutant (Lund et al., 2018). Importantly, the ectopic presynaptic protein accumulations of rab2−/−-deficient cells were not positive for standard markers of the lysosomal and autophagosomal system: Cathepsin L (Fig. 1, J and K; and Fig. S1, F and J), p62 (Fig. S2, I and J), or ATG8a (Fig. S2, K and L). Thus, presynaptic protein accumulations of rab2−/−-deficient neurons are apparently not representing generic degradative compartments.
Finally, we could rescue the BRP accumulation phenotype back to wild-type BRP levels by reexpression of Rab2 at endogenous protein levels into rab2−/− mutants (Fig. 1, L–N), proving the cell-autonomous nature of the rab2−/− phenotype.
Presynaptic biogenesis relies on the Rab2-dependent delivery of presynaptic material
Assuming that the accumulation of presynaptic proteins in the neuronal somata reflected deficits in assembly or transport of precursor cargo proteins, a consequent reduction of presynaptic proteins at motoneuron synaptic terminals is expected. Indeed, when analyzing neuromuscular terminals of rab2−/− mutants, we observed that the entire synaptic terminal of rab2−/− mutants appeared thinner with atypically small boutons, phenocopying larvae lacking kinesin Unc104/KIF1a (Pack-Chung et al., 2007; Zhang et al., 2016), or Arl8 (Vukoja et al., 2018). Overall BRP protein levels at rab2−/− mutant synaptic terminals were strongly reduced compared with wild-type (Fig. 2 A; Fig. 2 Bi; and Fig. 2 C), and AZ numbers were reduced (Fig. 2 D). However, AZ size remained unaltered (Fig. 2 E), indicating that presynapse assembly per se was not affected due to Rab2 absence. Apart from BRP, also other presynaptic proteins were reduced at rab2−/− mutant terminals: RIM-BP (Fig. 2 Bii; and Fig. 2, F and G), UNC13A (Fig. S4, A and E), VGlut (Fig. 2 Biii; and Fig. 2 H), Syt-1 (Fig. S4, B and F), and DAP160 (Fig. S4, C and G). Importantly, mitochondrial ATP synthase was again not altered at synaptic terminals (Fig. S4, D and H), indicating that axonal transport is not generically affected by absence of Rab2. Restricted Rab2 knockdown to motoneuron via RNAi caused a reduction of BRP levels (Fig. S4, I–K) and AZ number (Fig. S4 L) similar to rab2−/− mutants, while AZ size was again unaffected (Fig. S4 M). We further evaluated AZ architecture by superresolution stimulated emission depletion (STED) microscopy. The typical BRP C-terminal ring indicative of regular BRP incorporation and topology of the AZ scaffold (Fouquet et al., 2009) appeared normal in rab2−/− mutants (Fig. 2 I). BRP ring diameter showed a trend for a reduction but was not significantly altered (Fig. 2 J), and AZ number per bouton area was decreased at rab2−/− mutant (Fig. 2 K). We did not detect any diffuse BRP within the presynaptic cytoplasm, as would be expected in the case where BRP cargo would fail to incorporate into AZs (Fig. 2 L). Finally, we performed ultrastructural analyses by EM (Fig. 2 M). Bouton cross-sectional areas showed a trend for reduction in rab2−/− mutants (Fig. 2 N), and T-bar roof length and pedestal height were not significantly altered (Fig. 2, O and P). Moreover, SV numbers in proximity to the T-bar were not changed (Fig. 2 Q). Notably, overall numbers of SVs were reduced (Fig. 2 R), as well as numbers of SVs per bouton area (Fig. S4 N), consistent with a reduction of presynaptic SV proteins delivered (Fig. 2 Biii; and Fig. 2 H). Taken together, data from all resolution levels collectively indicate that Rab2 is primarily implied in the supply of biosynthetic presynaptic material to the synaptic terminal and not in the regulation of “acute” protein incorporation at assembling presynapses.
Two-electrode voltage clamp (TEVC) electrophysiological recordings showed that single action potential evoked response was reduced to half of the evoked excitatory junctional currents (eEJCs) in rab2−/− mutants (Fig. 2, S and T) and, correspondingly, evoked charge transfer (Fig. 2 U), consistent with the AZ numbers being reduced. However, eEJC kinetics (Fig. 2 V) and short-term plasticity (10 ms paired pulse [PP]; Fig. 2 W) were altered. Rab2-RNAi–mediated knockdown yielded similar results to the rab2−/− mutant: eEJC amplitude (Fig. S4, O and P) and eEJC charge (Fig. S4 Q) were significantly reduced, while both kinetics (Fig. S4 R) and short-term plasticity (Fig. S4 S) remained unaltered. Thus, reduction of presynaptic material available at the synaptic terminal in the absence of Rab2 clearly restricts neurotransmitter release. However, those AZs still forming remain apparently functional, in accordance with the overall normal AZ architecture confirmed by microscopy.
Rab2 is a component of anterogradely trafficking presynaptic precursors
Could Rab2, as a membrane anchored GTPase, localize to in vivo trafficking precursors? By use of a genomic Rab2-GFP construct (genRab2-GFP), with the GFP inserted in the Rab2 ORF (Lund et al., 2018), we colabeled larval axons for BRP and observed distinct Rab2-positive puncta, partially copositive for BRP (Fig. 3 A), also shown by line profile (Fig. 3 B) and Pearson’s correlation coefficient (Fig. 3 C). Similarly, we detected puncta copositive for genRab2-GFP and VGlut (Fig. 3, D–F), Syt-1 (Fig. S5, A–C), or the lysosomal marker Spinster (Rong et al., 2011; Sweeney and Davis, 2002; Fig. 3, G–I). Rab2 in a GTP locked constitutively active form (Rab2QL-YFP) expressed in motoneurons showed a similar colocalization profile (Fig. S5, D–L). Next, we directly monitored axonal trafficking in vivo in intact larvae (Andlauer and Sigrist, 2012a), therefore coexpressing strawberry-tagged “BRP short” (BRPD3-straw; Fouquet et al., 2009) and Rab2-YFP or constitutively active Rab2QL-YFP (Fig. 3, J and K; and Fig. S5, M and N). 50% of anterograde trafficking BRP-positive vesicles were copositive for Rab2 (wild-type) and 100% copositive for constitutively active Rab2QL (Fig. 3, J–L). Thus, Rab2 might be activated during transport. Velocities of Rab2/RabQL and BRP-positive vesicles were equal (∼0.79 µm/s for Rab2 and ∼0.69 µm/s for Rab2QL or only BRP-positive vesicles [∼0.78 µm/s; Fig. 3 M]). By contrast, trafficking vesicles devoid of BRP moved more slowly (∼0.55 µm/s Rab2/∼0.43 µm/s Rab2QL). Retrograde trafficking vesicles were in the majority copositive for BRP and either Rab2 construct (93% for Rab2 [wild-type], 72% for Rab2QL; Fig. S5 Q) and slower (∼0.57 µm/s wild-type Rab2/∼0.38 µm/s Rab2QL; Fig. S5 R). Also, coexpression of Spinster-mRFP and Rab2-YFP/Rab2QL-YFP showed that almost all Rab2-positive vesicles were copositive for Spinster, independent of the Rab2 activity status (96% wild-type/100% Rab2QL; Fig. 3, N–P) and moved with a velocity of ∼1.1 µm/s for Rab2 and ∼0.83 µm/s for constitutively active Rab2QL (Fig. 3 Q). Thus, all anterograde trafficking BRP- and Rab2-copositive vesicles are also Spinster-positive, confirming that Rab2, probably in its activated form, is a component of mature presynaptic precursors traveling to the synaptic terminal. Retrograde trafficking Rab2-positive vesicles were also mostly copositive for Spinster (74% for wild-type Rab2 and 80% for Rab2QL; Fig. S5 S) and slower (∼0.66 µm/s for wild-type Rab2 and ∼0.48 µm/s for Rab2QL; Fig. S5 T). As a membrane component of anterogradely trafficking precursor vesicles, we should expect Rab2 to reach the synaptic terminal. Indeed, we observed genRab2-GFP in a punctate pattern at the synaptic terminals, localizing mostly adjacent to the AZ (Fig. 3, R and S). Interestingly, Rab2 might here identify cargo-loaded precursors “waiting to be discharged” and/or vesicles destined for retrograde transport.
In summary, we provided evidence that Rab2 is associated to presynaptic precursors, supporting the hypothesis that ectopic presynaptic protein accumulations in rab2−/− mutants consist of immature, trafficking-incompetent precursors.
Presynaptic precursors derive from the trans-Golgi
Using the genRab2-GFP construct, we now thought to identify the site of Rab2 function during precursor biogenesis. Rab2 localized to large organelles in the cytoplasm of neuronal somata, which were frequently copositive for the Golgi marker GM130 (Fig. 4 A), also shown in the line profile (Fig. 4 B) and by a positive Pearson’s correlation coefficient (Fig. 4 C). We observed a similar Rab2 Golgi localization overexpressing Rab2-YFP (Fig. S6 A), while colabeling with markers for the endo-lysosomal system, such as the early endosomal marker Rabenosyn-5 (Rbsn-5), a RAB5-effector protein, (Fig. S6 B) or the late endosomal marker Rab7 (Fig. S6 C), as well as degradative organelles, using the multivesicular body marker deep orange (DOR/VPS18; Fig. S6 D), the lysosomal marker Cathepsin L (Fig. S6 E), or the autophagosomal marker ATG8a (Fig. S6 F), showed not a similarly strong signal overlap. We thus could validate the previously described Golgi localization of Rab2 (Liu and Storrie, 2012; Saraste, 2016; Tisdale and Balch, 1996) in neuronal somata of Drosophila motoneurons. However, this does not exclude Rab2 being present in other cellular compartments, as reported previously for other cell types (Ding et al., 2019; Fujita et al., 2017; Lőrincz et al., 2017; Lund et al., 2018).
Rab2 is described as a Golgi resident, and its depletion is reported to affect Golgi apparatus integrity (Aizawa and Fukuda, 2015; Maringer et al., 2016). Indeed, the number of Golgi complexes (measured by GM130) was increased by 20% (Fig. 4 E), while Golgi area was decreased by 35% (Fig. 4 F) in rab2−/− mutants, pointing toward a Golgi fragmentation defect. GM130 protein level was not altered (Fig. 4 G). As rab2−/− mutants showed an altered Golgi architecture, it is conceivable that generic protein biogenesis or general post-Golgi sorting defects could cause the ectopic presynaptic protein accumulation we observed at the trans-Golgi. To address this question, we first used RNAi lines targeting proteins of the coat protein II and I (COPII and COPI) complexes, implied in Golgi biogenesis, maintenance, and consequently function (Kondylis and Rabouille, 2009; Popoff et al., 2011; Ward et al., 2001), thus generically disrupting Golgi function. We restricted protein knockdown to the motoneurons to circumvent survival deficits. If generic Golgi disruption accounts for the rab2−/− mutant phenotype, we should observe presynaptic protein aggregation upon Golgi disruption. Depletion of the COPII complex protein Secretory 23 (Sec23), a component of the COPII vesicle coating (Fromme et al., 2008), caused a much stronger Golgi fragmentation than Rab2 depletion (Fig. 4 H), with an 87% increase in GM130 aggregate number (Fig. 4 I), 89% decrease of GM130 aggregate area (Fig. 4 J), and 43% decrease of GM130 intensity (Fig. 4 K). In clear contrast to rab2−/− mutants, we could not observe any ectopic accumulations of presynaptic proteins (BRP, VGlut; Fig. 4, L–O) and no depletion of synaptic proteins (BRP, VGlut) at the synaptic terminal (Fig. 4, P–U). Additionally, depletion of the COPI complex protein ADP-ribosylation factor 79 (Arf79; Donaldson et al., 2005; Popoff et al., 2011) affected Golgi architecture and caused a twofold reduction in GM130 intensity, while GM130 aggregate number and area were not affected (Fig. S6 G). Again, presynaptic proteins (BRP, VGlut) did not aggregate in the somata of Arf79 depleted neurons, although overall BRP and VGlut levels were slightly increased (Fig. S6, H and I). Second, we asked if other than presynaptic proteins were also ectopically aggregating in Rab2-depleted somata and hence analyzed the transsynaptic membrane glycoprotein Neuroglian (Nrg) and the early endosomal marker Rbsn-5. Nrg protein remained bound to the plasma membrane in rab2−/− mutant somata and did not aggregate into ectopic BRP accumulations, although the overall protein level was increased (Fig. 4, V and W). Nrg protein levels at the synaptic terminals were in fact increased, in contrast to the synaptic proteins being depleted in rab2−/− mutants (Fig. 4, X–Z). In rab2−/− somata labeled with Rbsn-5, uncommonly large Rbsn-5–positive organelles formed, which, however, were not copositive for BRP (Fig. S7, A–C). They likely represent enlarged early endosomes reflecting increased fluxes from the Golgi to early endosomes. Consistently, Rbsn-5 levels at Rab2-deficient terminals were not altered (Fig. S7, D and E). Thus, in the absence of Rab2, nonsynaptic proteins are still targeted to their membrane destination, and sorting from the Golgi is apparently not compromised, confirming that Rab2 specifically regulates presynaptic protein sorting from the trans-Golgi.
We next analyzed the subcellular localization of ectopic BRP-positive aggregates in relation to the Golgi in rab2−/−-deficient somata and colabeled with a cis-Golgi (GM130) or trans-Golgi (Syx16) marker. Interestingly, BRP accumulations localized in the majority adjacent to the cis-Golgi (Fig. 5 A, arrowheads), with peak-to-peak distances of typically only around 400 nm (Fig. 5 B), and overlapped with the trans-Golgi (Fig. 5, C and D). We correspondingly observed a similar distribution of RIM-BP– (Fig. S7, F–I) and VGlut-positive (Fig. 5, E–H) aggregates with relation to both Golgi markers and could confirm these findings by STED microscopy. BRP accumulations also, under elevated magnification, formed adjacent to the cis-Golgi consistently both in lateral and planar views (Fig. 5 I) and clearly localized to the trans-Golgi area (Fig. 5 J; quantified in Fig. 5 K).
The trans-Golgi is the site of protein and lipid sorting and export for post-Golgi destinations (Keller and Simons, 1997; Klumperman, 2011; Pols et al., 2013). Assuming that in an early step of precursor biogenesis, synaptic proteins could sort in a Rab2-dependent process from the trans-Golgi, we should be able to detect precursors at the Golgi of wild-type neurons. Indeed, we observed BRP-positive puncta at the Golgi of wild-type somata (Fig. 5 L, arrowheads). Their numbers were higher at the Golgi compared with random-shifted Golgi areas in the cytoplasm, as well as BRP mean pixel intensity (Fig. 5, M and N). Notably, we also observed BRP-positive puncta in the cytoplasm (Fig. 5 L), and it is tempting to speculate that they might represent post-Golgi precursors ready to be exported from the cell body.
Thus, accumulating presynaptic proteins at the trans-Golgi of rab2−/− mutant somata could represent immature precursors, interrupted in their biosynthetic maturation process due to loss of the membrane trafficking regulator Rab2.
Tubule-shaped vesicular membranes accumulate at the trans-Golgi in rab2−/− mutants
For a detailed evaluation of the presynaptic protein accumulations in rab2−/− mutants, we performed ultrastructural analysis by EM. Electron micrographs showed a striking accumulation of circular- to oval-shaped vesicular membranes at the trans-Golgi (Fig. 6, A and B; and Fig. S8, A and B), correlating with the position of presynaptic protein accumulations detected by light microscopy. No vesicle accumulations of this kind were observed at the Golgi in wild-type cells (Fig. 6, A and B; and Fig. S8, A and B), and the volume fraction of trans-Golgi network–associated vesicles of rab2−/− mutants compared with wild type was significantly increased (Fig. 6 C).
But how do these maturation-arrested, immature precursors compare with the PLVs we described previously to accumulate in the cell body of arl8−/− mutants (Vukoja et al., 2018)? The direct comparison of arl8−/− mutant and rab2−/− mutant electron micrographs revealed that PLVs detected in the absence of Arl8 were, with a ∼60-nm short diameter, significantly larger than rab2−/− mutant precursors at ∼40 nm, possessed a more electron dense core (Fig. 6, B and D), and had a clearly more uniform, circular shape. Long diameters of both mutants were similar at ∼75 nm (Fig. S8 C). Notably, the volume fraction taken by mitochondria (Fig. 6 E) was not altered in Rab2-deficient neurons, in agreement with the confocal data (Fig. S1, Q and R), again showing that absence of Rab2 does not affect anterograde transport generically. Interestingly, we detected a threefold increase in the incidence of somatic structures resembling autophagocytic organelles in rab2−/− mutants (Fig. S8 D), potentially interesting in the context of the previously described Rab2 function in the autophagic degradative pathway (Ding et al., 2019; Fujita et al., 2017; Lőrincz et al., 2017; Lund et al., 2018).
We performed 3D electron tomography to evaluate more deeply the apparent shape differences of precursor vesicles in rab2−/− and arl8−/− mutants. Clearly, immature precursors in rab2−/− mutants localized, in agreement with the confocal and STED light-microscopic data (Fig. 5, A–H), along the entire trans-Golgi, while the cis-Golgi facing the nucleus was devoid of such vesicular accumulations (Fig. 6 F). Immature precursors were detectable now in the 3D reconstruction as heterogeneous vesicular or short tubular structures (Fig. 6, G and H), with an average short/long diameter ratio of ∼1.8 (Fig. 6 I). By contrast, PLVs accumulating in arl8−/− mutants were rather homogenously shaped, circular vesicles with a short/long diameter ratio close to 1 (Fig. 6, H and I). Consequently, immature precursors of rab2−/− mutants were threefold smaller compared with those of arl8−/− mutants (Fig. 6 J) and had a twofold smaller surface area (Fig. 6 K).
In summary, we suggest that the 40-nm-size, tubule-shaped, clear-core vesicles forming at the trans-Golgi of rab2−/− mutants represent immature and biosynthetically early arrested presynaptic precursors.
Convergence of Golgi “trafficking routes” during maturation of presynaptic precursors
The confocal microscopy analysis of ectopic presynaptic protein aggregates forming in rab2−/− mutants had shown that not all investigated proteins overlapped with the BRP signal (Fig. 1), suggesting that variable Golgi-sorting routes for different presynaptic proteins might exist. Hence, we performed 3D confocal image reconstruction to scrutinize the heterogeneous presynaptic protein distribution in more detail. A 360° rotation of a 3D reconstructed precursor field costained for BRP and VGlut demonstrated that BRP- or VGlut-positive precursors occupied contiguous instead of overlapping areas (Fig. 7 A), clearly visible in the line profiles (Fig. 7 B) and indicated by a low Pearson’s correlation coefficient (Fig. 7 C). STED imaging further underlined this finding: BRP-positive or VGlut-positive precursors occupied neighboring areas (Fig. 7 Di; and Fig. 7 Ei, for line profile), and quantification of the area center distance of both precursor fields showed a clear shift of several hundred nanometers (Fig. 7 H). By contrast, BRP- and RIM-BP–positive areas overlapped precisely (Fig. 7 Dii; and Fig. 7 Eii), and center distances between them were nearly zero (Fig. 7 H). Interestingly, GFP-LAMP1 preferentially localized to the VGlut-positive rather than the BRP-positive aggregate area (Fig. 7, Diii and Div; and Fig. 7, Eiii and Eiv), also evident in the large center distance between BRP and LAMP1 areas compared with VGlut and LAMP1 distances (Fig. 7 H). Consistently, GFP-LAMP1 accumulations localized in close proximity to the cis-Golgi (GM130) in rab2−/− mutants (Fig. S8, E and F), similar to VGlut (Fig. 5 E). Notably, a subfraction of the ectopic LAMP1 accumulations in rab2−/− mutants was devoid of BRP signal and not associated with early presynaptic precursor fields (Fig. S8, G and H), potentially related to the function of Rab2 in autophagic and lysosomal degradation (Ding et al., 2019; Fujita et al., 2017; Lőrincz et al., 2017; Lund et al., 2018). These data suggest that presynaptic proteins sort in a Rab2-dependent step from the trans-Golgi on two independent routes, separating AZ scaffold from SV proteins.
Interestingly, this segregation of presynaptic cargo was not observed in arl8−/− mutants (Vukoja et al., 2018). We here applied STED microscopy on arl8−/− mutants for comparability. Also, upon higher resolution, both VGlut and LAMP1 localized to the same area as BRP (Fig. 7, Fi and Fii; and Fig. 7, Gi and Gii) with a small center distances (Fig. 7 H), clearly contrasting with rab2−/− mutants. Finally, we also expressed a GFP-tagged Arl8 (Vukoja et al., 2018) in rab2−/− mutant animals and found that Arl8 was already present on both BRP- and VGlut-positive immature precursors (Fig. S8, I–L), as BRP- and VGlut-positive aggregates showed a nearly 100% overlap with Arl8 (Fig. S8, M and N). Collectively, these data suggest that in a first biosynthesis step, presynaptic proteins sort from the trans-Golgi on at least two distinct routes, a process regulated by the small GTPase Rab2. Immature precursors converge to larger and circular PLVs compatible with Arl8-dependent transport in subsequent maturation steps, characterized by a uniform presynaptic cargo load and a lysosomal membrane identity. This would suggest that Rab2 functions upstream of Arl8 during precursor biogenesis.
Rab2 acts upstream of Arl8 in presynaptic precursor biogenesis
We created a rab2−/−/arl8−/− double mutant and compared the somatic presynaptic protein accumulation phenotype to the respective single mutants to access the epistatic relation of both proteins. Both number and intensity of BRP or VGlut accumulations were increased in arl8−/− compared with rab2−/− mutants (Fig. 8, A–E). Importantly, the rab2−/−/arl8−/− double-mutant larvae clearly resembled the phenotype of rab2−/− single mutants, with quantifications detecting no significant differences between single rab2−/− and rab2−/−/arl8−/− double mutants (Fig. 8, A–E). The segregation of BRP and VGlut within the ectopic accumulation of rab2−/− mutants was clearly detectable in rab2−/−/arl8−/− double mutants, while arl8−/− mutants again showed overlapping BRP and VGlut signals (Fig. 8 A). STED microscopy confirmed the BRP/VGlut segregation of rab2−/−/arl8−/− double mutants (Fig. 8, F–I), and quantification showed an equal center distance between rab2−/− single and rab2−/−/arl8−/− double mutants, while single arl8−/− mutants had a significantly reduced center distance (Fig. 8 J). Thus, disruption of precursor biogenesis in the absence of Rab2 prevents the entrance in a later Arl8-dependent biosynthetic maturation step, placing Rab2 upstream of Arl8. Consequently, when analyzing the synaptic terminals of rab2−/−/arl8−/− double mutants, no additive effect in the reduction of the synaptic material was observed (Fig. S8, O and P), consistent with our interpretation that both proteins indeed act in a sequential manner.
Although of fundamental importance, our understanding of the cell biological origin and molecular mechanics of presynaptic precursor biogenesis remains fragmentary. We here provide evidence that early precursor biogenesis is regulated by the conserved, Golgi-related small GTPase Rab2 (Chan et al., 2011; Lardong et al., 2015; Liu and Storrie, 2012; Lőrincz et al., 2017; Saraste, 2016; Tisdale and Balch, 1996). We suggest that precursor formation takes place in two consecutive steps: (1) sorting of precursor cargo proteins from the trans-Golgi on independent exit routes for AZ and SV proteins, (2) followed by a subsequent maturation step merging both exit routes. Convergence could occur either by direct fusion of early precursors or through conventional sorting at endo-lysosomal compartments, from which mature precursors are retrieved at a later time point (Fig. 8 K).
Several independent lines of evidence, including confocal, super-resolution light, and EM, live imaging, and genetic assays, support this model. Upon depletion of Rab2, biosynthetic presynaptic material accumulates ectopically at the trans-Golgi of motoneuronal somata, overlapping with the trans-Golgi marker Syx16 (Figs. 5 and 6). Interestingly, these accumulations did not show a uniform presynaptic material distribution but a segregation into two adjacent but not overlapping fields, positive for either presynaptic proteins (BRP, RIM-BP, UNC13) or SV proteins (VGlut, Syt-1; Figs. 1 and 7). These aggregates apparently do not represent degradative compartments as they were devoid of lysosomal and autophagosomal machinery markers, i.e., Cathepsin L, p62, and ATG8a (Figs. 1 and S1). Ultrastructural analysis revealed that ectopic accumulations consist of tubule-shaped, clear-core vesicles with a short diameter of ∼40 nm (Fig. 6). The elongated structures frequently appeared to consist of two to three attached spherical modules (Fig. 6 H), possibly consisting of immature precursors unable to proceed in a Rab2-dependent fusion or fission process. By contrast, the previously identified PLVs observed in arl8−/− mutants (Vukoja et al., 2018) were uniformly circular-shaped with ∼60-nm diameters, a consequently larger volume, and a dense core (Fig. 6), apparently representing more matured precursors. Moreover, the matured precursor field of arl8−/− mutants showed no segregation of AZ and SV proteins and is thus likely the result of a subsequent fusion step of the post-Golgi early precursors, either directly by mutual fusion or via the endo-lysosomal compartments (Fig. 8 K). As this fusion step evidently does not occur in the absence of Rab2, immature precursors here have lost the ability to fuse to mature vesicles. These findings would place Arl8 function downstream of Rab2 in the precursor biogenesis pathway and, indeed, we could show by genetic epistasis experiments that Rab2 acts upstream of Arl8 (Fig. 8).
Rab2 as a membrane-associated protein could function either directly at the trans-Golgi or on the arrested immature precursors by recruitment of Rab2 interaction partners to the nascent precursor membranes executing the subsequent fusion process, for example, SNARE proteins or proteins of the homotypic fusion and protein sorting complex, such as VPS39, already shown to be a direct interaction partner of Rab2 (Gillingham et al., 2014).
Interestingly, both the SV protein–positive early (rab2−/− mutants) and the matured (arl8−/− mutants) precursors are positive for the lysosomal marker LAMP1 (Fig. 7). Thus, the immature precursors already express some lysosomal membrane identity. LMPs traffic on nonclathrin-coated vesicles from the trans-Golgi network to the late endosome (Pols et al., 2013). It is tempting to speculate that LAMP1-positive early precursors might be related to these organelles. In contrast, Arl8, the lysosomal kinesin adaptor, localized to both types of immature precursors in rab2 mutants (Fig. S5). Arl8, as a regulator of lysosome and lysosome-related organelle mobility (Khatter et al., 2015; Rosa-Ferreira and Munro, 2011), could confer motility to the mature precursors in the presynaptic precursor pathway by hooking them to the microtubule motor kinesin, as previously suggested in fly, mouse (Vukoja et al., 2018), and Caenorhabditis elegans (Klassen et al., 2010).
Downstream of somatic precursor biogenesis, Rab2 apparently remains associated, probably in its activated form, to the mature precursors as it trafficked anterogradely on BRP/Spinster-positive vesicles and effectively reached the synaptic terminal (Figs. 3 and S3). Retention of presynaptic proteins in the neuronal soma of rab2−/− mutants consequentially caused a depletion of these proteins at the synaptic terminal and compromised presynapse formation and synaptic transmission (Fig. 2). The few synapses still forming remained roughly unaltered regarding architecture and SV number at the presynapse, suggesting a Rab2 function in the supply of presynaptic material rather than incorporation of precursor material into the assembling presynapse. However, as Rab2 is localized to the actively transported precursors (Fig. 3), it is conceivable that next to its upstream function in precursor biogenesis, Rab2 also directly contributes to the trafficking of presynaptic precursors.
Rab2 is described as a Golgi-related small GTPase implied in the bidirectional ER to Golgi trafficking (Liu and Storrie, 2012; Saraste, 2016; Tisdale and Balch, 1996), and indeed, we consistently observed subtle Golgi fragmentation defects in rab2−/− mutants (Fig. 4, D–G), raising the possibility that generic Golgi sorting or protein biogenesis defects might be the primary cause of the rab2−/− mutant phenotype. However, we could show that effective generic Golgi disruption did not phenocopy the rab2−/− mutant phenotype, and other nonpresynaptic proteins were not accumulating at the trans-Golgi (Figs. 4 and S5), confirming a specific role for Rab2 in presynaptic precursor biogenesis.
The here-described function of Rab2 in precursor biogenesis could share mechanistic traits with the role of Rab2 and its effectors (RIC-19, RUND-1, CCCP-1, TBC-8) in DCV maturation in C. elegans (Ailion et al., 2014; Edwards et al., 2009; Hannemann et al., 2012; Sumakovic et al., 2009). DCVs store, transport, and release neuropeptides and neurotrophins through exocytosis. They form through secretion from the trans-Golgi network (Kim et al., 2006; Morvan and Tooze, 2008), where inadequate secretory proteins become incorporated into the immature DCVs, which are subsequently removed by clathrin-mediated membrane trafficking into the endo-lysosomal system (Ailion et al., 2014; Edwards et al., 2009; Kim et al., 2006; Morvan and Tooze, 2008; Sumakovic et al., 2009). Elegant studies in C. elegans showed that Rab2 and its effectors are required for the retention of neuropeptides in immature DCVs. Although we can only speculate, Rab2 function during the early phase of precursor biogenesis could comprise a similar mechanism by the retention or sorting of presynaptic and/or lysosomal proteins. Furthermore, future studies also need to evaluate a potential role of Rab2 in DCV maturation in Drosophila.
Our study contributes to a comprehensive model of the biosynthetic pathway underlying presynaptic precursor biogenesis and could support future research by connecting Rab2-related neurodevelopmental defects, e.g., memory defects in human (Li et al., 2015), autism spectrum disorders, or schizophrenia (Kiral et al., 2018; Takata et al., 2016), with the Golgi pathway–related neurodegenerative diseases (Rasika et al., 2018).
Materials and methods
Drosophila melanogaster strains were reared under standard laboratory conditions and raised at 25°C and 70% humidity on semi-defined medium (Bloomington recipe). For RNAi experiments, flies were kept at 29°C. For electrophysiological recordings, only male larvae were used. For all other experiments, both male and female animals were used. For genotypes and fly strains used, see Tables 1 and 2.
We used the genomic rab2−/− mutant line (Chan et al., 2011), in which the Rab2 ORF was replaced by a Gal4, now transcriptionally regulated by the endogenous Rab2 transcription regulators, and thus expressed a UAS-Rab2-YFP at endogenous protein levels in the rab2−/− mutant background (see also Tables 1 and 2).
Immunostainings of larval neuromuscular junctions (NMJs) and brains
Immunohistochemistry for confocal and STED microscopy
For immunohistochemistry, dissections were performed in hemolymph-like solution 3 (HL3; Stewart et al., 1994; composition in millimolar: 70 NaCl, 5 KCl, 20 MgCl2, 10 NaHCO3, 5 trehalose, 115 sucrose, and 5 Hepes, pH adjusted to 7.2) by opening the third instar larvae dorsally along the midline and removing the entrails. Dissections were fixated with 4% paraformaldehyde in PBS (pH 7.2) for NMJ staining or Bouin’s fixative for staining of brain samples for 10 min. After fixation, the filets were washed with PBS plus 0.05% Triton X-100 and blocked for 60 min in 5% normal goat serum (Sigma-Aldrich; S2007). For immunostainings, the larvae were incubated with primary antibodies at 4°C overnight and subsequently washed in a PBS plus 0.05% Triton X-100 solution for 2 h at RT. Larvae were then incubated for 2–3 h with secondary antibodies at RT. Washing procedures were repeated. Larvae were finally mounted in Vectashield (Vector Laboratories), Mowiol (Sigma-Aldrich), or ProLong Gold (Thermo Fisher Scientific). Primary and secondary antibodies used in these studies are shown in Table 3.
Image acquisition and analysis
Conventional confocal and STED images were acquired with Leica DMI 6000 (SP8) and TCS SP8 gSTED 3× microscopes (Leica Microsystems), respectively. For confocal scans, an HC PL APO CS2 63× /1.40 NA oil objective (Leica Microsystems) was used; for STED, an HC PL APO CS2 100×/1.40 NA oil objective (Leica Microsystems) was used. Images were acquired at ∼20°C, and fluorochromes are indicated in the antibody section described above. Imaging medium was immersion oil (Sigma-Aldrich; 55822). Mounting medium was Mowiol 4-88 (Sigma-Aldrich). For signal detection, high-sensitive hybrid detectors with a 400–800-nm spectrum scanned for green and red channels and photomultiplier tubes with a 400–800-nm spectrum scanned for far red channels and were used for confocal scans. For STED, hybrid detector Sp GaAsP were used. Larval brain z-stacks had a step size of 0.3–0.5 µm between single optical slices. 40 z-slices of the central region along the dorsoventral axis were imaged. The NMJ z-stacks had a step size of 0.2–0.3 µm between single optical slices. All images were acquired using the LAS X software (Leica Microsystems). The ImageJ 1.52n software was used for the analyses of confocal and STED images. GraphPad Prism, v5.01, was used for statistical analyses. For STED microscopy, Huygens Deconvolution software was used applying a theoretical point spread function automatically computed based on pulsed- or continuous-wave STED-optimized function and the specific microscope parameters. Default deconvolution settings were applied. Images for figures were processed, if necessary, with ImageJ software to enhance brightness using the brightness/contrast function and smoothened (0.5-pixel Sigma radius) using the Gaussian blur function. Confocal stacks were processed with the image software Fiji (http://fiji.sc; Schindelin et al., 2012). Image analysis followed the standard protocol as described by Andlauer and Sigrist (2012b) and will be described in detail for the different analyses below.
NMJ and VNC quantifications
Type 1b NMJs on muscle 4 were analyzed. Zoom areas in the VNC were acquired of motorneuronal somata either lateral or central (dorsal) of the neutrophil. The HRP-Cy5 antibody signal was used as the template for a mask, restricting the quantified area to the shape of the NMJ. For larval brain aggregate quantification, a region of interest (ROI) was drawn via the freehand selection tool including only the cortex region of the larval VNC (sparing neuropil region and background). The original confocal stacks were converted to maximal projections (for brain, taking only the central 20 z-slices of the VNC without adjacent axons or weakly stained tissue in deeper layers), and after background subtraction, a mask of the synaptic area/brain aggregates was created by applying a manual threshold to remove the irrelevant lower-intensity pixels. The segmentation of single spots was done semi-automatically via the command “Find Maxima” embedded in the Fiji software and by hand with the pencil tool and a line thickness of 1 pixel. To remove high-frequency noise, a Gaussian blur filter (0.5-pixel Sigma radius) was applied. The processed picture was then transformed into a binary mask using the same lower threshold value as in the first step. This binary mask was then projected onto the original unmodified image using the “min” operation from the ImageJ image calculator.
The number of AZ or VNC aggregates in the graphs reflects the density of particles in both tissues and was normalized to either the NMJ (HRP) area or the freehand selection ROI of the cortex area. The AZ area (nm2) is the mean area value of all determined individual AZ area of one NMJ. For sum intensities “per NMJ area,” every single pixel intensity within the single spots/particles of the corresponding channel was first added up (calculated as “integrated density” per AZ/aggregate in Fiji), and subsequently all spot values per NMJ/cortex were summed up in order to estimate total protein levels of all particles within the respective tissue. Finally, this total sum was divided by the respective NMJ or cortex area to normalize these estimates for variations in imaged tissue size. n represents the number of NMJs/larval brains. For sum intensities “of NMJ,” protein mean pixel intensity within the HRP mask was calculated as “integrated density.” For protein mean pixel intensity “of AZ,” the mean pixel intensity of the individual AZs was determined and the mean value plotted. For protein mean intensity “of cortex area,” an ROI around the whole cortex area excluding the neuropile was created using the HRP signal as a reference. Mean pixel intensity of protein within this ROI was determined and the mean value plotted. The same procedure was applied for NMJ (“protein mean pixel intensity of NMJ”) using the HRP mask. All data were normalized to wild-type (100%). n for brain quantifications represents one VNC of one larva; four to eight larvae were analyzed. n for NMJs represents NMJs with one or two NMJs/animal and five or six animals/genotype.
Golgi (GM130) quantification of Golgi defects
Golgi quantifications were performed as described above regarding the basic image treatment. A stack of four to six slices from the VNC-zoom images was chosen and a maximum projection performed. Three to six clearly detectable cells were chosen for analysis and marked (ROI-area). Within this ROI, single GM130 spots were identified as described above, but without segmentation. The number of spots was divided by ROI-area and normalized to wild-type (100%). Area of the individual GM130 spot was determined, and the mean plotted, normalized to wild-type (100%). Mean pixel intensity of GM130 signal within the ROI-area was determined, and the mean plotted, normalized to wild-type (100%). n represents one VNC of one larva; four to eight larvae were analyzed.
Pearson’s correlation coefficient
For Pearson’s correlation coefficient, an image stack was acquired, and one focal plane was chosen for analysis after ROI selection within a single cortex soma or axonal bundle, respectively. ROIs were selected around aggregates with potential signal overlap of one to three aggregates in a range of 1.5–2.5 µm2. ImageJ Coloc_2 was used to determine the Pearson’s R value above threshold by testing the Rab2 channel against the other indicated second channel (BRP, VGlut, Spinster, etc.). As a control, the second channel was flipped horizontally (“inverted BRP, Spinster, VGlut, etc.”), and subsequently the Pearson’s R value above threshold was recalculated. This aims to show specific colocalization only in the original (not flipped/not inverted) overlapping channels and to exclude a general stochastic colocalization. Analyzed where overall 9–12 motoneuronal somata (one ROI per soma) of two or three VNCs, corresponding to two or three animals. For axonal colocalization, 10–15 axonal bundles (one ROI per axonal bundle) in three to five animals were selected.
BRP mean pixel intensity in boutons excluding AZs from STED analysis of NMJs
Bouton area was outlined manually (ROI-1), and outside bouton area was cleared. Gaussian blur with 2.0 was applied. A threshold of 30 was set and all AZs within this threshold marked (ROI-2). Subtraction of ROI-2 from ROI-1 (XOR function in Fiji ROI manager) and the new ROI was named ROI-3. In ROI-3 (bouton area without AZs), the mean pixel intensity was measured. Data were collected from two or three boutons of one NMJ per larva from three larvae.
Determination of the number of AZs per bouton area and AZ ring diameter from STED analysis of NMJs
Area overlap measurement of BRP/VGlut and Arl8-GFP in percentage from neuronal somata
Single cells were manually chosen from three to five images/brain, three brains in total. Threshold values were set manually, and ROIs in both channels were outlined and the area measured. The AND function of the Fiji ROI manager was used to determine area overlap. Only the Arl8-positive area was determined by subtraction of the copositive area (AND; see above) from the total Arl8 area and percentages calculated.
Frequency determination in percentage of peak-to-peak distance of BRP-positive aggregates in rab2−/− mutants and GM130-marked Golgi
We place line through BRP aggregates and Golgi, saved as ROI, used ROI manager MultiPlot tool, and determined peak-to-peak distance. For random shifted data, the BRP channel was shifted by 80 pixels. A cutoff at 1.2-µm distance was used, binned for 0.06 µm. We measured 1–14 distances from two VNC images per larva, with a total of three larvae analyzed.
Center distance measurements of STED-imaged presynaptic protein aggregates in rab2−/− mutant somata
In vivo live imaging and analysis
In vivo imaging of intact Drosophila larvae was performed as previously described (Andlauer and Sigrist, 2012a). Briefly, third instar larvae were put into a drop of Voltalef H10S oil (Arkema) within an airtight imaging chamber. The larvae were anesthetized before imaging with 10 short pulses of a desflurane (UAS) air mixture until the muscles relaxed. Axons were imaged using confocal microscopy immediately after exiting the VNC.
Kymographs were plotted using a custom-written Fiji script (based on tsp050706.txt TimeSpacePlot [Kymograph] from Jens Rietdorf, downloadable at http://www.embl.de/eamnet/html/kymograph.html). In brief, temporal recordings of proximal axons (pixel size, 100 nm; t resolution, 0.8–1.7 s), expressing YFP/mRFP/strawberry-labeled proteins of interest, were background-subtracted, filtered via a Gausian blur (Sigma radius 0.5), and stabilized via the Fiji image registration plugin MultiStackReg. A time stack limited to the time frame of a single specific cotransport event was cropped from the original file. Sequentially, a segmented line ROI was drawn along the trajectory of the cargo (minimum of 30 s stable visibility). The final kymograph was (automatically) created by adding the line ROI (average width, 3–10 pixels) per time frame on top of the preceding line ROI in a concatenated fashion.
Quantification of (co)transport events in proximal axons was performed as follows. Per indicated coexpression genotype, two or three different axonal nerve bundles (5-min imaging time per axon) were imaged in five or six animals. Every file was screened for transport events by quantifying events of sufficient intensity and traceability (stable identification of the unidirectional cargo trajectory for at least 4 µm in one channel) via drawing a linear ROI starting from the point of first discovery to the end point in the last sufficiently recognizable section. The length of this ROI was divided by the time (number of traceable images of the single event × t resolution per file), and the resulting speed value (µm/s) was listed according to the cargo direction (anterograde/retrograde). Every ROI (created in channel one) was compared with the corresponding coexpression channel and only counted as cotransport if this second channel indicated a similar unipolar movement event in various images (at least three or four different temporal sections within the same time frame while indicating a sufficient signal overlap).
TEVC recordings were performed as previously reported (Matkovic et al., 2013). In short, spontaneous recordings (miniature excitatory junction currents [mEJCs], 90 s), single-evoked recordings (eEJCs, 20 repetitions at 0.2 Hz), and high-frequency recordings (PP 10-ms or 30-ms interstimulus interval, PP10 or PP30, 10 repetitions at 0.2 Hz) were recorded. All experiments were performed on third instar larvae raised at 25°C on semi-defined medium (Bloomington recipe). Dissection and recording medium was extracellular HL3 (Stewart et al., 1994; in millimolar: 70 NaCl, 5 KCl, 20 MgCl2, 10 NaHCO3, 5 trehalose, 115 sucrose, and 5 Hepes, pH adjusted to 7.2). Dissections were done in ice-cold Ca2+-free HL3 medium, while mEJC, eEJC, and high-frequency recordings were performed in 1.5 mM Ca2+ HL3 at room temperature. For all physiological recordings, intracellular electrodes with a resistances of 15−25 MΩ (filled with 3 M KCl) were placed at muscle 6 of abdominal segment A2/A3. Acquired data were subjected to a low-pass filter at 1 kHz and sampled at 10 kHz. Traces of both miniature and evoked postsynaptic currents were recorded in TEVC mode (AxoClamp 2B; Axon Instruments). The command potential for mEJC recordings was −80 mV, and it was −60 mV for all other recordings. Only cells with an initial membrane potential between −50 and −70 mV and input resistances of ≥4 MΩ were used for further analysis. eEJC and PP traces were analyzed for standard parameters (amplitude, rise time, decay, charge flow, PP ratio) by using a semi-automatic custom-written MATLAB script (Mathworks; v. R2009a). Stimulation artifacts in eEJC recordings were removed for clarity. mEJC recordings were analyzed with pClamp 10 software (Molecular Devices).
Western blot analysis of larval brains
Larval central nervous system (CNS) protein extractions were performed as follows: 10 CNSs were dissected from third instar larvae. The resulting tissues were sheared manually in 5 µl of 2% SDS aqueous solution using a micropistil secured into a 1.5-ml cup. Subsequently, 0.5 µl of a 10% Triton X-100 aqueous solution and 5 µl of 2× Laemmli sample buffer were added. Samples were then heated at 95°C for 10 min. Centrifugation was performed for 2 min at 16,000 g, in order to pellet debris. For larval CNS, 10 µl (equivalent to 10 larval CNS) was subjected to denaturing SDS-PAGE using a 6% Tris-HCl gel.
Proteins of the PAGE gel were then transferred onto a nitrocellulose membrane blocked with 5% skim milk in 1× PBS supplemented with 0.1% Tween-20 and probed with rabbit anti-Rab2 (1:5,000; Santa Cruz Biotechnology; sc-285667) and mouse anti-tubulin (1:100,000; Thermo Fisher Scientific; 62204) diluted in 5% milk in 1× PBS, supplemented with 0.1% Tween-20. After washing, secondary anti-rabbit or anti-mouse IgG HRP-conjugated antibodies (Dianova) were used for detection (Dianova) in conjunction with an enhanced chemoluminescence (GE Healthcare ECL Prime; RPN 2232) detection system with Hyperfilm ECL (GE Healthcare).
Conventional embedding for transmission EM was performed as described previously (Matkovic et al., 2013). In brief, dissected third instar larvae were fixed with 4% PFA and 0.5% glutaraldehyde in 0.1 M dissolved in PBS for 10 min and subsequently in glutaraldehyde (2% glutaraldehyde in 0.1 M sodium cacodylate) for 60 min, then washed in sodium cacodylate buffer and postfixed with 1% osmium tetroxide and 0.8% KFeCn in 0.1 M sodium cacodylate buffer (1 h on ice). After washing with sodium cacodylate buffer and distilled water, the samples were stained with 1% uranyl acetate in distilled water. Samples were dehydrated, infiltrated in EPON resin, and subsequently muscles 6/7 of abdominal segment A2/3 were removed by dissection. Blocks, collected in an embedding mold, were polymerized and cut in 65–70-nm-thin serial sections. Sections were then postfixed and poststained with uranyl acetate/lead citrate. Micrographs were produced using a JEM 1011 microscope equipped with an Orius 1200A Gatan camera using the DigitalMicrograph software package (Gatan).
For quantification, the plasma membrane and the electron-dense T-bar were detected by eye and labeled manually. T-bar roof size was measured using a straight line connecting the furthest distance of the upmost T-bar dense material (in relation to the plasma membrane). T-bar area was obtained by surrounding the dense material and measuring the area of the created ROI.
For 2D morphometric analysis, images of neuronal soma from VNCs were imaged using a Zeiss 900 transmission EM. For 3D analyses, 250-nm sections were cut and collected on coated slotted grids with 10-nm gold fiducials. Once the Golgi was located, a series of images from +60° to −60° were taken with a 1° step using a Tecnai G20 microscope. Etomo/IMOD (Kremer et al., 1996; https://bio3d.colorado.edu/imod/) and Microscopy Image Browser (Belevich et al., 2016; http://mib.helsinki.fi/index.html) were used to work with 3D volumes and render 3D models of Golgi and Golgi-associated structures.
Quantification and statistical analysis
Data were analyzed using GraphPad Prism v5.01. Normality was tested with the D’Agostino and Pearson omnibus normality test. If the sample size was too low for D’Agostino and Pearson omnibus normality test (n > 7), Gaussian distribution was assumed. If data were normally distributed (or assumed to be normally distributed for n < 7) and two groups were compared, an unpaired, two-tailed t test was used (Fig. 1, B, C, E, G, I, and K; Fig. 2, C–H, K, O, Q, and T–W; Fig. 3, M and Q; Fig. 4, F, G, I–K, M, O, R, T, and W; Fig. 6 E; and Fig. 7 C). If data were normally distributed (or assumed to be normally distributed for n < 7) and more than two groups were compared, a one-way ANOVA with Tukey’s post hoc test was used (Fig. 1, M and N; and Fig. 8, B–E). If data were not normally distributed and two groups were compared, a Mann–Whitney U test was used (Fig. 2, J, L, N, P, and R; Fig. 3, C, F, and I; Fig. 4, C, E, S, U, and Z; Fig. 5, K, M, and N; Fig. 6, C, D, and I–K; and Fig. 7 H). If data were not normally distributed and more than two groups were compared, the Kruskal–Wallis test with Dunn’s post hoc test was used (Fig. 8 J). For supplemental figures, see supplemental figure legends. Significance is noted by *, P < 0.05; **, P < 0.01; ***, P < 0.001. No statistical methods were used to predetermine sample sizes, but sample sizes in these experiments were similar to those generally employed in the field. Data collection and analyses were not performed blind to the conditions of the experiments, nor was data collection randomly. For electrophysiological recordings, genotypes were measured in an alternating fashion on the same day and strictly analyzed in an unbiased manner.
Online supplemental material
Fig. S1 shows that presynaptic proteins accumulate in cell bodies of rab2−/−-deficient neurons. Fig. S2 shows presynaptic and nonsynaptic proteins in cell bodies of rab2−/−-deficient neurons. Fig. S3 shows that presynaptic Rab2 depletion mediated by RNAi knockdown caused accumulation of presynaptic proteins in neuronal cell bodies. Fig. S4 shows that presynaptic biogenesis relies on the Rab2-dependent delivery of presynaptic material. Fig. S5 shows that Rab2 is a component of anterogradely and retrogradely trafficking presynaptic precursors. Fig. S6 shows that Rab2 localizes to the Golgi in neuronal somata, and disrupted Golgi is not affecting presynaptic precursor biogenesis. Fig. S7 shows that Golgi sorting and targeting of nonsynaptic proteins are not affected, and RIM-BP accumulates at the trans-Golgi in rab2−/− mutants. Fig. S8 shows that tubule-shaped vesicular membranes accumulate at the trans-Golgi in rab2−/− mutants and shows LAMP1 and Arl8 distribution on presynaptic precursors.
We thank Marta Maglione, Niclas Gimber, Ina-Maria Schedina, Henrike Scherrer, and Luisa Fuchs for discussion and support.
T.W.B. Götz was financed by SFB 740, and SFB 1315; V. Lysiuk by SFB958; and A.G. Petzoldt by SFB958 and the Berlin Institute of Health gender equality fund.
The authors declare no competing financial interests.
Author contributions: T.W.B. Götz and A.G. Petzoldt performed immunofluorescence and genetic experiments, including analysis. T.W.B. Götz performed electrophysiological experiments. A.G. Petzoldt performed STED imaging and analysis. J. Lützkendorf performed Western blot analysis. D. Puchkov and C. Quentin performed EM experiments. D. Puchkov and A.G. Nikonenko performed EM analysis. M. Lehmann contributed to STED imaging and analysis. V. Lysiuk supported manuscript revision and performed immunofluorescence and genetic experiments. S.J. Sigrist and A.G. Petzoldt designed the study and wrote the manuscript.