Towards motion insensitive EEG-fMRI: Correcting motion-induced voltages and gradient artefact instability in EEG using an fMRI prospective motion correction (PMC) system

Camera-based Prospective Motion Correction in Paediatric Epilepsy Patients Enables EEG-fMRI Localization Even in High-motion States

BACKGROUND

EEG-fMRI is a useful additional test to localize the epileptogenic zone (EZ) particularly in MRI negative cases. However subject motion presents a particular challenge owing to its large effects on both MRI and EEG signal. Traditionally it is assumed that prospective motion correction (PMC) of fMRI precludes EEG artifact correction.

METHODS

Children undergoing presurgical assessment at Great Ormond Street Hospital were included into the study. PMC of fMRI was done using a commercial system with a Moiré Phase Tracking marker and MR-compatible camera. For retrospective EEG correction both a standard and a motion educated EEG artefact correction (REEGMAS) were compared to each other.

RESULTS

Ten children underwent simultaneous EEG-fMRI. Overall head movement was high (mean RMS velocity < 1.5 mm/s) and showed high inter- and intra-individual variability. Comparing motion measured by the PMC camera and the (uncorrected residual) motion detected by realignment of fMRI images, there was a five-fold reduction in motion from its prospective correction. Retrospective EEG correction using both standard approaches and REEGMAS allowed the visualization and identification of physiological noise and epileptiform discharges. Seven of 10 children had significant maps, which were concordant with the clinical EZ hypothesis in 6 of these 7.

CONCLUSION

To our knowledge this is the first application of camera-based PMC for MRI in a pediatric clinical setting. Despite large amount of movement PMC in combination with retrospective EEG correction recovered data and obtained clinically meaningful results during high levels of subject motion. Practical limitations may currently limit the widespread use of this technology.

Dynamic networks of P300-related process

P300 as an effective biomarker to index attention and memory has been widely used for brain-computer interface, cognitive evaluation, and clinical diagnosis. To evoke clear P300, an oddball paradigm consisting of two types of stimuli, i.e., infrequent target stimuli and frequent standard stimuli, is usually used. However, to simply and quickly explore the P300-related process, previous studies predominately focused on the target condition but ignored the fusion of target and standard conditions, as well as the difference of brain networks between them. Therefore, in this study, we used the hidden Markov model to investigate the fused multi-conditional electroencephalogram dataset of P300, aiming to effectively identify the underlying brain networks and explore the difference between conditions. Specifically, the inferred networks, including their transition sequences and spatial distributions, were scrutinized first. Then, we found that the difference between target and standard conditions was mainly concentrated in two phases. One was the stimulation phase that mainly related to the cortical activities of the postcentral gyrus and superior parietal lobule, and the other corresponded to the response phase that involved the activities of superior and medial frontal gyri. This might be attributed to distinct cognitive functions, as the stimulation phase is associated with visual information integration whereas the response phase involves stimulus discrimination and behavior control. Taken together, the current work explored dynamic networks underlying the P300-related process and provided a complementary understanding of distinct P300 conditions, which may contribute to the design of P300-related brain-machine systems.

Tracking of rigid head motion during MRI using an EEG system

Purpose

To demonstrate a novel method for tracking of head movements during MRI using electroencephalography (EEG) hardware for recording signals induced by native imaging gradients.

Theory and Methods

Gradient switching during simultaneous EEG–fMRI induces distortions in EEG signals, which depend on subject head position and orientation. When EEG electrodes are interconnected with high‐impedance carbon wire loops, the induced voltages are linear combinations of the temporal gradient waveform derivatives. We introduce head tracking based on these signals (CapTrack) involving 3 steps: (1) phantom scanning is used to characterize the target sequence and a fast calibration sequence; (2) a linear relation between changes of induced signals and head pose is established using the calibration sequence; and (3) induced signals recorded during target sequence scanning are used for tracking and retrospective correction of head movement without prolonging the scan time of the target sequence. Performance of CapTrack is compared directly to interleaved navigators.

Results

