Characterization of developmental and molecular factors underlying release heterogeneity at Drosophila synapses

Neurons communicate through neurotransmitter release at specialized synaptic regions known as active zones (AZs). Using biosensors to visualize single synaptic vesicle fusion events at Drosophila neuromuscular junctions, we analyzed the developmental and molecular determinants of release probability (Pr) for a defined connection with ~300 AZs. Pr was heterogeneous but represented a stable feature of each AZ. Pr remained stable during high frequency stimulation and retained heterogeneity in mutants lacking the Ca2+ sensor Synaptotagmin 1. Pr correlated with both presynaptic Ca2+ channel abundance and Ca2+ influx at individual release sites. Pr heterogeneity also correlated with glutamate receptor abundance, with high Pr connections developing receptor subtype segregation. Intravital imaging throughout development revealed that AZs acquire high Pr during a multi-day maturation period, with Pr heterogeneity largely reflecting AZ age. The rate of synapse maturation was activity-dependent, as both increases and decreases in neuronal activity modulated glutamate receptor field size and segregation.

The Drosophila neuromuscular junction (NMJ) has emerged as a useful system to study release heterogeneity. At this connection, motor neurons form glutamatergic synapses onto bodywall muscles in a stereotypical fashion, with the axon expanding to form~10-60 synaptic boutons that each contain many individual AZs (Harris and Littleton, 2015). Drosophila AZs contain a similar assortment of proteins to those identified at mammalian AZs (Bö hme et al., 2016;Bruckner et al., 2017Bruckner et al., , 2012Ehmann et al., 2014;Feeney et al., 1998;Fouquet et al., 2009;Graf et al., 2012;Jan and Jan, 1976;Kaufmann et al., 2002;Kittel et al., 2006;Liu et al., 2011;Owald et al., 2010;Wagh et al., 2006). Each AZ is specifically associated with a postsynaptic glutamate receptor field. Glutamate receptors at the Drosophila NMJ are excitatory inotropic non-NMDA receptors that exist as tetramers, with three obligatory subunits encoded by GluRIII, GluRIID and GluRIIE, and a variable fourth subunit encoded by either GluRIIA (A-type) or GluRIIB (B-type) Marrus et al., 2004;Petersen et al., 1997;Qin et al., 2005;Schuster et al., 1991). GluRIIA containing receptors generate a larger quantal size and display slower receptor desensitization than their GluRIIB counterparts (DiAntonio et al., 1999). The A-and B-subtypes compete for incorporation into the tetramer at individual postsynaptic densities (PSDs) in a developmental and activityregulated fashion (Chen and Featherstone, 2005;DiAntonio et al., 1999;Marrus and DiAntonio, 2004a;Rasse et al., 2005;Schmid et al., 2008). eLife digest To send a message to its neighbor, a neuron releases chemicals called neurotransmitters into the gap -or synapse -between them. The neurotransmitter molecules bind to proteins on the receiver neuron called receptors. But what causes the sender neuron to release neurotransmitter in the first place? The process starts when an electrical impulse called an action potential arrives at the sender cell. Its arrival causes channels in the membrane of the sender neuron to open, so that calcium ions flood into the cell. The calcium ions interact with packages of neurotransmitter molecules, known as synaptic vesicles. This causes some of the vesicles to empty their contents into the synapse.
But this process is not particularly reliable. Only a small fraction of action potentials cause vesicles to fuse with the synaptic membrane. How likely this is to occur varies greatly between neurons, and even between synapses formed by the same neuron. Synapses that are likely to release neurotransmitter are said to be strong. They are good at passing messages from the sender neuron to the receiver. Synapses with a low probability of release are said to be weak. But what exactly differs between strong and weak synapses?
Akbergenova et al. studied synapses between motor neurons and muscle cells in the fruit fly Drosophila. Each motor neuron forms several hundred synapses. Some of these synapses are 50 times more likely to release neurotransmitter than others. Using calcium imaging and genetics, Akbergenova et al. showed that sender cells at strong synapses have more calcium channels than sender cells at weak synapses. The subtypes and arrangement of receptor proteins also differ between the receiver neurons of strong versus weak synapses. Finally, studies in larvae revealed that newly formed synapses all start out weak and then gradually become stronger. How fast this strengthening occurs depends on how active the neuron at the synapse is.
This study has shown, in unprecedented detail, key molecular factors that make some fruit fly synapses more likely to release neurotransmitter than others. Many proteins at synapses of mammals resemble those at fruit fly synapses. This means that similar factors may also explain differences in synaptic strength in the mammalian brain. Changes in the strength of synapses underlie the ability to learn. Furthermore, many neurological and psychiatric disorders result from disruption of synapses. Understanding the molecular basis of synapses will thus provide clues to the origins of certain brain diseases.
The stereotypical alignment of individual AZs to distinct postsynaptic glutamate receptor fields in Drosophila allowed the generation of genetic tools to optically follow quantal fusion events at single release sites by visualizing glutamate receptor activation (Melom et al., 2013;Peled and Isacoff, 2011). Classically, studies of synaptic transmission have used electrophysiology to measure the postsynaptic effect of neurotransmitter release over a population of release sites (Katz andMiledi, 1969, 1967), precluding an analysis of how individual AZs contribute to the evoked response. By transgenically expressing GCaMP Ca 2+ sensors that target to the postsynaptic membrane, single vesicle fusion events at each individual AZ can be imaged by following spatially localized Ca 2+ influx induced upon glutamate receptor opening. This allows for the generation of P r maps for both evoked and spontaneous fusion for all AZs (Cho et al., 2015;Melom et al., 2013;Muhammad et al., 2015;Newman et al., 2017;Peled et al., 2014;Peled and Isacoff, 2011;Reddy-Alla et al., 2017). One surprising observation using this quantal imaging approach is that AZs formed by a single motor neuron have a heterogeneous distribution of P r , ranging from 0.01 to~0.5, with neighboring AZs often showing~50-fold differences in P r (Melom et al., 2013;Peled et al., 2014;Peled and Isacoff, 2011). These differences in P r result in AZs with distinct short-term plasticity properties, suggesting release heterogeneity has functional importance for synaptic transmission (Peled and Isacoff, 2011).
Key questions raised by these observations include how P r is uniquely set for individual AZs and how the heterogeneity in P r arises during development. P r variability is likely to be controlled in part by variable Ca 2+ channel abundance at release sites, consistent with the heterogeneity in VGCCs and other associated AZ proteins previously documented (Bö hme et al., 2016;Ehmann et al., 2014;Fulterer et al., 2018;Graf et al., 2012Graf et al., , 2009Guerrero et al., 2005;Peled and Isacoff, 2011). Indeed, P r has been shown to correlate with BRP abundance at AZs in Drosophila (Paul et al., 2015;Peled et al., 2014;Reddy-Alla et al., 2017) and BRP has a key function in clustering VGCCs . Furthermore, P r has previously been shown to correlate with the number of VGCCs at several vertebrate synapses (Chen et al., 2015;Nakamura et al., 2015;Sheng et al., 2012). P r could also be regulated by local synaptic vesicle pools and their number and/or state (i.e. phosphorylation status). Beyond the molecular factors that determine AZ P r , it is unclear how release heterogeneity at Drosophila NMJs arises during development. Intravital imaging at third instar larval NMJs has demonstrated that AZs are born small and gain pre-and postsynaptic components over time in a sequential manner (Andlauer and Sigrist, 2012;Fouquet et al., 2009;Füger et al., 2007;Rasse et al., 2005;Zhang et al., 2010). However, this approach has not been used during earlier stages of larval development to determine whether the release heterogeneity observed at third instar NMJs reflects AZ birth order. To characterize factors regulating P r at individual AZs, as well as the origin of P r diversity, we employed optical quantal analysis and intravital imaging to examine how P r heterogeneity arises during development.

