Deep brain activities can be detected with magnetoencephalography

Spatiotemporal whole-brain activity and functional connectivity of melodies recognition

Music is a non-verbal human language, built on logical, hierarchical structures, that offers excellent opportunities to explore how the brain processes complex spatiotemporal auditory sequences. Using the high temporal resolution of magnetoencephalography, we investigated the unfolding brain dynamics of 70 participants during the recognition of previously memorized musical sequences compared to novel sequences matched in terms of entropy and information content. Measures of both whole-brain activity and functional connectivity revealed a widespread brain network underlying the recognition of the memorized auditory sequences, which comprised primary auditory cortex, superior temporal gyrus, insula, frontal operculum, cingulate gyrus, orbitofrontal cortex, basal ganglia, thalamus, and hippocampus. Furthermore, while the auditory cortex responded mainly to the first tones of the sequences, the activity of higher-order brain areas such as the cingulate gyrus, frontal operculum, hippocampus, and orbitofrontal cortex largely increased over time during the recognition of the memorized versus novel musical sequences. In conclusion, using a wide range of analytical techniques spanning from decoding to functional connectivity and building on previous works, our study provided new insights into the spatiotemporal whole-brain mechanisms for conscious recognition of auditory sequences.

Simultaneous invasive and non-invasive recordings in humans: A novel Rosetta stone for deciphering brain activity

Simultaneous noninvasive and invasive electrophysiological recordings provide a unique opportunity to achieve a comprehensive understanding of human brain activity, much like a Rosetta stone for human neuroscience. In this review we focus on the increasingly-used powerful combination of intracranial electroencephalography (iEEG) with scalp electroencephalography (EEG) or magnetoencephalography (MEG). We first provide practical insight on how to achieve these technically challenging recordings. We then provide examples from clinical research on how simultaneous recordings are advancing our understanding of epilepsy. This is followed by the illustration of how human neuroscience and methodological advances could benefit from these simultaneous recordings. We conclude with a call for open data sharing and collaboration, while ensuring neuroethical approaches and argue that only with a true collaborative approach the promises of simultaneous recordings will be fulfilled.

Effect of Inverse Solutions, Connectivity Measures, and Node Sizes on EEG Source Network: A Simultaneous EEG Study

Brain network provides an essential perspective for studying normal and pathological brain activities. Reconstructing the brain network in the source space becomes more needed, for example, as a target in non-invasive neuromodulation. Precise estimating source activities from the scalp EEG is still challenging because it is an ill-posed question and because of the volume conduction effect. There is no consensus on how to reconstruct the EEG source network. This study uses simultaneous scalp EEG and stereo-EEG to investigate the effect of inverse solutions, connectivity measures, and node sizes on the reconstruction of the source network. We evaluated the performance of different methods on both source activity and network. Numerical simulation was also carried out for comparison. The weighted phase-lag index (wPLI) method achieved significantly better performance on the reconstructed networks in source space than five other connectivity measures (directed transfer function (DTF), partial directed coherence (PDC), efficient effective connectivity (EEC), Pearson correlation coefficient (PCC), and amplitude envelope correlation (AEC)). There is no significant difference between the inverse solutions (standardized low-resolution brain electromagnetic tomography (sLORETA), weighted minimum norm estimate (wMNE), and linearly constrained minimum variance (LCMV) beamforming) on the reconstructed source networks. The source network based on signal phases can fit intracranial activities better than signal waveform properties or causality. Our study provides a basis for reconstructing source space networks from scalp EEG, especially for future neuromodulation research.

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.

Hierarchical syntax model of music predicts theta power during music listening

Linguistic research showed that the depth of syntactic embedding is reflected in brain theta power. Here, we test whether this also extends to non-linguistic stimuli, specifically music. We used a hierarchical model of musical syntax to continuously quantify two types of expert-annotated harmonic dependencies throughout a piece of Western classical music: prolongation and preparation. Prolongations can roughly be understood as a musical analogue to linguistic coordination between constituents that share the same function (e.g., 'pizza' and 'pasta' in 'I ate pizza and pasta'). Preparation refers to the dependency between two harmonies whereby the first implies a resolution towards the second (e.g., dominant towards tonic; similar to how the adjective implies the presence of a noun in 'I like spicy … '). Source reconstructed MEG data of sixty-five participants listening to the musical piece was then analysed. We used Bayesian Mixed Effects models to predict theta envelope in the brain, using the number of open prolongation and preparation dependencies as predictors whilst controlling for audio envelope. We observed that prolongation and preparation both carry independent and distinguishable predictive value for theta band fluctuation in key linguistic areas such as the Angular, Superior Temporal, and Heschl's Gyri, or their right-lateralised homologues, with preparation showing additional predictive value for areas associated with the reward system and prediction. Musical expertise further mediated these effects in language-related brain areas. Results show that predictions of precisely formalised music-theoretical models are reflected in the brain activity of listeners which furthers our understanding of the perception and cognition of musical structure.

Beyond Stereo-EEG: Is It Worth Combining Stereo-EEG With Other Diagnostic Methods?

SUMMARY

Stereo-EEG is a widely used method to improve the diagnostic precision of presurgical workup in patients with refractory epilepsy. Its ability to detect epileptic activity and identify epileptic networks largely depends on the chosen implantation strategy. Even in an ideal situation, electrodes record activity generated in <10% of the brain and contacts only record from brain tissue in their immediate proximity. In this article, the authors discuss how recording stereo-EEG simultaneously with other diagnostic methods can improve its diagnostic value in clinical and research settings. It can help overcome the limited spatial coverage of intracranial recording and better understand the sources of epileptic activity. Simultaneous scalp EEG is the most widely available method, often used to understand large epileptic networks, seizure propagation, and EEG activity occurring on the contralateral hemisphere. Simultaneous magnetoencephalography allows for more precise source localization and identification of deep sources outside the stereo-EEG coverage. Finally, simultaneous functional MRI can highlight metabolic changes following epileptic activity and help understand the widespread network changes associated with interictal activity. This overview highlights advantages and methodological challenges for all these methods. Clinical use and research applications are presented for each approach.

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

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

Theta and alpha oscillations in human hippocampus and medial parietal cortex support the formation of location-based representations

Our ability to navigate in a new environment depends on learning new locations. Mental representations of locations are quickly accessible during navigation and allow us to know where we are regardless of our current viewpoint. Recent functional magnetic resonance imaging (fMRI) research using pattern classification has shown that these location-based representations emerge in the retrosplenial cortex and parahippocampal gyrus, regions theorized to be critically involved in spatial navigation. However, little is currently known about the oscillatory dynamics that support the formation of location-based representations. We used magnetoencephalogram (MEG) recordings to investigate region-specific oscillatory activity in a task where participants could form location-based representations. Participants viewed videos showing that two perceptually distinct scenes (180° apart) belonged to the same location. This "overlap" video allowed participants to bind the two distinct scenes together into a more coherent location-based representation. Participants also viewed control "non-overlap" videos where two distinct scenes from two different locations were shown, where no location-based representation could be formed. In a post-video behavioral task, participants successfully matched the two viewpoints shown in the overlap videos, but not the non-overlap videos, indicating they successfully learned the locations in the overlap condition. Comparing oscillatory activity between the overlap and non-overlap videos, we found greater theta and alpha/beta power during the overlap relative to non-overlap videos, specifically at time-points when we expected scene integration to occur. These oscillations localized to regions in the medial parietal cortex (precuneus and retrosplenial cortex) and the medial temporal lobe, including the hippocampus. Therefore, we find that theta and alpha/beta oscillations in the hippocampus and medial parietal cortex are likely involved in the formation of location-based representations.

Spontaneous HFO Sequences Reveal Propagation Pathways for Precise Delineation of Epileptogenic Networks

