Cerebral Gray Matter May Not Explain Sleep Slow-Wave Characteristics after Severe Brain Injury.

Sleep slow waves are the hallmark of deeper non-rapid eye movement sleep. It is generally assumed that gray matter properties predict slow-wave density, morphology, and spectral power in healthy adults. Here, we tested the association between gray matter volume (GMV) and slow-wave characteristics in 27 patients with moderate-to-severe traumatic brain injury (TBI, 32.0 ± 12.2 years old, eight women) and compared that with 32 healthy controls (29.2 ± 11.5 years old, nine women). Participants underwent overnight polysomnography and cerebral MRI with a 3 Tesla scanner. A whole-brain voxel-wise analysis was performed to compare GMV between groups. Slow-wave density, morphology, and spectral power (0.4-6 Hz) were computed, and GMV was extracted from the thalamus, cingulate, insula, precuneus, and orbitofrontal cortex to test the relationship between slow waves and gray matter in regions implicated in the generation and/or propagation of slow waves. Compared with controls, TBI patients had significantly lower frontal and temporal GMV and exhibited a subtle decrease in slow-wave frequency. Moreover, higher GMV in the orbitofrontal cortex, insula, cingulate cortex, and precuneus was associated with higher slow-wave frequency and slope, but only in healthy controls. Higher orbitofrontal GMV was also associated with higher slow-wave density in healthy participants. While we observed the expected associations between GMV and slow-wave characteristics in healthy controls, no such associations were observed in the TBI group despite lower GMV. This finding challenges the presumed role of GMV in slow-wave generation and morphology.

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.

Oscillatory and nonoscillatory sleep electroencephalographic biomarkers of the epileptic network.

OBJECTIVE: In addition to the oscillatory brain activity, the nonoscillatory (scale-free) components of the background electroencephalogram (EEG) may provide further information about the complexity of the underlying neuronal network. As epilepsy is considered a network disease, such scale-free metrics might help to delineate the epileptic network. Here, we performed an analysis of the sleep oscillatory (spindle, slow wave, and rhythmic spectral power) and nonoscillatory (H exponent) intracranial EEG using multiple interictal features to estimate whether and how they deviate from normalcy in 38 adults with drug-resistant epilepsy. METHODS: To quantify intracranial EEG abnormalities within and outside the seizure onset areas, patients' values were adjusted based on normative maps derived from the open-access Montreal Neurological Institute open iEEG Atlas. In a subset of 29 patients who underwent resective surgery, we estimated the predictive value of these features to identify the epileptogenic zone in those with a good postsurgical outcome. RESULTS: We found that distinct sleep oscillatory and nonoscillatory metrics behave differently across the epileptic network, with the strongest differences observed for (1) a reduction in spindle activity (spindle rates and rhythmic sigma power in the 10-16 Hz band), (2) a higher rhythmic gamma power (30-80 Hz), and (3) a higher H exponent (steeper 1/f slope). As expected, epileptic spikes were also highest in the seizure onset areas. Furthermore, in surgical patients, the H exponent achieved the highest performance (balanced accuracy of .76) for classifying resected versus nonresected channels in good outcome patients. SIGNIFICANCE: This work suggests that nonoscillatory components of the intracranial EEG signal could serve as promising interictal sleep candidates of epileptogenicity in patients with drug-resistant epilepsy. Our findings further advance the understanding of epilepsy as a disease, whereby absence or loss of sleep physiology may provide information complementary to pathological epileptic processes.

Technical challenges in REM sleep microstructure classification: A study of patients with REM sleep behaviour disorder.

