In this work we examine counts of SV release events measured in simple PF-MLI synapses when increasing the release probability with several methods. Our aim is to identify the source(s) of the increase of release probability in each case, taking advantage of the detailed information provided by SV counting in individual AZs. As further discussed below, we assume for the interpretation of our results that release occurs at a finite set of docking/release sites, and that SVs undergo exocytosis only once they are docked. Within this simple docking/release model, we attempt to distinguish increased probability of docking site occupancy by SVs from increased probability of exocytosis of an occupied docking site (Quastel, 1997).

Figure 1A illustrates the experimental protocol. Having established a recording in an MLI, a pipette was placed either in the granule cell layer or in the molecular layer to stimulate extracellularly a presynaptic granule cell or its associated PF axon. The pipette position and stimulation strength were carefully adjusted to restrict the stimulation to a single presynaptic PF (see Materials and methods; Malagon et al., 2016; Miki et al., 2017). Trains of 8 action potentials (APs) were applied at a frequency of 100 Hz inside a train, leaving 10 s intervals between trains. Typical traces in response to an 8-AP stimulation train (stimulation times: vertical red bars) are shown in Figure 1A. As illustrated in this example, responses to individual APs varied between single quantal responses (with unitary amplitudes on the order of 100 pA or more in this preparation: Llano and Gerschenfeld, 1993; examples in responses to 1 st and 2nd AP in 1 st trace), failures (as following 3rd and 4th AP of 1 st trace), and multivesicular responses (eg, following the 5th stimulation in the 1 st trace). Only results passing previously defined criteria for simple synapse recordings were accepted in this work (Malagon et al., 2016). Briefly, the criteria were: (i) stability of the recording as a function of time; (ii) amplitude occlusion for consecutive events at short intervals; (iii) low scatter of quantal current amplitudes (Malagon et al., 2016). In recordings that met these criteria, we identified single release events during AP trains by deconvolution, using a synapse-specific averaged quantal EPSC profile as template (Malagon et al., 2016). We then summed the numbers of SVs released in a 5 ms long time window following each AP stimulation. These numbers represent the output of individual AZs following AP stimulations. They are depicted in matrices as a function of AP number and train number in Figure 1B, where grey levels represent various SV numbers.

Figure 1 with 1 supplement see all

Download asset Open asset

A first set of experiments was performed under normal release conditions, using an external Ca²⁺ concentration ([Ca²⁺]_o) of 1.5 mM. We then observed a gradual increase in SV counts during a train, reflecting facilitation (Figure 1B, upper panel; group data in Figure 1D, pink). After obtaining a series of train responses, SV counts were calculated separately for each of the 8 successive AP stimulations, and we calculated means and variances of SV counts according to the AP stimulation number, i. This gave a variance vs. mean plot with 8 data points that could be fitted using a parabola (Figure 1C). To interpret these data we assumed that each simple synapse obeys a binomial release process, with N functional units (docking/release sites) acting in parallel (Malagon et al., 2016). The intersect of the parabola with the abscissa axis yielded the value of N. For stimulus number i, we calculated the probability P_i that a docking site releases a SV by dividing the mean SV count by N, with a mean N value of 3.6 across synapses. The mean of P₁ obtained in this manner under control conditions was 0.089 ± 0.015 while the mean of the maximum of P_i for i = 1–8 was P_max = 0.270 ± 0.027. These data show that facilitation markedly increases P under control conditions, with a ratio between P₁ and P_max of 3-fold. They also indicate that P is far from reaching the maximum value of 1 under these conditions.

Standard docking/release site models assume that a docking site can release only if it is occupied by a SV. Accordingly, these models predict that P_i is the product of δ_i and p_i, where δ_i is the mean docking site occupancy (sometimes noted p_occ) before stimulus number i, and p_i is the probability of release of a docked SV (sometimes noted p_ves) following stimulus number i (Vere-Jones, 1966; Zucker, 1973; Quastel, 1997; Scheuss and Neher, 2001). In the following we will accept the P_i = δ_i p_i relation as a guideline. Recent results calling for a qualification of this hypothesis (Miki et al., 2018) will be discussed below.

Increasing evidence indicates that both facilitation and depression mostly follow changes in docking site occupancy, reflecting a trade-off between SV loss by exocytosis and SV gain by docking site replenishment (Pulido et al., 2015; Miki et al., 2016; Miki et al., 2018). Accordingly, we assume in the following that p_i takes a common value, p, independently of i, so that P_i = δ_i p. For i = 1 we write δ₁ = δ, and P₁ = δ p, where δ is the resting docking site occupancy. In line with our previous work (Miki et al., 2016; Miki et al., 2017), we further assume that all release sites have the same release probability and that they are located at equal distances from presynaptic voltage-dependent calcium channels.

While the value of P₁ is readily accessible by performing variance-mean analysis of SV counts, as explained above, separating p from δ poses a challenge. An attractive way to determine δ is to increase the amount of Ca²⁺ entry until p reaches one because in this case, the relation P₁ = δ p becomes P₁ = δ. Therefore, measuring P₁ under very high release probability conditions provides a direct estimate of δ.

