Pulse artifact detection in simultaneous EEG-fMRI recording based on EEG map topography

Altered correlation of concurrently recorded EEG-fMRI connectomes in temporal lobe epilepsy

Whole-brain functional connectivity networks (connectomes) have been characterized at different scales in humans using EEG and fMRI. Multimodal epileptic networks have also been investigated, but the relationship between EEG and fMRI defined networks on a whole-brain scale is unclear. A unified multimodal connectome description, mapping healthy and pathological networks would close this knowledge gap. Here, we characterize the spatial correlation between the EEG and fMRI connectomes in right and left temporal lobe epilepsy (rTLE/lTLE). From two centers, we acquired resting-state concurrent EEG-fMRI of 35 healthy controls and 34 TLE patients. EEG-fMRI data was projected into the Desikan brain atlas, and functional connectomes from both modalities were correlated. EEG and fMRI connectomes were moderately correlated. This correlation was increased in rTLE when compared to controls for EEG-delta/theta/alpha/beta. Conversely, multimodal correlation in lTLE was decreased in respect to controls for EEG-beta. While the alteration was global in rTLE, in lTLE it was locally linked to the default mode network. The increased multimodal correlation in rTLE and decreased correlation in lTLE suggests a modality-specific lateralized differential reorganization in TLE, which needs to be considered when comparing results from different modalities. Each modality provides distinct information, highlighting the benefit of multimodal assessment in epilepsy.

Broadly applicable methods for the detection of artefacts in electroencephalography acquired simultaneously with hemodynamic recordings

Electroencephalography (EEG) data, acquired simultaneously with magnetic resonance imaging (MRI), must be corrected for artefacts related to MR gradient switches (GS) and the cardioballistic (CB) effect. Canonical approaches require additional signal acquisition for artefact detection (e.g., MR volume onsets, ECG), without which the EEG data would be rendered uncleanable from these artefacts.•We present two broadly applicable methods for artefact detection based on peak detection combined with temporal constraints with respect to periodicity directly from the EEG data itself; no additional signals are required. We validated the performance of our methods versus the two canonical approaches for detection of GS/CB artefact, respectively, on 26 healthy human EEG-functional MRI resting-state datasets. Utilising various performance metrics, we found our methods to perform as well as - and sometimes better than - the canonical standard approaches. With as little as one EEG channel recording, our methods can be applied to detect GS/CB artefacts in EEG data acquired simultaneously with MRI in the absence of MR volume onsets and/or an ECG recording. The detected artefact onsets can then be fed into the standard artefact correction software.

Self-Regulation of Attention in Children in a Virtual Classroom Environment: A Feasibility Study

Attention is a crucial cognitive function that enables us to selectively focus on relevant information from the surrounding world to achieve our goals. Impairments in sustained attention pose challenges, particularly in children with attention deficit hyperactivity disorder, a neurodevelopmental disorder characterized by impulsive and inattentive behavior. While psychostimulant medications are the most effective ADHD treatment, they often yield unwanted side effects, making it crucial to explore non-pharmacological treatments. We propose a groundbreaking protocol that combines electroencephalography-based neurofeedback with virtual reality (VR) as an innovative approach to address attention deficits. By integrating a virtual classroom environment, we aim to enhance the transferability of attentional control skills while simultaneously increasing motivation and interest among children. The present study demonstrates the feasibility of this approach through an initial assessment involving a small group of healthy children, showcasing its potential for future evaluation in ADHD children. Preliminary results indicate high engagement and positive feedback. Pre- and post-protocol assessments via EEG and fMRI recordings suggest changes in attentional function. Further validation is required, but this protocol is a significant advancement in neurofeedback therapy for ADHD. The integration of EEG-NFB and VR presents a novel avenue for enhancing attentional control and addressing behavioral challenges in children with ADHD.

Preservation of EEG spectral power features during simultaneous EEG-fMRI

Introduction

