NREM sleep brain networks modulate cognitive recovery from sleep deprivation.

Decrease in cognitive performance after sleep deprivation followed by recovery after sleep suggests its key role, and especially non-rapid eye movement (NREM) sleep, in the maintenance of cognition. It remains unknown whether brain network reorganization in NREM sleep stages N2 and N3 can uniquely be mapped onto individual differences in cognitive performance after a recovery nap following sleep deprivation. Using resting state functional magnetic resonance imaging (fMRI), we quantified the integration and segregation of brain networks during NREM sleep stages N2 and N3 while participants took a 1-hour nap following 24-hour sleep deprivation, compared to well-rested wakefulness. Here, we advance a new analytic framework called the hierarchical segregation index (HSI) to quantify network segregation across spatial scales, from whole-brain to the voxel level, by identifying spatio-temporally overlapping large-scale networks and the corresponding voxel-to-region hierarchy. Our results show that network segregation increased in the default mode, dorsal attention and somatomotor networks during NREM sleep compared to wakefulness. Segregation within the visual, limbic, and executive control networks exhibited N2 versus N3 sleep-specific voxel-level patterns. More segregation during N3 was associated with worse recovery of working memory, executive attention, and psychomotor vigilance after the nap. The level of spatial resolution of network segregation varied among brain regions and was associated with the recovery of performance in distinct cognitive tasks. We demonstrated the sensitivity and reliability of voxel-level HSI to provide key insights into within-region variation, suggesting a mechanistic understanding of how NREM sleep replenishes cognition after sleep deprivation.

EEG/MEG source imaging of deep brain activity within the maximum entropy on the mean framework: Simulations and validation in epilepsy.

Electro/Magneto-EncephaloGraphy (EEG/MEG) source imaging (EMSI) of epileptic activity from deep generators is often challenging due to the higher sensitivity of EEG/MEG to superficial regions and to the spatial configuration of subcortical structures. We previously demonstrated the ability of the coherent Maximum Entropy on the Mean (cMEM) method to accurately localize the superficial cortical generators and their spatial extent. Here, we propose a depth-weighted adaptation of cMEM to localize deep generators more accurately. These methods were evaluated using realistic MEG/high-density EEG (HD-EEG) simulations of epileptic activity and actual MEG/HD-EEG recordings from patients with focal epilepsy. We incorporated depth-weighting within the MEM framework to compensate for its preference for superficial generators. We also included a mesh of both hippocampi, as an additional deep structure in the source model. We generated 5400 realistic simulations of interictal epileptic discharges for MEG and HD-EEG involving a wide range of spatial extents and signal-to-noise ratio (SNR) levels, before investigating EMSI on clinical HD-EEG in 16 patients and MEG in 14 patients. Clinical interictal epileptic discharges were marked by visual inspection. We applied three EMSI methods: cMEM, depth-weighted cMEM and depth-weighted minimum norm estimate (MNE). The ground truth was defined as the true simulated generator or as a drawn region based on clinical information available for patients. For deep sources, depth-weighted cMEM improved the localization when compared to cMEM and depth-weighted MNE, whereas depth-weighted cMEM did not deteriorate localization accuracy for superficial regions. For patients' data, we observed improvement in localization for deep sources, especially for the patients with mesial temporal epilepsy, for which cMEM failed to reconstruct the initial generator in the hippocampus. Depth weighting was more crucial for MEG (gradiometers) than for HD-EEG. Similar findings were found when considering depth weighting for the wavelet extension of MEM. In conclusion, depth-weighted cMEM improved the localization of deep sources without or with minimal deterioration of the localization of the superficial sources. This was demonstrated using extensive simulations with MEG and HD-EEG and clinical MEG and HD-EEG for epilepsy patients.

A spatial perturbation framework to validate implantation of the epileptogenic zone.