We therefore explored ways to increase p in the present work. Previous studies indicate that addition of TEA increases p in many synapses including PF-MLI synapses by broadening AP duration and increasing the amount of Ca²⁺ entry per AP, mainly by blocking Kv3 channels (Sabatini and Regehr, 1997; Ishikawa et al., 2003). When applying TEA in granule cell recordings, we found in agreement with previous work (Sabatini and Regehr, 1997) a decrease in voltage-dependent K⁺ current amplitudes and an increase in AP duration (Figure 1—figure supplement 1). Furthermore, TEA increased the responses of PF-MLI synapses to AP stimulations (Figure 1B: compare first and second panel from above). Comparing released SV counts before and after the addition of TEA in the same experiment, as illustrated in the example of Figure 1B, shows a marked increase in the presence of TEA. Average SV counts across experiments for the first AP were n₁ = 0.36 ± 0.06 in control and n₁ = 0.85 ± 0.26 in TEA (n = 6; p<0.05, paired t-test). Variance-mean analysis of SV counts (Figure 1C, first and second graphs from top) showed no change of N upon application of TEA (individual experiments and corresponding means in Figure 1G). Given that N stays constant, the change in synaptic strength associated with TEA application is entirely mediated by an increase in P. After normalizing SV counts per docking site, we found that P₁ values increased from 0.103 to 0.213 in 1 mM TEA (Figure 1E).

Despite this marked increase in P₁, the responses to an 8-AP train displayed roughly the same amount of facilitation with and without TEA. Paired-pulse ratio (PPR) values, obtained by calculating the ratio of n₂/n₁, were similar with and without TEA (two first columns in Figure 1H). When examining the time course of n_i values during the train, the maximum P_max was attained earlier in TEA (for i = 3) than in control (for i = 5), but the shapes of the two curves were not strikingly different (compare pink and red curves in Figure 1D). These observations are in line with the fact that the value of P₁ obtained in TEA (0.213) is far from the maximum attainable value of 1, leaving room for a significant increase of P linked to facilitation.

Finally, the value of P_max in the presence of 1 mM TEA was 0.472, far below the maximum of 1. Therefore, increasing release probability by a combination of synaptic facilitation and of TEA application is not sufficient to bring the release probability to 1.

A possible interpretation of the results above is that addition of 1 mM TEA was insufficient to bring the release probability of a docked SV (p) to its maximum of 1. To investigate this possibility, we tested 5 mM TEA in a new series of experiments, starting again from the same 1.5 mM [Ca²⁺]_o control condition (Figure 1B, third row). As documented below, Ca²⁺ entry was substantially increased in 5 mM TEA compared to 1 mM TEA. Nevertheless, even with 5 mM TEA, failures were often observed in response to the 1 st AP (Figure 1B, third row), suggesting that P₁, the probability that a site releases a SV after the 1^st AP, was still far from reaching its maximum of 1. In these experiments, P₁ grew from 0.053 ± 0.018 without TEA to 0.146 ± 0.036 in 5 mM TEA. This represented a mean P₁ ratio of 3.0 in individual experiments, not significantly different from the 1.9-fold ratio obtained in the previous experiments with 1 mM TEA (unpaired t-test). Globally these results suggest that TEA increases P₁ only up to a certain point.

If p is close to 1 in both 1 mM and 5 mM TEA, then in view of the relation P₁ = p δ, P₁ is close to δ in both cases. Therefore, an obvious possible explanation for the limitation of P₁ near 0.2 is that this limit represents δ, the docking site occupancy. Starting with a common submaximal δ value leaves the same room for δ-driven facilitation in 1 mM TEA experiments and in 5 mM TEA experiments. In line with this suggestion, we found similar values for the PPR (2.47 ± 0.45 in 1 mM TEA vs. 1.80 ± 0.21 in 5 mM TEA; n. s.) and for P_max (0.472 ± 0.059 in 1 mM TEA vs. 0.351 ± 0.080 in 5 mM TEA; n. s.) in 1 mM TEA and in 5 mM TEA.

In summary, values of P₁, P_max and PPR are similar in 1 mM TEA and in 5 mM TEA, suggesting saturation of p as a function of TEA concentration. In 5 mM TEA as in 1 mM TEA, the variance-mean curve remains restricted to the ascending limb of the parabola, where the release probability is less than 0.5 (grey thick lines in Figure 1C, 2nd and 3rd panels from above).

In view of our failure to reach P = 1 by addition of TEA, we next studied SV release in elevated (3 mM) [Ca²⁺]_o conditions, both with and without TEA, and we compared the results with those previously obtained under control [Ca²⁺]_o (1.5 mM).

We first examined results in 3 mM [Ca²⁺]_o without TEA. In 3 mM [Ca²⁺]_o, P₁ was increased compared to 1.5 mM [Ca²⁺]_o (0.28 ± 0.03 vs. 0.089 ± 0.015; p<0.01, unpaired t-test), and it reached a mean value similar to that obtained in 1.5 mM [Ca²⁺]_o and TEA (0.28 ± 0.03 vs. 0.21 ± 0.06). Likewise P_max was larger in 3 mM [Ca²⁺]_o compared to 1.5 mM [Ca²⁺]_o (0.46 ± 0.04 vs. 0.27 ± 0.03; p<0.01, unpaired t-test), and was similar in 3 mM [Ca²⁺]_o compared to 1.5 mM [Ca²⁺]_o and TEA (0.46 ± 0.04 vs. 0.47 ± 0.06).

The finding that P₁ and P_max take similar values in 3 mM [Ca²⁺]_o and in 1.5 mM [Ca²⁺]_o plus TEA could reflect a common limit to the strength of the synapse. Alternatively, it could be a sheer numerical coincidence. To distinguish between these two possibilities, we next tested the effects of adding 1 mM TEA starting in 3 mM [Ca²⁺]_o. P₁ significantly increased from 0.28 ± 0.03 to 0.51 ± 0.09 (p<0.05, paired t-test; exemplar experiment in lowest 2 panels of Figure 1B; group results in Figure 1E). Likewise, P_max increased from 0.46 ± 0.04 to 0.66 ± 0.08 (p<0.05, paired t-test; Figure 1F). Therefore, whereas neither P₁ nor P_max increased significantly when increasing the TEA concentration from 1 to 5 mM in 1.5 mM [Ca²⁺]_o experiments, both P₁ and P_max were significantly higher in 3 mM [Ca²⁺]_o and 1 mM TEA compared to 3 mM [Ca²⁺]_o. These results indicate that combining [Ca²⁺]_o elevation with TEA application is an effective way to overcome the limitations in P apparent at high TEA concentrations with normal [Ca²⁺]_o.

