Frequency-dependent functional connectivity within resting-state networks: an atlas-based MEG beamformer solution

Functionally annotated electrophysiological neuromarkers of healthy ageing and memory function

The unprecedented increase in life expectancy presents a unique opportunity and the necessity to explore both healthy and pathological aspects of ageing. Electroencephalography (EEG) has been widely used to identify neuromarkers of cognitive ageing due to its affordability and richness in information. However, despite the growing volume of data and methodological advancements, the abundance of contradictory and non-reproducible findings has hindered clinical translation. To address these challenges, our study introduces a comprehensive workflow expanding on previous EEG studies and investigates various static and dynamic power and connectivity estimates as potential neuromarkers of cognitive ageing in a large dataset. We also assess the robustness of our findings by testing their susceptibility to band specification. Finally, we characterise our findings using functionally annotated brain networks to improve their interpretability and multi-modal integration. Our analysis demonstrates the effect of methodological choices on findings and that dynamic rather than static neuromarkers are not only more sensitive but also more robust. Consequently, they emerge as strong candidates for cognitive ageing neuromarkers. Moreover, we were able to replicate the most established EEG findings in cognitive ageing, such as alpha oscillation slowing, increased beta power, reduced reactivity across multiple bands, and decreased delta connectivity. Additionally, when considering individual variations in the alpha band, we clarified that alpha power is characteristic of memory performance rather than ageing, highlighting its potential as a neuromarker for cognitive ageing. Finally, our approach using functionally annotated source reconstruction allowed us to provide insights into domain-specific electrophysiological mechanisms underlying memory performance and ageing. HIGHLIGHTS: We provide an open and reproducible pipeline with a comprehensive workflow to investigate static and dynamic EEG neuromarkers. Neuromarkers related to neural dynamics are sensitive and robust. Individualised alpha power characterises cognitive performance rather than ageing. Functional annotation allows cross-modal interpretation of EEG findings.

Neurophysiological brain function predicts response to cognitive rehabilitation and mindfulness in multiple sclerosis: a randomized trial

BACKGROUND

Cognitive treatment response varies highly in people with multiple sclerosis (PwMS). Identification of mechanisms is essential for predicting response.

OBJECTIVES

This study aimed to investigate whether brain network function predicts response to cognitive rehabilitation therapy (CRT) and mindfulness-based cognitive therapy (MBCT).

METHODS

PwMS with cognitive complaints completed CRT, MBCT, or enhanced treatment as usual (ETAU) and performed three measurements (baseline, post-treatment, 6-month follow-up). Baseline magnetoencephalography (MEG) measures were used to predict treatment effects on cognitive complaints, personalized cognitive goals, and information processing speed (IPS) using mixed models (secondary analysis REMIND-MS study).

RESULTS

We included 105 PwMS (96 included in prediction analyses; 32 CRT, 31 MBCT, 33 ETAU), and 56 healthy controls with baseline MEG. MEG did not predict reductions in complaints. Higher connectivity predicted better goal achievement after MBCT (p = 0.010) and CRT (p = 0.018). Lower gamma power (p = 0.006) and higher connectivity (p = 0.020) predicted larger IPS benefits after MBCT. These MEG predictors indicated worse brain function compared to healthy controls (p < 0.05).

CONCLUSIONS

Brain network function predicted better cognitive goal achievement after MBCT and CRT, and IPS improvements after MBCT. PwMS with neuronal slowing and hyperconnectivity were most prone to show treatment response, making network function a promising tool for personalized treatment recommendations.

TRIAL REGISTRATION