Epilepsy, a neurological disorder affecting millions worldwide, poses great challenges in precisely delineating the epileptogenic zone - the brain region generating seizures - for effective treatment. High-frequency oscillations (HFOs) are emerging as promising biomarkers; however, the clinical utility is hindered by the difficulties in distinguishing pathological HFOs from non- epileptiform activities at single electrode and single patient resolution and understanding their dynamic role in epileptic networks. Here, we introduce an HFO-sequencing approach to analyze spontaneous HFOs traversing cortical regions in 40 drug-resistant epilepsy patients. This data- driven method automatically detected over 8.9 million HFOs, pinpointing pathological HFO- networks, and unveiled intricate millisecond-scale spatiotemporal dynamics, stability, and functional connectivity of HFOs in prolonged intracranial EEG recordings. These HFO sequences demonstrated a significant improvement in localization of epileptic tissue, with an 818.47% increase in concordance with seizure-onset zone (mean error: 2.92 mm), compared to conventional benchmarks. They also accurately predicted seizure outcomes for 90% AUC based on pre-surgical information using generalized linear models. Importantly, this mapping remained reliable even with short recordings (mean standard deviation: 3.23 mm for 30-minute segments). Furthermore, HFO sequences exhibited distinct yet highly repetitive spatiotemporal patterns, characterized by pronounced synchrony and predominant inward information flow from periphery towards areas involved in propagation, suggesting a crucial role for excitation-inhibition balance in HFO initiation and progression. Together, these findings shed light on the intricate organization of epileptic network and highlight the potential of HFO-sequencing as a translational tool for improved diagnosis, surgical targeting, and ultimately, better outcomes for vulnerable patients with drug-resistant epilepsy.

One Sentence Summary

Pathological fast brain oscillations travel like traffic along varied routes, outlining recurrently visited neural sites emerging as critical hotspots in epilepsy network.

Non-electrocardiogram-gated, free-breathing, off-resonance reduced, high-resolution, whole-heart myocardial T2 * mapping at 3 T within 5 min

PURPOSE

Widely used conventional 2D T2 * approaches that are based on breath-held, electrocardiogram (ECG)-gated, multi-gradient-echo sequences are prone to motion artifacts in the presence of incomplete breath holding or arrhythmias, which is common in cardiac patients. To address these limitations, a 3D, non-ECG-gated, free-breathing T2 * technique that enables rapid whole-heart coverage was developed and validated.

METHODS

A continuous random Gaussian 3D k-space sampling was implemented using a low-rank tensor framework for motion-resolved 3D T2 * imaging. This approach was tested in healthy human volunteers and in swine before and after intravenous administration of ferumoxytol.

RESULTS

Spatial-resolution matched T2 * images were acquired with 2-3-fold reduction in scan time using the proposed T2 * mapping approach relative to conventional T2 * mapping. Compared with the conventional approach, T2 * images acquired with the proposed method demonstrated reduced off-resonance and flow artifacts, leading to higher image quality and lower coefficient of variation in T2 *-weighted images of the myocardium of swine and humans. Mean myocardial T2 * values determined using the proposed and conventional approaches were highly correlated and showed minimal bias.

CONCLUSION

The proposed non-ECG-gated, free-breathing, 3D T2 * imaging approach can be performed within 5 min or less. It can overcome critical image artifacts from undesirable cardiac and respiratory motion and bulk off-resonance shifts at the heart-lung interface. The proposed approach is expected to facilitate faster and improved cardiac T2 * mapping in those with limited breath-holding capacity or arrhythmias.

Magnetoencephalographic detection of synchronized epileptic activity between the hippocampus and insular cortex

Most magnetoencephalographic signals are derived from synchronized activity in the brain surface cortex. By contrast, the contribution of synchronized activity in the deep brain to magnetoencephalography (MEG) has remained unclear. We compared stereotactic electroencephalography (sEEG) with simultaneous MEG findings in a patient with temporal lobe epilepsy to determine the conditions under which MEG could also detect sEEG findings. The synchrony and similarity of the waves were evaluated using visual inspection and wavelet coherence. A 45-year-old woman with intractable temporal lobe epilepsy underwent sEEG and MEG simultaneously to determine the laterality and precise location of the epileptic focus. When spike-and-waves were seen in the right hippocampal head alone, no distinct spike-and-waves were observed visually in the right temporal MEG. The seizure then spread to the right insula on sEEG with a rhythmic theta frequency while synchronous activity was observed in the right temporal MEG channels. When polyspikes appeared in the right hippocampus, the right temporal MEG showed electrical activity with relatively high similarity to that of the right hippocampal head and insular cortex but less similarity to that of the right lateral temporal lobe cortex. MEG might detect epileptic activity synchronized between the hippocampus and insular cortex.

Identification of Early Hippocampal Dynamics during Recognition Memory with Independent Component Analysis

The hippocampus is generally considered to have relatively late involvement in recognition memory, its main electrophysiological signature being between 400 and 800 ms after stimulus onset. However, most electrophysiological studies have analyzed the hippocampus as a single responsive area, selecting only a single-site signal exhibiting the strongest effect in terms of amplitude. These classical approaches may not capture all the dynamics of this structure, hindering the contribution of other hippocampal sources that are not located in the vicinity of the selected site. We combined intracerebral electroencephalogram recordings from epileptic patients with independent component analysis during a recognition memory task involving the recognition of old and new images. We identified two sources with different responses emerging from the hippocampus: a fast one (maximal amplitude at ∼250 ms) that could not be directly identified from raw recordings and a latter one, peaking at ∼400 ms. The former component presented different amplitudes between old and new items in 6 out of 10 patients. The latter component had different delays for each condition, with a faster activation (∼290 ms after stimulus onset) for recognized items. We hypothesize that both sources represent two steps of hippocampal recognition memory, the faster reflecting the input from other structures and the latter the hippocampal internal processing. Recognized images evoking early activations would facilitate neural computation in the hippocampus, accelerating memory retrieval of complementary information. Overall, our results suggest that the hippocampal activity is composed of several sources with an early activation related to recognition memory.

Towards the automated detection of interictal epileptiform discharges with magnetoencephalography

BACKGROUND

The analysis of clinical magnetoencephalography (MEG) in patients with epilepsy traditionally relies on visual identification of interictal epileptiform discharges (IEDs), which is time consuming and dependent on subjective criteria.

NEW METHOD

Here, we explore the ability of Independent Components Analysis (ICA) and Hidden Markov Modeling (HMM) to automatically detect and localize IEDs. We tested our pipelines on resting-state MEG recordings from 10 school-aged children with (multi)focal epilepsy.

RESULTS

In focal epilepsy patients, both pipelines successfully detected visually identified IEDs, but also revealed unidentified low-amplitude IEDs. Success was more mitigated in patients with multifocal epilepsy, as our automated pipeline missed IED activity associated with some foci-an issue that could be alleviated by post-hoc manual selection of epileptiform ICs or HMM states.

COMPARISON WITH EXISTING METHODS

We compared our results with visual IED detection by an experienced clinical magnetoencephalographer, getting heightened sensitivity and requiring minimal input from clinical practitioners.

CONCLUSIONS

IED detection based on ICA or HMM represents an efficient way to identify IED localization and timing. The development of these automatic IED detection algorithms provide a step forward in clinical MEG practice by decreasing the duration of MEG analysis and enhancing its sensitivity.

From cryogenic to on-scalp magnetoencephalography for the evaluation of paediatric epilepsy

