Spontaneous activity of auditory nerve fibers in the barn owl (Tyto alba): analyses of interspike interval distributions

Variability and Randomness of the Instantaneous Firing Rate

The apparent stochastic nature of neuronal activity significantly affects the reliability of neuronal coding. To quantify the encountered fluctuations, both in neural data and simulations, the notions of variability and randomness of inter-spike intervals have been proposed and studied. In this article we focus on the concept of the instantaneous firing rate, which is also based on the spike timing. We use several classical statistical models of neuronal activity and we study the corresponding probability distributions of the instantaneous firing rate. To characterize the firing rate variability and randomness under different spiking regimes, we use different indices of statistical dispersion. We find that the relationship between the variability of interspike intervals and the instantaneous firing rate is not straightforward in general. Counter-intuitively, an increase in the randomness (based on entropy) of spike times may either decrease or increase the randomness of instantaneous firing rate, in dependence on the neuronal firing model. Finally, we apply our methods to experimental data, establishing that instantaneous rate analysis can indeed provide additional information about the spiking activity.

A numerical method for computing interval distributions for an inhomogeneous Poisson point process modified by random dead times

The inhomogeneous Poisson point process is a common model for time series of discrete, stochastic events. When an event from a point process is detected, it may trigger a random dead time in the detector, during which subsequent events will fail to be detected. It can be difficult or impossible to obtain a closed-form expression for the distribution of intervals between detections, even when the rate function (often referred to as the intensity function) and the dead-time distribution are given. Here, a method is presented to numerically compute the interval distribution expected for any arbitrary inhomogeneous Poisson point process modified by dead times drawn from any arbitrary distribution. In neuroscience, such a point process is used to model trains of neuronal spikes triggered by the detection of excitatory events while the neuron is not refractory. The assumptions of the method are that the process is observed over a finite observation window and that the detector is not in a dead state at the start of the observation window. Simulations are used to verify the method for several example point processes. The method should be useful for modeling and understanding the relationships between the rate functions and interval distributions of the event and detection processes, and how these relationships depend on the dead-time distribution.

Mathematical Modeling and Analyses of Interspike-Intervals of Spontaneous Activity in Afferent Neurons of the Zebrafish Lateral Line

Without stimuli, hair cells spontaneously release neurotransmitter leading to spontaneous generation of action potentials (spikes) in innervating afferent neurons. We analyzed spontaneous spike patterns recorded from the lateral line of zebrafish and found that distributions of interspike intervals (ISIs) either have an exponential shape or an "L" shape that is characterized by a sharp decay but wide tail. ISI data were fitted to renewal-process models that accounted for the neuron refractory periods and hair-cell synaptic release. Modeling the timing of synaptic release using a mixture of two exponential distributions yielded the best fit for our ISI data. Additionally, lateral line ISIs displayed positive serial correlation and appeared to exhibit switching between faster and slower modes of spike generation. This pattern contrasts with previous findings from the auditory system where ISIs tended to have negative serial correlation due to synaptic depletion. We propose that afferent neuron innervation with multiple and heterogenous hair-cells synapses, each influenced by changes in calcium domains, can serve as a mechanism for the random switching behavior. Overall, our analyses provide evidence of how physiological similarities and differences between synapses and innervation patterns in the auditory, vestibular, and lateral line systems can lead to variations in spontaneous activity.

Contribution of action potentials to the extracellular field potential in the nucleus laminaris of barn owl

Extracellular field potentials (EFP) are widely used to evaluate in vivo neural activity, but identification of multiple sources and their relative contributions is often ambiguous, making the interpretation of the EFP difficult. We have therefore analyzed a model EFP from a simple brainstem circuit with separable pre- and postsynaptic components to determine whether we could isolate its sources. Our previous papers had shown that the barn owl neurophonic largely originates with spikes from input axons and synapses that terminate on the neurons in the nucleus laminaris (NL) (Kuokkanen PT, Wagner H, Ashida G, Carr CE, Kempter R. J Neurophysiol 104: 2274-2290, 2010; Kuokkanen PT, Ashida G, Carr CE, Wagner H, Kempter R. J Neurophysiol 110: 117-130, 2013; McColgan T, Liu J, Kuokkanen PT, Carr CE, Wagner H, Kempter R. eLife 6: e26106, 2017). To determine how much the postsynaptic NL neurons contributed to the neurophonic, we recorded EFP responses in NL in vivo. Power spectral analyses showed that a small spectral component of the evoked response, between 200 and 700 Hz, could be attributed to the NL neurons' spikes, while nucleus magnocellularis (NM) spikes dominate the EFP at frequencies ≳1 kHz. Thus, spikes of NL neurons and NM axons contribute to the EFP in NL in distinct frequency bands. We conclude that if the spectral components of source types are different and if their activities can be selectively modulated, the identification of EFP sources is possible. NEW & NOTEWORTHY Extracellular field potentials (EFPs) generate clinically important signals, but their sources are incompletely understood. As a model, we have analyzed the auditory neurophonic in the barn owl's nucleus laminaris. There the EFP originates predominantly from spiking in the afferent axons, with spectral power ≳1 kHz, while postsynaptic laminaris neurons contribute little. In conclusion, the identification of EFP sources is possible if they have different spectral components and if their activities can be modulated selectively.

A model of synaptic vesicle-pool depletion and replenishment can account for the interspike interval distributions and nonrenewal properties of spontaneous spike trains of auditory-nerve fibers