Electroencephalographic (EEG) data quality is severely compromised when recorded inside the magnetic resonance (MR) environment. Here we characterized the impact of the ballistocardiographic (BCG) artifact on resting-state EEG spectral properties and compared the effectiveness of seven common BCG correction methods to preserve EEG spectral features. We also assessed if these methods retained posterior alpha power reactivity to an eyes closure-opening (EC-EO) task and compared the results from EEG-informed fMRI analysis using different BCG correction approaches.

Method

Electroencephalographic data from 20 healthy young adults were recorded outside the MR environment and during simultaneous fMRI acquisition. The gradient artifact was effectively removed from EEG-fMRI acquisitions using Average Artifact Subtraction (AAS). The BCG artifact was corrected with seven methods: AAS, Optimal Basis Set (OBS), Independent Component Analysis (ICA), OBS followed by ICA, AAS followed by ICA, PROJIC-AAS and PROJIC-OBS. EEG signal preservation was assessed by comparing the spectral power of traditional frequency bands from the corrected rs-EEG-fMRI data with the data recorded outside the scanner. We then assessed the preservation of posterior alpha functional reactivity by computing the ratio between the EC and EO conditions during the EC-EO task. EEG-informed fMRI analysis of the EC-EO task was performed using alpha power-derived BOLD signal predictors obtained from the EEG signals corrected with different methods.

Results

The BCG artifact caused significant distortions (increased absolute power, altered relative power) across all frequency bands. Artifact residuals/signal losses were present after applying all correction methods. The EEG reactivity to the EC-EO task was better preserved with ICA-based correction approaches, particularly when using ICA feature extraction to isolate alpha power fluctuations, which allowed to accurately predict hemodynamic signal fluctuations during the EEG-informed fMRI analysis.

Discussion

Current software solutions for the BCG artifact problem offer limited efficiency to preserve the EEG spectral power properties using this particular EEG setup. The state-of-the-art approaches tested here can be further refined and should be combined with hardware implementations to better preserve EEG signal properties during simultaneous EEG-fMRI. Existing and novel BCG artifact correction methods should be validated by evaluating signal preservation of both ERPs and spontaneous EEG spectral power.

Data-driven beamforming technique to attenuate ballistocardiogram artefacts in electroencephalography-functional magnetic resonance imaging without detecting cardiac pulses in electrocardiography recordings

Simultaneous recording of electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) is a very promising non‐invasive neuroimaging technique. However, EEG data obtained from the simultaneous EEG–fMRI are strongly influenced by MRI‐related artefacts, namely gradient artefacts (GA) and ballistocardiogram (BCG) artefacts. When compared to the GA correction, the BCG correction is more challenging to remove due to its inherent variabilities and dynamic changes over time. The standard BCG correction (i.e., average artefact subtraction [AAS]), require detecting cardiac pulses from simultaneous electrocardiography (ECG) recording. However, ECG signals are also distorted and will become problematic for detecting reliable cardiac peaks. In this study, we focused on a beamforming spatial filtering technique to attenuate all unwanted source activities outside of the brain. Specifically, we applied the beamforming technique to attenuate the BCG artefact in EEG–fMRI, and also to recover meaningful task‐based neural signals during an attentional network task (ANT) which required participants to identify visual cues and respond accurately. We analysed EEG–fMRI data in 20 healthy participants during the ANT, and compared four different BCG corrections (non‐BCG corrected, AAS BCG corrected, beamforming + AAS BCG corrected, beamforming BCG corrected). We demonstrated that the beamforming approach did not only significantly reduce the BCG artefacts, but also significantly recovered the expected task‐based brain activity when compared to the standard AAS correction. This data‐driven beamforming technique appears promising especially for longer data acquisition of sleep and resting EEG–fMRI. Our findings extend previous work regarding the recovery of meaningful EEG signals by an optimized suppression of MRI‐related artefacts.

Artifact Reduction in Simultaneous EEG-fMRI: A Systematic Review of Methods and Contemporary Usage