Magnetoencephalography (MEG) is a neurophysiological technique based on the detection of brain magnetic fields. Whole-head MEG systems typically house a few hundred sensors requiring cryogenic cooling in a rigid one-size-fits-all (commonly adult-sized) helmet to keep a thermal insulation space. This leads to an increased brain-to-sensor distance in children, because of their smaller head circumference, and decreased signal-to-noise ratio. MEG allows detection and localization of interictal and ictal epileptiform discharges, and pathological high frequency oscillations, as a part of the presurgical assessment of children with refractory focal epilepsy, where electroencephalography is not contributive. MEG can also map the eloquent cortex before surgical resection. MEG also provides insights into the physiopathology of both generalized and focal epilepsy. On-scalp recordings based on cryogenic-free sensors have demonstrated their use in the field of childhood focal epilepsy and should become a reference technique for diagnosing epilepsy in the paediatric population. WHAT THIS PAPER ADDS: Magnetoencephalography (MEG) contributes to the diagnosis and understanding of paediatric epilepsy. On-scalp MEG recordings demonstrate some advantages over cryogenic MEG.

Combining independent component analysis and source localization for improving spatial sampling of stereoelectroencephalography in epilepsy

Stereoelectroencephalography is a powerful intracerebral EEG recording method for the presurgical evaluation of epilepsy. It consists in implanting depth electrodes in the patient's brain to record electrical activity and map the epileptogenic zone, which should be resected to render the patient seizure-free. Stereoelectroencephalography has high spatial accuracy and signal-to-noise ratio but remains limited in the coverage of the explored brain regions. Thus, the implantation might provide a suboptimal sampling of epileptogenic regions. We investigate the potential of improving a suboptimal stereoelectroencephalography recording by performing source localization on stereoelectroencephalography signals. We propose combining independent component analysis, connectivity measures to identify components of interest, and distributed source modelling. This approach was tested on two patients with two implantations each, the first failing to characterize the epileptogenic zone and the second giving a better diagnosis. We demonstrate that ictal and interictal source localization performed on the first stereoelectroencephalography recordings matches the findings of the second stereo-EEG exploration. Our findings suggest that independent component analysis followed by source localization on the topographies of interest is a promising method for retrieving the epileptogenic zone in case of suboptimal implantation.

Visualization of the Spatiotemporal Propagation of Interictal Spikes in Temporal Lobe Epilepsy: A MEG Pilot Study

Magnetoencephalography (MEG) is clinically used to localize interictal spikes in discrete brain areas of epilepsy patients through the equivalent current dipole (ECD) method, but does not account for the temporal dynamics of spike activity. Recent studies found that interictal spike propagation beyond the temporal lobe may be associated with worse postsurgical outcomes, but studies using whole-brain data such as in MEG remain limited. In this pilot study, we developed a tool that visualizes the spatiotemporal dynamics of interictal MEG spikes normalized to spike-free sleep activity to assess their onset and propagation patterns in patients with temporal lobe epilepsy (TLE). We extracted interictal source data containing focal epileptiform activity in awake and asleep states from seven patients whose MEG ECD clusters localized to the temporal lobe and normalized the data against spike-free sleep recordings. We calculated the normalized activity over time per cortical label, confirmed maximal activity at onset, and mapped the activity over a 10 ms interval onto each patient's brain using a custom-built Multi-Modal Visualization Tool. The onset of activity in all patients appeared near the clinically determined epileptogenic zone. By 10 ms, four of the patients had propagated source activity restricted to within the temporal lobe, and three had propagated source activity spread to extratemporal regions. Using this tool, we show that noninvasively identifying the onset and propagation of interictal spike activity in MEG can be achieved, which may help provide further insight into epileptic networks and guide surgical planning and interventions in patients with TLE.

Disentanglement of Resting State Brain Networks for Localizing Epileptogenic Zone in Focal Epilepsy

The objective of this study is to extract pathological brain networks from interictal period of E/MEG recordings to localize epileptic foci for presurgical evaluation. We proposed here a resting state E/MEG analysis framework, to disentangle brain functional networks represented by neural oscillations. By using an Embedded Hidden Markov Model, we constructed a state space for resting state recordings consisting of brain states with different spatiotemporal patterns. Functional connectivity analysis along with graph theory was applied on the extracted brain states to quantify the network features of the extracted brain states, based on which the source location of pathological states is determined. The method is evaluated by computer simulations and our simulation results revealed the proposed framework can extract brain states with high accuracy regarding both spatial and temporal profiles. We further evaluated the framework as compared with intracranial EEG defined seizure onset zone in 10 patients with drug-resistant focal epilepsy who underwent MEG recordings and were seizure free after surgical resection. The real patient data analysis showed very good localization results using the extracted pathological brain states in 6/10 patients, with localization error of about 15 mm as compared to the seizure onset zone. We show that the pathological brain networks can be disentangled from the resting-state electromagnetic recording and could be identified based on the connectivity features. The framework can serve as a useful tool in extracting brain functional networks from noninvasive resting state electromagnetic recordings, and promises to offer an alternative to aid presurgical evaluation guiding intracranial EEG electrodes implantation.

Eavesdropping on Tinnitus Using MEG: Lessons Learned and Future Perspectives

Tinnitus has been widely investigated in order to draw conclusions about the underlying causes and altered neural activity in various brain regions. Existing studies have based their work on different tinnitus frameworks, ranging from a more local perspective on the auditory cortex to the inclusion of broader networks and various approaches towards tinnitus perception and distress. Magnetoencephalography (MEG) provides a powerful tool for efficiently investigating tinnitus and aberrant neural activity both spatially and temporally. However, results are inconclusive, and studies are rarely mapped to theoretical frameworks. The purpose of this review was to firstly introduce MEG to interested researchers and secondly provide a synopsis of the current state. We divided recent tinnitus research in MEG into study designs using resting state measurements and studies implementing tone stimulation paradigms. The studies were categorized based on their theoretical foundation, and we outlined shortcomings as well as inconsistencies within the different approaches. Finally, we provided future perspectives on how to benefit more efficiently from the enormous potential of MEG. We suggested novel approaches from a theoretical, conceptual, and methodological point of view to allow future research to obtain a more comprehensive understanding of tinnitus and its underlying processes.

Helium Optically Pumped Magnetometers Can Detect Epileptic Abnormalities as Well as SQUIDs as Shown by Intracerebral Recordings

Magnetoencephalography based on superconducting quantum interference devices (SQUIDs) has been shown to improve the diagnosis and surgical treatment decision for presurgical evaluation of drug-resistant epilepsy. Still, its use remains limited because of several constraints such as cost, fixed helmet size, and the obligation of immobility. A new generation of sensors, optically pumped magnetometers (OPMs), could overcome these limitations. In this study, we validate the ability of helium-based OPM (4He-OPM) sensors to record epileptic brain activity thanks to simultaneous recordings with intracerebral EEG [stereotactic EEG (SEEG)]. We recorded simultaneous SQUIDs-SEEG and 4He-OPM-SEEG signals in one patient during two sessions. We show that epileptic activities on intracerebral EEG can be recorded by OPMs with a better signal-to noise ratio than classical SQUIDs. The OPM sensors open new venues for the widespread application of magnetoencephalography in the management of epilepsy and other neurologic diseases and fundamental neuroscience.

Connectivity alteration in thalamic nuclei and default mode network-related area in memory processes in mesial temporal lobe epilepsy using magnetoencephalography

This work aimed to investigate the involvement of the thalamic nuclei in mesial temporal lobe epilepsy (MTLE) and identify the influence of interictal epileptic discharges on the neural basis of memory processing by evaluating the functional connectivity (FC) between the thalamic nuclei and default mode network-related area (DMNRA) using magnetoencephalography. Preoperative datasets of nine patients with MTLE with seizure-free status after surgery and those of nine healthy controls were analyzed. The FC between the thalamic nuclei (anterior nucleus [ANT], mediodorsal nucleus [MD], intralaminar nuclei [IL]), hippocampus, and DMNRA was examined for each of the resting, pre-spike, spike, and post-spike periods in the delta to ripple bands using magnetoencephalography. The FC between the ANT, MD, hippocampus, and medial prefrontal cortex increased in the gamma to ripple bands, whereas the FC between the ANT, IL, and DMNRA decreased in the delta to beta bands, compared with that of the healthy controls at rest. Compared with the rest period, the pre-spike period had significantly decreased FC between the ANT, MD, and DMNRA in the ripple band. Different FC changes between the thalamic nuclei, hippocampus, and DMNRA of specific connections in a particular band may reflect impairment or compensation in the memory processes.

