Time-frequency analysis of neuronal populations with instantaneous resolution based on noise-assisted multivariate empirical mode decomposition

Empirical mode decomposition of local field potential data from optogenetic experiments

Introduction

This study investigated the effects of cocaine administration and parvalbumin-type interneuron stimulation on local field potentials (LFPs) recorded in vivo from the medial prefrontal cortex (mPFC) of six mice using optogenetic tools.

Methods

The local network was subject to a brief 10 ms laser pulse, and the response was recorded for 2 s over 100 trials for each of the six subjects who showed stable coupling between the mPFC and the optrode. Due to the strong non-stationary and nonlinearity of the LFP, we used the adaptive, data-driven, Empirical Mode Decomposition (EMD) method to decompose the signal into orthogonal Intrinsic Mode Functions (IMFs).

Results

Through trial and error, we found that seven is the optimum number of orthogonal IMFs that overlaps with known frequency bands of brain activity. We found that the Index of Orthogonality (IO) of IMF amplitudes was close to zero. The Index of Energy Conservation (IEC) for each decomposition was close to unity, as expected for orthogonal decompositions. We found that the power density distribution vs. frequency follows a power law with an average scaling exponent of ~1.4 over the entire range of IMF frequencies 2-2,000 Hz.

Discussion

The scaling exponent is slightly smaller for cocaine than the control, suggesting that neural activity avalanches under cocaine have longer life spans and sizes.

Medium spiny neurons activity reveals the discrete segregation of mouse dorsal striatum

Behavioral studies differentiate the rodent dorsal striatum (DS) into lateral and medial regions; however, anatomical evidence suggests that it is a unified structure. To understand striatal dynamics and basal ganglia functions, it is essential to clarify the circuitry that supports this behavioral-based segregation. Here, we show that the mouse DS is made of two non-overlapping functional circuits divided by a boundary. Combining in vivo optopatch-clamp and extracellular recordings of spontaneous and evoked sensory activity, we demonstrate different coupling of lateral and medial striatum to the cortex together with an independent integration of the spontaneous activity, due to particular corticostriatal connectivity and local attributes of each region. Additionally, we show differences in slow and fast oscillations and in the electrophysiological properties between striatonigral and striatopallidal neurons. In summary, these results demonstrate that the rodent DS is segregated in two neuronal circuits, in homology with the caudate and putamen nuclei of primates.

Comparative study of extracellular recording methods for analysis of afferent sensory information: Empirical modeling, data analysis and interpretation

BACKGROUND

Physiological studies of sensorial systems often require the acquisition and processing of data extracted from their multiple components to evaluate how the neural information changes in relation to the environment changes. In this work, a comparative study about methodological aspects of two electrophysiological approaches is described.

NEW METHOD

Extracellular recordings from deep vibrissal nerves were obtained by using a customized microelectrode Utah array during passive mechanical stimulation of rat´s whiskers. These recordings were compared with those obtained with bipolar electrodes. We also propose here a simplified empirical model of the electrophysiological activity obtained from a bundle of myelinated nerve fibers.

RESULTS

The peripheral activity of the vibrissal system was characterized through the temporal and spectral features obtained with both recording methods. The empirical model not only allows the correlation between anatomical structures and functional features, but also allows to predict changes in the CAPs morphology when the arrangement and the geometry of the electrodes changes.

COMPARISON WITH EXISTING METHOD(S)

This study compares two extracellular recording methods based on analysis techniques, empirical modeling and data processing of vibrissal sensory information.

CONCLUSIONS

This comparative study reveals a close relationship between the electrophysiological techniques and the processing methods necessary to extract sensory information. This relationship is the result of maximizing the extraction of information from recordings of sensory activity.

Unmixing Oscillatory Brain Activity by EEG Source Localization and Empirical Mode Decomposition