Stereo-electroencephalography (SEEG) is the gold standard to delineate surgical targets in focal drug-resistant epilepsy. SEEG uses electrodes placed directly into the brain to identify the seizure-onset zone (SOZ). However, its major constraint is limited brain coverage, potentially leading to misidentification of the 'true' SOZ. Here, we propose a framework to assess adequate SEEG sampling by coupling epileptic biomarkers with their spatial distribution and measuring the system's response to a perturbation of this coupling. We demonstrate that the system's response is strongest in well-sampled patients when virtually removing the measured SOZ. We then introduce the spatial perturbation map, a tool that enables qualitative assessment of the implantation coverage. Probability modelling reveals a higher likelihood of well-implanted SOZs in seizure-free patients or non-seizure free patients with incomplete SOZ resections, compared to non-seizure-free patients with complete resections. This highlights the framework's value in sparing patients from unsuccessful surgeries resulting from poor SEEG coverage.

Systematic review of seizure-onset patterns in stereo-electroencephalography: Current state and future directions.

OBJECTIVE: Increasing evidence suggests that the seizure-onset pattern (SOP) in stereo-electroencephalography (SEEG) is important for localizing the "true" seizure onset. Specifically, SOPs with low-voltage fast activity (LVFA) are associated with seizure-free outcome (Engel I). However, several classifications and various terms corresponding to the same pattern have been reported, challenging its use in clinical practice. METHOD: Following the Preferred Reporting Items of Systematic reviews and Meta-Analyses (PRISMA) guideline, we performed a systematic review of studies describing SOPs along with accompanying figures depicting the reported SOP in SEEG. RESULTS: Of 1799 studies, 22 met the selection criteria. Among the various SOPs, we observed that the terminology for low frequency periodic spikes exhibited the most variability, whereas LVFA is the most frequently used term of this pattern. Some SOP terms were inconsistent with standard EEG terminology. Finally, there was a significant but weak association between presence of LVFA and seizure-free outcome. CONCLUSION: Divergent terms were used to describe the same SOPs and some of these terms showed inconsistencies with the standard EEG terminology. Additionally, our results confirmed the link between patterns with LVFA and seizure-free outcomes. However, this association was not strong. SIGNIFICANCE: These results underline the need for standardization of SEEG terminology.

Consistency of electrical source imaging in presurgical evaluation of epilepsy across different vigilance states.

OBJECTIVE: The use of electrical source imaging (ESI) in assessing the source of interictal epileptic discharges (IEDs) is gaining increasing popularity in presurgical work-up of patients with drug-resistant focal epilepsy. While vigilance affects the ability to locate IEDs and identify the epileptogenic zone, we know little about its impact on ESI. METHODS: We studied overnight high-density electroencephalography recordings in focal drug-resistant epilepsy. IEDs were marked visually in each vigilance state, and examined in the sensor and source space. ESIs were calculated and compared between all vigilance states and the clinical ground truth. Two conditions were considered within each vigilance state, an unequalized and an equalized number of IEDs. RESULTS: The number, amplitude, and duration of IEDs were affected by the vigilance state, with N3 sleep presenting the highest number, amplitude, and duration for both conditions (P < 0.001), while signal-to-noise ratio only differed in the unequalized condition (P < 0.001). The vigilance state did not affect channel involvement (P > 0.05). ESI maps showed no differences in distance, quality, extent, or maxima distances compared to the clinical ground truth for both conditions (P > 0.05). Only when an absolute reference (wakefulness) was used, the channel involvement (P < 0.05) and ESI source extent (P < 0.01) were impacted during rapid-eye-movement (REM) sleep. Clustering of amplitude-sensitive and -insensitive ESI maps pointed to amplitude rather than the spatial profile as the driver (P < 0.05). INTERPRETATION: IED ESI results are stable across vigilance states, including REM sleep, if controlled for amplitude and IED number. ESI is thus stable and invariant to the vigilance state.

Validating MEG source imaging of resting state oscillatory patterns with an intracranial EEG atlas.