Simultaneous electroencephalography-functional MRI (EEG-fMRI) is a technique that combines temporal (largely from EEG) and spatial (largely from fMRI) indicators of brain dynamics. It is useful for understanding neuronal activity during many different event types, including spontaneous epileptic discharges, the activity of sleep stages, and activity evoked by external stimuli and decision-making tasks. However, EEG recorded during fMRI is subject to imaging, pulse, environment and motion artifact, causing noise many times greater than the neuronal signals of interest. Therefore, artifact removal methods are essential to ensure that artifacts are accurately removed, and EEG of interest is retained. This paper presents a systematic review of methods for artifact reduction in simultaneous EEG-fMRI from literature published since 1998, and an additional systematic review of EEG-fMRI studies published since 2016. The aim of the first review is to distill the literature into clear guidelines for use of simultaneous EEG-fMRI artifact reduction methods, and the aim of the second review is to determine the prevalence of artifact reduction method use in contemporary studies. We find that there are many published artifact reduction techniques available, including hardware, model based, and data-driven methods, but there are few studies published that adequately compare these methods. In contrast, recent EEG-fMRI studies show overwhelming use of just one or two artifact reduction methods based on literature published 15–20 years ago, with newer methods rarely gaining use outside the group that developed them. Surprisingly, almost 15% of EEG-fMRI studies published since 2016 fail to adequately describe the methods of artifact reduction utilized. We recommend minimum standards for reporting artifact reduction techniques in simultaneous EEG-fMRI studies and suggest that more needs to be done to make new artifact reduction techniques more accessible for the researchers and clinicians using simultaneous EEG-fMRI.

The relationship between EEG and fMRI connectomes is reproducible across simultaneous EEG-fMRI studies from 1.5T to 7T

Both electroencephalography (EEG) and functional Magnetic Resonance Imaging (fMRI) are non-invasive methods that show complementary aspects of human brain activity. Despite measuring different proxies of brain activity, both the measured blood-oxygenation (fMRI) and neurophysiological recordings (EEG) are indirectly coupled. The electrophysiological and BOLD signal can map the underlying functional connectivity structure at the whole brain scale at different timescales. Previous work demonstrated a moderate but significant correlation between resting-state functional connectivity of both modalities, however there is a wide range of technical setups to measure simultaneous EEG-fMRI and the reliability of those measures between different setups remains unknown. This is true notably with respect to different magnetic field strengths (low and high field) and different spatial sampling of EEG (medium to high-density electrode coverage). Here, we investigated the reproducibility of the bimodal EEG-fMRI functional connectome in the most comprehensive resting-state simultaneous EEG-fMRI dataset compiled to date including a total of 72 subjects from four different imaging centers. Data was acquired from 1.5T, 3T and 7T scanners with simultaneously recorded EEG using 64 or 256 electrodes. We demonstrate that the whole-brain monomodal connectivity reproducibly correlates across different datasets and that a moderate crossmodal correlation between EEG and fMRI connectivity of r ≈ 0.3 can be reproducibly extracted in low- and high-field scanners. The crossmodal correlation was strongest in the EEG-β frequency band but exists across all frequency bands. Both homotopic and within intrinsic connectivity network (ICN) connections contributed the most to the crossmodal relationship. This study confirms, using a considerably diverse range of recording setups, that simultaneous EEG-fMRI offers a consistent estimate of multimodal functional connectomes in healthy subjects that are dominantly linked through a functional core of ICNs across spanning across the different timescales measured by EEG and fMRI. This opens new avenues for estimating the dynamics of brain function and provides a better understanding of interactions between EEG and fMRI measures. This observed level of reproducibility also defines a baseline for the study of alterations of this coupling in pathological conditions and their role as potential clinical markers.

Modulation of epileptic networks by transient interictal epileptic activity: A dynamic approach to simultaneous EEG-fMRI

•

EEG-fMRI has been instrumental in characterizing brain networks in epilepsy.

•

