Exploring the dynamics of collective synchronizations in large ensembles of brain signals.

Developing methods characterizing the dynamics of synchronization in large ensemble of electromagnetic brain signals has become an important issue. In this article, we review a recently introduced method for analyzing multivariate phase synchronization in brain signals. The approach is based on the equivalence between phase locking and frequency locking in narrow band signals, which allows tracking multivariate phase synchronization in the time-frequency domain as periods of common frequency among multiple channels. The method is illustrated with simulations of multivariate phase dynamics in coupled oscillators and real multichannel electro- and magnetoencephalographic data recorded prior and during epileptic seizures. The reviewed results support the relevance of this method in the context of brain synchronization, in particular to track transient collective dynamics fluctuating in time, frequency and space.     

Functional Brain Connectivity as a New Feature for P300 Speller

The brain is a large-scale complex network often referred to as the “connectome”. Cognitive functions and information processing are mainly based on the interactions between distant brain regions. However, most of the ‘feature extraction’ methods used in the context of Brain Computer Interface (BCI) ignored the possible functional relationships between different signals recorded from distinct brain areas. In this paper, the functional connectivity quantified by the phase locking value (PLV) was introduced to characterize the evoked responses (ERPs) obtained in the case of target and non-targets visual stimuli. We also tested the possibility of using the functional connectivity in the context of ‘P300 speller’. The proposed approach was compared to the well-known methods proposed in the state of the art of “P300 Speller”, mainly the peak picking, the area, time/frequency based features, the xDAWN spatial filtering and the stepwise linear discriminant analysis (SWLDA). The electroencephalographic (EEG) signals recorded from ten subjects were analyzed offline. The results indicated that phase synchrony offers relevant information for the classification in a P300 speller. High synchronization between the brain regions was clearly observed during target trials, although no significant synchronization was detected for a non-target trial. The results showed also that phase synchrony provides higher performance than some existing methods for letter classification in a P300 speller principally when large number of trials is available. Finally, we tested the possible combination of both approaches (classical features and phase synchrony). Our findings showed an overall improvement of the performance of the P300-speller when using Peak picking, the area and frequency based features. Similar performances were obtained compared to xDAWN and SWLDA when using large number of trials.     

Impact of adaptation currents on synchronization of coupled exponential integrate-and-fire neurons

The ability of spiking neurons to synchronize their activity in a network depends on the response behavior of these neurons as quantified by the phase response curve (PRC) and on coupling properties. The PRC characterizes the effects of transient inputs on spike timing and can be measured experimentally. Here we use the adaptive exponential integrate-and-fire (aEIF) neuron model to determine how subthreshold and spike-triggered slow adaptation currents shape the PRC. Based on that, we predict how synchrony and phase locked states of coupled neurons change in presence of synaptic delays and unequal coupling strengths. We find that increased subthreshold adaptation currents cause a transition of the PRC from only phase advances to phase advances and delays in response to excitatory perturbations. Increased spike-triggered adaptation currents on the other hand predominantly skew the PRC to the right. Both adaptation induced changes of the PRC are modulated by spike frequency, being more prominent at lower frequencies. Applying phase reduction theory, we show that subthreshold adaptation stabilizes synchrony for pairs of coupled excitatory neurons, while spike-triggered adaptation causes locking with a small phase difference, as long as synaptic heterogeneities are negligible. For inhibitory pairs synchrony is stable and robust against conduction delays, and adaptation can mediate bistability of in-phase and anti-phase locking. We further demonstrate that stable synchrony and bistable in/anti-phase locking of pairs carry over to synchronization and clustering of larger networks. The effects of adaptation in aEIF neurons on PRCs and network dynamics qualitatively reflect those of biophysical adaptation currents in detailed Hodgkin-Huxley-based neurons, which underscores the utility of the aEIF model for investigating the dynamical behavior of networks. Our results suggest neuronal spike frequency adaptation as a mechanism synchronizing low frequency oscillations in local excitatory networks, but indicate that inhibition rather than excitation generates coherent rhythms at higher frequencies.     

To Perceive or Not Perceive: The Role of Gamma-band Activity in Signaling Object Percepts