Head‐pose tracking at 27.5 Hz during echo planar imaging (EPI) was demonstrated with close resemblance to rigid body alignment (mean absolute difference: [0.14 0.38 0.15]‐mm translation, [0.30 0.27 0.22]‐degree rotation). Retrospective correction of 3D gradient‐echo imaging shows an increase of average edge strength of 12%/−0.39% for instructed/uninstructed motion with CapTrack pose estimates, with a tracking interval of 1561 ms and high similarity to interleaved navigator estimates (mean absolute difference: [0.13 0.33 0.12] mm, [0.28 0.15 0.22] degrees).

Conclusion

Motion can be estimated from recordings of gradient switching with little or no sequence modification, optionally in real time at low computational burden and synchronized to image acquisition, using EEG equipment already found at many research institutions.

Simultaneous EEG-fMRI: What Have We Learned and What Does the Future Hold?

Simultaneous EEG-fMRI has developed into a mature measurement technique in the past 25 years. During this time considerable technical and analytical advances have been made, enabling valuable scientific contributions to a range of research fields. This review will begin with an introduction to the measurement principles involved in EEG and fMRI and the advantages of combining these methods. The challenges faced when combining the two techniques will then be considered. An overview of the leading application fields where EEG-fMRI has made a significant contribution to the scientific literature and emerging applications in EEG-fMRI research trends is then presented.

Application of Compressed Sensing 3D MR cholangiopancreatography (CS-MRCP) with Contact-Free Physiological Monitoring (CFPM) for Pancreaticobiliary Disorders

RATIONAL AND OBJECTIVES

To prospectively evaluate the clinical feasibility of the magnetic resonance cholangiopancreatography (MRCP) protocol using both contact-free physiological monitoring (CFPM) and compressed sensing (CS) (CS-CFPM-MRCP) and to compare its performance with that of the standard navigator-triggered (NT) CS-NT-MRCP and NT-MRCP.

MATERIALS AND METHODS

A total of 63 patients (36 males, 27 females, age range: 18-83 years, mean age: 52.30 ± 15.70 years) suspected with duct-related pathologies were prospectively enrolled and performed the three MRCP protocols randomly. The acquisition time was compared. The pancreaticobiliary system was divided into 12 segments and evaluated based on a five-point Likert scale and compared by the Friedman test with a post hoc test. The diagnostic performance of the 3 MRCP was evaluated by the AUC value and compared by Delong's test. The interobserver agreement was evaluated by Kendall's W test.

RESULTS

Compared to NT-MRCP, the acquisition time of CS-NT-MRCP and CS-CFPM-MRCP was significantly decreased (both p < 0.001). There is no significant difference in the overall imaging quality (p > 0.05) between the NT-MRCP and CS-CFPM-MRCP protocols. CS-CFPM-MRCP depicted pancreatic duct and intrahepatic ducts better than CS-NT-MRCP (all p < 0.05) and was comparable with that of the NT-MRCP (all p > 0.05). For identification of abnormalities and diseases associated with MPD anatomy, the mean AUC value for NT-MRCP and CS-CFPM-MRCP were 0.896 (95%CI: 0.834, 0.958) and 0.905 (95%CI: 0.846, 0.964), which were significantly higher when compared to that for CS-NT-MRCP (0.713 [95%CI:0.622, 0.805]) (p = 0.001 and < 0.001). All evaluations showed good to excellent agreement (0.619-0.897).

CONCLUSION

The combination of CS and CFPM is considered feasible for shortening the scan time of 3D free breath MRCP without impairing the imaging quality and CS-CFPM-MRCP is considered feasible for patients suspected with pancreaticobiliary diseases.

Unified Retrospective EEG Motion Educated Artefact Suppression for EEG-fMRI to Suppress Magnetic Field Gradient Artefacts During Motion

