Which spike train distance is most suitable for distinguishing rate and temporal coding?

BACKGROUND

It is commonly assumed in neuronal coding that repeated presentations of a stimulus to a coding neuron elicit similar responses. One common way to assess similarity are spike train distances. These can be divided into spike-resolved, such as the Victor-Purpura and the van Rossum distance, and time-resolved, e.g. the ISI-, the SPIKE- and the RI-SPIKE-distance.

NEW METHOD

We use independent steady-rate Poisson processes as surrogates for spike trains with fixed rate and no timing information to address two basic questions: How does the sensitivity of the different spike train distances to temporal coding depend on the rates of the two processes and how do the distances deal with very low rates?

RESULTS

Spike-resolved distances always contain rate information even for parameters indicating time coding. This is an issue for reasonably high rates but beneficial for very low rates. In contrast, the operational range for detecting time coding of time-resolved distances is superior at normal rates, but these measures produce artefacts at very low rates. The RI-SPIKE-distance is the only measure that is sensitive to timing information only.

COMPARISON WITH EXISTING METHODS

While our results on rate-dependent expectation values for the spike-resolved distances agree with Chicharro et al. (2011), we here go one step further and specifically investigate applicability for very low rates.

CONCLUSIONS

The most appropriate measure depends on the rates of the data being analysed. Accordingly, we summarize our results in one table that allows an easy selection of the preferred measure for any kind of data.

Measures of spike train synchrony for data with multiple time scales

BACKGROUND

Measures of spike train synchrony are widely used in both experimental and computational neuroscience. Time-scale independent and parameter-free measures, such as the ISI-distance, the SPIKE-distance and SPIKE-synchronization, are preferable to time scale parametric measures, since by adapting to the local firing rate they take into account all the time scales of a given dataset.

NEW METHOD

In data containing multiple time scales (e.g. regular spiking and bursts) one is typically less interested in the smallest time scales and a more adaptive approach is needed. Here we propose the A-ISI-distance, the A-SPIKE-distance and A-SPIKE-synchronization, which generalize the original measures by considering the local relative to the global time scales. For the A-SPIKE-distance we also introduce a rate-independent extension called the RIA-SPIKE-distance, which focuses specifically on spike timing.

RESULTS

The adaptive generalizations A-ISI-distance and A-SPIKE-distance allow to disregard spike time differences that are not relevant on a more global scale. A-SPIKE-synchronization does not any longer demand an unreasonably high accuracy for spike doublets and coinciding bursts. Finally, the RIA-SPIKE-distance proves to be independent of rate ratios between spike trains.

COMPARISON WITH EXISTING METHODS

We find that compared to the original versions the A-ISI-distance and the A-SPIKE-distance yield improvements for spike trains containing different time scales without exhibiting any unwanted side effects in other examples. A-SPIKE-synchronization matches spikes more efficiently than SPIKE-synchronization.

CONCLUSIONS

With these proposals we have completed the picture, since we now provide adaptive generalized measures that are sensitive to firing rate only (A-ISI-distance), to timing only (ARI-SPIKE-distance), and to both at the same time (A-SPIKE-distance).

Synchronization study in ring-like and grid-like neuronal networks

In this paper, we study the synchronization status of both two gap-junction coupled neurons and neuronal network with two different network connectivity patterns. One of the network connectivity patterns is a ring-like neuronal network, which only considers nearest-neighbor neurons. The other is a grid-like neuronal network, with all nearest neighbor couplings. We show that by varying some key parameters, such as the coupling strength and the external current injection, the neuronal network will exhibit various patterns of firing synchronization. Different types of firing synchronization are diagnosed by means of a mean field potential, a bifurcation diagram, a correlation coefficient and the ISI-distance method. Numerical simulations demonstrate that the synchronization status of multiple neurons is much dependent on the network patters, when the number of neurons is the same. It is also demonstrated that the synchronization status of two coupled neurons is similar with the grid-like neuronal network, but differs radically from that of the ring-like neuronal network. These results may be instructive in understanding synchronization transitions in neuronal systems.

Using interspike intervals to quantify noise effects on spike trains in temperature encoding neurons

This paper examines how noise interacts with the non-linear dynamical mechanisms of neuronal stimulus. We study the spike trains generated by a minimal Hodgkin-Huxley type model of a cold receptor neuron. The distributions of interspike intervals(ISIs) of purely deterministic simulations exhibit considerable differences compared to the noisy ones. We quantify the effect of noise using ISI return plots and the ISI-distance recently proposed by Kreuz et al. (J Neurosci Meth, 165:151-161, 2007). It is shown that the spike trains of a cold receptor neuron are more strongly affected by noise for low temperatures than for high temperatures. This trend is also observed in both regimes of cold receptors: tonic firing(which occurs for low and high temperatures) and bursting (which occurs for intermediate temperatures).