The REMIND-MS study was prospectively registered in the Dutch Trial registry (NL6285; https://trialsearch.who.int/Trial2.aspx?TrialID=NTR6459 ).

Data-driven MEG analysis to extract fMRI resting-state networks

The electrophysiological basis of resting-state networks (RSN) is still under debate. In particular, no principled mechanism has been determined that is capable of explaining all RSN equally well. While magnetoencephalography (MEG) and electroencephalography are the methods of choice to determine the electrophysiological basis of RSN, no standard analysis pipeline of RSN yet exists. In this article, we compare the two main existing data-driven analysis strategies for extracting RSNs from MEG data and introduce a third approach. The first approach uses phase-amplitude coupling to determine the RSN. The second approach extracts RSN through an independent component analysis of the Hilbert envelope in different frequency bands, while the third new approach uses a singular value decomposition instead. To evaluate these approaches, we compare the MEG-RSN to the functional magnetic resonance imaging (fMRI)-RSN from the same subjects. Overall, it was possible to extract RSN with MEG using all three techniques, which matched the group-specific fMRI-RSN. Interestingly the new approach based on SVD yielded significantly higher correspondence to five out of seven fMRI-RSN than the two existing approaches. Importantly, with this approach, all networks-except for the visual network-had the highest correspondence to the fMRI networks within one frequency band. Thereby we provide further insights into the electrophysiological underpinnings of the fMRI-RSNs. This knowledge will be important for the analysis of the electrophysiological connectome.

Deep brain stimulation of the central thalamus restores arousal and motivation in a zolpidem-responsive patient with akinetic mutism after severe brain injury

After severe brain injury, zolpidem is known to cause spectacular, often short-lived, restorations of brain functions in a small subgroup of patients. Previously, we showed that these zolpidem-induced neurological recoveries can be paralleled by significant changes in functional connectivity throughout the brain. Deep brain stimulation (DBS) is a neurosurgical intervention known to modulate functional connectivity in a wide variety of neurological disorders. In this study, we used DBS to restore arousal and motivation in a zolpidem-responsive patient with severe brain injury and a concomitant disorder of diminished motivation, more than 10 years after surviving hypoxic ischemia. We found that DBS of the central thalamus, targeted at the centromedian-parafascicular complex, immediately restored arousal and was able to transition the patient from a state of deep sleep to full wakefulness. Moreover, DBS was associated with temporary restoration of communication and ability to walk and eat in an otherwise wheelchair-bound and mute patient. With the use of magnetoencephalography (MEG), we revealed that DBS was generally associated with a marked decrease in aberrantly high levels of functional connectivity throughout the brain, mimicking the effects of zolpidem. These results imply that 'pathological hyperconnectivity' after severe brain injury can be associated with reduced arousal and behavioral performance and that DBS is able to modulate connectivity towards a 'healthier baseline' with lower synchronization, and, can restore functional brain networks long after severe brain injury. The presence of hyperconnectivity after brain injury may be a possible future marker for a patient's responsiveness for restorative interventions, such as DBS, and suggests that lower degrees of overall brain synchronization may be conducive to cognition and behavioral responsiveness.

The relationship between pathological brain activity and functional network connectivity in glioma patients

PURPOSE

Glioma is associated with pathologically high (peri)tumoral brain activity, which relates to faster progression. Functional connectivity is disturbed locally and throughout the entire brain, associating with symptomatology. We, therefore, investigated how local activity and network measures relate to better understand how the intricate relationship between the tumor and the rest of the brain may impact disease and symptom progression.

METHODS

We obtained magnetoencephalography in 84 de novo glioma patients and 61 matched healthy controls. The offset of the power spectrum, a proxy of neuronal activity, was calculated for 210 cortical regions. We calculated patients' regional deviations in delta, theta and lower alpha network connectivity as compared to controls, using two network measures: clustering coefficient (local connectivity) and eigenvector centrality (integrative connectivity). We then tested group differences in activity and connectivity between (peri)tumoral, contralateral homologue regions, and the rest of the brain. We also correlated regional offset to connectivity.

RESULTS

As expected, patients' (peri)tumoral activity was pathologically high, and patients showed higher clustering and lower centrality than controls. At the group-level, regionally high activity related to high clustering in controls and patients alike. However, within-patient analyses revealed negative associations between regional deviations in brain activity and clustering, such that pathologically high activity coincided with low network clustering, while regions with 'normal' activity levels showed high network clustering.

CONCLUSION

Our results indicate that pathological activity and connectivity co-localize in a complex manner in glioma. This insight is relevant to our understanding of disease progression and cognitive symptomatology.

Functional brain alterations in symptomatic dermographism patients-An exploratory magnetoencephalography study

Symptomatic dermographism (SD) is a common form of urticaria, which is triggered by stroking the skin. Brain involvement in its aetiology was investigated by means of magnetoencephalography (MEG) after provocation with histamine and dermography. Wheals were induced by histamine skin prick test and dermography in twelve SD patients and fourteen controls. Itch severity was scored on a Visual Analogue Scale (VAS). Relative power and functional connectivity (FC) were measured using a 306-channel whole-head MEG system at baseline and 10 min after histamine and dermography, and contrasted between groups and conditions. Furthermore, wheal diameter and itch scores after these procedures were correlated with the MEG values. SD patients had higher itch scores after histamine and dermography. No significant group-differences were observed in relative power or FC for any condition. In both groups, power decreases were mostly observed in the beta band, and power increases in the alpha bands, after provocation, with more regions involved in patients compared to controls. Increased FC was seen after histamine in patients, and after dermography in controls. In patients only, dermography and histamine wheal size correlated with the alpha2 power in the regions of interest that showed significant condition effects after these procedures. Our findings may be cautiously interpreted as aberrant itch processing, and suggest involvement of the central nervous system in the aetiology of SD.

Solving large-scale MEG/EEG source localisation and functional connectivity problems simultaneously using state-space models

State-space models are widely employed across various research disciplines to study unobserved dynamics. Conventional estimation techniques, such as Kalman filtering and expectation maximisation, offer valuable insights but incur high computational costs in large-scale analyses. Sparse inverse covariance estimators can mitigate these costs, but at the expense of a trade-off between enforced sparsity and increased estimation bias, necessitating careful assessment in low signal-to-noise ratio (SNR) situations. To address these challenges, we propose a three-fold solution: (1) Introducing multiple penalised state-space (MPSS) models that leverage data-driven regularisation; (2) Developing novel algorithms derived from backpropagation, gradient descent, and alternating least squares to solve MPSS models; (3) Presenting a K-fold cross-validation extension for evaluating regularisation parameters. We validate this MPSS regularisation framework through lower and more complex simulations under varying SNR conditions, including a large-scale synthetic magneto- and electro-encephalography (MEG/EEG) data analysis. In addition, we apply MPSS models to concurrently solve brain source localisation and functional connectivity problems for real event-related MEG/EEG data, encompassing thousands of sources on the cortical surface. The proposed methodology overcomes the limitations of existing approaches, such as constraints to small-scale and region-of-interest analyses. Thus, it may enable a more accurate and detailed exploration of cognitive brain functions.

Cortical Synchrony and Information Flow during Transition from Wakefulness to Light Non-Rapid Eye Movement Sleep

Sleep is a highly stereotyped phenomenon, requiring robust spatiotemporal coordination of neural activity. Understanding how the brain coordinates neural activity with sleep onset can provide insights into the physiological functions subserved by sleep and the pathologic phenomena associated with sleep onset. We quantified whole-brain network changes in synchrony and information flow during the transition from wakefulness to light non-rapid eye movement (NREM) sleep, using MEG imaging in a convenient sample of 14 healthy human participants (11 female; mean 63.4 years [SD 11.8 years]). We furthermore performed computational modeling to infer excitatory and inhibitory properties of local neural activity. The transition from wakefulness to light NREM was identified to be encoded in spatially and temporally specific patterns of long-range synchrony. Within the delta band, there was a global increase in connectivity from wakefulness to light NREM, which was highest in frontoparietal regions. Within the theta band, there was an increase in connectivity in fronto-parieto-occipital regions and a decrease in temporal regions from wakefulness to Stage 1 sleep. Patterns of information flow revealed that mesial frontal regions receive hierarchically organized inputs from broad cortical regions upon sleep onset, including direct inflow from occipital regions and indirect inflow via parieto-temporal regions within the delta frequency band. Finally, biophysical neural mass modeling demonstrated changes in the anterior-to-posterior distribution of cortical excitation-to-inhibition with increased excitation-to-inhibition model parameters in anterior regions in light NREM compared with wakefulness. Together, these findings uncover whole-brain corticocortical structure and the orchestration of local and long-range, frequency-specific cortical interactions in the sleep-wake transition.SIGNIFICANCE STATEMENT Our work uncovers spatiotemporal cortical structure of neural synchrony and information flow upon the transition from wakefulness to light non-rapid eye movement sleep. Mesial frontal regions were identified to receive hierarchically organized inputs from broad cortical regions, including both direct inputs from occipital regions and indirect inputs via the parieto-temporal regions within the delta frequency range. Biophysical neural mass modeling revealed a spatially heterogeneous, anterior-posterior distribution of cortical excitation-to-inhibition. Our findings shed light on the orchestration of local and long-range cortical neural structure that is fundamental to sleep onset, and support an emerging view of cortically driven regulation of sleep homeostasis.

Connectivity-based Meta-Bands: A new approach for automatic frequency band identification in connectivity analyses

The majority of electroencephalographic (EEG) and magnetoencephalographic (MEG) studies filter and analyse neural signals in specific frequency ranges, known as "canonical" frequency bands. However, this segmentation, is not exempt from limitations, mainly due to the lack of adaptation to the neural idiosyncrasies of each individual. In this study, we introduce a new data-driven method to automatically identify frequency ranges based on the topological similarity of the frequency-dependent functional neural network. The resting-state neural activity of 195 cognitively healthy subjects from three different databases (MEG: 123 subjects; EEG1: 27 subjects; EEG2: 45 subjects) was analysed. In a first step, MEG and EEG signals were filtered with a narrow-band filter bank (1 Hz bandwidth) from 1 to 70 Hz with a 0.5 Hz step. Next, the connectivity in each of these filtered signals was estimated using the orthogonalized version of the amplitude envelope correlation to obtain the frequency-dependent functional neural network. Finally, a community detection algorithm was used to identify communities in the frequency domain showing a similar network topology. We have called this approach the "Connectivity-based Meta-Bands" (CMB) algorithm. Additionally, two types of synthetic signals were used to configure the hyper-parameters of the CMB algorithm. We observed that the classical approaches to band segmentation are partially aligned with the underlying network topologies at group level for the MEG signals, but they are missing individual idiosyncrasies that may be biasing previous studies, as revealed by our methodology. On the other hand, the sensitivity of EEG signals to reflect this underlying frequency-dependent network structure is limited, revealing a simpler frequency parcellation, not aligned with that defined by the "canonical" frequency bands. To the best of our knowledge, this is the first study that proposes an unsupervised band segmentation method based on the topological similarity of functional neural network across frequencies. This methodology fully accounts for subject-specific patterns, providing more robust and personalized analyses, and paving the way for new studies focused on exploring the frequency-dependent structure of brain connectivity.

Tau protein spreads through functionally connected neurons in Alzheimer's disease: a combined MEG/PET study

Recent studies on Alzheimer's disease (AD) suggest that tau proteins spread through the brain following neuronal connections. Several mechanisms could be involved in this process: spreading between brain regions that interact strongly (functional connectivity); through the pattern of anatomical connections (structural connectivity); or simple diffusion. Using magnetoencephalography (MEG), we investigated which spreading pathways influence tau protein spreading by modelling the tau propagation process using an epidemic spreading model. We compared the modelled tau depositions with 18F-flortaucipir PET binding potentials at several stages of the AD continuum. In this cross-sectional study, we analysed source-reconstructed MEG data and dynamic 100-min 18F-flortaucipir PET from 57 subjects positive for amyloid-β pathology [preclinical AD (n = 16), mild cognitive impairment (MCI) due to AD (n = 16) and AD dementia (n = 25)]. Cognitively healthy subjects without amyloid-β pathology were included as controls (n = 25). Tau propagation was modelled as an epidemic process (susceptible-infected model) on MEG-based functional networks [in alpha (8-13 Hz) and beta (13-30 Hz) bands], a structural or diffusion network, starting from the middle and inferior temporal lobe. The group-level network of the control group was used as input for the model to predict tau deposition in three stages of the AD continuum. To assess performance, model output was compared to the group-specific tau deposition patterns as measured with 18F-flortaucipir PET. We repeated the analysis by using networks of the preceding disease stage and/or using regions with most observed tau deposition during the preceding stage as seeds. In the preclinical AD stage, the functional networks predicted most of the modelled tau-PET binding potential, with best correlations between model and tau-PET [corrected amplitude envelope correlation (AEC-c) alpha C = 0.584; AEC-c beta C = 0.569], followed by the structural network (C = 0.451) and simple diffusion (C = 0.451). Prediction accuracy declined for the MCI and AD dementia stages, although the correlation between modelled tau and tau-PET binding remained highest for the functional networks (C = 0.384; C = 0.376). Replacing the control-network with the network from the preceding disease stage and/or alternative seeds improved prediction accuracy in MCI but not in the dementia stage. These results suggest that in addition to structural connections, functional connections play an important role in tau spread, and highlight that neuronal dynamics play a key role in promoting this pathological process. Aberrant neuronal communication patterns should be taken into account when identifying targets for future therapy. Our results also suggest that this process is more important in earlier disease stages (preclinical AD/MCI); possibly, in later stages, other processes may be influential.

Comparing individual and group-level simulated neurophysiological brain connectivity using the Jansen and Rit neural mass model

Computational models are often used to assess how functional connectivity (FC) patterns emerge from neuronal population dynamics and anatomical brain connections. It remains unclear whether the commonly used group-averaged data can predict individual FC patterns. The Jansen and Rit neural mass model was employed, where masses were coupled using individual structural connectivity (SC). Simulated FC was correlated to individual magnetoencephalography-derived empirical FC. FC was estimated using phase-based (phase lag index (PLI), phase locking value (PLV)), and amplitude-based (amplitude envelope correlation (AEC)) metrics to analyze their goodness of fit for individual predictions. Individual FC predictions were compared against group-averaged FC predictions, and we tested whether SC of a different participant could equally well predict participants' FC patterns. The AEC provided a better match between individually simulated and empirical FC than phase-based metrics. Correlations between simulated and empirical FC were higher using individual SC compared to group-averaged SC. Using SC from other participants resulted in similar correlations between simulated and empirical FC compared to using participants' own SC. This work underlines the added value of FC simulations using individual instead of group-averaged SC for this particular computational model and could aid in a better understanding of mechanisms underlying individual functional network trajectories.

Neurophysiological alterations in mice and humans carrying mutations in APP and PSEN1 genes

BACKGROUND

Studies in animal models of Alzheimer's disease (AD) have provided valuable insights into the molecular and cellular processes underlying neuronal network dysfunction. Whether and how AD-related neurophysiological alterations translate between mice and humans remains however uncertain.

METHODS

We characterized neurophysiological alterations in mice and humans carrying AD mutations in the APP and/or PSEN1 genes, focusing on early pre-symptomatic changes. Longitudinal local field potential recordings were performed in APP/PS1 mice and cross-sectional magnetoencephalography recordings in human APP and/or PSEN1 mutation carriers. All recordings were acquired in the left frontal cortex, parietal cortex, and hippocampus. Spectral power and functional connectivity were analyzed and compared with wildtype control mice and healthy age-matched human subjects.

RESULTS

APP/PS1 mice showed increased absolute power, especially at higher frequencies (beta and gamma) and predominantly between 3 and 6 moa. Relative power showed an overall shift from lower to higher frequencies over almost the entire recording period and across all three brain regions. Human mutation carriers, on the other hand, did not show changes in power except for an increase in relative theta power in the hippocampus. Mouse parietal cortex and hippocampal power spectra showed a characteristic peak at around 8 Hz which was not significantly altered in transgenic mice. Human power spectra showed a characteristic peak at around 9 Hz, the frequency of which was significantly reduced in mutation carriers. Significant alterations in functional connectivity were detected in theta, alpha, beta, and gamma frequency bands, but the exact frequency range and direction of change differed for APP/PS1 mice and human mutation carriers.

CONCLUSIONS

Both mice and humans carrying APP and/or PSEN1 mutations show abnormal neurophysiological activity, but several measures do not translate one-to-one between species. Alterations in absolute and relative power in mice should be interpreted with care and may be due to overexpression of amyloid in combination with the absence of tau pathology and cholinergic degeneration. Future studies should explore whether changes in brain activity in other AD mouse models, for instance, those also including tau pathology, provide better translation to the human AD continuum.

Top-down reconfiguration of SMA cortical connectivity during action preparation

The Bereitschaftspotential (BP), a scalp potential recorded in humans during action preparation, is characterized by a slow amplitude increase over fronto-central regions as action execution approaches. We recorded TMS evoked-potentials (TEP) stimulating the supplementary motor area (SMA) at different time-points during a Go/No-Go task to assess whether and how cortical excitability and connectivity of this region change as the BP increases. When approaching BP peak, left SMA reactivity resulted greater. Concurrently, its effective connectivity increased with the left occipital areas, while it decreased with the right inferior frontal gyrus, indicating a fast reconfiguration of cortical networks during the preparation of the forthcoming action. Functional connectivity patterns supported these findings, suggesting a critical role of frequency-specific inter-areal interactions in implementing top-down mechanisms in the sensorimotor system prior to action. These findings reveal that BP time-course reflects quantitative and qualitative changes in SMA communication patterns that shape mechanisms involved in motor readiness.

Identifying good practices for detecting inter-regional linear functional connectivity from EEG

Aggregating voxel-level statistical dependencies between multivariate time series is an important intermediate step when characterising functional connectivity (FC) between larger brain regions. However, there are numerous ways in which voxel-level data can be aggregated into inter-regional FC, and the advantages of each of these approaches are currently unclear. In this study we generate ground-truth data and compare the performances of various pipelines that estimate directed and undirected linear phase-to-phase FC between regions. We test the ability of several existing and novel FC analysis pipelines to identify the true regions within which connectivity was simulated. We test various inverse modelling algorithms, strategies to aggregate time series within regions, and connectivity metrics. Furthermore, we investigate the influence of the number of interactions, the signal-to-noise ratio, the noise mix, the interaction time delay, and the number of active sources per region on the ability of detecting phase-to-phase FC. Throughout all simulated scenarios, lowest performance is obtained with pipelines involving the absolute value of coherency. Further, the combination of dynamic imaging of coherent sources (DICS) beamforming with directed FC metrics that aggregate information across multiple frequencies leads to unsatisfactory results. Pipelines that show promising results with our simulated pseudo-EEG data involve the following steps: (1) Source projection using the linearly-constrained minimum variance (LCMV) beamformer. (2) Principal component analysis (PCA) using the same fixed number of components within every region. (3) Calculation of the multivariate interaction measure (MIM) for every region pair to assess undirected phase-to-phase FC, or calculation of time-reversed Granger Causality (TRGC) to assess directed phase-to-phase FC. We formulate recommendations based on these results that may increase the validity of future experimental connectivity studies. We further introduce the free ROIconnect plugin for the EEGLAB toolbox that includes the recommended methods and pipelines that are presented here. We show an exemplary application of the best performing pipeline to the analysis of EEG data recorded during motor imagery.

The longitudinal relation between executive functioning and multilayer network topology in glioma patients

Many patients with glioma, primary brain tumors, suffer from poorly understood executive functioning deficits before and/or after tumor resection. We aimed to test whether frontoparietal network centrality of multilayer networks, allowing for integration across multiple frequencies, relates to and predicts executive functioning in glioma. Patients with glioma (n = 37) underwent resting-state magnetoencephalography and neuropsychological tests assessing word fluency, inhibition, and set shifting before (T1) and one year after tumor resection (T2). We constructed binary multilayer networks comprising six layers, with each layer representing frequency-specific functional connectivity between source-localized time series of 78 cortical regions. Average frontoparietal network multilayer eigenvector centrality, a measure for network integration, was calculated at both time points. Regression analyses were used to investigate associations with executive functioning. At T1, lower multilayer integration (p = 0.017) and epilepsy (p = 0.006) associated with poorer set shifting (adj. R2 = 0.269). Decreasing multilayer integration (p = 0.022) and not undergoing chemotherapy at T2 (p = 0.004) related to deteriorating set shifting over time (adj. R2 = 0.283). No significant associations were found for word fluency or inhibition, nor did T1 multilayer integration predict changes in executive functioning. As expected, our results establish multilayer integration of the frontoparietal network as a cross-sectional and longitudinal correlate of executive functioning in glioma patients. However, multilayer integration did not predict postoperative changes in executive functioning, which together with the fact that this correlate is also found in health and other diseases, limits its specific clinical relevance in glioma.

Non-reversibility outperforms functional connectivity in characterisation of brain states in MEG data

Characterising brain states during tasks is common practice for many neuroscientific experiments using electrophysiological modalities such as electroencephalography (EEG) and magnetoencephalography (MEG). Brain states are often described in terms of oscillatory power and correlated brain activity, i.e. functional connectivity. It is, however, not unusual to observe weak task induced functional connectivity alterations in the presence of strong task induced power modulations using classical time-frequency representation of the data. Here, we propose that non-reversibility, or the temporal asymmetry in functional interactions, may be more sensitive to characterise task induced brain states than functional connectivity. As a second step, we explore causal mechanisms of non-reversibility in MEG data using whole brain computational models. We include working memory, motor, language tasks and resting-state data from participants of the Human Connectome Project (HCP). Non-reversibility is derived from the lagged amplitude envelope correlation (LAEC), and is based on asymmetry of the forward and reversed cross-correlations of the amplitude envelopes. Using random forests, we find that non-reversibility outperforms functional connectivity in the identification of task induced brain states. Non-reversibility shows especially better sensitivity to capture bottom-up gamma induced brain states across all tasks, but also alpha band associated brain states. Using whole brain computational models we find that asymmetry in the effective connectivity and axonal conduction delays play a major role in shaping non-reversibility across the brain. Our work paves the way for better sensitivity in characterising brain states during both bottom-up as well as top-down modulation in future neuroscientific experiments.

Uncorrelated bilateral cortical input becomes timed across hippocampal subfields for long waves whereas gamma waves are largely ipsilateral

The role of interhemispheric connections along successive segments of cortico-hippocampal circuits is poorly understood. We aimed to obtain a global picture of spontaneous transfer of activity during non-theta states across several nodes of the bilateral circuit in anesthetized rats. Spatial discrimination techniques applied to bilateral laminar field potentials (FP) across the CA1/Dentate Gyrus provided simultaneous left and right readouts in five FP generators that reflect activity in specific hippocampal afferents and associative pathways. We used a battery of correlation and coherence analyses to extract complementary aspects at different time scales and frequency bands. FP generators exhibited varying bilateral correlation that was high in CA1 and low in the Dentate Gyrus. The submillisecond delays indicate coordination but not support for synaptic dependence of one side on another. The time and frequency characteristics of bilateral coupling were specific to each generator. The Schaffer generator was strongly bilaterally coherent for both sharp waves and gamma waves, although the latter maintained poor amplitude co-variation. The lacunosum-moleculare generator was composed of up to three spatially overlapping activities, and globally maintained high bilateral coherence for long but not short (gamma) waves. These two CA1 generators showed no ipsilateral relationship in any frequency band. In the Dentate Gyrus, strong bilateral coherence was observed only for input from the medial entorhinal areas, while those from the lateral entorhinal areas were largely asymmetric, for both alpha and gamma waves. Granger causality testing showed strong bidirectional relationships between all homonymous bilateral generators except the lateral entorhinal input and a local generator in the Dentate Gyrus. It also revealed few significant relationships between ipsilateral generators, most notably the anticipation of lateral entorhinal cortex toward all others. Thus, with the notable exception of the lateral entorhinal areas, there is a marked interhemispheric coherence primarily for slow envelopes of activity, but not for pulse-like gamma waves, except in the Schafer segment. The results are consistent with essentially different streams of activity entering from and returning to the cortex on each side, with slow waves reflecting times of increased activity exchange between hemispheres and fast waves generally reflecting ipsilateral processing.

Identifying oscillatory brain networks with hidden Gaussian graphical spectral models of MEEG

Identifying the functional networks underpinning indirectly observed processes poses an inverse problem for neurosciences or other fields. A solution of such inverse problems estimates as a first step the activity emerging within functional networks from EEG or MEG data. These EEG or MEG estimates are a direct reflection of functional brain network activity with a temporal resolution that no other in vivo neuroimage may provide. A second step estimating functional connectivity from such activity pseudodata unveil the oscillatory brain networks that strongly correlate with all cognition and behavior. Simulations of such MEG or EEG inverse problem also reveal estimation errors of the functional connectivity determined by any of the state-of-the-art inverse solutions. We disclose a significant cause of estimation errors originating from misspecification of the functional network model incorporated into either inverse solution steps. We introduce the Bayesian identification of a Hidden Gaussian Graphical Spectral (HIGGS) model specifying such oscillatory brain networks model. In human EEG alpha rhythm simulations, the estimation errors measured as ROC performance do not surpass 2% in our HIGGS inverse solution and reach 20% in state-of-the-art methods. Macaque simultaneous EEG/ECoG recordings provide experimental confirmation for our results with 1/3 times larger congruence according to Riemannian distances than state-of-the-art methods.

Network Hyperexcitability in Early Alzheimer's Disease: Is Functional Connectivity a Potential Biomarker?

Network hyperexcitability (NH) is an important feature of the pathophysiology of Alzheimer's disease. Functional connectivity (FC) of brain networks has been proposed as a potential biomarker for NH. Here we use a whole brain computational model and resting-state MEG recordings to investigate the relation between hyperexcitability and FC. Oscillatory brain activity was simulated with a Stuart Landau model on a network of 78 interconnected brain regions. FC was quantified with amplitude envelope correlation (AEC) and phase coherence (PC). MEG was recorded in 18 subjects with subjective cognitive decline (SCD) and 18 subjects with mild cognitive impairment (MCI). Functional connectivity was determined with the corrected AECc and phase lag index (PLI), in the 4-8 Hz and the 8-13 Hz bands. The excitation/inhibition balance in the model had a strong effect on both AEC and PC. This effect was different for AEC and PC, and was influenced by structural coupling strength and frequency band. Empirical FC matrices of SCD and MCI showed a good correlation with model FC for AEC, but less so for PC. For AEC the fit was best in the hyperexcitable range. We conclude that FC is sensitive to changes in E/I balance. The AEC was more sensitive than the PLI, and results were better for the thetaband than the alpha band. This conclusion was supported by fitting the model to empirical data. Our study justifies the use of functional connectivity measures as surrogate markers for E/I balance.

Magnetoencephalography to measure the effect of contact point-specific deep brain stimulation in Parkinson's disease: A proof of concept study

BACKGROUND

Deep brain stimulation (DBS) of the subthalamic nucleus (STN) is an effective treatment for disabling fluctuations in motor symptoms in Parkinson's disease (PD) patients. However, iterative exploration of all individual contact points (four in each STN) by the clinician for optimal clinical effects may take months.

OBJECTIVE

In this proof of concept study we explored whether magnetoencephalography (MEG) has the potential to noninvasively measure the effects of changing the active contact point of STN-DBS on spectral power and functional connectivity in PD patients, with the ultimate aim to aid in the process of selecting the optimal contact point, and perhaps reduce the time to achieve optimal stimulation settings.

METHODS

The study included 30 PD patients who had undergone bilateral DBS of the STN. MEG was recorded during stimulation of each of the eight contact points separately (four on each side). Each stimulation position was projected on a vector running through the longitudinal axis of the STN, leading to one scalar value indicating a more dorsolateral or ventromedial contact point position. Using linear mixed models, the stimulation positions were correlated with band-specific absolute spectral power and functional connectivity of i) the motor cortex ipsilateral tot the stimulated side, ii) the whole brain.