The data quality of simultaneously acquired electroencephalography and functional magnetic resonance imaging (EEG-fMRI) can be strongly affected by motion. Recent work has shown that the quality of fMRI data can be improved by using a Moiré-Phase-Tracker (MPT)-camera system for prospective motion correction. The use of the head position acquired by the MPT-camera-system has also been shown to correct motion-induced voltages, ballistocardiogram (BCG) and gradient artefact residuals separately. In this work we show the concept of an integrated framework based on the general linear model to provide a unified motion informed model of in-MRI artefacts. This model (retrospective EEG motion educated gradient artefact suppression, REEG-MEGAS) is capable of correcting voltage-induced, BCG and gradient artefact residuals of EEG data acquired simultaneously with prospective motion corrected fMRI. In our results, we have verified that applying REEG-MEGAS correction to EEG data acquired during subject motion improves the data quality in terms of motion induced voltages and also GA residuals in comparison to standard Artefact Averaging Subtraction and Retrospective EEG Motion Artefact Suppression. Besides that, we provide preliminary evidence that although adding more regressors to a model may slightly affect the power of physiological signals such as the alpha-rhythm, its application may increase the overall quality of a dataset, particularly when strongly affected by motion. This was verified by analysing the EEG traces, power spectra density and the topographic distribution from two healthy subjects. We also have verified that the correction by REEG-MEGAS improves higher frequency artefact correction by decreasing the power of Gradient Artefact harmonics. Our method showed promising results for decreasing the power of artefacts for frequencies up to 250 Hz. Additionally, REEG-MEGAS is a hybrid framework that can be implemented for real time prospective motion correction of EEG and fMRI data. Among other EEG-fMRI applications, the approach described here may benefit applications such as EEG-fMRI neurofeedback and brain computer interface, which strongly rely on the prospective acquisition and application of motion artefact removal.

Removal of physiological artifacts from simultaneous EEG and fMRI recordings

OBJECTIVE

Simultaneous recording of the electroencephalogram (EEG) and functional magnetic resonance imaging (fMRI) allows a combination of eletrophysiological and haemodynamic information to be used to form a more complete picture of cerebral dynamics. However, EEG recorded within the MRI scanner is contaminated by both imaging artifacts and physiological artifacts. The majority of the techniques used to pre-process such EEG focus on removal of the imaging and balistocardiogram artifacts, with some success, but don't remove all other physiological artifacts.

METHODS

We propose a new offline EEG artifact removal method based upon a combination of independent component analysis and fMRI-based head movement estimation to aid the removal of physiological artifacts from EEG recorded during EEG-fMRI recordings. Our method makes novel use of head movement trajectories estimated from the fMRI recording in order to assist with identifying physiological artifacts in the EEG and is designed to be used after removal of the fMRI imaging artifact from the EEG.

RESULTS

We evaluate our method on EEG recorded during a joint EEG-fMRI session from healthy adult participants. Our method significantly reduces the influence of all types of physiological artifacts on the EEG. We also compare our method with a state-of-the-art physiological artifact removal method and demonstrate superior performance removing physiological artifacts.

CONCLUSIONS

Our proposed method is able to remove significantly more physiological artifact components from the EEG, recorded during a joint EEG-fMRI session, than other state-of-the-art methods.

SIGNIFICANCE

Our proposed method represents a marked improvement over current processing pipelines for removing physiological noise from EEG recorded during a joint EEG-fMRI session.

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.

Retrospective Camera-Based Respiratory Gating in Clinical Whole-Heart 4D Flow MRI

Background

Respiratory gating is generally recommended in 4D flow MRI of the heart to avoid blurring and motion artifacts. Recently, a novel automated contact‐less camera‐based respiratory motion sensor has been introduced.

Purpose

To compare camera‐based respiratory gating (CAM) with liver‐lung‐navigator‐based gating (NAV) and no gating (NO) for whole‐heart 4D flow MRI.

Study Type

Retrospective.

Subjects

Thirty two patients with a spectrum of cardiovascular diseases.

Field Strength/Sequence

A 3T, 3D‐cine spoiled‐gradient‐echo‐T1‐weighted‐sequence with flow‐encoding in three spatial directions.

Assessment

Respiratory phases were derived and compared against each other by cross‐correlation. Three radiologists/cardiologist scored images reconstructed with camera‐based, navigator‐based, and no respiratory gating with a 4‐point Likert scale (qualitative analysis). Quantitative image quality analysis, in form of signal‐to‐noise ratio (SNR) and liver‐lung‐edge (LLE) for sharpness and quantitative flow analysis of the valves were performed semi‐automatically.

Statistical Tests

