Finding the event structure of neuronal spike trains

Asynchronous Spatiotemporal Spike Metric for Event Cameras

Event cameras as bioinspired vision sensors have shown great advantages in high dynamic range and high temporal resolution in vision tasks. Asynchronous spikes from event cameras can be depicted using the marked spatiotemporal point processes (MSTPPs). However, how to measure the distance between asynchronous spikes in the MSTPPs still remains an open issue. To address this problem, we propose a general asynchronous spatiotemporal spike metric considering both spatiotemporal structural properties and polarity attributes for event cameras. Technically, the conditional probability density function is first introduced to describe the spatiotemporal distribution and polarity prior in the MSTPPs. Besides, a spatiotemporal Gaussian kernel is defined to capture the spatiotemporal structure, which transforms discrete spikes into the continuous function in a reproducing kernel Hilbert space (RKHS). Finally, the distance between asynchronous spikes can be quantified by the inner product in the RKHS. The experimental results demonstrate that the proposed approach outperforms the state-of-the-art methods and achieves significant improvement in computational efficiency. Especially, it is able to better depict the changes involving spatiotemporal structural properties and polarity attributes.

Causal Network Inference for Neural Ensemble Activity

Interactions among cellular components forming a mesoscopic scale brain network (microcircuit) display characteristic neural dynamics. Analysis of microcircuits provides a system-level understanding of the neurobiology of health and disease. Causal discovery aims to detect causal relationships among variables based on observational data. A key barrier in causal discovery is the high dimensionality of the variable space. A method called Causal Inference for Microcircuits (CAIM) is proposed to reconstruct causal networks from calcium imaging or electrophysiology time series. CAIM combines neural recording, Bayesian network modeling, and neuron clustering. Validation experiments based on simulated data and a real-world reaching task dataset demonstrated that CAIM accurately revealed causal relationships among neural clusters.

Adaptive spike-artifact removal from local field potentials uncovers prominent beta and gamma band neuronal synchronization

BACKGROUND

Many neurons synchronize their action potentials to the phase of local field potential (LFP) fluctuations in one or more frequency bands. Analyzing this spike-to-LFP synchronization is challenging, however, when neural spikes and LFP are generated in the same local circuit, because the spike's action potential waveform leak into the LFP and distort phase synchrony estimates. Existing approaches to address this spike bleed-through artifact relied on removing the average action potential waveforms of neurons, but this leaves artifacts in the LFP and distorts synchrony estimates.

NEW METHOD

We describe a spike-removal method that surpasses these limitations by decomposing individual action potentials into their frequency components before their removal from the LFP. The adaptively estimated frequency components allow for variable spread, strength and temporal variation of the spike artifact.

RESULTS

This adaptive approach effectively removes spike bleed-through artifacts in synthetic data with known ground truth, and in single neuron and LFP recordings in nonhuman primate striatum. For a large population of neurons with both narrow and broad action potential waveforms, the use of adaptive artifact removal uncovered 20-35 Hz beta and 35-45 Hz gamma band spike-LFP synchronization that would have remained contaminated otherwise.

COMPARISON WITH EXISTING METHODS

We demonstrate that adaptive spike-artifact removal cleans LFP data that remained contaminated when applying existing Bayesian and non-Bayesian methods of average spike-artifact removal.

CONCLUSIONS

Applying adaptive spike-removal from field potentials allows to estimate the phase at which neurons synchronize and the consistency of their phase-locked firing for both beta and low gamma frequencies. These metrics may prove essential to understand cell-to-circuit neuronal interactions in multiple brain systems.

A Lower Bound for the First Passage Time Density of the Suprathreshold Ornstein-Uhlenbeck Process

We prove that the first passage time density ρ(t) for an Ornstein-Uhlenbeck process X(t) obeying dX = -β Xdt + σdW to reach a fixed threshold θ from a suprathreshold initial condition x0 > θ > 0 has a lower bound of the form ρ(t) > kexp[-pe6βt] for positive constants k and p for times t exceeding some positive value u. We obtain explicit expressions for k, p, and u in terms of β, σ, x0, and θ, and discuss the application of the results to the synchronization of periodically forced stochastic leaky integrate-and-fire model neurons.

Separating burst from background spikes in multichannel neuronal recordings using return map analysis