RESULTS

At group level, more dorsolateral stimulation was associated with lower low-beta absolute band power in the ipsilateral motor cortex (p = .019). More ventromedial stimulation was associated with higher whole-brain absolute delta (p = .001) and theta (p = .005) power, as well as higher whole-brain theta band functional connectivity (p = .040). At the level of the individual patient, switching the active contact point caused significant changes in spectral power, but the results were highly variable.

CONCLUSIONS

We demonstrate for the first time that stimulation of the dorsolateral (motor) STN in PD patients is associated with lower low-beta power values in the motor cortex. Furthermore, our group-level data show that the location of the active contact point correlates with whole-brain brain activity and connectivity. As results in individual patients were quite variable, it remains unclear if MEG is useful in the selection of the optimal DBS contact point.

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.

Resting-state oscillations reveal disturbed excitation-inhibition ratio in Alzheimer's disease patients

An early disruption of neuronal excitation-inhibition (E-I) balance in preclinical animal models of Alzheimer's disease (AD) has been frequently reported, but is difficult to measure directly and non-invasively in humans. Here, we examined known and novel neurophysiological measures sensitive to E-I in patients across the AD continuum. Resting-state magnetoencephalography (MEG) data of 86 amyloid-biomarker-confirmed subjects across the AD continuum (17 patients diagnosed with subjective cognitive decline, 18 with mild cognitive impairment (MCI) and 51 with dementia due to probable AD (AD dementia)), 46 healthy elderly and 20 young control subjects were reconstructed to source-space. E-I balance was investigated by detrended fluctuation analysis (DFA), a functional E/I (fE/I) algorithm, and the aperiodic exponent of the power spectrum. We found a disrupted E-I ratio in AD dementia patients specifically, by a lower DFA, and a shift towards higher excitation, by a higher fE/I and a lower aperiodic exponent. Healthy subjects showed lower fE/I ratios (< 1.0) than reported in previous literature, not explained by age or choice of an arbitrary threshold parameter, which warrants caution in interpretation of fE/I results. Correlation analyses showed that a lower DFA (E-I imbalance) and a lower aperiodic exponent (more excitation) was associated with a worse cognitive score in AD dementia patients. In contrast, a higher DFA in the hippocampi of MCI patients was associated with a worse cognitive score. This MEG-study showed E-I imbalance, likely due to increased excitation, in AD dementia, but not in early stage AD patients. To accurately determine the direction of shift in E-I balance, validations of the currently used markers and additional in vivo markers of E-I are required.

