The correlation grid: analysis of synchronous spiking in multi-dimensional spike train data and identification of feasible connection architectures

Coupling feature extraction method of resting state EEG Signals from amnestic mild cognitive impairment with type 2 diabetes mellitus based on weight permutation conditional mutual information

This study aimed to find a good coupling feature extraction method to effectively analyze resting state EEG signals (rsEEG) of amnestic mild cognitive impairment(aMCI) with type 2 diabetes mellitus(T2DM) and normal control (NC) with T2DM. A method of EEG signal coupling feature extraction based on weight permutation conditional mutual information (WPCMI) was proposed in this research. With the WPCMI method, coupling feature strength of two time series in Alpha1, Alpha2, Beta1, Beta2 and Gamma bands for aMCI with T2DM and NC with T2DM could be extracted respectively. Then selected three frequency bands coupling feature matrix with the help of multi-spectral image transformation method to map it as spectral image characteristics. And finally classified these characteristics through the convolution neural network method(CNN). For aMCI with T2DM and NC with T2DM, the highest classification accuracy of 96%, 95%, 95% could be achieved respectively in the combination of three frequency bands (Alpha1, Alpha2, Gamma), (Beta1, Beta2 and Gamma) and (Alpha2, Beta1, Beta2). This WPCMI method highlighted the coupling dynamic characteristics of EEG signals, and its classification performance was better than all previous methods in aMCI with T2DM diagnosis field. WPCMI method could be used as an effective biomarker to distinguish EEG signals of aMCI with T2DM and NC with T2DM. The coupling feature extraction method used in this paper provided a new perspective for the EEG analysis of aMCI with T2DM.

Estimating coupling strength between multivariate neural series with multivariate permutation conditional mutual information

Recently, coupling between groups of neurons or different brain regions has been widely studied to provide insights into underlying mechanisms of brain functions. To comprehensively understand the effect of such coupling, it is necessary to accurately extract the coupling strength information among multivariate neural signals from the whole brain. This study proposed a new method named multivariate permutation conditional mutual information (MPCMI) to quantitatively estimate the coupling strength of multivariate neural signals (MNS). The performance of the MPCMI method was validated on the simulated MNS generated by multi-channel neural mass model (MNMM). The coupling strength feature of simulated MNS extracted by MPCMI showed better performance compared with standard methods, such as permutation conditional mutual information (PCMI), multivariate Granger causality (MVGC), and Granger causality analysis (GCA). Furthermore, the MPCMI was applied to estimate the coupling strengths of two-channel resting-state electroencephalographic (rsEEG) signals from different brain regions of 19 patients with amnestic mild cognitive impairment (aMCI) with type 2 diabetes mellitus (T2DM) and 20 normal control (NC) with T2DM in Alpha1 and Alpha2 frequency bands. Empirical results showed that the MPCMI could effectively extract the coupling strength features that were significantly different between the aMCI and the NC. Hence, the proposed MPCMI method could be an effective estimate of coupling strengths of MNS, and might be a viable biomarker for clinical applications.

A global coupling index of multivariate neural series with application to the evaluation of mild cognitive impairment

Recently, the synchronization between neural signals has been widely used as a key indicator of brain function. To understand comprehensively the effect of synchronization on the brain function, accurate computation of the synchronization strength among multivariate neural series from the whole brain is necessary. In this study, we proposed a method named global coupling index (GCI) to estimate the synchronization strength of multiple neural signals. First of all, performance of the GCI method was evaluated by analyzing simulated EEG signals from a multi-channel neural mass model, including the effects of the frequency band, the coupling coefficient, and the signal noise ratio. Then, the GCI method was applied to analyze the EEG signals from 12 mild cognitive impairment (MCI) subjects and 12 normal controls (NC). The results showed that GCI method had two major advantages over the global synchronization index (GSI) or S-estimator. Firstly, simulation data showed that the GCI method provided both a more robust result on the frequency band and a better performance on the coupling coefficients. Secondly, the actual EEG data demonstrated that GCI method was more sensitive in differentiating the MCI from control subjects, in terms of the global synchronization strength of neural series of specific alpha, beta1 and beta2 frequency bands. Hence, it is suggested that GCI is a better method over GSI and S-estimator to estimate the synchronization strength of multivariate neural series for predicting the MCI from the whole brain EEG recordings.