In view of the relation P = δ p, we propose as an interpretation of the above results that TEA application leads to an increase of p without a change in δ, whereas an increase in [Ca²⁺]_o results in increases in both δ and p. In this interpretation, δ remains constant if [Ca²⁺]_o is constant, so that even if TEA increases p near its maximum p = 1, P remains limited by δ. This explains the limitation of P₁ values near 0.2 in 1.5 mM [Ca²⁺]_o (Figure 1E). By contrast if [Ca²⁺]_o is increased the value of δ increases so that the maximum attainable for P₁ increases according to the change in δ. This explains the higher P₁ values observed when combining TEA and high [Ca²⁺]_o, compared to applying TEA in normal [Ca²⁺]_o (Figure 1E–F).

The proposal of differential effects of TEA and [Ca²⁺]_o on p and δ also helps explaining short-term plasticity results shown in Figure 1. As detailed earlier (Pulido and Marty, 2018), models of synaptic function based on changes in docking site occupancy predict different effects of increasing p or δ on short-term synaptic plasticity. While all synapses become depressing if δ increases to the maximum δ = 1, some synapses with low δ values remain facilitating even if p increases to p = 1. This is because if the initial value of docking site occupancy is low, even if the first AP releases all docked SVs, replenishment of initially empty docking sites can overcompensate the loss of exocytosed SVs, resulting in facilitation. In the present case, while facilitation was prominent in 1.5 mM [Ca²⁺]_o, without and with TEA, it became weak or absent in 3 mM [Ca²⁺]_o (yellow and green curves in Figure 1D). Simultaneously depression increased, particularly in the presence of TEA (green curve in Figure 1D), and the value of the PPR decreased (Figure 1H). Overall, the finding that facilitation remains when applying TEA alone is consistent with a pure effect on p, while the more complex effects observed on short-term synaptic plasticity when increasing [Ca²⁺]_o are consistent with a mixed effect on p and δ.

In several synapses including MLI synapses onto Purkinje cells, increasing [Ca²⁺]_o leads to an increase of miniature current frequency, presumably following an increase in the presynaptic resting [Ca²⁺]_i (Llano et al., 2000). At PF-MLI synapses, we found that miniature EPSC frequency rose from 0.142 ± 0.033 Hz in 1.5 mM [Ca²⁺]_o to 0.364 ± 0.040 Hz in 3 mM [Ca²⁺]_o (p<0.01; n = 6), while mean peak miniature EPSC amplitudes did not change significantly (93.6 ± 9.1 pA in 1.5 mM [Ca²⁺]_o vs. 105.3 ± 20.4 pA in 3 mM [Ca²⁺]_o; p>0.05; n = 6). These data suggest that elevating [Ca²⁺]_o results in an elevation of presynaptic [Ca²⁺]_i at PF-MLI synapses. By contrast, we found that TEA (1 mM) did not change miniature EPSC frequency (control frequency: 0.123 ± 0.014 Hz; TEA frequency: 0.116 ± 0.023 Hz, p>0.05, n = 7). These results suggest that [Ca²⁺]_o elevation, but not TEA addition, leads to an increase in the presynaptic basal [Ca²⁺]_i. As SV replenishment is enhanced by elevations of [Ca²⁺]_i (Neher and Sakaba, 2008), it is plausible that presynaptic [Ca²⁺]_i elevation increases the docking site occupancy δ.

A striking outcome of the experiments described so far is that in all cases, variance-mean parabolas were far from complete, as illustrated by the thick grey curves in Figure 1C. Specifically, in spite of our efforts, we were unable to attain a value of P larger than 0.66 in the above experiments (i. e., the value of P_max in 3 mM [Ca²⁺]_o and 1 mM TEA). This contrasts with many studies reporting P values of 0.9 or larger when performing mean-variance analysis of peak EPSC amplitudes (Branco and Staras, 2009). One possible reason for the discrepancy is that P is limited by δ, as mentioned above, and that the value of δ is particularly low in PF-MLI synapses. Another possible reason is methodological. In the traditional method of peak EPSC mean-variance analysis, receptor saturation or desensitization, as well as synaptic jitter, can lead to substantial errors (Clements, 2003; Silver, 2003). These errors worsen at high release probability, and if they are not corrected, their effects tend to overestimate P and to underestimate N. While methods exist to correct for some of these errors (Silver, 2003; Taschenberger et al., 2005), corrections are difficult to implement, often requiring dedicated experiments, and they are at best partial. To examine the potential consequences of the errors involved, we next compared in the same recordings the P value obtained from peak EPSC analysis with the value obtained from SV counts.