Altered directional functional connectivity underlies post-stroke cognitive recovery

Cortical ischaemic strokes result in cognitive deficits depending on the area of the affected brain. However, we have demonstrated that difficulties with attention and processing speed can occur even with small subcortical infarcts. Symptoms appear independent of lesion location, suggesting they arise from generalized disruption of cognitive networks. Longitudinal studies evaluating directional measures of functional connectivity in this population are lacking. We evaluated six patients with minor stroke exhibiting cognitive impairment 6-8 weeks post-infarct and four age-similar controls. Resting-state magnetoencephalography data were collected. Clinical and imaging evaluations of both groups were repeated 6- and 12 months later. Network Localized Granger Causality was used to determine differences in directional connectivity between groups and across visits, which were correlated with clinical performance. Directional connectivity patterns remained stable across visits for controls. After the stroke, inter-hemispheric connectivity between the frontoparietal cortex and the non-frontoparietal cortex significantly increased between visits 1 and 2, corresponding to uniform improvement in reaction times and cognitive scores. Initially, the majority of functional links originated from non-frontal areas contralateral to the lesion, connecting to ipsilesional brain regions. By visit 2, inter-hemispheric connections, directed from the ipsilesional to the contralesional cortex significantly increased. At visit 3, patients demonstrating continued favourable cognitive recovery showed less reliance on these inter-hemispheric connections. These changes were not observed in those without continued improvement. Our findings provide supporting evidence that the neural basis of early post-stroke cognitive dysfunction occurs at the network level, and continued recovery correlates with the evolution of inter-hemispheric connectivity.

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.