We propose a preprocessing method to separate coherent neuronal network activity, referred to as “bursts”, from background spikes. High background activity in neuronal recordings reduces the effectiveness of currently available burst detection methods. For long-term, stationary recordings, burst and background spikes have a bimodal ISI distribution which makes it easy to select the threshold to separate burst and background spikes. Finite, nonstationary recordings lead to noisy ISIs for which the bimodality is not that clear. We introduce a preprocessing method to separate burst from background spikes to improve burst detection reliability because it efficiently uses both single and multichannel activity. The method is tested using a stochastic model constrained by data available in the literature and recordings from primary cortical neurons cultured on multielectrode arrays. The separation between burst and background spikes is obtained using the interspike interval return map. The cutoff threshold is the key parameter to separate the burst and background spikes. We compare two methods for selecting the threshold. The 2-step method, in which threshold selection is based on fixed heuristics. The iterative method, in which the optimal cutoff threshold is directly estimated from the data. The proposed preprocessing method significantly increases the reliability of several established burst detection algorithms, both for simulated and real recordings. The preprocessing method makes it possible to study the effects of diseases or pharmacological manipulations, because it can deal efficiently with nonstationarity in the data.

A transition to sharp timing in stochastic leaky integrate-and-fire neurons driven by frozen noisy input

The firing activity of intracellularly stimulated neurons in cortical slices has been demonstrated to be profoundly affected by the temporal structure of the injected current (Mainen & Sejnowski, 1995 ). This suggests that the timing features of the neural response may be controlled as much by its own biophysical characteristics as by how a neuron is wired within a circuit. Modeling studies have shown that the interplay between internal noise and the fluctuations of the driving input controls the reliability and the precision of neuronal spiking (Cecchi et al., 2000 ; Tiesinga, 2002 ; Fellous, Rudolph, Destexhe, & Sejnowski, 2003 ). In order to investigate this interplay, we focus on the stochastic leaky integrate-and-fire neuron and identify the Hölder exponent H of the integrated input as the key mathematical property dictating the regime of firing of a single-unit neuron. We have recently provided numerical evidence (Taillefumier & Magnasco, 2013 ) for the existence of a phase transition when [Formula: see text] becomes less than the statistical Hölder exponent associated with internal gaussian white noise (H=1/2). Here we describe the theoretical and numerical framework devised for the study of a neuron that is periodically driven by frozen noisy inputs with exponent H>0. In doing so, we account for the existence of a transition between two regimes of firing when H=1/2, and we show that spiking times have a continuous density when the Hölder exponent satisfies H>1/2. The transition at H=1/2 formally separates rate codes, for which the neural firing probability varies smoothly, from temporal codes, for which the neuron fires at sharply defined times regardless of the intensity of internal noise.

A phase transition in the first passage of a Brownian process through a fluctuating boundary with implications for neural coding

Finding the first time a fluctuating quantity reaches a given boundary is a deceptively simple-looking problem of vast practical importance in physics, biology, chemistry, neuroscience, economics, and industrial engineering. Problems in which the bound to be traversed is itself a fluctuating function of time include widely studied problems in neural coding, such as neuronal integrators with irregular inputs and internal noise. We show that the probability p(t) that a Gauss-Markov process will first exceed the boundary at time t suffers a phase transition as a function of the roughness of the boundary, as measured by its Hölder exponent H. The critical value occurs when the roughness of the boundary equals the roughness of the process, so for diffusive processes the critical value is Hc = 1/2. For smoother boundaries, H > 1/2, the probability density is a continuous function of time. For rougher boundaries, H < 1/2, the probability is concentrated on a Cantor-like set of zero measure: the probability density becomes divergent, almost everywhere either zero or infinity. The critical point Hc = 1/2 corresponds to a widely studied case in the theory of neural coding, in which the external input integrated by a model neuron is a white-noise process, as in the case of uncorrelated but precisely balanced excitatory and inhibitory inputs. We argue that this transition corresponds to a sharp boundary between rate codes, in which the neural firing probability varies smoothly, and temporal codes, in which the neuron fires at sharply defined times regardless of the intensity of internal noise.

Modifying spiking precision in conductance-based neuronal models

The temporal precision of a neuron's spiking can be characterized by calculating its "jitter," defined as the standard deviation of the timing of individual spikes in response to repeated presentations of a stimulus. Sub-millisecond jitters have been measured for neurons in a variety of experimental systems and appear to be functionally important in some instances. We have investigated how modifying a neuron's maximal conductances affects jitter using the leaky integrate-and-fire (LIF) model and an eight-conductance Hodgkin-Huxley type (HH8) model. We observed that jitter can be largely understood in the LIF model in terms of the neuron's filtering properties. In the HH8 model we found the role of individual conductances in determining jitter to be complicated and dependent on the model's spiking properties. Distinct behaviors were observed for populations with slow (<11.5 Hz) and fast (>11.5 Hz) spike rates and appear to be related to differences in a particular channel's activity at times just before spiking occurs.