Neuronal activity is composed of synchronous and asynchronous oscillatory activity at different frequencies. The neuronal oscillations occur at time scales well matched to the temporal resolution of electroencephalography (EEG); however, to derive meaning from the electrical brain activity as measured from the scalp, it is useful to decompose the EEG signal in space and time. In this study, we elaborate on the investigations into source-based signal decomposition of EEG. Using source localization, the electrical brain signal is spatially unmixed and the neuronal dynamics from a region of interest are analyzed using empirical mode decomposition (EMD), a technique aimed at detecting periodic signals. We demonstrate, first in simulations, that the EMD is more accurate when applied to the spatially unmixed signal compared to the scalp-level signal. Furthermore, on EEG data recorded simultaneously with transcranial magnetic stimulation (TMS) over the hand area of the primary motor cortex, we observe a link between the peak to peak amplitude of the motor-evoked potential (MEP) and the phase of the decomposed localized electrical activity before TMS onset. The results thus encourage combination of source localization and EMD in the pursuit of further insight into the mechanisms of the brain with respect to the phase and frequency of the electrical oscillations and their cortical origin.

A 3D Convolutional Neural Network to Model Retinal Ganglion Cell's Responses to Light Patterns in Mice

Deep Learning offers flexible powerful tools that have advanced our understanding of the neural coding of neurosensory systems. In this work, a 3D Convolutional Neural Network (3D CNN) is used to mimic the behavior of a population of mice retinal ganglion cells in response to different light patterns. For this purpose, we projected homogeneous RGB flashes and checkerboards stimuli with variable luminances and wavelength spectrum to mimic a more naturalistic stimuli environment onto the mouse retina. We also used white moving bars in order to localize the spatial position of the recorded cells. Then recorded spikes were smoothed with a Gaussian kernel and used as the output target when training a 3D CNN in a supervised way. To find a suitable model, two hyperparameter search stages were performed. In the first stage, a trial and error process allowed us to obtain a system that is able to fit the neurons firing rates. In the second stage, a systematic procedure was used to compare several gradient-based optimizers, loss functions and the model's convolutional layers number. We found that a three layered 3D CNN was able to predict the ganglion cells firing rates with high correlations and low prediction error, as measured with Mean Squared Error and Dynamic Time Warping in test sets. These models were either competitive or outperformed other models used already in neuroscience, as Feed Forward Neural Networks and Linear-Nonlinear models. This methodology allowed us to capture the temporal dynamic response patterns in a robust way, even for neurons with high trial-to-trial variable spontaneous firing rates, when providing the peristimulus time histogram as an output to our model.

Multiscale dynamics of interstimulus interval integration in visual cortex

Although the visual cortex receives information at multiple temporal patterns, much of the research in the field has focused only on intervals shorter than 1 second. Consequently, there is almost no information on what happens at longer temporal intervals. We have tried to address this question recording neuronal populations of the primary visual cortex during visual stimulation with repetitive grating stimuli and intervals ranging from 1 to 7 seconds. Our results showed that firing rate and response stability were dependent of interval duration. In addition, there were collective oscillations with different properties in response to changes in intervals duration. These results suggest that visual cortex could encode visual information at several time scales using oscillations at multiple frequencies.

Toward an Improvement of the Analysis of Neural Coding

Machine learning and artificial intelligence have strong roots on principles of neural computation. Some examples are the structure of the first perceptron, inspired in the retina, neuroprosthetics based on ganglion cell recordings or Hopfield networks. In addition, machine learning provides a powerful set of tools to analyze neural data, which has already proved its efficacy in so distant fields of research as speech recognition, behavioral states classification, or LFP recordings. However, despite the huge technological advances in neural data reduction of dimensionality, pattern selection, and clustering during the last years, there has not been a proportional development of the analytical tools used for Time–Frequency (T–F) analysis in neuroscience. Bearing this in mind, we introduce the convenience of using non-linear, non-stationary tools, EMD algorithms in particular, for the transformation of the oscillatory neural data (EEG, EMG, spike oscillations…) into the T–F domain prior to its analysis with machine learning tools. We support that to achieve meaningful conclusions, the transformed data we analyze has to be as faithful as possible to the original recording, so that the transformations forced into the data due to restrictions in the T–F computation are not extended to the results of the machine learning analysis. Moreover, bioinspired computation such as brain–machine interface may be enriched from a more precise definition of neuronal coding where non-linearities of the neuronal dynamics are considered.