The relation of gamma-band synchrony to holistic perception in which concerns the effects of sensory processing, high level perceptual gestalt formation, motor planning and response is still controversial. To provide a more direct link to emergent perceptual states we have used holistic EEG/ERP paradigms where the moment of perceptual “discovery” of a global pattern was variable. Using a rapid visual presentation of short-lived Mooney objects we found an increase of gamma-band activity locked to perceptual events. Additional experiments using dynamic Mooney stimuli showed that gamma activity increases well before the report of an emergent holistic percept. To confirm these findings in a data driven manner we have further used a support vector machine classification approach to distinguish between perceptual vs. non perceptual states, based on time-frequency features. Sensitivity, specificity and accuracy were all above 95%. Modulations in the 30–75 Hz range were larger for perception states. Interestingly, phase synchrony was larger for perception states for high frequency bands. By focusing on global gestalt mechanisms instead of local processing we conclude that gamma-band activity and synchrony provide a signature of holistic perceptual states of variable onset, which are separable from sensory and motor processing.     

Denoising of ECG signals based on noise reduction algorithms in EMD and wavelet domains

This paper presents a new ECG denoising approach based on noise reduction algorithms in empirical mode decomposition (EMD) and discrete wavelet transform (DWT) domains. Unlike the conventional EMD based ECG denoising approaches that neglect a number of initial intrinsic mode functions (IMFs) containing the QRS complex as well as noise, we propose to perform windowing in the EMD domain in order to reduce the noise from the initial IMFs instead of discarding them completely thus preserving the QRS complex and yielding a relatively cleaner ECG signal. The signal thus obtained is transformed in the DWT domain, where an adaptive soft thresholding based noise reduction algorithm is employed considering the advantageous properties of the DWT compared to that of the EMD in preserving the energy in the presence of noise and in reconstructing the original ECG signal with a better time resolution. Extensive simulations are carried out using the MIT-BIH arrythmia database and the performance of the proposed method is evaluated in terms of several standard metrics. The simulation results show that the proposed method is able to reduce noise from the noisy ECG signals more accurately and consistently in comparison to some of the stateof-the-art methods.     

Time-frequency filtering of MEG signals with matching pursuit.

Time-frequency signal analysis based on various decomposition techniques is widely used in biomedical applications. Matching Pursuit is a new adaptive approach for time-frequency decomposition of such biomedical signals. Its advantage is that it creates a concise signal approximation with the help of a small set of Gabor atoms chosen iteratively from a large and redundant set. In this paper, the usage of Matching Pursuit for time-frequency filtering of biomagnetic signals is proposed. The technique was validated on artificial signals and its performance was tested for varying signal-to-noise ratios using both simulated and real MEG somatic evoked magnetic field data.     

Quantifying Neural Oscillatory Synchronization: A Comparison between Spectral Coherence and Phase-Locking Value Approaches

Synchronization or phase-locking between oscillating neuronal groups is considered to be important for coordination of information among cortical networks. Spectral coherence is a commonly used approach to quantify phase locking between neural signals. We systematically explored the validity of spectral coherence measures for quantifying synchronization among neural oscillators. To that aim, we simulated coupled oscillatory signals that exhibited synchronization dynamics using an abstract phase-oscillator model as well as interacting gamma-generating spiking neural networks. We found that, within a large parameter range, the spectral coherence measure deviated substantially from the expected phase-locking. Moreover, spectral coherence did not converge to the expected value with increasing signal-to-noise ratio. We found that spectral coherence particularly failed when oscillators were in the partially (intermittent) synchronized state, which we expect to be the most likely state for neural synchronization. The failure was due to the fast frequency and amplitude changes induced by synchronization forces. We then investigated whether spectral coherence reflected the information flow among networks measured by transfer entropy (TE) of spike trains. We found that spectral coherence failed to robustly reflect changes in synchrony-mediated information flow between neural networks in many instances. As an alternative approach we explored a phase-locking value (PLV) method based on the reconstruction of the instantaneous phase. As one approach for reconstructing instantaneous phase, we used the Hilbert Transform (HT) preceded by Singular Spectrum Decomposition (SSD) of the signal. PLV estimates have broad applicability as they do not rely on stationarity, and, unlike spectral coherence, they enable more accurate estimations of oscillatory synchronization across a wide range of different synchronization regimes, and better tracking of synchronization-mediated information flow among networks.     

Electrical Brain Responses to an Auditory Illusion and the Impact of Musical Expertise