Its value is documented in the pre-surgical assessment of drug-resistant epilepsy.

•

The delineation of brain areas to resect is fundamental for the post-surgical outcome.

•

Standard EEG-fMRI in epilepsy assesses static functional connectivity of the network.

•

EEG-fMRI dynamic connectivity identifies transitory features of specific connections.

•

We integrate dynamic fMRI connectivity and dynamic patterns of simultaneous scalp EEG.

•

This allows to better characterize the spatiotemporal aspects of epileptic networks.

•

This may help in more efficiently target the surgical intervention.

Epileptic networks, defined as brain regions involved in epileptic brain activity, have been mapped by functional connectivity in simultaneous electroencephalography and functional magnetic resonance imaging (EEG-fMRI) recordings. This technique allows to define brain hemodynamic changes, measured by the Blood Oxygen Level Dependent (BOLD) signal, associated to the interictal epileptic discharges (IED), which together with ictal events constitute a signature of epileptic disease. Given the highly time-varying nature of epileptic activity, a dynamic functional connectivity (dFC) analysis of EEG-fMRI data appears particularly suitable, having the potential to identify transitory features of specific connections in epileptic networks. In the present study, we propose a novel method, defined dFC-EEG, that integrates dFC assessed by fMRI with the information recorded by simultaneous scalp EEG, in order to identify the connections characterised by a dynamic profile correlated with the occurrence of IED, forming the dynamic epileptic subnetwork. Ten patients with drug-resistant focal epilepsy were included, with different aetiology and showing a widespread (or multilobar) BOLD activation, defined as involving at least two distinct clusters, located in two different lobes and/or extended to the hemisphere contralateral to the epileptic focus. The epileptic focus was defined from the IED-related BOLD map. Regions involved in the occurrence of interictal epileptic activity; i.e., forming the epileptic network, were identified by a general linear model considering the timecourse of the fMRI-defined focus as main regressor. dFC between these regions was assessed with a sliding-window approach. dFC timecourses were then correlated with the sliding-window variance of the IED signal (VarIED), to identify connections whose dynamics related to the epileptic activity; i.e., the dynamic epileptic subnetwork. As expected, given the very different clinical picture of each individual, the extent of this subnetwork was highly variable across patients, but was but was reduced of at least 30% with respect to the initially identified epileptic network in 9/10 patients. The connections of the dynamic subnetwork were most commonly close to the epileptic focus, as reflected by the laterality index of the subnetwork connections, reported higher than the one within the original epileptic network. Moreover, the correlation between dFC timecourses and VarIED was predominantly positive, suggesting a strengthening of the dynamic subnetwork associated to the occurrence of IED. The integration of dFC and scalp IED offers a more specific description of the epileptic network, identifying connections strongly influenced by IED. These findings could be relevant in the pre-surgical evaluation for the resection or disconnection of the epileptogenic zone and help in reaching a better post-surgical outcome. This would be particularly important for patients characterised by a widespread pathological brain activity which challenges the surgical intervention.

Concurrent EEG- and fMRI-derived functional connectomes exhibit linked dynamics

Long-range connectivity has become the most studied feature of human functional Magnetic Resonance Imaging (fMRI), yet the spatial and temporal relationship between its whole-brain dynamics and electrophysiological connectivity remains largely unknown. FMRI-derived functional connectivity exhibits spatial reconfigurations or time-varying dynamics at infraslow (<0.1Hz) speeds. Conversely, electrophysiological connectivity is based on cross-region coupling of fast oscillations (~1-100Hz). It is unclear whether such fast oscillation-based coupling varies at infraslow speeds, temporally coinciding with infraslow dynamics across the fMRI-based connectome. If so, does the association of fMRI-derived and electrophysiological dynamics spatially vary over the connectome across the functionally distinct electrophysiological oscillation bands? In two concurrent electroencephalography (EEG)-fMRI resting-state datasets, oscillation-based coherence in all canonical bands (delta through gamma) indeed reconfigured at infraslow speeds in tandem with fMRI-derived connectivity changes in corresponding region-pairs. Interestingly, irrespective of EEG frequency-band the cross-modal tie of connectivity dynamics comprised a large proportion of connections distributed across the entire connectome. However, there were frequency-specific differences in the relative strength of the cross-modal association. This association was strongest in visual to somatomotor connections for slower EEG-bands, and in connections involving the Default Mode Network for faster EEG-bands. Methodologically, the findings imply that neural connectivity dynamics can be reliably measured by fMRI despite heavy susceptibility to noise, and by EEG despite shortcomings of source reconstruction. Biologically, the findings provide evidence that contrast with known territories of oscillation power, oscillation coupling in all bands slowly reconfigures in a highly distributed manner across the whole-brain connectome.