Effects of gamma frequency binaural beats on attention and anxiety

Binaural beats (BB) are an auditory phenomenon produced from a combination of two sine waves with slightly different frequencies presented to each ear. Previous research has implicated the role of BBs through brainwave entrainment in potentially giving rise to benefits ranging from enhanced memory and attention to reduced anxiety and stress. Here, we investigated the effect of gamma (40-Hz) BBs on attention using the attention network test (ANT), a previously unused task that assesses three subtypes of attention: Alerting, Orienting, and Executive Control. Fifty-eight healthy adults performed the ANT remotely under the exposure of 340-Hz BBs and a 380-Hz control tone. All completed a rating scale for levels of anxiety before and after each exposure. Performance on the ANT task (reaction time and error rates) between BB and control groups was evaluated using Wilcoxon signed-rank tests. We found no significant differences in Reaction Time (RT), Error Rate (ER), or the efficacy of the Attention Networks (AN) between the experimental and control conditions (p > 0.05). We found no effect of BB on self-rated measures of anxiety. Our findings do not provide evidence for improvement in attention with gamma BB.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12144-023-04681-3.

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.

Brain-wide network analysis of resting-state neuromagnetic data

The present study performed a brain-wide network analysis of resting-state magnetoencephalograms recorded from 53 healthy participants to visualize elaborate brain maps of phase- and amplitude-derived graph-theory metrics at different frequencies. To achieve this, we conducted a vertex-wise computation of threshold-independent graph metrics by combining proportional thresholding and a conjunction analysis and applied them to a correlation analysis of age and brain networks. Source power showed a frequency-dependent cortical distribution. Threshold-independent graph metrics derived from phase- and amplitude-based connectivity showed similar or different distributions depending on frequency. Vertex-wise age-brain correlation maps revealed that source power at the beta band and the amplitude-based degree at the alpha band changed with age in local regions. The present results indicate that a brain-wide analysis of neuromagnetic data has the potential to reveal neurophysiological network features in the human brain in a resting state.

Non-invasive measurements of ictal and interictal epileptiform activity using optically pumped magnetometers

Magneto- and electroencephalography (MEG/EEG) are important techniques for the diagnosis and pre-surgical evaluation of epilepsy. Yet, in current cryogen-based MEG systems the sensors are offset from the scalp, which limits the signal-to-noise ratio (SNR) and thereby the sensitivity to activity from deep structures such as the hippocampus. This effect is amplified in children, for whom adult-sized fixed-helmet systems are typically too big. Moreover, ictal recordings with fixed-helmet systems are problematic because of limited movement tolerance and/or logistical considerations. Optically Pumped Magnetometers (OPMs) can be placed directly on the scalp, thereby improving SNR and enabling recordings during seizures. We aimed to demonstrate the performance of OPMs in a clinical population. Seven patients with challenging cases of epilepsy underwent MEG recordings using a 12-channel OPM-system and a 306-channel cryogen-based whole-head system: three adults with known deep or weak (low SNR) sources of interictal epileptiform discharges (IEDs), along with three children with focal epilepsy and one adult with frequent seizures. The consistency of the recorded IEDs across the two systems was assessed. In one patient the OPMs detected IEDs that were not found with the SQUID-system, and in two patients no IEDs were found with either system. For the other patients the OPM data were remarkably consistent with the data from the cryogenic system, noting that these were recorded in different sessions, with comparable SNRs and IED-yields overall. Importantly, the wearability of OPMs enabled the recording of seizure activity in a patient with hyperkinetic movements during the seizure. The observed ictal onset and semiology were in agreement with previous video- and stereo-EEG recordings. The relatively affordable technology, in combination with reduced running and maintenance costs, means that OPM-based MEG could be used more widely than current MEG systems, and may become an affordable alternative to scalp EEG, with the potential benefits of increased spatial accuracy, reduced sensitivity to volume conduction/field spread, and increased sensitivity to deep sources. Wearable MEG thus provides an unprecedented opportunity for epilepsy, and given its patient-friendliness, we envisage that it will not only be used for presurgical evaluation of epilepsy patients, but also for diagnosis after a first seizure.

Alpha and gamma EEG coherence during on-task and mind wandering states in schizophrenia

OBJECTIVE

Electroencephalographic (EEG) coherence is one of the most relevant physiological measures used to detect abnormalities in patients with schizophrenia. The present study applies a task-related EEG coherence approach to understand cognitive processing in patients with schizophrenia and healthy controls.

METHODS

EEG coherence for alpha and gamma frequency bands was analyzed in a group of patients with schizophrenia and a group of healthy controls during the performance of an ecological task of sustained attention. We compared EEG coherence when participants presented externally directed cognitive states (On-Task) and when they presented cognitive distraction episodes (Mind-Wandering).

RESULTS

Results reflect cortical differences between groups (higher coherence for schizophrenia in the frontocentral and fronto-temporal regions, and higher coherence for healthy-controls in the postero-central regions), especially in the On-Task condition for the alpha band, compared to Mind-Wandering episodes. Few individual differences in gamma coherence were found.

CONCLUSIONS

The current study provides evidence of neurophysiological differences underlying different cognitive states in schizophrenia and healthy controls.

SIGNIFICANCE

Differences between groups may reflect inhibitory processes necessary for the successful processing of information, especially in the alpha band, given its role in cortical inhibition processes. Patients may activate compensatory inhibitory mechanisms when performing the task, reflected in increased coherence in fronto-temporal regions.

Association between cognitive impairments and aberrant dynamism of overlapping brain sub-networks in unmedicated major depressive disorder: A resting-state MEG study

OBJECTIVE

Little is known about the pathogenesis underlying cognitive impairment in major depressive disorder (MDD). We aimed to explore the mechanisms of cognitive impairments among patients with MDD by investigating the dynamics of overlapping brain sub-networks.

METHODS

Forty unmedicated patients with MDD and 28 healthy controls (HC) were enrolled in this study. Cognitive function was measured using the Chinese versions of MATRICS Consensus Cognitive Battery (MCCB). All participants were scanned using a whole-head resting-state magnetoencephalography (MEG) machine. The dynamism of neural sub-networks was analyzed based on the detection of overlapping communities in five frequency bands of oscillatory brain signals.

RESULTS

MDD demonstrated poorer cognitive performance in six domains compared to HC. The difference in community detection (functional integration mode) in MDD was frequency-dependent. MDD showed significantly decreased community dynamics in all frequency bands compared to HC. Specifically, differences in the visual network (VN) and default mode network (DMN) were detected in all frequency bands, differences in the cognitive control network (CCN) were detected in the alpha2 and beta frequency bands, and differences in the bilateral limbic network (BLN) were only detected in the beta frequency band. Moreover, community dynamics in the alpha2 frequency band were positively correlated with verbal learning and reasoning problem solving abilities in MDD.

CONCLUSIONS

Our study found that decreasing in the dynamics of overlapping sub-networks may differ by frequency bands. The aberrant dynamics of overlapping neural sub-networks revealed by frequency-specific MEG signals may provide new information on the mechanism of cognitive impairments that result from MDD.

Multimodal multilayer network centrality relates to executive functioning