One‐way repeated measured analysis of variance (ANOVA) with Wilks's lambda testing and follow‐up pairwise comparisons. Significance level of P ≤ 0.05. Krippendorff's‐alpha‐test for inter‐rater reliability.

Results

The respiratory signal analysis revealed that CAM and NAV phases were highly correlated (C = 0.93 ± 0.09, P < 0.01). Image scoring showed poor inter‐rater reliability and no significant differences were observed (P ≥ 0.16). The image quality comparison showed that NAV and CAM were superior to NO with higher SNR (P = 0.02) and smaller LLE (P < 0.01). The quantitative flow analysis showed significant differences between the three respiratory‐gated reconstructions in the tricuspid and pulmonary valves (P ≤ 0.05), but not in the mitral and aortic valves (P > 0.05). Pairwise comparisons showed that reconstructions without respiratory gating were different in flow measurements to either CAM or NAV or both, but no differences were found between CAM and NAV reconstructions.

Data Conclusion

Camera‐based respiratory gating performed as well as conventional liver‐lung‐navigator‐based respiratory gating. Quantitative image quality analysis showed that both techniques were equivalent and superior to no‐gating‐reconstructions. Quantitative flow analysis revealed local flow differences (tricuspid/pulmonary valves) in images of no‐gating‐reconstructions, but no differences were found between images reconstructed with camera‐based and navigator‐based respiratory gating.

Level of Evidence

3

Technical Efficacy

Stage 2

Schrödinger filtering: a precise EEG despiking technique for EEG-fMRI gradient artifact

In EEG data acquired in the presence of fMRI, gradient-related spike artifacts contaminate the signal following the common preprocessing step of average artifact subtraction. Spike artifacts compromise EEG data quality since they overlap with the EEG signal in frequency, thereby confounding frequency-based inferences on activity. As well, spike artifacts can inflate or deflate correlations among time series, thereby confounding inferences on functional connectivity. We present Schrödinger filtering, which uses the Schrödinger equation to decompose the spike-containing input. The basis functions of the decomposition are localized and pulse-shaped, and selectively capture the various input peaks, with the spike components clustered at the beginning of the spectrum. Schrödinger filtering automatically subtracts the spike components from the data. On real and simulated data, we show that Schrödinger filtering (1) simultaneously accomplishes high spike removal and high signal preservation without affecting evoked activity, and (2) reduces spurious pairwise correlations in spontaneous activity. In these regards, Schrödinger filtering was significantly better than three other despiking techniques: median filtering, amplitude thresholding, and wavelet denoising. These results encourage the use of Schrödinger filtering in future EEG-fMRI pipelines, as well as in other spike-related applications (e.g., fMRI motion artifact removal or action potential extraction).

Concurrent electrophysiological and hemodynamic measurements of evoked neural oscillations in human visual cortex using sparsely interleaved fast fMRI and EEG

Electroencephalography (EEG) concurrently collected with functional magnetic resonance imaging (fMRI) is heavily distorted by the repetitive gradient coil switching during the fMRI acquisition. The performance of the typical template-based gradient artifact suppression method can be suboptimal because the artifact changes over time. Gradient artifact residuals also impede the subsequent suppression of ballistocardiography artifacts. Here we propose recording continuous EEG with temporally sparse fast fMRI (fast fMRI-EEG) to minimize the EEG artifacts caused by MRI gradient coil switching without significantly compromising the field-of-view and spatiotemporal resolution of fMRI. Using simultaneous multi-slice inverse imaging to achieve whole-brain fMRI with isotropic 5-mm resolution in 0.1 ​s, and performing these acquisitions once every 2 ​s, we have 95% of the duty cycle available to record EEG with substantially less gradient artifact. We found that the standard deviation of EEG signals over the entire acquisition period in fast fMRI-EEG was reduced to 54% of that in conventional concurrent echo-planar imaging (EPI) and EEG recordings (EPI-EEG) across participants. When measuring 15-Hz steady-state visual evoked potentials (SSVEPs), the baseline-normalized oscillatory neural response in fast fMRI-EEG was 2.5-fold of that in EPI-EEG. The functional MRI responses associated with the SSVEP delineated by EPI and fast fMRI were similar in the spatial distribution, the elicited waveform, and detection power. Sparsely interleaved fast fMRI-EEG provides high-quality EEG without substantially compromising the quality of fMRI in evoked response measurements, and has the potential utility for applications where the onset of the target stimulus cannot be precisely determined, such as epilepsy.