BACKGROUND: Magnetoencephalography (MEG) is a widely used non-invasive tool to estimate brain activity with high temporal resolution. However, due to the ill-posed nature of the MEG source imaging (MSI) problem, the ability of MSI to identify accurately underlying brain sources along the cortical surface is still uncertain and requires validation. METHOD: We validated the ability of MSI to estimate the background resting state activity of 45 healthy participants by comparing it to the intracranial EEG (iEEG) atlas (https://mni-open-ieegatlas. RESEARCH: mcgill.ca/). First, we applied wavelet-based Maximum Entropy on the Mean (wMEM) as an MSI technique. Next, we converted MEG source maps into intracranial space by applying a forward model to the MEG-reconstructed source maps, and estimated virtual iEEG (ViEEG) potentials on each iEEG channel location; we finally quantitatively compared those with actual iEEG signals from the atlas for 38 regions of interest in the canonical frequency bands. RESULTS: The MEG spectra were more accurately estimated in the lateral regions compared to the medial regions. The regions with higher amplitude in the ViEEG than in the iEEG were more accurately recovered. In the deep regions, MEG-estimated amplitudes were largely underestimated and the spectra were poorly recovered. Overall, our wMEM results were similar to those obtained with minimum norm or beamformer source localization. Moreover, the MEG largely overestimated oscillatory peaks in the alpha band, especially in the anterior and deep regions. This is possibly due to higher phase synchronization of alpha oscillations over extended regions, exceeding the spatial sensitivity of iEEG but detected by MEG. Importantly, we found that MEG-estimated spectra were more comparable to spectra from the iEEG atlas after the aperiodic components were removed. CONCLUSION: This study identifies brain regions and frequencies for which MEG source analysis is likely to be reliable, a promising step towards resolving the uncertainty in recovering intracerebral activity from non-invasive MEG studies.

Hierarchical Bayesian modeling of the relationship between task-related hemodynamic responses and cortical excitability.

Investigating the relationship between task-related hemodynamic responses and cortical excitability is challenging because it requires simultaneous measurement of hemodynamic responses while applying noninvasive brain stimulation. Moreover, cortical excitability and task-related hemodynamic responses are both associated with inter-/intra-subject variability. To reliably assess such a relationship, we applied hierarchical Bayesian modeling. This study involved 16 healthy subjects who underwent simultaneous Paired Associative Stimulation (PAS10, PAS25, Sham) while monitoring brain activity using functional Near-Infrared Spectroscopy (fNIRS), targeting the primary motor cortex (M1). Cortical excitability was measured by Motor Evoked Potentials (MEPs), and the motor task-related hemodynamic responses were measured using fNIRS 3D reconstructions. We constructed three models to investigate: (1) PAS effects on the M1 excitability, (2) PAS effects on fNIRS hemodynamic responses to a finger tapping task, and (3) the correlation between PAS effects on M1 excitability and PAS effects on task-related hemodynamic responses. Significant increase in cortical excitability was found following PAS25, whereas a small reduction of the cortical excitability was shown after PAS10 and a subtle increase occurred after sham. Both HbO and HbR absolute amplitudes increased after PAS25 and decreased after PAS10. The probability of the positive correlation between modulation of cortical excitability and hemodynamic activity was 0.77 for HbO and 0.79 for HbR. We demonstrated that PAS stimulation modulates task-related cortical hemodynamic responses in addition to M1 excitability. Moreover, the positive correlation between PAS modulations of excitability and hemodynamics brought insight into understanding the fundamental properties of cortical function and cortical excitability.

Clinical Yield of Electromagnetic Source Imaging and Hemodynamic Responses in Epilepsy: Validation With Intracerebral Data.

OBJECTIVES: Accurate delineation of the seizure-onset zone (SOZ) in focal drug-resistant epilepsy often requires stereo-electroencephalography (SEEG) recordings. We aimed at: (1) proposing a truly objective and quantitative comparison between electro-encephalography/magnetoencephalography (EEG/MEG) source-imaging (EMSI), EEG/functional MRI (EEG/fMRI) responses for similar spikes with primary-irritative zone (PIZ) and SOZ defined by SEEG and (2) evaluating the value of EMSI and EEG/fMRI to predict postsurgical outcome. METHODS: We identified patients with drug-resistant epilepsy who underwent EEG/MEG, EEG/fMRI, and subsequent SEEG at the Epilepsy Service from the Montreal Neurological Institute and Hospital. We quantified multimodal concordance within the SEEG channel-space, as spatial overlap with PIZ/SOZ and distances to the Spike-onset, Spike-maximum-amplitude and Seizure-core intracerebral channels, by applying a new methodology consisting of converting EMSI results into SEEG electrical potentials (EMSIe-SEEG) and projecting the most significant fMRI response on the SEEG channels (fMRIp-SEEG). Spatial overlaps with PIZ/SOZ (AUCPIZ, AUCSOZ) were assessed by using the area under the receiver operating characteristic curve (AUC). Here, AUC represents the probability that a randomly picked active contact exhibited higher amplitude when located inside the spatial reference than outside. RESULTS: Seventeen patients were included. Mean spatial overlaps with the primary-irritative zone and seizure-onset zone were 0.71 and 0.65 for EMSIe-SEEG, and 0.57 and 0.62 for fMRIp-SEEG. Good EMSIe-SEEG  spatial overlap with the primary-irritative zone was associated with smaller distance from the maximum EMSIe-SEEG contact to the Spike-maximum-amplitude channel (median distance 14 mm). Conversely, good fMRIp-SEEG spatial overlap with the seizure-onset zone was associated with smaller distances from the maximum  fMRIp-SEEG contact to the Spike-onset and Seizure-core channels (median distances 10 mm and 5mm respectively). Surgical outcomes were correctly predicted by EEG/MEG in 12/15 (80%) patients and EEG/fMRI in 6/11(54%) patients. CONCLUSIONS: Using a unique quantitative approach estimating EMSI and fMRI results in the reference SEEG channel-space, EEG/MEG and EEG/fMRI accurately localized the seizure-onset zone as well as the primary-irritative zone. Precisely, EEG/MEG more accurately localized the primary-irritative zone, whereas EEG/fMRI was more sensitive to the seizure-onset zone. Both neuro-imaging techniques provide complementary localization that can help guiding SEEG implantation and selecting good candidates for surgery.

Diffuse optical reconstructions of functional near infrared spectroscopy data using maximum entropy on the mean.

Functional near-infrared spectroscopy (fNIRS) measures the hemoglobin concentration changes associated with neuronal activity. Diffuse optical tomography (DOT) consists of reconstructing the optical density changes measured from scalp channels to the oxy-/deoxy-hemoglobin concentration changes within the cortical regions. In the present study, we adapted a nonlinear source localization method developed and validated in the context of Electro- and Magneto-Encephalography (EEG/MEG): the Maximum Entropy on the Mean (MEM), to solve the inverse problem of DOT reconstruction. We first introduced depth weighting strategy within the MEM framework for DOT reconstruction to avoid biasing the reconstruction results of DOT towards superficial regions. We also proposed a new initialization of the MEM model improving the temporal accuracy of the original MEM framework. To evaluate MEM performance and compare with widely used depth weighted Minimum Norm Estimate (MNE) inverse solution, we applied a realistic simulation scheme which contained 4000 simulations generated by 250 different seeds at different locations and 4 spatial extents ranging from 3 to 40[Formula: see text] along the cortical surface. Our results showed that overall MEM provided more accurate DOT reconstructions than MNE. Moreover, we found that MEM was remained particularly robust in low signal-to-noise ratio (SNR) conditions. The proposed method was further illustrated by comparing to functional Magnetic Resonance Imaging (fMRI) activation maps, on real data involving finger tapping tasks with two different montages. The results showed that MEM provided more accurate HbO and HbR reconstructions in spatial agreement with the main fMRI cluster, when compared to MNE.

An altered balance of integrated and segregated brain activity is a marker of cognitive deficits following sleep deprivation.

Sleep deprivation (SD) leads to impairments in cognitive function. Here, we tested the hypothesis that cognitive changes in the sleep-deprived brain can be explained by information processing within and between large-scale cortical networks. We acquired functional magnetic resonance imaging (fMRI) scans of 20 healthy volunteers during attention and executive tasks following a regular night of sleep, a night of SD, and a recovery nap containing nonrapid eye movement (NREM) sleep. Overall, SD was associated with increased cortex-wide functional integration, driven by a rise of integration within cortical networks. The ratio of within versus between network integration in the cortex increased further in the recovery nap, suggesting that prolonged wakefulness drives the cortex towards a state resembling sleep. This balance of integration and segregation in the sleep-deprived state was tightly associated with deficits in cognitive performance. This was a distinct and better marker of cognitive impairment than conventional indicators of homeostatic sleep pressure, as well as the pronounced thalamocortical connectivity changes that occurs towards falling asleep. Importantly, restoration of the balance between segregation and integration of cortical activity was also related to performance recovery after the nap, demonstrating a bidirectional effect. These results demonstrate that intra- and interindividual differences in cortical network integration and segregation during task performance may play a critical role in vulnerability to cognitive impairment in the sleep-deprived state.

Evaluation of a personalized functional near infra-red optical tomography workflow using maximum entropy on the mean.

In the present study, we proposed and evaluated a workflow of personalized near infra-red optical tomography (NIROT) using functional near-infrared spectroscopy (fNIRS) for spatiotemporal imaging of cortical hemodynamic fluctuations. The proposed workflow from fNIRS data acquisition to local 3D reconstruction consists of: (a) the personalized optimal montage maximizing fNIRS channel sensitivity to a predefined targeted brain region; (b) the optimized fNIRS data acquisition involving installation of optodes and digitalization of their positions using a neuronavigation system; and (c) the 3D local reconstruction using maximum entropy on the mean (MEM) to accurately estimate the location and spatial extent of fNIRS hemodynamic fluctuations along the cortical surface. The workflow was evaluated on finger-tapping fNIRS data acquired from 10 healthy subjects for whom we estimated the reconstructed NIROT spatiotemporal images and compared with functional magnetic resonance imaging (fMRI) results from the same individuals. Using the fMRI activation maps as our reference, we quantitatively compared the performance of two NIROT approaches, the MEM framework and the conventional minimum norm estimation (MNE) method. Quantitative comparisons were performed at both single subject and group-level. Overall, our results suggested that MEM provided better spatial accuracy than MNE, while both methods offered similar temporal accuracy when reconstructing oxygenated (HbO) and deoxygenated hemoglobin (HbR) concentration changes evoked by finger-tapping. Our proposed complete workflow was made available in the brainstorm fNIRS processing plugin-NIRSTORM, thus providing the opportunity for other researchers to further apply it to other tasks and on larger populations.

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.

How cerebral cortex protects itself from interictal spikes: The alpha/beta inhibition mechanism.

Interactions between interictal epileptiform discharges (IEDs) and distant cortical regions subserve potential effects on cognition of patients with focal epilepsy. We hypothesize that "healthy" brain areas at a distance from the epileptic focus may respond to the interference of IEDs by generating inhibitory alpha and beta oscillations. We predict that more prominent alpha-beta oscillations can be found in patients with less impaired neurocognitive profile. We performed a source imaging magnetoencephalography study, including 41 focal epilepsy patients: 21 with frontal lobe epilepsy (FLE) and 20 with mesial temporal lobe epilepsy. We investigated the effect of anterior (i.e., frontal and temporal) IEDs on the oscillatory pattern over posterior head regions. We compared cortical oscillations (5-80 Hz) temporally linked to 3,749 IEDs (1,945 frontal and 1,803 temporal) versus an equal number of IED-free segments. We correlated results from IED triggered oscillations to global neurocognitive performance. Only frontal IEDs triggered alpha-beta oscillations over posterior head regions. IEDs with higher amplitude triggered alpha-beta oscillations of higher magnitude. The intensity of posterior head region alpha-beta oscillations significantly correlated with a better neuropsychological profile. Our study demonstrated that cerebral cortex protects itself from IEDs with generation of inhibitory alpha-beta oscillations at distant cortical regions. The association of more prominent oscillations with a better cognitive status suggests that this mechanism might play a role in determining the cognitive resilience in patients with FLE.

Errata: Best practices for fNIRS publications.

[This corrects the article DOI: 10.1117/1.NPh.8.1.012101.].

Fast oscillations >40 Hz localize the epileptogenic zone: An electrical source imaging study using high-density electroencephalography.

OBJECTIVE: Fast Oscillations (FO) >40 Hz are a promising biomarker of the epileptogenic zone (EZ). Evidence using scalp electroencephalography (EEG) remains scarce. We assessed if electrical source imaging of FO using 256-channel high-density EEG (HD-EEG) is useful for EZ identification. METHODS: We analyzed HD-EEG recordings of 10 focal drug-resistant epilepsy patients with seizure-free postsurgical outcome. We marked FO candidate events at the time of epileptic spikes and verified them by screening for an isolated peak in the time-frequency plot. We performed electrical source imaging of spikes and FO within the Maximum Entropy of the Mean framework. Source localization maps were validated against the surgical cavity. RESULTS: We identified FO in five out of 10 patients who had a superficial or intermediate deep generator. The maximum of the FO maps was localized inside the cavity in all patients (100%). Analysis with a reduced electrode coverage using the 10-10 and 10-20 system showed a decreased localization accuracy of 60% and 40% respectively. CONCLUSIONS: FO recorded with HD-EEG localize the EZ. HD-EEG is better suited to detect and localize FO than conventional EEG approaches. SIGNIFICANCE: This study acts as proof-of-concept that FO localization using 256-channel HD-EEG is a viable marker of the EZ.

Best practices for fNIRS publications.

The application of functional near-infrared spectroscopy (fNIRS) in the neurosciences has been expanding over the last 40 years. Today, it is addressing a wide range of applications within different populations and utilizes a great variety of experimental paradigms. With the rapid growth and the diversification of research methods, some inconsistencies are appearing in the way in which methods are presented, which can make the interpretation and replication of studies unnecessarily challenging. The Society for Functional Near-Infrared Spectroscopy has thus been motivated to organize a representative (but not exhaustive) group of leaders in the field to build a consensus on the best practices for describing the methods utilized in fNIRS studies. Our paper has been designed to provide guidelines to help enhance the reliability, repeatability, and traceability of reported fNIRS studies and encourage best practices throughout the community. A checklist is provided to guide authors in the preparation of their manuscripts and to assist reviewers when evaluating fNIRS papers.

Cortical gradients of functional connectivity are robust to state-dependent changes following sleep deprivation.

Sleep deprivation leads to significant impairments in cognitive performance and changes to the interactions between large scale cortical networks, yet the hierarchical organization of cortical activity across states is still being explored. We used functional magnetic resonance imaging to assess activations and connectivity during cognitive tasks in 20 healthy young adults, during three states: (i) following a normal night of sleep, (ii) following 24hr of total sleep deprivation, and (iii) after a morning recovery nap. Situating cortical activity during cognitive tasks along hierarchical organizing gradients based upon similarity of functional connectivity patterns, we found that regional variations in task-activations were captured by an axis differentiating areas involved in executive control from default mode regions and paralimbic cortex. After global signal regression, the range of functional differentiation along this axis at baseline was significantly related to decline in working memory performance (2-back task) following sleep deprivation, as well as the extent of recovery in performance following a nap. The relative positions of cortical regions within gradients did not significantly change across states, except for a lesser differentiation of the visual system and increased coupling of the posterior cingulate cortex with executive control areas after sleep deprivation. This was despite a widespread increase in the magnitude of functional connectivity across the cortex following sleep deprivation. Cortical gradients of functional differentiation thus appear relatively insensitive to state-dependent changes following sleep deprivation and recovery, suggesting that there are no large-scale changes in cortical functional organization across vigilance states. Certain features of particular gradient axes may be informative for the extent of decline in performance on more complex tasks following sleep deprivation, and could be beneficial over traditional voxel- or parcel-based approaches in identifying realtionships between state-dependent brain activity and behavior.

Effects of Independent Component Analysis on Magnetoencephalography Source Localization in Pre-surgical Frontal Lobe Epilepsy Patients.

Objective: Magnetoencephalography source imaging (MSI) of interictal epileptiform discharges (IED) is a useful presurgical tool in the evaluation of drug-resistant frontal lobe epilepsy (FLE) patients. Yet, failures in MSI can arise related to artifacts and to interference of background activity. Independent component analysis (ICA) is a popular denoising procedure but its clinical application remains challenging, as the selection of multiple independent components (IC) is controversial, operator dependent, and time consuming. We evaluated whether selecting only one IC of interest based on its similarity with the average IED field improves MSI in FLE. Methods: MSI was performed with the equivalent current dipole (ECD) technique and two distributed magnetic source imaging (dMSI) approaches: minimum norm estimate (MNE) and coherent Maximum Entropy on the Mean (cMEM). MSI accuracy was evaluated under three conditions: (1) ICA of continuous data (Cont_ICA), (2) ICA at the time of IED (IED_ICA), and (3) without ICA (No_ICA). Localization performance was quantitatively measured as actual distance of the source maximum in relation to the focus (Dmin), and spatial dispersion (SD) for dMSI. Results: After ICA, ECD Dmin did not change significantly (p > 0.200). For both dMSI techniques, ICA application worsened the source localization accuracy. We observed a worsening of both MNE Dmin (p < 0.05, consistently) and MNE SD (p < 0.001, consistently) for both ICA approaches. A similar behaviour was observed for cMEM, for which, however, Cont_ICA seemed less detrimental. Conclusion: We demonstrated that a simplified ICA approach selecting one IC of interest in combination with distributed magnetic source imaging can be detrimental. More complex approaches may provide better results but would be rather difficult to apply in real-world clinical setting. In a broader perspective, caution should be taken in applying ICA for source localization of interictal activity. To ensure optimal and useful results, effort should focus on acquiring good quality data, minimizing artifacts, and determining optimal candidacy for MEG, rather than counting on data cleaning techniques.

Accuracy and spatial properties of distributed magnetic source imaging techniques in the investigation of focal epilepsy patients.

Source localization of interictal epileptiform discharges (IEDs) is clinically useful in the presurgical workup of epilepsy patients. We aimed to compare the performance of four different distributed magnetic source imaging (dMSI) approaches: Minimum norm estimate (MNE), dynamic statistical parametric mapping (dSPM), standardized low-resolution electromagnetic tomography (sLORETA), and coherent maximum entropy on the mean (cMEM). We also evaluated whether a simple average of maps obtained from multiple inverse solutions (Ave) can improve localization accuracy. We analyzed dMSI of 206 IEDs derived from magnetoencephalography recordings in 28 focal epilepsy patients who had a well-defined focus determined through intracranial EEG (iEEG), epileptogenic MRI lesions or surgical resection. dMSI accuracy and spatial properties were quantitatively estimated as: (a) distance from the epilepsy focus, (b) reproducibility, (c) spatial dispersion (SD), (d) map extension, and (e) effect of thresholding on map properties. Clinical performance was excellent for all methods (median distance from the focus MNE = 2.4 mm; sLORETA = 3.5 mm; cMEM = 3.5 mm; dSPM = 6.8 mm, Ave = 0 mm). Ave showed the lowest distance between the map maximum and epilepsy focus (Dmin lower than cMEM, MNE, and dSPM, p = .021, p = .008, p < .001, respectively). cMEM showed the best spatial features, with lowest SD outside the focus (SD lower than all other methods, p < .001 consistently) and high contrast between the generator and surrounding regions. The average map Ave provided the best localization accuracy, whereas cMEM exhibited the lowest amount of spurious distant activity. dMSI techniques have the potential to significantly improve identification of iEEG targets and to guide surgical planning, especially when multiple methods are combined.

Magnetoencephalography resting state connectivity patterns as indicatives of surgical outcome in epilepsy patients.

OBJECTIVE: Focal epilepsy is a disorder affecting several brain networks; however, epilepsy surgery usually targets a restricted region, the so-called epileptic focus. There is a growing interest in embedding resting state (RS) connectivity analysis into pre-surgical workup. APPROACH: In this retrospective study, we analyzed Magnetoencephalography (MEG) long-range RS functional connectivity patterns in patients with drug-resistant focal epilepsy. MEG recorded prior to surgery from seven seizure-free (Engel Ia) and five non seizure-free (Engel III or IV) patients were analyzed (minimum 2-years post-surgical follow-up). MEG segments without any detectable epileptic activity were source localized using wavelet-based Maximum Entropy on the Mean method. Amplitude envelope correlation in the theta (4-8 Hz), alpha (8-13 Hz), and beta (13-26 Hz) bands were used for assessing connectivity. MAIN RESULTS: For seizure-free patients, we found an isolated epileptic network characterized by weaker connections between the brain region where interictal epileptic discharges (IED) are generated and the rest of the cortex, when compared to connectivity between the corresponding contralateral homologous region and the rest of the cortex. Contrarily, non seizure-free patients exhibited a widespread RS epileptic network characterized by stronger connectivity between the IED generator and the rest of the cortex, in comparison to the contralateral region and the cortex. Differences between the two seizure outcome groups concerned mainly distant long-range connections and were found in the alpha-band. SIGNIFICANCE: Importantly, these connectivity patterns suggest specific mechanisms describing the underlying organization of the epileptic network and were detectable at the individual patient level, supporting the prospect use of MEG connectivity patterns in epilepsy to predict post-surgical seizure outcome.