Results
Drosophila NMJ synapses display heterogeneity in P r , ranging from functionally silent sites to high P r AZs Recent studies indicate that release sites possess structural and functional heterogeneity (Éltes et al., 2017;Holderith et al., 2012;Maschi and Klyachko, 2017;Melom et al., 2013;Peled et al., 2014;Peled and Isacoff, 2011;Reddy-Alla et al., 2017;Sugie et al., 2015). Using the Drosophila NMJ, we previously observed that evoked P r is non-uniform across a population of~300 AZs formed by motor neuron MN4-Ib onto muscle 4, ranging from 0.01 to~0.5 in HL3 saline containing 1.3 mM extracellular Ca 2+ and 20 mM Mg 2+ (Melom et al., Video 1. Representative movie showing evoked and spontaneous GCaMP6s events (green) in larvae expressing GluRIIA-RFP (red) that were stimulated at 0.3 Hz. DOI: https://doi.org/10.7554/eLife.38268.003 2013). In our original study, each AZ was identified by the location of postsynaptic Ca 2+ flashes, but AZs were not directly labeled in the live preparation. To more precisely map AZ P r heterogeneity, we identified the position of each corresponding PSD by co-expressing the RFP-tagged glutamate receptor subunit GluRIIA under the control of its endogenous promoter  along with a newer version of our previous biosensor, N-terminal myristoylated GCaMP6s, expressed in muscles using Mef2-GAL4. We monitored postsynaptic Ca 2+ influx from activation of glutamate receptors after either spontaneous release or nerve stimulation (0.3 Hz for 5 min) in muscle 4 of early stage third instar larvae (Video 1).
Using this approach, we mapped all myrGCaMP6s visualized release events to the position of in vivo GluRIIA-RFP labeled PSDs ( Figure 1A). Consistent with previous data, we observed a heterogeneous distribution of AZ P r , with an average P r of 0.073 ± 0.002 (n = 1933 AZs from 16 NMJs from 16 animals). However, there was a~50-fold difference in P r between the highest and lowest releasing sites. The AZ P r dataset did not fit a normal distribution (D'Agostino K 2 test (p<0.0001), Shapiro-Wilk test (p<0.0001), Kolmogorov-Smirnov test (p<0.0001)) and instead was skewed to the right, with a majority of AZs rarely releasing a synaptic vesicle following an action potential (P r in the range of 0.01 to 0.2) and a small number of AZs consistently showing high release rates (75% percentile of P r was 0.1, with a maximum P r of 0.73; Figure 1B). 9.7% of all release sites defined by their apposed GluRIIA receptors displayed only spontaneous fusion events, and another 14.6% of the AZ population was silent for both spontaneous and evoked release during the recording period ( Figure 1B). We categorized all AZs with a release rate greater than two standard deviations above average as 'high P r ', and the remaining AZs that showed evoked release as 'low P r '. Using these criteria, 65.8% of all AZs fell in the low P r category with an average P r of 0.049 ± 0.004. In contrast, 9.9% of AZs were classified as high P r sites, with an average P r of 0.277 ± 0.015 ( Figure 1C). High P r AZs displayed on average a 5.7-fold higher chance of vesicle fusion following an action potential compared to low P r AZs.
High P r AZs correspond to single release sites with enhanced levels of the AZ protein Bruchpilot One potential caveat to the interpretation of P r heterogeneity is the possibility that multiple closelypositioned release sites could be falsely characterized as single high P r AZs using conventional light microscopy. Presynaptic AZ position can be precisely identified at the NMJ by labeling the core AZ T-bar component Bruchpilot (BRP), the homolog of mammalian ELKS/CAST (Fouquet et al., 2009;Wagh et al., 2006). To determine if the high P r sites we observed were actually due to release from closely clustered AZs, we used high-resolution structured illumination microscopy (SIM) on fixed tissue stained with anti-BRP ( Figure 1D) following dual color (myrGCaMP6s/GluRIIA-RFP) quantal imaging. SIM provides a lateral resolution of~110 nm (Wegel et al., 2016), providing clear resolution of the BRP ring structure, which has a diameter of~200 nm (Owald et al., 2012), smaller than the resolution limit of conventional light microscopy. The presence of GluRIIA-RFP allowed us to precisely match individual high P r sites observed during live imaging with their position in fixed and stained tissue during SIM imaging ( Figure 1D). Using an automated detection algorithm in the Volocity 3D image analysis software, we were able to identify all AZs labeled with BRP ( Figure 1D, far right panel), and to resolve individual AZ clusters that were not separated using conventional spinning disk microscopy ( Figure 1D, white box). SIM analysis of distances between neighboring AZs indicated that 2.45 ± 0.4% (n = 9 NMJs from nine animals) of all AZs were located close enough to each other (within 280 nm) such that they would not be resolvable during live imaging. In contrast, 9.9% (n = 16 NMJs from 16 animals) of AZs were functionally classified as high P r sites, suggesting that the majority of these sites are not likely to be explained by release from multiple closely clustered AZs.
SIM visualization of BRP at AZs following P r mapping revealed that a majority of high P r sites were represented by a single BRP-positive AZ that was not further resolvable after SIM ( Figure 1D, red circles). We identified high P r sites (n = 42 AZs from 5 NMJs from five animals) and then determined what fraction of these sites were truly clusters of multiple neighboring AZs based on their SIM profiles. Only 16 ± 3% of these high P r sites were resolved into multiple AZs upon SIM analysis, indicating that most high P r AZs correspond to single release sites. These single BRP clusters at high P r sites had larger sum fluorescence intensities than most other BRP positive puncta ( Figure 1E). The average sum fluorescence of single BRP puncta from high P r AZs (3.95 Â 10 6 ± 2.67 x 10 5 , n = 24 AZs from 9 NMJs from nine animals) was 1.7-fold greater than the fluorescence of randomly selected low P r BRP clusters (2.33 Â 10 6 ± 0.98 x 10 5 , n = 60 AZs from 9 NMJs from nine animals, p<0.0001).
To further examine large single BRP clusters that could not be resolved using conventional spinning disk microscopy, all BRP clusters larger than 280 nm were automatically detected and assigned their P r measured during live imaging. We then determined whether these sites were represented by single or multiple AZs using SIM microscopy. Clusters > 280 nm in diameter that could be resolved to multiple BRP positive AZs after SIM imaging had a lower P r (0.10 ± 0.02, n = 35 AZs from 5 NMJs from five animals) than those comprised of a single large BRP positive AZ (0.19 ± 0.02, n = 42 AZs from 5 NMJs from five animals; Figure 1F). As such, high resolution SIM analysis confirms that most high P r sites correspond to single AZs with more intense BRP labeling, consistent with previous data supporting the positive role of BRP in regulating P r (Paul et al., 2015;Peled et al., 2014;Reddy-Alla et al., 2017).

Release heterogeneity is retained in Synaptotagmin one mutants
Heterogeneous release rates between AZs could solely reflect stable differences in protein content of the AZs themselves. However, the accumulation of different synaptic vesicle populations with variable levels of Ca 2+ sensitivity or fusogenicity might also contribute to release heterogeneity. The synchronous Ca 2+ sensor Synaptotagmin 1 (Syt1) resides on synaptic vesicles and plays a major role in P r determination at Drosophila NMJs (DiAntonio and Schwarz, 1994;Guan et al., 2017;Lee et al., 2013;Littleton et al., 1994Littleton et al., , 1993Yoshihara et al., 2003;Yoshihara and Littleton, 2002). We hypothesized that if differential synaptic vesicle Ca 2+ sensitivity is a major determinant of release heterogeneity in addition to its established role in determining overall P r , then elimination of Syt1 would disrupt P r heterogeneity. Consistent with electrophysiological findings, quantal imaging in syt1 null mutants expressing GluRIIA-RFP and myrGCaMP6s revealed a dramatic reduction in evoked release, a shift from synchronous to highly asynchronous The percentage of AZs that were low P r (65.8%), high P r (9.9%), spontaneous-only (9.7%) and silent (14.6%) is displayed in the inset. (C) Average P r determined for each individual experiment for the AZ population categorized based on low and high activity sites (>2 standard deviations above the mean). Each point represents the average for all AZs (classified as either high or low P r ) from a single animal. (D) Individual BRP puncta for three NMJs from three different animals imaged with high resolution structured illumination microscopy (SIM, left panel) or confocal microscopy (middle panel). The right panel displays the heat map for evoked P r from the same NMJs determined by GCaMP6s imaging prior to fixation. Representative high P r sites are circled with red in the middle panels. Representative example of a large BRP puncta that would be classified as a single AZ using conventional light microscopy but resolved into multiple clustered AZs using SIM analysis is boxed with the dotted white line. The far right top panel displays the results from the automated detection algorithm that outlines individual AZs. (E) AZs were separated into high and low P r based on their activity and the sum fluorescence intensity of the corresponding BRP puncta is shown. (F) AZs with high BRP intensity (two standard deviations above average) were preselected from conventional confocal images and identified on corresponding SIM images. In cases where the BRP signal was resolvable into more than one AZ by SIM microscopy, it was assigned to the AZ cluster group. In cases where the BRP signal mapped to a single BRP puncta by SIM imaging, it was assigned to the single AZ group. P r is plotted for each group. Student's t-test was used for statistical analysis (***=p 0.001 fusion, and an increase in spontaneous release rates (Video 2). To estimate AZ release heterogeneity in syt1 nulls, preparations were stimulated at 5 Hz and release events were normalized to the number of stimuli ( Figure 2A). The average release rate per AZ per second in syt1 nulls during 5 Hz stimulation was 0.03 ± 0.001 (n = 719 AZs from 7 NMJs from six animals; Figure 2B). In contrast, spontaneous release rate per AZ in the absence of stimulation was 0.018 ± 0.001 per second in syt1 nulls (n = 719 AZs from 7 NMJs from six animals) compared to 0.011 ± 0.001 in controls (n = 559 AZs from 6 NMJs from four animals, p<0.0001; Figure 2B).
Although release rate is dramatically reduced in syt1 nulls, AZs still maintained overall heterogeneity in P r distribution. Comparing the distribution of AZ release rates for syt1 nulls and controls, release was proportionally decreased across all AZs in syt1 ( Figure 2C); frequency distribution analysis of AZs with normalized release rates (from 0 to maximum release) confirmed that there was no significant change in the heterogeneity of release between syt1 mutants and controls ( Figure 2D,E). Given that AZ P r remains highly heterogeneous in the absence of Syt1, these data indicate that heterogeneity in synaptic vesicle Ca 2+ sensitivity between AZs is unlikely to play a major role in P r distribution.
Individual AZ P r remains stable through extensive vesicle cycling Release heterogeneity in syt1 null animals suggests that synaptic vesicle Ca 2+ sensitivity is not a major determinant of P r heterogeneity at this synapse; however, it is possible that other synaptic vesicle components influence P r on an AZ-specific level. If unique local synaptic vesicle pools contribute to the distribution of P r between AZs, we predicted that P r would be highly dynamic at individual AZs over time. In contrast, stability of P r would argue that release heterogeneity is more likely associated with stably resident proteins at individual AZs. To assess P r stability over time, we conducted a 3 min imaging session using 0.3 Hz stimulation to generate an initial P r map, and then allowed the preparation to rest for 5 min without stimulation or imaging before re-mapping P r in a final 3 min imaging session. We were constrained in our ability to examine P r continuously over longer time intervals due to bleaching of GCaMP6s from the high frequency capture rate. P r at individual AZs was very stable between the two sessions ( Figure 3A). This was especially evident for high P r sites, which sustained high levels of activity during both imaging sessions. Plotting release rate for all AZs revealed a strong correlation for P r across the two imaging sessions (Pearson r = 0.77, R 2 = 0.59, p<0.0001, n = 988 AZs from 8 NMJs from seven animals; Figure 3B).
We next used a strong stimulation paradigm to drive vesicle cycling to promote intermixing of synaptic vesicles. NMJ preparations were imaged during two low frequency 0.3 Hz stimulation periods separated by a 5 min 5 Hz stimulation session ( Figure 3C). Release maps were not dramatically altered by 5 Hz stimulation, with the overall correlation of P r between the two imaging sessions similar with and without stimulation (Pearson r = 0.78, R 2 = 0.61, p<0.0001, n = 613 AZs from 6 NMJs from six animals; Figure 3D). Thus, inducing vesicle recycling with 5 Hz stimulation does not dramatically change P r across the AZ population, arguing that stably resident AZ components, rather than AZ-specific synaptic vesicle populations, are likely to represent the major driver of P r heterogeneity at this synapse. P r is correlated with Ca 2+ channel abundance at AZs Given that vesicle fusion is highly sensitive to Ca 2+ and most effective in close proximity to VGCCs (Augustine et al., 1985;Bö hme et al., 2016;Chen et al., 2015;Heidelberger et al., 1994;Katz and Miledi, 1967;Katz, 1969Zito et al., 1999Keller et al., 2015;Meinrenken et al., 2003Meinrenken et al., , 2002Nakamura et al., 2015;Sheng et al., 2012;Stanley, 2016;Wang et al., 2008), Ca 2+ channel abundance at individual AZs is predicted to be a key variable for P r heterogeneity at Drosophila NMJs as well. Cacophony (cac) encodes the Drosophila voltage-activated Ca 2+ channel a1 subunit required for neurotransmitter release (Fouquet et al., 2009;Kawasaki et al., 2004;Littleton and Ganetzky, 2000;Liu et al., 2011;Rieckhof et al., 2003;Smith et al., 1996). Transgenic animals expressing fluorescently tagged Cac channels have been previously generated, demonstrating that Cac localizes specifically to AZs and its abundance appears heterogenous across release sites (Kawasaki et al., 2004;Matkovic et al., 2013;Yu et al., 2011).
To examine the heterogeneity of Cac abundance across AZs, we used SIM to measure the distribution of Cac-GFP and BRP at muscle 4 NMJs. Using this approach, we observed a heterogeneous distribution of mean Cac-GFP fluorescence at AZs, similar to the variable levels of BRP intensity described earlier ( Figure 1D, Figure 4-figure supplement 1A,B). 5.72% of AZs displayed Cac-GFP fluorescence greater than two standard deviations above average (n = 2011 AZs from 11 NMJs from three animals). The mean Cac-GFP fluorescence for these bright AZs was 2.1-fold greater than that observed for the remaining sites (p<0.0001; Figure 4-figure supplement 1C). Mean AZ intensities of Cac-GFP and BRP were positively correlated (Pearson r = 0.46, R 2 = 0.21, p<0.0001, n = 730 AZs from 6 NMJs from three animals; Figure 4-figure supplement 1D,E).
To determine whether Cac distribution correlates with P r heterogeneity, we used dual color imaging experiments where vesicle fusion events were detected by myrGCaMP6s driven in muscles using mef2-GAL4 and Ca 2+ channel distribution was visualized using red-labeled Cac-TdTomato expressed pan-neuronally using elav-GAL4 ( Figure 4A). We observed a strong positive correlation (average Pearson r = 0.61, R 2 = 0.38, p<0.0001, n = 483 AZs from 7 NMJs from seven animals) between AZ Cac fluorescence intensity and evoked AZ P r ( Figure 4B). In contrast, single AZ release rates for spontaneous events showed only a mild correlation with Cac intensity (average Pearson r = 0.19, R 2 = 0.036, p<0.0001, n = 483 AZs from 7 NMJs from seven animals; Figure 4C). These results are The corresponding P r heatmap is shown on the right. The arrows denote several high P r sites opposed by bright GluRIIA positive PSDs. (B) AZ release events per second for spontaneous release and evoked by 5 Hz stimulation are shown for syt1 nulls mutants, and for spontaneous release in controls. (C) Frequency distribution of P r is shown for syt1 nulls and controls. (D) Plot of normalized P r frequency distribution (from 0 to 1 (max)) for syt1 nulls and controls. (E) Cumulative frequency distribution for normalized release rates for syt1 nulls and controls is shown. One-way ANOVA followed by Dunnett's multiple comparisons test used for statistical analysis (***=p 0.001). Error bars represent SEM. DOI: https://doi.org/10.7554/eLife.38268.005 Figure 3. Stability of release maps at the NMJ. (A) P r heatmaps for the same muscle 4 NMJ were generated for two individual imaging sessions, separated by a 5 min resting period. High P r AZs were numbered and re-identified in each heatmap. Representative high P r AZs that sustain release rates during the second imaging session are noted with arrows. (B) Correlation of AZ P r between two imaging sessions separated by a 5 min resting Figure 3 continued on next page consistent with previous observations that release rates for evoked and spontaneous fusion are not correlated at Drosophila AZs (Melom et al., 2013;Peled et al., 2014), and that spontaneous fusion is largely independent of extracellular Ca 2+ at this synapse (Jorquera et al., 2012;Lee et al., 2013).
To increase confidence that the observed Cac-TdTomato intensity accurately reflects Cac distribution, we also measured the correlation between P r and Cac channels transgenically tagged with GFP. We generated transgenic lines expressing myristoylated red Ca 2+ indicators previously characterized in the field, including RCaMP1h, R-GECO1 and jRGECO1a. Although RCaMP1h and R-GECO1 were too dim to visualize localized Ca 2+ transients at PSDs, transgenic lines expressing the myristoylated Ca 2+ sensor jRGECO1a (Dana et al., 2016) in muscle four allowed detection of Ca 2+ influx following vesicle fusion at single AZs (Video 3). In contrast to the more robust GCaMP6s, jRGECO1a has a shorter fluorescent lifetime and the signal amplitude decays more rapidly. We observed that quantal events imaged with myr-jRGECO1a were dimmer and fully bleached within 7-10 min of imaging. Therefore, preparations were stimulated at 1 Hz for shorter two-minute imaging sessions to generate P r maps in larvae expressing myr-jRGECO1a ( Figure 4D). We observed a strong correlation between AZ P r detected by myr-jRGECO1a and Cac-GFP intensity (average Pearson r = 0.54, R 2 = 0.29, p<0.0001, n = 651 AZs from 7 NMJs from seven animals, correlation from a representative experiment shown in Figure 4E). Again, we found a weaker correlation between spontaneous fusion and Cac-GFP intensity (Pearson r = 0.17, R 2 = 0.03, p<0.0001, n = 651 AZs from 6 NMJs from six animals). Hence, P r for action-potential evoked fusion is strongly correlated with the local abundance of Cac channels at individual Drosophila AZs.
We next compared Cac-GFP fluorescence at AZs that were functionally classified as either low or high P r sites by quantal imaging using myr-jRGECO1a ( Figure 4F). The average fluorescence of single Cac-GFP puncta from high P r AZs (normalized intensity = 0.6 ± 0.04, n = 38 AZs from 7 NMJs from seven animals) was 2.09-fold greater than the average fluorescence of low P r AZs (normalized intensity = 0.29 ± 0.01, n = 638 AZs from 7 NMJs from seven animals, p<0.0001). We also examined the P r of AZs classified by Cac-GFP fluorescence. The average P r for AZs displaying high Cac-GFP fluorescence (>2 standard deviations above average) was 0.2 ± 0.016 (n = 7 NMJs from seven animals) compared to 0.06 ± 0.003 (n = 7 NMJs from seven animals, p<0.0001) for the remaining AZs with lower levels of Cac-GFP. Although the absolute number of Cac channels at single Drosophila AZs is unknown, these data suggest that a ! 2 fold difference in channel number exists between low and high P r AZs. Given the steep third to fourth order non-linear dependence of synaptic vesicle fusion on Ca 2+ (Dodge and Rahamimoff, 1967;Heidelberger et al., 1994;Jan and Jan, 1976), a small change in channel number is likely to have a large effect on P r .

P r correlates with the level of presynaptic Ca 2+ influx at individual AZs
VGCCs are extensively modulated by second messenger pathways that can alter channel conductivity (Catterall and Few, 2008;Dolphin et al., 1991;Evans and Zamponi, 2006;Reid et al., 2003;Tedford and Zamponi, 2006;Zamponi and Snutch, 1998). Although Cac channel abundance correlates with AZ P r , an important readout of channel activity is the local Ca 2+ influx occurring at each AZ. Assaying presynaptic Ca 2+ influx directly would also be useful to bypass any unknown effects on P r generated by expressing fluorescently tagged Ca 2+ channels. As such, we generated transgenic animals expressing GCaMP6m fused to the N-terminus of BRP, which localizes directly to the base of the AZ where Ca 2+ channels cluster (Fouquet et al., 2009;Kittel et al., 2006). These GCaMP6m fusions were made to a BRP fragment (BRP short ) corresponding to amino acids 473-1226 of the full 1740 amino acid protein (Schmid et al., 2008). At rest, N-terminal GCaMP-BRP short was dim, consistent with low levels of resting Ca 2+ ( Figure 5A). Stimulation at 10 Hz resulted in a robust increase in discrete punctated presynaptic fluorescence that remained confined to single AZs during stimulation ( Figure 5A). During multiple rounds of 5 s 10 Hz stimulation, GCaMP-BRP short fluorescence increase period. (C) P r heatmaps for the same NMJ separated by a 5 min 5 Hz stimulation. Representative high P r AZs that did not change activity levels are noted with arrows. (D) Correlation of AZ P r between two imaging sessions separated by a 5 min 5 Hz stimulation period. DOI: https://doi.org/10.7554/eLife.38268.007 (DF) varied between AZs, but was stable at the same AZ for each independent stimulation (n = 205 AZs from 6 NMJs from three animals; Figure 5B,C).
Given that GCaMP-BRP short abundance at an AZ likely reflects the absolute amount of BRP at that AZ, we assayed if heterogeneity in GCaMP-BRP short fluorescence during 10 Hz stimulation could be solely explained by differences in sensor distribution across AZs. We applied 200 nM of the Ca 2+ ionophore ionomycin to elevate Ca 2+ concentrations uniformly throughout the terminal independent of Cac abundance. In the presence of ionomycin, differences in fluorescent signals between AZs should be entirely due to heterogeneity in sensor abundance. We observed a rightward shift in the GCaMP-BRP short intensity distribution among AZs upon ionomycin application compared to 10 Hz stimulation (average AZ fluorescence during 10 Hz stimulation was 1924 ± 63, and following ionomycin addition was 6105 ± 175; Figure 5-figure supplement 1A,B), indicating that during 10 Hz stimulation, detection of Ca 2+ by GCaMP-BRP short is not limited by sensor abundance. Furthermore, we observed a significant difference in the shape of the distribution during 10 Hz stimulation compared to both before stimulation and after ionomycin. The distribution of fluorescence intensities is narrower both at rest and upon ionomycin application; these two distributions should primarily reflect sensor distribution. In contrast, the distribution of GCaMP-BRP short fluorescence during 10 Hz stimulation is wider, indicating that the sensor is reporting local changes in Ca 2+ influx and not just sensor distribution (Figure 5, Figure 5-figure supplement 1C,D). Thus, although GCaMP-BRP short abundance is likely to contribute to the levels of Ca 2+ influx detected, these results are consistent with heterogeneity in Ca 2+ influx across individual AZs.
We next assayed if Ca 2+ influx detected by GCaMP-BRP short is correlated with Cac channel abundance. Animals expressing both Cac-TdTomato and GCaMP-BRP short transgenes in the presynaptic compartment displayed a strong correlation between Ca 2+ dependent excitation of GCaMP-BRP short (DF) and Cac-TdTomato intensity at individual AZs during stimulation (Pearson r = 0.73, R 2 = 0.53, p<0.0001, n = 176 AZs from 7 NMJs from six animals; Figure 5D, E). We observed a weaker correlation between Cac intensity and GCaMP-BRP short DF at rest (Pearson r = 0.18, R 2 = 0.03, p<0.001, n = 338 AZs from 8 NMJs from six animals). A few Evoked release visualized at the same NMJ by myr-jRGECO1a is shown in the second panel. The P r heatmap for evoked release is shown in the third panel. A heatmap distribution of Cac-GFP fluorescence intensities, based on same two standard deviation criteria as color-coding of P r , is shown in the right panel. The arrows denote several higher P r sites that correlated with bright Cac-GFP puncta. (E) Correlation between AZ P r and Cac-GFP fluorescence intensity for evoked release. (F) Cac-GFP fluorescence for AZs functionally classified as either low or high P r (>2 standard deviations above mean) by quantal imaging with myr-jRGECO1a. Student's t-test was used for statistical analysis (***=p 0.001). Error bars represent SEM. DOI: https://doi.org/10.7554/eLife.38268.008 The following figure supplement is available for figure 4:  instances were noted where specific AZs experienced a disproportionally low GCaMP-BRP short DF signal relative to their Cac-TdTomato intensity ( Figure 5D, arrows), suggesting Ca 2+ influx may be fine-tuned at certain release sites. We next analyzed the correlation between GCaMP-BRP short DF induced by 10 Hz stimulation and release rate visualized by postsynaptic myr-jRGECO1a during 1 Hz stimulation ( Figure 6A,B). AZs that experienced stronger Ca 2+ influx displayed higher P r during stimulation. Overall, there was a strong correlation between Ca 2+ influx and AZ P r (Pearson r = 0.56, R 2 = 0.31, p<0.0001, n = 492 AZs from 6 NMJs from six animals; Figure 6B). In contrast, the frequency of spontaneous vesicle fusion per AZ was only mildly correlated with GCaMP-BRP short DF (Pearson r = 0.23, R 2 = 0.07, n = 492 AZs from 6 NMJs from six animals; correlations from a representative experiment shown in Figure 6C). It is worth noting that although a strong correlation between Ca 2+ influx and evoked P r was observed at most AZs, a minority population of release sites that displayed robust Ca 2+ influx had very low P r ( Figure 6B).
To functionally test if the level of Ca 2+ influx rather than the structural presence of the Ca 2+ channel is responsible for determining P r , we generated P r maps in the cac NT27 mutant using dual color quantal imaging with GluRIIA-RFP and myr-GCamP6s ( Figure 6-figure supplement 1A). Cac NT27 channels have reduced Ca 2+ conductance due to a point mutation in the S4 voltage sensor (Rieckhof et al., 2003). We observed that cac NT27 results in a global decrease in P r across AZs; evoked P r in controls ranged from 0 to 0.73 with an average of 0.073 ± 0.0021 (n = 1933 AZs from 16 NMJs from 16 animals), while cac NT27 P r ranged from 0 to 0.47 with a significantly lower average P r of 0.049 ± 0.0045 (n = 275 AZs from 5 NMJs from five animals; Figure 6-figure supplement 1B). While the entire P r distribution was shifted lower compared to controls (Figure 6-figure supplement 1C), the normalized distribution of P r was nearly identical to controls (Figure 6-figure supplement 1D). These findings indicate that the levels of Ca 2+ influx through Cac channels, rather than the physical presence of the channels, is a primary determinant of P r .

Segregation of postsynaptic glutamate receptor subunits at high P r AZs
We next examined if postsynaptic glutamate receptor composition varied at low P r versus high P r AZs. Drosophila glutamate receptors at the NMJ assemble as heteromeric tetramers containing three essential subunits (GluRIII, IID and IIE) and a variable fourth subunit of GluRIIA or GluRIIB, with the GluRIIA subtypes having a higher conductance than GluRIIB Marrus et al., 2004;Petersen et al., 1997;Qin et al., 2005;Schuster et al., 1991). To determine if the GluR subtypes differentially accumulate at AZs in a manner that correlates with presynaptic P r , we visualized GluRIIA-RFP and GluRIIB-GFP expressed under the control of their endogenous promoters . To image myrGCaMP6s activity without obscuring GluRIIB-GFP, myrG-CaMP6s was expressed at low levels using the LexA/LexOP system in muscle four with Mef2-LexA. LexA driven myrGCaMP6s signal is dimmer than UAS-myrGCaMP6s at rest and does not obscure the brighter GluRIIB-GFP PSD puncta ( Figure 7A). However, upon Ca 2+ binding to myrGCaMP6s, the fluorescence dramatically increases compared to the level of endogenous GluRIIB-GFP signal, allowing simultaneous imaging of baseline GluRIIB-GFP levels and synaptic activity detected by myrGCaMP6s (Video 4).
Simultaneous expression of GluRIIA-RFP and GluRIIB-GFP revealed a heterogeneous distribution of each subunit across the AZ population, with GluRIIA levels appearing more variable than GluRIIB ( Figure 7A). Similar to the relatively sparse localization of high P r AZs (Figure 1), a similar sparse distribution of AZs apposed by very bright GluRIIA fields was observed ( Figure 7A). To determine if AZs that preferentially accumulate high levels of GluRIIA corresponded to high P r release sites, we mapped P r across the AZ population in GluRIIA-RFP/GluRIIB-GFP expressing animals. Analysis of the P r map revealed a strong positive correlation between mean GluRIIA-RFP intensity and P r (Pearson r = 0.56, R 2 = 0.32, p<0.0001, n = 756 AZs from 8 NMJs from four animals; Figure 7B). In contrast, correlation with the levels of mean GluRIIB-GFP intensity was weaker (Pearson r = 0.32, R 2 = 0.1, p<0.0001, n = 756 AZs from 8 NMJs from four animals; Figure 7C). Heterogeneity in GluRIIA PSD brightness and AZ P r was also observed in syt1 mutants (Figure 2A, arrows). These findings are consistent with previous observations that glutamate receptors preferentially cluster at sites with high P r based on electrophysiological studies in a Drosophila GluRIII hypomorphic mutant (Marrus and DiAntonio, 2004b). As observed in controls (Figure 1), a positive correlation between AZ P r and BRP levels was also observed in larvae expressing labeled glutamate receptors (Pearson r = 0.44, R 2 = 0.2, p<0.0001, n = 399 AZs from 6 NMJs from four animals; Figure 7D). In summary, these data indicate that GluRIIA accumulates more at PSDs apposing high P r AZs. Beyond the preferential GluRIIA accumulation at high P r sites, a change in GluRIIA/B distribution within single PSDs was also observed. PSDs associated with the highest P r AZs showed a segregated distribution of the receptor subtypes, with GluRIIA concentrating in the center of the receptor field immediately apposing the presynaptic BRP cluster ( Figure 7E, arrows). At these sites, GluRIIB occupied a more peripheral position around the central GluRIIA cluster. A similar localization pattern with a ring of GluRIIB surrounding a central GluRIIA patch was previously noted with antibody staining for the two receptors at a population of AZs in wildtype late third instar larvae (Marrus et al., 2004). To analyze this receptor segregation in greater detail, GluRIIA-RFP and GluRIIB-GFP were examined in the absence of co-expressed myrGCaMP6s. Prior analysis ( Figure 7B) indicated the brightest GluRIIA PSDs corresponded to high P r sites. Bright PSDs were selected based on their GluRIIA intensity (two standard deviations above average) and line profiles were drawn across each PSD. The intensity of pixels along that line for each fluorophore was then analyzed. Average pixel intensity revealed drastically distinct profiles for GluRIIB distribution between 'bright' and 'dim' PSDs classified based on their GluRIIA intensity. GluRIIB was more evenly distributed across the entire PSD at dim GluRIIA sites, but was segregated outward, forming a ring around central GluRIIA puncta at bright GluRIIA sites ( Figure 7F). In addition, presynaptic BRP intensity was more strongly correlated with postsynaptic GluRIIA levels (Pearson r = 0.53, R 2 = 0.28, p<0.0001, n = 2496 AZs from 19 NMJs from seven animals; Figure 7G) than with GluRIIB (Pearson r = 0.24, R 2 = 0.05, p<0.0001, n = 2496 AZs from 19 NMJs from seven animals; Figure 7H). These findings suggest that the postsynaptic cell accumulates GluRIIA and redistributes GluRIIB to the PSD periphery at high P r sites. Intravital imaging of glutamate receptors and cac channels throughout synapse development The Drosophila larval NMJ is a highly dynamic structure, with new synaptic boutons and AZs undergoing continuous addition throughout development (Harris and Littleton, 2015;Rasse et al., 2005;Schuster et al., 1996;Zito et al., 1999). Given the correlation between Ca 2+ channel abundance, GluRIIA/B segregation and high P r , we were interested in determining how AZs acquire a specific P r during a larval developmental period that lasts 6-7 days. One model is that certain AZs gain a higher P r status during development through preferential accumulation of key AZ components compared to their neighbors. Alternatively, high P r AZs might simply be more mature than their low P r neighbors, having an earlier birthdate and a longer timeframe to accumulate AZ material. To test these models for release heterogeneity, it would be desirable to follow P r development from the embryonic through larval stages. However, this is not technically feasible due to the small size of AZs and the rapid locomotion that larvae undergo, preventing generation of P r maps in moving animals. Instead, we employed an alternative intravital approach to repeatedly image the same NMJ at muscle 26 directly through the cuticle of intact larvae during anesthesia (Andlauer and Sigrist, 2012;Fouquet et al., 2009;Füger et al., 2007;Rasse et al., 2005;Zhang et al., 2010). During anesthesia, endogenous action potential-induced release and the associated GCaMP signals were eliminated, preventing direct P r measurements in anesthetized larvae. We instead focused on imaging GluRIIA accumulation and GluRIIA/B segregation, which was strongly correlated with P r (Figure 7), as a proxy for the emergence of high P r sites.
Previously described in vivo imaging approaches with anesthesia at Drosophila NMJs employed early third instar larvae as the starting time point, and followed the distribution of fluorescentlylabeled synaptic proteins during the final~36 hr of development prior to pupation (Andlauer and Sigrist, 2012;Fouquet et al., 2009;Füger et al., 2007;Rasse et al., 2005;Zhang et al., 2010). To follow AZ P r development beginning soon after synapse formation, we modified these techniques to allow imaging of glutamate receptors at earlier stages of development (see Materials and methods). This allowed successful birth dating and successive imaging of the same AZ over a 6 day period beginning shortly after synapse formation in the early first instar period through the late third instar stage (Figure 7-figure supplement 1). In early first instar larvae (within 12 hr of hatching) GluRIIA and GluRIIB were largely co-localized at postsynaptic puncta (Figure 7-figure supplement 1). One exception was the presence of diffuse GluRIIA that accumulated around unusually long axonal extensions that emerged from presynaptic boutons (Figure 7-figure supplement 2). These structures were devoid of any detectable GluRIIB or the bright GluRIIA puncta that are associated with AZs, and may be remnants of previously described muscle filopodial structures, termed myopodia, that interact with presynaptic filopodia to dynamically regulate early synaptic target recognition (Kohsaka and Nose, 2009;Ritzenthaler et al., 2000;Ritzenthaler and Chiba, 2003). GluRIIA was diffusely present on these structures, as has been observed for the leucine-rich repeat cell adhesion protein Capricious (Kohsaka and Nose, 2009). Repeated imaging of these thinner GluRIIA-positive processes revealed that they were capable of developing into mature synaptic boutons with concentrated GluRIIA and GluRIIB synaptic puncta (Figure 7-figure supplement 3). By 24 hr of larval growth, GluRIIA rich extensions were no longer observed, indicating these structures are restricted to early developmental stages.
Prior studies indicated GluRIIA PSD levels closely track with Cac accumulation at corresponding AZs of third instar larvae (Fouquet et al., 2009;Rasse et al., 2005), indicating the two compartments are likely to mature at similar rates. To examine this directly, we assayed whether GluRIIA and Cac accumulation were correlated during earlier stages of development (Figure 7-figure supplement 3). Indeed, the intensity of Cac-GFP and GluRIIA-RFP puncta were strongly correlated at individual AZs during both early and late larval development (first instar; Pearson r = 0.82, R 2 = 0.6771, p<0.0001, n = 441 AZs from 8 NMJs from eight animals; third instar; Pearson r = 0.63, R 2 = 0.395, n = 874 AZs from 8 NMJs from five animals; Figure 7-figure supplement 3D,E). One exception was observed in very early first instar larvae, where a few Cac-GFP puncta accumulated along the previously described GluRIIA positive axonal extensions prior to the specific accumulation of GluRIIA at PSDs (Figure 7-figure supplement 3A). In contrast, postsynaptic appearance of GluRIIA at the PSD could be observed to slightly precede the accumulation of Cac-GFP at later developmental stages (Figure 7-figure supplement 3B,C), similar to previous observations at mature third instar NMJs . In summary, the intensity of Cac-GFP and GluRIIA-RFP puncta are strongly correlated at AZs during larval development, indicating that GluRIIA provides a robust marker that reflects the corresponding levels of presynaptic Cac at individual AZs.

Analysis of Pr acquisition during AZ development using glutamate receptor segregation as a proxy
We examined how GluRIIA accumulation and GluRIIA/B segregation emerged during development of the NMJ. Live imaging of GluRIIA and GluRIIB distribution at early PSDs in anesthetized first instar larvae demonstrated that the receptors were co-localized and lacked the segregation where GluRIIB clustered around central GluRIIA puncta that was observed at high P r sites later in development ( Figure 8A). The first emergence of GluRIIA/B segregation was observed after 36 hr of imaging starting from the first instar period. The GluRIIA/B segregation always emerged first at the oldest synapses that existed previously on the first day of imaging ( Figure 8A,B). The most mature PSDs also contained more GluRIIA fluorescent signal (17430 ± 634.0, n = 86 AZs from 8 NMJs from five animals) compared to younger synapses that emerged during the 48 hr imaging session (8909 ± 289.8, n = 210 AZs from 8 NMJs from five animals). During later larval development, the cuticle thickness changed dramatically and prevented reliable comparison of absolute receptor density with earlier stages. However, GluRIIA intensities that were more uniform at the first instar larval stage became more heterogeneous at the third instar stage (Figure 8-figure supplement 1). Indeed, histograms of normalized fluorescence intensity (relative intensity scaled from 0 to 1) revealed that GluRIIA and GluRIIB were relatively uniformly distributed at first instar larval PSDs, with GluRIIA distribution becoming more skewed at later stages ( Figure 9A,B). To determine whether muscle 26 exhibits P r heterogeneity similar to muscle 4, we performed P r mapping using GluRIIA-RFP and myr-GCaMP6s in dissected non-anesthetized third larvae and confirmed that P r heterogeneity is extremely similar between muscle 26 (mean P r = 0.068 ± 0.0048, skewness = 2.27, 365 AZs from 5 NMJs from five animals) and muscle 4 (mean P r = 0.073 ± 0.0021, skewness = 2.23, 1933 AZs from 16 NMJs from 16 animals; Figure 9-figure supplement 1B-D). These results indicate that GluRIIA/B fluorescence and P r distribution are both highly heterogeneous by the early third instar stage in muscle 26, with the brightest GluRIIA PSDs, and by extension their corresponding high P r AZs, representing those that appeared earliest in development.
Over what time frame do synapses developmentally acquire markers of high P r sites? To estimate the average time required for AZ maturation, we calculated the time interval from the first emergence of a PSD in an imaging session to the time point when segregation of GluRIIB around GluRIIA central puncta occurred. This analysis was restricted to newly formed AZs that appeared during the imaging sessions and excluded AZs that were present in the first imaging session performed in first instar larvae. The average time from the first emergence of a PSD to when it acquired the segregated GluRIIA/B pattern observed at high P r AZs was 3.20 ± 0.08 days (n = 41 AZs from 7 NMJs from three animals; Figure 9C). In a small subset of PSDs (5%), a slightly faster accumulation of GluR-IIA and the formation of GluRIIB peripheral rings was observed, but never faster than 2 days. In addition, glutamate receptor fields increased in size throughout the first 72 hr of development ( Figure 9D). GluRIIA diameter increased by 1.47-fold during the first 24 hr of development but then plateaued, while GluRIIB field diameter continued to grow over 72 hr, increasing in size by 2.8-fold over that time. At 72 hr after initial glutamate receptor field formation, the average GluRIIA field diameter was 0.59 ± 0.058 mm (range 0.45-0.80 mm, 12 PSDs from three animals) and the average GluRIIB field diameter was 1.01 ± 0.035 mm (range 0.89-1.12 mm, 12 PSDs from three animals), consistent with the formation of segregated GluRIIA/B fields at mature AZs.
To directly assess whether AZ age corresponds to P r , we followed PSDs in animals expressing GluRIIA-RFP, GluRIIB-GFP and myr-GCaMP for 24 hr using in vivo imaging, and then dissected the animals and mapped P r directly ( Figure 9F). We observed that PSDs that emerged on the second day of imaging (less than 24 hr old) were consistently associated with very low P r AZs. The mean P r of these newly-formed AZs was 0.035 ± 0.0047 (n = 77 AZs from 4 NMJs from four animals), with a range of 0-0.14. In contrast, the average P r for AZs older than 24 hr was 0.14 ± 0.0093 with a range of 0-0.61 (n = 188 AZs from 4 NMJs from four animals; Figure 9G). These findings support a model in which the vast majority of newly formed AZs are very weak, with increased P r requiring more than 24 hr to develop.
Given the developing NMJ is adding AZs at a rapid rate Schuster et al., 1996), we estimated whether AZ maturation time identified over the course of our live imaging experiments could lead to the~10% of high P r sites observed at the early third instar stage. We quantified the number of synapses present at the same NMJ from the first instar through the early third instar stage from live imaging experiments ( Figure 9E). AZ number roughly doubled each day, such that the average number of AZs found at the first instar stage (day 1) represented 14.7 ± 1.4% (n = 8 NMJs from three animals) of all AZs present by day 4 (3 days after initial imaging in first  (0) to brightest (1). GluRIIA shows a more skewed distribution of fluorescence at day 4, consistent with its accumulation at high P r AZs. (C) Representative muscle 26 NMJ image sequence showing appearance and maturation of two Figure 9 continued on next page instars). Overall, these data are consistent with the hypothesis that AZ maturation is a key factor in regulating P r , leading to increased accumulation of Ca 2+ channels and GluRIIA/B segregation at high P r sites compared to AZs that are newly formed (<2 days).
If heterogeneity in AZ age underlies the majority of diversity in P r , we would expect to see a different P r distribution at the earlier second instar stage when there is a smaller range of AZ ages. Since the number of AZs nearly doubles each day and newly formed AZs display low P r , we hypothesized that third instar NMJs would have a greater proportion of low P r sites compared to second instar NMJs. To test this hypothesis, we mapped P r at muscle 4 NMJs in the second instar and compared the distribution to that seen at the third instar stage (Figure 9-figure supplement 1A,B). A rightward shift in the distribution of P r was observed in the earlier second instar stage, with a greater proportion of AZs in the higher P r category, a smaller population of low P r AZs, and a significant increase in mean P r in second instars compared to third instars (third instar: 0.07 ± 0.002, n = 1933 AZs from 16 NMJs from 16 animals; second instar: 0.13 ± 0.008, n = 282 AZs from 6 NMJs from six animals; Figure 9-figure supplement 1C,D). We were unable to map P r at muscle 26 in earlier stages in dissected animals because it is covered by muscles 6 and 7 until the third instar stage. These data support the hypothesis that P r heterogeneity reflects AZ age and maturation time.