Intrinsic connectome organization across temporal scales: New insights from cross-modal approaches

The discovery of a stable, whole-brain functional connectivity organization that is largely independent of external events has drastically extended our view of human brain function. However, this discovery has been primarily based on functional magnetic resonance imaging (fMRI). The role of this whole-brain organization in fast oscillation-based connectivity as measured, for example, by electroencephalography (EEG) and magnetoencephalography (MEG) is only beginning to emerge. Here, we review studies of intrinsic connectivity and its whole-brain organization in EEG, MEG, and intracranial electrophysiology with a particular focus on direct comparisons to connectome studies in fMRI. Synthesizing this literature, we conclude that irrespective of temporal scale over four orders of magnitude, intrinsic neurophysiological connectivity shows spatial similarity to the connectivity organization commonly observed in fMRI. A shared structural connectivity basis and cross-frequency coupling are possible mechanisms contributing to this similarity. Acknowledging that a stable whole-brain organization governs long-range coupling across all timescales of neural processing motivates researchers to take "baseline" intrinsic connectivity into account when investigating brain-behavior associations, and further encourages more widespread exploration of functional connectomics approaches beyond fMRI by using EEG and MEG modalities.

Investigating the variability of cardiac pulse artifacts across heartbeats in simultaneous EEG-fMRI recordings: A 7T study

Electroencephalography (EEG) recordings performed in magnetic resonance imaging (MRI) scanners are affected by complex artifacts caused by heart function, often termed pulse artifacts (PAs). PAs can strongly compromise EEG data quality, and remain an open problem for EEG-fMRI. This study investigated the properties and mechanisms of PA variability across heartbeats, which has remained largely unaddressed to date, and evaluated its impact on PA correction approaches. Simultaneous EEG-fMRI was performed at 7T on healthy participants at rest or under visual stimulation, with concurrent recordings of breathing and cardiac activity. PA variability was found to contribute to EEG variance with more than 500 μV² at 7T, which extrapolates to 92 μV² at 3T. Clustering analyses revealed that PA variability not only is linked to variations in head position/orientation, as previously hypothesized, but also, and more importantly, to the respiratory cycle and to heart rate fluctuations. The latter mechanisms are associated to short-timescale variability (even across consecutive heartbeats), and their importance varied across EEG channels. In light of this PA variability, three PA correction techniques were compared: average artifact subtraction (AAS), optimal basis sets (OBS), and an approach based on K-means clustering. All methods allowed the recovery of visual evoked potentials from the EEG data; nonetheless, OBS and K-means tended to outperform AAS, likely due to the inability of the latter in modeling short-timescale variability. Altogether, these results offer novel insights into the dynamics and underlying mechanisms of the pulse artifact, with important consequences for its correction, relevant to most EEG-fMRI applications.

Exploring the relative efficacy of motion artefact correction techniques for EEG data acquired during simultaneous fMRI

