Mahmoud Hassan

EEG source connectivity analysis: from dense array recordings to brain networks.

The recent past years have seen a noticeable increase of interest for electroencephalography (EEG) to analyze functional connectivity through brain sources reconstructed from scalp signals. Although considerable advances have been done both on the recording and analysis of EEG signals, a number of methodological questions are still open regarding the optimal way to process the data in order to identify brain networks. In this paper, we analyze the impact of three factors that intervene in this processing: i) the number of scalp electrodes, ii) the combination between the algorithm used to solve the EEG inverse problem and the algorithm used to measure the functional connectivity and iii) the frequency bands retained to estimate the functional connectivity among neocortical sources. Using High-Resolution (hr) EEG recordings in healthy volunteers, we evaluated these factors on evoked responses during picture recognition and naming task. The main reason for selection this task is that a solid literature background is available about involved brain networks (ground truth). From this a priori information, we propose a performance criterion based on the number of connections identified in the regions of interest (ROI) that belong to potentially activated networks. Our results show that the three studied factors have a dramatic impact on the final result (the identified network in the source space) as strong discrepancies were evidenced depending on the methods used. They also suggest that the combination of weighted Minimum Norm Estimator (wMNE) and the Phase Synchronization (PS) methods applied on High-Resolution EEG in beta/gamma bands provides the best performance in term of topological distance between the identified network and the expected network in the above-mentioned cognitive task.

Combination of Canonical Correlation Analysis and Empirical Mode Decomposition Applied to Denoising the Labor Electrohysterogram

The electrohysterogram (EHG) is often corrupted by electronic and electromagnetic noise as well as movement artifacts, skeletal electromyogram, and ECGs from both mother and fetus. The interfering signals are sporadic and/or have spectra overlapping the spectra of the signals of interest rendering classical filtering ineffective. In the absence of efficient methods for denoising the monopolar EHG signal, bipolar methods are usually used. In this paper, we propose a novel combination of blind source separation using canonical correlation analysis (BSS_CCA) and empirical mode decomposition (EMD) methods to denoise monopolar EHG. We first extract the uterine bursts by using BSS_CCA then the biggest part of any residual noise is removed from the bursts by EMD. Our algorithm, called CCA_EMD, was compared with wavelet filtering and independent component analysis. We also compared CCA_EMD with the corresponding bipolar signals to demonstrate that the new method gives signals that have not been degraded by the new method. The proposed method successfully removed artifacts from the signal without altering the underlying uterine activity as observed by bipolar methods. The CCA_EMD algorithm performed considerably better than the comparison methods.

A new algorithm for spatiotemporal analysis of brain functional connectivity.

Specific networks of interacting neuronal assemblies distributed within and across distinct brain regions underlie brain functions. In most cognitive tasks, these interactions are dynamic and take place at the millisecond time scale. Among neuroimaging techniques, magneto/electroencephalography - M/EEG - allows for detection of very short-duration events and offers the single opportunity to follow, in time, the dynamic properties of cognitive processes (sub-millisecond temporal resolution). In this paper, we propose a new algorithm to track the functional brain connectivity dynamics. During a picture naming task, this algorithm aims at segmenting high-resolution EEG signals (hr-EEG) into functional connectivity microstates. The proposed algorithm is based on the K-means clustering of the connectivity graphs obtained from the phase locking value (PLV) method applied on hr-EEG. Results show that the analyzed evoked responses can be divided into six clusters representing distinct networks sequentially involved during the cognitive task, from the picture presentation and recognition to the motor response.

Dynamic reorganization of functional brain networks during picture naming.

For efficient information processing during cognitive activity, functional brain networks have to rapidly and dynamically reorganize on a sub-second time scale. Tracking the spatiotemporal dynamics of large scale networks over this short time duration is a very challenging issue. Here, we tackle this problem by using dense electroencephalography (EEG) recorded during a picture naming task. We found that (i) the picture naming task can be divided into six brain network states (BNSs) characterized by significantly high synchronization of gamma (30-45 Hz) oscillations, (ii) fast transitions occur between these BNSs that last from 30 msec to 160 msec, (iii) based on the state of the art of the picture naming task, we consider that the spatial location of their nodes and edges, as well as the timing of transitions, indicate that each network can be associated with one or several specific function (from visual processing to articulation) and (iv) the comparison with previously-used approach aimed at localizing the sources showed that the network-based approach reveals networks that are more specific to the performed task. We speculate that the persistence of several brain regions in successive BNSs participates to fast and efficient information processing in the brain.

Identification of Interictal Epileptic Networks from Dense-EEG.

Epilepsy is a network disease. The epileptic network usually involves spatially distributed brain regions. In this context, noninvasive M/EEG source connectivity is an emerging technique to identify functional brain networks at cortical level from noninvasive recordings. In this paper, we analyze the effect of the two key factors involved in EEG source connectivity processing: (i) the algorithm used in the solution of the EEG inverse problem and (ii) the method used in the estimation of the functional connectivity. We evaluate four inverse solutions algorithms (dSPM, wMNE, sLORETA and cMEM) and four connectivity measures (r2, h2, PLV, and MI) on data simulated from a combined biophysical/physiological model to generate realistic interictal epileptic spikes reflected in scalp EEG. We use a new network-based similarity index to compare between the network identified by each of the inverse/connectivity combination and the original network generated in the model. The method will be also applied on real data recorded from one epileptic patient who underwent a full presurgical evaluation for drug-resistant focal epilepsy. In simulated data, results revealed that the selection of the inverse/connectivity combination has a significant impact on the identified networks. Results suggested that nonlinear methods (nonlinear correlation coefficient, phase synchronization and mutual information) for measuring the connectivity are more efficient than the linear one (the cross correlation coefficient). The wMNE inverse solution showed higher performance than dSPM, cMEM and sLORETA. In real data, the combination (wMNE/PLV) led to a very good matching between the interictal epileptic network identified from noninvasive EEG recordings and the network obtained from connectivity analysis of intracerebral EEG recordings. These results suggest that source connectivity method, when appropriately configured, is able to extract highly relevant diagnostic information about networks involved in interictal epileptic spikes from non-invasive dense-EEG data.