Executive functioning (EF) is a higher order cognitive process that is thought to depend on a network organization facilitating integration across subnetworks, in the context of which the central role of the fronto-parietal network (FPN) has been described across imaging and neurophysiological modalities. However, the potentially complementary unimodal information on the relevance of the FPN for EF has not yet been integrated. We employ a multilayer framework to allow for integration of different modalities into one 'network of networks.' We used diffusion MRI, resting-state functional MRI, MEG, and neuropsychological data obtained from 33 healthy adults to construct modality-specific single-layer networks as well as a single multilayer network per participant. We computed single-layer and multilayer eigenvector centrality of the FPN as a measure of integration in this network and examined their associations with EF. We found that higher multilayer FPN centrality, but not single-layer FPN centrality, was related to better EF. We did not find a statistically significant change in explained variance in EF when using the multilayer approach as compared to the single-layer measures. Overall, our results show the importance of FPN integration for EF and underline the promise of the multilayer framework toward better understanding cognitive functioning.

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.

Network alterations in eating epilepsy during resting state and during eating using Magnetoencephalography

OBJECTIVE

Eating epilepsy presents various imaging and electrophysiological features along with various seizure triggers. As such, network changes in eating epilepsy have not been comprehensively explored. This study was conducted to illustrate resting state network changes in eating epilepsy and to study the changes in network configurations during eating.

METHODS

Magnetoencephalography recordings of nineteen patients with drug-resistant eating epilepsy were compared with healthy controls during resting state. A subgroup of nine patients and 12 controls had MEG recordings during eating. Network changes were analyzed using phase lag index across 5 frequency bands [delta, theta, alpha, beta, and gamma] using clustering coefficient (CC), betweenness centrality (BC), path length (PL), modularity (Q), and small worldness (SW).

RESULTS

During the resting state, PL was decreased in patients with epilepsy in the delta, theta, and gamma band. Q was lower in patients with epilepsy in the beta and gamma bands. During eating, in patients with epilepsy, PL and SW were increased in all frequency bands, and Q was decreased in the beta band and increased in the rest of the frequency bands. Patients with mixed types of seizures showed higher PL in all bands except alpha, higher Q in all bands, and higher SW in the alpha and beta bands. Node-wise changes in CC and BC implicated changes in DMN and 'eating' networks.

CONCLUSION

Reflex Eating epilepsy presents with a hyperconnected network that exacerbates during eating. The cause of seizure onset and loss of consciousness in eating epilepsy might be due to aberrant network interaction between the regions of the brain involved with eating, such as the sensorimotor cortex, lateral parietal cortex, and insula with the limbic cortex and default mode network across multiple frequency bands.

A multiscale brain network model links Alzheimer’s disease-mediated neuronal hyperactivity to large-scale oscillatory slowing

Background

Neuronal hyperexcitability and inhibitory interneuron dysfunction are frequently observed in preclinical animal models of Alzheimer’s disease (AD). This study investigates whether these microscale abnormalities explain characteristic large-scale magnetoencephalography (MEG) activity in human early-stage AD patients.

Methods

To simulate spontaneous electrophysiological activity, we used a whole-brain computational network model comprised of 78 neural masses coupled according to human structural brain topology. We modified relevant model parameters to simulate six literature-based cellular scenarios of AD and compare them to one healthy and six contrast (non-AD-like) scenarios. The parameters include excitability, postsynaptic potentials, and coupling strength of excitatory and inhibitory neuronal populations. Whole-brain spike density and spectral power analyses of the simulated data reveal mechanisms of neuronal hyperactivity that lead to oscillatory changes similar to those observed in MEG data of 18 human prodromal AD patients compared to 18 age-matched subjects with subjective cognitive decline.

Results

All but one of the AD-like scenarios showed higher spike density levels, and all but one of these scenarios had a lower peak frequency, higher spectral power in slower (theta, 4–8Hz) frequencies, and greater total power. Non-AD-like scenarios showed opposite patterns mainly, including reduced spike density and faster oscillatory activity. Human AD patients showed oscillatory slowing (i.e., higher relative power in the theta band mainly), a trend for lower peak frequency and higher total power compared to controls. Combining model and human data, the findings indicate that neuronal hyperactivity can lead to oscillatory slowing, likely due to hyperexcitation (by hyperexcitability of pyramidal neurons or greater long-range excitatory coupling) and/or disinhibition (by reduced excitability of inhibitory interneurons or weaker local inhibitory coupling strength) in early AD.

Conclusions

Using a computational brain network model, we link findings from different scales and models and support the hypothesis of early-stage neuronal hyperactivity underlying E/I imbalance and whole-brain network dysfunction in prodromal AD.

Supplementary Information

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

Modulation of cortical beta oscillations influences motor vigor: A rhythmic TMS-EEG study

Previous electro- or magnetoencephalography (Electro/Magneto EncephaloGraphic; E/MEG) studies using a correlative approach have shown that β (13-30 Hz) oscillations emerging in the primary motor cortex (M1) are implicated in regulating motor response vigor and associated with an anti-kinetic role, that is, slowness of movement. However, the functional role of M1 β oscillations in regulation of motor responses remains unclear. To address this gap, we combined EEG with rhythmic TMS (rhTMS) delivered to M1 at the β (20 Hz) frequency shortly before subjects performed an isometric ramp-and-hold finger force production task at three force levels. rhTMS is a novel approach that can modulate rhythmic patterns of neural activity. β-rhTMS over M1 induced a modulation of neural oscillations to β frequency in the sensorimotor area and reduced peak force rate during the ramp-up period relative to sham and catch trials. Interestingly, this rhTMS effect occurred only in the large force production condition. To distinguish whether the effects of rhTMS on EEG and behavior stemmed from phase-resetting by each magnetic pulse or neural entrainment by the periodicity of rhTMS, we performed a control experiment using arrhythmic TMS (arTMS). arTMS did not induce changes in EEG oscillations nor peak force rate during the rump-up period. Our results provide novel evidence that β neural oscillations emerging the sensorimotor area influence the regulation of motor response vigor. Furthermore, our findings further demonstrate that rhTMS is a promising tool for tuning neural oscillations to the target frequency.

Modulations of local synchrony over time lead to resting-state functional connectivity in a parsimonious large-scale brain model

Biophysical models of large-scale brain activity are a fundamental tool for understanding the mechanisms underlying the patterns observed with neuroimaging. These models combine a macroscopic description of the within- and between-ensemble dynamics of neurons within a single architecture. A challenge for these models is accounting for modulations of within-ensemble synchrony over time. Such modulations in local synchrony are fundamental for modeling behavioral tasks and resting-state activity. Another challenge comes from the difficulty in parametrizing large scale brain models which hinders researching principles related with between-ensembles differences. Here we derive a parsimonious large scale brain model that can describe fluctuations of local synchrony. Crucially, we do not reduce within-ensemble dynamics to macroscopic variables first, instead we consider within and between-ensemble interactions similarly while preserving their physiological differences. The dynamics of within-ensemble synchrony can be tuned with a parameter which manipulates local connectivity strength. We simulated resting-state static and time-resolved functional connectivity of alpha band envelopes in models with identical and dissimilar local connectivities. We show that functional connectivity emerges when there are high fluctuations of local and global synchrony simultaneously (i.e. metastable dynamics). We also show that for most ensembles, leaning towards local asynchrony or synchrony correlates with the functional connectivity with other ensembles, with the exception of some regions belonging to the default-mode network.

A systematic data-driven approach to analyze sensor-level EEG connectivity: Identifying robust phase-synchronized network components using surface Laplacian with spectral-spatial PCA

Although conventional averaging across predefined frequency bands reduces the complexity of EEG functional connectivity (FC), it obscures the identification of resting-state brain networks (RSN) and impedes accurate estimation of FC reliability. Extending prior work, we combined scalp current source density (CSD; spherical spline surface Laplacian) and spectral-spatial PCA to identify FC components. Phase-based FC was estimated via debiased-weighted phase-locking index from CSD-transformed resting EEGs (71 sensors, 8 min, eyes open/closed, 35 healthy adults, 1-week retest). Spectral PCA extracted six robust alpha and theta components (86.6% variance). Subsequent spatial PCA for each spectral component revealed seven robust regionally focused (posterior, central, and frontal) and long-range (posterior-anterior) alpha components (peaks at 8, 10, and 13 Hz) and a midfrontal theta (6 Hz) component, accounting for 37.0% of FC variance. These spatial FC components were consistent with well-known networks (e.g., default mode, visual, and sensorimotor), and four were sensitive to eyes open/closed conditions. Most FC components had good-to-excellent internal consistency (odd/even epochs, eyes open/closed) and test-retest reliability (ICCs ≥ .8). Moreover, the FC component structure was generally present in subsamples (session × odd/even epoch, or smaller subgroups [n = 7-10]), as indicated by high similarity of component loadings across PCA solutions. Apart from systematically reducing FC dimensionality, our approach avoids arbitrary thresholds and allows quantification of meaningful and reliable network components that may prove to be of high relevance for basic and clinical research applications.

