Open source tools for the information theoretic analysis of neural data

Variable but not random: temporal pattern coding in a songbird brain area necessary for song modification

Practice of a complex motor gesture involves motor exploration to attain a better match to target, but little is known about the neural code for such exploration. We examine spiking in a premotor area of the songbird brain critical for song modification and quantify correlations between spiking and time in the motor sequence. While isolated spikes code for time in song during performance of song to a female bird, extended strings of spiking and silence, particularly bursts, code for time in song during undirected (solo) singing, or "practice." Bursts code for particular times in song with more information than individual spikes, and this spike-spike synergy is significantly higher during undirected singing. The observed pattern information cannot be accounted for by a Poisson model with a matched time-varying rate, indicating that the precise timing of spikes in both bursts in undirected singing and isolated spikes in directed singing code for song with a temporal code. Temporal coding during practice supports the hypothesis that lateral magnocellular nucleus of the anterior nidopallium neurons actively guide song modification at local instances in time.NEW & NOTEWORTHY This paper shows that bursts of spikes in the songbird brain during practice carry information about the output motor pattern. The brain's code for song changes with social context, in performance versus practice. Synergistic combinations of spiking and silence code for time in the bird's song. This is one of the first uses of information theory to quantify neural information about a motor output. This activity may guide changes to the song.

NeuroChaT: A toolbox to analyse the dynamics of neuronal encoding in freely-behaving rodents in vivo

There is a dearth of freely-available, standardised open source analysis tools available for the analysis of neuronal signals recorded in vivo in the freely-behaving animal. In response, we have developed a freely-available, open-source toolbox, NeuroChaT ( Neuron Characterisation Toolbox), specifically addressing this lacuna. Although we have particularly emphasised single unit analyses for spatial coding, NeuroChaT also characterises rhythmic properties of units and their dynamics associated with local field potential signals. NeuroChaT was developed using Python and facilitates a complete pipeline from automation of analysis to producing and managing publication-quality figures. Additionally, we have adopted a platform-independent format (Hierarchical Data Format version 5) for storing recorded and analysed data. By providing an easy-to-use software package, we aim to simplify the adoption of standardised analyses for behavioural neurophysiology and facilitate open data sharing and collaboration between laboratories.

Comparing Open-Source Toolboxes for Processing and Analysis of Spike and Local Field Potentials Data

Analysis of spike and local field potential (LFP) data is an essential part of neuroscientific research. Today there exist many open-source toolboxes for spike and LFP data analysis implementing various functionality. Here we aim to provide a practical guidance for neuroscientists in the choice of an open-source toolbox best satisfying their needs. We overview major open-source toolboxes for spike and LFP data analysis as well as toolboxes with tools for connectivity analysis, dimensionality reduction and generalized linear modeling. We focus on comparing toolboxes functionality, statistical and visualization tools, documentation and support quality. To give a better insight, we compare and illustrate functionality of the toolboxes on open-access dataset or simulated data and make corresponding MATLAB scripts publicly available.

High-density multi-fiber photometry for studying large-scale brain circuit dynamics

Animal behavior originates from neuronal activity distributed across brain-wide networks. However, techniques available to assess large-scale neural dynamics in behaving animals remain limited. Here we present compact, chronically implantable, high-density arrays of optical fibers that enable multi-fiber photometry and optogenetic perturbations across many regions in the mammalian brain. In mice engaged in a texture discrimination task, we achieved simultaneous photometric calcium recordings from networks of 12-48 brain regions, including striatal, thalamic, hippocampal and cortical areas. Furthermore, we optically perturbed subsets of regions in VGAT-ChR2 mice by targeting specific fiber channels with a spatial light modulator. Perturbation of ventral thalamic nuclei caused distributed network modulation and behavioral deficits. Finally, we demonstrate multi-fiber photometry in freely moving animals, including simultaneous recordings from two mice during social interaction. High-density multi-fiber arrays are versatile tools for the investigation of large-scale brain dynamics during behavior.