Scalable Orthonormal Projective NMF via Diversified Stochastic Optimization

The increasing availability of large-scale neuroimaging initiatives opens exciting opportunities for discovery science of human brain structure and function. Data-driven techniques, such as Orthonormal Projective Non-negative Matrix Factorization (opNMF), are well positioned to explore multivariate relationships in big data towards uncovering brain organization. opNMF enjoys advantageous interpretability and reproducibility compared to commonly used matrix factorization methods like Principal Component Analysis (PCA) and Independent Component Analysis (ICA), which led to its wide adoption in clinical computational neuroscience. However, applying opNMF in large-scale cohort studies is hindered by its limited scalability caused by its accompanying computational complexity. In this work, we address the computational challenges of opNMF using a stochastic optimization approach that learns over mini-batches of the data. Additionally, we diversify the stochastic batches via repulsive point processes, which reduce redundancy in the mini-batches and in turn lead to lower variance in the updates. We validated our framework on gray matter tissue density maps estimated from 1000 subjects part of the Open Access Series of Imaging (OASIS) dataset. We demonstrated that operations over mini-batches of data yield significant reduction in computational cost. Importantly, we showed that our novel optimization does not compromise the accuracy or interpretability of factors when compared to standard opNMF. The proposed model enables new investigations of brain structure using big neuroimaging data that could improve our understanding of brain structure in health and disease.

Test and Analysis of High-Permeability Material's Microstructure in Magnetic Shielding Device

The magnetic shielding device is used to provide an extreme weak magnetic field, which plays a key role in variety of fields. Since the high-permeability material constituting the magnetic shielding device determines the magnetic shielding performance, it is important to evaluate the property of the high-permeability material. In this paper, the relationship between the microstructure and the magnetic properties of the high-permeability material is analyzed using minimum free energy principle based on magnetic domain theory, and the test method of the material's microstructure including the material composition, the texture and the grain structure to reflect the magnetic properties is put forward. The test result shows that the grain structure is closely related to the initial permeability and the coercivity, which is highly consistent with the theory. As a result, it provides a more efficient way to evaluate the property of the high-permeability material. The test method proposed in the paper has important significance in the high efficiency sampling inspection of the high-permeability material.

Cortical and Subcortical Changes in MEG Activity Reflect Parkinson's Progression over a Period of 7 Years

In this study of early functional changes in Parkinson's disease (PD), we aimed to provide a comprehensive assessment of the development of changes in both cortical and subcortical neurophysiological brain activity, including their association with clinical measures of disease severity. Repeated resting-state MEG recordings and clinical assessments were obtained in the context of a unique longitudinal cohort study over a seven-year period using a multiple longitudinal design. We used linear mixed-models to analyze the relationship between neurophysiological (spectral power and functional connectivity) and clinical data. At baseline, early-stage (drug-naïve) PD patients demonstrated spectral slowing compared to healthy controls in both subcortical and cortical brain regions, most outspoken in the latter. Over time, spectral slowing progressed in strong association with clinical measures of disease progression (cognitive and motor). Global functional connectivity was not different between groups at baseline and hardly changed over time. Therefore, investigation of associations with clinical measures of disease progression were not deemed useful. An analysis of individual connections demonstrated differences between groups at baseline (higher frontal theta, lower parieto-occipital alpha2 band functional connectivity) and over time in PD patients (increase in frontal delta and theta band functional connectivity). Our results suggest that spectral measures are promising candidates in the search for non-invasive markers of both early-stage PD and of the ongoing disease process.

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.

Neurophysiology, Neuropsychology, and Epilepsy, in 2022: Hills We Have Climbed and Hills Ahead. Neurophysiology in epilepsy

Since the discovery of the human electroencephalogram (EEG), neurophysiology techniques have become indispensable tools in our armamentarium to localize epileptic seizures. New signal analysis techniques and the prospects of artificial intelligence and big data will offer unprecedented opportunities to further advance the field in the near future, ultimately resulting in improved quality of life for many patients with drug-resistant epilepsy. This article summarizes selected presentations from Day 1 of the two-day symposium "Neurophysiology, Neuropsychology, Epilepsy, 2022: Hills We Have Climbed and the Hills Ahead". Day 1 was dedicated to highlighting and honoring the work of Dr. Jean Gotman, a pioneer in EEG, intracranial EEG, simultaneous EEG/ functional magnetic resonance imaging, and signal analysis of epilepsy. The program focused on two main research directions of Dr. Gotman, and was dedicated to "High-frequency oscillations, a new biomarker of epilepsy" and "Probing the epileptic focus from inside and outside". All talks were presented by colleagues and former trainees of Dr. Gotman. The extended summaries provide an overview of historical and current work in the neurophysiology of epilepsy with emphasis on novel EEG biomarkers of epilepsy and source imaging and concluded with an outlook on the future of epilepsy research, and what is needed to bring the field to the next level.

Biomagnetism: The First Sixty Years

Biomagnetism is the measurement of the weak magnetic fields produced by nerves and muscle. The magnetic field of the heart-the magnetocardiogram (MCG)-is the largest biomagnetic signal generated by the body and was the first measured. Magnetic fields have been detected from isolated tissue, such as a peripheral nerve or cardiac muscle, and these studies have provided insights into the fundamental properties of biomagnetism. The magnetic field of the brain-the magnetoencephalogram (MEG)-has generated much interest and has potential clinical applications to epilepsy, migraine, and psychiatric disorders. The biomagnetic inverse problem, calculating the electrical sources inside the brain from magnetic field recordings made outside the head, is difficult, but several techniques have been introduced to solve it. Traditionally, biomagnetic fields are recorded using superconducting quantum interference device (SQUID) magnetometers, but recently, new sensors have been developed that allow magnetic measurements without the cryogenic technology required for SQUIDs.

Measurement and Analysis of Magnetic Properties of Permalloy for Magnetic Shielding Devices under Different Temperature Environments

The relative permeability, coercivity, and remanence of permalloy are closely related to the performance of magnetic shielding devices. In this paper, the relationship between the magnetic properties of permalloy and the working temperature of magnetic shielding devices is measured. Firstly, the measurement method of permalloy properties based on the simulated impact method is analyzed. What is more, a magnetic property test system consisting of a soft magnetic material tester and a high-low temperature chamber for permalloy ring samples at different temperatures was established to measure DC and AC (0.01 Hz to 1 kHz) magnetic properties at different temperatures (-60 °C to 140 °C). Finally, the results show that compared with room temperature (25 °C), the initial permeability (μi) decreases by 69.64% at -60 °C and increases by 38.23% at 140 °C, and the coercivity (hc) decreases by 34.81% at -60 °C and increases by 8.93% at 140 °C, which are the key parameters in the magnetic shielding device. It can be concluded that the relative permeability and remanence of permalloy are positively correlated with temperature, while the saturation magnetic flux density and coercivity are negatively correlated with temperature. This paper is of great significance to the magnetic analysis and design of magnetic shielding devices.

Studying memory processes at different levels with simultaneous depth and surface EEG recordings

