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 (I_A(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.

To study voltage-dependent vesicular glutamate release from rods, we recorded I_A(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 I_Ca 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 Ca²⁺ influx and derived from a pool of ribbon-associated vesicles shared with evoked release.

Figure 1

Download asset Open asset

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 Ca²⁺ independent fusion events or from chance openings of voltage-gated Ca²⁺ channels (Cork et al., 2016; Kavalali, 2015; Kavalali, 2020). The persistence of spontaneous events in rods held at −70 mV, even after blocking Ca²⁺ channels with extracellular Cd²⁺ or strongly buffering intracellular Ca²⁺ with 10 mM BAPTA, showed evidence for Ca²⁺-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 (R² = 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 (R² = 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).

Figure 2

Download asset Open asset

These results indicate that the rate of Ca²⁺-independent spontaneous release from rods at hyperpolarized membrane potentials is a Poisson process. Does Poisson release persist when Ca²⁺-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, R² = 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 I_A(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 R² 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 R² = 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:

${\displaystyle \mathrm{C}\mathrm{V}=1/(\mathrm{k}{)}^{1/2}}$

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.

Figure 3

Download asset Open asset

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.

Figure 4

Download asset Open asset

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).

Figure 5

Download asset Open asset

Figure 5—source data 1

Start times for multivesicular events during simulated single-photon light responses.

https://cdn.elifesciences.org/articles/67446/elife-67446-fig5-data1-v2.xlsx

Download elife-67446-fig5-data1-v2.xlsx

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.

Figure 6

Download asset Open asset

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 $\tau$ (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, E_f. When the counter reached E_f, 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(E_f). 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).

Figure 7

Download asset Open asset

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 (R_s), the Erlang factor (E_f), the rebound amplitude (A), the suppression time (T), the decay time ($\tau$), 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.

Figure 8

Download asset Open asset

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 (E_f), 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.

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 I_A(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 I_Ca, 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 I_Ca. 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).

All data generated or analyzed during this study are included in the manuscript and supporting file. Source data for Figure 5B are provided.

1. 1. Arriza JL
   2. Eliasof S
   3. Kavanaugh MP
   4. Amara SG

   (1997) Excitatory amino acid transporter 5, a retinal glutamate transporter coupled to a chloride conductance

   PNAS 94:4155–4160.

   https://doi.org/10.1073/pnas.94.8.4155

   • PubMed
   • Google Scholar
2. 1. Barlow HB

   (1956) Retinal noise and absolute threshold

   Journal of the Optical Society of America 46:634–639.

   https://doi.org/10.1364/JOSA.46.000634

   • PubMed
   • Google Scholar
3. 1. Baylor DA
   2. Nunn BJ
   3. Schnapf JL

   (1984) The photocurrent, noise and spectral sensitivity of rods of the monkey Macaca fascicularis

   The Journal of Physiology 357:575–607.

   https://doi.org/10.1113/jphysiol.1984.sp015518

   • PubMed
   • Google Scholar
4. 1. Behrens C
   2. Schubert T
   3. Haverkamp S
   4. Euler T
   5. Berens P

   (2016) Connectivity map of bipolar cells and photoreceptors in the mouse retina

   eLife 5:e20041.

   https://doi.org/10.7554/eLife.20041

   • PubMed
   • Google Scholar
5. 1. Berntson A
   2. Smith RG
   3. Taylor WR

   (2004) Transmission of single photon signals through a binary synapse in the mammalian retina

   Visual Neuroscience 21:693–702.

   https://doi.org/10.1017/S0952523804215048

   • PubMed
   • Google Scholar
6. 1. Burns ME
   2. Mendez A
   3. Chen J
   4. Baylor DA

   (2002) Dynamics of cyclic GMP synthesis in retinal rods

   Neuron 36:81–91.

   https://doi.org/10.1016/S0896-6273(02)00911-X

   • PubMed
   • Google Scholar
7. 1. Calkins DJ
   2. Sterling P

   (1999) Evidence that circuits for spatial and color vision segregate at the first retinal synapse

   Neuron 24:313–321.

   https://doi.org/10.1016/S0896-6273(00)80846-6

   • PubMed
   • Google Scholar
8. 1. Cangiano L
   2. Asteriti S
   3. Cervetto L
   4. Gargini C

   (2012) The photovoltage of rods and cones in the dark-adapted mouse retina

   The Journal of Physiology 590:3841–3855.

   https://doi.org/10.1113/jphysiol.2011.226878

   • PubMed
   • Google Scholar
9. 1. Chichilnisky EJ
   2. Rieke F

   (2005) Detection sensitivity and temporal resolution of visual signals near absolute threshold in the salamander retina

   Journal of Neuroscience 25:318–330.

   https://doi.org/10.1523/JNEUROSCI.2339-04.2005

   • PubMed
   • Google Scholar
10. 1. Cork KM
    2. Van Hook MJ
    3. Thoreson WB

    (2016) Mechanisms, pools, and sites of spontaneous vesicle release at synapses of rod and cone photoreceptors

    European Journal of Neuroscience 44:2015–2027.

    https://doi.org/10.1111/ejn.13288

    • PubMed
    • Google Scholar
11. 1. Dembla E
    2. Dembla M
    3. Maxeiner S
    4. Schmitz F

    (2020) Synaptic ribbons foster active zone stability and illumination-dependent active zone enrichment of RIM2 and Cav1.4 in photoreceptor synapses

    Scientific Reports 10:5957.

    https://doi.org/10.1038/s41598-020-62734-0

    • PubMed
    • Google Scholar
12. Book

    1. Duda R
    2. Hart P
    3. Stork D

    (2001)

    Pattern Classification

    New York: Wiley & Sons.

    • Google Scholar
13. 1. Fairless R
    2. Williams SK
    3. Katiyar R
    4. Maxeiner S
    5. Schmitz F
    6. Diem R

    (2020) ERG responses in mice with deletion of the synaptic ribbon component RIBEYE

    Investigative Opthalmology & Visual Science 61:37.

    https://doi.org/10.1167/iovs.61.5.37

    • Google Scholar
14. 1. Field GD
    2. Uzzell V
    3. Chichilnisky EJ
    4. Rieke F

    (2019) Temporal resolution of single-photon responses in primate rod photoreceptors and limits imposed by cellular noise

    Journal of Neurophysiology 121:255–268.

    https://doi.org/10.1152/jn.00683.2018

    • PubMed
    • Google Scholar
15. 1. Field GD
    2. Rieke F

    (2002a) Mechanisms regulating variability of the single photon responses of mammalian rod photoreceptors

    Neuron 35:733–747.

    https://doi.org/10.1016/S0896-6273(02)00822-X

    • PubMed
    • Google Scholar
16. 1. Field GD
    2. Rieke F

    (2002b) Nonlinear signal transfer from mouse rods to bipolar cells and implications for visual sensitivity

    Neuron 34:773–785.

    https://doi.org/10.1016/S0896-6273(02)00700-6

    • PubMed
    • Google Scholar
17. 1. Fu Y
    2. Kefalov V
    3. Luo DG
    4. Xue T
    5. Yau KW

    (2008) Quantal noise from human red cone pigment

    Nature Neuroscience 11:565–571.

    https://doi.org/10.1038/nn.2110

    • PubMed
    • Google Scholar
18. 1. Goodchild AK
    2. Ghosh KK
    3. Martin PR

    (1996) Comparison of photoreceptor spatial density and ganglion cell morphology in the retina of human, Macaque Monkey, cat, and the marmoset callithrix jacchus

    The Journal of Comparative Neurology 366:55–75.

    https://doi.org/10.1002/(SICI)1096-9861(19960226)366:1<55::AID-CNE5>3.0.CO;2-J

    • PubMed
    • Google Scholar
19. 1. Grassmeyer JJ
    2. Cahill AL
    3. Hays CL
    4. Barta C
    5. Quadros RM
    6. Gurumurthy CB
    7. Thoreson WB

    (2019) Ca²⁺ sensor synaptotagmin-1 mediates exocytosis in mammalian photoreceptors

    eLife 8:e45946.

    https://doi.org/10.7554/eLife.45946

    • PubMed
    • Google Scholar
20. 1. Grimes WN
    2. Hoon M
    3. Briggman KL
    4. Wong RO
    5. Rieke F

    (2014) Cross-synaptic synchrony and transmission of signal and noise across the mouse retina

    eLife 3:e03892.

    https://doi.org/10.7554/eLife.03892

    • PubMed
    • Google Scholar
21. 1. Gross OP
    2. Pugh EN
    3. Burns ME

    (2015) cGMP in mouse rods: the spatiotemporal dynamics underlying single photon responses

    Frontiers in Molecular Neuroscience 8:6.

    https://doi.org/10.3389/fnmol.2015.00006

    • PubMed
    • Google Scholar
22. 1. Hasegawa J
    2. Obara T
    3. Tanaka K
    4. Tachibana M

    (2006) High-density presynaptic transporters are required for glutamate removal from the first visual synapse

    Neuron 50:63–74.

    https://doi.org/10.1016/j.neuron.2006.02.022

    • PubMed
    • Google Scholar
23. 1. Hays CL
    2. Grassmeyer JJ
    3. Wen X
    4. Janz R
    5. Heidelberger R
    6. Thoreson WB

    (2020a) Simultaneous release of multiple vesicles from rods involves synaptic ribbons and syntaxin 3B

    Biophysical Journal 118:967–979.

    https://doi.org/10.1016/j.bpj.2019.10.006

    • PubMed
    • Google Scholar
24. 1. Hays CL
    2. Sladek AL
    3. Thoreson WB

    (2020b) Resting and stimulated mouse rod photoreceptors show distinct patterns of vesicle release at ribbon synapses

    Journal of General Physiology 152:e202012716.

    https://doi.org/10.1085/jgp.202012716

    • Google Scholar
25. 1. Hecht S
    2. Shlaer S
    3. Pirenne MH

    (1942) Energy, quanta, and vision

    Journal of General Physiology 25:819–840.

    https://doi.org/10.1085/jgp.25.6.819

    • Google Scholar
26. 1. Hornstein EP
    2. Verweij J
    3. Li PH
    4. Schnapf JL

    (2005) Gap-junctional coupling and absolute sensitivity of photoreceptors in macaque retina

    Journal of Neuroscience 25:11201–11209.

    https://doi.org/10.1523/JNEUROSCI.3416-05.2005

    • PubMed
    • Google Scholar
27. 1. Jin N
    2. Zhang Z
    3. Keung J
    4. Youn SB
    5. Ishibashi M
    6. Tian LM
    7. Marshak DW
    8. Solessio E
    9. Umino Y
    10. Fahrenfort I
    11. Kiyama T
    12. Mao CA
    13. You Y
    14. Wei H
    15. Wu J
    16. Postma F
    17. Paul DL
    18. Massey SC
    19. Ribelayga CP

    (2020) Molecular and functional architecture of the mouse photoreceptor network

    Science Advances 6:eaba7232.

    https://doi.org/10.1126/sciadv.aba7232

    • PubMed
    • Google Scholar
28. 1. Kavalali ET

    (2015) The mechanisms and functions of spontaneous neurotransmitter release

    Nature Reviews Neuroscience 16:5–16.

    https://doi.org/10.1038/nrn3875

    • PubMed
    • Google Scholar
29. 1. Kavalali ET

    (2018) Spontaneous neurotransmission: a form of neural communication comes of age

    Journal of Neuroscience Research 96:331–334.

    https://doi.org/10.1002/jnr.24207

    • PubMed
    • Google Scholar
30. 1. Kavalali ET

    (2020) Neuronal Ca²⁺ signalling at rest and during spontaneous neurotransmission

    The Journal of Physiology 598:1649–1654.

    https://doi.org/10.1113/JP276541

    • PubMed
    • Google Scholar
31. 1. Malagon G
    2. Miki T
    3. Llano I
    4. Neher E
    5. Marty A

    (2016) Counting vesicular release events reveals binomial release statistics at single glutamatergic synapses

    The Journal of Neuroscience 36:4010–4025.

    https://doi.org/10.1523/JNEUROSCI.4352-15.2016

    • PubMed
    • Google Scholar
32. 1. Miki T

    (2019) What we can learn from cumulative numbers of vesicular release events

    Frontiers in Cellular Neuroscience 13:257.

    https://doi.org/10.3389/fncel.2019.00257

    • PubMed
    • Google Scholar
33. 1. Okawa H
    2. Miyagishima KJ
    3. Arman AC
    4. Hurley JB
    5. Field GD
    6. Sampath AP

    (2010) Optimal processing of photoreceptor signals is required to maximize behavioural sensitivity

    The Journal of Physiology 588:1947–1960.

    https://doi.org/10.1113/jphysiol.2010.188573

    • PubMed
    • Google Scholar
34. 1. Okawa H
    2. Yu WQ
    3. Matti U
    4. Schwarz K
    5. Odermatt B
    6. Zhong H
    7. Tsukamoto Y
    8. Lagnado L
    9. Rieke F
    10. Schmitz F
    11. Wong ROL

    (2019) Dynamic assembly of ribbon synapses and circuit maintenance in a vertebrate sensory system

    Nature Communications 10:2167.

    https://doi.org/10.1038/s41467-019-10123-1

    • PubMed
    • Google Scholar
35. 1. Okawa H
    2. Sampath AP

    (2007) Optimization of single-photon response transmission at the rod-to-rod bipolar synapse

    Physiology 22:279–286.

    https://doi.org/10.1152/physiol.00007.2007

    • PubMed
    • Google Scholar
36. 1. Otis TS
    2. Jahr CE

    (1998) Anion currents and predicted glutamate flux through a neuronal glutamate transporter

    The Journal of Neuroscience 18:7099–7110.

    https://doi.org/10.1523/JNEUROSCI.18-18-07099.1998

    • PubMed
    • Google Scholar
37. 1. Rao R
    2. Buchsbaum G
    3. Sterling P

    (1994) Rate of quantal transmitter release at the mammalian rod synapse

    Biophysical Journal 67:57–63.

    https://doi.org/10.1016/S0006-3495(94)80454-0

    • PubMed
    • Google Scholar
38. 1. Rao-Mirotznik R
    2. Buchsbaum G
    3. Sterling P

    (1998) Transmitter concentration at a Three-Dimensional synapse

    Journal of Neurophysiology 80:3163–3172.

    https://doi.org/10.1152/jn.1998.80.6.3163

    • Google Scholar
39. 1. Reingruber J
    2. Holcman D
    3. Fain GL

    (2015) How rods respond to single photons: key adaptations of a G-protein cascade that enable vision at the physical limit of perception

    BioEssays 37:1243–1252.

    https://doi.org/10.1002/bies.201500081

    • PubMed
    • Google Scholar
40. 1. Rieke F
    2. Baylor DA

    (1998) Origin of reproducibility in the responses of retinal rods to single photons

    Biophysical Journal 75:1836–1857.

    https://doi.org/10.1016/S0006-3495(98)77625-8

    • PubMed
    • Google Scholar
41. 1. Sakitt B

    (1972) Counting every quantum

    The Journal of Physiology 223:131–150.

    https://doi.org/10.1113/jphysiol.1972.sp009838

    • PubMed
    • Google Scholar
42. 1. Sampath AP
    2. Rieke F

    (2004) Selective transmission of single photon responses by saturation at the rod-to-rod bipolar synapse

    Neuron 41:431–443.

    https://doi.org/10.1016/S0896-6273(04)00005-4

    • PubMed
    • Google Scholar
43. 1. Schein S
    2. Ahmad KM

    (2005) A clockwork hypothesis: synaptic release by rod photoreceptors must be regular

    Biophysical Journal 89:3931–3949.

    https://doi.org/10.1529/biophysj.105.070623

    • PubMed
    • Google Scholar
44. 1. Schein S
    2. Ahmad KM

    (2006) Efficiency of synaptic transmission of single-photon events from rod photoreceptor to rod bipolar dendrite

    Biophysical Journal 91:3257–3267.

    https://doi.org/10.1529/biophysj.106.091744

    • PubMed
    • Google Scholar
45. 1. Schneeweis DM
    2. Schnapf JL

    (1995) Photovoltage of rods and cones in the macaque retina

    Science 268:1053–1056.

    https://doi.org/10.1126/science.7754386

    • PubMed
    • Google Scholar
46. 1. Schneider N
    2. Cordeiro S
    3. Machtens J-P
    4. Braams S
    5. Rauen T
    6. Fahlke C

    (2014) Functional properties of the retinal glutamate transporters GLT-1c and EAAT5

    Journal of Biological Chemistry 289:1815–1824.

    https://doi.org/10.1074/jbc.M113.517177

    • Google Scholar
47. 1. Sheng Z
    2. Choi SY
    3. Dharia A
    4. Li J
    5. Sterling P
    6. Kramer RH

    (2007) Synaptic Ca2+ in darkness is lower in rods than cones, causing slower tonic release of vesicles

    Journal of Neuroscience 27:5033–5042.

    https://doi.org/10.1523/JNEUROSCI.5386-06.2007

    • PubMed
    • Google Scholar
48. 1. Smith RG
    2. Dhingra NK

    (2009) Ideal observer analysis of signal quality in retinal circuits

    Progress in Retinal and Eye Research 28:263–288.

    https://doi.org/10.1016/j.preteyeres.2009.05.001

    • PubMed
    • Google Scholar
49. Book

    1. Sterling P
    2. Laughlin S

    (2015)

    Principles of Neural Design

    Cambridge MIT Press.

    • Google Scholar
50. 1. Takeshita D
    2. Smeds L
    3. Ala-Laurila P

    (2017) Processing of single-photon responses in the mammalian on and off retinal pathways at the sensitivity limit of vision

    Philosophical Transactions of the Royal Society B: Biological Sciences 372:20160073.

    https://doi.org/10.1098/rstb.2016.0073

    • Google Scholar
51. 1. Taylor WR
    2. Smith RG

    (2004) Transmission of scotopic signals from the rod to rod-bipolar cell in the mammalian retina

    Vision Research 44:3269–3276.

    https://doi.org/10.1016/j.visres.2004.07.043

    • PubMed
    • Google Scholar
52. 1. Tinsley JN
    2. Molodtsov MI
    3. Prevedel R
    4. Wartmann D
    5. Espigulé-Pons J
    6. Lauwers M
    7. Vaziri A

    (2016) Direct detection of a single photon by humans

    Nature Communications 7:12172.

    https://doi.org/10.1038/ncomms12172

    • PubMed
    • Google Scholar
53. 1. Trexler EB
    2. Casti AR
    3. Zhang Y

    (2011) Nonlinearity and noise at the rod-rod bipolar cell synapse

    Visual Neuroscience 28:61–68.

    https://doi.org/10.1017/S0952523810000301

    • PubMed
    • Google Scholar
54. 1. Tsukamoto Y
    2. Omi N

    (2013) Functional allocation of synaptic contacts in microcircuits from rods via rod bipolar to AII amacrine cells in the mouse retina

    Journal of Comparative Neurology 521:3541–3555.

    https://doi.org/10.1002/cne.23370

    • Google Scholar
55. 1. van Rossum MC
    2. Smith RG

    (1998) Noise removal at the rod synapse of mammalian retina

    Visual Neuroscience 15:809–821.

    https://doi.org/10.1017/S0952523898155037

    • PubMed
    • Google Scholar
56. 1. Yue WW
    2. Frederiksen R
    3. Ren X
    4. Luo DG
    5. Yamashita T
    6. Shichida Y
    7. Cornwall MC
    8. Yau KW

    (2017) Spontaneous activation of visual pigments in relation to openness/closedness of chromophore-binding pocket

    eLife 6:e18492.

    https://doi.org/10.7554/eLife.18492

    • PubMed
    • Google Scholar
57. 1. Zampighi GA
    2. Schietroma C
    3. Zampighi LM
    4. Woodruff M
    5. Wright EM
    6. Brecha NC

    (2011) Conical tomography of a ribbon synapse: structural evidence for vesicle fusion

    PLOS ONE 6:e16944.

    https://doi.org/10.1371/journal.pone.0016944

    • PubMed
    • Google Scholar
58. 1. Zhang C
    2. Peskin CS

    (2015) Improved signaling as a result of randomness in synaptic vesicle release

    PNAS 112:14954–14959.

    https://doi.org/10.1073/pnas.1513160112

    • PubMed
    • Google Scholar

Author details

1. Cassandra L Hays

   1. Truhlsen Eye Institute and Department of Ophthalmology and Visual Sciences, College of Medicine, University of Nebraska Medical Center, Omaha, United States
   2. Cellular and Integrative Physiology, Omaha, United States

   Contribution

   Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing - original draft, Writing - review and editing

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0003-4481-3495

2. Asia L Sladek

   Pharmacology and Experimental Neuroscience, Omaha, United States

   Contribution

   Investigation, Writing - review and editing

   Competing interests

   No competing interests declared

3. Greg D Field

   Department of Neurobiology, Duke University School of Medicine, Durham, United States

   Contribution

   Conceptualization, Software, Formal analysis, Investigation, Methodology, Writing - original draft, Writing - review and editing

   For correspondence

   field@neuro.duke.edu

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0001-5942-2679

4. Wallace B Thoreson

   1. Truhlsen Eye Institute and Department of Ophthalmology and Visual Sciences, College of Medicine, University of Nebraska Medical Center, Omaha, United States
   2. Pharmacology and Experimental Neuroscience, Omaha, United States

   Contribution

   Conceptualization, Data curation, Formal analysis, Supervision, Funding acquisition, Investigation, Methodology, Writing - original draft, Project administration, Writing - review and editing

   For correspondence

   wbthores@unmc.edu

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0001-7104-042X

• 897

  Page views

• 148

  Downloads

• 9

  Citations

Article citation count generated by polling the highest count across the following sources: PubMed Central, Crossref, Scopus.