Abstract
Hair cell ribbon synapses exhibit several distinguishing features. Structurally, a dense body, or ribbon, is anchored to the presynaptic membrane and tethers synaptic vesicles; functionally, neurotransmitter release is dominated by large EPSC events produced by seemingly synchronous multivesicular release. However, the specific role of the synaptic ribbon in promoting this form of release remains elusive. Using complete ultrastructural reconstructions and capacitance measurements of bullfrog amphibian papilla hair cells dialyzed with high concentrations of a slow Ca2+ buffer (10 mm EGTA), we found that the number of synaptic vesicles at the base of the ribbon correlated closely to those vesicles that released most rapidly and efficiently, while the rest of the ribbon-tethered vesicles correlated to a second, slower pool of vesicles. Combined with the persistence of multivesicular release in extreme Ca2+ buffering conditions (10 mm BAPTA), our data argue against the Ca2+-dependent compound fusion of ribbon-tethered vesicles at hair cell synapses. Moreover, during hair cell depolarization, our results suggest that elevated Ca2+ levels enhance vesicle pool replenishment rates. Finally, using Ca2+ diffusion simulations, we propose that the ribbon and its vesicles define a small cytoplasmic volume where Ca2+ buffer is saturated, despite 10 mm BAPTA conditions. This local buffer saturation permits fast and large Ca2+ rises near release sites beneath the synaptic ribbon that can trigger multiquantal EPSCs. We conclude that, by restricting the available presynaptic volume, the ribbon may be creating conditions for the synchronous release of a small cohort of docked vesicles.
Introduction
At hair cell afferent synapses, a synaptic ribbon is affiliated with the sites of Ca2+ channel clustering and exocytosis (Zenisek et al., 2004; Frank et al., 2010). The ribbon tethers a shell of synaptic vesicles, a small subset of which dock onto, or reside very close to, the presynaptic membrane. Remarkably, hair cell ribbon synapses can sustain high release rates of several hundred vesicles per second over prolonged periods (Parsons et al., 1994). This seemingly inexhaustible supply of vesicles is presumed to result from a vigorous replenishment of the ribbon-associated vesicles by the multitude of vesicles that populate the surrounding cytoplasm. Mechanisms for signaling, transporting, and mobilizing these vesicles remain unclear, although vesicle pool replenishment is influenced by Ca2+ and depends on a Ca2+-sensor protein unique to hair cells called otoferlin (Roux et al., 2006; Johnson et al., 2008; Pangrsic et al., 2010).
Another hallmark of ribbon function is multivesicular release, in which more than one vesicle fuses at a single synapse to produce EPSCs that are much larger than mEPSCs (Glowatzki and Fuchs, 2002; Grant et al., 2010). In frog hair cells, a transition from mEPSCs to predominantly large EPSCs occurs with a small difference in presynaptic membrane potential of only a few millivolts (Li et al., 2009), suggesting that this transition occurs with the incremental opening of a few Ca2+ channels per synapse (Brandt et al., 2005; Jarsky et al., 2010). Two leading hypotheses have been proposed to explain multivesicular release. One involves the homotypic fusion of vesicles (so called, compound fusion) (Matthews and Sterling, 2008); the other, the coordinated fusions of several docked vesicles at the plasma membrane (Singer et al., 2004). However, both explanations are confounded by the fact that multivesicular release remains robust even upon the introduction of strong presynaptic Ca2+ buffering (Goutman and Glowatzki, 2007), which is estimated to restrict the mean distance of free Ca2+ diffusion to <10 nm from the pore of the Ca2+ channel for 10 mm BAPTA (Bucurenciu et al., 2008; Parekh, 2008). In the context of these tight spatial constraints, the large size of some multivesicular EPSCs (composed of five or more vesicles) is baffling, especially considering the low Ca2+ sensitivity of the release machinery of the hair cell (Beutner et al., 2001).
The complex relationship between vesicle populations and Ca2+ influx into the crowded presynaptic space under a synaptic ribbon has yet to be fully characterized. Here, we combine electrophysiology with a thorough ultrastructural analysis to determine the size of morphologically and functionally defined vesicle pools at an adult hair cell synapse. We find that a slow Ca2+ buffer (10 mm EGTA) can decrease vesicle recruitment rates and thus uncovers the size of two distinct vesicle pools. Furthermore, by incorporating our ultrastructural findings into Monte Carlo simulations of Ca2+ influx at ribbon-type active zones, we suggest that volume exclusion by the presence of the ribbon and its associated vesicles can produce a sharp local Ca2+ rise and account for the observed characteristics of multivesicular release at auditory hair cell synapses.
Materials and Methods
Tissue and hair cell preparation.
All procedures followed the Oregon Health & Sciences University or National Institutes of Health (protocol 1215-08) approved animal care protocols and guidelines. For physiology studies, adult bullfrogs (Rana catesbeiana) of either sex were sedated by placing them in an ice bath for ∼20 min, and then double-pithed and decapitated. Amphibian papillae were carefully dissected from the frog's ear. To expose hair cells and afferent fibers for patch clamping, amphibian papilla were stretched out and split with fine microtools in a recording chamber with oxygenated artificial perilymph (Keen and Hudspeth, 2006; Li et al., 2009). Artificial perilymph contained the following (in mm): 95 NaCl, 2 KCl, 2 CaCl2, 1 MgCl2, 25 NaHCO3, 3 glucose, 1 creatine, 1 sodium pyruvate, pH adjusted to 7.3, and bubbled with 95% O2 and 5% CO2. Salts and chemicals were bought from Sigma-Aldrich or from Ascent. During the recordings, the preparation was perfused continuously (2–3 ml/min) with oxygenated artificial perilymph.
Electrophysiology.
Preparations were viewed with an Olympus BX51WI or Zeiss Axioscope II FS microscope equipped with a 40 or 60× water-immersion objective and a CCD camera (XC-75; Sony) with contrast enhancement (C2741; Hamamatsu). Whole-cell recordings were performed with a double EPC-9/2 or EPC-10/2 (HEKA Elektronik) patch-clamp amplifier and Pulse software (HEKA) at room temperature. Patch pipettes of borosilicate glass (World Precision Instruments) were pulled to resistances of 5–6 MΩ for hair cells and 8–10 MΩ for afferent fibers with a Narishige model PP-830 puller or a Sutter model P-97 puller. The standard internal solution contained the following (in mm): 77 Cs-gluconate, 20 CsCl, 1 MgCl2, 10 TEA (tetraethylammonium)-Cl, 10 HEPES, 2 EGTA, 3 Mg-ATP, 1 Na-GTP, and 5 Na2-phosphocreatine, adjusted to pH 7.3 with CsOH. The averaged uncompensated series resistances (Rs) in whole-cell recordings were 13.2 ± 0.4 MΩ for 27 hair cell recordings and 31 ± 9 MΩ for 15 afferent-fiber recordings. We did not use electronic series resistance compensation for most of the fiber recordings because this added noise and often resulted in ringing and the loss of the recording. Errors in the voltage command due to uncompensated series resistance were not corrected off-line. In five afferent fibers held at −90 mV, we recorded spontaneous EPSCs before and after 50–70% series resistance compensation. No significant changes were observed for peak amplitudes (124 ± 34 vs 159 ± 54 pA), charge transferred (152 ± 58 vs 170 ± 69 fC), or for 10–90% rise time (0.27 ± 0.01 vs 0.26 ± 0.02 ms; n = 5). Hair cells and afferent fibers were held at −90 mV after being corrected for liquid junction potentials. The current signal was low-pass filtered at 2 kHz and sampled at 20 kHz or higher.
Capacitance measurements.
Whole-cell membrane capacitance (Cm) was measured from hair cells under voltage clamp with the “Sine + DC” method (Lindau and Neher, 1988) using a double EPC9/2 amplifier and Pulse software. Patch pipettes were coated with dental wax (Cavex) to minimize noise and pipette capacitance. A 1 kHz sine wave (50 mV peak-to-peak) was superposed on the holding potential of the hair cells of −90 mV, and the resulting current was used to calculate Cm via the Pulse software emulator of a lock-in amplifier (Gillis, 2000). The change of capacitance, ΔCm = Cm (response) − Cm (baseline), was used as a proxy of total synaptic vesicle exocytosis from a single hair cell. An average of Cm (response) and Cm (baseline) was obtained by averaging Cm data points before and after the depolarizing pulse. Data analysis was performed with IGOR Pro software (Wavemetrics) and Prism (GraphPad Software).
Resonant frequencies of hair cells.
To measure the resonant frequency, we injected step current into the hair cells and measured the resulting membrane potential oscillations. These membrane potential oscillations decayed exponentially in amplitude and reached a steady-state membrane potential (Vss). To calculate the resonant frequency (fc), we fit the current-clamp data points using the following equation: V(t) = A · sin(2π · fc · (t − t0) + ϕ) · exp[−(t − t0)/τ] + Vss, where A is the amplitude of voltage oscillation, ϕ is the phase, and τ is the single exponential decay time constant (Crawford and Fettiplace, 1981).
Noise analysis on Ca2+ tail currents.
Nonstationary noise analysis of Ca2+ currents was used to estimate the number of Ca2+ channels per hair cell and their single-channel current (Roberts et al., 1990). To obtain a complete relationship between the mean and variance of the Ca2+ currents for low- and high-amplitude Ca2+currents, we included 5 μm 1,4-dihydro-2,6-dimethyl-5-nitro-4-[2-(trifluoromethyl)phenyl]-3-pyridinecarboxylic acid, methyl ester (BayK 8644), in the bath solution to maximize the open probability of the L-type Ca2+ channels of the hair cell (Fox et al., 1987). To keep resting Ca2+ levels very low, hair cells were first held at −90 mV, and then stepped to −110 mV to relieve Ca2+ channels from any Ca2+ or voltage-dependent inactivation. This was followed by a step depolarization to +30 mV, opening all (or nearly all) of the Ca2+ channels. A membrane voltage step back to −90 mV produced a large and sudden increase of driving force. As a result, a large Ca2+ tail current was elicited through the open channels, which decreased quickly as the Ca2+ channels reclosed (see Fig. 6B). For each hair cell, we collected a total of 100 leak-subtracted current traces (elicited every 100 ms), and the mean and variance of these currents were calculated point-by-point using a custom-made IGOR Pro analysis program. After plotting the variance against the mean current, the data were fit to a parabolic function: Var(I) = i · I − I2/NCa + Enoise, where I is the mean current, i is the single-channel current, Enoise is the electrical noise, and NCa is the number of Ca2+ channels.
Serial section reconstructions.
For electron microscopy, bullfrogs were anesthetized by intraperitoneal injection of a buffered solution of MS-222 (ethyl 3-aminobenzoate methanesulfonic acid) followed by double-pithing and decapitation. Under a dissecting scope, the entire bullfrog inner ear was excised from the surrounding bone and fixed immediately by immersion in 2.5% glutaraldehyde/2% paraformaldehyde/3 mm CaCl2/0.1 m cacodylate buffer, pH 7.3. Alternatively, some bullfrog tissue (“zero-calcium” condition) was incubated for 10 min in a Ca2+-free frog Ringer's (in mm: 110 Na, 2 K, 3 Mg, 1 EGTA, 5 HEPES, 118 Cl, 3 d-glucose, 1 pyruvate, 1 creatine; pH 7.3) before fixation in 2.5% glutaraldehyde/2% paraformaldehyde/0.1 m cacodylate buffer, pH 7.3. After fixing, the tissue was postfixed in 1% osmium tetroxide in 0.1 m cacodylate buffer, pH 7.3, for 1 h. For some preparations, the tissue was then stained en bloc with 2% aqueous uranyl acetate for 1 h. The tissue was washed with dH2O, dehydrated in a graded series of acetone or PPO [poly(p-phenylene oxide)], and embedded in EMbed 812 resin (Electron Microscopy Sciences).
For whole-cell reconstructions, sets of 150-nm-thick serial sections (typically 150–200 sections per set) were collected on 2 × 1 mm single-slot grids coated with Formvar. The sets were then imaged in a JEOL 1010 electron microscope operating at 80 kV at a magnification of 2500×. Overlapping images were taken of each section and later montaged in either Reconstruct (Fiala, 2005) or Adobe Photoshop CS3 (Adobe Systems). From this data set, we quantified the number, diameter, and spacing of all the synaptic ribbons and mapped their afferent fiber connectivity. The number of ribbons per hair cell was either counted directly from four hair cells that were entirely reconstructed, or extrapolated for the hair cells that were reconstructed halfway through their longitudinal diameter (using ∼53 consecutive sections, or ∼8 μm).
Determination of vesicle pool sizes from serial sections.
For single-ribbon synapse reconstructions, serial 50-nm-thick sections were collected on grids and imaged at 40,000× as described above. Alignment, reconstructions, and morphometric analysis were made in Reconstruct software. We used second-order polynomial fits to describe the relationship between ribbon diameters and the number of associated synaptic vesicles. For each ribbon, we counted the number of synaptic vesicles that were within 30 nm from the ribbon surface—this represented the ribbon-attached vesicle pool (RAP). A subset of the ribbon-attached vesicles, where the distance of the vesicle was within <20 nm from the presynaptic membrane, was considered the docked pool or the immediately releasable pool (IRP)—the presumed pool of vesicles that release first during a depolarization (Lenzi et al., 1999; Schnee et al., 2011). The value of 20 nm was chosen because t-SNARE and v-SNARE proteins have a cytosolic length of ∼10 nm, so SNARE interactions may occur from a distance of 20 nm from the presynaptic plasma membrane (Zenisek et al., 2000; Castorph et al., 2010). By combining our whole-cell reconstruction data with our individual ribbon reconstruction data and taking into account the ribbon size distribution, we estimated the total size of the RAP and IRP per hair cell.
Synaptic vesicle size distributions.
The diameters of vesicles in the RAP were calculated in our serial sections by drawing a best-fit ellipse over each synaptic vesicle membrane through the center of the lipid bilayer. The diameters of each vesicle profile, approximated as a circle, were then calculated from the area of each ellipse.
Tomographic reconstructions.
Plastic sections (200 nm) taken from the preparations described above were collected on hex mesh grids. Colloidal gold (10 nm; Sigma-Aldrich) was applied to both sides of the sections to act as fiducial markers before imaging with a Zeiss 922 transmission electron microscope operating at 160 kV and equipped with a 2k Gatan CCD camera and Digital Micrograph (Gatan) tomography acquisition software. Tilt series (from −65 to +65°) images were acquired at magnifications producing image pixel sizes ranging from 0.7 to 1.4 nm. Images were aligned with IMOD software (Kremer et al., 1996), reconstructed using EM3D software (Ress et al., 1999), and then rendered, segmented, and analyzed using Reconstruct software.
Determination of vesicle pool sizes from tomographic reconstructions.
To estimate vesicle pools in our tomographic reconstructions of ribbon synapses, we used only data sets that included the center of the ribbon. Because ribbons were approximately spherical, for each data set we measured the largest cross-sectional area of the ribbon to estimate a corresponding ribbon diameter value. Ribbon-attached vesicles formed a shell enclosing the ribbon, so vesicle pool estimates were made based on the percentage of the vesicle shell reconstructed. The diameter of this shell incorporated the diameter of the ribbon, one-half the diameter of attached vesicles, and the length (mean length ± SD, 30.2 ± 6.0; n = 10) of the filamentous tether attaching the vesicle to the ribbon measured from the tomographic reconstructions. For these reconstructions, we defined the IRP in the same way as in the serial section reconstructions. We also counted vesicles associated with the ribbon and touching the presynaptic membrane, a feature strictly discernible only by using electron tomography (Lenzi et al., 2002).
Monte Carlo modeling.
All simulations were done using MCell 3 software (Stiles et al., 1996; Stiles and Bartol, 2001) and included 3D models generated in the open-source program Blender (www.blender.org) (Czech et al., 2009). Results from these simulations were analyzed using Excel (Microsoft) and IGOR Pro software. Model parameters were sourced from our ribbon synapse reconstructions. Specifically, a model of a 200 nm, spherical ribbon was placed opposing a presynaptic membrane that contoured the base of the ribbon. A shell of 78 total vesicles (35 nm outer diameter) was added 30 nm from the edge of the ribbon to account for the length of the tether. The ribbon and vesicles were considered to be reflective to Ca2+, unbound buffer, and bound buffer. At the base of the ribbon, 14 (see Fig. 7) or 15 (see Fig. 8) of the vesicles were clustered in the 37 nm gap between the ribbon and presynaptic membrane. A total of 45 Ca2+ channels were randomly tiled underneath the ribbon on two stripes of presynaptic membrane (33 × 100 nm) spaced 50 nm apart, consistent with previous accounts of channel arrangement at ribbon synapses (Roberts et al., 1990; Frank et al., 2010). Each open channel was given a Ca2+ current consistent with our noise analysis estimates (see Fig. 6). Because local free Ca2+ concentrations stabilized within 1–2 ms after channels opened, stochastic timing of channel openings and closings was simplified: channels opened together and remained open for the duration of the simulation. A box (400 × 400 × 300 nm) enclosed the entire model. The walls of the box released unbound buffer in proportion to Ca2+ influx from the channels, effectively “clamping” its concentration in the simulation. The box walls absorbed both free Ca2+ and bound buffer. The diffusion coefficients for free Ca2+ (DCa = 200 μm2 s−1) and for Ca2+ buffer (bound or unbound, DBu = 20 μm2 s−1) are the same as those used in a previous study (Roberts, 1994). The time step used in the simulations (30 ns) resulted in average diffusion lengths per time step of 5.5 nm for Ca2+ and 1.7 nm for buffer. Our buffer Ca2+ binding rates (kON = 9.6 μm−1 s−1 and kOFF = 0.86 s−1; or kON = 500 μm−1 s−1 and kOFF = 96 s−1) were chosen to reflect EGTA and BAPTA, respectively (Bortolozzi et al., 2008). Unbound buffer was allowed to equilibrate to a uniform concentration before Ca2+ channels were opened. At eight locations, progressively distant from the base of the ribbon, 10 nm boxes were incorporated into the simulation to monitor (in 50 μs increments) the concentrations of buffer and free Ca2+ following the opening of channels (see Fig. 7A).
Results
Multiquantal EPSCs evoked under strong Ca2+ buffering are highly synchronous
For this comparative physiology and ultrastructure study of auditory hair cell ribbon synapses, we used hair cells and afferent fibers in the midfrequency (350–550 Hz) region of bullfrog amphibian papillae (Fig. 1A, dashed line) (Lewis et al., 1982). In vivo recordings of turtle and frog auditory hair cells have described resting membrane potentials in the range of −55 to −50 mV (Crawford and Fettiplace, 1980; Pitchford and Ashmore, 1987). In five paired recordings in which our hair cells were depolarized to −54 mV, large amplitude EPSC events were observed in the afferent fiber even though the hair cell was dialyzed with a high concentration of a fast Ca2+ buffer (10 mm BAPTA; Fig. 1B). These large events occurred in a stochastic manner during the depolarization, but were absent when the hair cell potential was held at −90 mV. In previous studies, this concentration of BAPTA completely blocked release from the ribbon-type synapses of retinal bipolar cells even after strong depolarizations (von Gersdorff and Matthews, 1994; Coggins and Zenisek, 2009). How does the hair cell release machinery overcome this strong Ca2+ buffering?
By drastically limiting the spread of Ca2+ into the terminal, fast and strong Ca2+ buffering might change the exact timing of exocytosis and, as a result, alter the kinetics of individual EPSC events. However, the individual EPSC events of Figure 1B were monophasic with a sharp rise and decay time course, having the same kinetics as small mEPSCs (Fig. 1C), and the majority of the release events had large amplitudes (mean amplitude, −203 pA; Fig. 1D). This indicated that the EPSCs were multiquantal, with a quantal content of four to six vesicles, given that the average quantal size is −57 pA at an afferent fiber holding potential of −90 mV (Li et al., 2009). To determine whether Ca2+ buffering alters hair cell release kinetics by desynchronizing multiquantal release, we analyzed the kinetics of individual EPSCs obtained under several different conditions: first, with 2 mm EGTA or 10 mm BAPTA (both done with paired recordings), and second, when the hair cell was unperturbed (afferent fiber recordings of spontaneous EPSCs). With an internal solution containing 2 mm EGTA, multivesicular release began to occur for hair cell potentials more depolarized than −70 mV (Li et al., 2009). Small events (uniquantal: 0 ∼ −100 pA) and large events (multiquantal: −100 ∼ −200 pA) had similar 10∼90% rise time (0.26 ± 0.06 vs 0.27 ± 0.03 ms; n = 5 pairs; p > 0.05 for paired Student's t test) and similar exponential decay time constants (0.62 ± 0.17 vs 0.63 ± 0.16 ms; n = 5 pairs; p > 0.05 for paired Student's t test). For an internal solution containing 10 mm BAPTA multivesicular release began to occur only for hair cell potentials more depolarized than −60 mV (Li et al., 2009). Once again, the small events and large events had similar waveforms (0.27 ± 0.04 vs 0.26 ± 0.02 ms for 10∼90% rise time; 0.50 ± 0.05 vs 0.53 ± 0.06 ms for decay time constant; n = 5; p > 0.05 for paired Student's t test). Statistical analysis for EPSC events with 2 mm EGTA versus 10 mm BAPTA revealed no significant change of 10∼90% rise time or decay time constant for either small or large events (p > 0.05 for unpaired Student's t test). For spontaneous EPSCs, in which the hair cell was not patch clamped, small events had a 10–90% rise time of 0.30 ± 0.06 ms and a decay time constant of 0.58 ± 0.05 ms, which was not significantly different from large events (0.27 ± 0.04 ms rise time, and 0.56 ± 0.05 ms decay; n = 6; p > 0.05 for paired Student's t test). Therefore, the kinetics of individual afferent fiber EPSCs appears to be remarkably insensitive to the strength of Ca2+ buffering in the hair cell.
Vesicle pool sizes determined with 10 mm EGTA
Because strong Ca2+ buffering conditions in the hair cell are unable to desynchronize the large multivesicular EPSC events, multivesicular release may be mediated by a subpopulation of vesicles tightly coupled (within a few nanometers) to Ca2+ channels. To determine the size of synaptic vesicle pools and the maximal rates of exocytosis, we used strongly depolarizing pulses from −90 to −30 mV (the peak of Ca2+ current) and measured changes in the membrane capacitance (ΔCm) of the hair cell (Parsons et al., 1994; von Gersdorff and Matthews, 1994; Hallermann et al., 2003). The durations of the depolarizing pulses were varied so that we could determine vesicle pool sizes by the saturation of ΔCm increases. Different intracellular Ca2+ buffers were also used to better isolate plateaus in Cm changes. The ΔCm jumps obtained with either 2 mm EGTA (Fig. 2A) or 10 mm BAPTA (Fig. 2B) were dramatically different in size (note the fivefold difference in vertical scale). For example, the ΔCm evoked by a short 100 ms pulse with 2 mm EGTA was 60 ± 4 fF (n = 26), which was not significantly different from the ΔCm = 74 ± 11 fF (n = 16; p = 0.16) evoked by a 3-s-long pulse with 10 mm BAPTA. Despite the effectiveness of BAPTA at blocking most of the exocytosis evoked by prolonged 1- to 3-s-long depolarizing pulses, a small amount of exocytosis (∼10%) was resistant, consistent with the EPSC data of Figure 1B. Fast Ca2+ current inactivation was consistently observed immediately after the initial current peak for hair cells held at −90 mV (Cho et al., 2011), and was followed by a slow form of Ca2+-dependent inactivation, more pronounced with EGTA than BAPTA, as has been described at turtle auditory hair cells (Schnee and Ricci, 2003) (Fig. 2A).
Averaging ΔCm across cell populations might introduce systematic errors due to the variation in cell sizes because cells of different sizes may contain different numbers and types of synaptic ribbons (von Gersdorff et al., 1996; Martinez-Dunst et al., 1997). We next checked whether the size of the hair cell was correlated with the total amount of exocytosis. We plotted ΔCm evoked by 20-ms-long [n = 16; coefficient of variation (CV) = 0.34)], 100-ms-long (n = 17; CV = 0.20), and 500-ms-long (n = 22; CV = 0.22) pulses against the resting Cm of each cell, which varied from 10 to 14.8 pF. However, regardless of pulse duration, the resting Cm was very weakly correlated with ΔCm (data not shown), justifying our calculations of average ΔCm across hair cells in the midfrequency range of amphibian papilla (Smotherman and Narins, 2000).
A plot of average ΔCm against pulse duration for the 2 mm EGTA data revealed three distinct kinetic subgroups (Fig. 3A,B): the first from 2 to 50 ms (Fig. 3B, left), the second from 100 ms to 1 s (Fig. 3B, middle), and the third from 1 to 3 s (Fig. 3B, right). The data from each of these groups were well described by a linear function (R2 ≥ 0.97; Fig. 3B, black dashed lines). The slope of the linear fit in the first subgroup (slope, 0.8 fF/ms) was the steepest, the second subgroup had a smaller slope (=0.2 fF/ms), and the third subgroup showed the least slope (=0.1 fF/ms), suggesting that transmitter release kinetics evoked by short pulses were eightfold faster than those by longer pulses. However, the capacity for fast exocytosis lasted only for ∼50 ms, whereas the slower forms of exocytosis could be sustained for longer periods. In contrast, plotting ΔCm against total Ca2+ charge (the integral of Ca2+ currents after leak subtraction) revealed only two distinct subgroups (Fig. 3C), although the relationship between ΔCm and Ca2+ charge in both subgroups was again well described by a linear fit (R2 ≥ 0.98; Fig. 3C, lines). The first group ranged from 2 to 50 ms (Fig. 3C, no. 1), while the second group ranged from 100 ms to 3 s (Fig. 3C, no. 2). Assuming that the rate of exocytosis is limited by the rate of Ca2+ influx, the slow Ca2+-dependent inactivation component of the Ca2+ current (ICa) during long depolarizing stimuli (Fig. 2A) may partially account for the lack of a third subgroup. However, the fast ICa inactivation was too rapid to account for the distinction between subgroups 1 and 2. Distinct kinetic subgroups with different rates of release suggested that there were functionally distinct vesicle pools. However, determining the size of these putative vesicle pools directly was difficult because the ΔCm did not plateau (Fig. 3A, black data points).
We next attempted to isolate the vesicle pools by restricting the spread of the synaptic Ca2+ microdomains with higher concentrations of a slow internal Ca2+ buffer (10 mm EGTA; Fig. 3A,B, red), which may slow down the recruitment rates of vesicles at ribbon synapses (Mennerick and Matthews, 1996; Babai et al., 2010). The first two subgroups then showed clear single-exponential relationships with plateaus in ΔCm plotted against pulse duration (Fig. 3B, red in left and middle panels). The asymptote (t → +∞) of the first subgroup (5–50 ms) was 29 fF, and the data were well fit by ΔCm = 29 fF · (1 − exp(−t/τ)) with τ = 23.4 ms (R2 = 0.99). The asymptote of the second subgroup (100 ms to 1 s) was 138.7 fF and the average ΔCm = 138.7 fF · (1 − exp(−t/τ)) with τ = 323 ms (R2 = 0.99). The third subgroup (1–3 s; Fig. 1B, red in right panel) showed a nonsaturating linear relationship with the same slope as in 2 mm EGTA. Surprisingly, we found that 10 mm EGTA did not block the vesicle release significantly except for four data points (50, 100, 500, and 1000 ms; Fig. 3B, green asterisks). This contrasts starkly with results from goldfish bipolar cells and mouse inner hair cells in which 5–10 mm EGTA was very effective at blocking the sustained component of exocytosis evoked by long depolarizing pulses (Moser and Beutner, 2000; Coggins and Zenisek, 2009). With 10 mm BAPTA, however, ΔCm decreased significantly for all data points by 80–90% (Fig. 3A,B, blue). To confirm that the ΔCm plateaus with 10 mm EGTA were not due to differences in Ca2+ influx, we compared the Ca2+ influx charges (integral of ICa) in 2 and 10 mm EGTA (Fig. 3D). The Ca2+ charges in 10 mm EGTA were not significantly different from those in 2 mm EGTA for all the depolarizing pulse durations.
Next, we determined the efficacy of exocytosis among the three subgroups identified in Figure 3B, calculated as a ratio of the average ΔCm per Ca2+ charge influx (or per number of Ca2+ ions) during the depolarizing pulse (Mansvelder and Kits, 1998; Hull et al., 2006). The efficacy for exocytosis was initially high for both 2 mm EGTA (Fig. 3E, black) and 10 mm EGTA (Fig. 3E,F, red) internal solutions, but decreased as depolarization duration increased, presumably due to vesicle pool depletion. With 10 mm EGTA, the efficacy decreased sharply between different subgroups compared with changes within one subgroup (Fig. 3F, arrows), suggesting greatly reduced efficacy after depletion of the first pool and a further reduction in efficacy after depletion of the second pool.
Morphological characterization of synaptic vesicle populations
To provide an anatomical correlate for our physiological estimates of vesicle pools, we inventoried synaptic vesicle populations using ultrastructural methods, assuming that the preferential site of exocytosis is located near synaptic ribbons (Zenisek et al., 2000; Zampighi et al., 2011). We first quantified the number and size of ribbon synapses per hair cell by serial section electron microscopy (Fig. 4). Four hair cells from two bullfrogs were entirely sectioned (two examples shown in Fig. 4A), and five additional hair cells were partially sectioned from another bullfrog. For our entirely sectioned cells, we counted a mean of 55.5 ribbon synapses per hair cell, a mean ribbon diameter of 208.8 nm, and the CV (a ratio of the SD to the mean) of the mean ribbon diameter at 9.8% (Fig. 4C). Additional estimates from the five partially sectioned cells were consistent with these numbers (mean ribbon number, 51.4; mean diameter, 218.7 nm; mean CV, 9.9%). Auditory nerve fibers were also inventoried for three of the entirely sectioned hair cells. These hair cells had a total of 53, 54, and 65 ribbons, and each hair cell had one large claw-like afferent receiving input from 44, 31, and 27 ribbons, respectively (Fig. 4A,B, brown fiber). Other ribbons appeared to be presynaptic to additional fibers (Fig. 4B, green fiber); however, the limited spatial dimensions of our reconstructions made it impossible to determine whether these fibers are distinct or arise from the branching of a shared fiber outside of the reconstructed volume.
In addition, for the two segmented hair cells (shown in Fig. 4A), the nearest-neighbor distance between ribbon synapses was calculated in three-dimensional space. A minority of synapses exhibited two ribbons sharing the same active zone (7 of 111 synapses). Treating these ribbon doublets as single ribbon synapses, we found the average distance between ribbon synapses to be 1.71 μm (SD = 0.65) and 1.34 μm (SD = 0.38), similar to the distance measured in mammalian inner hair cells (Meyer et al., 2009). By considering the basolateral membrane surface area at the synaptic pole of the two hair cells where ribbons were counted, we calculated the density of ribbon synapses at 0.158 ribbon synapses/μm2 and 0.117 ribbon synapses/μm2. Individual ribbon synapses were thus well separated from each other and may operate as independent units under strong Ca2+ buffering conditions (Issa and Hudspeth, 1996).
We next defined the relationship between the diameter of a ribbon and its vesicle pools using serial sections (Fig. 5A) and electron tomography (Fig. 5B). A total of 34 synapses (22 normal-external Ca2+ fixed, 12 zero-external Ca2+ fixed) with ribbons of varying diameters were entirely reconstructed and quantified (Fig. 5C). In addition to serial section reconstructions, nine ribbon synapses (four normal-Ca2+ fixed, five zero-Ca2+ fixed) were reconstructed by electron tomography. Each electron tomography reconstruction contained approximately one-half of the associated vesicle pools of the ribbon (mean, 49.2% reconstructed; n = 9), allowing us to estimate the full vesicle pools based on the percentage reconstructed.
Polynomial fits of the relationship between synaptic ribbon diameter and the size of vesicle pools (Fig. 5C, fit lines) allowed us to estimate whole-cell vesicle pools for our serial sectioned hair cells. For our four fully reconstructed hair cells, the IRP is estimated to be 702 ± 100 (mean ± SD), or 12.6 vesicles per ribbon, and the RAP is estimated to be 4308 ± 613, or 77.6 vesicles per ribbon. Similarly, for our tomography data sets, the average number of vesicles in the IRP was 14.8, with an average of 6.8 vesicles touching the membrane. These results agree well with synaptic ribbon reconstructions from the caudal (high frequency) amphibian papilla of Rana pipiens frogs where synaptic ribbons varied from 100 to 300 nm in diameter and had 59–81 attached vesicles (Simmons et al., 1995).
The correlation between morphologically and functionally defined vesicle pools
Determining vesicle pool sizes from whole-cell capacitance recordings required an estimate of the capacitance of a single synaptic vesicle. The Cm of a single synaptic vesicle (Cv) in frog saccular hair cells has been reported as 37 aF (vesicle diameter, 34.3) (Lenzi et al., 1999) and has been determined to be 45 aF for bullfrog hair cells from a small sample of vesicles (n = 109 vesicles) (Li et al., 2009). Here, we calculated a Cv of 34 aF (Cv =πd2Cs) using our mean vesicle diameter (d) of 32.9 nm (n = 1467 vesicles; Fig. 5D) and a specific membrane capacitance (Cs) of 9 mF/m2 (Fettiplace et al., 1971; Albillos et al., 1997). Accordingly, the first pool of rapidly releasing vesicles in our physiological recordings, with a size of 29 fF, consisted of ∼725–850 vesicles depending on whether Cv = 40 aF (average from previous reports) or 34 aF. This functional estimate of the IRP correlated well with the morphological IRP, which had an average size of 702. The rate of release of this pool was ∼0.8 fF/s, or 360–420 vesicles · s−1 · ribbon−1. The second pool of vesicles in our recordings (139 fF) consisted of ∼3500–4090 vesicles and was released at a slower pace. It correlated well with the morphological RAP, which had an average size of 4308 vesicles. To summarize, our average synaptic ribbon had ∼12.6 docked vesicles and a fast releasing pool of 13–15 vesicles that could be completely released by a strong depolarization within ∼20–50 ms (Fig. 3B). In comparison, frog saccular hair cells have been reported as having a larger synaptic ribbon (diameter of ∼400 nm) and a larger pool of 43–66 docked vesicles (Rutherford and Roberts, 2006).
A central difficulty in relating EPSCs to ultrastructure concerns the partially overlapping concepts of the IRP (immediately releasable pool), “docked” vesicles, and vesicles physically touching the presynaptic membrane. Our capacitance recordings (Fig. 3) showed a first wave of rapid release lasting 20–50 ms, presumably allowing IRP vesicles not initially docked (but very close to the membrane) enough time to dock and then release. For the nine tomography data sets used, the average number of vesicles in the <20 nm IRP group was 14.8 (similar in size to the 12.6 vesicles reported from serial section data sets), while the “touching the membrane” IRP group was 6.8. The number of vesicles touching the membrane was therefore sufficient to account for the size of the large multiquantal EPSCs that we observed.
Variability in vesicle size and mEPSC distributions
Our measured vesicle diameter distribution was normally distributed with a mean of 32.9 nm and a CV of 0.110 (Fig. 5D). Using EM tomography, Lenzi et al. (1999) reported a CV of 0.261 for synaptic vesicle volumes, noting that for normal vesicle diameter distributions, the CV of vesicle volume: CVvol = 3 · CVdiam. So our CVvol = 0.328, consistent with the CV of single quanta mEPSC amplitudes reported previously for our hair cell synapse (CV = 0.29) (Li et al., 2009). The variability of mEPSC amplitudes may thus originate mostly from the variability of vesicle volumes.
Vesicle mobility and sustained exocytosis
We observed sustained rates of exocytosis of 100 fF/s after 1- to 3-s-long depolarizing pulses (Fig. 3). Assuming a single vesicle capacitance of ∼34 aF, this rate of sustained release corresponds to 53 vesicles per second for each of the 55.5 ribbons of the hair cell. Can the rate at which freely diffusing cytoplasmic vesicles hit the ribbon account for this exocytosis rate? If we assume that only diffusion transports vesicles to the synaptic ribbon, where they are immediately released, and we ignore cytoskeletal and plasma membrane barrier effects, the hit rate of vesicles (F) bombarding a spherical ribbon is as follows: F = 4πDRC, where D = 1.5 × 10−2 μm2 s−1 (synaptic vesicle diffusion coefficient) (Holt et al., 2004), R = 147 nm (sum of synaptic ribbon radius, 30 nm tethers, and vesicle radius), and C is the concentration of neighboring cytoplasmic vesicles (Jackson, 2006; Rutherford and Roberts, 2008). To estimate the reserve pool of cytoplasmic vesicles at “actively exocytosing” ribbon synapses, five serial section data sets from the “Ca2+-fixed” condition were arbitrarily chosen (see Materials and Methods). The concentration of cytoplasmic vesicles (within ∼1 μm of the ribbon) ranged between 2.55 and 4.78% of maximum vesicle packing, averaging 3.65%, so that C ≈ 851 vesicles/μm3. We calculate a collision rate of 24 vesicles per second for each ribbon. Diffusion alone, therefore, appears to be inadequate to maintain even the slowest rate of sustained release observed at these synapses. Consequently, an active ATP-dependent, and probably Ca2+-dependent, replenishment mechanism seems to be required to enhance vesicle mobility and concentrate vesicles near synaptic ribbons (Heidelberger et al., 2002).
Physiological estimate of the number of Ca2+ channels
To determine the total number of Ca2+ channels per hair cell, we used nonstationary noise analysis of Ca2+ tail currents (Roberts et al., 1990). Figure 6A shows the Ca2+ current I–V relationship (n = 5 cells) (Edmonds et al., 2004). Figure 6B shows the voltage-clamp protocol we used to elicited Ca2+ tail currents. By plotting the variance of this tail current against its mean (Fig. 6C,D), the number of Ca2+ channels and their single-channel current was extrapolated by fitting the data with a parabolic function (see Materials and Methods) (Alvarez et al., 2002). We estimated a total number of Ca2+ channels NCa = 1998 ± 355 (n = 6 hair cells) and a single-channel current i = −0.62 ± 0.19 pA at −90 mV. Assuming the reversal potential ECa = +90 mV, approximately equal to the thermodynamic ECa, this yields a single-channel conductance of 3.4 pS. At −30 mV, the peak of the Ca2+ current IV curve, we calculated that i = −0.4 pA. This single-channel current value is very similar to that found in ribbon synapses of inner hair cells and retinal rod bipolar cells (Brandt et al., 2005; Jarsky et al., 2010). Estimates of number of Ca2+ channels in mouse inner hair cells and frog saccular hair cells are 1700 and 1800, respectively (Roberts et al., 1990; Brandt et al., 2005).
An alternative method to estimate the number of Ca2+ channels uses the resonant frequency of nonmammalian hair cells (Wu and Fettiplace, 1996; Martinez-Dunst et al., 1997). For our amphibian papilla hair cells, at the average steady-state membrane potential Vss = −49.6 ± 0.4 mV, the average resonant or tuning frequency (fc) was 410 ± 2 Hz (n = 15 hair cells), and when Vss = −54.8 ± 0.3 mV, the fc = 353 ± 2 Hz (n = 12). The number of Ca2+ channels NCa = 2 · (NKCa) and the number of Ca2+-activated K+ channels NKCa = 3.75 · fc (Wu and Fettiplace, 1996). Using this method, we calculate that NCa varies from 2648 to 3075. We thus estimate from two independent data sets that the number of functional Ca2+ channels contained in our hair cells ranges from ∼2000 to 3000. Given our average number of 55.5 ribbons per hair cell, we estimate that the average number of Ca2+ channels per ribbon is ∼40–50, assuming that all hair cell Ca2+ channels are clustered near to ribbons (Brandt et al., 2005; Sheets et al., 2011).
Monte Carlo modeling of Ca2+ influx and propagation under the synaptic ribbon space
One hypothesis of how multiquantal EPSC events can still occur in 10 mm BAPTA conditions (Fig. 1B) is that the restricted, small space available under a synaptic ribbon permits the rapid saturation or depletion of Ca2+ buffers near open Ca2+ channels (Roberts, 1994). We therefore asked how the kinetics of Ca2+ influx and Ca2+ buffering influence the local concentration of free Ca2+ ions in the crowded space under a synaptic ribbon.
To examine the changes in free Ca2+ concentration following the onset of Ca2+ channel opening, we used Monte Carlo diffusion simulations incorporating the detailed spatial parameters of our ribbon synapse reconstructions as well as the Ca2+ buffer concentrations, channel numbers, and Ca2+ currents of our physiology experiments (Fig. 7A). To simulate Ca2+ currents resulting from a depolarizing pulse from −90 to −30 mV, we opened 22 or 23 Ca2+ channels of a total of 45 located under the ribbon (for a Ca2+ current of approximately −9 pA, using −0.4 pA as the single Ca2+ channel current). Simulations were run in either 2 mm EGTA (Fig. 7B) or 10 mm BAPTA (Fig. 7C) buffering conditions while monitoring free Ca2+ concentration at eight different locations (Fig. 7A, red boxes). In 2 mm EGTA conditions, free Ca2+ reached stable concentrations rapidly (∼100–200 μs; Fig. 7B), with the region underneath the ribbon experiencing uniformly large (millimolar) concentrations. However, with increasing distance from the base of the ribbon (and into less encumbered cytoplasmic space), free Ca2+ concentrations diminished sharply. Simulations in 10 mm BAPTA conditions (Fig. 7C) revealed dramatic reductions in free Ca2+ concentrations, particularly around the sides of the ribbon. Although the slower rise of Ca2+ with stronger buffer suggests a longer EPSC latency for our physiological recordings, variability between synapses in the number of open Ca2+ channels and size of vesicle pools, along with variable numbers of ribbon synapses per afferent fiber, made analysis of EPSC latency difficult. These simple simulations likely overestimated the local Ca2+ concentration near the channels, as we did not take into account the dramatic alterations of the Ca2+ electrochemical driving force that such high local Ca2+ concentration would cause. However, despite this overestimation, rises in local Ca2+ concentrations were still confined to the bottom row of vesicles in 10 mm BAPTA conditions.
Since multivesicular release persisted despite the 10 mm BAPTA presynaptic Ca2+ buffering and at relatively low hair cell membrane potentials (Fig. 1B), it was not clear how the spatial spread of free Ca2+ related to synaptic vesicles when fewer channels were open (Fig. 8). In 2 mm EGTA, large Ca2+ microdomains encompassing several vesicles arose from the opening of a single channel. Comparable peak Ca2+ concentration could be obtained in 10 mm BAPTA by the opening of six channels. For both of these conditions, we found a dramatic enhancement of the Ca2+ microdomain when the ribbon was considered to be an opaque diffusion barrier and restrictor of presynaptic space.
Previously, we described a sharp transition from sparse mEPSC events that are uniquantal (quantal size, −57 pA) to abundant multiquantal EPSC events (amplitudes >100 pA) as hair cells were progressively depolarized during paired recordings [Li et al. (2009), their Fig. 5]. An example of these previous data is shown in Figure 9, A and B, for hair cells dialyzed with 2 mm EGTA and 10 mm BAPTA, respectively. Note that a small 3 mV increase in hair cell depolarization caused the amplitude of the EPSCs to suddenly increase in size. With 2 mm EGTA, this transition point occurred at approximately −70 mV, where the macroscopic ICa < −10 pA (Fig. 6A, inset). From ICa = i · Po · NCa (where i = −0.54 pA at −70 mV, NCa = 2000, and Po is the single-channel open probability), we estimated that Po < 0.009 at −70 mV. With 40–50 Ca2+ channels per ribbon, the probability of two or more opening simultaneously to trigger a multiquantal release event was very low. However, with 10 mm BAPTA, the transition to multiquantal EPSC events occurred only at potentials more depolarized than −60 mV, when Ca2+ channel open probability was considerably higher. These results suggest that one open Ca2+ channel may be sufficient to trigger multivesicular release under weak Ca2+ buffering conditions (2 mm EGTA), whereas several Ca2+ channels are required to open before a multivesicular release event can occur under strong Ca2+ buffering conditions (10 mm BAPTA).
To test the hypothesis that a local depletion of Ca2+ buffers by open Ca2+ channels underlies the transition to multiquantal EPSC events in 10 mm BAPTA conditions (Fig. 1B), we resorted again to our Monte Carlo simulations. We monitored the concentrations of free Ca2+ and BAPTA in a 100-nm-diameter region under the ribbon (Fig. 9C,D), while varying the number of open Ca2+ channels. In addition, the same simulations were repeated for three different ribbon models: (1) a model incorporating vesicles numbers consistent with our reconstructions (Fig. 9E, black inset); (2) a model incorporating maximum vesicle packing (Fig. 9E, red inset); and (3) a model incorporating vesicles numbers consistent with our reconstructions, but where the ribbon was not a barrier to diffusion (Fig. 9E, blue inset). For the model closely resembling our reconstructed ribbons, opening less than ∼11 channels resulted in relatively modest increase in free Ca2+ concentration (Fig. 9E, black trace). However, large increases in free Ca2+ concentration resulted from opening >11 channels, suggesting a transition point reminiscent of EPSC recordings. Consistent with the idea that tighter packing of ribbon-associated vesicles excludes more presynaptic free volume, thereby reducing local Ca2+ buffering capability, the model incorporating maximum vesicle packing showed a shift in the transition point to approximately five channels (Fig. 9E, red trace). In contrast, the model lacking a ribbon diffusion barrier showed relatively modest increases in free Ca2+ across the entire range of channel openings that were simulated (Fig. 9E, blue trace). Our Monte Carlo simulations thus suggest that, by acting as a “diffusion barrier” to incoming Ca2+ ions, the synaptic ribbon may function as a local Ca2+ concentration booster that promotes synchronous multivesicular release.
Discussion
How does one account for the unusually large amplitudes of individual EPSC events at hair cell synapses? One potential explanation is that presynaptic membrane-docked vesicles are released in a coordinated, synchronous fashion. Whether this coordination is mediated by tight colocalization of vesicles with Ca2+ channels, by very large presynaptic Ca2+ microdomains that might encompass multiple vesicles, or even by mechanical coupling of vesicles remains unclear. Alternatively, large EPSC events could be mediated by the compound fusion and coalescence of several vesicles tethered to the ribbon before the final release event (Parsons and Sterling, 2003; Edmonds et al., 2004). Our results favor a synchronous release mechanism over a two-step compound fusion mechanism in three main ways. First, the kinetics of EPSCs remained fast across presynaptic Ca2+ buffering conditions, and even with concentrations of BAPTA that severely restrict the spread of free Ca2+ away from the presynaptic membrane. Second, our anatomical and physiological inventories of synaptic vesicle populations in hair cells revealed a close correlation between the population of ribbon-attached vesicles most proximal to the presynaptic membrane and the size of the first, ultrafast wave of vesicle release (Mennerick and Matthews, 1996). Third, our Ca2+ diffusion simulations, incorporating the volume-excluding effect of the ribbon and its associated vesicles, demonstrate that small differences in the number of open Ca2+ channels may result in large and rapid increases in free Ca2+ concentration under the ribbon. As such, the influence of presynaptic structures on Ca2+ dynamics at the base of the ribbon may facilitate a synchronous release mechanism and account for the sharp transition to, and fast kinetics of, the large EPSCs.
Compound fusion and EPSCs at hair cell synapses
Recently, compelling EM evidence for compound fusion was presented for strongly stimulated bipolar cell ribbon synapses (Matthews and Sterling, 2008). However, compound fusion at other synapses exhibits large EPSCs with slower rise times compared with mEPSCs (He et al., 2009). Assuming both single and compound vesicles form the same narrow fusion pore with the plasma membrane, the slower rise time of large EPSCs likely results from a compound vesicle requiring more time to fully void its transmitter content than a single synaptic vesicle. Indeed, the amount of glutamate released from a vesicle, Cr in time t is described by the following: Cr = Co · [1 − exp(−t/τv)], where Co is the amount initially in the vesicle, and τv = VL/(πro2D), where V is the vesicle volume, L is the pore length, ro is the pore radius, and D is the diffusion coefficient of glutamate in the pore (Almers et al., 1991). To account for the fast rise times (<0.2–0.3 ms) of large EPSCs in our experiments, compound fusion would require at least two processes: first, homotypic fusion of vesicles along the upper rows of a ribbon should occur in a Ca2+-independent manner (since these vesicles are located >40–100 nm from Ca2+ channels and are bathed by 10 mm BAPTA; Fig. 7); second, to avoid a slow EPSC rise time, Ca2+-dependent exocytosis of a large, compound vesicle should produce a rapidly expanding fusion pore (Stiles et al., 1996). Since synaptic vesicles are not known to fuse with each other in a Ca2+-independent manner, and previous evidence for compound exocytosis does not imply fast fusion pore expansion rates (Van der Kloot and Naves, 1996; He et al., 2009), the classically defined compound fusion mechanism appears implausible for hair cell synapses.
A crowded, dynamic presynaptic space
That the ribbon may serve as a diffusion barrier to Ca2+ ions is supported by estimates that a large fraction (>67%) of the volume of the ribbon is occupied by a tightly packed protein called RIBEYE (Zenisek et al., 2004), which constitutes only one of several different protein fractions present in the ribbon (Magupalli et al., 2008). This suggests that the dense matrix of proteins forming the ribbon has little, if any, open space to permit free Ca2+ diffusion in its interior. Computer modeling of Ca2+ influx at frog saccular hair cells demonstrated that a large ribbon (diameter, ∼400 nm) that acts as a diffusion barrier causes a dramatic enhancement in the steady-state Ca2+ concentration and local depletion of unbound Ca2+ buffer in the space between the ribbon and the presynaptic membrane (Roberts, 1994). In addition, docked synaptic vesicles can also act as diffusion barriers that distort the single Ca2+ channel nanodomains, resulting in large concentration differences across the surface of the vesicle (Shahrezaei and Delaney, 2004). Combining these two putative diffusion barriers into a single model based on our ultrastructural characterizations, we reach similar conclusions.
Electron tomography of frog sacculus ribbon synapses following inhibition of release in low external Ca2+ conditions showed an increase in the vesicle packing density around and underneath the ribbon (Lenzi et al., 2002). For our hair cells, which are held at −90 mV before release, vesicles are likely at higher packing densities as well. Under these conditions, ribbons fully packed with vesicles presumably maximize the release probability of large, fast multivesicular events in several ways: first, minimizing the available volume for Ca2+ diffusion serves to maximize the increase in Ca2+ concentration for a given Ca2+ current; second, multiple docked vesicles in the vicinity of a channel experience this maximized concentration (Fig. 8); and third, the fusion of a vesicle cohort suddenly creates extra diffusional volume at the active zone, temporarily decreasing the release probability until new docked vesicles can repopulate the available space next to the Ca2+ channels. Importantly, simulations with 10 mm BAPTA in the cytoplasm reveal that high Ca2+ levels exist only near docked vesicles, where unbound BAPTA is severely depleted near open Ca2+ channels. Once the BAPTA is depleted, incremental increases in the number of open Ca2+ channels can result in large increases in local Ca2+, consistent with the sharp transition from mEPSCs to large EPSCs that we observe across a range of only a few millivolts (Fig. 9A). Further studies of the voltage-dependent flicker open probability of single Ca2+ channels are needed to better understand this transition from single quanta to multiquantal release (Zampini et al., 2010).
At the cellular level, enhancement of local Ca2+ concentrations by ribbons with full pools of docked vesicles might account for our observations of initially high efficacy of exocytosis, decaying dramatically with depletion of the releasable pool (Fig. 3F). Accordingly, experiments using the flash photolysis of caged-Ca2+ in bipolar and inner hair cells show that the fast component (rise time, <1 ms) of exocytosis requires large elevations of presynaptic Ca2+ to ∼100 μm (Heidelberger et al., 1994; Beutner et al., 2001). Our in vitro recordings and in silico simulations thus suggest that docked vesicle fusions mediate the large and fast EPSC events. Furthermore, the rate-limiting step for ongoing release during continuous hair cell stimulation may not be exocytosis itself, but the Ca2+-dependent recovery rate from depletion of the vesicle pool docked underneath the synaptic ribbon (Furukawa et al., 1982; Cho et al., 2011; Schnee et al., 2011).
Vesicles pools and synaptic ribbon geometry
Recent work has begun to reveal some major differences between ribbon synapses of different sensory neurons (LoGiudice and Matthews, 2009). For instance, the large sheet-like ribbons of cone photoreceptor terminals are capable of acting as “capacitors” that release a large bout of vesicles after being fully “charged” with vesicles in hyperpolarized conditions (Jackman et al., 2009). The synaptic ribbon of retinal bipolar cells is a relatively small plate-like structure, so vesicles do not dock underneath the ribbon, but along its lateral ridge. Accordingly, only 30% of spontaneous EPSCs are multivesicular, and these have a quantal content of only two vesicles (Singer et al., 2004; Jarsky et al., 2010). In contrast, the majority of the spontaneous EPSCs at bullfrog and rat inner hair cells are multiquantal (Li et al., 2009; Grant et al., 2010). Unlike mammalian inner hair cells where a single ribbon synapse apposes a single, dedicated afferent fiber, our afferent fiber reconstructions indicate that dozens of ribbon synapses often feed into a single claw-like afferent fiber (Fig. 4B). Is there a functional advantage to this arrangement? One idea originates from the finding that a pool of docked vesicles tightly coupled to the Ca2+ channel cluster beneath the ribbon reduces first spike latency and jitter in the auditory nerve (Buran et al., 2010; Frank et al., 2010). Likewise, a large pool of docked vesicles distributed among dozens of ribbon synapses with perhaps different release probabilities may account for the extremely short spike latency and low jitter of some frog afferent fibers (Wittig and Parsons, 2008). Moreover, we argue that, after a large EPSC event, the synaptic ribbon may become refractory to release until it is fully loaded again, so an afferent fiber with multiple ribbon synapses may be better able to sustain a prolonged high rate of large EPSC events. We thus propose that the coordination of several docked vesicles underneath a spherical ribbon endows hair cell synapses with an unmatched ability for large and fast-rising EPSC events that are well suited to trigger phase-locked spikes at high stimulation frequencies.
Footnotes
This work was supported by NIDCD Grant DC04274 (H.v.G.), a Deafness Research Foundation Grant and a Tartar Fellowship (S.C.), a K99 Research Award from NIDCD (G.-L.L.), and the NIDCD Division of Intramural Research (C.W.G., B.K.). We thank Will Grimes and Gary Matthews for discussions.
- Correspondence should be addressed to either of the following: Bechara Kachar, National Institute on Deafness and Other Communication Disorders, National Institutes of Health, 50 South Drive, MSC 8027, Bethesda, MD 20892-8027, kacharb{at}nidcd.nih.gov; or Henrique von Gersdorff, The Vollum Institute, Oregon Health & Science University, 3181 Southwest Sam Jackson Park Road, Portland, OR 97239, vongersd{at}ohsu.edu