The presentation of two sinusoidal tones, one to each ear, with a slight frequency mismatch yields an auditory illusion of a beating frequency equal to the frequency difference between the two tones; this is known as binaural beat (BB). The effect of brief BB stimulation on scalp EEG is not conclusively demonstrated. Further, no studies have examined the impact of musical training associated with BB stimulation, yet musicians'' brains are often associated with enhanced auditory processing. In this study, we analysed EEG brain responses from two groups, musicians and non-musicians, when stimulated by short presentation (1 min) of binaural beats with beat frequency varying from 1 Hz to 48 Hz. We focused our analysis on alpha and gamma band EEG signals, and they were analysed in terms of spectral power, and functional connectivity as measured by two phase synchrony based measures, phase locking value and phase lag index. Finally, these measures were used to characterize the degree of centrality, segregation and integration of the functional brain network. We found that beat frequencies belonging to alpha band produced the most significant steady-state responses across groups. Further, processing of low frequency (delta, theta, alpha) binaural beats had significant impact on cortical network patterns in the alpha band oscillations. Altogether these results provide a neurophysiological account of cortical responses to BB stimulation at varying frequencies, and demonstrate a modulation of cortico-cortical connectivity in musicians'' brains, and further suggest a kind of neuronal entrainment of a linear and nonlinear relationship to the beating frequencies.     

Using eigenstructure decompositions of time-varying autoregressions in common spatial patterns-based EEG signal classification

Brain–computer interfaces based on common spatial patterns (CSP) depend on the operational frequency bands of the events to be discriminated. This problem has been addressed through sub-band decompositions of the electroencephalographic signals using filter banks, then the performance relies on the number of filters that are stacked and the criteria to select their bandwidths. Here, we propose an alternative approach based on an eigenstructure decomposition of the signals’ time-varying autoregressions (TVAR). The eigen-based decomposition of the TVAR allows for subject-specific estimation of the principal time-varying frequencies, then such principal eigencomponents can be used in the traditional CSP-based classification. We show through a series of numerical experiments that the proposed classification scheme can achieve a performance which is comparable with the one obtained through the filter bank-based approach. However, our method does not rely on a preliminary selection of a frequency band, yet good performance is achieved under realistic conditions (such as reduced number of sensors and small amount of training data) independently of the time interval selected.     

Empirical Mode Decomposition and Wavelet Transform Based ECG Data Compression Scheme

ObjectiveIn health-care systems, compression is an essential tool to solve the storage and transmission problems. In this regard, this paper reports a new electrocardiogram (ECG) data compression scheme which employs sifting function based empirical mode decomposition (EMD) and discrete wavelet transform.MethodEMD based on sifting function is utilized to get the first intrinsic mode function (IMF). After EMD, the first IMF and four significant sifting functions are combined together. This combination is free from many irrelevant components of the signal. Discrete wavelet transform (DWT) with mother wavelet ‘bior4.4’ is applied to this combination. The transform coefficients obtained after DWT are passed through dead-zone quantization. It discards small transform coefficients lying around zero. Further, integer conversion of coefficients and run-length encoding are utilized to achieve a compressed form of ECG data.ResultsCompression performance of the proposed scheme is evaluated using 48 ECG records of the MIT-BIH arrhythmia database. In the comparison of compression results, it is observed that the proposed method exhibits better performance than many recent ECG compressors. A mean opinion score test is also conducted to evaluate the true quality of the reconstructed ECG signals.ConclusionThe proposed scheme offers better compression performance with preserving the key features of the signal very well.     

Reduced Theta-Band Power and Phase Synchrony during Explicit Verbal Memory Tasks in Female,Non-Clinical Individuals with Schizotypal Traits

The study of non-clinical individuals with schizotypal traits has been considered to provide a promising endophenotypic approach to understanding schizophrenia, because schizophrenia is highly heterogeneous, and a number of confounding factors may affect neuropsychological performance. Here, we investigated whether deficits in explicit verbal memory in individuals with schizotypal traits are associated with abnormalities in the local and inter-regional synchrony of brain activity. Memory deficits have been recognized as a core problem in schizophrenia, and previous studies have consistently shown explicit verbal memory impairment in schizophrenic patients. However, the mechanism of this impairment has not been fully revealed. Seventeen individuals with schizotypal traits and 17 age-matched, normal controls participated. Multichannel event-related electroencephalograms (EEGs) were recorded while the subjects performed a continuous recognition task. Event-related spectral perturbations (ERSPs) and inter-regional theta-band phase locking values (TPLVs) were investigated to determine the differences in local and global neural synchrony between the two subject groups. Additionally, the connection patterns of the TPLVs were quantitatively analyzed using graph theory measures. An old/new effect was found in the induced theta-band ERSP in both groups. However, the difference between the old and new was larger in normal controls than in schizotypal trait group. The tendency of elevated old/new effect in normal controls was observed in anterior-posterior theta-band phase synchrony as well. Our results suggest that explicit memory deficits observed in schizophrenia patients can also be found in non-clinical individuals with psychometrically defined schizotypal traits.     