Postsynaptic maturation rate depends on presynaptic activity
To determine whether the developmental time-course of synapse maturation could be modulated by changes in presynaptic activity, we measured the rate of PSD growth (fold-increase in GluRIIB area over 24 hr) and the percent of PSDs displaying GluRIIA/B segregation in mutants with altered presynaptic activity ( Figure 10). We first measured postsynaptic maturation in BRP 69/def null animals, which have a dramatic reduction in evoked synaptic transmission . Consistent with previous findings, Cac-GFP intensity in BRP 69/def animals was reduced to 25% (mean fluorescence intensity 3549 ± 23, 4 NMJs from two animals) of control levels (mean intensity 14177 ± 220, 4 NMJs from four animals; Figure 10-figure supplement 1A,B). We expressed GluRIIA-RFP and GluRIIB-GFP in BRP 69/def mutants and imaged muscle 26 NMJs in anesthetized larvae intravitally over 24 hr. A significant reduction in postsynaptic maturation rate was observed; newly formed BRP 69/def PSDs only increased in GluRIIB area by 1.16-fold (±0.11, 4 NMJs from three animals) over the first 24 hr of development, compared to a 1.61-fold (±0.11, 5 NMJs from five animals) increase in controls ( Figure 10A,C). Furthermore, a significant reduction in the percent of PSDs with GluRIIA/B rings (defined by a > 10% central dip in the GluRIIB intensity profile) was observed in the second day of larval development; only 4.9% of BRP 69/def PSDs (±1.3%, 10 NMJs from five animals) showed receptor segregation compared to 22% (±5%, 11 NMJs from five animals) in age-matched controls ( Figure 10D). We were unable to observe any GluRIIB rings at the third instar stage in BRP 69/def , in contrast to the clear rings seen in controls at this stage. Instead, BRP 69/def mutants displayed highly disorganized GluR clusters ( Figure 10B). new synapses (#3 and 4) that were not present in the initial imaging session. Several preexisting synapses (#1 and 2) that developed the typical GluRIIB donut structure later in development are also labeled. The dashed box surrounds the actual imaged segment of the NMJ. New GluRIIA and GluRIIB clusters appear initially as small puncta (day 2, arrows) that become brighter on day 3. By day four they begin to display the donut like GluRIIB profile. At day 5, GluRIIB distribution to the periphery around a bright GluRIIA PSD representative of high P r sites becomes prominent. The bottom panels show normalized GluRIIA and GluRIIB fluorescent line profiles for the newly identified PSDs (#3, 4) throughout the 5 day imaging series. (D) Diameter of glutamate receptor fields during the first 72 hr of PSD development. Error bars represent SEM. (E) Changes in AZ number during larval maturation at muscle 26 presented as a ratio of AZs observed during the first day of imaging (day 1). (F) Representative serial time points of GluRIIA-RFP in vivo imaging over 24 hr (left two panels). Newly formed PSDs are marked with white arrows; these displayed uniformly low P r during release mapping (right panel). Several older PSDs with bright GluRIIA intensity are denoted with white numbers. (G) Average P r of old AZs (those present in the 0 hr time point initial imaging session) and new AZs (AZs under 24 hr old that were first seen at the 24 hr time point). Student's t-test was used for statistical analysis (***=p 0.001). Error bars represent SEM. DOI: https://doi.org/10.7554/eLife.38268.022 The following figure supplement is available for figure 9: To examine the consequences of reductions in presynaptic activity independent of structural alterations in the AZ that might occur at NMJs in BRP mutants, we measured PSD maturation rate in nap TS and syt1 null mutants that reduce release through separate mechanisms. Loss of Syt1 causes a dramatic decrease in synchronous fusion (DiAntonio and Schwarz, 1994;Guan et al., 2017;Lee et al., 2013;Littleton et al., 1994Littleton et al., , 1993Yoshihara et al., 2003;Yoshihara and Littleton, 2002), while nap TS results in constitutively reduced neuronal excitability due to decreased sodium channel activity (Kernan et al., 1991;Wu et al., 1978). In both syt1 and nap TS mutants, significantly reduced PSD growth rate and GluRIIA/B segregation was observed compared to controls; GluRIIB field area in nap TS and syt1 null mutants only increased by 1.24-fold (±0.05 4 NMJs from three animals) and 1.15-fold (±0.05, 4 NMJs from four animals) respectively during the first 24 hr of PSD development compared to 1.61-fold (±0.11, 5 NMJs from five animals) in controls ( Figure 10C). At the second instar stage, when 22% of control PSDs displayed GluRIIA/B segregation, only 5.4% of syt1 (±1.5%, 7 NMJs from three animals) and 6.0% of nap TS (±2.0%, 13 NMJs from four animals) PSDs showed GluR segregation ( Figure 10D). Both nap TS and syt1 null mutants also formed less segregated GluRIIA/B fields compared to controls at the third instar stage ( Figure 10B). To determine whether increasing presynaptic release would accelerate maturation, we examined PSD size and GluR segregation in shaker 120 , eag 1 (sh, eag) double mutants that display increased excitability due to loss of several K + currents Salkoff and Wyman, 1981;Wu et al., 1983). Sh, eag showed a significantly increased rate of GluRIIB field size increase (2.37-fold ± 0.11, 5 NMJs from four animals) over 24 hr compared to controls (1.61 fold (±0.11, 5 NMJs from five animals; Figure 10A,C). Furthermore, 48.7% of sh, eag mutant PSDs (±4.4%, 13 NMJs from five animals) displayed GluRIIB rings at the early second instar stage, compared to only 22% of control PSDs ( Figure 10D).
These observations indicate that increasing (sh, eag) or decreasing (nap TS , BRP 69/def , syt1) presynaptic activity results in a respective increase or decrease in PSD maturation rate. One factor that might contribute to these changes in PSD maturation rate is altered growth rate of the entire NMJ secondary to changes in neuronal activity. Indeed, we observed that the rate of new PSD addition over 24 hr in BRP 69/def , nap TS , and sh, eag mutants showed a similar trend to PSD maturation, suggesting that overall NMJ growth rate is altered in these mutants (control animals increase AZ number by 1.8 ± 0.1 fold over 24 hr compared to sh, eag (2.2 ± 0.1 fold), nap TS (1.3 ± 0.04 fold), syt1 (1.5 ± 0.08 fold), and BRP 69/def (1.2 ± 0.06 fold); Figure 10E). Previous work on activity-dependent NMJ growth in Drosophila focused primarily on bouton number, not AZ addition rate, and it will be interesting to explore the molecular mechanisms mediating activity-dependent seeding of new AZs. To measure PSD maturation rate in the context of an otherwise healthy animal, we took advantage of the rab3 null mutant, where roughly half of the PSDs are apposed by presynaptic areas lacking BRP and Cac, and the remaining AZs contain higher than normal levels of these presynaptic components (Graf et al., 2009). Rab3 AZs containing BRP and Cac have been shown to have higher P r , whereas AZs lacking these presynaptic components are functionally silent (Peled and Isacoff, 2011). This redistribution of P r provides an opportunity to compare PSD maturation with less influence of overall animal health and NMJ growth rate.
The rab3 phenotype has been previously analyzed at the third instar stage of larval development. To characterize the dynamics of pre-and postsynaptic maturation during earlier development, we followed rab3 animals expressing Cac-GFP and GluRIIA-RFP for 48 hr starting in the early first instar larval stage. AZs populated with Cac-GFP early in development were stable over time; Cac-GFP was never lost at these AZs over 48 hr of imaging ( Figure 10-figure supplement 1E). However, a divergence in the rate of PSD versus Cac-GFP addition was observed after the first instar stage; new GluRIIA-containing PSDs were added to the NMJ at a rate similar to control, with total PSD number nearly doubling every 24 hr. However, the rate of addition of new Cac-GFP-containing AZs was much lower than in controls, resulting in a decrease in Cac/GluRIIA apposition in rab3 mutants over development ( Figure 10-figure supplement 1D,F). This is consistent with the observation that PSD growth rate was decreased in the rab3 null (1.23-fold ± 0.05, 6 NMJs from six animals) when we measured PSDs born in the late first instar or early second instar stage ( Figure 10C), while the percent of GluRIIA/B rings seen in the second instar stage was elevated (34.7 ± 3.3% of PSDs had rings, 14 NMJs from six animals) compared to age matched controls (22%; Figure 10D). To avoid complications from the atypical dynamics of AZ apposition during rab3 NMJ development, we focused the remainder of our analysis on the early first instar stage where the age distribution of AZs is much narrower.
We next examined first instar rab3 mutants to determine whether AZs with higher presynaptic activity had more mature postsynaptic receptor fields than their silent neighbors. We imaged GluR-IIA-RFP and GluRIIB-GFP in rab3 1 st instars and then dissected and stained for BRP. In muscle 26 of rab3 mutants, roughly 50% of first instar AZs were populated with BRP (the ratio of GluRIIA-RFP puncta to Cac-GFP puncta was 1.9 ± 0.15); these AZs were opposed by large PSDs, many of which already showed GluRIIA/IIB segregation ( Figure 10D, Figure 10-figure supplement 1C). In contrast, PSDs lacking BRP were much smaller and lacked receptor segregation. All of the PSDs with GluRIIB rings were BRP-positive, suggesting that PSDs opposite AZs with reduced presynaptic release are less mature at this stage of larval development. Though the mechanism by which rab3 regulates AZ assembly remains unknown, the larger PSD size and enhanced GluRIIA/B segregation opposite BRP-positive AZs and the small size and lack of receptor segregation opposite BRP-negative AZs within the same terminal suggests that PSD maturation depends on presynaptic activity at the level of individual AZs.

Discussion
In the current study we used quantal imaging, super resolution SIM, and intravital imaging to examine the development of heterogeneity in evoked P r across the AZ population at Drosophila NMJs. We first confirmed that release heterogeneity was not caused by summation of fusion events from multiple unresolvable AZs. Indeed, high P r sites corresponded to single AZs with enhanced levels of BRP. These findings are consistent with previous observations using conventional light microscopy that indicate P r correlates with BRP levels (Muhammad et al., 2015;Peled et al., 2014;Reddy-Alla et al., 2017). By monitoring release over intervals of extensive vesicle fusion during strong stimulation, we also observed that P r is a stable feature of each AZ. In addition, loss of the synaptic vesicle Ca 2+ sensor Syt1 globally reduced P r without altering the heterogeneous distribution of P r across AZs, indicating that AZ-local synaptic vesicle pools with differential Ca 2+ sensitivity are not likely to account for P r heterogeneity. Since VGCC abundance, gating, and organization within the AZ are well established regulators of P r across synapses (Borst and Sakmann, 1996;Chen et al., 2015;Meinrenken et al., 2002Zito et al., 1999Nakamura et al., 2015;Sheng et al., 2012;Südhof, 2012;Wang et al., 2008), heterogeneity in presynaptic Ca 2+ channel abundance was a clear candidate for the generation of P r heterogeneity at Drosophila NMJs. Indeed, the Cac Ca 2+ channel responsible for neurotransmitter release is heterogeneously distributed across the NMJ (Fouquet et al., 2009;Kawasaki et al., 2004;Littleton and Ganetzky, 2000;Liu et al., 2011;Rieckhof et al., 2003;Smith et al., 1996). Using transgenically labeled Cac lines, we observed that Cac density at AZs is indeed strongly correlated with P r . To directly visualize presynaptic Ca 2+ influx at single AZs, we generated GCaMP fusions to the core AZ component BRP. Ca 2+ influx at single AZs was highly correlated with both Cac density and P r . The cac NT27 mutant with decreased conductance also resulted in a global reduction in P r without disrupting heterogeneity, further confirming that Ca 2+ influx regulates P r across the range of release heterogeneity.
Postsynaptically, high P r AZs were enriched in GluRIIA-containing receptors and displayed a distinct pattern of glutamate receptor clustering. While most synapses showed GluRIIA and GluRIIB spread over the entire PSD, high P r AZs were apposed by PSDs where GluRIIA concentrated at the center of the AZ, with GluRIIB forming a ring at the PSD periphery. Indeed, anti-glutamate receptor antibody staining of wildtype larvae lacking tagged glutamate receptors had previously identified a GluRIIB ring around the GluRIIA core in some mature third instar NMJ AZs (Marrus et al., 2004). In addition, activity-dependent segregation of GluRIIA and a GluRIIA gating mutant has been observed at individual AZs in Drosophila (Petzoldt et al., 2014). The correlation of P r with GluRIIA accumulation is especially intriguing considering that this subunit has been implicated in homeostatic and activity-dependent plasticity (Davis, 2006;Frank, 2014;Petersen et al., 1997;Sigrist et al., 2003). By following the developmental acquisition of this postsynaptic property as a proxy for P r from the first through third instar larval stages via intravital imaging in control and mutant backgrounds, we observed that the earliest formed AZs are the first to acquire this high P r signature over a time course of~3 days, and that PSD maturation rate can be modulated by changes in presynaptic activity.
Similar to prior observations (Melom et al., 2013;Peled and Isacoff, 2011), we found that most AZs at the Drosophila NMJ have a low P r . For the current study, the AZ pool was artificially segregated into low and high release sites, with high releasing sites defined based on having a release rate greater than two standard deviations above the mean. Given that birthdate is a key predictor of glutamate receptor segregation, and by proxy P r , we expect the AZ pool to actually reflect a continuum of P r values based on their developmental history. However, using the two standard deviation criteria, 9.9% of AZs fell into the high P r category, with an average P r of 0.28. We also observed that 9.7% of the AZs analyzed displayed only spontaneous release. We could detect no fusion events for either evoked or spontaneous release for another 14.6% of AZs that were defined by a GluRIIA-positive PSD in live imaging. Future investigation will be required to determine whether these cases represent immature AZs with extremely low evoked P r , or distinct categories reflective of differences in AZ content. The remaining AZs that participated in evoked release had an average P r of 0.05. Ca 2+ channel density and Ca 2+ influx at individual AZs was a key determinant of evoked P r heterogeneity, as P r and the intensity of Cac channels tagged with either TdTomato or GFP displayed a strong positive correlation. Spontaneous fusion showed a much weaker correlation with both Cac density and Ca 2+ influx at individual AZs, consistent with prior studies indicating spontaneous release rates are poorly correlated with external Ca 2+ levels at this synapse (Jorquera et al., 2012;Lee et al., 2013). With synaptic vesicle fusion showing a steep non-linear dependence upon external Ca 2+ with a slope of~3-4 (Dodge and Rahamimoff, 1967;Heidelberger et al., 1994;Jan and Jan, 1976), a robust change in P r could occur secondary to a relatively modest increase in Ca 2+ channel abundance over development. Although the number of VGCCs at a Drosophila NMJ AZ is unknown, estimates of Cac-GFP fluorescence during quantal imaging indicate a~2 fold increase in channel number would be necessary to move a low P r AZ into the high P r category. Similar correlations between evoked P r and Ca 2+ channel abundance have been found at mammalian synapses (Holderith et al., 2012;Nakamura et al., 2015;Sheng et al., 2012), suggesting this represents a common evolutionarily conserved mechanism for determining release strength at synapses.
We did not test the correlation of P r with other AZ proteins besides Cac and BRP, but it would not be surprising to see a positive correlation with the abundance of many AZ proteins based on the observation that maturation time is a key determinant for P r . Indeed, recent studies have begun to correlate P r with specific AZ proteins at Drosophila NMJs (Reddy-Alla et al., 2017). We also observed that PSD size was robustly increased by 1.6-fold over a 24 hr period of AZ development during the early larval period in control animals. AZ maturation is likely to promote increased synaptic vesicle docking and availability, consistent with observations that correlate AZ size with either P r or the readily releasable pool (Han et al., 2011;Holderith et al., 2012;Matkovic et al., 2013;Matz et al., 2010;Nakamura et al., 2015;Schikorski and Stevens, 1997;Wadel et al., 2007).
We considered several models for how AZs acquire this heterogeneous nature of P r distribution during a developmental period lasting several days. One possibility is that unique AZs gain high P r status through a mechanism that would result in preferential accumulation of key AZ components compared to their neighbors. Given that retrograde signaling from the muscle is known to drive synaptic development at Drosophila NMJs (Ball et al., 2010;Berke et al., 2013;Harris and Littleton, 2015;Keshishian and Kim, 2004;McCabe et al., 2003;Piccioli and Littleton, 2014;Yoshihara et al., 2005), certain AZ populations might have preferential access to specific signaling factors that would alter their P r state. Another model is that AZs compete for key presynaptic P r regulators through an activity-dependent process. High P r AZs might also be more mature than their low P r neighbors, having a longer timeframe to accumulate AZ components. Given the Drosophila NMJ is constantly forming new AZs at a rapid pace during development Schuster et al., 1996), newly formed AZs would be less mature compared to a smaller population of 'older' high P r AZs.
To examine if the release heterogeneity observed at the third instar stage reflects AZ birth order over several days of development, we needed to extend the intravital imaging through a longer time period beginning in the first instar larval stage. GCaMP imaging indicated high P r sites segregate GluRIIA and GluRIIB differently from low P r sites, with the IIA isoform preferentially localizing at the center of PSDs apposing high P r AZs. As such, we used developmental acquisition of this property as an indicator of high P r sites. Although segregation of glutamate receptors may not perfectly replicate the timing of P r acquisition during development, it is currently the best tool for estimating P r during sequential live imaging. Based on the acquisition of GluRIIA/B segregation, the data support the hypothesis that increases in P r reflect a time-dependent maturation process at the NMJ. The continuous addition of new AZs, which double in number during each day of development, ensures that the overall ratio of high to low P r sites represents a low percentage as the NMJ grows ( Figure 11). We further established confidence in the age-dependency model of P r by mapping release after following NMJs intravitally for 24 hr; using this approach, we found that newly formed AZs are consistently very low P r . Finally, we mapped P r in the second instar stage and observed that the heterogeneity at this stage is shifted towards higher releasing sites, with a reduction in the fraction of low-releasing sites.
These data indicate AZs are born with low P r and gain pre-and postsynaptic material over 3 days on an upward trajectory toward higher P r status ( Figure 11). Is AZ age a static determinant of P r , or can growth rate be regulated to allow faster or slower acquisition of high P r AZs? To answer this question, we assayed whether mutants that alter presynaptic activity influence the rate of PSD growth and GluRIIA/B segregation. In BRP 69/def , syt1 null , and nap TS mutants with decreased Figure 11. Model depicting the acquisition of release heterogeneity over three days of NMJ development. (A) Presynaptically, Ca 2+ channels (yellow) and BRP (blue) are added over time to increase release probability. Postsynaptically, GluRIIA and GluRIIB are added during maturation and eventually segregate at high releasing sites. (B) The number of AZs roughly doubles every 24 hr, with all new AZs having low release probability (blue dot). AZs increase in release probability over time and eventually reach high P r status (red dot). DOI: https://doi.org/10.7554/eLife.38268.026 presynaptic release, a significant reduction in postsynaptic maturation rate and in the percentage of PSDs with GluRIIA/IIB rings was observed. Conversely, in the shaker, eag double mutant with increased presynaptic excitability, a significantly increased rate of GluRIIB field size and a significant increase in GluRIIA/IIB rings was found compared to controls. To investigate whether differences in synapse growth and maturation rate could be seen between AZs enriched in BRP and Cac versus neighboring AZs that are deficient in these components, we imaged development in the rab3 null mutant. At the first instar stage when AZs are roughly age matched, release sites enriched in presynaptic components had developed large and mature PSDs in stark contrast with non-enriched AZs, whose PSDs appeared immature. These results indicate that PSD maturation can be influenced by presynaptic activity at the resolution of single AZs.
Although how AZs are assembled during development is still being established, our data do not support a model where AZs are fully preassembled during transport and then deposited as a single 'quantal' entity onto the presynaptic membrane. Rather, these data support a model of seeding of AZ material that increases developmentally over time as AZs matures, consistent with previous studies of AZ development in Drosophila (Bö hme et al., 2016;Fouquet et al., 2009). Although no evidence for rapid changes in P r were detected in the steady-state conditions used in the current study, homeostatic plasticity is known to alter P r over a rapid time frame (~10 min) at the NMJ (Davis and Müller, 2015;Frank, 2014;Frank et al., 2006). It will be interesting to determine if Cac abundance can change over such a rapid window, or whether the enhanced release is mediated solely through changes in Cac function and Ca 2+ influx (Müller and Davis, 2012). Changes in the temporal order of P r development could also occur secondary to altered transport or capture of AZ material. For example, the large NMJ on muscle fibers 6 and 7 displays a gradient in synaptic transmission, with terminal branch boutons often showing a larger population of higher P r AZs (Guerrero et al., 2005;Peled and Isacoff, 2011). If AZ material is not captured by earlier synapses along the arbor, it would be predicted to accumulate in terminal boutons, potentially allowing these AZs greater access to key components, and subsequently increasing their rate of P r acquisition. Alternatively, the gradient of P r along the axon could be due to terminal boutons being slightly older than the rest of the arbor.
In summary, our data indicate that heterogeneity in release correlates highly with Ca 2+ channel abundance and Ca 2+ influx at AZs. Postsynaptically, PSDs apposed to high-releasing AZs display increased GluRIIA abundance and form segregated receptor fields, with GluRIIB forming a ring around a central core of GluRIIA. Release sites accumulate these high P r markers during a synapse maturation process in which newly formed AZs are consistently low P r , with AZs gaining signatures of high releasing sites over several days. Finally, mutations that increase or decrease presynaptic activity result in faster or slower rates of PSD maturation, respectively. These data add to our understanding of the molecular and developmental features associated with high versus low P r AZs.