Investigating cognitive brain functions using non-invasive electrophysiology can be challenging due to the particularities of the task-related EEG activity, the depth of the activated brain areas, and the extent of the networks involved. Stereoelectroencephalographic (SEEG) investigations in patients with drug-resistant epilepsy offer an extraordinary opportunity to validate information derived from non-invasive recordings at macro-scales. The SEEG approach can provide brain activity with high spatial specificity during tasks that target specific cognitive processes (e.g., memory). Full validation is possible only when performing simultaneous scalp SEEG recordings, which allows recording signals in the exact same brain state. This is the approach we have taken in 12 subjects performing a visual memory task that requires the recognition of previously viewed objects. The intracranial signals on 965 contact pairs have been compared to 391 simultaneously recorded scalp signals at a regional and whole-brain level, using multivariate pattern analysis. The results show that the task conditions are best captured by intracranial sensors, despite the limited spatial coverage of SEEG electrodes, compared to the whole-brain non-invasive recordings. Applying beamformer source reconstruction or independent component analysis does not result in an improvement of the multivariate task decoding performance using surface sensor data. By analyzing a joint scalp and SEEG dataset, we investigated whether the two types of signals carry complementary information that might improve the machine-learning classifier performance. This joint analysis revealed that the results are driven by the modality exhibiting best individual performance, namely SEEG.

Reconstruction and localization of auditory sources from intracerebral SEEG using independent component analysis

Stereo-electroencephalography (SEEG) is the surgical implantation of electrodes in the brain to better localize the epileptic network in pharmaco-resistant epileptic patients. This technique has exquisite spatial and temporal resolution. Still, the number and the position of the electrodes in the brain is limited and determined by the semiology and/or preliminary non-invasive examinations, leading to a large number of unexplored brain structures in each patient. Here, we propose a new approach to reconstruct the activity of non-sampled structures in SEEG, based on independent component analysis (ICA) and dipole source localization. We have tested this approach with an auditory stimulation dataset in ten patients. The activity directly recorded from the auditory cortex served as ground truth and was compared to the ICA applied on all non-auditory electrodes. Our results show that the activity from the auditory cortex can be reconstructed at the single trial level from contacts as far as ∼40 mm from the source. Importantly, this reconstructed activity is localized via dipole fitting in the proximity of the original source. In addition, we show that the size of the confidence interval of the dipole fitting is a good indicator of the reliability of the result, which depends on the geometry of the SEEG implantation. Overall, our approach allows reconstructing the activity of structures far from the electrode locations, partially overcoming the spatial sampling limitation of intracerebral recordings.

Task matters: Individual MEG signatures from naturalistic and neurophysiological brain states

The discovery that human brain connectivity data can be used as a "fingerprint" to identify a given individual from a population, has become a burgeoning research area in the neuroscience field. Recent studies have identified the possibility to extract these brain signatures from the temporal rich dynamics of resting-state magneto encephalography (MEG) recordings. Nevertheless, it is still uncertain to what extent MEG signatures can serve as an indicator of human identifiability during task-related conduct. Here, using MEG data from naturalistic and neurophysiological tasks, we show that identification improves in tasks relative to resting-state, providing compelling evidence for a task dependent axis of MEG signatures. Notably, improvements in identifiability were more prominent in strictly controlled tasks. Lastly, the brain regions contributing most towards individual identification were also modified when engaged in task activities. We hope that this investigation advances our understanding of the driving factors behind brain identification from MEG signals.

Extracting the Invisible: Mesial Temporal Source Detection in Simultaneous EEG and SEEG Recordings

Epileptic source detection relies mainly on visual expertise of scalp EEG signals, but it is recognised that epileptic discharges can escape to this expertise due to a deep localization of the brain sources that induce a very low, even negative, signal to noise ratio. In this methodological study, we aimed to investigate the feasibility of extracting deep mesial temporal sources that were invisible in scalp EEG signals using blind source separation (BSS) methods (infomax ICA, extended infomax ICA, and JADE) combined with a statistical measure (kurtosis). We estimated the effect of different methodological and physiological parameters that could alter or improve the extraction. Using nine well-defined mesial epileptic networks (1949 spikes) obtained from seven patients and simultaneous EEG-SEEG recordings, the first independent component extracted from the scalp EEG signals was validated in mean from 46 to 80% according to the different parameters. The three BSS methods equally performed (no significant difference) and no influence of the number of scalp electrodes used was found. At the opposite, the number and amplitude of spikes included in the averaging before the extraction modified the performance. Anyway, despite their invisibility in scalp EEG signals, this study demonstrates that deep source extraction is feasible under certain conditions and with the use of common signal analysis toolboxes. This finding confirms the crucial need to continue the signal analysis of scalp EEG recordings which contains subcortical signals that escape to expert visual analysis but could be found by signal processing.

Novel noninvasive identification of patient-specific epileptic networks in focal epilepsies: Linking single-photon emission computed tomography perfusion during seizures with resting-state magnetoencephalography dynamics

Single-photon emission computed tomography (SPECT) during seizures and magnetoencephalography (MEG) during the interictal state are noninvasive modalities employed in the localization of the epileptogenic zone in patients with drug-resistant focal epilepsy (DRFE). The present study aims to investigate whether there exists a preferentially high MEG functional connectivity (FC) among those regions of the brain that exhibit hyperperfusion or hypoperfusion during seizures. We studied MEG and SPECT data in 30 consecutive DRFE patients who had resective epilepsy surgery. We parcellated each ictal perfusion map into 200 regions of interest (ROIs) and generated ROI time series using source modeling of MEG data. FC between ROIs was quantified using coherence and phase-locking value. We defined a generalized linear model to relate the connectivity of each ROI, ictal perfusion z score, and distance between ROIs. We compared the coefficients relating perfusion z score to FC of each ROI and estimated the connectivity within and between resected and unresected ROIs. We found that perfusion z scores were strongly correlated with the FC of hyper-, and separately, hypoperfused ROIs across patients. High interictal connectivity was observed between hyperperfused brain regions inside and outside the resected area. High connectivity was also observed between regions of ictal hypoperfusion. Importantly, the ictally hypoperfused regions had a low interictal connectivity to regions that became hyperperfused during seizures. We conclude that brain regions exhibiting hyperperfusion during seizures highlight a preferentially connected interictal network, whereas regions of ictal hypoperfusion highlight a separate, discrete and interconnected, interictal network.

Manifestation of Hippocampal Interictal Discharges on Clinical Scalp EEG Recordings

PURPOSE

Epileptiform activity limited to deep sources such as the hippocampus currently lacks reliable scalp correlates. Recent studies, however, have found that a subset of hippocampal interictal discharges may be associated with visible scalp signals, suggesting that some types of hippocampal activity may be monitored noninvasively. The purpose of this study is to characterize the relationship between these scalp waveforms and the underlying intracranial activity.

METHODS

Paired intracranial and scalp EEG recordings obtained from 16 patients were used to identify hippocampal interictal discharges. Discharges were grouped by waveform shape, and spike-triggered averages of the intracranial and scalp signals were calculated for each group. Cross-correlation of intracranial and scalp spike-triggered averages was used to determine their temporal relationship, and topographic maps of the scalp were generated for each group.

RESULTS

Cross-correlation of intracranial and scalp correlates resulted in two classes of scalp waveforms-those with and without time delays from the associated hippocampal discharges. Scalp signals with no delay showed topographies with a broad field with higher amplitudes on the side ipsilateral to the discharges and a left-right flip in polarity-observations consistent with the volume conduction of a single unilateral deep source. In contrast, scalp correlates with time lags showed rotational dynamics, suggesting synaptic propagation mechanisms.

CONCLUSIONS

The temporal relationship between the intracranial and scalp signals suggests that both volume conduction and synaptic propagation contribute to these scalp manifestations. Furthermore, the topographic evolution of these scalp waveforms may be used to distinguish spikes that are limited to the hippocampus from those that travel to or engage other brain areas.

Hampered gamma oscillations induced by sad emotion underlying suicide attempt in major depressive disorder

AIM

Major depressive disorder (MDD) is associated with high suicidality, especially for those with suicide attempt (SA). Although impaired oscillatory activity has been previously reported in patients with SA, little is known about precise temporal-spatial variability of its neural dynamics. To solve this, the current study probed the spectral power and network interactions underlying SA in MDD.