Prospective motion correction of fMRI: Improving the quality of resting state data affected by large head motion

The quality of functional MRI (fMRI) data is affected by head motion. It has been shown that fMRI data quality can be improved by prospectively updating the gradients and radio-frequency pulses in response to head motion during image acquisition by using an MR-compatible optical tracking system (prospective motion correction, or PMC). Recent studies showed that PMC improves the temporal Signal to Noise Ratio (tSNR) of resting state fMRI data (rs-fMRI) acquired from subjects not moving intentionally. Besides that, the time courses of Independent Components (ICs), resulting from Independent Component Analysis (ICA), were found to present significant temporal correlation with the motion parameters recorded by the camera. However, the benefits of applying PMC for improving the quality of rs-fMRI acquired under large head movements and its effects on resting state networks (RSN) and connectivity matrices are still unknown. In this study, subjects were instructed to cross their legs at will while rs-fMRI data with and without PMC were acquired, which generated head motion velocities ranging from 4 to 30 ​mm/s. We also acquired fMRI data without intentional motion. Independent component analysis of rs-fMRI was performed to evaluate IC maps and time courses of RSNs. We also calculated the temporal correlation among different brain regions and generated connectivity matrices for the different motion and PMC conditions. In our results we verified that the crossing leg movements reduced the tSNR of sessions without and with PMC by 45 and 20%, respectively, when compared to sessions without intentional movements. We have verified an interaction between head motion speed and PMC status, showing stronger attenuation of tSNR for acquisitions without PMC than for those with PMC. Additionally, the spatial definition of major RSNs, such as default mode, visual, left and right central executive networks, was improved when PMC was enabled. Furthermore, motion altered IC-time courses by decreasing power at low frequencies and increasing power at higher frequencies (typically associated with artefacts). PMC partially reversed these alterations of the power spectra. Finally, we showed that PMC provides temporal correlation matrices for data acquired under motion conditions more comparable to those obtained by fMRI sessions where subjects were instructed not to move.

Changing concepts in presurgical assessment for epilepsy surgery

Candidates for epilepsy surgery must undergo presurgical evaluation to establish whether and how surgical treatment can stop seizures without causing neurological deficits. Various techniques, including MRI, PET, single-photon emission CT, video-EEG, magnetoencephalography and invasive EEG, aim to identify the diseased brain tissue and the involved network. Recent technical and methodological developments, encompassing both advances in existing techniques and new combinations of technologies, are enhancing the ability to define the optimal resection strategy. Multimodal interpretation and predictive computer models are expected to aid surgical planning and patient counselling, and multimodal intraoperative guidance is likely to increase surgical precision. In this Review, we discuss how the knowledge derived from these new approaches is challenging our way of thinking about surgery to stop focal seizures. In particular, we highlight the importance of looking beyond the EEG seizure onset zone and considering focal epilepsy as a brain network disease in which long-range connections need to be taken into account. We also explore how new diagnostic techniques are revealing essential information in the brain that was previously hidden from view.

Direct modelling of gradient artifacts for EEG-fMRI denoising and motion tracking

OBJECTIVE

Simultaneous electroencephalography and functional magnetic resonance imaging recording (EEG-fMRI) has been widely used in neuroscientific and clinical research. The artifacts in the recorded EEG resulting from rapidly switching magnetic field gradients are usually corrected by average-artifact subtraction (AAS) due to their repetitive nature. But the performance of AAS is often disrupted by altered artifact waveforms across epochs, notably due to head motion.

APPROACH

Here, a method is proposed to make use of the known MR sequence gradient waveforms for a direct modelling of gradient artifacts. After accounting for filtering effects on the gradient artifacts, a continuous modulation of the gradient waveforms superimposed on the EEG signal is obtained.

MAIN RESULTS

