Multivariate Autoregressive Model Constrained by Anatomical Connectivity to Reconstruct Focal Sources

International audienceIn this paper, we present a framework to reconstruct spatially localized sources from Magnetoencephalogra-phy (MEG)/Electroencephalography (EEG) using spatiotempo-ral constraint. The source dynamics are represented by a Mul-tivariate Autoregressive (MAR) model whose matrix elements are constrained by the anatomical connectivity obtained from diffusion Magnetic Resonance Imaging (dMRI). The framework assumes that the whole brain dynamic follows a constant MAR model in a time window of interest. The source activations and the MAR model parameters are estimated iteratively. We could confirm the accuracy of the framework using simulation experiments in both high and low noise levels. The proposed framework outperforms the two-stage approach

Estimating EEG Source Dipole Orientation Based on Singular-value Decomposition for Connectivity Analysis.

In the last decade, the use of high-density electrode arrays for EEG recordings combined with the improvements of source reconstruction algorithms has allowed the investigation of brain networks dynamics at a sub-second scale. One powerful tool for investigating large-scale functional brain networks with EEG is time-varying effective connectivity applied to source signals obtained from electric source imaging. Due to computational and interpretation limitations, the brain is usually parcelled into a limited number of regions of interests (ROIs) before computing EEG connectivity. One specific need and still open problem is how to represent the time- and frequency-content carried by hundreds of dipoles with diverging orientation in each ROI with one unique representative time-series. The main aim of this paper is to provide a method to compute a signal that explains most of the variability of the data contained in each ROI before computing, for instance, time-varying connectivity. As the representative time-series for a ROI, we propose to use the first singular vector computed by a singular-value decomposition of all dipoles belonging to the same ROI. We applied this method to two real datasets (visual evoked potentials and epileptic spikes) and evaluated the time-course and the frequency content of the obtained signals. For each ROI, both the time-course and the frequency content of the proposed method reflected the expected time-course and the scalp-EEG frequency content, representing most of the variability of the sources (~ 80%) and improving connectivity results in comparison to other procedures used so far. We also confirm these results in a simulated dataset with a known ground truth

Large brain effective network from EEG/MEG data and dMR information

International audienceOver the past 30 years, neuroimaging has become a predominant technique. One might envision that over the next years it will play a major role in disclosing the brain's functional interactions. In this work, we use information coming from diffusion magnetic resonance imaging (dMRI) to reconstruct effective brain network from two functional modalities: electroencephalography (EEG) and magnetoen-cephalography (MEG)

EEG-Based Quantification of Cortical Current Density and Dynamic Causal Connectivity Generalized across Subjects Performing BCI-Monitored Cognitive Tasks.

Quantification of dynamic causal interactions among brain regions constitutes an important component of conducting research and developing applications in experimental and translational neuroscience. Furthermore, cortical networks with dynamic causal connectivity in brain-computer interface (BCI) applications offer a more comprehensive view of brain states implicated in behavior than do individual brain regions. However, models of cortical network dynamics are difficult to generalize across subjects because current electroencephalography (EEG) signal analysis techniques are limited in their ability to reliably localize sources across subjects. We propose an algorithmic and computational framework for identifying cortical networks across subjects in which dynamic causal connectivity is modeled among user-selected cortical regions of interest (ROIs). We demonstrate the strength of the proposed framework using a "reach/saccade to spatial target" cognitive task performed by 10 right-handed individuals. Modeling of causal cortical interactions was accomplished through measurement of cortical activity using (EEG), application of independent component clustering to identify cortical ROIs as network nodes, estimation of cortical current density using cortically constrained low resolution electromagnetic brain tomography (cLORETA), multivariate autoregressive (MVAR) modeling of representative cortical activity signals from each ROI, and quantification of the dynamic causal interaction among the identified ROIs using the Short-time direct Directed Transfer function (SdDTF). The resulting cortical network and the computed causal dynamics among its nodes exhibited physiologically plausible behavior, consistent with past results reported in the literature. This physiological plausibility of the results strengthens the framework's applicability in reliably capturing complex brain functionality, which is required by applications, such as diagnostics and BCI

Structural connectivity to reconstruct brain activation and effective connectivity between brain regions