A Tutorial for Information Theory in Neuroscience

Understanding how neural systems integrate, encode, and compute information is central to understanding brain function. Frequently, data from neuroscience experiments are multivariate, the interactions between the variables are nonlinear, and the landscape of hypothesized or possible interactions between variables is extremely broad. Information theory is well suited to address these types of data, as it possesses multivariate analysis tools, it can be applied to many different types of data, it can capture nonlinear interactions, and it does not require assumptions about the structure of the underlying data (i.e., it is model independent). In this article, we walk through the mathematics of information theory along with common logistical problems associated with data type, data binning, data quantity requirements, bias, and significance testing. Next, we analyze models inspired by canonical neuroscience experiments to improve understanding and demonstrate the strengths of information theory analyses. To facilitate the use of information theory analyses, and an understanding of how these analyses are implemented, we also provide a free MATLAB software package that can be applied to a wide range of data from neuroscience experiments, as well as from other fields of study.

Bilateral Tactile Input Patterns Decoded at Comparable Levels But Different Time Scales in Neocortical Neurons

The presence of contralateral tactile input can profoundly affect ipsilateral tactile perception, and unilateral stroke in somatosensory areas can result in bilateral tactile deficits, suggesting that bilateral tactile integration is an important part of brain function. Although previous studies have shown that bilateral tactile inputs exist and that there are neural interactions between inputs from the two sides, no previous study explored to what extent the local neuronal circuitry processing contains detailed information about the nature of the tactile input from the two sides. To address this question, we used a recently introduced approach to deliver a set of electrical, reproducible, tactile afferent, spatiotemporal activation patterns, which permits a high-resolution analysis of the neuronal decoding capacity, to the skin of the second forepaw digits of the anesthetized male rat. Surprisingly, we found that individual neurons of the primary somatosensory can decode contralateral and ipsilateral input patterns to comparable extents. Although the contralateral input was stronger and more rapidly decoded, given sufficient poststimulus processing time, ipsilateral decoding levels essentially caught up to contralateral levels. Moreover, there was a weak but significant correlation for neurons with high decoding performance for contralateral tactile input to also perform well on decoding ipsilateral input. Our findings shed new light on the brain mechanisms underlying bimanual haptic integration.SIGNIFICANCE STATEMENT Here we demonstrate that the spiking activity of single neocortical neurons in the somatosensory cortex of the rat can be used to decode patterned tactile stimuli delivered to the distal ventral skin of the second forepaw digits on both sides of the body. Even though comparable levels of decoding of the tactile input were achieved faster for contralateral input, given sufficient integration time each neuron was found to decode ipsilateral input with a comparable level of accuracy. Given that the neocortical neurons could decode ipsilateral inputs with such small differences between the patterns suggests that S1 cortex has access to very precise information about ipsilateral events. The findings shed new light on possible network mechanisms underlying bimanual haptic processing.

NeoAnalysis: a Python-based toolbox for quick electrophysiological data processing and analysis

Background

In a typical electrophysiological experiment, especially one that includes studying animal behavior, the data collected normally contain spikes, local field potentials, behavioral responses and other associated data. In order to obtain informative results, the data must be analyzed simultaneously with the experimental settings. However, most open-source toolboxes currently available for data analysis were developed to handle only a portion of the data and did not take into account the sorting of experimental conditions. Additionally, these toolboxes require that the input data be in a specific format, which can be inconvenient to users. Therefore, the development of a highly integrated toolbox that can process multiple types of data regardless of input data format and perform basic analysis for general electrophysiological experiments is incredibly useful.

Results