Although a moving AAS template can adjust to slow drifts in gradient artifact variation, it fails to adapt to abrupt motion, resulting in residual noise. We demonstrate how this modelling approach can reduce motion-affected gradient artifacts without distorting the underlying neuronal signals. Moreover, the method provides useful head motion information highly correlated with motion tracked by an optical camera.

SIGNIFICANCE

Our work provides a novel way to improve gradient artifact removal in EEG-fMRI, and shows a potential to detect head motion without requiring additional hardware.

Conductive gel bridge sensor for motion tracking in simultaneous EEG-fMRI recordings

EEG-fMRI allows the localization of the hemodynamic correlates of neural activity and has been shown to be useful as a diagnostic tool in pre-surgical evaluation of refractory epilepsy. However, EEG recordings may be highly contaminated by artifacts induced by movements inside the magnetic field thus rendering the scan difficult for interpretation. Existing methods for motion correction require additional equipment or hardware modification. We introduce a simple method for motion artifact detection, the conductive gel bridge sensor (CGBS), easily applicable using the standard setup. We report examples of CGBS use in two patients with epilepsy and demonstrate the method's ability to successfully differentiate between epochs of brain activity and those of movement.

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.

Annual Research Review: Not just a small adult brain: understanding later neurodevelopment through imaging the neonatal brain

BACKGROUND

There has been a recent proliferation in neuroimaging research focusing on brain development in the prenatal, neonatal and very early childhood brain. Early brain injury and preterm birth are associated with increased risk of neurodevelopmental disorders, indicating the importance of this early period for later outcome.

SCOPE AND METHODOLOGY

Although using a wide range of different methodologies and investigating diverse samples, the common aim of many of these studies has been to both track normative development and investigate deviations in this development to predict behavioural, cognitive and neurological function in childhood. Here we review structural and functional neuroimaging studies investigating the developing brain. We focus on practical and technical complexities of studying this early age range and discuss how neuroimaging techniques have been successfully applied to investigate later neurodevelopmental outcome.

CONCLUSIONS

Neuroimaging markers of later outcome still have surprisingly low predictive power and their specificity to individual neurodevelopmental disorders is still under question. However, the field is still young, and substantial challenges to both acquiring and modeling neonatal data are being met.

Exploring the origins of EEG motion artefacts during simultaneous fMRI acquisition: Implications for motion artefact correction

Motion artefacts (MAs) are induced within EEG data collected simultaneously with fMRI when the subject's head rotates relative to the magnetic field. The effects of these artefacts have generally been ameliorated by removing periods of data during which large artefact voltages appear in the EEG traces. However, even when combined with other standard post-processing methods, this strategy does not remove smaller MAs which can dominate the neuronal signals of interest. A number of methods are therefore being developed to characterise the MA by measuring reference signals and then using these in artefact correction. These methods generally assume that the head and EEG cap, plus any attached sensors, form a rigid body which can be characterised by a standard set of six motion parameters. Here we investigate the motion of the head/EEG cap system to provide a better understanding of MAs. We focus on the reference layer artefact subtraction (RLAS) approach, as this allows measurement of a separate reference signal for each electrode that is being used to measure brain activity.

Through a series of experiments on phantoms and subjects, we find that movement of the EEG cap relative to the phantom and skin on the forehead is relatively small and that this non-rigid body movement does not appear to cause considerable discrepancy in artefacts between the scalp and reference signals. However, differences in the amplitude of these signals is observed which may be due to differences in geometry of the system from which the reference signals are measured compared with the brain signals. In addition, we find that there is non-rigid body movement of the skull and skin which produces an additional MA component for a head shake, which is not present for a head nod. This results in a large discrepancy in the amplitude and temporal profile of the MA measured on the scalp and reference layer, reducing the efficacy of MA correction based on the reference signals.

Together our data suggest that the efficacy of the correction of MA using any reference-based system is likely to differ for different types of head movement with head shake being the hardest to correct. This provides new information to inform the development of hardware and post-processing methods for removing MAs from EEG data acquired simultaneously with fMRI data.

EEG-Informed fMRI: A Review of Data Analysis Methods