Statistical technique for analysing functional connectivity of multiple spike trains

A new statistical technique, the Cox method, used for analysing functional connectivity of simultaneously recorded multiple spike trains is presented. This method is based on the theory of modulated renewal processes and it estimates a vector of influence strengths from multiple spike trains (called reference trains) to the selected (target) spike train. Selecting another target spike train and repeating the calculation of the influence strengths from the reference spike trains enables researchers to find all functional connections among multiple spike trains. In order to study functional connectivity an "influence function" is identified. This function recognises the specificity of neuronal interactions and reflects the dynamics of postsynaptic potential. In comparison to existing techniques, the Cox method has the following advantages: it does not use bins (binless method); it is applicable to cases where the sample size is small; it is sufficiently sensitive such that it estimates weak influences; it supports the simultaneous analysis of multiple influences; it is able to identify a correct connectivity scheme in difficult cases of "common source" or "indirect" connectivity. The Cox method has been thoroughly tested using multiple sets of data generated by the neural network model of the leaky integrate and fire neurons with a prescribed architecture of connections. The results suggest that this method is highly successful for analysing functional connectivity of simultaneously recorded multiple spike trains.

iRaster: a novel information visualization tool to explore spatiotemporal patterns in multiple spike trains

Over the last few years, simultaneous recordings of multiple spike trains have become widely used by neuroscientists. Therefore, it is important to develop new tools for analysing multiple spike trains in order to gain new insight into the function of neural systems. This paper describes how techniques from the field of visual analytics can be used to reveal specific patterns of neural activity. An interactive raster plot called iRaster has been developed. This software incorporates a selection of statistical procedures for visualization and flexible manipulations with multiple spike trains. For example, there are several procedures for the re-ordering of spike trains which can be used to unmask activity propagation, spiking synchronization, and many other important features of multiple spike train activity. Additionally, iRaster includes a rate representation of neural activity, a combined representation of rate and spikes, spike train removal and time interval removal. Furthermore, it provides multiple coordinated views, time and spike train zooming windows, a fisheye lens distortion, and dissemination facilities. iRaster is a user friendly, interactive, flexible tool which supports a broad range of visual representations. This tool has been successfully used to analyse both synthetic and experimentally recorded datasets. In this paper, the main features of iRaster are described and its performance and effectiveness are demonstrated using various types of data including experimental multi-electrode array recordings from the ganglion cell layer in mouse retina. iRaster is part of an ongoing research project called VISA (Visualization of Inter-Spike Associations) at the Visualization Lab in the University of Plymouth. The overall aim of the VISA project is to provide neuroscientists with the ability to freely explore and analyse their data. The software is freely available from the Visualization Lab website (see www.plymouth.ac.uk/infovis).

Scale-free topology of the CA3 hippocampal network: a novel method to analyze functional neuronal assemblies

Cognitive mapping functions of the hippocampus critically depend on the recurrent network of the CA3 pyramidal cells. However, it is still not known in detail how network activity patterns emerge, or how they encode information. By using functional multineuron calcium imaging, we simultaneously recorded the activity of >100 neurons in the CA3 region of hippocampal slice cultures. We utilized a novel computational method to analyze the multichannel spike trains and to depict functional neuronal assemblies. By means of event synchronization and the correlation matrix analysis method, we found that: 1), the average functional neuronal cluster consists of 23 neurons, and neurons could be part of multiple assemblies; 2), the clustering strength, size, and mean distance among cells in neuronal assemblies follow a power-law-like distribution; 3), the clustering strength and size of neuronal assemblies are not correlated with the total number of neurons and their physical distance; and 4), the clustering distance of neuronal assemblies is weakly correlated with the total number of neurons and their physical distance. These findings suggest that the functional organization of the spontaneously firing CA3 hippocampal network is a scale-free structure in slice culture.

Global synchronization index as a biological correlate of cognitive decline in Alzheimer's disease

OBJECTIVE

The recently developed global synchronization index (GSI) quantifies synchronization between neuronal signals at multiple sites. This study explored the clinical significance of the GSI in Alzheimer's disease (AD) patients.

METHODS

Electroencephalograms were recorded from 25 AD patients and 22 age-matched healthy normal controls (NC). GSI values were computed both across the entire frequency band and separately in the delta, theta, alpha, beta1, beta2, beta3, and gamma bands. The Mini-Mental Status Examination (MMSE) and Clinical Dementia Rating scale (CDR) were used to assess the symptom severity.

