From EEG to BOLD: Brain mapping and estimating transfer functions in simultaneous EEG-fMRI acquisitions

Simultaneous acquisition of electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) aims to disentangle the description of brain processes by exploiting the advantages of each technique. Most studies in this field focus on exploring the relationships between fMRI signals and the power spectrum at some specific frequency bands (alpha, beta, etc.). On the other hand, brain mapping of EEG signals (e.g., interictal spikes in epileptic patients) usually assumes an haemodynamic response function for a parametric analysis applying the GLM, as a rough approximation. The integration of the information provided by the high spatial resolution of MR images and the high temporal resolution of EEG may be improved by referencing them by transfer functions, which allows the identification of neural driven areas without strong assumptions about haemodynamic response shapes or brain haemodynamic`s homogeneity. The difference on sampling rate is the first obstacle for a full integration of EEG and fMRI information. Moreover, a parametric specification of a function representing the commonalities of both signals is not established. In this study, we introduce a new data-driven method for estimating the transfer function from EEG signal to fMRI signal at EEG sampling rate. This approach avoids EEG subsampling to fMRI time resolution and naturally provides a test for EEG predictive power over BOLD signal fluctuations, in a well-established statistical framework. We illustrate this concept in resting state (eyes closed) and visual simultaneous fMRI-EEG experiments. The results point out that it is possible to predict the BOLD fluctuations in occipital cortex by using EEG measurements. (C) 2010 Elsevier Inc. All rights reserved.

Optimization of fMRI Processing Parameters for Simutaneous Acquisition of EEG/fMRI in Focal Epilepsy

In the context of focal epilepsy, the simultaneous combination of electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) holds a great promise as a technique by which the hemodynamic correlates of interictal spikes detected on scalp EEG can be identified. The fact that traditional EEG recordings have not been able to overcome the difficulty in correlating the ictal clinical symptoms to the onset in particular areas of the lobes, brings the need of mapping with more precision the epileptogenic cortical regions. On the other hand, fMRI suggested localizations more consistent with the ictal clinical manifestations detected. This study was developed in order to improve the knowledge about the way parameters involved in the physical and mathematical data, produced by the EEG/fMRI technique processing, would influence the final results. The evaluation of the accuracy was made by comparing the BOLD results with: the high resolution EEG maps; the malformative lesions detected in the T1 weighted MR images; and the anatomical localizations of the diagnosed symptomatology of each studied patient. The optimization of the set of parameters used, will provide an important contribution to the diagnosis of epileptogenic focuses, in patients included on an epilepsy surgery evaluation program. The results obtained allowed us to conclude that: by associating the BOLD effect with interictal spikes, the epileptogenic areas are mapped to localizations different from those obtained by the EEG maps representing the electrical potential distribution across the scalp (EEG); there is an important and solid bond between the variation of particular parameters (manipulated during the fMRI data processing) and the optimization of the final results, from which smoothing, deleted volumes, HRF (used to convolve with the activation design), and the shape of the Gamma function can be certainly emphasized.

Simultaneous EEG-fMRI at ultra-high field: Artifact prevention and safety assessment.

The simultaneous recording of scalp electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) can provide unique insights into the dynamics of human brain function, and the increased functional sensitivity offered by ultra-high field fMRI opens exciting perspectives for the future of this multimodal approach. However, simultaneous recordings are susceptible to various types of artifacts, many of which scale with magnetic field strength and can seriously compromise both EEG and fMRI data quality in recordings above 3T. The aim of the present study was to implement and characterize an optimized setup for simultaneous EEG-fMRI in humans at 7T. The effects of EEG cable length and geometry for signal transmission between the cap and amplifiers were assessed in a phantom model, with specific attention to noise contributions from the MR scanner coldheads. Cable shortening (down to 12cm from cap to amplifiers) and bundling effectively reduced environment noise by up to 84% in average power and 91% in inter-channel power variability. Subject safety was assessed and confirmed via numerical simulations of RF power distribution and temperature measurements on a phantom model, building on the limited existing literature at ultra-high field. MRI data degradation effects due to the EEG system were characterized via B0 and B1(+) field mapping on a human volunteer, demonstrating important, although not prohibitive, B1 disruption effects. With the optimized setup, simultaneous EEG-fMRI acquisitions were performed on 5 healthy volunteers undergoing two visual paradigms: an eyes-open/eyes-closed task, and a visual evoked potential (VEP) paradigm using reversing-checkerboard stimulation. EEG data exhibited clear occipital alpha modulation and average VEPs, respectively, with concomitant BOLD signal changes. On a single-trial level, alpha power variations could be observed with relative confidence on all trials; VEP detection was more limited, although statistically significant responses could be detected in more than 50% of trials for every subject. Overall, we conclude that the proposed setup is well suited for simultaneous EEG-fMRI at 7T.

Towards high-quality simultaneous EEG-fMRI at 7T: Detection and reduction of EEG artifacts due to head motion.

The enhanced functional sensitivity offered by ultra-high field imaging may significantly benefit simultaneous EEG-fMRI studies, but the concurrent increases in artifact contamination can strongly compromise EEG data quality. In the present study, we focus on EEG artifacts created by head motion in the static B0 field. A novel approach for motion artifact detection is proposed, based on a simple modification of a commercial EEG cap, in which four electrodes are non-permanently adapted to record only magnetic induction effects. Simultaneous EEG-fMRI data were acquired with this setup, at 7T, from healthy volunteers undergoing a reversing-checkerboard visual stimulation paradigm. Data analysis assisted by the motion sensors revealed that, after gradient artifact correction, EEG signal variance was largely dominated by pulse artifacts (81-93%), but contributions from spontaneous motion (4-13%) were still comparable to or even larger than those of actual neuronal activity (3-9%). Multiple approaches were tested to determine the most effective procedure for denoising EEG data incorporating motion sensor information. Optimal results were obtained by applying an initial pulse artifact correction step (AAS-based), followed by motion artifact correction (based on the motion sensors) and ICA denoising. On average, motion artifact correction (after AAS) yielded a 61% reduction in signal power and a 62% increase in VEP trial-by-trial consistency. Combined with ICA, these improvements rose to a 74% power reduction and an 86% increase in trial consistency. Overall, the improvements achieved were well appreciable at single-subject and single-trial levels, and set an encouraging quality mark for simultaneous EEG-fMRI at ultra-high field.

EEG spike source localization before and after surgery for temporal lobe epilepsy: a BOLD EEG-fMRI and independent component analysis study

Simultaneous measurements of EEG-functional magnetic resonance imaging (fMRI) combine the high temporal resolution of EEG with the distinctive spatial resolution of fMRI. The purpose of this EEG-fMRI study was to search for hemodynamic responses (blood oxygen level-dependent - BOLD responses) associated with interictal activity in a case of right mesial temporal lobe epilepsy before and after a successful selective amygdalohippocampectomy. Therefore, the study found the epileptogenic source by this noninvasive imaging technique and compared the results after removing the atrophied hippocampus. Additionally, the present study investigated the effectiveness of two different ways of localizing epileptiform spike sources, i.e., BOLD contrast and independent component analysis dipole model, by comparing their respective outcomes to the resected epileptogenic region. Our findings suggested a right hippocampus induction of the large interictal activity in the left hemisphere. Although almost a quarter of the dipoles were found near the right hippocampus region, dipole modeling resulted in a widespread distribution, making EEG analysis too weak to precisely determine by itself the source localization even by a sophisticated method of analysis such as independent component analysis. On the other hand, the combined EEG-fMRI technique made it possible to highlight the epileptogenic foci quite efficiently.

Realistic mean field forward predictions for the integration of co-registered EEG/fMRI

Brain activity can be measured non-invasively with functional imaging techniques. Each pixel in such an image represents a neural mass of about 105 to 107 neurons. Mean field models (MFMs) approximate their activity by averaging out neural variability while retaining salient underlying features, like neurotransmitter kinetics. However, MFMs incorporating the regional variability, realistic geometry and connectivity of cortex have so far appeared intractable. This lack of biological realism has led to a focus on gross temporal features of the EEG. We address these impediments and showcase a "proof of principle" forward prediction of co-registered EEG/fMRI for a full-size human cortex in a realistic head model with anatomical connectivity, see figure 1. MFMs usually assume homogeneous neural masses, isotropic long-range connectivity and simplistic signal expression to allow rapid computation with partial differential equations. But these approximations are insufficient in particular for the high spatial resolution obtained with fMRI, since different cortical areas vary in their architectonic and dynamical properties, have complex connectivity, and can contribute non-trivially to the measured signal. Our code instead supports the local variation of model parameters and freely chosen connectivity for many thousand triangulation nodes spanning a cortical surface extracted from structural MRI. This allows the introduction of realistic anatomical and physiological parameters for cortical areas and their connectivity, including both intra- and inter-area connections. Proper cortical folding and conduction through a realistic head model is then added to obtain accurate signal expression for a comparison to experimental data. To showcase the synergy of these computational developments, we predict simultaneously EEG and fMRI BOLD responses by adding an established model for neurovascular coupling and convolving "Balloon-Windkessel" hemodynamics. We also incorporate regional connectivity extracted from the CoCoMac database [1]. Importantly, these extensions can be easily adapted according to future insights and data. Furthermore, while our own simulation is based on one specific MFM [2], the computational framework is general and can be applied to models favored by the user. Finally, we provide a brief outlook on improving the integration of multi-modal imaging data through iterative fits of a single underlying MFM in this realistic simulation framework.

Multimodal (EEG-fMRI) functional connectivity study of levodopa effect in Parkinson’s disease

Aim: To assess if the intake of levodopa in patients with Parkinson’s Disease (PD) changes cerebral connectivity, as revealed by simultaneous recording of hemodynamic (functional MRI, or fMRI) and electric (electroencephalogram, EEG) signals. Particularly, we hypothesize that the strongest changes in FC will involve the motor network, which is the most impaired in PD. Methods: Eight patients with diagnosis of PD “probable”, therapy with levodopa exclusively, normal cognitive and affective status, were included. Exclusion criteria were: moderate-severe rest tremor, levodopa induced dyskinesia, evidence of gray or white matter abnormalities on structural MRI. Scalp EEG (64 channels) were acquired inside the scanner (1.5 Tesla) before and after the intake of levodopa. fMRI functional connectivity was computed from four regions of interest: right and left supplementary motor area (SMA) and right and left precentral gyrus (primary motor cortex). Weighted partial directed coherence (w-PDC) was computed in the inverse space after the removal of EEG gradient and cardioballistic artifacts. Results and discussion: fMRI group analysis shows that the intake of levodopa increases hemodynamic functional connectivity among the SMAs / primary motor cortex and: sensory-motor network itself, attention network and default mode network. w-PDC analysis shows that EEG connectivity among regions of the motor network has the tendency to decrease after the intake the levodopa; furthermore, regions belonging to the DMN have the tendency to increase their outflow toward the rest of the brain. These findings, even if in a small sample of patients, suggest that other resting state physiological functional networks, beyond the motor one, are affected in patients with PD. The behavioral and cognitive tasks corresponding to the affected networks could benefit from the intake of levodopa.

Epileptogenic developmental venous anomaly: insights from simultaneous EEG/fMRI

Developmental venous anomalies (DVAs) are associated with epileptic seizures; however, the role of DVA in the epileptogenesis is still not established. Simultaneous interictal electroencephalogram/functional magnetic resonance imaging (EEG/fMRI) recordings provide supplementary information to electroclinical data about the epileptic generators, and thus aid in the differentiation of clinically equivocal epilepsy syndromes. The main objective of our study was to characterize the epileptic network in a patient with DVA and epilepsy by simultaneous EEG/fMRI recordings. A 17-year-old woman with recently emerging generalized tonic-clonic seizures, and atypical generalized discharges, was investigated using simultaneous EEG/fMRI at the university hospital. Previous high-resolution MRI showed no structural abnormalities, except a DVA in the right frontal operculum. Interictal EEG recordings showed atypical generalized discharges, corresponding to positive focal blood oxygen level dependent (BOLD) correlates in the right frontal operculum, a region drained by the DVA. Additionally, widespread cortical bilateral negative BOLD correlates in the frontal and parietal lobes were delineated, resembling a generalized epileptic network. The EEG/fMRI recordings support a right frontal lobe epilepsy, originating in the vicinity of the DVA, propagating rapidly to both frontal and parietal lobes, as expressed on the scalp EEG by secondary bilateral synchrony. The DVA may be causative of focal epilepsies in cases where no concomitant epileptogenic lesions can be detected. Advanced imaging techniques, such as simultaneous EEG/fMRI, may thus aid in the differentiation of clinically equivocal epilepsy syndromes.

Ordinal patterns in epileptic brains: Analysis of intracranial EEG and simultaneous EEG-fMRI

Epileptic seizures are associated with high behavioral stereotypy of the patients. In the EEG of epilepsy patients characteristic signal patterns can be found during and between seizures. Here we use ordinal patterns to analyze EEGs of epilepsy patients and quantify the degree of signal determinism. Besides relative signal redundancy and the fraction of forbidden patterns we introduce the fraction of under-represented patterns as a new measure. Using the logistic map, parameter scans are performed to explore the sensitivity of the measures to signal determinism. Thereafter, application is made to two types of EEGs recorded in two epilepsy patients. Intracranial EEG shows pronounced determinism peaks during seizures. Finally, we demonstrate that ordinal patterns may be useful for improving analysis of non-invasive simultaneous EEG-fMRI.

Brain rhythms alterations and their effect in functional networks: A simultaneous EEG/fMRI study in Fixation-off epilepsy

Fixation-off sensitivity (FOS) denotes the forms of EEG abnormalities, which are elicited by elimination of central vision or fixation. The phenomenon seems to depend on variables that modulate the alpha rhythm, however, the cerebral mechanisms underlying FOS remain unclear [1]. The scarce previous fMRI findings related to FOS have shown activation in extrastriate cortex [2] and also in frontal areas [3][4]. On the other hand, simultaneous EEG-fMRI technique has been used to assess the relationship between spontaneous power fluctuations of electrical rhythms and associated fMRI signal modulations. These studies have identified that lateral frontoparietal networks show a negative correlation with alpha band in healthy subjects. This neuroanatomical pattern is related to attentional processes and cognitive resources. Moreover, a sub-beta band (17-23 Hz) has been identified with posterior cingulate, temporoparietal junction and dorso-medial prefrontal cortex activations, which correspond to the DMN [5][6].

Comparison of Methods and Co-Registration Maps of EEG and fMRI in Occipital Lobe Epilepsy

Clinically childhood occipital lobe epilepsy (OLE) manifests itself with distinct syndromes. The traditional EEG recordings have not been able to overcome the difficulty in correlating the ictal clinical symptoms to the onset in particular areas of the occipital lobes. To understand these syndromes it is important to map with more precision the epileptogenic cortical regions in OLE. Experimentally, we studied three idiopathic childhood OLE patients with EEG source analysis and with the simultaneous acquisition of EEG and fMRI, to map the BOLD effect associated with EEG spikes. The spatial overlap between the EEG and BOLD results was not very good, but the fMRI suggested localizations more consistent with the ictal clinical manifestations of each type of epileptic syndrome. Since our first results show that by associating the BOLD effect with interictal spikes the epileptogenic areas are mapped to localizations different from those calculated from EEG sources and that by using different EEG/fMRI processing methods our results differ to some extent, it is very important to compare the different methods of processing the localization of activation and develop a good methodology for obtaining co-registration maps of high resolution EEG with BOLD localizations.

Ambiti applicativi dell'EEG e della fMRI nello studio dei processi cognitivi

L’obiettivo della mia tesi è quello di presentare e confrontare due tipologie di tecniche di indagine cerebrale, l’EEG (Elettroencefalogramma) e la fMRI (Risonanza Magnetica funzionale), evidenziandone i vantaggi e gli svantaggi, e le loro applicazioni in campo medico. Successivamente è presentato lo sviluppo di un modello sperimentale volto allo studio del fenomeno della sinestesia, a partire da dati estratti mediante le tecniche precedenti.

BOLD correlates of EEG alpha phase-locking and the fMRI default mode network

Phase locking or synchronization of brain areas is a key concept of information processing in the brain. Synchronous oscillations have been observed and investigated extensively in EEG during the past decades. EEG oscillations occur over a wide frequency range. In EEG, a prominent type of oscillations is alpha-band activity, present typically when a subject is awake, but at rest with closed eyes. The spectral power of alpha rhythms has recently been investigated in simultaneous EEG/fMRI recordings, establishing a wide-range cortico-thalamic network. However, spectral power and synchronization are different measures and little is known about the correlations between BOLD effects and EEG synchronization. Interestingly, the fMRI BOLD signal also displays synchronous oscillations across different brain regions. These oscillations delineate so-called resting state networks (RSNs) that resemble the correlation patterns of simultaneous EEG/fMRI recordings. However, the nature of these BOLD oscillations and their relations to EEG activity is still poorly understood. One hypothesis is that the subunits constituting a specific RSN may be coordinated by different EEG rhythms. In this study we report on evidence for this hypothesis. The BOLD correlates of global EEG synchronization (GFS) in the alpha frequency band are located in brain areas involved in specific RSNs, e.g. the 'default mode network'. Furthermore, our results confirm the hypothesis that specific RSNs are organized by long-range synchronization at least in the alpha frequency band. Finally, we could localize specific areas where the GFS BOLD correlates and the associated RSN overlap. Thus, we claim that not only the spectral dynamics of EEG are important, but also their spatio-temporal organization.

Shifted Coupling of EEG Driving Frequencies and fMRI Resting State Networks in Schizophrenia Spectrum Disorders

INTRODUCTION: The cerebral resting state in schizophrenia is altered, as has been demonstrated separately by electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) resting state networks (RSNs). Previous simultaneous EEG/fMRI findings in healthy controls suggest that a consistent spatiotemporal coupling between neural oscillations (EEG frequency correlates) and RSN activity is necessary to organize cognitive processes optimally. We hypothesized that this coupling is disorganized in schizophrenia and related psychotic disorders, in particular regarding higher cognitive RSNs such as the default-mode (DMN) and left-working-memory network (LWMN). METHODS: Resting state was investigated in eleven patients with a schizophrenia spectrum disorder (n = 11) and matched healthy controls (n = 11) using simultaneous EEG/fMRI. The temporal association of each RSN to topographic spectral changes in the EEG was assessed by creating Covariance Maps. Group differences within, and group similarities across frequencies were estimated for the Covariance Maps. RESULTS: The coupling of EEG frequency bands to the DMN and the LWMN respectively, displayed significant similarities that were shifted towards lower EEG frequencies in patients compared to healthy controls. CONCLUSIONS: By combining EEG and fMRI, each measuring different properties of the same pathophysiology, an aberrant relationship between EEG frequencies and altered RSNs was observed in patients. RSNs of patients were related to lower EEG frequencies, indicating functional alterations of the spatiotemporal coupling. SIGNIFICANCE: The finding of a deviant and shifted coupling between RSNs and related EEG frequencies in patients with a schizophrenia spectrum disorder is significant, as it might indicate how failures in the processing of internal and external stimuli, as commonly seen during this symptomatology (i.e. thought disorders, hallucinations), arise.

Towards the utilization of EEG as a brain imaging tool.

Recent advances in signal analysis have engendered EEG with the status of a true brain mapping and brain imaging method capable of providing spatio-temporal information regarding brain (dys)function. Because of the increasing interest in the temporal dynamics of brain networks, and because of the straightforward compatibility of the EEG with other brain imaging techniques, EEG is increasingly used in the neuroimaging community. However, the full capability of EEG is highly underestimated. Many combined EEG-fMRI studies use the EEG only as a spike-counter or an oscilloscope. Many cognitive and clinical EEG studies use the EEG still in its traditional way and analyze grapho-elements at certain electrodes and latencies. We here show that this way of using the EEG is not only dangerous because it leads to misinterpretations, but it is also largely ignoring the spatial aspects of the signals. In fact, EEG primarily measures the electric potential field at the scalp surface in the same way as MEG measures the magnetic field. By properly sampling and correctly analyzing this electric field, EEG can provide reliable information about the neuronal activity in the brain and the temporal dynamics of this activity in the millisecond range. This review explains some of these analysis methods and illustrates their potential in clinical and experimental applications.