METHODS

The present study recruited 104 subjects including 56 subjects with MDD (30 with SA and 26 without SA) and 48 healthy controls, who performed sad expressions recognition task during magnetoencephalography (MEG) recording. By investigating source-reconstructed MEG-data, brain states representing different task stages were estimated from a Hidden Markov model. Spectrum power and network connectivity were compared via Gaussian Mixture Models, and fractional occupancy (FO) of states were compared via an independent F-test.

RESULTS

Brain states were corresponding to various frequencies (theta/beta/low gamma/ high gamma). In low gamma band (35-45 Hz), the early visual state exhibited increased activation and hyper inter-network connectivity between visual regions and the limbic system, while the middle fronto-parietal state exhibited attenuated activation and decreased intra-network connectivity within fronto-parietal regions in SA group. Crucially, FO values of these two states were significantly correlated with the suicide risks.

CONCLUSIONS

Suicide behavior of patients with MDD was significantly associated with aberrant oscillations in low gamma band. Elevated oscillations in occipital cortices and attenuated oscillations in fronto-parietal cortices were significantly associated with SA. Manifesting sadness indulging and reckless decision-making, the hampered temporal characteristics could help explain the neural-electric basis of SA.

Comparison of beamformer and ICA for dynamic connectivity analysis: A simultaneous MEG-SEEG study

Magnetoencephalography (MEG) is a powerful tool for estimating brain connectivity with both good spatial and temporal resolution. It is particularly helpful in epilepsy to characterize non-invasively the epileptic networks. However, using MEG to map brain networks requires solving a difficult inverse problem that introduces uncertainty in the activity localization and connectivity measures. Our goal here was to compare independent component analysis (ICA) followed by dipole source localization and the linearly constrained minimum-variance beamformer (LCMV-BF) for characterizing regions with interictal epileptic activity and their dynamic connectivity. After a simulation study, we compared ICA and LCMV-BF results with intracerebral EEG (stereotaxic EEG, SEEG) recorded simultaneously in 8 epileptic patients, which provide a unique 'ground truth' to which non-invasive results can be confronted. We compared the signal time courses extracted applying ICA and LCMV-BF on MEG data to that of SEEG, both for the actual signals and the dynamic connectivity computed using cross-correlation (evolution of links in time). With our simulations, we illustrated the different effect of the temporal and spatial correlation among sources on the two methods. While ICA was more affected by the temporal correlation but robust against spatial configurations, LCMV-BF showed opposite behavior. Moreover, ICA seems more suited to retrieve the simulated networks. In case of real patient data, good MEG/SEEG correlation and good localization were obtained in 6 out of 8 patients. In 4 of them ICA had the best performance (higher correlation, lower localization distance). In terms of dynamic connectivity, the evolution in time of the cross-correlation links could be retrieved in 5 patients out of 6, however, with more variable results in terms of correlation and distance. In two patients LCMV-BF had better results than ICA. In one patient the two methods showed equally good outcomes, and in the remaining two patients ICA performed best. In conclusion, our results obtained by exploiting simultaneous MEG/SEEG recordings suggest that ICA and LCMV-BF have complementary qualities for retrieving the dynamics of interictal sources and their network interactions.

Preliminary Observations of Resting-State Magnetoencephalography in Nonmedicated Children with Obsessive-Compulsive Disorder

Background: Cortico-striato-thalamo-cortical (CSTC) network alterations are hypothesized to contribute to symptoms of obsessive-compulsive disorder (OCD). To date, very few studies have examined whether CSTC network alterations are present in children with OCD, who are medication naive. Medication-naive pediatric imaging samples may be optimal to study neural correlates of illness and identify brain-based markers, given the proximity to illness onset. Methods: Magnetoencephalography (MEG) data were analyzed at rest, in 18 medication-naive children with OCD (M = 12.1 years ±2.0 standard deviation [SD]; 10 M/8 F) and 13 typically developing children (M = 12.3 years ±2.2 SD; 6 M/7 F). Whole-brain MEG-derived resting-state functional connectivity (rs-fc), for alpha- and gamma-band frequencies were compared between OCD and typically developing (control) groups. Results: Increased MEG-derived rs-fc across alpha- and gamma-band frequencies was found in the OCD group compared to the control group. Increased MEG-derived rs-fc at alpha-band frequencies was evident across a number of regions within the CSTC circuitry and beyond, including the cerebellum and limbic regions. Increased MEG-derived rs-fc at gamma-band frequencies was restricted to the frontal and temporal cortices. Conclusions: This MEG study provides preliminary evidence of altered alpha and gamma networks, at rest, in medication-naive children with OCD. These results support prior findings pointing to the relevance of CSTC circuitry in pediatric OCD and further support accumulating evidence of altered connectivity between regions that extend beyond this network, including the cerebellum and limbic regions. Given the substantial portion of children and youth whose OCD symptoms do not respond to conventional treatments, our findings have implications for future treatment innovation research aiming to target and track whether brain patterns associated with having OCD may change with treatment and/or predict treatment response.

Virtual MEG sensors based on beamformer and independent component analysis can reconstruct epileptic activity as measured on simultaneous intracerebral recordings

The prevailing gold standard for presurgical determination of epileptogenic brain networks is intracerebral EEG, a potent yet invasive approach. Magnetoencephalography (MEG) is a state-of-the art non-invasive method for investigating epileptiform discharges. However, it is not clear at what level the precision offered by MEG can reach that of SEEG. Here, we present a strategy for non-invasively retrieving the constituents of the interictal network, with high spatial and temporal precision. Our method is based on MEG and a combination of spatial filtering and independent component analysis (ICA). We validated this approach in twelve patients with drug-resistant focal epilepsy, thanks to the unprecedented ground truth provided by simultaneous recordings of MEG and SEEG. A minimum variance adaptive beamformer estimated the source time series and ICA was used to further decompose these time series into network constituents (MEG-ICs), each having a time series (virtual electrode) and a topography (spatial distribution of amplitudes in the brain). We show that MEG has a considerable sensitivity of 0.80 and 0.84 and a specificity of 0.93 and 0.91 for reconstructing deep and superficial sources, respectively, when compared to the ground truth (SEEG). For each epileptic MEG-IC (n = 131), we found at least one significantly correlating SEEG contact close to zero lag after correcting for multiple comparisons. All the patients except one had at least one epileptic component that was highly correlated (Spearman rho>0.3) with that of SEEG traces. MEG-ICs correlated well with SEEG traces. The strength of correlation coefficients did not depend on the depth of the SEEG contacts or the clinical outcome of the patient. A significant proportion of the MEG-ICs (n = 83/131) were localized in proximity with their maximally correlating SEEG, within a mean distance of 20±12.18mm. Our research is the first to validate the MEG-retrieved beamformer IC sources against SEEG-derived ground truth in a simultaneous MEG-SEEG framework. Observations from the present study suggest that non-invasive MEG source components may potentially provide additional information, comparable to SEEG in a number of instances.

The hearing hippocampus

The hippocampus has a well-established role in spatial and episodic memory but a broader function has been proposed including aspects of perception and relational processing. Neural bases of sound analysis have been described in the pathway to auditory cortex, but wider networks supporting auditory cognition are still being established. We review what is known about the role of the hippocampus in processing auditory information, and how the hippocampus itself is shaped by sound. In examining imaging, recording, and lesion studies in species from rodents to humans, we uncover a hierarchy of hippocampal responses to sound including during passive exposure, active listening, and the learning of associations between sounds and other stimuli. We describe how the hippocampus' connectivity and computational architecture allow it to track and manipulate auditory information - whether in the form of speech, music, or environmental, emotional, or phantom sounds. Functional and structural correlates of auditory experience are also identified. The extent of auditory-hippocampal interactions is consistent with the view that the hippocampus makes broad contributions to perception and cognition, beyond spatial and episodic memory. More deeply understanding these interactions may unlock applications including entraining hippocampal rhythms to support cognition, and intervening in links between hearing loss and dementia.