Regional healthy brain activity, glioma occurrence and symptomatology

It is unclear why exactly gliomas show preferential occurrence in certain brain areas. Increased spiking activity around gliomas leads to faster tumour growth in animal models, while higher non-invasively measured brain activity is related to shorter survival in patients. However, it is unknown how regional intrinsic brain activity, as measured in healthy controls, relates to glioma occurrence. We first investigated whether gliomas occur more frequently in regions with intrinsically higher brain activity. Second, we explored whether intrinsic cortical activity at individual patients’ tumour locations relates to tumour and patient characteristics.

Across three cross-sectional cohorts, 413 patients were included. Individual tumour masks were created. Intrinsic regional brain activity was assessed through resting-state magnetoencephalography acquired in healthy controls and source-localized to 210 cortical brain regions. Brain activity was operationalized as: (i) broadband power; and (ii) offset of the aperiodic component of the power spectrum, which both reflect neuronal spiking of the underlying neuronal population. We additionally assessed (iii) the slope of the aperiodic component of the power spectrum, which is thought to reflect the neuronal excitation/inhibition ratio. First, correlation coefficients were calculated between group-level regional glioma occurrence, as obtained by concatenating tumour masks across patients, and group-averaged regional intrinsic brain activity. Second, intrinsic brain activity at specific tumour locations was calculated by overlaying patients’ individual tumour masks with regional intrinsic brain activity of the controls and was associated with tumour and patient characteristics.

As proposed, glioma preferentially occurred in brain regions characterized by higher intrinsic brain activity in controls as reflected by higher offset. Second, intrinsic brain activity at patients’ individual tumour locations differed according to glioma subtype and performance status: the most malignant isocitrate dehydrogenase-wild-type glioblastoma patients had the lowest excitation/inhibition ratio at their individual tumour locations as compared to isocitrate dehydrogenase-mutant, 1p/19q-codeleted glioma patients, while a lower excitation/inhibition ratio related to poorer Karnofsky Performance Status, particularly in codeleted glioma patients.

In conclusion, gliomas more frequently occur in cortical brain regions with intrinsically higher activity levels, suggesting that more active regions are more vulnerable to glioma development. Moreover, indices of healthy, intrinsic excitation/inhibition ratio at patients’ individual tumour locations may capture both tumour biology and patients’ performance status. These findings contribute to our understanding of the complex and bidirectional relationship between normal brain functioning and glioma growth, which is at the core of the relatively new field of ‘cancer neuroscience’.

Abnormal information interaction in multilayer directed network based on cross-frequency integration of mild cognitive impairment and Alzheimer’s disease

Mild cognitive impairment (MCI) and Alzheimer’s disease (AD) have been reported to result in abnormal cross-frequency integration. However, previous studies have failed to consider specific abnormalities in receiving and outputting information among frequency bands during integration. Here, we investigated heterogeneity in receiving and outputting information during cross-frequency integration in patients. The results showed that during cross-frequency integration, information interaction first increased and then decreased, manifesting in the heterogeneous distribution of inter-frequency nodes for receiving information. A possible explanation was that due to damage to some inter-frequency hub nodes, intra-frequency nodes gradually became new inter-frequency nodes, whereas original inter-frequency nodes gradually became new inter-frequency hub nodes. Notably, damage to the brain regions that receive information between layers was often accompanied by a strengthened ability to output information and the emergence of hub nodes for outputting information. Moreover, an important compensatory mechanism assisted in the reception of information in the cingulo-opercular and auditory networks and in the outputting of information in the visual network. This study revealed specific abnormalities in information interaction and compensatory mechanism during cross-frequency integration, providing important evidence for understanding cross-frequency integration in patients with MCI and AD.

The expanding horizons of network neuroscience: From description to prediction and control

The field of network neuroscience has emerged as a natural framework for the study of the brain and has been increasingly applied across divergent problems in neuroscience. From a disciplinary perspective, network neuroscience originally emerged as a formal integration of graph theory (from mathematics) and neuroscience (from biology). This early integration afforded marked utility in describing the interconnected nature of neural units, both structurally and functionally, and underscored the relevance of that interconnection for cognition and behavior. But since its inception, the field has not remained static in its methodological composition. Instead, it has grown to use increasingly advanced graph-theoretic tools and to bring in several other disciplinary perspectives-including machine learning and systems engineering-that have proven complementary. In doing so, the problem space amenable to the discipline has expanded markedly. In this review, we discuss three distinct flavors of investigation in state-of-the-art network neuroscience: (i) descriptive network neuroscience, (ii) predictive network neuroscience, and (iii) a perturbative network neuroscience that draws on recent advances in network control theory. In considering each area, we provide a brief summary of the approaches, discuss the nature of the insights obtained, and highlight future directions.

Space-time resolved inference-based neurophysiological process imaging: application to resting-state alpha rhythm

Neural processes are complex and difficult to image. This paper presents a new space-time resolved brain imaging framework, called Neurophysiological Process Imaging (NPI), that identifies neurophysiological processes within cerebral cortex at the macroscopic scale. By fitting uncoupled neural mass models to each electromagnetic source time-series using a novel nonlinear inference method, population averaged membrane potentials and synaptic connection strengths are efficiently and accurately inferred and imaged across the whole cerebral cortex at a resolution afforded by source imaging. The efficiency of the framework enables return of the augmented source imaging results overnight using high performance computing. This suggests it can be used as a practical and novel imaging tool. To demonstrate the framework, it has been applied to resting-state magnetoencephalographic source estimates. The results suggest that endogenous inputs to cingulate, occipital, and inferior frontal cortex are essential modulators of resting-state alpha power. Moreover, endogenous input and inhibitory and excitatory neural populations play varied roles in mediating alpha power in different resting-state sub-networks. The framework can be applied to arbitrary neural mass models and has broad applicability to image neural processes in different brain states.

EEG functional connectivity analysis in the source space

There is a growing interest in using electroencephalography (EEG) and source modeling to investigate functional interactions among cortical processes, particularly when dealing with pediatric populations. This paper introduces two pipelines that have been recently used to conduct EEG FC analysis in the cortical source space. The analytic streams of these pipelines can be summarized into the following steps: 1) cortical source reconstruction of high-density EEG data using realistic magnetic resonance imaging (MRI) models created with age-appropriate MRI templates; 2) segmentation of reconstructed source activities into brain regions of interest; and 3) estimation of FC in age-related frequency bands using robust EEG FC measures, such as weighted phase lag index and orthogonalized power envelope correlation. In this paper we demonstrate the two pipelines with resting-state EEG data collected from children at 12 and 36 months of age. We also discuss the advantages and limitations of the methods/techniques integrated into the pipelines. Given there is a need in the research community for open-access analytic toolkits that can be used for pediatric EEG data, programs and codes used for the current analysis are made available to the public.

Exploring sex differences in alpha brain activity as a potential neuromarker associated with neuropathic pain

ABSTRACT

Alpha oscillatory activity (8-13 Hz) is the dominant rhythm in the awake brain and is known to play an important role in pain states. Previous studies have identified alpha band slowing and increased power in the dynamic pain connectome (DPC) of people with chronic neuropathic pain. However, a link between alpha-band abnormalities and sex differences in brain organization in healthy individuals and those with chronic pain is not known. Here, we used resting-state magnetoencephalography to test the hypothesis that peak alpha frequency (PAF) abnormalities are general features across chronic central and peripheral conditions causing neuropathic pain but exhibit sex-specific differences in networks of the DPC (ascending nociceptive pathway [ANP], default mode network, salience network [SN], and subgenual anterior cingulate cortex). We found that neuropathic pain (N = 25 men and 25 women) was associated with increased PAF power in the DPC compared with 50 age- and sex-matched healthy controls, whereas slower PAF in nodes of the SN (temporoparietal junction) and the ANP (posterior insula) was associated with higher trait pain intensity. In the neuropathic pain group, women exhibited lower PAF power in the subgenual anterior cingulate cortex and faster PAF in the ANP and SN than men. The within-sex analyses indicated that women had neuropathic pain-related increased PAF power in the ANP, SN, and default mode network, whereas men with neuropathic pain had increased PAF power restricted to the ANP. These findings highlight neuropathic pain-related and sex-specific abnormalities in alpha oscillations across the DPC that could underlie aberrant neuronal communication in nociceptive processing and modulation.