While commonly treated as a uniform state in practice, rapid eye movement sleep contains two distinct microstructures-phasic (presence of rapid eye movement) and tonic (no rapid eye movement). This study aims to identify technical challenges during rapid eye movement sleep microstructure visual classification in patients with rapid eye movement sleep behaviour disorder, and to propose solutions to enhance reliability between scorers. Fifty-seven sleep recordings were randomly allocated into three subsequent batches (n = 10, 13 and 34) for scoring. To reduce single-centre bias, we recruited three raters/scorers, with each trained from a different institution. Two raters independently scored each 30-s rapid eye movement sleep into 10â€…× fSEM3-s phasic/tonic microstructures based on the AASM guidelines. The third rater acted as an "arbitrator" to resolve opposite opinions persisting during the revision between batches. Besides interrater differences in artefact rejection rate, interrater variance frequently occurred due to transitioning between microstructures and moderate-to-severe muscular/electrode artefact interference. To enhance interrater agreement, a rapid eye movement scoring schematic graph was developed, incorporating proxy electrode use, filters and cut-offs for microstructure transitioning. To assess potential effectiveness of the schematic graph proposed, raters were instructed to systematically apply it in scoring for the third batch. Of the 34 recordings, 27 reached a Cohen's kappa score above 0.8 (i.e. almost perfect agreement between raters), significantly improved from the prior batches (p = 0.0003, Kruskal-Wallis test). Our study illustrated potential solutions and guidance for challenges that may be encountered during rapid eye movement sleep microstructure classification.

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.

Optoelectronic Pressure Sensor Based on the Bending Loss of Plastic Optical Fibers Embedded in Stretchable Polydimethylsiloxane.

We report the design and testing of a sensor pad based on optical and flexible materials for the development of pressure monitoring devices. This project aims to create a flexible and low-cost pressure sensor based on a two-dimensional grid of plastic optical fibers embedded in a pad of flexible and stretchable polydimethylsiloxane (PDMS). The opposite ends of each fiber are connected to an LED and a photodiode, respectively, to excite and measure light intensity changes due to the local bending of the pressure points on the PDMS pad. Tests were performed in order to study the sensitivity and repeatability of the designed flexible pressure sensor.

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.

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.

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.

White matter information flow mapping from diffusion MRI and EEG.

The human brain can be described as a network of specialized and spatially distributed regions. The activity of individual regions can be estimated using electroencephalography and the structure of the network can be measured using diffusion magnetic resonance imaging. However, the communication between the different cortical regions occurring through the white matter, coined information flow, cannot be observed by either modalities independently. Here, we present a new method to infer information flow in the white matter of the brain from joint diffusion MRI and EEG measurements. This is made possible by the millisecond resolution of EEG which makes the transfer of information from one region to another observable. A subject specific Bayesian network is built which captures the possible interactions between brain regions at different times. This network encodes the connections between brain regions detected using diffusion MRI tractography derived white matter bundles and their associated delays. By injecting the EEG measurements as evidence into this model, we are able to estimate the directed dynamical functional connectivity whose delays are supported by the diffusion MRI derived structural connectivity. We present our results in the form of information flow diagrams that trace transient communication between cortical regions over a functional data window. The performance of our algorithm under different noise levels is assessed using receiver operating characteristic curves on simulated data. In addition, using the well-characterized visual motor network as grounds to test our model, we present the information flow obtained during a reaching task following left or right visual stimuli. These promising results present the transfer of information from the eyes to the primary motor cortex. The information flow obtained using our technique can also be projected back to the anatomy and animated to produce videos of the information path through the white matter, opening a new window into multi-modal dynamic brain connectivity.

The effects of bandpass filtering on the oscillatory potentials of the electroretinogram.

PURPOSE: In order to study the OPs, the ERG signal must be filtered to eliminate the low-frequency waves known as the a- and b-waves. Unfortunately, the ISCEV ERG standard does not give clear guidelines on how to proceed apart from indicating that frequencies below 75 Hz should be filtered out when recording scotopic OPs, while no suggestions are offered for the photopic OPs. The purpose of this study was thus to characterize more extensively the effects of various digital filters on the photopic OP waveforms in order to suggest the most appropriate filtering method to record them. METHODS: Filtered OPs (N = 9600 tracings) were extracted from a photopic ERG databank of 40 normal subjects [intensity: 4.4 cd s m-2; background: 30 cd m-2] using 240 different combinations of five digital filters types (Bessel; Butterworth; Elliptic; Chebyshev type 1 and 2), eight bandwidth ranges (50-300; 75-300; 100-300; 125-300; 50-1000; 75-1000; 100-1000; 125-1000 Hz), three filter orders (1, 2 and 5) and with/without phase lag corrections that were generated using MATLAB 2015b. The peak time and the percentage of OPs (sum of OP amplitudes on the b-wave amplitude) were calculated in the time domain (TD%OP). RESULTS: The timing of the OPs was less affected than the amplitude by the different filters used. Depending on the filter used, the resulting OPs were either severely depressed (16.16% of broadband OP content) or slightly reduced (93.63%). The filters that most successfully eliminated the slow components of the ERG (i.e., < 12% of broadband value) were the Bessel, the Butterworth and the Chebyshev type 1 filters and out of the latter, the Butterworth filter was that which most faithfully reproduced the high-frequency OPs (i.e., > 96%). CONCLUSION: Our results vividly demonstrate the need to better define the characteristics of the filter that is used to record the OPs as it does have a significant impact on the resulting waveform.