Here, we report the development of a Python based open-source toolbox, referred to as NeoAnalysis, to be used for quick electrophysiological data processing and analysis. The toolbox can import data from different data acquisition systems regardless of their formats and automatically combine different types of data into a single file with a standardized format. In cases where additional spike sorting is needed, NeoAnalysis provides a module to perform efficient offline sorting with a user-friendly interface. Then, NeoAnalysis can perform regular analog signal processing, spike train, and local field potentials analysis, behavioral response (e.g. saccade) detection and extraction, with several options available for data plotting and statistics. Particularly, it can automatically generate sorted results without requiring users to manually sort data beforehand. In addition, NeoAnalysis can organize all of the relevant data into an informative table on a trial-by-trial basis for data visualization. Finally, NeoAnalysis supports analysis at the population level.

Conclusions

With the multitude of general-purpose functions provided by NeoAnalysis, users can easily obtain publication-quality figures without writing complex codes. NeoAnalysis is a powerful and valuable toolbox for users doing electrophysiological experiments.

Electronic supplementary material

The online version of this article (10.1186/s12938-017-0419-7) contains supplementary material, which is available to authorized users.

PRANAS: A New Platform for Retinal Analysis and Simulation

The retina encodes visual scenes by trains of action potentials that are sent to the brain via the optic nerve. In this paper, we describe a new free access user-end software allowing to better understand this coding. It is called PRANAS (https://pranas.inria.fr), standing for Platform for Retinal ANalysis And Simulation. PRANAS targets neuroscientists and modelers by providing a unique set of retina-related tools. PRANAS integrates a retina simulator allowing large scale simulations while keeping a strong biological plausibility and a toolbox for the analysis of spike train population statistics. The statistical method (entropy maximization under constraints) takes into account both spatial and temporal correlations as constraints, allowing to analyze the effects of memory on statistics. PRANAS also integrates a tool computing and representing in 3D (time-space) receptive fields. All these tools are accessible through a friendly graphical user interface. The most CPU-costly of them have been implemented to run in parallel.

The structure function as new integral measure of spatial and temporal properties of multichannel EEG

The first-order temporal structure functions (SFs), i.e., the first-order statistical moment of absolute increments of scaled multichannel resting state EEG signals in healthy children and teenagers over a wide range of temporal separation (time lags) are computed. Our research shows that the sill level (asymptote) of the SF is mainly defined by a determinant of EEG correlation matrix reflecting the EEG spatial structure. The temporal structure of EEG is found to be characterized by power-law scaling or statistical-scale invariance over time scales less than 0.028 s and at least by two dominant frequencies differing by less than 0.3 Hz. These frequencies define the oscillation behavior of the SF and are mainly distributed within the range of 7.5-12.0 Hz. In this paper, we propose the combined Bessel and exponential model that fits well the empirical SF. It provides a good fit with the mean relative error fitting of 2.8 % over the time lag range of 1 s, using a sampling interval of 4 ms, for all cases under analysis. We also show that the hyper gamma distribution (HGD) fits to the empirical probability density functions (PDFs) of absolute increments of scaled multichannel resting state EEG signals at any given time lag. It means that only two parameters (sample mean of absolute increments and relevant coefficient of variation) may approximately define the empirical PDFs for a given number of channels. A three-dimensional feature vector constructed from the shape and scale parameters of the HGD and the sill level may be used to estimate the closeness of the real EEG to the "random" EEG characterized by the absence of temporal and spatial correlation.

Temporal Code-Driven Stimulation: Definition and Application to Electric Fish Signaling

Closed-loop activity-dependent stimulation is a powerful methodology to assess information processing in biological systems. In this context, the development of novel protocols, their implementation in bioinformatics toolboxes and their application to different description levels open up a wide range of possibilities in the study of biological systems. We developed a methodology for studying biological signals representing them as temporal sequences of binary events. A specific sequence of these events (code) is chosen to deliver a predefined stimulation in a closed-loop manner. The response to this code-driven stimulation can be used to characterize the system. This methodology was implemented in a real time toolbox and tested in the context of electric fish signaling. We show that while there are codes that evoke a response that cannot be distinguished from a control recording without stimulation, other codes evoke a characteristic distinct response. We also compare the code-driven response to open-loop stimulation. The discussed experiments validate the proposed methodology and the software toolbox.

Transition between Functional Regimes in an Integrate-And-Fire Network Model of the Thalamus

The thalamus is a key brain element in the processing of sensory information. During the sleep and awake states, this brain area is characterized by the presence of two distinct dynamical regimes: in the sleep state activity is dominated by spindle oscillations (7 − 15 Hz) weakly affected by external stimuli, while in the awake state the activity is primarily driven by external stimuli. Here we develop a simple and computationally efficient model of the thalamus that exhibits two dynamical regimes with different information-processing capabilities, and study the transition between them. The network model includes glutamatergic thalamocortical (TC) relay neurons and GABAergic reticular (RE) neurons described by adaptive integrate-and-fire models in which spikes are induced by either depolarization or hyperpolarization rebound. We found a range of connectivity conditions under which the thalamic network composed by these neurons displays the two aforementioned dynamical regimes. Our results show that TC-RE loops generate spindle-like oscillations and that a minimum level of clustering (i.e. local connectivity density) in the RE-RE connections is necessary for the coexistence of the two regimes. We also observe that the transition between the two regimes occurs when the external excitatory input on TC neurons (mimicking sensory stimulation) is large enough to cause a significant fraction of them to switch from hyperpolarization-rebound-driven firing to depolarization-driven firing. Overall, our model gives a novel and clear description of the role that the two types of neurons and their connectivity play in the dynamical regimes observed in the thalamus, and in the transition between them. These results pave the way for the development of efficient models of the transmission of sensory information from periphery to cortex.

Dynamics of Population Activity in Rat Sensory Cortex: Network Correlations Predict Anatomical Arrangement and Information Content

To study the spatiotemporal dynamics of neural activity in a cortical population, we implanted a 10 × 10 microelectrode array in the vibrissal cortex of urethane-anesthetized rats. We recorded spontaneous neuronal activity as well as activity evoked in response to sustained and brief sensory stimulation. To quantify the temporal dynamics of activity, we computed the probability distribution function (PDF) of spiking on one electrode given the observation of a spike on another. The spike-triggered PDFs quantified the strength, temporal delay, and temporal precision of correlated activity across electrodes. Nearby cells showed higher levels of correlation at short delays, whereas distant cells showed lower levels of correlation, which tended to occur at longer delays. We found that functional space built based on the strength of pairwise correlations predicted the anatomical arrangement of electrodes. Moreover, the correlation profile of electrode pairs during spontaneous activity predicted the "signal" and "noise" correlations during sensory stimulation. Finally, mutual information analyses revealed that neurons with stronger correlations to the network during spontaneous activity, conveyed higher information about the sensory stimuli in their evoked response. Given the 400-μm-distance between adjacent electrodes, our functional quantifications unravel the spatiotemporal dynamics of activity among nearby and distant cortical columns.

Processing and Analysis of Multichannel Extracellular Neuronal Signals: State-of-the-Art and Challenges

In recent years multichannel neuronal signal acquisition systems have allowed scientists to focus on research questions which were otherwise impossible. They act as a powerful means to study brain (dys)functions in in-vivo and in in-vitro animal models. Typically, each session of electrophysiological experiments with multichannel data acquisition systems generate large amount of raw data. For example, a 128 channel signal acquisition system with 16 bits A/D conversion and 20 kHz sampling rate will generate approximately 17 GB data per hour (uncompressed). This poses an important and challenging problem of inferring conclusions from the large amounts of acquired data. Thus, automated signal processing and analysis tools are becoming a key component in neuroscience research, facilitating extraction of relevant information from neuronal recordings in a reasonable time. The purpose of this review is to introduce the reader to the current state-of-the-art of open-source packages for (semi)automated processing and analysis of multichannel extracellular neuronal signals (i.e., neuronal spikes, local field potentials, electroencephalogram, etc.), and the existing Neuroinformatics infrastructure for tool and data sharing. The review is concluded by pinpointing some major challenges that are being faced, which include the development of novel benchmarking techniques, cloud-based distributed processing and analysis tools, as well as defining novel means to share and standardize data.

A robust and representative lower bound on object processing speed in humans

How early does the brain decode object categories? Addressing this question is critical to constrain the type of neuronal architecture supporting object categorization. In this context, much effort has been devoted to estimating face processing speed. With onsets estimated from 50 to 150 ms, the timing of the first face-sensitive responses in humans remains controversial. This controversy is due partially to the susceptibility of dynamic brain measurements to filtering distortions and analysis issues. Here, using distributions of single-trial event-related potentials (ERPs), causal filtering, statistical analyses at all electrodes and time points, and effective correction for multiple comparisons, we present evidence that the earliest categorical differences start around 90 ms following stimulus presentation. These results were obtained from a representative group of 120 participants, aged 18-81, who categorized images of faces and noise textures. The results were reliable across testing days, as determined by test-retest assessment in 74 of the participants. Furthermore, a control experiment showed similar ERP onsets for contrasts involving images of houses or white noise. Face onsets did not change with age, suggesting that face sensitivity occurs within 100 ms across the adult lifespan. Finally, the simplicity of the face-texture contrast, and the dominant midline distribution of the effects, suggest the face responses were evoked by relatively simple image properties and are not face specific. Our results provide a new lower benchmark for the earliest neuronal responses to complex objects in the human visual system.

Phase-coherence transitions and communication in the gamma range between delay-coupled neuronal populations

Synchronization between neuronal populations plays an important role in information transmission between brain areas. In particular, collective oscillations emerging from the synchronized activity of thousands of neurons can increase the functional connectivity between neural assemblies by coherently coordinating their phases. This synchrony of neuronal activity can take place within a cortical patch or between different cortical regions. While short-range interactions between neurons involve just a few milliseconds, communication through long-range projections between different regions could take up to tens of milliseconds. How these heterogeneous transmission delays affect communication between neuronal populations is not well known. To address this question, we have studied the dynamics of two bidirectionally delayed-coupled neuronal populations using conductance-based spiking models, examining how different synaptic delays give rise to in-phase/anti-phase transitions at particular frequencies within the gamma range, and how this behavior is related to the phase coherence between the two populations at different frequencies. We have used spectral analysis and information theory to quantify the information exchanged between the two networks. For different transmission delays between the two coupled populations, we analyze how the local field potential and multi-unit activity calculated from one population convey information in response to a set of external inputs applied to the other population. The results confirm that zero-lag synchronization maximizes information transmission, although out-of-phase synchronization allows for efficient communication provided the coupling delay, the phase lag between the populations, and the frequency of the oscillations are properly matched.

Estimating Neuronal Information: Logarithmic Binning of Neuronal Inter-Spike Intervals

Neurons communicate via the relative timing of all-or-none biophysical signals called spikes. For statistical analysis, the time between spikes can be accumulated into inter-spike interval histograms. Information theoretic measures have been estimated from these histograms to assess how information varies across organisms, neural systems, and disease conditions. Because neurons are computational units that, to the extent they process time, work not by discrete clock ticks but by the exponential decays of numerous intrinsic variables, we propose that neuronal information measures scale more naturally with the logarithm of time. For the types of inter-spike interval distributions that best describe neuronal activity, the logarithm of time enables fewer bins to capture the salient features of the distributions. Thus, discretizing the logarithm of inter-spike intervals, as compared to the inter-spike intervals themselves, yields histograms that enable more accurate entropy and information estimates for fewer bins and less data. Additionally, as distribution parameters vary, the entropy and information calculated from the logarithm of the inter-spike intervals are substantially better behaved, e.g., entropy is independent of mean rate, and information is equally affected by rate gains and divisions. Thus, when compiling neuronal data for subsequent information analysis, the logarithm of the inter-spike intervals is preferred, over the untransformed inter-spike intervals, because it yields better information estimates and is likely more similar to the construction used by nature herself.