RESULTS

GSI values in the beta1, beta2, beta3, and gamma bands were significantly lower in AD patients than in NC. GSI values in the beta and gamma bands were positively correlated with the MMSE scores in all participants (AD and NC). In AD patients, GSI values were negatively correlated with MMSE scores in the delta bands, but positively correlated in the beta1 and gamma band. Also, GSI values were positively correlated with CDR scores in the delta bands, but negatively correlated in the gamma band.

CONCLUSIONS

GSI values of mainly high-frequency bands were significantly lower in AD patients than in NC, they were significantly correlated with scores on symptom severity scales.

SIGNIFICANCE

Our results suggest that GSI values are a useful biological correlate of cognitive decline in AD patients.

A color-based visualization technique for multielectrode spike trains

Multi electrode recordings of neuronal activity provide an overwhelming amount of data that is often difficult to analyze and interpret. Although various methods exist for treating multielectrode datasets quantitatively, there is a particularly prominent lack of techniques that enable a quick visual exploration of such datasets. Here, by using Kohonen self-organizing maps, we propose a simple technique that allows for the representation of multiple spike trains through a sequence of color-coded population activity vectors. When multiple color sequences are grouped according to a certain criterion, e.g., by stimulation condition or recording time, one can inspect an entire dataset visually and extract quickly information about the identity, stimulus-locking and temporal distribution of multi-neuron activity patterns. Color sequences can be computed on various time scales revealing different aspects of the temporal dynamics and can emphasize high-order correlation patterns that are not detectable with pairwise techniques. Furthermore, this technique is useful for determining the stability of neuronal responses during a recording session. Due to its simplicity and reliance on perceptual grouping, the method is useful for both quick on-line visualization of incoming data and for more detailed post hoc analyses.

Visual perception of ambiguous figures: synchronization based neural models

We develop and study two neural network models of perceptual alternations. Both models have a star-like architecture of connections with a central element connected to a set of peripheral elements. A particular perception is simulated in terms of partial synchronization between the central element and some sub-group of peripheral elements. The first model is constructed from phase oscillators and the mechanism of perceptual alternations is based on chaotic intermittency under fixed parameter values. Similar to experimental evidence, the distribution of times between perceptual alternations is represented by the gamma distribution. The second model is built of spiking neurons of the Hodgkin-Huxley type. The mechanism of perceptual alternations is based on plasticity of inhibitory synapses which increases the inhibition from the central unit to the neural assembly representing the current percept. As a result another perception is formed. Simulations show that the second model is in good agreement with behavioural data on switching times between percepts of ambiguous figures and with experimental results on binocular rivalry of two and four percepts.

A stimulus-dependent connectivity analysis of neuronal networks

We present an analysis of interactions among neurons in stimulus-driven networks that is designed to control for effects from unmeasured neurons. This work builds on previous connectivity analyses that assumed connectivity strength to be constant with respect to the stimulus. Since unmeasured neuron activity can modulate with the stimulus, the effective strength of common input connections from such hidden neurons can also modulate with the stimulus. By explicitly accounting for the resulting stimulus-dependence of effective interactions among measured neurons, we are able to remove ambiguity in the classification of causal interactions that resulted from classification errors in the previous analyses. In this way, we can more reliably distinguish causal connections among measured neurons from common input connections that arise from hidden network nodes. The approach is derived in a general mathematical framework that can be applied to other types of networks. We illustrate the effects of stimulus-dependent connectivity estimates with simulations of neurons responding to a visual stimulus.

Synchronization Measurement of Multiple Neuronal Populations

The purpose of the present paper is to develop a method, based on equal-time correlation, correlation matrix analysis and surrogate resampling, that is able to quantify and describe properties of synchronization of population neuronal activity recorded simultaneously from multiple sites. Initially, Lorenz-type oscillators were used to model multiple time series with different patterns of synchronization. Eigenvalue and eigenvector decomposition was then applied to identify “clusters” of locally synchronized activity and to calculate a “global synchronization index.” This method was then applied to multichannel data recorded from an in vitro model of epileptic seizures. The results demonstrate that this novel method can be successfully used to analyze synchronization between multiple neuronal population series.