A comparison between the two kinds of variance-mean analyses is presented in Figure 2A in an experiment using 3 mM [Ca²⁺]_o. The variance-mean parabola indicated a P₁ value of 0.307 with the SV count method (left, red triangles) and of 0.487 with the peak EPSC method (center, open purple rectangles); P_max values were 0.436 with the SV count method and 0.710 with the peak EPSC method. As illustrated in this example, P estimates were invariably higher with the EPSC amplitude analysis than with the SV count analysis. The discrepancy became larger when correcting the variance for EPSC amplitudes for the variations in q values in the same experiment (filled purple rectangles, right panel: P₁ = 0.530, and P_max = 0.779), as is usually done in peak amplitude variance-mean analysis (Silver, 2003). To perform the correction, we assumed a CV value of 0.30 for quantal amplitude variance within one synapse, based on previous estimates (Malagon et al., 2016). Average results revealed significant discrepancies in P₁ estimates in 3 mM [Ca²⁺]_o (0.294 ± 0.032 with the SV count method vs. 0.612 ± 0.058 with the peak EPSC method, p<0.001; Figure 2B; note that the present values obtained with the peak EPSC method are in agreement with previous estimates in the same preparation using the same method: Ishiyama et al., 2014). Similarly, P_max estimates were smaller with the SV count method than with the peak EPSC method (0.456 ± 0.025 vs. 0.806 ± 0.058, p<0.01; Figure 2C). Estimated N values were larger with the SV counting method than with the peak EPSC method (3.94 ± 0.46 vs. 2.19 ± 0.29, p<0.05; Figure 2D), suggesting that the peak EPSC method leads to an underestimate of N. Finally, direct measurements of quantal current amplitudes (217 ± 22 pA) were significantly different from q estimates derived from EPSC amplitude variance-mean analysis (147 ± 24 pA, p<0.01; green vs. purple means in Figure 2E). These numbers indicate a significant (about 1.5-fold) under-estimate of q values with the EPSC amplitude variance-mean analysis compared to direct measurements. Altogether these results indicate that the peak EPSC amplitude mean-variance analysis leads to over-estimates of P₁ and P_max, and to under-estimates of N and q, likely due to a combination of receptor saturation and synaptic jitter. Consequently, the difference between the moderate maximum P values found in this work, and higher values found in earlier studies, may partially arise from uncorrected errors linked to receptor saturation and synaptic jitter in these earlier studies.

Figure 2

Download asset Open asset

The above results show that values of P per release site are markedly smaller than 1. If the value of δ were 1, as is commonly assumed, then the low range of P values would predict that P is highly sensitive to Ca²⁺ entry, because the dose-response curve of release probability as a function of [Ca²⁺]_i is cooperative at the foot of the curve (Dodge and Rahamimoff, 1967; Schneggenburger and Neher, 2000). Therefore we next investigated the changes in Ca²⁺ entry associated with [Ca²⁺]_o changes and with TEA applications. For this purpose, we performed two-photon imaging of single PF varicosities (Figure 3A) to determine Ca²⁺ entry elicited by AP trains. We started fluorescence measurements about 30 min after establishing somatic whole-cell recording to allow dye equilibration in the imaged axonal area. Ca²⁺ entry was evaluated as the relative increase of the peak fluorescence measurement of OGB-6F (ΔF/F₀: Figure 3B). Figure 3B shows an experiment where we compared ΔF/F₀ signals under control conditions and during TEA application, finding a reversible increase. In conformity with previous results (Sabatini and Regehr, 1997; Brenowitz and Regehr, 2007; Miki et al., 2016), we found that the amount of Ca²⁺ entry per AP did not change significantly during an AP train (Figure 3C). Three series of experiments were performed: adding 1 mM TEA starting from a 1.5 mM [Ca²⁺]_o solution, as in panel B (n = 5); increasing the TEA concentration from 1 mM to 5 mM, starting in a 1.5 mM [Ca²⁺]_o plus 1 mM TEA solution (n = 4); and adding 1 mM TEA in 3 mM [Ca²⁺]_o (n = 5). Group results showed significant increases of Ca²⁺ entry upon application of 1 mM TEA, both when starting from a 1.5 mM [Ca²⁺]_o solution (control: 26.7 ± 3.6%; 1 mM TEA: 40.5 ± 3.6%; p<0.01, paired t-test) and when starting from a 3 mM [Ca²⁺]_o solution (control: 43.8 ± 11.6%; 1 mM TEA: 86.4 ± 25.2%; p<0.05; paired t-test; Figure 3C, lower panel). Likewise, we found a significant increase of fluorescence measurements in 1.5 mM [Ca²⁺]_o when increasing the TEA concentration from 1 to 5 mM (1 mM TEA: 70.8 ± 10.6%; 5 mM TEA: 100.7 ± 11.3%; p<0.05; paired t-test; Figure 3C, lower panel). Therefore, even though our previous results indicate that P₁ does not increase when elevating the TEA concentration from 1 to 5 mM, the present imaging results show that Ca²⁺ entry is substantially larger at the higher TEA concentration. The opposite situation is found when comparing results in the presence of 1 mM TEA and 3 mM [Ca²⁺]_o (green bar) with those in the presence of 5 mM TEA and 1.5 mM [Ca²⁺]_o (purple bar). Whereas P₁ is larger in the first solution than in the second (Figure 1E), ΔF/F₀ signals appear similar (Figure 3C). Therefore altogether, comparison of the results in Figure 1E and Figure 3C reveals discrepancies between P₁ results and ΔF/F₀ results, indicating that p, the release probability of docked SVs, is not the only determinant of P₁.

Figure 3

Download asset Open asset