In mammalian auditory systems, the spiking characteristics of each primary afferent (type I auditory-nerve fiber; ANF) are mainly determined by a single ribbon synapse in a single receptor cell (inner hair cell; IHC). ANF spike trains therefore provide a window into the operation of these synapses and cells. It was demonstrated previously (Heil et al., 2007) that the distribution of interspike intervals (ISIs) of cat ANFs during spontaneous activity can be modeled as resulting from refractoriness operating on a non-Poisson stochastic point process of excitation (transmitter release events from the IHC). Here, we investigate nonrenewal properties of these cat-ANF spontaneous spike trains, manifest as negative serial ISI correlations and reduced spike-count variability over short timescales. A previously discussed excitatory process, the constrained failure of events from a homogeneous Poisson point process, can account for these properties, but does not offer a parsimonious explanation for certain trends in the data. We then investigate a three-parameter model of vesicle-pool depletion and replenishment and find that it accounts for all experimental observations, including the ISI distributions, with only the release probability varying between spike trains. The maximum number of units (single vesicles or groups of simultaneously released vesicles) in the readily releasable pool and their replenishment time constant can be assumed to be constant (∼4 and 13.5 ms, respectively). We suggest that the organization of the IHC ribbon synapses not only enables sustained release of neurotransmitter but also imposes temporal regularity on the release process, particularly when operating at high rates.

Refractoriness enhances temporal coding by auditory nerve fibers

A universal property of spiking neurons is refractoriness, a transient decrease in discharge probability immediately following an action potential (spike). The refractory period lasts only one to a few milliseconds, but has the potential to affect temporal coding of acoustic stimuli by auditory neurons, which are capable of submillisecond spike-time precision. Here this possibility was investigated systematically by recording spike times from chicken auditory nerve fibers in vivo while stimulating with repeated pure tones at characteristic frequency. Refractory periods were tightly distributed, with a mean of 1.58 ms. A statistical model was developed to recapitulate each fiber's responses and then used to predict the effect of removing the refractory period on a cell-by-cell basis for two largely independent facets of temporal coding: faithful entrainment of interspike intervals to the stimulus frequency and precise synchronization of spike times to the stimulus phase. The ratio of the refractory period to the stimulus period predicted the impact of refractoriness on entrainment and synchronization. For ratios less than ∼0.9, refractoriness enhanced entrainment and this enhancement was often accompanied by an increase in spike-time precision. At higher ratios, little or no change in entrainment or synchronization was observed. Given the tight distribution of refractory periods, the ability of refractoriness to improve temporal coding is restricted to neurons responding to low-frequency stimuli. Enhanced encoding of low frequencies likely affects sound localization and pitch perception in the auditory system, as well as perception in nonauditory sensory modalities, because all spiking neurons exhibit refractoriness.

A new method to infer higher-order spike correlations from membrane potentials

What is the role of higher-order spike correlations for neuronal information processing? Common data analysis methods to address this question are devised for the application to spike recordings from multiple single neurons. Here, we present a new method which evaluates the subthreshold membrane potential fluctuations of one neuron, and infers higher-order correlations among the neurons that constitute its presynaptic population. This has two important advantages: Very large populations of up to several thousands of neurons can be studied, and the spike sorting is obsolete. Moreover, this new approach truly emphasizes the functional aspects of higher-order statistics, since we infer exactly those correlations which are seen by a neuron. Our approach is to represent the subthreshold membrane potential fluctuations as presynaptic activity filtered with a fixed kernel, as it would be the case for a leaky integrator neuron model. This allows us to adapt the recently proposed method CuBIC (cumulant based inference of higher-order correlations from the population spike count; Staude et al., J Comput Neurosci 29(1-2):327-350, 2010c) with which the maximal order of correlation can be inferred. By numerical simulation we show that our new method is reasonably sensitive to weak higher-order correlations, and that only short stretches of membrane potential are required for their reliable inference. Finally, we demonstrate its remarkable robustness against violations of the simplifying assumptions made for its construction, and discuss how it can be employed to analyze in vivo intracellular recordings of membrane potentials.

Mechanism of spontaneous activity in afferent neurons of the zebrafish lateral-line organ

Many auditory, vestibular, and lateral-line afferent neurons display spontaneous action potentials. This spontaneous spiking is thought to result from hair-cell glutamate release in the absence of stimuli. Spontaneous release at hair-cell resting potentials presumably results from Ca(V)1.3 L-type calcium channel activity. Here, using intact zebrafish larvae, we recorded robust spontaneous spiking from lateral-line afferent neurons in the absence of external stimuli. Consistent with the above assumptions, spiking was absent in mutants that lacked either Vesicular glutamate transporter 3 (Vglut3) or Ca(V)1.3. We then tested the hypothesis that spontaneous spiking resulted from sustained Ca(V)1.3 activity due to depolarizing currents that are active at rest. Mechanotransduction currents (I(MET)) provide a depolarizing influence to the resting potential. However, following block of I(MET), spontaneous spiking persisted and was characterized by longer interspike intervals and increased periods of inactivity. These results suggest that an additional depolarizing influence maintains the resting potential within the activation range of Ca(V)1.3. To test whether the hyperpolarization-activated cation current, I(h) participates in setting the resting potential, we applied I(h) antagonists. Both ZD7288 and DK-AH 269 reduced spontaneous activity. Finally, concomitant block of I(MET) and I(h) essentially abolished spontaneous activity, ostensibly by hyperpolarization outside of the activation range for Ca(V)1.3. Together, our data support a mechanism for spontaneous spiking that results from Ca(2+)-dependent neurotransmitter release at hair-cell resting potentials that are maintained within the activation range of Ca(V)1.3 channels through active I(MET) and I(h).