SPIKY: a graphical user interface for monitoring spike train synchrony.

Techniques for recording large-scale neuronal spiking activity are developing very fast. This leads to an increasing demand for algorithms capable of analyzing large amounts of experimental spike train data. One of the most crucial and demanding tasks is the identification of similarity patterns with a very high temporal resolution and across different spatial scales. To address this task, in recent years three time-resolved measures of spike train synchrony have been proposed, the ISI-distance, the SPIKE-distance, and event synchronization. The Matlab source codes for calculating and visualizing these measures have been made publicly available. However, due to the many different possible representations of the results the use of these codes is rather complicated and their application requires some basic knowledge of Matlab. Thus it became desirable to provide a more user-friendly and interactive interface. Here we address this need and present SPIKY, a graphical user interface that facilitates the application of time-resolved measures of spike train synchrony to both simulated and real data. SPIKY includes implementations of the ISI-distance, the SPIKE-distance, and the SPIKE-synchronization (an improved and simplified extension of event synchronization) that have been optimized with respect to computation speed and memory demand. It also comprises a spike train generator and an event detector that makes it capable of analyzing continuous data. Finally, the SPIKY package includes additional complementary programs aimed at the analysis of large numbers of datasets and the estimation of significance levels.

A Measure of Concurrent Neural Firing Activity Based on Mutual Information.

Multiple methods have been developed in an attempt to quantify stimulus-induced neural coordination and to understand internal coordination of neuronal responses by examining the synchronization phenomena in neural discharge patterns. In this work we propose a novel approach to estimate the degree of concomitant firing between two neural units, based on a modified form of mutual information (MI) applied to a two-state representation of the firing activity. The binary profile of each single unit unfolds its discharge activity in time by decomposition into the state of neural quiescence/low activity and state of moderate firing/bursting. Then, the MI computed between the two binary streams is normalized by their minimum entropy and is taken as positive or negative depending on the prevalence of identical or opposite concomitant states. The resulting measure, denoted as Concurrent Firing Index based on MI (CFIMI), relies on a single input parameter and is otherwise assumption-free and symmetric. Exhaustive validation was carried out through controlled experiments in three simulation scenarios, showing that CFIMI is independent on firing rate and recording duration, and is sensitive to correlated and anti-correlated firing patterns. Its ability to detect non-correlated activity was assessed using ad-hoc surrogate data. Moreover, the evaluation of CFIMI on experimental recordings of spiking activity in retinal ganglion cells brought insights into the changes of neural synchrony over time. The proposed measure offers a novel perspective on the estimation of neural synchrony, providing information on the co-occurrence of firing states in the two analyzed trains over longer temporal scales compared to existing measures.

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

Time-frequency characterization of local field potential in a decision making task.

This study seeks to characterize the neuronal mechanisms underlying voluntary decisions to check/verify. In order to describe and potentially decode decisions from brain signals we analyzed intracortical recordings from monkey prefrontal regions obtained during a cognitive task requiring self-initiated as well as cue-instructed decisions. Using local field potentials (LFP) and single units, we analyzed power spectral density, oscillatory modes, power profiles in time, single unit firing rate, and spike-phase relationships in the ß band. Our results point toward specific but variable activation patterns of oscillations in ß band from separate recordings, with task-dependent frequency preference and amplitude modulation of power. The results suggest relationships between particular LFP oscillations and functions engaged at specific time in the task.