The simultaneous acquisition of electroencephalography (EEG) with functional magnetic resonance imaging (fMRI) is a very promising non-invasive technique for the study of human brain function. Despite continuous improvements, it remains a challenging technique, and a standard methodology for data analysis is yet to be established. Here we review the methodologies that are currently available to address the challenges at each step of the data analysis pipeline. We start by surveying methods for pre-processing both EEG and fMRI data. On the EEG side, we focus on the correction for several MR-induced artifacts, particularly the gradient and pulse artifacts, as well as other sources of EEG artifacts. On the fMRI side, we consider image artifacts induced by the presence of EEG hardware inside the MR scanner, and the contamination of the fMRI signal by physiological noise of non-neuronal origin, including a review of several approaches to model and remove it. We then provide an overview of the approaches specifically employed for the integration of EEG and fMRI when using EEG to predict the blood oxygenation level dependent (BOLD) fMRI signal, the so-called EEG-informed fMRI integration strategy, the most commonly used strategy in EEG-fMRI research. Finally, we systematically review methods used for the extraction of EEG features reflecting neuronal phenomena of interest.

Exploring the advantages of multiband fMRI with simultaneous EEG to investigate coupling between gamma frequency neural activity and the BOLD response in humans

We established an optimal combination of EEG recording during sparse multiband (MB) fMRI that preserves high-resolution, whole-brain fMRI coverage while enabling broad-band EEG recordings which are uncorrupted by MRI gradient artefacts (GAs). We first determined the safety of simultaneous EEG recording during MB fMRI. Application of MB factor = 4 produced <1°C peak heating of electrode/hardware during 20 min of GE-EPI data acquisition. However, higher SAR sequences require specific safety testing, with greater heating observed using PCASL with MB factor = 4. Heating was greatest in the electrocardiogram channel, likely due to it possessing longest lead length. We investigated the effect of MB factor on the temporal signal-to-noise ratio for a range of GE-EPI sequences (varying MB factor and temporal interval between slice acquisitions). We found that, for our experimental purpose, the optimal acquisition was achieved with MB factor = 3, 3mm isotropic voxels, and 33 slices providing whole head coverage. This sequence afforded a 2.25 s duration quiet period (without GAs) in every 3 s TR. Using this sequence, we demonstrated the ability to record gamma frequency (55-80 Hz) EEG oscillations, in response to right index finger abduction, that are usually obscured by GAs during continuous fMRI data acquisition. In this novel application of EEG-MB fMRI to a motor task, we observed a positive correlation between gamma and BOLD responses in bilateral motor regions. These findings support and extend previous work regarding coupling between neural and hemodynamic measures of brain activity in humans and showcase the utility of EEG-MB fMRI for future investigations.

DeepIED: An epileptic discharge detector for EEG-fMRI based on deep learning

Presurgical evaluation that can precisely delineate the epileptogenic zone (EZ) is one important step for successful surgical resection treatment of refractory epilepsy patients. The noninvasive EEG-fMRI recording technique combined with general linear model (GLM) analysis is considered an important tool for estimating the EZ. However, the manual marking of interictal epileptic discharges (IEDs) needed in this analysis is challenging and time-consuming because the quality of the EEG recorded inside the scanner is greatly deteriorated compared to the usual EEG obtained outside the scanner. This is one of main impediments to the widespread use of EEG-fMRI in epilepsy. We propose a deep learning based semi-automatic IED detector that can find the candidate IEDs in the EEG recorded inside the scanner which resemble sample IEDs marked in the EEG recorded outside the scanner. The manual marking burden is greatly reduced as the expert need only edit candidate IEDs. The model is trained on data from 30 patients. Validation of IEDs detection accuracy on another 37 consecutive patients shows our method can improve the median sensitivity from 50.0% for the previously proposed template-based method to 84.2%, with false positive rate as 5 events/min. Reproducibility validation on 15 patients is applied to evaluate if our method can produce similar hemodynamic response maps compared with the manual marking ground truth results. We explore the concordance between the maximum hemodynamic response and the intracerebral EEG defined EZ and find that both methods produce similar percentage of concordance (76.9%, 10 out of 13 patients, electrode was absent in the maximum hemodynamic response in two patients). This tool will make EEG-fMRI analysis more practical for clinical usage.