Information quantification of empirical mode decomposition and applications to field potentials

The empirical mode decomposition (EMD) method can adaptively decompose a non-stationary time series into a number of amplitude or frequency modulated functions known as intrinsic mode functions. This paper combines the EMD method with information analysis and presents a framework of information-preserving EMD. The enhanced EMD method has been exploited in the analysis of neural recordings. It decomposes a signal and extracts only the most informative oscillations contained in the non-stationary signal. Information analysis has shown that the extracted components retain the information content of the signal. More importantly, a limited number of components reveal the main oscillations presented in the signal and their instantaneous frequencies, which are not often obvious from the original signal. This information-coupled EMD method has been tested on several field potential datasets for the analysis of stimulus coding in visual cortex, from single and multiple channels, and for finding information connectivity among channels. The results demonstrate the usefulness of the method in extracting relevant responses from the recorded signals. An investigation is also conducted on utilizing the Hilbert phase for cases where phase information can further improve information analysis and stimulus discrimination. The components of the proposed method have been integrated into a toolbox and the initial implementation is also described.     

Empirical mode decomposition based filtering techniques for power line interference reduction in electrocardiogram using various adaptive structures and subtraction methods

Electrocardiogram (ECG) is a vital sign monitoring measurement of the cardiac activity. One of the main problems in biomedical signals like electrocardiogram is the separation of the desired signal from noises caused by power line interference, muscle artifacts, baseline wandering and electrode artifacts. Different types of digital filters are used to separate signal components from unwanted frequency ranges. Adaptive filter is one of the primary methods to filter, because it does not need the signal statistic characteristics. In contrast with Fourier analysis and wavelet methods, a new technique called EMD, a fully data-driven technique is used. It is an adaptive method well suited to analyze biomedical signals. This paper foregrounds an empirical mode decomposition based two-weight adaptive filter structure to eliminate the power line interference in ECG signals. This paper proposes four possible methods and each have less computational complexity compared to other methods. These methods of filtering are fully a signal-dependent approach with adaptive nature, and hence it is best suited for denoising applications. Compared to other proposed methods, EMD based direct subtraction method gives better SNR irrespective of the level of noises.     

癫痫发作间期alpha波的窄带相位同步分析

神经元同步化放电是癫痫发作的一个重要特征,作者提出并运用窄带相位同步技术对比分析了54个癫痫病人和10个正常成年人的脑电信号(EEG)数据.结果表明,相对于对照组,癫痫组的Alpha波的平均窄带相位同步值有显著下降(P=0.02058).为了更准确地刻画和衡量癫痫组和其对照组在同步模式上的差异,提出了一个新的具体的量化指标,即alpha波的窄带相位同步发散率.分析结果显示,对照组的窄带相位同步发散率明显低于癫痫组(P=0.003060),说明对照组的alpha波振子间互诱导强度更高.这可能反映了癫痫组窄带相位同步发散率的升高及所代表的alpha振子间互诱导强度的减弱与患者在癫痫发作间期的状态有密切的联系.     

Partial phase synchronization of neural populations due to random Poisson inputs

We show that populations of identical uncoupled neurons exhibit partial phase synchronization when stimulated with independent, random unidirectional current spikes with interspike time intervals drawn from a Poisson distribution. We characterize this partial synchronization using the phase distribution of the population, and consider analytical approximations and numerical simulations of phase-reduced models and the corresponding conductance-based models of typical Type I (Hindmarsh-Rose) and Type II (Hodgkin-Huxley) neurons, showing quantitatively how the extent of the partial phase synchronization depends on the magnitude and mean interspike frequency of the stimulus. Furthermore, we present several simple examples that disprove the notion that phase synchrony must be strongly related to spike synchrony. Instead, the importance of partial phase synchrony is shown to lie in its influence on the response of the population to stimulation, which we illustrate using first spike time histograms.     

Statistical Detection of EEG Synchrony Using Empirical Bayesian Inference