Good scientific practice in EEG and MEG research: Progress and perspectives

Good scientific practice (GSP) refers to both explicit and implicit rules, recommendations, and guidelines that help scientists to produce work that is of the highest quality at any given time, and to efficiently share that work with the community for further scrutiny or utilization. For experimental research using magneto- and electroencephalography (MEEG), GSP includes specific standards and guidelines for technical competence, which are periodically updated and adapted to new findings. However, GSP also needs to be regularly revisited in a broader light. At the LiveMEEG 2020 conference, a reflection on GSP was fostered that included explicitly documented guidelines and technical advances, but also emphasized intangible GSP: a general awareness of personal, organizational, and societal realities and how they can influence MEEG research. This article provides an extensive report on most of the LiveMEEG contributions and new literature, with the additional aim to synthesize ongoing cultural changes in GSP. It first covers GSP with respect to cognitive biases and logical fallacies, pre-registration as a tool to avoid those and other early pitfalls, and a number of resources to enable collaborative and reproducible research as a general approach to minimize misconceptions. Second, it covers GSP with respect to data acquisition, analysis, reporting, and sharing, including new tools and frameworks to support collaborative work. Finally, GSP is considered in light of ethical implications of MEEG research and the resulting responsibility that scientists have to engage with societal challenges. Considering among other things the benefits of peer review and open access at all stages, the need to coordinate larger international projects, the complexity of MEEG subject matter, and today's prioritization of fairness, privacy, and the environment, we find that current GSP tends to favor collective and cooperative work, for both scientific and for societal reasons.

Deep neural networks constrained by neural mass models improve electrophysiological source imaging of spatiotemporal brain dynamics

Many efforts have been made to image the spatiotemporal electrical activity of the brain with the purpose of mapping its function and dysfunction as well as aiding the management of brain disorders. Here, we propose a non-conventional deep learning-based source imaging framework (DeepSIF) that provides robust and precise spatiotemporal estimates of underlying brain dynamics from noninvasive high-density electroencephalography (EEG) recordings. DeepSIF employs synthetic training data generated by biophysical models capable of modeling mesoscale brain dynamics. The rich characteristics of underlying brain sources are embedded in the realistic training data and implicitly learned by DeepSIF networks, avoiding complications associated with explicitly formulating and tuning priors in an optimization problem, as often is the case in conventional source imaging approaches. The performance of DeepSIF is evaluated by 1) a series of numerical experiments, 2) imaging sensory and cognitive brain responses in a total of 20 healthy subjects from three public datasets, and 3) rigorously validating DeepSIF's capability in identifying epileptogenic regions in a cohort of 20 drug-resistant epilepsy patients by comparing DeepSIF results with invasive measurements and surgical resection outcomes. DeepSIF demonstrates robust and excellent performance, producing results that are concordant with common neuroscience knowledge about sensory and cognitive information processing as well as clinical findings about the location and extent of the epileptogenic tissue and outperforming conventional source imaging methods. The DeepSIF method, as a data-driven imaging framework, enables efficient and effective high-resolution functional imaging of spatiotemporal brain dynamics, suggesting its wide applicability and value to neuroscience research and clinical applications.

Advances in human intracranial electroencephalography research, guidelines and good practices

Since the second-half of the twentieth century, intracranial electroencephalography (iEEG), including both electrocorticography (ECoG) and stereo-electroencephalography (sEEG), has provided an intimate view into the human brain. At the interface between fundamental research and the clinic, iEEG provides both high temporal resolution and high spatial specificity but comes with constraints, such as the individual's tailored sparsity of electrode sampling. Over the years, researchers in neuroscience developed their practices to make the most of the iEEG approach. Here we offer a critical review of iEEG research practices in a didactic framework for newcomers, as well addressing issues encountered by proficient researchers. The scope is threefold: (i) review common practices in iEEG research, (ii) suggest potential guidelines for working with iEEG data and answer frequently asked questions based on the most widespread practices, and (iii) based on current neurophysiological knowledge and methodologies, pave the way to good practice standards in iEEG research. The organization of this paper follows the steps of iEEG data processing. The first section contextualizes iEEG data collection. The second section focuses on localization of intracranial electrodes. The third section highlights the main pre-processing steps. The fourth section presents iEEG signal analysis methods. The fifth section discusses statistical approaches. The sixth section draws some unique perspectives on iEEG research. Finally, to ensure a consistent nomenclature throughout the manuscript and to align with other guidelines, e.g., Brain Imaging Data Structure (BIDS) and the OHBM Committee on Best Practices in Data Analysis and Sharing (COBIDAS), we provide a glossary to disambiguate terms related to iEEG research.

Magnetoencephalography can reveal deep brain network activities linked to memory processes

Recording from deep neural structures such as hippocampus noninvasively and yet with high temporal resolution remains a major challenge for human neuroscience. Although it has been proposed that deep neuronal activity might be recordable during cognitive tasks using magnetoencephalography (MEG), this remains to be demonstrated as the contribution of deep structures to MEG recordings may be too small to be detected or might be eclipsed by the activity of large‐scale neocortical networks. In the present study, we disentangled mesial activity and large‐scale networks from the MEG signals thanks to blind source separation (BSS). We then validated the MEG BSS components using intracerebral EEG signals recorded simultaneously in patients during their presurgical evaluation of epilepsy. In the MEG signals obtained during a memory task involving the recognition of old and new images, we identified with BSS a putative mesial component, which was present in all patients and all control subjects. The time course of the component selectively correlated with stereo‐electroencephalography signals recorded from hippocampus and rhinal cortex, thus confirming its mesial origin. This finding complements previous studies with epileptic activity and opens new possibilities for using MEG to study deep brain structures in cognition and in brain disorders.

Cortical and subcortical changes in resting-state neuronal activity and connectivity in early symptomatic ALS and advanced frontotemporal dementia

The objective of this study was to examine if patterns of resting-state brain activity and functional connectivity in cortical and subcortical regions in patients with early symptomatic amyotrophic lateral sclerosis (ALS) resemble those of behavioural variant frontotemporal dementia (bvFTD). In a cross-sectional design, eyes-closed resting-state magnetoencephalography (MEG) data of 34 ALS patients, 18 bvFTD patients and 18 age- and gender-matched healthy controls (HCs) were projected to source-space using an atlas-based beamformer. Group differences in peak frequency, band-specific oscillatory activity and functional connectivity (corrected amplitude envelope correlation) in 78 cortical regions and 12 subcortical regions were determined. False discovery rate was used to correct for multiple comparisons. BvFTD patients, as compared to ALS and HCs, showed lower relative beta power in parietal, occipital, temporal and nearly all subcortical regions. Compared to HCs, patients with ALS and patients with bvFTD had a higher delta (0.5-4 Hz) and gamma (30-48 Hz) band resting-state functional connectivity in a high number of overlapping regions in the frontal lobe and in limbic and subcortical regions. Higher delta band connectivity was widespread in the bvFTD patients compared to HCs. ALS showed a more widespread higher gamma band functional connectivity compared to bvFTD. In conclusion, MEG in early symptomatic ALS patients shows resting-state functional connectivity changes in frontal, limbic and subcortical regions that overlap considerably with bvFTD. The findings show the potential of MEG to detect brain changes in early symptomatic phases of ALS and contribute to our understanding of the disease spectrum, with ALS and bvFTD at the two extreme ends.

Simulation Study of Different OPM-MEG Measurement Components