Current models of SV release assume that a presynaptic AP induces a local [Ca²⁺]_i transient at the level of vesicular Ca²⁺ sensors with a sub-ms duration and an amplitude of several 10 s of μM (Bollmann et al., 2000; Schneggenburger and Neher, 2000). The Ca²⁺ sensors respond to the [Ca²⁺]_i transient with a sequence of fast Ca²⁺ binding to a series of 4–5 cooperative binding sites, followed by an irreversible exocytosis step (Dodge and Rahamimoff, 1967; Wu and Saggau, 1994; Bollmann et al., 2000; Schneggenburger and Neher, 2000; Lou et al., 2005). We found that depending on the mode of stimulation, such models predict different dose-response curves representing the synapse output as a function of peak global [Ca²⁺]_i. In Figure 4A–B, we compare the dose-response curve for step [Ca²⁺]_i increases (Figure 4A) or for AP stimulations (Figure 4B), using [Ca²⁺]_Ii-sensitive reaction steps and [Ca²⁺]_i profiles previously developed to simulate release at PF-MLI synapses (Miki et al., 2018). Even though Ca²⁺-sensitive reaction steps are the same for the two sets of simulations both using a modified version of the allosteric model of Lou et al. (2005), as explained in Miki et al. (2018), the dose-response curve is markedly steeper in response to step [Ca²⁺]_i increases (with an apparent Hill coefficient n = 3; Figure 4A) than in response to AP stimulations (with an apparent Hill coefficient n = 1.7; Figure 4B). The difference reflects the fact that in AP-induced responses, effective calcium levels vary in a complex manner with calcium entry, due to the diffusion of entering calcium ions and of their interaction with a variety of calcium buffers. These simulations indicate that the Hill coefficient of the dose-response curve does not depend only on the number of Ca²⁺ binding steps in the kinetic model (5 in the present case), but also on the mode of stimulation (Ca²⁺ uncaging vs. AP).

Figure 4 with 1 supplement see all

Download asset Open asset

Next, we asked how the dose-response curve of Figure 4B could be used to estimate δ. For this purpose, we combined together P₁ values obtained under the various experimental conditions of Figure 1, together with peak global [Ca²⁺]_I data obtained from Figure 3, creating an experimental plot of P₁ as a function of peak global [Ca²⁺]_I (Figure 4C). In this plot, P₁ values are overall means derived from the paired control/TEA experiments of Figure 1 as well as from other experiments including the TEA/(TEA + 4-AP) experiments to be described below. In addition, P₁ values were corrected to only incorporate release events coming directly from docking sites ('one-step release'), as opposed to indirect release events coming from replacement sites ('two-step release': Miki et al., 2018). Because 2-step release events do not originate from docked SVs, they should not be considered in the estimate of δ. Fortunately, in ordinary experimental conditions, 2-step release appears only starting from the 2^nd AP in a train, and is not detectable after the 1^st AP. In these cases, no correction for P₁ is needed to account for 2-step release. Only in the most extreme condition of 3 mM [Ca²⁺]_o + TEA was a correction necessary (by 18%). Also, peak global values from Figure 3C were corrected for partial dye saturation. To relate the corrected experimental results in Figure 4C to the simulation in Figure 4B, we note that in Figure 4B, the ordinate is p, the release probability of docked SVs, whereas in Figure 4C, the ordinate is P₁, the global release probability in response to the first stimulus in a train. As already pointed out, P₁ data points incorporate both δ and p factors, according to the relation P₁ = δ p. Clearly, P₁ data differ from the curve in Figure 4B and furthermore, they cannot be fit from this curve by using a single scaling factor. However, acceptable fits are obtained when using two different scalings depending on [Ca²⁺]_o (Figure 4C). In these plots, δ appears as the asymptotic value of P₁ at high global [Ca²⁺]_i, yielding δ = 0.22 at 1.5 mM [Ca²⁺]_o and δ = 0.47 at 3 mM [Ca²⁺]_o (Figure 4C). Altogether this analysis suggests that δ values increase as a function of [Ca²⁺]_o and therefore, that changes in both δ and p contribute to the sensitivity of P as a function of [Ca²⁺]_o.

It may be asked whether the above results may be extended to other blockers of potassium channels than TEA. Both in MLI terminals and in the calyx of Held, 4-AP increases synaptic transmission more efficiently than TEA (Tan and Llano, 1999; Ishikawa et al., 2003). In MLI terminals, 4-aminopyridine (4-AP) markedly increases calcium entry whereas TEA alone is almost ineffective; combining TEA and 4-AP together leads to an even larger AP-induced calcium elevation than when applying 4-AP alone (Tan and Llano, 1999). In view of these earlier results, we expected that combining TEA and 4-AP in our preparation should lead to a very large presynaptic calcium entry. Nevertheless, supposing that the model of Figure 4 is valid, P₁ should not be able to increase beyond the value of δ. Group results showed P₁ values of 0.140 ± 0.053 in 1.5 mM [Ca²⁺]_o and 1 mM TEA without 4-AP, and of 0.187 ± 0.048 after further addition of 15 μM 4-AP (Figure 4—figure supplement 1). After adding 4-AP, both 2-step release and facilitation were very pronounced, probably due to a marked increase in calcium entry. After correction for 2-step release, P₁ in TEA + 4-AP was 0.125 ± 0.032, a value similar to that in TEA alone. Meanwhile, P_max increased from 0.263 ± 0.046 in TEA alone to 0.550 ± 0.149 in TEA + 4-AP (Figure 4—figure supplement 1), a value similar to that obtained in 3 mM [Ca²⁺]_o and 1 mM TEA (Figure 1). Thus, while P₁ was limited by the resting docking site occupancy, the synapse could deliver more SVs upon repetitive stimulation as SV replenishment was increased. In conclusion, whether calcium entry is increased with TEA alone or with a combination of TEA and 4-AP, P₁ remains limited by the value of δ.