There is growing interest in understanding how the brain utilizes synchronized oscillatory activity to integrate information across functionally connected regions. Computing phase-locking values (PLV) between EEG signals is a popular method for quantifying such synchronizations and elucidating their role in cognitive tasks. However, high-dimensionality in PLV data incurs a serious multiple testing problem. Standard multiple testing methods in neuroimaging research (e.g., false discovery rate, FDR) suffer severe loss of power, because they fail to exploit complex dependence structure between hypotheses that vary in spectral, temporal and spatial dimension. Previously, we showed that a hierarchical FDR and optimal discovery procedures could be effectively applied for PLV analysis to provide better power than FDR. In this article, we revisit the multiple comparison problem from a new Empirical Bayes perspective and propose the application of the local FDR method (locFDR; Efron, 2001) for PLV synchrony analysis to compute FDR as a posterior probability that an observed statistic belongs to a null hypothesis. We demonstrate the application of Efron''s Empirical Bayes approach for PLV synchrony analysis for the first time. We use simulations to validate the specificity and sensitivity of locFDR and a real EEG dataset from a visual search study for experimental validation. We also compare locFDR with hierarchical FDR and optimal discovery procedures in both simulation and experimental analyses. Our simulation results showed that the locFDR can effectively control false positives without compromising on the power of PLV synchrony inference. Our results from the application locFDR on experiment data detected more significant discoveries than our previously proposed methods whereas the standard FDR method failed to detect any significant discoveries.     

基于EMD的左右手运动想象脑电信号分析研究

基于经验模态分解（EMD）理论,提出一种左右手运动想象脑电信号分析方法。首先利用时间窗对脑电信号数据进行划分,对每段数据通过经验模态分解法将其分解为一组固有模态函数IMF,提取主要信号所在的IMF层去除信号中的噪声。对含有主要信号的几层IMF进行Hilbert变换,得到瞬时频率与对应的瞬时幅值。再提取左右手想象的特定频段mu节律和beta节律的能量信号作为特征,分别利用支持向量机（SVM）和Fisher进行了分类比较。对EMD和小波包在去噪和特征提取进行了比较。结果表明,EMD是一种很有效的去噪方法,经过EMD分解后提取的能量信号在区分左右手想象上更具有优势,识别率高。     

Classification of Motor Imagery EEG Based on Time-Domain and Frequency-Domain Dual-Stream Convolutional Neural Network

Background and objectiveAn important task of the brain-computer interface (BCI) of motor imagery is to extract effective time-domain features, frequency-domain features or time-frequency domain features from the raw electroencephalogram (EEG) signals for classification of motor imagery. However, choosing an appropriate method to combine time domain and frequency domain features to improve the performance of motor imagery recognition is still a research hotspot.MethodsIn order to fully extract and utilize the time-domain and frequency-domain features of EEG in classification tasks, this paper proposed a novel dual-stream convolutional neural network (DCNN), which can use time domain signal and frequency domain signal as the inputs, and the extracted time-domain features and frequency-domain features are fused by linear weighting for classification training. Furthermore, the weight can be learned by the DCNN automatically.ResultsThe experiments based on BCI competition II dataset III and BCI competition IV dataset 2a showed that the model proposed by this study has better performance than other conventional methods. The model used time-frequency signal as the inputs had better performance than the model only used time-domain signals or frequency-domain signals. The accuracy of classification was improved for each subject compared with the models only used one signals as the inputs.ConclusionsFurther analysis shown that the fusion weight of different subject is specifically, adjusting the weight coefficient automatically is helpful to improve the classification accuracy.     

Frequency Dependent Topological Patterns of Resting-State Brain Networks

The topological organization underlying brain networks has been extensively investigated using resting-state fMRI, focusing on the low frequency band from 0.01 to 0.1 Hz. However, the frequency specificities regarding the corresponding brain networks remain largely unclear. In the current study, a data-driven method named complementary ensemble empirical mode decomposition (CEEMD) was introduced to separate the time series of each voxel into several intrinsic oscillation rhythms with distinct frequency bands. Our data indicated that the whole brain BOLD signals could be automatically divided into five specific frequency bands. After applying the CEEMD method, the topological patterns of these five temporally correlated networks were analyzed. The results showed that global topological properties, including the network weighted degree, network efficiency, mean characteristic path length and clustering coefficient, were observed to be most prominent in the ultra-low frequency bands from 0 to 0.015 Hz. Moreover, the saliency of small-world architecture demonstrated frequency-density dependency. Compared to the empirical mode decomposition method (EMD), CEEMD could effectively eliminate the mode-mixing effects. Additionally, the robustness of CEEMD was validated by the similar results derived from a split-half analysis and a conventional frequency division method using the rectangular window band-pass filter. Our findings suggest that CEEMD is a more effective method for extracting the intrinsic oscillation rhythms embedded in the BOLD signals than EMD. The application of CEEMD in fMRI data analysis will provide in-depth insight in investigations of frequency specific topological patterns of the dynamic brain networks.