Clinical yield of magnetoencephalography distributed source imaging in epilepsy: A comparison with equivalent current dipole method.

OBJECTIVE: Source localization of interictal epileptic discharges (IEDs) is clinically useful in the presurgical workup of epilepsy patients. It is usually obtained by equivalent current dipole (ECD) which localizes a point source and is the only inverse solution approved by clinical guidelines. In contrast, magnetic source imaging using distributed methods (dMSI) provides maps of the location and the extent of the generators, but its yield has not been clinically validated. We systematically compared ECD versus dMSI performed using coherent Maximum Entropy on the Mean (cMEM), a method sensitive to the spatial extent of the generators. METHODS: 340 source localizations of IEDs derived from 49 focal epilepsy patients with foci well-defined through intracranial EEG, MRI lesions, and surgery were analyzed. The comparison was based on the assessment of the sublobar concordance with the focus and of the distance between the source and the focus. RESULTS: dMSI sublobar concordance was significantly higher than ECD (81% vs 69%, P < 0.001), especially for extratemporal lobe sources (dMSI = 84%; ECD = 67%, P < 0.001) and for seizure free patients (dMSI = 83%; ECD = 70%, P < 0.001). The median distance from the focus was 4.88 mm for ECD and 3.44 mm for dMSI (P < 0.001). ECD dipoles were often wrongly localized in deep brain regions. CONCLUSIONS: dMSI using cMEM exhibited better accuracy. dMSI also offered the advantage of recovering more realistic maps of the generator, which could be exploited for neuronavigation aimed at targeting invasive EEG and surgical resection. Therefore, dMSI may be preferred to ECD in clinical practice. Hum Brain Mapp 39:218-231, 2018. © 2017 Wiley Periodicals, Inc.

Reproducibility of EEG-MEG fusion source analysis of interictal spikes: Relevance in presurgical evaluation of epilepsy.

Fusion of electroencephalography (EEG) and magnetoencephalography (MEG) data using maximum entropy on the mean method (MEM-fusion) takes advantage of the complementarities between EEG and MEG to improve localization accuracy. Simulation studies demonstrated MEM-fusion to be robust especially in noisy conditions such as single spike source localizations (SSSL). Our objective was to assess the reliability of SSSL using MEM-fusion on clinical data. We proposed to cluster SSSL results to find the most reliable and consistent source map from the reconstructed sources, the so-called consensus map. Thirty-four types of interictal epileptic discharges (IEDs) were analyzed from 26 patients with well-defined epileptogenic focus. SSSLs were performed on EEG, MEG, and fusion data and consensus maps were estimated using hierarchical clustering. Qualitative (spike-to-spike reproducibility rate, SSR) and quantitative (localization error and spatial dispersion) assessments were performed using the epileptogenic focus as clinical reference. Fusion SSSL provided significantly better results than EEG or MEG alone. Fusion found at least one cluster concordant with the clinical reference in all cases. This concordant cluster was always the one involving the highest number of spikes. Fusion yielded highest reproducibility (SSR EEG = 55%, MEG = 71%, fusion = 90%) and lowest localization error. Also, using only few channels from either modality (21EEG + 272MEG or 54EEG + 25MEG) was sufficient to reach accurate fusion. MEM-fusion with consensus map approach provides an objective way of finding the most reliable and concordant generators of IEDs. We, therefore, suggest the pertinence of SSSL using MEM-fusion as a valuable clinical tool for presurgical evaluation of epilepsy.

Cortical thinning explains changes in sleep slow waves during adulthood.