In the present work we take advantage of the recently developed method of SV counting (Malagon et al., 2016) to examine variations of the release probability of single docking sites (P) during short term synaptic plasticity, following addition of TEA, and when changing [Ca²⁺]_o. Our main finding is that the synaptic strength at the beginning of an AP train (P₁) has a maximum attainable value that depends on the value of [Ca²⁺]_o. With 1.5 mM [Ca²⁺]_o, the maximum of P₁ is 0.22 SV per AP and per release/docking site, whereas with 3 mM [Ca²⁺]_o, the maximum of P₁ is 0.47.

Our interpretation of this finding is that when the release probability of docked SVs, p, is close to 1, P₁ is limited by the resting docking site probability δ, and that δ is lower in 1.5 mM [Ca²⁺]_o compared to 3 mM [Ca²⁺]_o. The maximum output per site gives an estimate of resting docking site occupancy, a basic parameter of synaptic function that has been difficult to evaluate.

Our results confirm earlier suggestions that δ is less than 1 at PF-MLI synapses (Miki et al., 2016; Miki et al., 2018). In addition, we now provide reliable δ estimates. The present finding that δ = 0.22 under physiological conditions (1.5 mM [Ca²⁺]_o) implies that there is a large discrepancy between the release probability per docking site, P₁, and the release probability per occupied docking site, p. By inverting the P₁ = δ p relation, we obtain from P₁ = 0.089 and δ = 0.22 that p = 0.39 in 1.5 mM [Ca²⁺]_o. At 3 mM [Ca²⁺]_o, we have P₁ = 0.28 and δ = 0.47, so that p = 0.60. Therefore both at 1.5 mM [Ca²⁺]_o and at 3 mM [Ca²⁺]_o, the actual release probability of docked SVs is substantially higher than it would appear without the correction. In the presence of TEA, p values calculated from the pooled data of Figure 4 approach 1 (1.5 mM [Ca²⁺]_o: p = 0.85 in 1 mM TEA and p = 0.95 in 5 mM TEA; 3 mM [Ca²⁺]_o: p = 0.89 in 1 mM TEA, value corrected for 2-step release). Therefore in the present preparation, P₁ values in the presence of TEA give reasonably good approximations of δ.

This correction is important when relating release probability and local calcium profiles. For example, at synapses between PFs and Purkinje cells, P values are comparatively small, yet EGTA application experiments suggest tight coupling between SVs and Ca²⁺ entry, indicating high p values (Schmidt et al., 2013). A low δ value at these synapses, similar to the value found here for PF-MLI synapses, is a plausible explanation for these apparently discrepant results.

Increasing evidence suggests that part of inter-synaptic variability may reside in differences in docking site occupancy. Specifically, synapses follow a gradient from low release probability, facilitating synapses to high release probability, depressing synapses, presumably depending on the proportion of occupancy of release sites by SVs (Pan and Zucker, 2009; Pulido and Marty, 2018; Neher and Brose, 2018). PF-MLI synapses have a low δ value in normal physiological conditions (0.22 at 1.5 mM [Ca²⁺]_o according to the present study) and are facilitating, indicating that they belong to the first category, called 'tonic synapses' (Pan and Zucker, 2009; Neher and Brose, 2018). By contrast, SV release statistics at simple MLI-MLI synapses indicate a higher resting δ value near 0.5–0.7 in 2 mM [Ca²⁺]_o (Trigo et al., 2012; Pulido et al., 2015). Furthermore MLI-MLI synapses are depressing, and altogether MLI-MLI synapses appear closer to the 'phasic synapse' category (Pan and Zucker, 2009; Neher and Brose, 2018).

Because of the difficulty of separating δ from p, estimates of δ in most preparations are not available. Nevertheless, earlier variance/mean analysis studies, often carried out on peak synaptic current amplitudes, give constraints on δ. Maximum P values reported in these studies are usually in the range 0.5–1 (Branco and Staras, 2009). In view of the relation P = δ p, since p<1, such large P values suggest resting δ values above 0.5. However these maximum P values are almost always obtained in high [Ca²⁺]_o conditions (calyx of Held: Taschenberger et al., 2005; PF-MLI synapses: Ishiyama et al., 2014; fly neuromuscular junction: Reddy-Alla et al., 2017; hippocampal cultures: Sakamoto et al., 2018). In view of the present evidence that δ increases as a function of [Ca²⁺]_o, the high P values at high [Ca²⁺]_o are compatible with δ values < 0.5 under normal [Ca²⁺]_o conditions, at least in some (tonic) synapses. Apart from using high [Ca²⁺]_o conditions, another factor that may have contributed in some studies to high apparent P values is uncorrected synaptic jitter and receptor saturation, as documented above (Figure 2). One notable exception is the climbing fibre-Purkinje cell synapse, that may represent an extreme case of phasic synapses, and that exhibits P values above 0.9 at normal [Ca²⁺]_o (Silver et al., 1998). Overall our conclusions appear well compatible with earlier P estimates obtained with variance/mean analysis of peak synaptic current amplitudes.

An increase in [Ca²⁺]_o likely acts on the docking site occupancy by inducing an increase in [Ca²⁺]_i after reequilibration of the cytosolic Ca²⁺ concentration. It is widely accepted that following docking site depletion, docking site replenishment is sensitive to the bulk [Ca²⁺]_i (Neher and Sakaba, 2008). This suggests that the proportion of occupied docking sites at rest may be enhanced by a [Ca²⁺]_i increase, as the equilibrium is shifted in favor of SV binding to the docking site. Recent electron microscopy studies are in accord with this proposal (Chang et al., 2018; Kusick et al., 2018). Several AZ molecules could be responsible for Ca²⁺-dependent SV docking, including synaptotagmin-1 (Chang et al., 2018) and Munc13 (Shin et al., 2010). Another possibility is that Ca²⁺-dependent SV docking is due to activation of actomyosin-driven SV movement by [Ca²⁺]_i (Lee et al., 2012; Miki et al., 2016).