Exact Event-Driven Implementation for Recurrent Networks of Stochastic Perfect Integrate-and-Fire Neurons

In vivo cortical recording reveals that indirectly driven neural assemblies can produce reliable and temporally precise spiking patterns in response to stereotyped stimulation. This suggests that despite being fundamentally noisy, the collective activity of neurons conveys information through temporal coding. Stochastic integrate-and-fire models delineate a natural theoretical framework to study the interplay of intrinsic neural noise and spike timing precision. However, there are inherent difficulties in simulating their networks’ dynamics in silico with standard numerical discretization schemes. Indeed, the well-posedness of the evolution of such networks requires temporally ordering every neuronal interaction, whereas the order of interactions is highly sensitive to the random variability of spiking times. Here, we answer these issues for perfect stochastic integrate-and-fire neurons by designing an exact event-driven algorithm for the simulation of recurrent networks, with delayed Dirac-like interactions. In addition to being exact from the mathematical standpoint, our proposed method is highly efficient numerically. We envision that our algorithm is especially indicated for studying the emergence of polychronized motifs in networks evolving under spike-timing-dependent plasticity with intrinsic noise.

Fluctuating inhibitory inputs promote reliable spiking at theta frequencies in hippocampal interneurons

Theta-frequency (4–12 Hz) rhythms in the hippocampus play important roles in learning and memory. CA1 interneurons located at the stratum lacunosum-moleculare and radiatum junction (LM/RAD) are thought to contribute to hippocampal theta population activities by rhythmically pacing pyramidal cells with inhibitory postsynaptic potentials. This implies that LM/RAD cells need to fire reliably at theta frequencies in vivo. To determine whether this could occur, we use biophysically based LM/RAD model cells and apply different cholinergic and synaptic inputs to simulate in vivo-like network environments. We assess spike reliabilities and spiking frequencies, identifying biophysical properties and network conditions that best promote reliable theta spiking. We find that synaptic background activities that feature large inhibitory, but not excitatory, fluctuations are essential. This suggests that strong inhibitory input to these cells is vital for them to be able to contribute to population theta activities. Furthermore, we find that Type I-like oscillator models produced by augmented persistent sodium currents (I_NaP) or diminished A-type potassium currents (I_A) enhance reliable spiking at lower theta frequencies. These Type I-like models are also the most responsive to large inhibitory fluctuations and can fire more reliably under such conditions. In previous work, we showed that I_NaP and I_A are largely responsible for establishing LM/RAD cells’ subthreshold activities. Taken together with this study, we see that while both these currents are important for subthreshold theta fluctuations and reliable theta spiking, they contribute in different ways – I_NaP to reliable theta spiking and subthreshold activity generation, and I_A to subthreshold activities at theta frequencies. This suggests that linking subthreshold and suprathreshold activities should be done with consideration of both in vivo contexts and biophysical specifics.

On the efficient calculation of van Rossum distances

The van Rossum metric measures the distance between two spike trains. Measuring a single van Rossum distance between one pair of spike trains is not a computationally expensive task, however, many applications require a matrix of distances between all the spike trains in a set or the calculation of a multi-neuron distance between two populations of spike trains. Moreover, often these calculations need to be repeated for many different parameter values. An algorithm is presented here to render these calculation less computationally expensive, making the complexity linear in the number of spikes rather than quadratic.

Quantifying statistical interdependence, part III: N > 2 point processes

Stochastic event synchrony (SES) is a recently proposed family of similarity measures. First, "events" are extracted from the given signals; next, one tries to align events across the different time series. The better the alignment, the more similar the N time series are considered to be. The similarity measures quantify the reliability of the events (the fraction of "nonaligned" events) and the timing precision. So far, SES has been developed for pairs of one-dimensional (Part I) and multidimensional (Part II) point processes. In this letter (Part III), SES is extended from pairs of signals to N > 2 signals. The alignment and SES parameters are again determined through statistical inference, more specifically, by alternating two steps: (1) estimating the SES parameters from a given alignment and (2), with the resulting estimates, refining the alignment. The SES parameters are computed by maximum a posteriori (MAP) estimation (step 1), in analogy to the pairwise case. The alignment (step 2) is solved by linear integer programming. In order to test the robustness and reliability of the proposed N-variate SES method, it is first applied to synthetic data. We show that N-variate SES results in more reliable estimates than bivariate SES. Next N-variate SES is applied to two problems in neuroscience: to quantify the firing reliability of Morris-Lecar neurons and to detect anomalies in EEG synchrony of patients with mild cognitive impairment. Those problems were also considered in Parts I and II, respectively. In both cases, the N-variate SES approach yields a more detailed analysis.