Magnetoencephalography (MEG) is a neuroimaging technique that measures the magnetic fields of the brain outside of the head. In the past, the most suitable magnetometer for MEG was the superconducting quantum interference device (SQUID), but in recent years, a new type has also been used, the optically pumped magnetometer (OPM). OPMs can be configured to measure multiple directions of magnetic field simultaneously. This work explored whether combining multiple directions of the magnetic field lowers the source localization error of brain sources under various conditions of noise. We simulated dipolar-like sources for multiple configurations of both SQUID- and OPM-MEG systems. To test the performance of a given layout, we calculated the average signal-to-noise ratio and the root mean square of the simulated magnetic field; furthermore, we evaluated the performance of the dipole fit. The results showed that the field direction normal to the scalp yields a higher signal-to-noise ratio and that ambient noise has a much lower impact on its localization error; therefore, this is the optimal choice for source localization when only one direction of magnetic field can be measured. For a low number of OPMs, combining multiple field directions greatly improves the source localization results. Lastly, we showed that MEG sensors that can be placed closer to the brain are more suitable for localizing deeper sources.

Oscillatory Activity of the Hippocampus in Prodromal Alzheimer’s Disease: A Source-Space Magnetoencephalography Study

Background: In Alzheimer’s disease (AD), oscillatory activity of the human brain slows down. However, oscillatory slowing varies between individuals, particularly in prodromal AD. Cortical oscillatory changes have shown suboptimal accuracy as diagnostic markers. We speculated that focusing on the hippocampus might prove more successful, particularly using magnetoencephalography (MEG) for capturing subcortical oscillatory activity. Objective: We explored MEG-based detection of hippocampal oscillatory abnormalities in prodromal AD patients. Methods: We acquired resting-state MEG data of 18 AD dementia patients, 18 amyloid-β-positive amnestic mild cognitive impairment (MCI, prodromal AD) patients, and 18 amyloid-β-negative persons with subjective cognitive decline (SCD). Oscillatory activity in 78 cortical regions and both hippocampi was reconstructed using beamforming. Between-group and hippocampal-cortical differences in spectral power were assessed. Classification accuracy was explored using ROC curves. Results: The MCI group showed intermediate power values between SCD and AD, except for the alpha range, where it was higher than both (p < 0.05 and p < 0.001). The largest differences between MCI and SCD were in the theta band, with higher power in MCI (p < 0.01). The hippocampi showed several unique group differences, such as higher power in the higher alpha band in MCI compared to SCD (p < 0.05). Classification accuracy (MCI versus SCD) was best for absolute theta band power in the right hippocampus (AUC = 0.87). Conclusion: In this MEG study, we detected oscillatory abnormalities of the hippocampi in prodromal AD patients. Moreover, hippocampus-based classification performed better than cortex-based classification. We conclude that a focus on hippocampal MEG may improve early detection of AD-related neuronal dysfunction.

Improving Detection of Hippocampal Epileptiform Activity Using Magnetoencephalography

PURPOSE

Magnetoencephalography (MEG) defines the spike-generating zone and provides targets for invasive monitoring with stereotactic electroencephalography. This retrospective, blinded, cross-sectional study determined whether MEG virtual sensors could identify hippocampal epileptiform activity.

METHODS

Using MEG beamformer analysis, virtual sensors were manually placed in bilateral hippocampi and corresponding virtual sensor waveforms were analyzed for the presence of epileptiform activity. These findings were compared with hippocampal stereotactic electroencephalography in the same patients. Concordance was determined using sensitivity and specificity.

RESULTS

Thirty patients (mean age 12.5 ± 5.9 years) and 35 hippocampi were included. Patients were also placed into subgroups based on conventional MEG analysis: temporal (n = 19), extratemporal (n = 10), and normal (n = 1). Overall, sensitivity and specificity were 57.9% and 50.0%, respectively (n = 35). Patients with temporal sources based on conventional MEG analysis had sensitivity and specificity of 80.0% and 36.4%, respectively (n = 21). Those with extratemporal sources based on conventional MEG had sensitivity and specificity of 42.9% and 80.0%, respectively (n = 12).

CONCLUSIONS

When grouped by conventional MEG analysis, virtual sensors can be useful to confirm mesial temporal dipoles seen with conventional analysis.

SIGNIFICANCE

This work may help support the use of MEG for the detection of epileptiform activity in the hippocampus and influence the planning of invasive electrode placement.

Sensitive and reproducible MEG resting-state metrics of functional connectivity in Alzheimer's disease

Background

Analysis of functional brain networks in Alzheimer’s disease (AD) has been hampered by a lack of reproducible, yet valid metrics of functional connectivity (FC). This study aimed to assess both the sensitivity and reproducibility of the corrected amplitude envelope correlation (AEC-c) and phase lag index (PLI), two metrics of FC that are insensitive to the effects of volume conduction and field spread, in two separate cohorts of patients with dementia due to AD versus healthy elderly controls.

Methods

Subjects with a clinical diagnosis of AD dementia with biomarker proof, and a control group of subjective cognitive decline (SCD), underwent two 5-min resting-state MEG recordings. Data consisted of a test (AD = 28; SCD = 29) and validation (AD = 29; SCD = 27) cohort. Time-series were estimated for 90 regions of interest (ROIs) in the automated anatomical labelling (AAL) atlas. For each of five canonical frequency bands, the AEC-c and PLI were calculated between all 90 ROIs, and connections were averaged per ROI. General linear models were constructed to compare the global FC differences between the groups, assess the reproducibility, and evaluate the effects of age and relative power. Reproducibility of the regional FC differences was assessed using the Mann-Whitney U tests, with correction for multiple testing using the false discovery rate (FDR).

Results

The AEC-c showed significantly and reproducibly lower global FC for the AD group compared to SCD, in the alpha (8–13 Hz) and beta (13–30 Hz) bands, while the PLI revealed reproducibly lower FC for the AD group in the delta (0.5–4 Hz) band and higher FC for the theta (4–8 Hz) band. Regionally, the beta band AEC-c showed reproducibility for almost all ROIs (except for 13 ROIs in the frontal and temporal lobes). For the other bands, the AEC-c and PLI did not show regional reproducibility after FDR correction. The theta band PLI was susceptible to the effect of relative power.

Conclusion

For MEG, the AEC-c is a sensitive and reproducible metric, able to distinguish FC differences between patients with AD dementia and cognitively healthy controls. These two measures likely reflect different aspects of neural activity and show differential sensitivity to changes in neural dynamics.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13195-022-00970-4.

EEG effective connectivity during the first year of life mirrors brain synaptogenesis, myelination, and early right hemisphere predominance

INTRODUCTION

The maturation of electroencephalogram (EEG) effective connectivity in healthy infants during the first year of life is described.

METHODS

Participants: A cross-sectional sample of 125 healthy at-term infants, from 0 to 12 months of age, underwent EEG in a state of quiet sleep.

PROCEDURES

The EEG primary currents at the source were described with the sLoreta method. An unmixing algorithm was applied to reduce the leakage, and the isolated effective coherence, a direct and directed measurement of information flow, was calculated.

RESULTS AND DISCUSSION

Initially, the highest indices of connectivity are at the subcortical nuclei, continuing to the parietal lobe, predominantly the right hemisphere, then expanding to temporal, occipital, and finally the frontal areas, which is consistent with the myelination process. Age-related connectivity changes were mostly long-range and bilateral. Connections increased with age, mainly in the right hemisphere, while they mainly decreased in the left hemisphere. Increased connectivity from 20 to 30 Hz, mostly at the right hemisphere. These findings were consistent with right hemisphere predominance during the first three years of life. Theta and alpha connections showed the greatest changes with age. Strong connectivity was found between the parietal, temporal, and occipital regions to the frontal lobes, responsible for executive functions and consistent with behavioral development during the first year. The thalamus exchanges information bidirectionally with all cortical regions and frequency bands.

CONCLUSIONS

The maturation of EEG connectivity during the first year in healthy infants is very consistent with synaptogenesis, reductions in synaptogenesis, myelination, and functional and behavioral development.