A [Ca²⁺]_i –induced change in resting docking site occupancy may not only apply when changing [Ca²⁺]_o, as proposed in the present work, but it may also occur following a number of other manipulations leading to a change in the global presynaptic [Ca²⁺]_i concentration. In the calyx of Held, as well as in MLI terminals, subthreshold presynaptic depolarization activates low-threshold voltage-gated Ca²⁺ channels, thus increasing release probability (Awatramani et al., 2005; Bouhours et al., 2011). Activation of Cl^--permeant presynaptic glycinergic or GABAergic receptors in these synapses engages this pathway, leading to a presynaptic [Ca²⁺]_i increase and to an enhancement of synaptic strength (Turecek and Trussell, 2001; Trigo et al., 2007; Zorrilla de San Martin et al., 2017). Activation of cationic channels by presynaptic nicotinic receptors, AMPA receptors, or NMDA receptors, is also able to increase presynaptic [Ca²⁺]_i, in various preparations including MLI terminals, and these actions have been shown to increase release (Sharma and Vijayaraghavan, 2003; Rossi et al., 2008; Rossi et al., 2012). Finally, enhancement of neurotransmitter release by passive spread of subthreshold somatic depolarization in hippocampal and cortical synapses ('analog signaling') depends on an increase in presynaptic [Ca²⁺]_i concentration (Shu et al., 2006; Alle and Geiger, 2006). In view of the present results, it is possible that these effects are mediated at least in part by an increase in resting docking site occupancy. Altogether, [Ca²⁺]_i-driven changes of docking site occupancy appears as a simple and powerful mechanism of gain control of synaptic strength that may be used by various forms of neuromodulation.

Sagittal slices (200 μm thick) were prepared from the cerebellar vermis of Sprague-Dawley rats (P12‒P17) following the animal care guidelines of Paris Descartes University (approval no. A-750607). Recordings were obtained from either basket or stellate cells in the molecular layer; these cells were collectively called molecular layer interneurons (MLIs). The composition of the extracellular solution was (in mM): 130 NaCl, 2.5 KCl, 26 NaHCO₃, 1.3 NaH₂PO₄, 10 glucose, 2 CaCl₂, and 1 MgCl₂ (osmolarity: 300 mosm). This solution was equilibrated with 95% O₂ and 5% CO₂ (pH 7.4). The internal recording solution contained (in mM): 144 K-gluconate, 6 KCl, 4.6 MgCl₂, 1 EGTA, 0.1 CaCl₂, 10 HEPES, 4 ATP-Na, 0.4 GTA-Na; pH 7.3 (osmolarity: 300 mosm). Recordings were at 30‒34°C.

To study simple PF-MLI synapses, potential postsynaptic MLIs were recorded under voltage clamp at ‒60 mV. For extracellular stimulation, a monopolar pipette was filled with internal solution. NMDA receptors and GABA_A receptors were blocked by inclusion of D(‒)−2-amino-5-phosphonopentanoic acid (APV, 50 μM) and gabazine (15 μM). Two alternative procedures were used to find an appropriate location for electrical stimulation. In the first procedure, an electrical stimulation pipette was located in the molecular layer, above the dendritic field of the recorded MLI, and the pipette position and stimulation strength were adjusted to obtain minimal stimulation (Miki et al., 2016). In the second procedure, we puff-applied the internal solution including 150 mM K⁺ from a pipette using small pressure steps while moving the pipette in the granule cell layer (Miki et al., 2017). When we found a burst-like EPSC response in the postsynaptic cell, we reduced the pressure of puffing to better define the spot for stimulation. Then we switched to electrical stimulation, using the same pipette, and we adjusted stimulation intensity to fire a connected granule cell under minimal stimulation conditions. With either procedure, definitive acceptance of the experiment as a usable simple synapse recording occurred after analysis and depended on three criteria (Malagon et al., 2016): (i) a decrement of EPSC amplitudes of second events in a pair, reflecting activation of a common set of receptors belonging to one PSD; (ii) a Gaussian distribution of EPSC amplitudes with a coefficient of variation (CV) less than 0.5; and (iii) stability of the overall responsiveness over time. Single stimulations and trains of two or eight stimulation pulses were applied repetitively with intervals of 10 s between sweeps. Statistical data were derived from sequences of 10‒30 trains.

We determined occurrence times of individual EPSCs based on deconvolution analysis, as detailed in Malagon et al. (2016), and we built latency distributions by averaging the occurrence times across experiments. We briefly describe the analysis here. Firstly, we made an average of single EPSCs obtained during asynchronous release to obtain a template in a given synapse. Then the average mEPSC was fit by triple-exponential function with five free parameters (rise time, amplitude, fast decay time constant, slow decay time constant, and amplitude fraction of slow decay). Next, mEPSC and individual data traces were deconvolved using the five parameters. The deconvolved mEPSC resulted in a narrow spike, and the deconvolved data traces resulted in sequences of spikes. Finally, we fit a given deconvolved trace by a sum of scaled narrow spikes in order to obtain the timing of each event. The amplitude parameter was free because the peak EPSC amplitude varied during a train due to receptor saturation and desensitization. The above procedure had a detection limit that caused a failure of separation of two events occurring within 0.2 ms. To correct for missed events, we split into two the events having amplitudes at least 1.7 times larger than the average amplitude obtained during asynchronous release.

In TEA experiments, we checked somatic potassium currents by applying 3 ms voltage step to 0 mV until TEA reduced potassium current amplitudes to a stable level, and then we started to collect data.