Sleep slow waves (SWs) change considerably throughout normal aging. In humans, SWs are generated and propagate on a structural backbone of highly interconnected cortical regions that form most of the default mode network, such as the insula, cingulate cortices, temporal lobe, parietal lobe, and medial frontal lobe. Regions in this network undergo cortical thinning and breakdown in structural and functional connectivity over the course of normal aging. In this study, we investigated how changes in cortical thickness (CT), a measure of gray matter integrity, are involved in modifications of sleep SWs during adulthood in humans. Thirty young (mean age = 23.49 years; SD = 2.79) and 33 older (mean age = 60.35 years; SD = 5.71) healthy subjects underwent a nocturnal polysomnography and T1 MRI. We show that, when controlling for age, higher SW density (nb/min of nonrapid eye movement sleep) was associated with higher CT in cortical regions involved in SW generation surrounding the lateral fissure (insula, superior temporal, parietal, middle frontal), whereas higher SW amplitude was associated with higher CT in middle frontal, medial prefrontal, and medial posterior regions. Mediation analyses demonstrated that thinning in a network of cortical regions involved in SW generation and propagation, but also in cognitive functions, explained the age-related decrease in SW density and amplitude. Altogether, our results suggest that microstructural degradation of specific cortical regions compromise SW generation and propagation in older subjects, critically contributing to age-related changes in SW oscillations.

SPARK: Sparsity-based analysis of reliable k-hubness and overlapping network structure in brain functional connectivity.

Functional hubs are defined as the specific brain regions with dense connections to other regions in a functional brain network. Among them, connector hubs are of great interests, as they are assumed to promote global and hierarchical communications between functionally specialized networks. Damage to connector hubs may have a more crucial effect on the system than does damage to other hubs. Hubs in graph theory are often identified from a correlation matrix, and classified as connector hubs when the hubs are more connected to regions in other networks than within the networks to which they belong. However, the identification of hubs from functional data is more complex than that from structural data, notably because of the inherent problem of multicollinearity between temporal dynamics within a functional network. In this context, we developed and validated a method to reliably identify connectors and corresponding overlapping network structure from resting-state fMRI. This new method is actually handling the multicollinearity issue, since it does not rely on counting the number of connections from a thresholded correlation matrix. The novelty of the proposed method is that besides counting the number of networks involved in each voxel, it allows us to identify which networks are actually involved in each voxel, using a data-driven sparse general linear model in order to identify brain regions involved in more than one network. Moreover, we added a bootstrap resampling strategy to assess statistically the reproducibility of our results at the single subject level. The unified framework is called SPARK, i.e. SParsity-based Analysis of Reliable k-hubness, where k-hubness denotes the number of networks overlapping in each voxel. The accuracy and robustness of SPARK were evaluated using two dimensional box simulations and realistic simulations that examined detection of artificial hubs generated on real data. Then, test/retest reliability of the method was assessed using the 1000 Functional Connectome Project database, which includes data obtained from 25 healthy subjects at three different occasions with long and short intervals between sessions. We demonstrated that SPARK provides an accurate and reliable estimation of k-hubness, suggesting a promising tool for understanding hub organization in resting-state fMRI.

Intracranial EEG potentials estimated from MEG sources: A new approach to correlate MEG and iEEG data in epilepsy.

Detection of epileptic spikes in MagnetoEncephaloGraphy (MEG) requires synchronized neuronal activity over a minimum of 4cm2. We previously validated the Maximum Entropy on the Mean (MEM) as a source localization able to recover the spatial extent of the epileptic spike generators. The purpose of this study was to evaluate quantitatively, using intracranial EEG (iEEG), the spatial extent recovered from MEG sources by estimating iEEG potentials generated by these MEG sources. We evaluated five patients with focal epilepsy who had a pre-operative MEG acquisition and iEEG with MRI-compatible electrodes. Individual MEG epileptic spikes were localized along the cortical surface segmented from a pre-operative MRI, which was co-registered with the MRI obtained with iEEG electrodes in place for identification of iEEG contacts. An iEEG forward model estimated the influence of every dipolar source of the cortical surface on each iEEG contact. This iEEG forward model was applied to MEG sources to estimate iEEG potentials that would have been generated by these sources. MEG-estimated iEEG potentials were compared with measured iEEG potentials using four source localization methods: two variants of MEM and two standard methods equivalent to minimum norm and LORETA estimates. Our results demonstrated an excellent MEG/iEEG correspondence in the presumed focus for four out of five patients. In one patient, the deep generator identified in iEEG could not be localized in MEG. MEG-estimated iEEG potentials is a promising method to evaluate which MEG sources could be retrieved and validated with iEEG data, providing accurate results especially when applied to MEM localizations. Hum Brain Mapp 37:1661-1683, 2016. © 2016 Wiley Periodicals, Inc.

