A new similarity measure for spike trains: sensitivity to bursts and periods of inhibition

High-resolution raindrop counting via instantaneous frequency sensing on hydrophobic elastic membranes

In this paper, we introduce a novel approach that paves the way for the creation of affordable, high-precision rainfall sensors utilizing microphone data. The cornerstone of this methodology is an innovative algorithm capable of converting audio recordings into distinctive features, which are subsequently processed by a compact machine learning model. Our findings demonstrate that this technique can attain a temporal resolution of 10 milliseconds with an accuracy of 80%, underscoring its potential to overcome the limitations imposed by the necessity for power infrastructure and specialized expertise in traditional rain sensing methods.

STSimM: A new tool for evaluating neuron model performance and detecting spike trains similarity

BACKGROUND

In computational neuroscience, performance measures are essential for quantitatively assessing the predictive power of neuron models, while similarity measures are used to estimate the level of synchrony between two or more spike trains. Most of the measures proposed in the literature require setting an appropriate time-scale and often neglect silent periods.

NEW METHOD

Four time-scale adaptive performance and similarity measures are proposed and implemented in the STSimM (Spike Trains Similarity Measures) Python tool. These measures are designed to accurately capture both the precise timing of individual spikes and shared periods of inactivity among spike trains.

RESULTS

The proposed ST-measures demonstrate enhanced sensitivity over Spike-contrast and SPIKE-distance in detecting spike train similarity, aligning closely with SPIKE-synchronization. Correlations among all similarity measures were observed in Poisson datasets, whereas in vivo-like synaptic stimulations showed correlations only between ST-measures and SPIKE-synchronization.

COMPARISON OF EXISTING METHOD

The STSimM measures are compared with SPIKE-distance, SPIKE-synchronization and Spike-contrast using four spike train datasets with varying similarity levels.

CONCLUSION

ST-measures appear more suitable for detecting both the precise timing of single spikes and shared periods of inactivity among spike trains compared to those considered in this work. Their flexibility originates from two primary factors: firstly, the inclusion of four key measures - ST-Accuracy, ST-Precision, ST-Recall, ST-Fscore - capable of discerning similarity levels across neuronal activity, whether interleaved with silent periods or solely focusing on spike timing accuracy; secondly, the integration of three model parameters that govern both precise spike detection and the weighting of silent periods.

Neural burst codes disguised as rate codes

The burst coding hypothesis posits that the occurrence of sudden high-frequency patterns of action potentials constitutes a salient syllable of the neural code. Many neurons, however, do not produce clearly demarcated bursts, an observation invoked to rule out the pervasiveness of this coding scheme across brain areas and cell types. Here we ask how detrimental ambiguous spike patterns, those that are neither clearly bursts nor isolated spikes, are for neuronal information transfer. We addressed this question using information theory and computational simulations. By quantifying how information transmission depends on firing statistics, we found that the information transmitted is not strongly influenced by the presence of clearly demarcated modes in the interspike interval distribution, a feature often used to identify the presence of burst coding. Instead, we found that neurons having unimodal interval distributions were still able to ascribe different meanings to bursts and isolated spikes. In this regime, information transmission depends on dynamical properties of the synapses as well as the length and relative frequency of bursts. Furthermore, we found that common metrics used to quantify burstiness were unable to predict the degree with which bursts could be used to carry information. Our results provide guiding principles for the implementation of coding strategies based on spike-timing patterns, and show that even unimodal firing statistics can be consistent with a bivariate neural code.

T-pattern analysis and spike train dissimilarity for the analysis of structure in blinking behavior

Spontaneous eye-blinks are a ubiquitous behavior. However, blink timing is not random, nor does it always follow physiological demands. Research rather suggests that blink timing, and thus the structure of blinking behavior, is influenced by cognitive processes, such as attention. Since attention is regarded a necessary precursor of media use phenomena, the present study investigates the relation between the structure of blinking behavior and the media use phenomenon of spatial presence. To this end, spontaneous eye-blinks have been observed in an experiment during the reception of a video story. The methods of T-pattern analysis, ISI distance, and IBI variability have been used to quantify stimulus-dependent blink structure, which has then been related to self-reports of spatial presence experiences. While the T-pattern analysis and ISI distance showed converging results for behavior structure, a hypothesized relation between more stimulus-dependent blink structure and stronger presence experiences was not found. On the contrary, blink data suggested a difference in attention allocation, whereas self-report data indicated no difference in presence experiences. This demonstrates that beyond self-report and the analysis of event frequencies, the analysis of behavior structure offers insights into behavior synchronization between participants, allowing for new inferences on internal processing of media stimuli.

Dynamical networks: Finding, measuring, and tracking neural population activity using network science