EEGNET: An Open Source Tool for Analyzing and Visualizing M/EEG Connectome.

The brain is a large-scale complex network often referred to as the “connectome”. Exploring the dynamic behavior of the connectome is a challenging issue as both excellent time and space resolution is required. In this context Magneto/Electroencephalography (M/EEG) are effective neuroimaging techniques allowing for analysis of the dynamics of functional brain networks at scalp level and/or at reconstructed sources. However, a tool that can cover all the processing steps of identifying brain networks from M/EEG data is still missing. In this paper, we report a novel software package, called EEGNET, running under MATLAB (Math works, inc), and allowing for analysis and visualization of functional brain networks from M/EEG recordings. EEGNET is developed to analyze networks either at the level of scalp electrodes or at the level of reconstructed cortical sources. It includes i) Basic steps in preprocessing M/EEG signals, ii) the solution of the inverse problem to localize / reconstruct the cortical sources, iii) the computation of functional connectivity among signals collected at surface electrodes or/and time courses of reconstructed sources and iv) the computation of the network measures based on graph theory analysis. EEGNET is the unique tool that combines the M/EEG functional connectivity analysis and the computation of network measures derived from the graph theory. The first version of EEGNET is easy to use, flexible and user friendly. EEGNET is an open source tool and can be freely downloaded from this webpage: https://sites.google.com/site/eegnetworks/.

Performance analysis of four nonlinearity analysis methods using a model with variable complexity and application to uterine EMG signals

Several measures have been proposed to detect nonlinear characteristics in time series. Results on time series, multiple surrogates and their z-score are used to statistically test for the presence or absence of non-linearity. The z-score itself has sometimes been used as a measure of nonlinearity. The sensitivity of nonlinear methods to the nonlinearity level and their robustness to noise have rarely been evaluated in the past. While surrogates are important tools to rigorously detect nonlinearity, their usefulness for evaluating the level of nonlinearity is not clear. In this paper we investigate the performance of four methods arising from three families that are widely used in non-linearity detection: statistics (time reversibility), predictability (sample entropy, delay vector variance) and chaos theory (Lyapunov exponents). We used sensitivity to increasing complexity and the mean square error (MSE) of Monte Carlo instances for quantitative comparison of their performances. These methods were applied to a Henon nonlinear synthetic model in which we can vary the complexity degree (CD). This was done first by applying the methods directly to the signal and then using the z-score (surrogates) with and without added noise. The methods were then applied to real uterine EMG signals and used to distinguish between pregnancy and labor contraction bursts. The discrimination performances were compared to linear frequency based methods classically used for the same purpose such as mean power frequency (MPF), peak frequency (PF) and median frequency (MF). The results show noticeable difference between different methods, with a clear superiority of some of the nonlinear methods (time reversibility, Lyapunov exponents) over the linear methods. Applying the methods directly to the signals gave better results than using the z-score, except for sample entropy.

Nonlinear estimation of coupling and directionality between signals: Application to uterine EMG propagation

Understanding the direction and quantity of information flowing in a complex system is a fundamental task in signal processing. Several measures have been proposed to detect the quantity of synchronization and the directionality between time series and in physiological data. In this paper we use two methods that are widely used in synchronization and directionality analysis: Nonlinear correlation coefficient (h2) and the general synchronization (H). The performances of both methods were tested on four dimensional coupled synthetic nonlinear Rössler models. They were then applied to a single real labor contraction uterine EMG burst with the aim of using them to detect synchronization and to plot the map of direction of information flow between the whole signal channels. The results on synthetic signal show a slight superiority of H over h2. The results obtained on a single contraction are encouraging for the future use of these tools for resolving the open question of the directionality of uterine contractions and may provide a way of finding their source loci.

The dynamic functional core network of the human brain at rest

The human brain is an inherently complex and dynamic system. Even at rest, functional brain networks dynamically reconfigure in a well-organized way to warrant an efficient communication between brain regions. However, a precise characterization of this reconfiguration at very fast time-scale (hundreds of millisecond) during rest remains elusive. In this study, we used dense electroencephalography data recorded during task-free paradigm to track the fast temporal dynamics of spontaneous brain networks. Results obtained from network-based analysis methods revealed the existence of a functional dynamic core network formed of a set of key brain regions that ensure segregation and integration functions. Brain regions within this functional core share high betweenness centrality, strength and vulnerability (high impact on the network global efficiency) and low clustering coefficient. These regions are mainly located in the cingulate and the medial frontal cortex. In particular, most of the identified hubs were found to belong to the Default Mode Network. Results also revealed that the same central regions may dynamically alternate and play the role of either provincial (local) or connector (global) hubs.

Spatiotemporal Analysis of Brain Functional Connectivity

Brain functions are based on interactions between neural assemblies distributed within and across distinct cerebral regions. During cognitive tasks, these interactions are dynamic and take place at the millisecond time scale. In this context, the excellent temporal resolution (<1 ms) of the Electroencephalographic –EEG- signals allows for detection of very short-duration events and therefore, offers the unique opportunity to follow, over time, the dynamic properties of cognitive processes.