Our understanding of endocytic pathway dynamics is severely restricted by the diffraction limit of light microscopy. To address this, we implemented a novel technique based on the subdiffractional tracking of internalized molecules (sdTIM). This allowed us to image anti–green fluorescent protein Atto647N-tagged nanobodies trapped in synaptic vesicles (SVs) from live hippocampal nerve terminals expressing vesicle-associated membrane protein 2 (VAMP2)–pHluorin with 36-nm localization precision. Our results showed that, once internalized, VAMP2–pHluorin/Atto647N–tagged nanobodies exhibited a markedly lower mobility than on the plasma membrane, an effect that was reversed upon restimulation in presynapses but not in neighboring axons. Using Bayesian model selection applied to hidden Markov modeling, we found that SVs oscillated between diffusive states or a combination of diffusive and transport states with opposite directionality. Importantly, SVs exhibiting diffusive motion were relatively less likely to switch to the transport motion. These results highlight the potential of the sdTIM technique to provide new insights into the dynamics of endocytic pathways in a wide variety of cellular settings.
Introduction
Although live fluorescence microscopy has played a key part in deciphering the different mechanisms underpinning the endocytic pathway, one of its limitations is the diffraction of light, which restricts imaging to cellular structures >200 nm in diameter (Li et al., 2015). A large fraction of endocytic structures are subdiffractional and, in some specialized cells such as neurons, these can constitute a considerable proportion of the cell volume. Developing adequate superresolution techniques to unravel the dynamic nature of subdiffractional endocytic structures is therefore warranted.
Synaptic vesicles (SVs) are 45-nm endocytic structures that are typically located in nerve terminals. They contain, and are able to release, neurotransmitters upon exocytic fusion with the plasma membrane at specific sites called active zones, thereby mediating fast neuronal communication (Couteaux and Pécot-Dechavassine, 1974; Südhof, 2012). Several SVs undergo fast recycling (Ryan et al., 1993; Pyle et al., 2000; Sankaranarayanan and Ryan, 2000; Gandhi and Stevens, 2003) and are capable of undergoing endocytosis, docking, and priming, thereby regaining fusion competence in a short period of time. Although the development of pH-sensitive markers has provided a wealth of information on SV recycling (Kavalali and Jorgensen, 2014), the sequence of the molecular interactions that control the recycling process has remained unclear because of the scarcity of methods that allow direct visualization of SVs in the crowded environment of the presynapse. Advances in superresolution methods have allowed the visualization and tracking of individual recycling SVs (Lemke and Klingauf, 2005; Shtrahman et al., 2005; Yeung et al., 2007; Westphal et al., 2008); however, the results obtained in these studies have relied on a limited number of trajectories from spatially isolated SVs, excluding discrete diffusional and transport states. Defining these mobility states is critical to our understanding of the intra- and intermolecular interactions that control both docking and priming, as changes in their dynamics underpin these essential processes.
In this study, we describe a method based on a pulse-chase of fluorescently tagged ligands destined to undergo endocytic transport, which we have termed subdiffractional tracking of internalized molecules (sdTIM). Using this technique, we were able to image a large number of SVs simultaneously in active zones of live hippocampal neurons with unprecedented 36-nm precision. This in turn allowed us to image the activity-dependent internalization of anti-GFP Atto647N-tagged nanobodies bound to pHluorin-tagged vesicle-associated membrane protein 2 (VAMP2) to study the mobility of recycling SVs in live hippocampal nerve endings. As a proof of concept, we demonstrated that the mobility of internalized VAMP2–pHluorin–bound anti-GFP Atto647N-tagged nanobodies was significantly lower than that of those transiting the plasma membrane (Giannone et al., 2010). We also investigated the changes in SV mobility elicited upon restimulation and showed that the mobility of SVs in presynapses, but not in the adjacent axons, increased significantly. In addition to classical mean square displacement (MSD) and diffusion coefficient characterization of SV mobility, we also accounted for nonlinear diffusion by using a combination of single-particle tracking, the moment scaling spectrum (MSS), and Bayesian model selection applied to hidden Markov modeling (HMM). By investigating these anomalous and subdiffusive events, we were able to annotate heterogeneous mobility along a single SV trajectory and discovered that, in most nerve terminals, SVs stochastically switch between either purely diffusive and/or transport mobility states. We detected similar mobility patterns upon restimulation of internalized VAMP2–pHluorin–bound anti-GFP Atto647N-tagged nanobodies. Our results highlight the power of sdTIM in studying the entire SV recycling process and the potential of the technique to be applied to investigate other subdiffractional endocytic events.
Results
Subdiffractional imaging of internalized single molecules in live hippocampal neurons
The sdTIM technique uses fluorescently labeled single molecules to track the internalization of extracellular ligands, allowing high-density single-particle tracking of relatively long trajectories. Implementation of this technique requires the use of an adjusted oblique illumination laser angle, just slightly smaller than the critical angle required for total internal reflection fluorescence microscopy to optimize synaptic access (Fig. 1 A). We first applied sdTIM to live hippocampal neurons from rats to assess the mobility of internalized vesicular proteins after their activity-dependent translocation to the plasma membrane and subsequent endocytic internalization. Like many other vesicular proteins, VAMP2 is rapidly internalized in recycling SVs (Harper et al., 2016) and has been shown to exhibit similar mobility to that of other major SV proteins, such as synaptophysin and synaptotagmin 1 (Gimber et al., 2015). Importantly, the kinetics of SV recycling are not affected by VAMP2–pHluorin overexpression in neurons, as demonstrated by comparing the kinetics of antibody-labeled endogenous SV proteins and pH-sensitive SV markers (Hua et al., 2011). To image single SVs using sdTIM (Fig. 1 B), hippocampal neurons growing on glass-bottom dishes were transfected at 14–17 d in vitro (DIV) with VAMP2–pHluorin for 24 h (Fig. 1 C) and then simulated with high K+ buffer in the presence of anti-GFP Atto-647N–tagged nanobodies (Atto647N nanobodies) for 5 min. pHluorin is a variant of GFP, the fluorescence of which is quenched in the acidic environment of the SV and becomes unquenched upon fusion with the plasma membrane (Miesenböck et al., 1998). The 5-min stimulation (pulse) provides a sufficient time frame for the uptake of Atto647N nanobodies in the recycling pool of SVs (Gaffield and Betz, 2006). The neurons were then washed extensively with low K+ buffer to remove unbound Atto647N nanobodies and incubated for 10 min in low K+ buffer (chase). To detect Atto647N nanobodies trapped in single SVs, they were then partially bleached to discriminate single SVs in the highly crowded presynaptic environment, after which we acquired time-lapse videos (50 Hz) of single-molecule fluorescence datasets of Atto647N nanobodies (Video 1), which were used to obtain high-resolution mean intensity images (Fig. 1 D), diffusion coefficients (Fig. 1 E), and single-molecule trajectory maps (Fig. 1 F). To minimize the nonspecific background, we omitted analyzing single-molecule tracks shorter than eight time points. To visualize the confined (Fig. 1 G) and mobile (Fig. 1 H) fractions of SVs, the single-particle tracking data were filtered based on the diffusion coefficients. Internalized single molecules exhibited various mobility patterns (Fig. 1 I): single (i) or dual confinement (ii) and free diffusion (iii) and directed movement (iv). SVs were also observed to move in groups where multiple vesicles visited a confinement site (Fig. 1 I, v) or moved in a circular directed pattern at the perimeter of an active nerve terminal (Fig. 1 I, vi).
Mobility of VAMP2–pHluorin–bound Atto647N nanobodies decreases when internalized in SVs
To validate the sdTIM technique, we compared the mobility of VAMP2–pHluorin/Atto647N nanobodies on the plasma membrane with that of internalized nanobodies. For this, we used a technique called Universal Point Accumulation Imaging in Nanoscale Topography (uPAINT; Giannone et al., 2010, 2013), which is ideal to establish the mobility of plasma membrane proteins (Fig. S1 A) and has previously been used to track single glutamate receptors using anti-GFP nanobodies (Nair et al., 2013). Hippocampal neurons were transfected with VAMP2–pHluorin (Fig. S1 B) for 24 h and imaged at 50 Hz during high K+ stimulation in the presence of Atto647N nanobodies. After single-particle tracking (Nair et al., 2013), we generated superresolved mean intensity images (Fig. S1 C), diffusion coefficients (Fig. S1 D), and single-molecule trajectory maps (Fig. S1 E) of VAMP2–pHluorin/Atto647N nanobodies on the plasma membrane.
We then compared the mobility of VAMP2 on the plasma membrane (uPAINT) with that of internalized VAMP2 (sdTIM) by calculating the MSD (µm2) using the PALMtracer plug-in (Kechkar et al., 2013; Nair et al., 2013) for MetaMorph. As expected, we observed a significant reduction in the mobility of internalized VAMP2 compared with that on the plasma membrane (Fig. 2 A). The area under the MSD curve (AUC) for each cell in each condition also showed a significant decrease in mobility upon VAMP2 internalization (Fig. 2 B). We next analyzed the diffusion coefficient (D) of the single-molecule trajectories. Based on our observations, the frequency distribution of the Log10D values of VAMP2–pHluorin/Atto647N nanobodies on the plasma membrane and internalized in SVs was bimodal. The threshold for the immobile and mobile fractions (indicated by the dashed line in Fig. 2 C) was calculated as described earlier (Constals et al., 2015): the displacement threshold was 0.03 µm2 s−1 (Log10D = −1.45, if [D] = µm2 s−1). The immobile fraction was composed of VAMP2–pHluorin/Atto647N nanobodies for which displacement within four frames was below the spatial detection limit of our method (106 nm). Internalized VAMP2–pHluorin/Atto647N nanobodies displayed a significantly higher frequency of immobile molecules (Log10D ≤ −1.45) and a lower frequency of mobile molecules (Log10D > −1.45) than those detected on the plasma membrane (Fig. 2 C and Fig. S2 A). Correspondingly, the mobile-to-immobile ratio of VAMP2–pHluorin/Atto647N nanobodies was significantly higher for molecules on the plasma membrane than for internalized molecules (Fig. 2 C, inset). The mobile-to-immobile ratio of internalized VAMP2–pHluorin/Atto647N nanobodies was similar to that reported earlier for other major SV proteins, synaptophysin– and synaptotagmin–pHluorin (Gimber et al., 2015), as well as for postsynaptic AMPA receptors (Groc et al., 2007), validating the use of VAMP2–pHluorin as an SV marker. The average diffusion coefficient for the mobile fraction of VAMP2–pHluorin on the plasma membrane was 0.14 µm2 s−1, which is similar to that reported earlier for VAMP2–pHluorin (Gimber et al., 2015). The average diffusion coefficient for the immobile fraction of internalized VAMP2–pHluorin was 0.009 µm2 s−1, which is similar to the mobility reported in a recent study on SV mobility in pools (Rothman et al., 2016). The localization precision (σ; Michalet, 2010) was 36.0 ± 1.6 nm for uPAINT and 36.2 ± 1.3 nm for sdTIM, similar to the value reported previously for SV tracking (Westphal et al., 2008).
To further characterize the mobility of SVs and to account for nonlinear diffusion, we analyzed their MSS and compared the results of sdTIM experiments to uPAINT controls. Previously recorded movies were analyzed using the utrack program (Jaqaman et al., 2008). In contrast to MSD measurements, MSS (Ferrari et al., 2001; Ewers et al., 2005) allows examination of cases other than the second moment of moving particles (i.e., anomalous and subdiffusive events). The utrack program also allowed us to recover diffusion coefficients for each detected trajectory and plot the individual MSS slopes versus their respective diffusion coefficients. The distributions of the diffusion coefficient versus MSS slope values were significantly different for the sdTIM and the uPAINT experiments (Fig. S1 F). We also fitted the distributions to 2D Gaussian functions that further highlighted the different states and the shift toward smaller diffusion coefficient/MSS slope pairs for the sdTIM experiments (Fig. S1, G and H). These results demonstrate the presence of a distinct subset of VAMP2–pHluorin/Atto647N nanobody mobility in the sdTIM experiments, further validating that these trajectories are derived from Atto647N nanobodies internalized in SVs.
sdTIM enables simultaneous detection of multiple SVs in live hippocampal nerve terminals
Despite exhibiting some variability (Mutch et al., 2011), a mean of 70 copies of VAMP2 can be found in SVs (Wilhelm et al., 2014). Given this large number of VAMP2 copies per SV, there was a high probability of multiple nanobodies being trapped in a single SV, even when used at very low concentrations. The stochastic number of nanobodies per vesicle and the random number of Atto647N molecules per nanobody made it challenging to determine the number of vesicles for which movement was visualized in a crowded presynaptic terminal. To overcome this issue, we used single-molecule detection and stepwise photobleaching to calculate the number of conjugated Atto647N molecules per nanobody and the number of nanobodies in the presynaptic terminal. A low concentration of Atto647N nanobodies was deposited on a clean glass surface to minimize the chance of detecting two or more nanobodies within the same diffraction-limited spot. We then monitored the stepwise emission of fluorescence to directly assess the number of Atto647N molecules per nanobody. At a neutral pH, 28% of emission traces had one step, 36% showed two steps, and 28% had three steps (Fig. 3, A and B). A small number of traces (8%) had more than three steps, which is consistent with a random distribution of the protein on the glass surface. When we monitored the fluorescence emission of Atto647N nanobodies at pH 5.0 to mimic the low intravesicular pH (Mel’nik et al., 1985; Fig. 3 C), no change in the average number of Atto647N emission steps per nanobody was detected (Fig. 3 D). However, we noticed a slight increase in the signal to noise ratio of individual emission steps in the intensity-time traces. The fluorescence stability of Atto647N allowed us to reliably count more than five emission steps (example trace in Fig. S2 B).
In our sdTIM experiments, hippocampal neurons were exposed to partial bleaching so that single moving vesicles could be followed. The number of bleaching steps helped us to determine the number of vesicles in the nerve terminals. This was achieved by imaging fixed neurons after uptake of Atto647N nanobodies (sdTIM). The bleaching time, before imaging, was adjusted to maximize the detection of single step emission from single SVs in presynapses. This ensured that most of the fluorophores were bleached and that the large majority of SVs contained only one nanobody, with only one or two emitting Atto647N fluorophores. Representative fluorescence-time traces obtained from synaptic regions in fixed hippocampal neurons are shown in Fig. 3 E, together with the distribution of Atto647N nanobody emission steps after the initial bleaching (Fig. 3 F). The average number of emission steps originating from a single SV was 1.5.
In live hippocampal neurons, internalized Atto647N nanobodies could easily be traced after a similar initial bleaching step before imaging in sdTIM. We detected single SVs entering the presynaptic nerve terminals, followed by their confined motion within this region and exit into the perisynaptic region. To determine the number of SVs that resided within the active zone at any given time, two methods were used. First, to obtain intensity-time traces, we recorded the fluorescence intensity of each detected spot within the single particle trajectories. The fluorescence of the trajectories did not change as long as only one vesicle was followed in the presynaptic region. When a second vesicle entered the region and was found within the same diffraction-limited spot, a steplike increase in fluorescence was observed. We counted the total number of emission steps from all trajectories recovered from one nerve terminal to estimate the maximum number of vesicles for which positions could be detected within the same diffraction-limited spot. The average number of steps was 5.4 ± 1.3 (SD). In the second approach, a region of interest was chosen to cover the whole presynaptic terminal. The fluorescence intensity over time exhibited similar step-like emission to that in fixed cells, with the number of steps reflecting the number of nanobodies within that region (Fig. 3 G). We found a mean of 6.0 ± 1.6 (SD) steps, which was in good agreement with the results from the first method. Given that after bleaching there was on average, 1.5 fluorescently active Atto647N molecules per nanobody and estimating one nanobody per vesicle, we calculated that the number of SVs we could detect within each presynaptic active zone, at any given time, was 4 ± 1 (Fig. 3 H). These results highlight the potential of simultaneous detection of several SVs in the crowded environment of the hippocampal presynapse.
SV movement in active presynapses is highly confined
The assessment of SV mobility in nerve terminals has previously been based on morphological features of the neuron (Rizzoli and Betz, 2004). To unequivocally discriminate active nerve terminals from surrounding axonal segments, VAMP2–pHluorin fluorescence (Fig. 4, A and B) was imaged during high K+ stimulation and concomitant uptake of Atto647N nanobodies. As expected, after SV fusion and VAMP2–pHluorin unquenching, the fluorescence intensity increased significantly in presynaptic sites compared with axons (Fig. 4, C and D), as previously established (Gimber et al., 2015). Nerve terminals (Fig. 4, E and F) and intersynaptic axonal segments (Fig. 4, E and I) were then partially bleached (Fig. 4, G and J) and imaged at 50 Hz to detect internalized Atto647N nanobodies. Single-molecule tracking of Atto647N nanobodies trapped in SVs in active presynapses was achieved by delineating regions that exhibited activity-dependent fluorescence unquenching. The mobility of internalized Atto647N nanobodies in presynapses was mostly confined (Fig. 4 H), with more mobile tracks being found in axonal segments (Fig. 4 K). As previously established, we could detect SVs moving in and out of active presynapses (Fig. 4 H, dark blue track); this most likely reflects intersynaptic vesicle exchange (Kamin et al., 2010), a process that has been proposed to be involved in synaptic homeostasis and plasticity in vivo (Herzog et al., 2011).
SVs were significantly less mobile in active presynapses than in axons, as demonstrated by comparing the MSD curves and AUC of VAMP2–pHluorin/Atto647N nanobodies in the respective regions (Fig. 4, L and M, respectively). Accordingly, the Log10D showed a significantly larger immobile fraction of SVs in presynapses compared with axons (Fig. 4 N and Fig. S2 C). Consequently, we detected a significantly smaller mobile-to-immobile ratio in presynapses than in axons (Fig. 4 O). These results are in good agreement with previous studies performed by SV immunolabeling (Kamin et al., 2010) using QuantumDot-based presynaptic probes (Lee et al., 2012) or surface-resident VAMP2–pHluorin (Sankaranarayanan and Ryan, 2000). Our results confirm that intersynaptic SV mobility (which occurs in axons) is greater than in presynapses. To analyze the mobility of internalized VAMP2–pHluorin/Atto647N nanobodies in live hippocampal neurons (Fig. S3 A) in more detail, we calculated the vector displacement for each SV trajectory and the cumulative magnitude of these vectors in 100-nm2 grids to obtain a quiver plot of the detected single molecules (Fig. S3 B). The resulting vector map clearly exhibited areas of vector convergence and lower displacement consistent with the expected SV confinement in presynapses (Fig. S3 C) as compared with axonal regions (Fig. S3 D).
To confirm the correct targeting of the Atto647N nanobodies in SVs, hippocampal neurons were subjected to sdTIM, fixed, and immunostained against another SV protein, Synapsin-1. As expected, the unquenched signal of VAMP2–pHluorin in active presynapses after stimulation colocalized with labeling against endogenous Synapsin-1 (Fig. 5 A), confirming the identity of the active presynapses. Importantly, we only detected Atto647N nanobodies internalized in neurons expressing VAMP2–pHluorin (Fig. 5, B and C), confirming the anti-GFP specificity of the nanobodies (Kubala et al., 2010). To assess the correct localization of the nanobodies in SVs, we performed correlative light microscopy and EM. Hippocampal neurons were grown on gridded cell culture dishes and transfected with VAMP2–pHluorin for 24 h, after which the localization of the neurons expressing VAMP2–pHluorin on the grids was recorded (Fig. 5 D). Cells were then processed for sdTIM except for the use of anti-GFP HRP-mCherry-GFP nanotrap (GNT) nanobodies and omitting the live-cell imaging of the SVs. After internalization, the neurons were fixed, subjected to cytochemical staining, and processed for EM. The VAMP2–pHluorin–positive neurons were then located based on their recorded position on the grids. This allowed us to identify the precise localization of the nanobodies in these neurons based on the electron-dense precipitate produced by the HRP substrate (Fig. 5 E). Although some of the HRP precipitate was detected in endocytic structures >45 nm in diameter, the great majority (88.5 ± 3.0%) was found in SVs (Fig. 5 F). We also observed some residual staining on the plasma membrane, indicating that not all nanobodies were internalized. Importantly, the neurons not expressing VAMP2–pHluorin did not contain internalized HRP precipitate. These results provide validation of the specific internalization of the GFP-bound nanobodies into SVs at VAMP2–pHluorin–positive nerve terminals and further validate the use of sdTIM to track SVs in active presynapses.
Restimulation increases presynaptic SV mobility
Neurotransmission is mediated by SV exocytosis at presynapses. Although considerable efforts have been dedicated to studying the events leading to SV exocytosis (Betz et al., 1992; Lemke and Klingauf, 2005; Kamin et al., 2010), the details of how SVs move toward the plasma membrane to undergo fusion (Kamin et al., 2010) remain unclear.
In agreement with a previous study demonstrating that recently endocytosed SVs are more mobile and become gradually less mobile when integrated into vesicle clusters (Kamin et al., 2010), our MSD measurements of the recently endocytosed SVs (t0–t120s) compared with resting SVs (t0–t200s) indicated a significant decrease in mobility over time (Fig. S4, A–C). To study the dynamic steps leading to exocytosis and neurotransmitter release, we introduced a second round of high K+ stimulation (restimulation) of live VAMP2–pHluorin–expressing hippocampal neurons after the initial uptake of Atto647N nanobodies and compared the mobility with the corresponding time frame of the sdTIM (t0–t120s). Restimulation with high K+ after the internalization of Atto647N nanobodies led to the unquenching of VAMP2–pHluorin and an increase in fluorescence similar to that detected during the first round of stimulation (Fig. 6, A–C), thereby confirming the localization of functional release sites. Restimulation increased the SV mobility significantly when compared with unstimulated conditions in both whole neurons (Fig. 6, D–F; and Fig. S2 D) and presynaptic regions of interest (Fig. 6, G–I; and Fig. S2 E). Interestingly, the restimulation had no effect on the mobility of intersynaptic SVs (Fig. 6, J–L) imaged in axonal segments. These results are different from those of a previous study that showed that stimulation had no effect on the mobility of recently endocytosed SVs (Kamin et al., 2010). However, the authors could only assess SV mobility minutes after the stimulus, suggesting that they may have missed effects at the onset or offset of stimulation. The discrepancy in results might therefore be explained by our ability to follow the restimulation effects on SV mobility immediately after introduction of the second stimulus with high K+ (within 10 s).
When comparing the effects of restimulation on SV mobility between axons and presynapses, we observed similar, although less profound, differences in the MSD and Log10D measurements (Fig. 7, A–C; and Fig. S2 F) than in our previous sdTIM results (Fig. 4, L–O), with SVs being more mobile in axons than in presynapses. Based on our earlier observations, the mobility of SVs increased upon restimulation in presynapses (Fig. 6, G–I; and Fig. S2 E) but remained unaltered in axons (Fig. 6, J–L). Interestingly, although the SV mobility increased in presynapses, it did not rise to the level observed in axons (Fig. 7, A–C; and Fig. S2 F). Together, these results indicate that the mechanisms responsible for the SV mobility detected in axons are most likely different from those in presynapses.
Interestingly, our data demonstrate that the increase in SV mobility elicited by restimulation is indistinguishable from that of VAMP2–pHluorin–bound Atto647N nanobodies on the plasma membrane (Fig. 7, D–F). These results indicate that, upon restimulation, SVs relocate to the plasma membrane and undergo exocytosis, after which the SV contents, including VAMP2–pHluorin–bound Atto647N nanobodies, are released to move on the plasma membrane. These results are consistent with those of a previous study showing that the mobility of vesicular proteins increases after their translocation to the plasma membrane (Kamin et al., 2010). Together, our results suggest that the VAMP2–pHluorin/Atto647N nanobody mobility detected after restimulation contained information about the SV mobility inside presynapses and the dynamic events of exocytosis. Dissecting the exact time points of these events will provide a wealth of information regarding these crucial events and should prove an interesting subject for future studies. Furthermore, these results indicate that the sdTIM technique can be used to study the mobility of SVs after restimulation and, upon fusion, to investigate the dynamics of the VAMP2–pHluorin/Atto647N nanobodies on the plasma membrane. In contrast, although the mobility of single molecules on the plasma membrane can be addressed using styryl dyes, the dye repartitions into the extracellular solution shortly after exocytosis (Zenisek et al., 2002; Richards et al., 2005).
SVs stochastically switch between distinct motion states in presynapses and axons
Single-molecule trajectories are classically characterized using MSD and diffusion coefficient analyses (Gimber et al., 2015; Maschi and Klyachko, 2015; Wang et al., 2015). The main disadvantage of these approaches is the time averaging, which results in loss of the alternating motion states along a single trajectory. HMMs have been successfully used to infer the number of motion states and the state transition probabilities from experimental data and to annotate each step of the trajectories with the respective motion state (Monnier et al., 2015). By using Bayesian model selection in the inference process, the simplest mobility model describing a set of trajectories can then be selected in an objective manner (Persson et al., 2013).
We therefore applied Bayesian model selection to HMMs to investigate the heterogeneous mobility patterns of internalized VAMP2–pHluorin/Atto647N nanobodies in live hippocampal neurons (see Materials and methods). It is common practice to analyze single-particle trajectories with a fixed number of hidden states (Das et al., 2009, 2015; Chung et al., 2010; Thurner et al., 2014; Freeman et al., 2015; Katz et al., 2016). When we assumed a maximum of four hidden states, we found that the optimal number of states was three for almost all presynapses (86%). Thus, for simplicity, we therefore assumed a maximum of three hidden states in our analysis. We observed heterogeneous transport dynamics across presynapses and axons. In 30 out of 31 presynapses and 28 out of 29 axonal segments, we detected three motion states, with the rest displaying two states. In both presynapses and axons, the mobility states displaying transport motion (65 and 93%, respectively) were more common than purely diffusive states (35 and 7%, respectively; Fig. 8 A). In 26% of the presynapses and 83% of the axonal segments, one state exhibited pure diffusive motion (D) with a low apparent diffusion coefficient (0.06 ± 0.01 and 0.06 ± 0.02 µm2 s−1 on average in presynapses and axons, respectively) representing the immobile subpopulation, and the other two were active transport states (DV).
The angle between the velocity vectors of the two transport states in these D-DV-DV mobility states was close to 180° (173.7 ± 4.1°), indicating bidirectional movement. Interestingly, we found that the SVs underwent stochastic switching between these states (Fig. 8 B) and that the transition probabilities between the states differed (Fig. 8 C). In 32% of the presynapses, all three states were purely diffusive, with average apparent diffusion coefficients of 0.03 ± 0.01, 0.11 ± 0.01, and 0.32 ± 0.01 µm2 s−1 (Fig. 8 D). The values for the other motion model types are given in Table S1.
In addition, our analysis provided population distributions of SV motion parameters in presynapses and axons (Fig. 8, E–G). In agreement with our earlier results (Fig. 4 L), we found that the average velocity of transport states (Fig. 8 F) was significantly higher (P < 0.0001) in axons (5.1 µm s−1) than in presynapses (1.9 µm s−1). These results are in agreement with a previous study (Shtrahman et al., 2005) showing that recycling SVs switch between bound and unbound states and spend more time in the immobile state. Furthermore, based on our results, the probability of remaining in the same state was higher for the diffusive state than the transport state in both presynapses and axons (Fig. 8 G), indicating that once an SV entered a diffusive state, it was relatively more likely to stay in that state than change to a transport state.
We next applied Bayesian model selection to HMMs for the restimulation data collected earlier. In 24 out of 34 presynapses and 32 out of 34 axonal segments, we detected three motion states, with the rest displaying two states. In both presynapses and axons, the mobility states displaying transport (74 and 97%, respectively) were more common than purely diffusive states (26 and 3%, respectively) (Fig. 9 A). The values for the other motion model are given in Table S1. The population distributions of SV motion parameters in presynapses and axons are shown in Fig. 9 (B–D). The average velocity magnitude of SVs was significantly higher (P < 0.0001) in axons than in presynapses (5.0 and 3.3 µm s−1, respectively). These results reveal that synaptic vesicles are actively transported within presynapses rather than relying on passive diffusion alone, as suggested by previous studies (Holt et al., 2004; Rea et al., 2004; Tokuoka and Goda, 2006; LoGiudice et al., 2008; Rothman et al., 2016). HMM-Bayes model selection applied to SVs in presynapses and axons demonstrated both a heterogeneous combination of diffusive and transport states and stochastic switching between these states.
Discussion
Single-particle tracking allows the characterization of the mobility of a single protein and can provide unprecedented insights into the motion heterogeneities associated with crucial cellular events (Freeman et al., 2016). Most insights regarding SV mobility have been gained by indirect methods, such as fluorescence recovery after photobleaching (Henkel and Betz, 1996; Kraszewski et al., 1996; Gaffield et al., 2006; Rothman et al., 2016) and correlation spectroscopy (Jordan et al., 2005; Shtrahman et al., 2005; Yeung et al., 2007). Other approaches have used a variety of styryl dyes to study endosomal trafficking in general (Lemke and Klingauf, 2005; Gaffield and Betz, 2006; Peng et al., 2012) or viral delivery (Gandhi and Stevens, 2003) to accomplish sparse labeling of SVs and the subsequent tracking of a relatively low number of SVs per presynapse. A hallmark study in the field provided high-resolution imaging of isolated vesicles in live hippocampal terminals and throughout the entire axons using stimulated emission depletion microscopy (Westphal et al., 2008).
Commonly, only two diffusive states of SV mobility (bound and free) are assumed, and single-particle tracking data have been fitted based on this assumption (Bakshi et al., 2012). Most current single-particle tracking techniques, such as the ones using photoconvertible fusion proteins, generate large numbers of relatively short trajectories, demonstrating the need for methods that can track and analyze a large number of trajectories simultaneously over longer periods (Manley et al., 2008; Persson et al., 2013). Several other techniques have previously been used to characterize SV mobility (Lemke and Klingauf, 2005; Westphal et al., 2008; Hoopmann et al., 2010; Kamin et al., 2010; Hua et al., 2011; Gimber et al., 2015). The main advantage of the sdTIM technique compared with these methods is that it takes advantage of several technological advances, such as the use of smaller nanobodies conjugated with stable Atto647-fluorophores with increased fluorescence lifetime, which are ideal for superresolution microscopy. This allowed us to follow several SVs simultaneously for longer timeframes. Furthermore, sdTIM is a relatively fast and simple technique with a high spatiotemporal resolution (36-nm localization precision and 20-ms imaging rate), and is capable of generating high-density output data. Monitoring pHluorin fluorescence during nanobody loading or unloading has allowed us to be highly confident that the trajectories we retrieved stemmed from either nerve terminals or axonal segments. The resulting high-density output data enabled us, in turn, not only to determine the average diffusion but also to isolate the heterogeneous stages of the individual SVs, providing novel insights into SV mobility in live hippocampal neurons. One of the major advantages of sdTIM is that it is also likely to prove valuable in the context of genetic or pharmacological manipulations.
The internalization of Atto647N nanobodies was only detected in VAMP2–pHluorin–positive neurons, indicating a selective endocytic uptake in SVs, a result that was confirmed by correlative light microscopy and EM. It is worth noting, however, that upon pulse-chase uptake of VAMP2–pHluorin–bound Atto647N nabobodies, we also detected a fraction of molecules that retained the same mobility as that of molecules detected on the plasma membrane with uPAINT (Gordon et al., 2011), as well as a residual HRP precipitate on the plasma membrane and in endosomes in the electron micrographs. Although the diffusion coefficient threshold between mobile and immobile fractions has been previously used as a criterion to identify the plasma membrane-resident population from internalized SVs (Lee et al., 2012), our results reveal a remarkable heterogeneity in SV mobility, even along a single trajectory. We therefore cannot exclude the possibility that a small fraction of our single-molecule detections originate from structures other than SVs. Although providing a powerful tool to study SV mobility, sdTIM uses oblique illumination and is therefore limited to imaging a thin optical section. Furthermore, our tracking is performed in 2D, and further advances in the field such as the use of astigmatic lenses will be required to assess the mobility of single molecules in 3D. Moreover, although sdTIM can be used to study the mobility steps of SV recycling, it cannot be used to detect all pools of SVs at the simultaneously. The use of dual fluorescence imaging should address this limitation.
Tracking single vesicles revealed mobile and immobile fractions that exhibited various mobility patterns and confinement. Directed movement at the perimeter of presynapses suggests that cytoskeletal elements could be involved in driving active motion. Another nonexclusive possibility is that the directed motion of a vesicle or other organelle could generate nanofluidic effects capable of mobilizing a trail of freely moving SVs before their recruitment. Such behavior is compatible with previous findings indicating that decreased mobility of SVs results from their incorporation into vesicle clusters and docking sites (Kamin et al., 2010) and that vesicle collisions and hydrodynamic interactions are among the major determinants of SV mobility in crowded presynaptic terminals (Rothman et al., 2016). Our data showed that the mobility of internalized VAMP2–pHluorin/Atto647N nanobodies was significantly lower in active nerve terminals than on the plasma membrane. Furthermore, our results demonstrating that the internalized VAMP2–pHluorin/Atto647N nanobody mobility decreased gradually and was significantly lower than that found in axonal segments are in good agreement with previous studies (Westphal et al., 2008; Kamin et al., 2010). Restimulation reversed these effects in the presynapse, but not in the axonal segments. Ultimately, the SVs resumed mobility that was indistinguishable from that of VAMP2–pHluorin/Atto647N nanobodies on the plasma membrane, indicating that upon SV exocytosis, these nanobodies could be detected transiting on the plasma membrane. Dissecting the precise temporal nature of these highly dynamic SV recycling events will allow further investigations of the key steps of SV mobility during neurotransmission. Although these questions remain to be addressed in future studies, our results validate the power of the sdTIM technique for the study of the entire SV cycle.
The HMM-Bayes analysis allowed us to unveil a range of apparent diffusion coefficients for SVs in presynapses. When averaged, the values obtained were consistent with those of some earlier studies (Holt et al., 2004; Rea et al., 2004) but differed from those of other studies (Jordan et al., 2005; Lemke and Klingauf, 2005; Shtrahman et al., 2005; Gaffield and Betz, 2006). One of the features of sdTIM is that this method is capable for producing hundreds of trajectories for each presynapse and thousands for each neuron. This allowed us to identify a wide range of apparent diffusion coefficients without the caveat inherent in time averaging, which leads to the loss of hidden mobility states and state transitions. The commonly presented average values therefore do not accurately represent the fluctuations in SV mobility. Furthermore, the use of different markers, cell types, methods, and time scales in the experiments might also lead to differing results. Finally, it is important to stress that the values reported in our study were not corrected for localization precision and motion blur and therefore only present apparent diffusion coefficients.
By applying the HMM-Bayes analysis to single-molecule data, we have provided a comprehensive analysis of the heterogeneous motion states of SVs in presynapses and axons in resting and stimulated neurons. In contrast to the current view, our results revealed a stochastic switching between the diffusive and transport motion of SVs. Our analysis also showed that SVs are relatively less likely to switch from diffusive states to transport states in either resting or stimulated neurons. It is noteworthy that the main conclusion from this analysis, that SVs stochastically switch between different motion states, does not depend on the maximum number of hidden states allowed. The sdTIM technique described in this study can potentially be applied to study other subdiffractional endocytic structures, enabling measurements of changes in SV mobility in response to genetic or pharmacological manipulations. Although further experiments will be required to address the mechanisms underlying these diverse diffusive states and transport dynamics, our study provides the first quantitative account of the heterogeneous nature of SV mobility.
Materials and methods
Neuronal cultures
Hippocampal neurons were obtained from Sprague Dawley rats at embryonic day 18 as previously described (Kaech and Banker, 2006; Heine et al., 2008; Harper et al., 2011) and cultured on poly-d-lysine–coated glass-bottom dishes (D29-20-1.5N; Cellvis) at 50,000–100,000 cell densities. Experiments were conducted with the approval of The University of Queensland Animal Ethics Committee. Neurons were transfected with Lipofectamine 2000 (11668–019; Invitrogen) at 14–17 DIV, according to the manufacturer’s protocol, and imaged at 15–18 DIV.
Constructs and antibodies
VAMP2–pHluorin was provided by J. Rothman (Yale University, New Haven, CT; Miesenböck et al., 1998). The primary Synapsin-1 antibody (106011; Synaptic Systems) and the secondary antibody Alexa Fluor 564 (A11030; Invitrogen) were obtained from the respective manufacturers. The anti-GFP nanobody (gift from K. Alexandrov, Institute for Molecular Bioscience, The University of Queensland, Brisbane, Australia) was labeled with Atto647N according to the manufacturer’s protocol (ATTO-TEC GmbH) as previously described (Harper et al., 2016). Nanobodies were diluted to 1 mg ml−1 and dialyzed against 20:1 PBS, pH 7.4/0.2 M NaHCO3, pH 9.0. A two- to threefold molar excess of Atto-647N normal human serum–ester, which was reconstituted in DMSO, was then added to the solution of nanobodies and incubated at room temperature for 1 h. Excess dye was removed using a Sephadex gel filtration column (17-0851-01; GE Healthcare), and a Bradford protein assay (B6916; Sigma-Aldrich) was performed to determine the concentration of the labeled protein. HRP-mCherry-GNT was constructed as follows: a prokaryotic mCherry-GNT expression vector yielding N-terminal fusion of the 6His-Ubiquitin tag (Catanzariti et al., 2004) was constructed using circular polymerase extension cloning (Quan and Tian, 2011). In brief, mCherry (Takara Bio Inc.) and GNT (Kirchhofer et al., 2010) were inserted into the pHUE backbone, using a set of appropriate primers (pHUE forward: 5′-AGCTTAGATCCGGCTGCTAACAAAGCCCG-3′; and pHUE reverse: 5′-TCCACCGCGGAGGCGCAACACCAGGTGCAG-3′; mCherry forward: 5′-TGTTGCGCCTCCGCGGTGGATGCAGCAGCGGTGTGAGCAAGGGCGA-3′; and mCherry reverse: 5′-CCGCTTTCCACCAGTTGAACCTGCATCTTGTACAGCTCGTCCATGCCGCCGGT-3′; and GNT forward: 5′-ATGCAGGTTCAACTGGTGGAAAGCGG-3′; and GNT reverse: 5′-CGGGCTTTGTTAGCAGCCGGATCTAAGCTTCAAGAGCTCACCGTCACCTGAGTCCCCTG-3′). Sequences encoding a cysteine residue (indicated in italics) were introduced between the ubiquitin tag and mCherry. The engineered cysteine was exposed to the N-terminal of the mCherry-GNT fragment because of proteolytic cleavage using deubiquitylating enzymes (Catanzariti et al., 2004). The cysteine residue appended to an N-terminal extension of the mCherry-GNT protein was reduced to free sulfhydryl with mercaptoethylamine hydrochloride (2-MEA; Thermo Fisher Scientific) at a concentration of 20 mM in PBS at room temperature for 1 h. Reduced protein was separated from excess reducing agent by gel filtration over a desalting gel (PD-10; GE Healthcare). The protein was eluted with ice-cold PBS, concentrated to 5 mg ml−1, and then applied for labeling. Peroxidase maleimide activated (P1709; Sigma-Aldrich) was reconstituted to 5 mg in ice-cold PBS and then combined with mCherry-GNT in a 1:2 ratio. The reaction mixture was incubated at 4°C with gentle agitation overnight. After reaction, the unreacted free HRP and mCherry-GNT were separated using gel filtration chromatography (Superdex-75; GE Healthcare).
Reagents
Hepes (H3375), ascorbic acid (A5960), CaCl2 (C5080), and BSA (A8022) were purchased from Sigma-Aldrich. NaCl2 (X190) and d-glucose (0188) were purchased from AMRESCO. KCl (1206119; Ajax Finechem Pty Limited) and MgCl2 (MA029; Chem-Supply) were purchased from the respective manufacturers.
sdTIM
Cultured hippocampal neurons expressing VAMP2–pHluorin were washed once with a low K+ buffer (0.5 mM MgCl2, 2.2 mM CaCl2, 5.6 mM KCl, 145 mM NaCl, 5.6 mM d-glucose, 0.5 mM ascorbic acid, 0.1% BSA, and 15 mM Hepes, pH 7.4), after which the neurons were incubated with high K+ buffer (as in low K+ with the exception of 95 mM NaCl and 56 mM KCl) containing 3.19 pg µl−1 Atto647N-tagged nanobodies for 5 min at 37°C. After stimulation, unbound Atto647N nanobodies were washed off with low K+ buffer, and neurons were then chased for 10 min at 37°C to maximize the internalization of VAMP2–pHluorin–bound Atto647N-tagged nanobodies. To detect single molecules of the internalized nanobodies, the nerve terminals were exposed to high laser illumination to partially bleach the majority of the Atto647N fluorophores, and the remaining internalized Atto647N nanobodies were imaged using oblique illumination microscopy. Localization precision, σ (36.2 ± 1.3 nm), was calculated according to ε = 4σ2 (Michalet, 2010), where ε is the y-intercept of the average of the MSD values.
uPAINT
Experiments were performed as per Giannone et al. (2010). Hippocampal neurons from rats were transfected with VAMP2–pHluorin. To track VAMP2–pHluorin, we used anti-GFP Atto647N nanobodies (Kubala et al., 2010) at 3.19 pg µl−1 concentration. σ (36.0 ± 1.6 nm) was calculated as described in the sdTIM section.
Oblique illumination microscopy
For live-cell oblique illumination, neurons were visualized at 50-Hz (10,000–20,000 frames by image streaming) and 20-ms exposure, at 37°C on a microscope equipped with an ILas2 double-laser illuminator (Roper Technologies), a CFI Apo TIRF 100× 1.49 NA objective (Nikon), and an Evolve512 delta EMCCD camera (Photometrics). Image acquisition was performed using Metamorph software (MetaMorph Microscopy Automation and Image Analysis Software, version 7.7.8; Molecular Devices). A quadruple beam splitter (LF 405/488/561/635-A-000-ZHE; Semrock) and a QUAD band emitter (FF01-446/510/581/703-25; Semrock) were used.
Single-molecule detection and tracking
The single-molecule localization and dynamics were extracted from the 6,000–10,000 frame oblique illumination acquisitions as described earlier (Nair et al., 2013). The VAMP2–pHluorin–bound Atto647N-tagged nanobodies were detected and tracked using a combination of wavelet segmentation (Izeddin et al., 2012) and optimization of multi-frame object correspondence by simulated annealing (Racine et al., 2006). PALM-Tracer, a custom-written program that runs as a plug-in within the Metamorph software, was used for molecule localization and tracking of the VAMP2 dynamics (MSD and diffusion coefficient, D). The Log10D distribution was sorted into immobile and mobile fractions as previously described (Constals et al., 2015). We only included tracks longer than eight frames in the analysis to minimize nonspecific background. The color coding for the superresolved images was performed using ImageJ (2.0.0-rc-43/1.50e; National Institutes of Health). Diffusion coefficients were calculated for each trajectory and presented in a color-coded pixel at the site of localization. In the mean intensity maps, each pixel indicates an individual molecule. The area with the highest density is represented in black, whereas white represents regions with no detection. The color-coding of the trajectory maps represents the detection time point along the acquisition: cold colors indicate early detections.
Nanobody fluorescence intensity fluctuation analysis
Single Atto647N-conjugated nanobodies were imaged in a 22 × 10 × 1-mm chamber, the bottom of which was made of microscope coverglass (Menzel Glaser, 1.5; Thermo Fisher Scientific). Before assembling the chamber, the cover glass was cleaned by sonication for 20 min in 1M KOH, followed by extensive rinsing in Milli-Q water (EMD Millipore) and drying by airflow. 200-µl of nanobodies diluted in PBS (1:104) were incubated in the chamber for 15 min before washing with PBS only. For experiments in an acidic environment, PBS was titrated with hydrochloric acid until it reached pH 5.0 and then used immediately. The density of nanobodies was 50–100 per 27 × 27 µm2. For the characterization of nanobodies, 5 × 5 pixels around the center of each fluorescent spot were chosen to obtain intensity-time traces using a custom-written routine described previously (Durisic et al., 2014). The background was subtracted locally and calculated from the mean intensity of the region two pixels larger than the edge of the Atto647N nanobody signal. Partially overlapping nanobodies and fluorophores that photobleached in less than two frames were excluded from the analysis. Emission steps were counted manually and defined as events in which the mean intensity of the smallest step was at least twice the SD of the background fluorescence. The steps within each single trace were similar in amplitude. Missed events were found in <2% of traces, and an appropriate correction was made. It was imperative that the signal-to-noise ratio was at least 1.5 so that the number of counted steps accurately reflected the number of fluorophores. Intensity-time traces from synaptic regions were extracted and characterized in the same way, except that the region of interest was adjusted to cover the whole synapse. The average number of steps in presynaptic vesicles was 6.0 (±1.6 SD). Given that the ratio of Atto647N fluorophore dyes per nanobody was 1.5:1 and assuming one Atto647N nanobody per vesicle, we estimate that we were able to detect, on average, 4 ± 1 labeled SVs in the presynapse at any given time.
Correlative light and EM
Rat hippocampal neurons (embryonic day 18) were grown on gridded glass-bottom dishes (P35G-1.5-14-CGRD-D; MatTek Corporation) and transfected with VAMP2–pHluorin on 14–16 DIV. On 15–17 DIV, the neurons were imaged live at 37°C on an inverted point-scanning laser confocal microscope (ZIS LMS 510 META; ZEISS), with a Plan-APO 10× 0.45 NA objective (ZEISS) using Zen 2009 software and an argon laser, to record the location of VAMP2–pHluorin–positive neurons on the cover grids. Neurons were then processed as described in the sdTIM section, with the exception of using anti-GFP HRP-mCherry-GNT and omitting single-particle imaging of SVs, and then fixed with 2.5% glutaraldehyde in PBS for 1 h at room temperature. After fixation, neurons were processed for DAB cytochemistry using standard protocols. Fixed cells were contrasted with 1% osmium tetroxide and 4% uranyl acetate before dehydration and embedded in LX-112 resin (Harper et al., 2011). Ultrathin sections (∼60 nm) were cut using an ultramicrotome (UC6FCS; Leica Biosystems). Neurons corresponding to those identified by fluorescence microscopy were visualized on a transmission electron microscope (model 1011; JEOL USA, Inc.) equipped with a cooled charge-coupled device camera (Morada; Olympus). The percentage of HRP-positive (cytochemically stained) SVs versus endocytic structures >45 nm in diameter was calculated from correlated electron micrographs of 25 VAMP2–pHluorin–expressing hippocampal neurons.
Vector quiver plot
To analyze the mobility of internalized VAMP2 in presynapses and axonal areas, previously recorded videos of sdTIM were analyzed with a custom python Qt5 GUI-based application called TrackerApp, which is available online as open-source software at https://github.com/QBI-Software/Tracking. The 2D displacement of each particle was used to calculate the polar coordinates for each trajectory. The xy coordinates of each particle position were then used to calculate the cumulative magnitude of the vectors in 100-nm2 grids, as presented in the quiver plot, generated in the open source python library at matplotlib (http://matplotlib.org; based on MATLAB; MathWorks, Inc.). The vector plot is presented as arrows on a 2D Cartesian plane with reference to the point of origin at t0.
Restimulation
After the pulse and chase uptake of VAMP2–pHluorin–bound Atto647N nanobodies, hippocampal neurons were restimulated with high K+, and the mobility of internalized Atto647N nanobodies was recorded immediately. To detect the release sites, we followed the unquenching of VAMP2–pHluorin upon restimulation as described earlier. We compared the mobility of restimulated SVs to resting SVs (sdTIM) by collecting 6,000 frames at 50 Hz in both conditions using the imaging setup described earlier. We estimate a 10-s delay in imaging after the second addition of high K+.
MSS
To obtain trajectories and the MSS, previously recorded videos of sdTIM and uPAINT were analyzed using utrack (Jaqaman et al., 2008). utrack uses Gaussian mixture models to detect particles with overlapping signals. The displacement of a particle is given by <rv> ∼ tγv and the MSS plot shows γv versus v. We only analyzed tracks that were at least eight frames in length and obtained MSS and diffusion data from these trajectories.
Applying Bayesian model selection to HMM
We marked regions of interest around the presynapses and axonal segment and analyzed the trajectories from each region of interest separately using HMM-Bayes software (Monnier et al., 2015). For comparison, we also allowed a maximum of 4 hidden states and analyzed 29 presynapses from sdTIM experiments. However, in the great majority of cases (86% of presynapses), we found three or less motion states and thus for simplicity present only the results assuming a maximum of three hidden states in Figs. 8 and 9. The sdTIM results contain quantification from presynapses and axons, which resulted in an average of 386 ± 69 and 576 ± 132 trajectories per region of interest, respectively. The corresponding numbers in restimulation analysis were 230 ± 20 and 372 ± 48, respectively. All of the analyses were performed using MATLAB (MathWorks, Inc.). The average angle between transport states was calculated from all D-DV-DV models in axons and presynapses. The apparent diffusion coefficients were not corrected for localization precision and motion blur.
Image analysis and processing
The brightness and contrast of fluorescence images were adjusted with Photoshop CC 2015.1.1 (Adobe Systems Inc.). Color scaling of trajectories and fluorescence intensity analysis of time-lapse videos were done using ImageJ (National Institutes of Health), and the change in fluorescence was normalized to the first value for each experiment. Presynapses were selected based on their punctate morphology and the increased VAMP2 fluorescence intensity upon high K+ stimulation. Axonal areas were selected based on their morphology and the absence of fluorescence variations elicited by stimulation.
Statistical analysis
All of the frequency-distributed data were tested in bins using two-tailed (unequal variances) Student’s t tests in Excel (Microsoft). We used the nonparametric Mann–Whitney U test for the AUC and the velocity magnitudes in the HMM-Bayes results and the Kolmogorov-Smirnoff test for cumulative frequency of Log10D Prism (both with Prism C for Mac OS X Version 6.0g or Version 7.0a; GraphPad Software). To compare the D/MSS distributions, a 2D Kolmogorov-Smirnoff test was performed using MATLAB (MathWorks, Inc.). The AUC and frequency distribution were calculated with Prism (GraphPad Software). Error bars are ±SEM for independent experiments unless otherwise stated. The level of significance was set to P < 0.05.
Online supplemental material
Fig. S1 describes the principle of uPAINT to track VAMP2–pHluorin–bound Atto647N nanobodies on the plasma membrane and shows the results of the MSS analysis comparing sdTIM with uPAINT methods. Fig. S2 shows the cumulative frequencies of VAMP2–pHluorin/Atto647N nanobody mobility and an example of a six-step Atto647N nanobody emission curve. Fig. S3 shows the vector analysis. Fig. S4 shows the mobility differences of recently endocytosed SVs versus those at a later time point.
Acknowledgments
The superresolution microscopy was carried out at the Queensland Brain Institute’s (QBI) Advanced Microimaging and Analysis Facility with the help of L. Hammond and R. Amor. We thank R. Martínez-Mármol (QBI) for help with the imaging, R. Gormal (QBI) for help with nanobody preparation, J. S. Borbely and M. Lakadamyali (Advanced Fluorescence Imaging and Biophysics group at The Institute of Photonic Sciences) for help with nanobody fluorescence intensity analysis, S. Manley, D. Choquet, J.-B. Sibarita, E. Hosy, and L. Cognet for insightful discussions, and R. Tweedale for critical appraisal of the manuscript.
This work was supported by an Australian Research Council Discovery Project grant (DP120104057) and an Australian Research Council Linkage Infrastructure, Equipment, and Facilities grant (LE130100078). F.A. Meunier is a Senior Research Fellow of the National Health and Medical Research Council and was the recipient of a Queensland International Fellowship. R.G. Parton is a Principal Research Fellow of the National Health and Medical Research Council and supported by National Health and Medical Research Council Program grant 1037320. P. Padmanabhan was supported by the University of Queensland Postdoctoral Research Fellowship.
The authors declare no competing financial interests.
Author contributions: M. Joensuu performed and analyzed the uPAINT and restimulation experiments. M. Joensuu and A. Papadopulos performed the sdTIM experiments. M. Joensuu analyzed the data, prepared the figures, and wrote the paper with P. Padmanabhan and F.A. Meunier. N. Durisic and M. Joensuu performed the nanobody fluorescence intensity fluctuation analysis. P. Padmanabhan, M. Joensuu, and A.T.D. Bademosi performed the Bayesian model selection to HMM, which was supervised by G.J. Goodhill. L. Cooper-Williams wrote the TrackerApp and designed the vector quiver plot analysis with A. Papadopulos. M. Joensuu and I.C. Morrow prepared the correlative EM samples. M. Joensuu, N. Durisic, and C. Harper prepared the Atto647N-tagged anti-GFP nanobodies. W. Jung prepared and characterized the peroxidase-labeled nanobodies, which was supervised by R.G. Parton. A. Papadopulos performed the MSS analysis and codesigned the study with F.A. Meunier. F.A. Meunier supervised the study.
References
- AUC
area under the mean square displacement curve
- DIV
days in vitro
- GNT
GFP nanotrap
- HMM
hidden Markov model
- MSD
mean square displacement
- MSS
moment scaling spectrum
- sdTIM
subdiffractional tracking of internalized molecules
- SV
synaptic vesicle
- uPAINT
Universal Point Accumulation Imaging in Nanoscale Topography
- VAMP2
vesicle-associated membrane protein 2
Author notes
C.B. Harper’s present address is Centre for Integrative Physiology, University of Edinburgh, Edinburgh EH8 9XD, Scotland, UK.