For Ca²⁺ imaging experiments, sagittal (200 μm) slices were prepared using a modified extracellular saline, as detailed in Brenowitz and Regehr (2007). Experiments were conducted in the same conditions as in the electrophysiology experiments (32‒34°C, with the 3 mM extracellular Ca²⁺ saline including APV and gabazine). Granule cells were loaded under whole-cell recording with a solution containing (in mM): 140 K-gluconate, 5.4 KCl, 4.1 MgCl₂, 9.9 HEPES, 0.36 Na-GTP, 3.6 Na-ATP, 500 μM of Oregon green 488 BAPTA-6F (OGB-6F: K_d for calcium of 5.1 μM; Invitrogen) and 20 μM Alexa-594 (Invitrogen). Imaging was performed with a custom-built 2-photon system, with 820 nm excitation provided by a MaiTai-Sapphire laser (Spectra Physics, USA). In order to visualize the granule cell axon, large raster scans were performed while acquiring the Alexa 594 fluorescence with a red channel photomultiplier (Hamamatsu H7422 PA-sel, bandpass emission filter 635 ± 65 nm, Chroma Technology; or avalanche photodiode Perkin Elmer, SPCM-AQR-13). Single varicosity imaging was performed using raster scans of 5 by 2 μm dimensions at dwell times of 2 ms. The granule cells were kept under current clamp conditions. We found that keeping a hyperpolarized holding potential improved recording stability, so that resting membrane potential was kept around ‒90 mV. APs were evoked by 1 ms steps of 350‒500 pA. Stimulation protocols were 4 or 8 APs at 100 Hz and were repeated every 1 min. Calcium signaling was analyzed in the pixels encompassing the varicosity in terms of fluorescence changes relative to pre-stimulus values (ΔF/F_o, expressed in %) with software written in the Igor Pro programming environment (Wavemetric, Lake Oswego, OR, USA).

To fit the data in Figure 3C, we assumed gradual saturation of OGB-6F as a function of Ca²⁺ entry, with equal contributions for each AP (Miki et al., 2016). Estimated K_d values are indicated in Figure 3C, upper panel. As before, dye saturation was set at 6.5 times the baseline level (Miki et al., 2016). The data in Figure 4C are based on the pooled data in the lower panel of Figure 3C, after correcting for dye saturation.

We calculated the release probabilities for various Ca²⁺ waveforms by numerical integration of differential equations using Igor Pro with a time interval of 0.001 ms as described (Miki et al., 2018). We performed the simulations using a simple one-step release model without replenishment at docking sites (Miki et al., 2018), Figure 1) with the allosteric model of Lou et al. (2005). According to the fit of data in Miki et al. (2018), we used the following parameter values for the allosteric model: K_on = 5×10⁸ M⁻¹s⁻¹, K_off = 5000 s⁻¹, b = 0.75, γ = 1500 s⁻¹, and f = 31.1. γ is the fusion rate of the 5Ca-binding state (V_5Ca), which is identical to l₊ x f⁵. f is a factor determining the increase in vesicle fusion rate upon Ca²⁺ binding. The fusion rate at V_0Ca-4Ca is l₊ x f^0-4. b is a cooperativity factor (Lou et al., 2005). In Figure 4A-B, the probability of occupancy at docking site was set to 1. For the Ca²⁺ uncaging simulation, we assumed step [Ca²⁺] increases from a basal [Ca²⁺] of 50 nM to 0.5, 1, 2, 4, 8, 16, or 24 μM, and we calculated the release probabilities by dividing the total number of vesicle release between 0 and 5 ms by the docking site number. For the AP simulation, we used the local [Ca²⁺] obtained from Ca²⁺ simulation at 40 nm distance from the nearest Ca²⁺ channels (Miki et al., 2018). Global [Ca²⁺] was obtained at the center of a bouton in Ca²⁺ simulation. The Ca²⁺ simulation was described in detail in a previous report (Miki et al., 2018). Briefly, we calculated the spatiotemporal distribution of [Ca²⁺] at a PF bouton. We used the bouton size of 0.9 × 0.5×0.5 μm³ cuboid, and the Ca²⁺ channel distributions based on the electron microscopy observations (Miki et al., 2017; Miki et al., 2018). The amplitude of Ca²⁺ influx through each channel during an AP was set at 0.2 pA and a fraction of channel open probability during an AP of ~0.7. We assumed a Gaussian-shaped Ca²⁺ influx with a half-amplitude duration of 0.34 ms. Other model parameters were set following Miki et al. (2018). To create Ca²⁺ waveforms having various amplitudes for the AP simulation, we multiplied the simulated local and global [Ca²⁺] by a factor of 0.25, 0.5, 1, 1.5, 2 and 3 while keeping the resting [Ca²⁺] of 50 nM. We calculated the release probabilities by dividing the total number of vesicle release between 0 and 5 ms after the beginning of the local [Ca²⁺] increase by the number of docking sites.

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

We thank Isabel Llano for her advice and support in two-photon imaging experiments. This work was supported by CNRS (UMR 8118, and UMR 8003), by the European Community (ERC Advanced Grant ‘Single Site’ to AM, nb. 294509), by JSPS (KAKENHI Grant JP18K06472 to TM, and Core-to-Core Program A, Advanced Research Networks), and by Fondation pour la Recherche Médicale (grant SPF201809007190).

1. Received: September 23, 2019
2. Accepted: March 23, 2020
3. Accepted Manuscript published: March 31, 2020 (version 1)
4. Version of Record published: April 6, 2020 (version 2)

© 2020, Malagon et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.