Source localization of the seizure onset zone from ictal EEG/MEG data.

INTRODUCTION: Surgical treatment of drug-resistant epilepsy relies on the identification of the seizure onset zone (SOZ) and often requires intracranial EEG (iEEG). We have developed a new approach for non-invasive magnetic and electric source imaging of the SOZ (MSI-SOZ and ESI-SOZ) from ictal magnetoencephalography (MEG) and EEG recordings, using wavelet-based Maximum Entropy on the Mean (wMEM) method. We compared the performance of MSI-SOZ and ESI-SOZ with interictal spike source localization (MSI-spikes and ESI-spikes) and clinical localization of the SOZ (i.e., based on iEEG or lesion topography, denoted as clinical-SOZ). METHODS: A total of 46 MEG or EEG seizures from 13 patients were analyzed. wMEM was applied around seizure onset, centered on the frequency band showing the strongest power change. Principal component analysis applied to spatiotemporal reconstructed wMEM sources (0.4-1 s around seizure onset) identified the main spatial pattern of ictal oscillations. Qualitative sublobar concordance and quantitative measures of distance and spatial overlaps were estimated to compare MSI/ESI-SOZ with MSI/ESI-Spikes and clinical-SOZ. RESULTS: MSI/ESI-SOZ were concordant with clinical-SOZ in 81% of seizures (MSI 90%, ESI 64%). MSI-SOZ was more accurate and identified sources closer to the clinical-SOZ (P = 0.012) and to MSI-Spikes (P = 0.040) as compared with ESI-SOZ. MSI/ESI-SOZ and MSI/ESI-Spikes did not differ in terms of concordance and distance from the clinical-SOZ. CONCLUSIONS: wMEM allows non-invasive localization of the SOZ from ictal MEG and EEG. MSI-SOZ performs better than ESI-SOZ. MSI/ESI-SOZ can provide important additional information to MSI/ESI-Spikes during presurgical evaluation. Hum Brain Mapp 37:2528-2546, 2016. © 2016 Wiley Periodicals, Inc.

Detection and Magnetic Source Imaging of Fast Oscillations (40-160 Hz) Recorded with Magnetoencephalography in Focal Epilepsy Patients.

We present a framework to detect fast oscillations (FOs) in magnetoencephalography (MEG) and to perform magnetic source imaging (MSI) to determine the location and extent of their generators in the cortex. FOs can be of physiologic origin associated to sensory processing and memory consolidation. In epilepsy, FOs are of pathologic origin and biomarkers of the epileptogenic zone. Seventeen patients with focal epilepsy previously confirmed with identified FOs in scalp electroencephalography (EEG) were evaluated. To handle data deriving from large number of sensors (275 axial gradiometers) we used an automatic detector with high sensitivity. False positives were discarded by two human experts. MSI of the FOs was performed with the wavelet based maximum entropy on the mean method. We found FOs in 11/17 patients, in only one patient the channel with highest FO rate was not concordant with the epileptogenic region and might correspond to physiologic oscillations. MEG FOs rates were very low: 0.02-4.55 per minute. Compared to scalp EEG, detection sensitivity was lower, but the specificity higher in MEG. MSI of FOs showed concordance or partial concordance with proven generators of seizures and epileptiform activity in 10/11 patients. We have validated the proposed framework for the non-invasive study of FOs with MEG. The excellent overall concordance with other clinical gold standard evaluation tools indicates that MEG FOs can provide relevant information to guide implantation for intracranial EEG pre-surgical evaluation and for surgical treatment, and demonstrates the important added value of choosing appropriate FOs detection and source localization methods.