Systems neuroscience is in a headlong rush to record from as many neurons at the same time as possible. As the brain computes and codes using neuron populations, it is hoped these data will uncover the fundamentals of neural computation. But with hundreds, thousands, or more simultaneously recorded neurons come the inescapable problems of visualizing, describing, and quantifying their interactions. Here I argue that network science provides a set of scalable, analytical tools that already solve these problems. By treating neurons as nodes and their interactions as links, a single network can visualize and describe an arbitrarily large recording. I show that with this description we can quantify the effects of manipulating a neural circuit, track changes in population dynamics over time, and quantitatively define theoretical concepts of neural populations such as cell assemblies. Using network science as a core part of analyzing population recordings will thus provide both qualitative and quantitative advances to our understanding of neural computation.

Detecting joint pausiness in parallel spike trains

BACKGROUND

Transient periods with reduced neuronal discharge - called 'pauses' - have recently gained increasing attention. In dopamine neurons, pauses are considered important teaching signals, encoding negative reward prediction errors. Particularly simultaneous pauses are likely to have increased impact on information processing.

COMPARISON WITH EXISTING METHODS

Available methods for detecting joint pausing analyze temporal overlap of pauses across spike trains. Such techniques are threshold dependent and can fail to identify joint pauses that are easily detectable by eye, particularly in spike trains with different firing rates.

NEW METHOD

We introduce a new statistic called pausiness that measures the degree of synchronous pausing in spike train pairs and avoids threshold-dependent identification of specific pauses. A new graphic termed the cross-pauseogram compares the joint pausiness of two spike trains with its time shifted analogue, such that a (pausiness) peak indicates joint pausing. When assessing significance of pausiness peaks, we use a stochastic model with synchronous spikes to disentangle joint pausiness arising from synchronous spikes from additional 'joint excess pausiness' (JEP). Parameter estimates are obtained from auto- and cross-correlograms, and statistical significance is assessed by comparison to simulated cross-pauseograms.

RESULTS

Our new method was applied to dopamine neuron pairs recorded in the ventral tegmental area of awake behaving mice. Significant JEP was detected in about 20% of the pairs.

CONCLUSION

Given the neurophysiological importance of pauses and the fact that neurons integrate multiple inputs, our findings suggest that the analysis of JEP can reveal interesting aspects in the activity of simultaneously recorded neurons.

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 of action potentials during low-magnesium-induced bursting

The relationship between mono- and polysynaptic strength and action potential synchronization was explored using a reduced external Mg(2+) model. Single and dual whole cell patch-clamp recordings were performed in hippocampal cultures in three concentrations of external Mg(2+). In decreased Mg(2+) medium, the individual cells transitioned to spontaneous bursting behavior. In lowered Mg(2+) media the larger excitatory synaptic events were observed more frequently and fewer transmission failures occurred, suggesting strengthened synaptic transmission. The event synchronization was calculated for the neural action potentials of the cell pairs, and it increased in media where Mg(2+) concentration was lowered. Analysis of surrogate data where bursting was present, but no direct or indirect connections existed between the neurons, showed minimal action potential synchronization. This suggests the synchronization of action potentials is a product of the strengthening synaptic connections within neuronal networks.

Monitoring Spiking Activity of Many Individual Neurons in Invertebrate Ganglia

Optical recording with fast voltage sensitive dyes makes it possible, in suitable preparations, to simultaneously monitor the action potentials of large numbers of individual neurons. Here we describe methods for doing this, including considerations of different dyes and imaging systems, methods for correlating the optical signals with their source neurons, procedures for getting good signals, and the use of Independent Component Analysis for spike-sorting raw optical data into single neuron traces. These combined tools represent a powerful approach for large-scale recording of neural networks with high temporal and spatial resolution.

Recent developments in VSD imaging of small neuronal networks

Voltage-sensitive dye (VSD) imaging is a powerful technique that can provide, in single experiments, a large-scale view of network activity unobtainable with traditional sharp electrode recording methods. Here we review recent work using VSDs to study small networks and highlight several results from this approach. Topics covered include circuit mapping, network multifunctionality, the network basis of decision making, and the presence of variably participating neurons in networks. Analytical tools being developed and applied to large-scale VSD imaging data sets are discussed, and the future prospects for this exciting field are considered.

Revealing cell assemblies at multiple levels of granularity

BACKGROUND

Current neuronal monitoring techniques, such as calcium imaging and multi-electrode arrays, enable recordings of spiking activity from hundreds of neurons simultaneously. Of primary importance in systems neuroscience is the identification of cell assemblies: groups of neurons that cooperate in some form within the recorded population.

NEW METHOD

We introduce a simple, integrated framework for the detection of cell-assemblies from spiking data without a priori assumptions about the size or number of groups present. We define a biophysically-inspired measure to extract a directed functional connectivity matrix between both excitatory and inhibitory neurons based on their spiking history. The resulting network representation is analyzed using the Markov Stability framework, a graph theoretical method for community detection across scales, to reveal groups of neurons that are significantly related in the recorded time-series at different levels of granularity.

RESULTS AND COMPARISON WITH EXISTING METHODS

Using synthetic spike-trains, including simulated data from leaky-integrate-and-fire networks, our method is able to identify important patterns in the data such as hierarchical structure that are missed by other standard methods. We further apply the method to experimental data from retinal ganglion cells of mouse and salamander, in which we identify cell-groups that correspond to known functional types, and to hippocampal recordings from rats exploring a linear track, where we detect place cells with high fidelity.