Simultaneous EEG‐fMRI allows multiparametric characterisation of brain function, in principle enabling a more complete understanding of brain responses; unfortunately the hostile MRI environment severely reduces EEG data quality. Simply eliminating data segments containing gross motion artefacts [MAs] (generated by movement of the EEG system and head in the MRI scanner's static magnetic field) was previously believed sufficient. However recently the importance of removal of all MAs has been highlighted and new methods developed. A systematic comparison of the ability to remove MAs and retain underlying neuronal activity using different methods of MA detection and post‐processing algorithms is needed to guide the neuroscience community. Using a head phantom, we recorded MAs while simultaneously monitoring the motion using three different approaches: Reference Layer Artefact Subtraction (RLAS), Moiré Phase Tracker (MPT) markers and Wire Loop Motion Sensors (WLMS). These EEG recordings were combined with EEG responses to simple visual tasks acquired on a subject outside the MRI environment. MAs were then corrected using the motion information collected with each of the methods combined with different analysis pipelines. All tested methods retained the neuronal signal. However, often the MA was not removed sufficiently to allow accurate detection of the underlying neuronal signal. We show that the MA is best corrected using the RLAS combined with post‐processing using a multichannel, recursive least squares (M‐RLS) algorithm. This method needs to be developed further to enable practical utility; thus, WLMS combined with M‐RLS currently provides the best compromise between EEG data quality and practicalities of motion detection.

Adaptive optimal basis set for BCG artifact removal in simultaneous EEG-fMRI

Electroencephalography (EEG) signals recorded during simultaneous functional magnetic resonance imaging (fMRI) are contaminated by strong artifacts. Among these, the ballistocardiographic (BCG) artifact is the most challenging, due to its complex spatio-temporal dynamics associated with ongoing cardiac activity. The presence of BCG residuals in EEG data may hide true, or generate spurious correlations between EEG and fMRI time-courses. Here, we propose an adaptive Optimal Basis Set (aOBS) method for BCG artifact removal. Our method is adaptive, as it can estimate the delay between cardiac activity and BCG occurrence on a beat-to-beat basis. The effective creation of an optimal basis set by principal component analysis (PCA) is therefore ensured by a more accurate alignment of BCG occurrences. Furthermore, aOBS can automatically estimate which components produced by PCA are likely to be BCG artifact-related and therefore need to be removed. The aOBS performance was evaluated on high-density EEG data acquired with simultaneous fMRI in healthy subjects during visual stimulation. As aOBS enables effective reduction of BCG residuals while preserving brain signals, we suggest it may find wide application in simultaneous EEG-fMRI studies.

Heart-Brain Interactions in the MR Environment: Characterization of the Ballistocardiogram in EEG Signals Collected During Simultaneous fMRI

The ballistocardiographic (BCG) artifact is linked to cardiac activity and occurs in electroencephalographic (EEG) recordings acquired inside the magnetic resonance (MR) environment. Its variability in terms of amplitude, waveform shape and spatial distribution over subject's scalp makes its attenuation a challenging task. In this study, we aimed to provide a detailed characterization of the BCG properties, including its temporal dependency on cardiac events and its spatio-temporal dynamics. To this end, we used high-density EEG data acquired during simultaneous functional MR imaging in six healthy volunteers. First, we investigated the relationship between cardiac activity and BCG occurrences in the EEG recordings. We observed large variability in the delay between ECG and subsequent BCG events (ECG-BCG delay) across subjects and non-negligible epoch-by-epoch variations at the single subject level. The inspection of spatial-temporal variations revealed a prominent non-stationarity of the BCG signal. We identified five main BCG waves, which were common across subjects. Principal component analysis revealed two spatially distinct patterns to explain most of the variance (85% in total). These components are possibly related to head rotation and pulse-driven scalp expansion, respectively. Our results may inspire the development of novel, more effective methods for the removal of the BCG, capable of isolating and attenuating artifact occurrences while preserving true neuronal activity.

Evaluating the impact of fast-fMRI on dynamic functional connectivity in an event-based paradigm

The human brain is known to contain several functional networks that interact dynamically. Therefore, it is desirable to analyze the temporal features of these networks by dynamic functional connectivity (dFC). A sliding window approach was used in an event-related fMRI (visual stimulation using checkerboards) to assess the impact of repetition time (TR) and window size on the temporal features of BOLD dFC. In addition, we also examined the spatial distribution of dFC and tested the feasibility of this approach for the analysis of interictal epileptiforme discharges. 15 healthy controls (visual stimulation paradigm) and three patients with epilepsy (EEG-fMRI) were measured with EPI-fMRI. We calculated the functional connectivity degree (FCD) by determining the total number of connections of a given voxel above a predefined threshold based on Pearson correlation. FCD could capture hemodynamic changes relative to stimulus onset in controls. A significant effect of TR and window size was observed on FCD estimates. At a conventional TR of 2.6 s, FCD values were marginal compared to FCD values using sub-seconds TRs achievable with multiband (MB) fMRI. Concerning window sizes, a specific maximum of FCD values (inverted u-shape behavior) was found for each TR, indicating a limit to the possible gain in FCD for increasing window size. In patients, a dynamic FCD change was found relative to the onset of epileptiform EEG patterns, which was compatible with their clinical semiology. Our findings indicate that dynamic FCD transients are better detectable with sub-second TR than conventional TR. This approach was capable of capturing neuronal connectivity across various regions of the brain, indicating a potential to study the temporal characteristics of interictal epileptiform discharges and seizures in epilepsy patients or other brain diseases with brief events.

Effectiveness of Reference Signal-Based Methods for Removal of EEG Artifacts Due to Subtle Movements During fMRI Scanning

OBJECTIVE

Subtle motion of an epileptic patient examined with co-registered EEG and functional MRI (EEG-fMRI) may often lead to spurious fMRI activation patterns when true epileptic spikes are contaminated with motion artefacts. In recent years, methods relying on reference signals for correcting these subtle movements in the EEG have emerged. In this study, the performance of two reference-based devices are compared to the template-based method with regard to their ability to remove movement-related artifacts in EEG measured during scanning.

METHODS

Measurements were performed with a novel double layer cap consisting of 29 EEG and 29 reference electrodes, and with a current loop cap consisting of 60 electrodes and three current loop wires attached to the cap. EEG was acquired inside the scanner during resting state, as well as when the subject was performing a cued movement task. For the double layer cap recordings, newly developed artifact removal algorithms are introduced and both reference signal-based methods are compared to a template-based correction method.

RESULTS

The BCG artifacts occurring at resting state could be removed successfully by both the reference signal-based methods as well as by the template-based method. However, the reference signal-based methods were also capable of removing EEG artifacts induced by subtle movements, whereas the template-based method failed to remove these artifacts.

CONCLUSION

Reference signal-based methods enable to correct for artifacts due to subtle movements, which are not removed by commonly used template-based removal algorithms.

SIGNIFICANCE

Sensitivity of EEG-fMRI analysis in patients with focal epilepsy is improved by avoiding erroneous detections of subtle movements as epileptic spikes in the EEG.

Epileptic networks are strongly connected with and without the effects of interictal discharges

OBJECTIVE

Epilepsy is increasingly considered as the dysfunction of a pathologic neuronal network (epileptic network) rather than a single focal source. We aimed to assess the interactions between the regions that comprise the epileptic network and to investigate their dependence on the occurrence of interictal epileptiform discharges (IEDs).

METHODS

We analyzed resting state simultaneous electroencephalography-functional magnetic resonance imaging (EEG-fMRI) recordings in 10 patients with drug-resistant focal epilepsy with multifocal IED-related blood oxygen level-dependent (BOLD) responses and a maximum t-value in the IED field. We computed functional connectivity (FC) maps of the epileptic network using two types of seed: (1) a 10-mm diameter sphere centered in the global maximum of IED-related BOLD map, and (2) the independent component with highest correlation to the IED-related BOLD map, named epileptic component. For both approaches, we compared FC maps before and after regressing out the effect of IEDs in terms of maximum and mean t-values and percentage of map overlap.

RESULTS

Maximum and mean FC maps t-values were significantly lower after regressing out IEDs at the group level (p < 0.01). Overlap extent was 85% ± 12% and 87% ± 12% when the seed was the 10-mm diameter sphere and the epileptic component, respectively.

SIGNIFICANCE

Regions involved in a specific epileptic network show coherent BOLD fluctuations independent of scalp EEG IEDs. FC topography and strength is largely preserved by removing the IED effect. This could represent a signature of a sustained pathologic network with contribution from epileptic activity invisible to the scalp EEG.

Ballistocardiogram artifact correction taking into account physiological signal preservation in simultaneous EEG-fMRI

The ballistocardiogram (BCG) artifact is currently one of the most challenging in the EEG acquired concurrently with fMRI, with correction invariably yielding residual artifacts and/or deterioration of the physiological signals of interest. In this paper, we propose a family of methods whereby the EEG is decomposed using Independent Component Analysis (ICA) and a novel approach for the selection of BCG-related independent components (ICs) is used (PROJection onto Independent Components, PROJIC). Three ICA-based strategies for BCG artifact correction are then explored: 1) BCG-related ICs are removed from the back-reconstruction of the EEG (PROJIC); and 2-3) BCG-related ICs are corrected for the artifact occurrences using an Optimal Basis Set (OBS) or Average Artifact Subtraction (AAS) framework, before back-projecting all ICs onto EEG space (PROJIC-OBS and PROJIC-AAS, respectively). A novel evaluation pipeline is also proposed to assess the methods performance, which takes into account not only artifact but also physiological signal removal, allowing for a flexible weighting of the importance given to physiological signal preservation. This evaluation is used for the group-level parameter optimization of each algorithm on simultaneous EEG-fMRI data acquired using two different setups at 3T and 7T. Comparison with state-of-the-art BCG correction methods showed that PROJIC-OBS and PROJIC-AAS outperformed the others when priority was given to artifact removal or physiological signal preservation, respectively, while both PROJIC-AAS and AAS were in general the best choices for intermediate trade-offs. The impact of the BCG correction on the quality of event-related potentials (ERPs) of interest was assessed in terms of the relative reduction of the standard error (SE) across trials: 26/66%, 32/62% and 18/61% were achieved by, respectively, PROJIC, PROJIC-OBS and PROJIC-AAS, for data collected at 3T/7T. Although more significant improvements were achieved at 7T, the results were qualitatively comparable for both setups, which indicate the wide applicability of the proposed methodologies and recommendations.