State Changes During Resting-State (Magneto)encephalographic Studies: The Effect of Drowsiness on Spectral, Connectivity, and Network Analyses

Background

A common problem in resting-state neuroimaging studies is that subjects become drowsy or fall asleep. Although this could drastically affect neurophysiological measurements, such as magnetoencephalography (MEG), its specific impact remains understudied. We aimed to systematically investigate how often drowsiness is present during resting-state MEG recordings, and how the state changes alter quantitative estimates of oscillatory activity, functional connectivity, and network topology.

Methods

About 8-min MEG recordings of 19 healthy subjects, split into ~13-s epochs, were scored for the presence of eyes-open (EO), alert eyes-closed (A-EC), or drowsy eyes-closed (D-EC) states. After projection to source-space, results of spectral, functional connectivity, and network analyses in 6 canonical frequency bands were compared between these states on a global and regional levels. Functional connectivity was analyzed using the phase lag index (PLI) and corrected amplitude envelope correlation (AECc), and network topology was analyzed using the minimum spanning tree (MST).

Results

Drowsiness was present in >55% of all epochs that did not fulfill the AASM criteria for sleep. There were clear differences in spectral results between the states (A-EC vs. D-EC) and conditions (EO vs. A-EC). The influence of state and condition was far less pronounced for connectivity analyses, with only minimal differences between D-EC and EO in the AECc in the delta band. There were no effects of drowsiness on any of the MST measures.

Conclusions

Drowsiness during eyes-closed resting-state MEG recordings is present in the majority of epochs, despite the instructions to stay awake. This has considerable influence on spectral properties, but much less so on functional connectivity and network topology. These findings are important for interpreting the results of EEG/MEG studies using spectral analyses in neurological disease, where recordings should be evaluated for the presence of drowsiness. For connectivity analyses or studies on network topology, this seems of far less importance.

The natural frequencies of the resting human brain: An MEG-based atlas

Brain oscillations are considered to play a pivotal role in neural communication. However, detailed information regarding the typical oscillatory patterns of individual brain regions is surprisingly scarce. In this study we applied a multivariate data-driven approach to create an atlas of the natural frequencies of the resting human brain on a voxel-by-voxel basis. We analysed resting-state magnetoencephalography (MEG) data from 128 healthy adult volunteers obtained from the Open MEG Archive (OMEGA). Spectral power was computed in source space in 500 ms steps for 82 frequency bins logarithmically spaced from 1.7 to 99.5 Hz. We then applied k-means clustering to detect characteristic spectral profiles and to eventually identify the natural frequency of each voxel. Our results provided empirical confirmation of the canonical frequency bands and revealed a region-specific organisation of intrinsic oscillatory activity, following both a medial-to-lateral and a posterior-to-anterior gradient of increasing frequency. In particular, medial fronto-temporal regions were characterised by slow rhythms (delta/theta). Posterior regions presented natural frequencies in the alpha band, although with differentiated generators in the precuneus and in sensory-specific cortices (i.e., visual and auditory). Somatomotor regions were distinguished by the mu rhythm, while the lateral prefrontal cortex was characterised by oscillations in the high beta range (>20 Hz). Importantly, the brain map of natural frequencies was highly replicable in two independent subsamples of individuals. To the best of our knowledge, this is the most comprehensive atlas of ongoing oscillatory activity performed to date. Critically, the identification of natural frequencies is a fundamental step towards a better understanding of the functional architecture of the human brain.

Network-level permutation entropy of resting-state MEG recordings: A novel biomarker for early-stage Alzheimer's disease?

Increasing evidence suggests that measures of signal variability and complexity could present promising biomarkers for Alzheimer's disease (AD). Earlier studies have however been limited to the characterization of local activity. Here, we investigate whether a network version of permutation entropy could serve as a novel biomarker for early-stage AD. Resting-state source-space magnetoencephalography was recorded in 18 subjects with subjective cognitive decline (SCD) and 18 subjects with mild cognitive impairment (MCI). Local activity was characterized by permutation entropy (PE). Network-level interactions were studied using the inverted joint permutation entropy (JPEinv), corrected for volume conduction. The JPEinv showed a reduction of nonlinear connectivity in MCI subjects in the theta and alpha band. Local PE showed increased theta band entropy. Between-group differences were widespread across brain regions. Receiver operating characteristic (ROC) analysis of classification of MCI versus SCD subjects revealed that a logistic regression model trained on JPEinv features (78.4% [62.5-93.3%]) slightly outperformed PE (76.9% [60.3-93.4%]) and relative theta power-based models (76.9% [60.4-93.3%]). Classification performance of theta JPEinv was at least as good as the relative theta power benchmark. The JPEinv is therefore a potential biomarker for early-stage AD that should be explored in larger studies.

Predicting time-resolved electrophysiological brain networks from structural eigenmodes

How temporal modulations in functional interactions are shaped by the underlying anatomical connections remains an open question. Here, we analyse the role of structural eigenmodes, in the formation and dissolution of temporally evolving functional brain networks using resting-state magnetoencephalography and diffusion magnetic resonance imaging data at the individual subject level. Our results show that even at short timescales, phase and amplitude connectivity can partly be expressed by structural eigenmodes, but hardly by direct structural connections. Albeit a stronger relationship was found between structural eigenmodes and time-resolved amplitude connectivity. Time-resolved connectivity for both phase and amplitude was mostly characterised by a stationary process, superimposed with very brief periods that showed deviations from this stationary process. For these brief periods, dynamic network states were extracted that showed different expressions of eigenmodes. Furthermore, the eigenmode expression was related to overall cognitive performance and co-occurred with fluctuations in community structure of functional networks. These results implicate that ongoing time-resolved resting-state networks, even at short timescales, can to some extent be understood in terms of activation and deactivation of structural eigenmodes and that these eigenmodes play a role in the dynamic integration and segregation of information across the cortex, subserving cognitive functions.

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.

Magnetoencephalography Imaging Reveals Abnormal Information Flow in Temporal Lobe Epilepsy

Background/Introduction: Widespread network disruption has been hypothesized to be an important predictor of outcomes in patients with refractory temporal lobe epilepsy (TLE). Most studies examining functional network disruption in epilepsy have largely focused on the symmetric bidirectional metrics of the strength of network connections. However, a more complete description of network dysfunction impacts in epilepsy requires an investigation of the potentially more sensitive directional metrics of information flow. Methods: This study describes a whole-brain magnetoencephalography-imaging approach to examine resting-state directional information flow networks, quantified by phase-transfer entropy (PTE), in patients with TLE compared with healthy controls (HCs). Associations between PTE and clinical characteristics of epilepsy syndrome are also investigated. Results: Deficits of information flow were specific to alpha-band frequencies. In alpha band, while HCs exhibit a clear posterior-to-anterior directionality of information flow, in patients with TLE, this pattern of regional information outflow and inflow was significantly altered in the frontal and occipital regions. The changes in information flow within the alpha band in selected brain regions were correlated with interictal spike frequency and duration of epilepsy. Conclusions: Impaired information flow is an important dimension of network dysfunction associated with the pathophysiological mechanisms of TLE.

Normative brain mapping of interictal intracranial EEG to localize epileptogenic tissue

The identification of abnormal electrographic activity is important in a wide range of neurological disorders, including epilepsy for localizing epileptogenic tissue. However, this identification may be challenging during non-seizure (interictal) periods, especially if abnormalities are subtle compared to the repertoire of possible healthy brain dynamics. Here, we investigate if such interictal abnormalities become more salient by quantitatively accounting for the range of healthy brain dynamics in a location-specific manner. To this end, we constructed a normative map of brain dynamics, in terms of relative band power, from interictal intracranial recordings from 234 participants (21 598 electrode contacts). We then compared interictal recordings from 62 patients with epilepsy to the normative map to identify abnormal regions. We proposed that if the most abnormal regions were spared by surgery, then patients would be more likely to experience continued seizures postoperatively. We first confirmed that the spatial variations of band power in the normative map across brain regions were consistent with healthy variations reported in the literature. Second, when accounting for the normative variations, regions that were spared by surgery were more abnormal than those resected only in patients with persistent postoperative seizures (t = -3.6, P = 0.0003), confirming our hypothesis. Third, we found that this effect discriminated patient outcomes (area under curve 0.75 P = 0.0003). Normative mapping is a well-established practice in neuroscientific research. Our study suggests that this approach is feasible to detect interictal abnormalities in intracranial EEG, and of potential clinical value to identify pathological tissue in epilepsy. Finally, we make our normative intracranial map publicly available to facilitate future investigations in epilepsy and beyond.