CONCLUSIONS

We present a versatile method to detect neural assemblies in spiking data applicable across a spectrum of relevant scales that contributes to understanding spatio-temporal information gathered from systems neuroscience experiments.

Respiratory cycle entrainment of septal neurons mediates the fast coupling of sniffing rate and hippocampal theta rhythm

Memory for odour information may result from temporal coupling between the olfactory and hippocampal systems. Respiration defines the frequency of olfactory perception, but how the respiratory rate affects hippocampal oscillations remains poorly understood. The afferent connectivity of the medial septum/diagonal band of Broca complex (MS/DB) proposes this region as a crossroads between respiratory and limbic pathways. Here we investigate if the firing rates of septal neurons integrate respiratory rate signals. We demonstrate that approximately 50% of MS/DB neurons are temporally correlated with sniffing frequency. Moreover, a group of slow-spiking septal neurons are phase-locked to the sniffing cycle. We show that inter-burst intervals of MS/DB theta cells relate to the sniff rate. Intranasal odour infusion evokes sniff phase preference for the activity of fast-spiking MS/DB neurons. Concurrently, the infusion augments the correlation between sniffing and limbic theta oscillations. During periods of sniffing–theta correlation, CA1 place cells fired preferentially during the inhalation phase, suggesting the theta cycle as a coherent time frame for central olfactory processing. Furthermore, injection of the GABAergic agonist muscimol into medial septum induces a parallel decrease of sniffing and theta frequencies. Our findings provide experimental evidence that MS/DB does not merely generate theta rhythm, but actively integrates sensorimotor stimuli that reflect sniffing rate. Such integration may provide temporal oscillatory synchronisation of MS/DB-innervated limbic structures with the sniffing cycle.

Synchronous firing of antennal-lobe projection neurons encodes the behaviorally effective ratio of sex-pheromone components in male Manduca sexta

Olfactory stimuli that are essential to an animal's survival and reproduction are often complex mixtures of volatile organic compounds in characteristic proportions. Here, we investigated how these proportions are encoded in the primary olfactory processing center, the antennal lobe, of male Manduca sexta moths. Two key components of the female's sex pheromone, present in an approximately 2:1 ratio, are processed in each of two neighboring glomeruli in the macroglomerular complex (MGC) of males of this species. In wind-tunnel flight experiments, males exhibited behavioral selectivity for ratios approximating the ratio released by conspecific females. The ratio between components was poorly represented, however, in the firing-rate output of uniglomerular MGC projection neurons (PNs). PN firing rate was mostly insensitive to the ratio between components, and individual PNs did not exhibit a preference for a particular ratio. Recording simultaneously from pairs of PNs in the same glomerulus, we found that the natural ratio between components elicited the most synchronous spikes, and altering the proportion of either component decreased the proportion of synchronous spikes. The degree of synchronous firing between PNs in the same glomerulus thus selectively encodes the natural ratio that most effectively evokes the natural behavioral response to pheromone.

Reliability of spike and burst firing in thalamocortical relay cells

The reliability and precision of the timing of spikes in a spike train is an important aspect of neuronal coding. We investigated reliability in thalamocortical relay (TCR) cells in the acute slice and also in a Morris-Lecar model with several extensions. A frozen Gaussian noise current, superimposed on a DC current, was injected into the TCR cell soma. The neuron responded with spike trains that showed trial-to-trial variability, due to amongst others slow changes in its internal state and the experimental setup. The DC current allowed to bring the neuron in different states, characterized by a well defined membrane voltage (between -80 and -50 mV) and by a specific firing regime that on depolarization gradually shifted from a predominantly bursting regime to a tonic spiking regime. The filtered frozen white noise generated a spike pattern output with a broad spike interval distribution. The coincidence factor and the Hunter and Milton measure were used as reliability measures of the output spike train. In the experimental TCR cell as well as the Morris-Lecar model cell the reliability depends on the shape (steepness) of the current input versus spike frequency output curve. The model also allowed to study the contribution of three relevant ionic membrane currents to reliability: a T-type calcium current, a cation selective h-current and a calcium dependent potassium current in order to allow bursting, investigate the consequences of a more complex current-frequency relation and produce realistic firing rates. The reliability of the output of the TCR cell increases with depolarization. In hyperpolarized states bursts are more reliable than single spikes. The analytically derived relations were capable to predict several of the experimentally recorded spike features.

Detecting connectivity changes in neuronal networks

We develop a method from semiparametric statistics (Cox, 1972) for the purpose of tracking links and connection strengths over time in a neuronal network from spike train data. We consider application of the method as implemented in Masud and Borisyuk (2011), and evaluate its use on data generated independently of the Cox model hypothesis, in particular from a spiking model of Izhikevich in four different dynamical regimes. Then, we show how the Cox method can be used to determine statistically significant changes in network connectivity over time. Our methodology is demonstrated using spike trains from multi-electrode array measurements of networks of cultured mammalian spinal cord cells.