Presurgical brain mapping in epilepsy using simultaneous EEG and functional MRI at ultra-high field: feasibility and first results

OBJECTIVES

The aim of this study was to demonstrate that eloquent cortex and epileptic-related hemodynamic changes can be safely and reliably detected using simultaneous electroencephalography (EEG)-functional magnetic resonance imaging (fMRI) recordings at ultra-high field (UHF) for clinical evaluation of patients with epilepsy.

MATERIALS AND METHODS

Simultaneous EEG-fMRI was acquired at 7 T using an optimized setup in nine patients with lesional epilepsy. According to the localization of the lesion, mapping of eloquent cortex (language and motor) was also performed in two patients.

RESULTS

Despite strong artifacts, efficient correction of intra-MRI EEG could be achieved with optimized artifact removal algorithms, allowing robust identification of interictal epileptiform discharges. Noise-sensitive topography-related analyses and electrical source localization were also performed successfully. Localization of epilepsy-related hemodynamic changes compatible with the lesion were detected in three patients and concordant with findings obtained at 3 T. Local loss of signal in specific regions, essentially due to B 1 inhomogeneities were found to depend on the geometric arrangement of EEG leads over the cap.

CONCLUSION

These results demonstrate that presurgical mapping of epileptic networks and eloquent cortex is both safe and feasible at UHF, with the benefits of greater spatial resolution and higher blood-oxygenation-level-dependent sensitivity compared with the more traditional field strength of 3 T.