•

A deep learning based epileptic discharge detector for EEG-fMRI is proposed.

•

The burden of manually marking epileptic discharges is greatly reduced.

•

Our method can produce similar EEG-fMRI results compared with traditional method.

Two-Dimensional Temporal Clustering Analysis for Patients with Epilepsy: Detecting Epilepsy-Related Information in EEG-fMRI Concordant, Discordant and Spike-Less Patients

EEG acquired simultaneously with fMRI (EEG-fMRI) is a multimodal method that has shown promise in mapping the seizure onset zone in patients with focal epilepsy. However, there are many instances when this method is unsuccessful or not applicable, and other data driven fMRI methods may be utilized. One such method is the two-dimensional temporal clustering analysis (2dTCA). In this study we compared the classic EEG-fMRI and 2dTCA performance in mapping regions related to the seizure onset region in 18 focal epilepsy patients (12 presenting interictal epileptiform discharges (IEDs), during EEG-fMRI acquisition) with Engel I or II surgical outcome. Activation maps of both 2dTCA timing outputs (positive and negative histograms) and EEG detected IEDs were computed and compared to the region of epilepsy surgical resection. Patients were evaluated in three categories based on frequency of EEG detected spiking during the MRI. EEG-fMRI maps were concordant to the epilepsy region in 5/12 subjects, four with frequent IEDs on EEG. The 2dTCA was successful in mapping 13/18 patients including 3/6 with no IEDs detected (10/12 with IEDs detected). The epilepsy-related activities were successfully mapped by both methods in only 4/12 patients. This work suggests that the epilepsy-related information detected by each method may be different: while EEG-fMRI is more accurate in patients with high rather than lower numbers of EEG detected IEDs; 2dTCA can be useful in evaluating patients even when no concurrent EEG spikes are detected or EEG-fMRI is not effective. Therefore, our results support that 2dTCA might be an alternative for mapping epilepsy-related BOLD activity in negative EEG-fMRI (6/7 patients) and spike-less patients.

Minimizing noise in pediatric task-based functional MRI; Adolescents with developmental disabilities and typical development

Functional Magnetic Resonance Imaging (fMRI) represents a powerful tool with which to examine brain functioning and development in typically developing pediatric groups as well as children and adolescents with clinical disorders. However, fMRI data can be highly susceptible to misinterpretation due to the effects of excessive levels of noise, often related to head motion. Imaging children, especially with developmental disorders, requires extra considerations related to hyperactivity, anxiety and the ability to perform and maintain attention to the fMRI paradigm. We discuss a number of methods that can be employed to minimize noise, in particular movement-related noise. To this end we focus on strategies prior to, during and following the data acquisition phase employed primarily within our own laboratory. We discuss the impact of factors such as experimental design, screening of potential participants and pre-scan training on head motion in our adolescents with developmental disorders and typical development. We make some suggestions that may minimize noise during data acquisition itself and finally we briefly discuss some current processing techniques that may help to identify and remove noise in the data. Many advances have been made in the field of pediatric imaging, particularly with regard to research involving children with developmental disorders. Mindfulness of issues such as those discussed here will ensure continued progress and greater consistency across studies.

Prospective motion correction in functional MRI

Due to the intrinsic low sensitivity of BOLD-fMRI long scanning is required. Subject motion during fMRI scans reduces statistical significance of the activation maps and increases the prevalence of false activations. Motion correction is therefore an essential tool for a successful fMRI data analysis. Retrospective motion correction techniques are now commonplace and are incorporated into a wide range of fMRI analysis toolboxes. These techniques are advantageous due to robustness, sequence independence and have minimal impact on the fMRI study setup. Retrospective techniques however, do not provide an accurate intra-volume correction, nor can these techniques correct for the spin-history effects. The application of prospective motion correction in fMRI appears to be effective in reducing false positives and increasing sensitivity when compared to retrospective techniques, particularly in the cases of substantial motion. Especially advantageous in this regard is the combination of prospective motion correction with dynamic distortion correction. Nevertheless, none of the recent methods are able to recover activations in presence of motion that are comparable to no-motion conditions, which motivates further research in the area of adaptive dynamic imaging.