International audienceUnderstanding how brain regions interact to perform a specific task is very challenging. EEG and MEG are two non-invasive imaging modalities that allow the measurement of brain activation with high temporal resolution. Several works in EEG/MEG source reconstruction show that estimating brain activation can be improved by considering spatio-temporal constraints but only few of them use structural information to do so. In this work, we present a source reconstruction algorithm that uses brain structural connectivity, estimated from diffusion MRI (dMRI), to constrain the EEG/MEG source reconstruction. Contrarily to most source reconstruction methods which reconstruct activation for each time instant, the proposed method estimates an initial reconstruction for the first time instants and a multivariate autoregressive model that explains the data in further time instants. This au-toregressive model can be thought as an estimation of the effective connectivity between brain regions. We called this algorithm iterative Source and Dynamics reconstruction (iSDR). This paper presents the overall iSDR approach and how the proposed model is optimized to obtain both brain activation and brain region interactions. The accuracy of our method is demonstrated using synthetic data in which it shows a good capability to reconstruct both activation and connectivity. iSDR is also tested with real data obtained from [dataset] (face recognition task). The results are in phase with other works published with the same data and others that used different imaging modalities with the same task showing that the choice of using an autoregressive model gives relevant results

Influence of the head model on EEG and MEG source connectivity analysis

The results of brain connectivity analysis using reconstructed source time courses derived from EEG and MEG data depend on a number of algorithmic choices. While previous studies have investigated the influence of the choice of source estimation method or connectivity measure, the effects of the head modeling errors or simplifications have not been studied sufficiently. In the present simulation study, we investigated the influence of particular properties of the head model on the reconstructed source time courses as well as on source connectivity analysis in EEG and MEG. Therefore, we constructed a realistic head model and applied the finite element method to solve the EEG and MEG forward problem. We considered the distinction between white and gray matter, the distinction between compact and spongy bone, the inclusion of a cerebrospinal fluid (CSF) compartment, and the reduction to a simple 3-layer model comprising only skin, skull, and brain. Source time courses were reconstructed using a beamforming approach and the source connectivity was estimated by the imaginary coherence (ICoh) and the generalized partial directed coherence (GPDC). Our results show that in both EEG and MEG, neglecting the white and gray matter distinction or the CSF causes considerable errors in reconstructed source time courses and connectivity analysis, while the distinction between spongy and compact bone is just of minor relevance, provided that an adequate skull conductivity value is used. Large inverse and connectivity errors are found in the same regions that show large topography errors in the forward solution. Moreover, we demonstrate that the very conservative ICoh is relatively safe from the crosstalk effects caused by imperfect head models, as opposed to the GPDC

Iterative two-stage approach to estimate sources and their interactions

International audienceNon-iterative two-stage approaches have been used to estimate source interactions. They first reconstruct sources and then compute the MAR model for the localized sources. They showed good results when working in high signal-to-noise ratio (SNR) settings, but fail in detecting the true interactions when working in low SNR. Our framework is based on two steps. First, we estimate sources activations for a given MAR model. Then, we estimate the MAR model. We repeat the two steps until a stopping criterion is achieved

Quantifying the effect of demixing approaches on directed connectivity estimated between reconstructed EEG sources

Electrical activity recorded on the scalp using electroencephalography (EEG) results from the mixing of signals originating from different regions of the brain as well as from artifactual sources. In order to investigate the role of distinct brain areas in a given experiment, the signal recorded on the sensors is typically projected back into the brain (source reconstruction) using algorithms that address the so-called EEG inverse problem. Once the activity of sources located inside of the brain has been reconstructed, it is often desirable to study the statistical dependencies among them, in particular to quantify directional dynamical interactions between brain areas. Unfortunately, even when performing source reconstruction, the superposition of signals that is due to the propagation of activity from sources to sensors cannot be completely undone, resulting in potentially biased estimates of directional functional connectivity. Here we perform a set of simulations involving interacting sources to quantify source connectivity estimation performance as a function of the location of the sources, their distance to each other, the noise level, the source reconstruction algorithm, and the connectivity estimator. The generated source activity was projected onto the scalp and projected back to the cortical level using two source reconstruction algorithms, linearly constrained minimum variance beamforming and Exact' low-resolution tomography (eLORETA). In source space, directed connectivity was estimated using multi-variate Granger causality and time-reversed Granger causality, and compared with the imposed ground truth. Our results demonstrate that all considered factors significantly affect the connectivity estimation performance