Transgenic constructs
The fluorescent Ca 2+ sensor GCaMP6s was tethered to the plasma membrane with an N-terminal myristoylation (myr) sequence. The first 90 amino acids of Src64b, containing a myristoylation target sequence, were subcloned into pBID-UASc with EcoRI and BglII (creating pBID-UASc-myr). GCaMP6s cDNA (Addgene plasmid 40753) was cloned into pBID-UASc-myr with BglII and XbaI. To generate the UAS-GCaMP6m-Brp-short line, GCaMP6m (Addgene plasmid 40754) cDNA and Brpshort (gift from Dr. Tobias Rasse) were PCR amplified and double digested with EcoRI/BglII and BglII/XbaI, respectively. The two cDNA fragments were ligated and the product was used to PCR amplify the fused GCaMP6m-Brp-short cDNA. The PCR product was inserted into the vector backbone pBID-UASc after digestion with EcoRI and XbaI to generate the final plasmid pBID-UASc-GCaMP6m-Brp-short. To create UAS-myrjRGECO, the vector backbone pBID-UASc-myr was digested with BglII and XbaI. jRGECO sequence was amplified from plasmid pGP-CMV-NES-jRGE-CO1a (gift from Dr. Douglas Kim, Addgene plasmid #61563). The digested backbone and insert were fused according to the Gibson assembly protocol using NEBuilder HighFidelity DNA Assembly Cloning Kit (E5520). To generate pBid-lexAop-myrGcaMP6s, myrGCaMP6s was amplified by PCR and inserted into pBiD-lexAop-DSCP (gift from Brian McCabe) between NotI and XbaI sites. All transgenic Drosophila strains were generated by BestGene.

Confocal imaging and data analysis
Confocal images were obtained on a Zeiss Axio Imager two equipped with a spinning-disk confocal head (CSU-X1; Yokagawa) and ImagEM X2 EM-CCD camera (Hammamatsu). An Olympus LUMFL N 60X objective with a 1.10 NA was used to acquire GCaMP6s imaging data at 7 to 8 Hz. A Zeiss pan-APOCHROMAT 63X objective with 1.40 NA was used for imaging stained or live animals. third instar larvae were dissected in Ca 2+ -free HL3 containing 20 mM MgCl 2 . After dissection, preparations were maintained in HL3 with 20 mM MgCl 2 and 1.3 mM Ca 2+ for 5 min. To stimulate the NMJ, motor nerves were cut close to the ventral ganglion and sucked into a pipette. Single pulses of current were delivered every one second for myr-jRGECO mapping or every three seconds for GCaMP6s mapping with an AMPI Master-8 stimulator using a stimulus strength just above the threshold for evoking EJPs. A 3D image stack was taken before the GCaMP imaging session to generate a full map of GluRIIA or Cac channel distribution. Later, single focal planes were imaged continuously for 4-5 min to collect GCaMP signals. Volocity 3D Image Analysis software (PerkinElmer) was used to analyze images. All images were Gaussian filtered (fine) to reduce noise and a movement-correction algorithm was applied. To enhance identification of myrGCaMP6 flashes, background myrGCaMP fluorescence was subtracted by creating a composite stack of 5-6 images during intervals when no synaptic release occurred. To identify the position of GluRIIA receptors and corresponding Ca 2+ events, a 3D stack image of GluRIIA was merged to create a single plane. AZ position was identified using the 'find spot' algorithm in Volocity 3.2 software that detects fluorescent peaks. ROIs with identical 5-pixel size (0.138 mm/pixel) were automatically generated by the software from identified GluRIIA spots. All GCaMP flashes were detected using the intensity threshold tool and assigned to specific ROIs based on proximity of their centroids. The time and location of Ca 2+ events were imported into Excel or Matlab for further analysis. The number of observed GCaMP events per AZ was divided by the number of delivered stimuli to calculate AZ P r . Analysis of Cac, BRP, GluRIIA or GluRIIB intensities was performed similarly, identifying AZ fluorescence peaks and defining three pixel square ROIs around each peak to calculate average fluorescence. Average AZ fluorescence intensities of 3-pixel square ROIs was also used for correlation analysis.

SIM and airyscan imaging
SIM microscopy was performed on an Applied Precision DeltaVision-OMX BLAZE-3D-Structural Illumination Microscope equipped with 60X, 1.4 NA oil objective and 3 sCMOS cameras. 3D-SIM images were obtained with 125 nm z-steps. Samples were illuminated by three central diffraction orders with 488, 562, and 640 nm lasers. For initial identification of specific NMJs, larvae were imaged in conventional confocal mode using a 20X oil objective. The positions of NMJ were marked and recorded to provide transition between objectives. A ZEISS LSM 800 microscope with Airyscan was also used to image anesthetized animals. Fluorescence was detected by a concentricallyarranged hexagonal detector array consisting of 32 single detector elements.

Live imaging
Larvae were anesthetized with SUPRANE (desflurane, USP) from Amerinet Choice (Zhang et al., 2010). Larvae were incubated in a petri dish with a small paper towel containing Suprane for 1-2 min in a fume hood. Anesthetized larvae were positioned ventral side up on a glass slide between spacers made by transparent tape, which prevented extreme compression of the larvae. Different size spacers were required for the various larval stages. Larvae were covered with a thin film of halocarbon oil and then with a cover glass. NMJ synapses on muscle 26 in hemi-segment 2 or three were imaged. After an imaging session, larvae were placed in numbered chambers with food in a 25˚C incubator. The same data acquisition settings where used to visualize NMJs at different larval stages. Larvae were imaged with either 6, 24 and 36 hr intervals for one data set ( Figure 8A-C), or for 24 hr intervals for the remaining datasets. To keep the size consistent between different time periods, images of the corresponding NMJ area at younger stages were cut (dashed areas in figures) and placed onto a black background. This presentation generated a similar orientation of the different size NMJs for easier comparison for Figure 8, Figure 9 and

Ionomycin application
Ionomycin (Sigma Aldrich) was dissolved in ethanol to make a 10 mM stock solution and was stored at 4˚C. Ionomycin was used at a working concentration of 200 nM dissolved in HL3 with 1.3 mM Ca 2+ ; this solution was applied to dissected preparations and NMJs were imaged one minute after application.

Statistical analysis
Statistical analysis was performed with GraphPad Prism using one-way ANOVA followed by Dunnett's Multiple Comparisons test for comparison of samples within an experimental group, or Student's t-test for comparing two groups. Asterisks denote p values of: *p 0.05; **p 0.01; and ***p 0.001. All histograms and measurements are shown as mean ± SEM. Pearson coefficient of correlation was calculated in GraphPad Prism using the following parameters: -two-tailed P value and 95% confidence interval.