Properties of multivesicular release from mouse rod photoreceptors support transmission of single-photon responses
Abstract
Vision under starlight requires rod photoreceptors to transduce and transmit single-photon responses to the visual system. Small single-photon voltage changes must therefore cause detectable reductions in glutamate release. We found that rods achieve this by employing mechanisms that enhance release regularity and its sensitivity to small voltage changes. At the resting membrane potential in darkness, mouse rods exhibit coordinated and regularly timed multivesicular release events, each consisting of ~17 vesicles and occurring two to three times more regularly than predicted by Poisson statistics. Hyperpolarizing rods to mimic the voltage change produced by a single photon abruptly reduced the probability of multivesicular release nearly to zero with a rebound increase at stimulus offset. Simulations of these release dynamics indicate that this regularly timed, multivesicular release promotes transmission of single-photon responses to post-synaptic rod-bipolar cells. Furthermore, the mechanism is efficient, requiring lower overall release rates than uniquantal release governed by Poisson statistics.
Introduction
One of the most impressive features of the visual system is the ability to detect single photons. Pioneering psychophysical studies showed that humans can detect flashes consisting of a few photons hitting the retina, indicating that individual rods can respond to the absorption of single photons and then reliably signal these events to post-synaptic neurons (Barlow, 1956; Hecht et al., 1942; Sakitt, 1972; Tinsley et al., 2016). The ability of rods to respond to single photons was subsequently confirmed by electrophysiological recordings (Baylor et al., 1984; Gross et al., 2015; Reingruber et al., 2015; Rieke and Baylor, 1998). In darkness, rods maintain a relatively depolarized membrane potential of around −40 mV, and the absorption of a photon produces a small hyperpolarization, between 1 and 3.5 mV (Cangiano et al., 2012; Hornstein et al., 2005; Schneeweis and Schnapf, 1995). Rods signal the absorption of a photon to post-synaptic rod-bipolar cells by decreasing the rate of glutamate release. However, synaptic vesicle release is an intrinsically noisy process that is typically described by Poisson statistics (Malagon et al., 2016; Miki, 2019; Zhang and Peskin, 2015), posing a problem for post-synaptic bipolar cells to distinguish a genuine decrease in vesicle release caused by the absorption of a photon from a stochastic pause in vesicle release. How does the rod transform a small voltage change into a sufficiently large and reliable change in vesicle release?
One proposed solution is for a rod to maintain a high basal rate of vesicle release. Under the assumption of Poisson release, a rod would need to maintain a basal release rate of ~100 vesicles/s to generate a sufficiently reliable signal to distinguish single-photon responses from stochastic variability in vesicle release (Rao-Mirotznik et al., 1998; Rao et al., 1994; Schein and Ahmad, 2005; van Rossum and Smith, 1998). This solution requires an extremely high turnover of vesicles that is energetically and physically demanding. Furthermore, measurements from salamander, gecko, and mouse rods show basal release rates of only 10–20 vesicles/s (Hays et al., 2020b; Sheng et al., 2007).
A second proposed solution is to make release more regular than a Poisson process. For example, implementing an Erlang process that releases vesicles at regular intervals would reduce the likelihood of mistaking a stochastic decrease in vesicle release for the absorption of a photon (Schein and Ahmad, 2005). Erlang events are generated by accumulating a certain number of Poisson events, which reduces interevent variability to regularize release. While Poisson statistics accurately describe vesicle release at many central synapses, this may not necessarily be true at the ribbon-style synapses in rods. Ribbons are plate-like, pre-synaptic structures that tether multiple vesicles along the surface and it has been suggested that delivery of vesicles down a ribbon may help to make release more regular (Schein and Ahmad, 2005).
To identify how single-photon responses can be signaled by changes in vesicle release from rods, we recorded pre-synaptic anion currents associated with glutamate transporter activity in mouse rod terminals (Hays et al., 2020b). The glutamate transporters in rods (largely EAAT5) are linked to an uncoupled anion conductance, so that an anion channel opens as glutamate is retrieved (Arriza et al., 1997; Schneider et al., 2014). Glutamate reuptake into rods continues even after post-synaptic mGluR6 receptors in rod-bipolar cells become saturated, suggesting that transporters are not saturated during rod release (Hasegawa et al., 2006). In that range, glutamate anion currents (IA(glu)) vary linearly with glutamate levels in the synaptic cleft and can thus provide a pre-synaptic measure of glutamate release (Hasegawa et al., 2006; Otis and Jahr, 1998). Using these techniques, we observed stochastic release of individual vesicles that obeyed Poisson statistics when rods were voltage-clamped at −70 mV. However, at the typical resting membrane potential in darkness of −40 mV, release transitioned to regularly spaced multivesicular release events of ~17 vesicles apiece (Hays et al., 2020b). We found that the probability of initiating multivesicular release events was nearly abolished by a small membrane hyperpolarization similar to single-photon responses. A model of the empirically determined, voltage-dependent vesicle release dynamics combined with an ideal observer analysis reproduced the detection sensitivity of mouse rod-bipolar cells (Field and Rieke, 2002a). These results show that regularly timed multivesicular release events in rods can enhance post-synaptic detection of small voltage changes evoked by the absorption of individual photons. Furthermore, this strategy is efficient, achieving greater sensitivity than Poisson release while requiring fewer vesicles per second.
Results
To study voltage-dependent vesicular glutamate release from rods, we recorded IA(glu) from rods using flatmount mouse retina preparations (Figure 1). When rods were held for many seconds at −60 or −70 mV, similar to the membrane potential achieved in bright light, we observed occasional inward currents arising from activation of glutamate transporters on the rod terminal (Figure 1; Grassmeyer et al., 2019). As described earlier (Hays et al., 2020b), these inward currents had waveforms typical of quantal post-synaptic currents but with a slow time course (10–90% rise time of 9–10 ms with decay time constants of ~40 ms). Spontaneous inward currents in rods voltage-clamped at −70 mV showed a unimodal amplitude distribution suggesting they consisted entirely of uniquantal events (Figure 1C). The overall rate of release increased at depolarized membrane potentials but achieved a rate at room temperature of only 4–5 vesicles/s at −40 mV, similar to the membrane potential in darkness. As shown previously, release rates and ICa were doubled by increasing the bath temperature to 35°C, attaining an average rate of ~11 v/s/ribbon at −40 mV (Hays et al., 2020b). Accompanying this increased rate of vesicle release was a switch to more coordinated forms of release with semi-regular multivesicular release events dominating release in rods held at −40 mV (Figure 1A,B). Occasional unitary events were interspersed among these multivesicular events and so the overall distribution of event amplitudes measured at −40 mV was bimodal (Figure 1C). As described previously (Hays et al., 2020b), multivesicular release events involved the rapid sequential release of 10–20 vesicles triggered by Ca2+ influx and derived from a pool of ribbon-associated vesicles shared with evoked release.
Multivesicular release was not modeled by Poisson statistics
It is generally accepted that under conditions of low release probability (e.g., when a neuron is hyperpolarized), vesicle exocytosis occurs stochastically, obeying Poisson statistics (Malagon et al., 2016; Miki, 2019; Zhang and Peskin, 2015). Spontaneous release can arise from Ca2+ independent fusion events or from chance openings of voltage-gated Ca2+ channels (Cork et al., 2016; Kavalali, 2015; Kavalali, 2020). The persistence of spontaneous events in rods held at −70 mV, even after blocking Ca2+ channels with extracellular Cd2+ or strongly buffering intracellular Ca2+ with 10 mM BAPTA, showed evidence for Ca2+-independent release (Hays et al., 2020b). To analyze frequency distributions, only rods in which we observed >50 events were included in the data set. The frequency distributions of interevent intervals for release events detected in rods voltage-clamped at −70 mV (e.g., Figure 2A) consistently exhibited a single exponential decay as predicted for a Poisson release process (R2 = 0.97 ± 0.04, mean ± SD, n = 9 rods; Figure 2C). The probability of the number of release events that occur in a 1 s interval was also well fit using a Poisson probability distribution (Equation 6; Figure 2E). The Poisson fit for spontaneous events was robust (R2 = 0.94 ± 0.07, n = 9 rods), and the mean number of events per second predicted from the Poisson fit, λ = 1.06 ± 0.55 (n = 9), was very close to the actual mean number of events observed in the same sample (λ = 1.0 ± 0.6, p=0.83, n = 9).
These results indicate that the rate of Ca2+-independent spontaneous release from rods at hyperpolarized membrane potentials is a Poisson process. Does Poisson release persist when Ca2+-dependent release is engaged at more depolarized potentials? At the typical resting potential in darkness (−40 mV), we observed periodic multivesicular events along with occasional unitary events (Figure 2B). Like those observed at −70 mV, the intervals between unitary events at −40 mV were well described by a single exponential (Figure 2D, filled squares, R2 = 0.90 ± 0.07, n = 9 rods), consistent with an underlying Poisson release process. However, if we treated each multivesicular event as an individual release event, histograms of interevent intervals could not be fit by single exponentials, suggesting that multivesicular events are not released by a Poisson process (Figure 2D, bars). It therefore appears that single and multivesicular release events emerge independently.
We estimated the number of vesicles released in each event by dividing the total charge transfer during each multivesicular event by the average charge transfer for unitary events measured at the same potential in the same cell. Glutamate transporter currents continue to increase even after saturation of post-synaptic mGluR6 receptors, indicating that transporters are not normally saturated during release from rods (Hasegawa et al., 2006). This implies that the charge transfer for each IA(glu) vesicle release event should sum linearly. A Poisson distribution could not describe the frequency distribution of individual vesicle release events counted in this way. During each 1 s interval, there was either a multivesicular release event or no release event and so interevent intervals were distributed bimodally (Figure 2F). The frequency distribution could not be fit with a Poisson model (Equation 6), yielding negative R2 values (−0.7 ± 0.6, n = 9) both when λ was unconstrained (λ = 2.5 in Figure 2F) and when λ was constrained to the actual mean number of events/s (λ = 11 in Figure 2F). Thus, unlike unitary event statistics, the statistics of multivesicular release were not well described by a Poisson distribution.
Fitting the interevent intervals for multivesicular events with a Gaussian function yielded better fits (mean R2 = 0.84 ± 0.12) than fitting with an exponential function. Interevent intervals averaged 2.34 ± 0.64 s (N = 16) with a coefficient of variation (CV) among these cells of 0.43 ± 0.09. For a Poisson distribution where variance = mean, the average mean intervals between multivesicular events in this sample of cells predicted a CV of 0.67. Narrowing the interval distribution is a characteristic of an Erlang distribution that waits a certain number of Poisson intervals before each Erlang event. The number of Poisson intervals, k, that occur before each Erlang event is related to the CV by the formula:
The intervals between multivesicular events were consistent with Erlang factors ranging from 3 to 12 and averaged k = 6.2 ± 3.12 (n = 16).
Impact of experimental conditions on statistics of multivesicular release
The conditions of our experiments differ from dark-adapted conditions in vivo. We therefore investigated how different aspects of these experiments might have impacted the observed release dynamics. We began by examining the impact of dark adaptation. During dark adaptation, ribbons can increase in length in mouse rods (Dembla et al., 2020) and bursts involve release from the pool of vesicles at the base of the ribbon (Hays et al., 2020b). Thus, dark adaptation might influence the statistics of vesicle release. However, when holding dark-adapted rods at −40 mV, we found that the number of vesicles in each multivesicular release event (16.7 ± 5.8, n = 7) and the intervals between multivesicular event start times (2.0 ± 0.54 s, n = 6, room temperature) did not differ significantly from those measured in retinas prepared in the light (17 ± 7 v/burst, n = 22, p=0.74, unpaired t-test; 2.3 ± 0.64 s intervals, n = 16, p=0.34, unpaired t-test) (Hays et al., 2020b). Thus, recording from rods in light-adapted versus dark-adapted states did not strongly alter release statistics. One caveat is that in both light- and dark-adapted retinas, we mechanically removed outer segments to expose the inner segments and cell bodies prior to recording.
Next, we examined the potential impact of physiological rundown following membrane rupture by the patch pipette for voltage-clamping the rods. As illustrated in Figure 3, multivesicular events increased in amplitude during the first few minutes of recording as SCN− diffused into the rod from the patch pipette. There was also a decrease in multivesicular event frequency after patch rupture, presumably due to rundown of the release mechanism. When the intervals between multivesicular events were plotted against time after patch rupture, the slopes of the linear regressions were all significantly non-zero (n = 5 rods, room temperature). The y-intercept values averaged 1.23 ± 0.45 s, whereas the overall mean interval between multivesicular events in this sample averaged 2.26 s ± 0.60 s (n = 5, room temperature). This suggests that the interevent intervals present at patch rupture were roughly half the overall average interval between multivesicular release events yielding approximately twice the release rate.
Finally, we examine the impact of temperature. Due to more rapid rundown at higher temperatures, most of our recordings were performed at room temperature (~22°C). Elevating the bath temperature also significantly shortened release intervals, from 2.3 ± 0.64 s (N = 16) at room temperature to 1.13 ± 0.67 s (n = 5) at 35°C (p<0.0001, unpaired t-test). While rundown was more rapid at 35°C, if we simply assume a similar amount of rundown as seen at room temperature, this suggests interevent intervals at 35°C as short as 615 ms (i.e., a frequency of 1.6 Hz). Alternatively, if we assume that the most rapid rates measured at 35°C come from the healthiest cells and are thus more representative of genuine rates in vivo, then interevent intervals one standard deviation below the mean of ~465 ms suggest a rate of 2.2 Hz. These various lines of reasoning converge on a spontaneous rate of multivesicular release events in darkness at physiological temperatures from intact mouse rods of ~2 Hz.
Multivesicular release is sensitive to small voltage changes
We examined the impact of small voltage changes on the probability of multivesicular release events from rods. We hypothesized that multivesicular release may play a role in transmitting single-photon responses from rods because multivesicular release occurred most prominently in rods near the typical membrane potential in darkness and was not observed in cones. In mouse rods, the voltage changes evoked by single-photon absorption average ~3.5 mV (Cangiano et al., 2012). To test the sensitivity of multivesicular release to small voltage changes, we voltage-clamped the rod at –40 mV for 30 s and then applied a hyperpolarizing step to −43.5 mV. This small hyperpolarization consistently produced an appreciable pause in release (Figure 4A). The average interevent intervals measured over a 30 s period at −43.5 mV (3744 ± 2243, 129 intervals) were significantly longer than those measured in the same cells at −40 mV (2405 ± 1430 ms, 279 intervals, p<0.0001, unpaired t-test, n = 9 rods). In addition, the interval to the first multivesicular event observed after the step to −43.5 mV (red triangles in Figure 4B show measurements from individual trials) was consistently greater than 1 standard deviation above the mean interevent intervals measured at −40 mV. Single-photon responses of primate rods appear be smaller than those of mouse rods, averaging only ~1 mV (Schneeweis and Schnapf, 1995). We found that steps of −1 (n = 7; filled circles in Figure 4B) and −2 mV (n = 7; stars in Figure 4B) also produced appreciable pauses in multivesicular events immediately after the step. This lengthening of release intervals could potentially improve the signaling of single-photon responses to post-synaptic bipolar cells, allowing genuine reductions in release caused by a photon to be distinguished from a stochastic fluctuation in release.
To further investigate the capability of a single-photon absorption event to reduce multivesicular release, we clamped rods at –40 mV and applied a voltage waveform that mimicked the amplitude and kinetics of a single-photon voltage response in a mouse rod (Figure 5). This voltage waveform was constructed from a two exponential function that attained a peak amplitude of 3.4 mV after 245 ms and then declined with a time constant of 2.5 s (Figure 5A; Cangiano et al., 2012). In rods voltage-clamped at −40 mV, we repeated this single-photon voltage stimulus for 20–60 trials (Figure 5A). Binning the multivesicular release events in time revealed that release probability declined almost to zero at the peak of the simulated light response, falling in parallel with membrane hyperpolarization (Figure 5B; n = 1892 events from 28 rods). Multivesicular release then rebounded with an overshoot as the membrane potential recovered, showing ~50% more events initiated during the first few hundred milliseconds after termination of the simulated light response (between 3000 and 4000 ms, Figure 5B). The voltage stimulus did not truncate multivesicular events that had started prior to the stimulus; events that continued during the stimulus had the same number of quanta, measured by charge transfer, as prior events (p=0.53, t-test, n = 5 rods, not shown).
White noise analysis
To identify voltage changes that preferentially evoke multivesicular release in rods, we held rods at −40 mV and applied a white noise voltage stimulus with a standard deviation of 3.5 mV and high frequency cutoff of 5 Hz, consistent with the rod power spectrum (Chichilnisky and Rieke, 2005). The average stimulus waveform that preceded each multivesicular event, that is, the event triggered average (ETA) waveform, consisted of a small hyperpolarizing excursion followed by a larger depolarizing excursion with a time course roughly similar to termination of a single-photon response (n = 4 rods; Figure 6A). Using a white noise voltage stimulus with a smaller standard deviation (0.8 mV), the ETA showed a waveform with a similar shape and time course (Figure 6B). When we tested a similarly small stimulus but with 20 Hz cutoff, the ETA appeared as noise (Figure 6C), suggesting the rod terminal preferentially responds to slow voltage changes with frequency response characteristics of rods in light and dark, filtering out higher frequency membrane voltage noise. These data reinforce the idea that transmitter release from rods is very sensitive to small, slow voltage changes similar to those experienced during a single-photon response, consistent with a role of these multivesicular events in transmitting small light responses at the rod ribbon synapse.
Simulation of vesicle release dynamics indicates efficient transmission of the single-photon response
The previous experiments suggest the multivesicular release dynamics observed in rods held at −40 mV could be used for transmitting single-photon responses to rod-bipolar cells. We examined the consequences of these release statistics on transmitting single-photon responses by generating a model of rod release dynamics matched to those we measured and by utilizing an ideal observer analysis to examine the resulting signal fidelity (Chichilnisky and Rieke, 2005; Field et al., 2019; Smith and Dhingra, 2009).
Multivesicular release was modeled as a time-dependent Erlang process across a population of simulated rods forming the receptive field of a rod-bipolar cell (Figure 7A; Erlang release model, see Materials and methods). In darkness, the release probability was constant in time (Figure 7F); the absorption of a photon caused a suppression of release lasting T seconds (Figure 7B). Note, this will generate the same integrated change in release probability as exponential decay with a time constant of T. Following suppression, release probability rebounded to an amplitude of A above baseline and then decayed exponentially back to baseline with a time constant (Figure 7B). To generate Erlang release statistics, we first sampled an inhomogeneous Poisson process (Figure 7C). Poisson events were counted until reaching the Erlang factor, Ef. When the counter reached Ef, this generated a release event (Figure 7D) and the counter was reset. This yielded sub-Poisson variability with a coefficient of variation less than a Poisson process by 1/sqrt(Ef). Release events were convolved with a template, and independent Gaussian noise was added to approximate recorded rod release events (compare Figure 7E with Figure 2B).
The fidelity of rod output with these vesicle release dynamics was assayed using a two-alternative force-choice (2AFC) task and an ideal observer analysis (Chichilnisky and Rieke, 2005; Field et al., 2019; Smith and Dhingra, 2009). For the 2AFC task, one response was generated by simulating a flash of strength S delivered to a population of N rods. A second response was generated with no flash. The task of the ideal observer was to identify the response containing the flash. The ideal observer was constructed by calculating the optimal linear discriminant separating flash from no-flash trials in each individual rod. The dot product between this template, and each individual response was calculated and then weighted by the likelihood it was generated by the given flash strength (see Materials and methods). This procedure instantiates an optimal linear filtering of the rod output followed by an optimal non-linear weighting that approximates the non-linear threshold performed at the rod-to-rod-bipolar synapse (Field and Rieke, 2002b; Field et al., 2019). The fraction of trials in which the flash response was correctly identified was calculated over 400 trials at each flash strength; flash strengths ranged from 0.001 to 1.0 photo-isomerizations (Rh*) per rod.
For a simulated population of 24 rods, the probability of detecting flashes varied smoothly with flash strength, reaching 75% correct in the 2AFC task between 0.1 and 0.2 Rh*/rod (Figure 7J, blue curve). This performance fell short of an optimal photon detector that is limited only by Poisson variability in the number of photons absorbed from trial-to-trial (Figure 7J, dashed curve) (Field et al., 2019). However, flash detection performance compared well to previous measurements of rod-bipolar cell sensitivity, which place their detection threshold (75% correct in the 2AFC task) between 0.1 and 0.3 Rh*/rod (Figure 7J, vertical gray bar) (Field and Rieke, 2002b; Okawa et al., 2010). Thus, the performance of the simulated rod output with temporally regular multivesicular release falls within the range of flash sensitivities demanded by rod-bipolar cell measurements. This indicates that the observed rod vesicle release dynamics can account for rod-bipolar cell flash sensitivity.
We next examined how this result depended on different choices of model parameters: while the chosen parameters were taken from our measurement of rod vesicle release (biased toward the earliest portion of the recordings and the highest temperatures), our recordings allow for a range of potential values under physiological conditions. We examined how detection performance depended on the multivesicular event rate (Rs), the Erlang factor (Ef), the rebound amplitude (A), the suppression time (T), the decay time (), and the number of rods in the pool (N) (Figure 7J–O). In general, the detection threshold was not strongly impacted over a reasonable range of these model parameters: the flash strength generating 75% correct performance was consistent with the detection threshold of rod-bipolar cells (Figure 7J–O, gray bars).
Two features of this multivesicular release model distinguish it from previous examinations of rod vesicle release. The first feature is multivesicular release. Each multivesicular event contains ~17 vesicles. With a release rate of ~2 multivesicular events/s, this means the rod terminal will release ~34 vesicles/s. Thus, we analyzed how performance of the multivesicular release model compared with asynchronous Poisson release of 34 vesicles/s (Figure 8A; Poisson release model, see Materials and methods). For comparison, we also included performance of an ideal detector and that predicted from the rod photocurrent (Figure 8A; rod photocurrent model, see Materials and methods). The temporally regular multivesicular release model exhibited higher sensitivity than asynchronous Poisson release, given a fixed vesicle budget of 34/s. Thus, the observed statistics of vesicle release from the rod terminal appear to be efficient relative to asynchronous release, given a fixed vesicle budget.
The second feature distinguishing this model is the temporal regularity of vesicle release. The idea of temporally regular release has been postulated previously (Schein and Ahmad, 2005), but not shown experimentally or considered in the context of multivesicular release. The temporal regularity of release is modeled here by the Erlang factor (Ef), which reduces variability in the intervals between multivesicular events. Our data showed a reduction in interevent variability relative to a Poisson process, consistent with Erlang factors of 3–12 (average 6; see Figure 2). Examining detection thresholds as a function of the Erlang factor revealed that over a range of multivesicular release rates, detection thresholds fell considerably as the Erlang factor increased to four and Erlang factors greater than nine exhibited rapidly diminishing returns on detection thresholds (Figure 8B). Thus, further increasing the regularity of interevent intervals (increasing the Erlang factor above 9) yields relatively small improvements to the fidelity of transmitting single-photon responses to rod-bipolar cells.
Discussion
To better understand how responses of rods to single photons can be transmitted reliably to downstream neurons, we examined the mechanisms by which ongoing release from rods can be regulated by small changes in membrane potential. Using IA(glu) as a presynaptic assay for release, we found that at a membrane potential similar to that in darkness, rods released glutamate-filled vesicles in coordinated multivesicular release events at rates that were not well modeled by Poisson statistics. Multivesicular events occurred at more regular intervals consistent with an Erlang distribution in which events occurred after waiting 3–12 Poisson intervals. Multivesicular release events were also extremely sensitive to small voltage changes. Simulating the impact of these properties on release showed that they can improve the ability of rod-bipolar cells to discriminate suppressed release caused by absorption of a single photon from stochastic pauses in release.
The ability of the vertebrate visual system to detect single photons arises from high photon capture rates combined with high gain in the phototransduction cascade of rods. Next, the relatively small changes in membrane voltage produced by absorption of a single photon must produce an appreciable change in the rate of release at the rod synapse. At the synapse from rods to rod-bipolar cell dendrites, there is a non-linear thresholding mechanism in which only responses that exceed a certain size are transmitted, removing baseline noise along with small responses that fall below that threshold (Berntson et al., 2004; Field and Rieke, 2002b; Sampath and Rieke, 2004). In mouse retina, where 20–35 rods converge onto each rod-bipolar cell (Behrens et al., 2016; Calkins and Sterling, 1999; Rao-Mirotznik et al., 1998; Tsukamoto and Omi, 2013), comparisons between rods and rod-bipolar cells showed thresholds that excluded 40–85% of the single-photon events in rods (Berntson et al., 2004; Field and Rieke, 2002b; Schein and Ahmad, 2006). In rabbit retina, where 100 rods may converge onto each rod-bipolar cell, as many as 90% of the single-photon responses are removed by the thresholding mechanism (Trexler et al., 2011). Finally, as many as 100,000 rods eventually converge onto a single ganglion cell, allowing the capture of scarce single photons to produce detectable changes in retinal output (Field et al., 2019; Goodchild et al., 1996; Okawa and Sampath, 2007; Takeshita et al., 2017; Taylor and Smith, 2004).
The modeling presented here (Figure 7) does not strongly constrain the fraction of single-photon responses that fail to traverse the synapse with rod-bipolar cells. This is because these simulations did not incorporate variability in the voltage of the single-photon voltage response or fluctuations in the membrane potential due to continuous noise in phototransduction. Instead, the model began with the estimated dynamics of transmitter release during a simulated light flash (Figure 7B), bypassing these upstream noise sources. To understand how this mode of transmitter release impacts the optimal location of the non-linear threshold at the rod to rod-bipolar synapse would require incorporating noise in phototransduction and determining how this noise (both single-photon response variability and continuous noise) influences the observed temporal variability in multivesicular release events.
The non-linear thresholding mechanism used at rod-bipolar synapses is thought to occur before summation of rod responses in the rod-bipolar cell soma (Field and Rieke, 2002b; Sampath and Rieke, 2004). A key requirement for this non-linear thresholding is that glutamate levels in the cleft must be maintained at levels high enough to keep the rod-bipolar cell glutamate receptors saturated. By relieving saturation in the receptor signaling cascade, larger reductions in glutamate release have disproportionately larger effects (Sampath and Rieke, 2004). In addition to improving regularity, multivesicular release by rods may help ensure that glutamate levels remain high enough at individual synapses to maintain saturation in darkness.
AII amacrine cells show synchronized membrane currents in darkness that arise from synchronous rod inputs into rod-bipolar cells (Grimes et al., 2014). Two rod-bipolar cell dendrites terminate beneath each rod ribbon in mouse retina (Behrens et al., 2016). Together with electrical coupling among rods via cones (Jin et al., 2020), shared responses of two such bipolar cells evoked by large multivesicular release events might contribute to the synchrony among AII amacrine cells.
Our previous recordings of release from a sample of mouse rods suggested an overall rate of release at −40 mV of ~11 v/s at 35°C (Hays et al., 2020b). However, assessments of rundown in the present study suggest that the true rate in vivo prior to patch rupture may be twice this average rate. Earlier computer simulations suggested that, assuming a purely Poisson release process, post-synaptic detection by rod-bipolar cells would require rod release rates of 80–100 v/s (Rao et al., 1994; van Rossum and Smith, 1998). Consistent with this conclusion, our model indicated that Poisson release at a rate of 34 v/s would not provide the level of performance observed in rod-bipolar cells unless we employed an unusually long suppression duration of 500 ms.
A second strategy suggested to overcome noise in release rates is to make release more regular. At hyperpolarized membrane potentials, we found that glutamate release from rods followed a Poisson distribution similar to that seen at many other neurons (Malagon et al., 2016; Miki, 2019). This nearly ubiquitous feature of neurons was established in the 1950s by Bernard Katz who modeled spontaneous release events to develop the quantal hypothesis of neurotransmission (Kavalali, 2015; Kavalali, 2018). The occasional unitary events observed in rods held at −40 mV also showed a Poisson distribution of intervals. However, the dominant feature of release at −40 mV was the appearance of multivesicular release events occurring at fairly regular intervals. The distribution of intervals between multivesicular events could not be explained by Poisson statistics and were instead more consistent with an Erlang process, having a CV narrower than expected from a Poisson process. By waiting a specific number of Poisson intervals between each Erlang event, Erlang processes effectively average the underlying Poisson intervals and thereby reduce variability. It was suggested previously that an Erlang process might play a role in transmitting single-photon responses from rods (Schein and Ahmad, 2005; Schein and Ahmad, 2006). Our experimental data showed a narrower distribution than expected for a Poisson process consistent with the accumulation of as many as 12 Poisson events before each multivesicular release event. The cumulative wait times for replenishment of multiple vesicles to the releasable pool at the base of the synaptic ribbon provides a possible mechanism for achieving an Erlang-like process to enhance regularity in release rates. Predictions of Schein and Ahmad suggested that achieving reliable single-photon detection with a release rate of 100 v/s/ribbon would require an Erlang process with an Erlang factor of 66 (i.e., an accumulation of 66 Poisson intervals) (Schein and Ahmad, 2006). Our model, which among other things differed by including a thresholding non-linearity at the rod-to-rod bipolar cell synapse, showed that even a modest increase in regularity combined with multivesicular release can improve the post-synaptic detection of single-photon events.
Another surprising characteristic of the multivesicular release in rods was its sensitivity to small voltage changes. While the overall rate of release parallels changes in ICa, a small hyperpolarizing voltage change caused appreciable pauses in release. The average interevent interval dropped close to zero during a simulated single-photon response, rather than declining only ~20–25% as predicted from the decline in ICa. Hyperpolarizing the membrane by only 1 mV, similar to the amplitude of single-photon responses in primate rods (Schneeweis and Schnapf, 1995), also caused an appreciable pause in release. The mechanisms that coordinate release among vesicles to produce this exquisite sensitivity remain unclear.
Recent studies on mice in which synaptic ribbons have been eliminated by deleting the key ribbon protein, Ribeye, showed surprisingly modest deficits (Fairless et al., 2020; Okawa et al., 2019) Multivesicular release in both salamander and mouse rods involves ribbon-associated vesicles (Hays et al., 2020a; Hays et al., 2020b). One role of the ribbon may thus be to help organize vesicles to promote multivesicular release, implementing an Erlang process to improve the encoding of light responses near scotopic threshold (Schein and Ahmad, 2005).
The need to rapidly recycle vesicles is a barrier to rods maintaining high release rates. Rod spherules in mouse retina possess 5800–7500 vesicles apiece (Zampighi et al., 2011) and so the entire pool of vesicles would be depleted in roughly a minute if release were maintained at 100 vesicles/s. Our results suggest that the temporally regular multivesicular release strategy employed by rods provides an efficient means of encoding small voltage changes to meet both constraints on high release rates and requirements for reliable signaling at rod ribbon synapses (Sterling and Laughlin, 2015).
Materials and methods
Animals
Control C57/Bl6 mice (RRID:IMSR_JAX:000664) were kept on 12 hr dark–light cycles, and animal handling protocols were approved by the University of Nebraska Medical Center Institutional Animal Care and Use Committee. Mice of both sexes aged 4–12 weeks were euthanized in accordance with AVMA Guidelines for the Euthanasia of Animals by CO2 asphyxiation followed by cervical dislocation (protocol 18-077-06).
Electrophysiology
Request a detailed protocolRecordings from rods were performed using a flatmount preparation of isolated retina. Immediately after euthanizing the mouse, it was enucleated and eyes were placed in Ames’ medium (US Biological; RRID:SCR_013653) bubbled with 95% O2/5% CO2. The cornea of an eye was punctured with a scalpel and the anterior segment removed. The retina was isolated by cutting optic nerve attachments. After making four fine cuts at opposite poles of the retina, it was flattened photoreceptors face up onto a glass slide in the perfusion chamber. The retina was anchored in place using a brain slice harp (Warner Instruments, cat. no. 64–0250). The perfusion chamber was placed on an upright fixed-stage microscope (Nikon E600FN) with a 60× water immersion, long-working distance objective (1.0 NA). Unless otherwise noted, the retina was superfused with room temperature Ames’ solution bubbled with 95% O2/5% CO2 at ~1 mL/min. Prior to recording, outer segments from a region of the retina were removed by gentle suction applied through a broken patch pipette. Rod cell bodies were identified morphologically and targeted with positive pressure using recording electrodes mounted on Huxley–Wall micromanipulators (Sutter Instruments).
Rod recordings were performed in whole-cell voltage-clamp using a Multiclamp 700A amplifier (Axon Instruments/Molecular Devices) and signals digitized with a DigiData 1550 (Axon Instruments/Molecular Devices). Data acquisition and analysis was performed using pClamp 10 Software (Molecular Devices; RRID:SCR_011323). Voltages were not corrected for liquid junction potential (KSCN pipette solution: 3.9 mV). Most experiments were performed in room light, but in some experiments, retinas were prepared and electrodes positioned in darkness under infrared illumination using Gen III night vision goggles (Nitemate NavIII, Litton Industries).
Intracellular pipette solutions for IA(glu) measurements contained (in mM): 120 KSCN, 10 TEA-Cl, 10 HEPES, 1 CaCl2, 1 MgCl2, 0.5 Na-GTP, 5 Mg-ATP, 5 EGTA, 5 phosphocreatine, pH 7.3. Rods were distinguished from cones by their smaller membrane capacitance. Passive membrane properties of rods averaged Cm = 3.2 ± 0.2 pF; Rm = 2.3 ± 0.04 GΩ (mean ± SD, n = 20).
IA(glu) event frequency, kinetics, and charge transfer were measured from events identified with the event finder function in pClamp directed by a template that included a sampling of ~10 manually curated events.
Rod vesicular release and photocurrent models
Erlang release model (model 1)
Request a detailed protocolMultivesicular release from the rod terminal was modeled as an Erlang process (Figure 7). To generate stochastic Erlang events with a spontaneous event rate of Rs and a coefficient of variation , a Poisson process was first generated with an event rate of Rs * Ef. For every Ef Poisson events that were generated, an Erlang event was generated producing stochastic events that were more regularly spaced in time than a Poisson process (see Figure 7B–D and F–H). These discrete events were modeled as delta-functions and then convolved with an alpha waveform:
Finally, independent Gaussian noise was added (Figure 7E,I). To model the release dynamics following a single-photon response, the Poisson process was made time dependent (Figure 7B), such that following the absorption of a simulated photon, the Poisson rate was set to zero for a T seconds, then set to A and exponentially relaxed back to the baseline rate with a time constant of . In these simulations, time was sampled at 0.1 ms to ensure the probability that the Poisson process generated a value other than 0 or 1 was miniscule (<0.0001). In the unlikely event that a value > 1 was sampled, the value was set to 1.
Poisson release model (model 2)
Request a detailed protocolTo compare the flash detection performance of multivesicular release governed by an Erlang process to that of asynchronous Poisson release (Figure 8), a similar procedure was followed as that described in the preceding paragraph. However, Poisson events were not accumulated to generate an Erlang event and the suppression of the vesicular release following a simulated single-photon response was set to 25% of the spontaneous release rate to match previous estimates given the membrane potential of rods and given the voltage dependence of calcium channel activation (Schein and Ahmad, 2005; Schein and Ahmad, 2006).
Rod photocurrent model (model 3)
Request a detailed protocolTo model the rod photocurrent and the detection sensitivity limits set by photocurrent noise sources (Figure 8), we used a generative model of the rod response (Field et al., 2019). A rod response, r(t), was generated from the following equation:
Continuous dark noise, was generated by sampling from a Gaussian distribution and filtered in time to match the power spectrum of measured continuous noise (Field et al., 2019). N photon responses in an individual rod were generated by the first term in Equation 2 by sampling from a Poisson distribution with a mean given by the flash strength on a given trial. The mean single-photon response is given by , and the covariance of the single-photon response is captured by summing over a weighted set, wi, of eigenvectors, , derived from the covariance matrix of the single-photon response (Field et al., 2019). Finally, discrete noise events caused by the thermal activation of rhodopsin were captured by the second term of Equation 2: Nd thermal isomerizations were generated using the same formulation as for the single-photon response, but Nd was determined by sampling from a Poisson distribution with a mean given by the thermal isomerization rate at 35°C and each isomerization event occurred at random and independent times given by td.
In the original formulation of this model (Field et al., 2019), the noise parameters were derived from measurements of the photocurrent from primate rods. We adjusted these parameters to account for the greater amounts of continuous noise and variability in the single-photon response present in mouse rods (both relative to the amplitude of the mean single-photon response). Specifically, the continuous noise was increased by 22% and single-photon response variability was increased by 37.5% (Field and Rieke, 2002a; Field and Rieke, 2002b). The kinetics of the single-photon response and the shape of the power spectrum of the continuous noise are both similar between mouse and primate rods, so no adjustments were made to these quantities. The thermal rate in mammalian and mouse rods is somewhat uncertain, with the literature allowing for a wide range of values between 0.001 and 0.015 Rh*/rod/s, so we chose an intermediate value of 0.005 Rh*/rod/s (Burns et al., 2002; Field et al., 2019; Fu et al., 2008; Yue et al., 2017). Importantly, this model for simulating rod responses reproduces rod detection and temporal sensitivities in a 2AFC task (Field et al., 2019), supporting the application of this model for analyzing task performance of pools of rods.
2AFC task and ideal observer analysis
Request a detailed protocolUsing each of the three models described in the previous section, detection threshold was determined using a 2AFC task and an ideal observer analysis. For each pair of trials, the task was to determine which contained the ‘flash’ and which did not. The ideal observer was constructed from a training set of simulated flash and no-flash (null) responses (1000 of each). Note, in the training phase, the simulated flash responses consisted only of single-photon responses without Poisson fluctuations in the number of photons absorbed. From these simulated responses, a difference-of-means linear discriminant was constructed (Duda et al., 2001):
where D is the linear discriminant, and are the mean of the flash and no-flash responses, respectively. Then the dot products of the flash and no-flash responses with the discriminant were calculated; each set of projections yielded an approximately Gaussian distribution associated with the flash and no-flash trials. Thus, their means and standard deviation were used to summarize the distributions: µA, µB, σA, and σB, where A and B denote the distributions of flash and no-flash responses, respectively, after computing their dot products with the discriminant. Then a new ‘test’ set of flash and no-flash responses were simulated and discriminated by the ideal observer using a Bayes-optimal non-linear weighting.
where P is the number of rods in the pool, r(t) is a given response (flash or no flash), D is the discriminant from Equation 3, and wi is the non-linear weight given by:
G(X|μ,σ) is the probability of sampling X from a Gaussian distribution with mean, μ, and standard deviation σ; P(Y|λ) is the probability of sampling Y from a Poisson distribution with mean λ; and f in Equation 5 is the flash strength used to generate the ‘flash’ responses, which acts like a prior in the Bayesian sense. Note, in the ‘test’ set, the simulated flash responses included Poisson fluctuations in the number of photons absorbed from trial-to-trial.
Responses were simulated for flash strength f and compared to responses of flash strength 0. If Rf (flash response) > R0 (no-flash response), the trial was scored as correct; if Rf < R0, the trial was scored as incorrect; and if Rf = R0, the trial was randomly scored as correct with 50% probability.
Statistical analysis
Request a detailed protocolStatistical analysis was performed on GraphPad Prism versions 7–9 (RRID:SCR_002798), and all data are represented as mean ± SD unless otherwise noted. Recordings were performed using rods from at least three different retinas.
To test whether release rates reflected a Poisson process, interevent interval histograms were fit with a one-phase exponential decay. In a second approach, the number of release events was tallied over a long period (30 s) to calculate the average number of release events per second (λ). A fractional frequency distribution of events was then made using 1 s bins and fit with a Poisson model:
where P is the probability of observing x number of events and e is Euler’s constant.
Data availability
All data generated or analyzed during this study are included in the manuscript and supporting file. Source data for Figure 5B are provided.
References
-
Retinal noise and absolute thresholdJournal of the Optical Society of America 46:634–639.https://doi.org/10.1364/JOSA.46.000634
-
The photocurrent, noise and spectral sensitivity of rods of the monkey Macaca fascicularisThe Journal of Physiology 357:575–607.https://doi.org/10.1113/jphysiol.1984.sp015518
-
The photovoltage of rods and cones in the dark-adapted mouse retinaThe Journal of Physiology 590:3841–3855.https://doi.org/10.1113/jphysiol.2011.226878
-
Mechanisms, pools, and sites of spontaneous vesicle release at synapses of rod and cone photoreceptorsEuropean Journal of Neuroscience 44:2015–2027.https://doi.org/10.1111/ejn.13288
-
ERG responses in mice with deletion of the synaptic ribbon component RIBEYEInvestigative Opthalmology & Visual Science 61:37.https://doi.org/10.1167/iovs.61.5.37
-
Temporal resolution of single-photon responses in primate rod photoreceptors and limits imposed by cellular noiseJournal of Neurophysiology 121:255–268.https://doi.org/10.1152/jn.00683.2018
-
Quantal noise from human red cone pigmentNature Neuroscience 11:565–571.https://doi.org/10.1038/nn.2110
-
cGMP in mouse rods: the spatiotemporal dynamics underlying single photon responsesFrontiers in Molecular Neuroscience 8:6.https://doi.org/10.3389/fnmol.2015.00006
-
Resting and stimulated mouse rod photoreceptors show distinct patterns of vesicle release at ribbon synapsesJournal of General Physiology 152:e202012716.https://doi.org/10.1085/jgp.202012716
-
Energy, quanta, and visionJournal of General Physiology 25:819–840.https://doi.org/10.1085/jgp.25.6.819
-
Gap-junctional coupling and absolute sensitivity of photoreceptors in macaque retinaJournal of Neuroscience 25:11201–11209.https://doi.org/10.1523/JNEUROSCI.3416-05.2005
-
The mechanisms and functions of spontaneous neurotransmitter releaseNature Reviews Neuroscience 16:5–16.https://doi.org/10.1038/nrn3875
-
Spontaneous neurotransmission: a form of neural communication comes of ageJournal of Neuroscience Research 96:331–334.https://doi.org/10.1002/jnr.24207
-
Neuronal Ca2+ signalling at rest and during spontaneous neurotransmissionThe Journal of Physiology 598:1649–1654.https://doi.org/10.1113/JP276541
-
Counting vesicular release events reveals binomial release statistics at single glutamatergic synapsesThe Journal of Neuroscience 36:4010–4025.https://doi.org/10.1523/JNEUROSCI.4352-15.2016
-
What we can learn from cumulative numbers of vesicular release eventsFrontiers in Cellular Neuroscience 13:257.https://doi.org/10.3389/fncel.2019.00257
-
Optimal processing of photoreceptor signals is required to maximize behavioural sensitivityThe Journal of Physiology 588:1947–1960.https://doi.org/10.1113/jphysiol.2010.188573
-
Anion currents and predicted glutamate flux through a neuronal glutamate transporterThe Journal of Neuroscience 18:7099–7110.https://doi.org/10.1523/JNEUROSCI.18-18-07099.1998
-
Rate of quantal transmitter release at the mammalian rod synapseBiophysical Journal 67:57–63.https://doi.org/10.1016/S0006-3495(94)80454-0
-
Transmitter concentration at a Three-Dimensional synapseJournal of Neurophysiology 80:3163–3172.https://doi.org/10.1152/jn.1998.80.6.3163
-
Origin of reproducibility in the responses of retinal rods to single photonsBiophysical Journal 75:1836–1857.https://doi.org/10.1016/S0006-3495(98)77625-8
-
Counting every quantumThe Journal of Physiology 223:131–150.https://doi.org/10.1113/jphysiol.1972.sp009838
-
A clockwork hypothesis: synaptic release by rod photoreceptors must be regularBiophysical Journal 89:3931–3949.https://doi.org/10.1529/biophysj.105.070623
-
Functional properties of the retinal glutamate transporters GLT-1c and EAAT5Journal of Biological Chemistry 289:1815–1824.https://doi.org/10.1074/jbc.M113.517177
-
Synaptic Ca2+ in darkness is lower in rods than cones, causing slower tonic release of vesiclesJournal of Neuroscience 27:5033–5042.https://doi.org/10.1523/JNEUROSCI.5386-06.2007
-
Ideal observer analysis of signal quality in retinal circuitsProgress in Retinal and Eye Research 28:263–288.https://doi.org/10.1016/j.preteyeres.2009.05.001
-
Processing of single-photon responses in the mammalian on and off retinal pathways at the sensitivity limit of visionPhilosophical Transactions of the Royal Society B: Biological Sciences 372:20160073.https://doi.org/10.1098/rstb.2016.0073
-
Direct detection of a single photon by humansNature Communications 7:12172.https://doi.org/10.1038/ncomms12172
-
Nonlinearity and noise at the rod-rod bipolar cell synapseVisual Neuroscience 28:61–68.https://doi.org/10.1017/S0952523810000301
-
Functional allocation of synaptic contacts in microcircuits from rods via rod bipolar to AII amacrine cells in the mouse retinaJournal of Comparative Neurology 521:3541–3555.https://doi.org/10.1002/cne.23370
-
Noise removal at the rod synapse of mammalian retinaVisual Neuroscience 15:809–821.https://doi.org/10.1017/S0952523898155037
Article and author information
Author details
Funding
National Eye Institute (EY10542)
- Wallace B Thoreson
National Eye Institute (EY32396)
- Wallace B Thoreson
UNMC College of Medicine (Graduate Fellowship)
- Cassandra L Hays
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
Funding was provided by UNMC Graduate Fellowship (CLH) and NIH grants EY10542 and EY32396 (WBT).
Ethics
Animal experimentation: Animal handling protocols were approved by the University of Nebraska Medical Center Institutional Animal Care and Use Committee. Mice of both sexes aged 4-12 weeks were euthanized in accordance with AVMA Guidelines for the Euthanasia of Animals by CO2 asphyxiation followed by cervical dislocation (protocol 18-077-06).
Copyright
© 2021, Hays 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.
Metrics
-
- 958
- views
-
- 157
- downloads
-
- 13
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
In congenital stationary night blindness, type 2 (CSNB2)—a disorder involving the Cav1.4 (L-type) Ca2+ channel—visual impairment is mild considering that Cav1.4 mediates synaptic release from rod and cone photoreceptors. Here, we addressed this conundrum using a Cav1.4 knockout (KO) mouse and a knock-in (G369i KI) mouse expressing a non-conducting Cav1.4. Surprisingly, Cav3 (T-type) Ca2+ currents were detected in cones of G369i KI mice and Cav1.4 KO mice but not in cones of wild-type mouse, ground squirrels, and macaque retina. Whereas Cav1.4 KO mice are blind, G369i KI mice exhibit normal photopic (i.e. cone-mediated) visual behavior. Cone synapses, which fail to form in Cav1.4 KO mice, are present, albeit enlarged, and with some errors in postsynaptic wiring in G369i KI mice. While Cav1.4 KO mice lack evidence of cone synaptic responses, electrophysiological recordings in G369i KI mice revealed nominal transmission from cones to horizontal cells and bipolar cells. In CSNB2, we propose that Cav3 channels maintain cone synaptic output provided that the nonconducting role of Cav1.4 in cone synaptogenesis remains intact. Our findings reveal an unexpected form of homeostatic plasticity that relies on a non-canonical role of an ion channel.
-
- Developmental Biology
- Neuroscience
Otolith organs in the inner ear and neuromasts in the fish lateral-line harbor two populations of hair cells oriented to detect stimuli in opposing directions. The underlying mechanism is highly conserved: the transcription factor EMX2 is regionally expressed in just one hair cell population and acts through the receptor GPR156 to reverse cell orientation relative to the other population. In mouse and zebrafish, loss of Emx2 results in sensory organs that harbor only one hair cell orientation and are not innervated properly. In zebrafish, Emx2 also confers hair cells with reduced mechanosensory properties. Here, we leverage mouse and zebrafish models lacking GPR156 to determine how detecting stimuli of opposing directions serves vestibular function, and whether GPR156 has other roles besides orienting hair cells. We find that otolith organs in Gpr156 mouse mutants have normal zonal organization and normal type I-II hair cell distribution and mechano-electrical transduction properties. In contrast, gpr156 zebrafish mutants lack the smaller mechanically evoked signals that characterize Emx2-positive hair cells. Loss of GPR156 does not affect orientation-selectivity of afferents in mouse utricle or zebrafish neuromasts. Consistent with normal otolith organ anatomy and afferent selectivity, Gpr156 mutant mice do not show overt vestibular dysfunction. Instead, performance on two tests that engage otolith organs is significantly altered – swimming and off-vertical-axis rotation. We conclude that GPR156 relays hair cell orientation and transduction information downstream of EMX2, but not selectivity for direction-specific afferents. These results clarify how molecular mechanisms that confer bi-directionality to sensory organs contribute to function, from single hair cell physiology to animal behavior.