Spatiotemporal signal space separation for regions of interest: Application for extracting neuromagnetic responses evoked by deep brain stimulation.

Magnetoencephalography (MEG) recordings are often contaminated by interference that can exceed the amplitude of physiological brain activity by several orders of magnitude. Furthermore, the activity of interference sources may spatially extend (known as source leakage) into the activity of brain signals of interest, resulting in source estimation inaccuracies. This problem is particularly apparent when using MEG to interrogate the effects of brain stimulation on large-scale cortical networks. In this technical report, we develop a novel denoising approach for suppressing the leakage of interference source activity into the activity representing a brain region of interest. This approach leverages spatial and temporal domain projectors for signal arising from prespecified anatomical regions of interest. We apply this denoising approach to reconstruct simulated evoked response topographies to deep brain stimulation (DBS) in a phantom recording. We highlight the advantages of our approach compared to the benchmark-spatiotemporal signal space separation-and show that it can more accurately reveal brain stimulation-evoked response topographies. Finally, we apply our method to MEG recordings from a single patient with Parkinson's disease, to reveal early cortical-evoked responses to DBS of the subthalamic nucleus.

Multi-modal and multi-model interrogation of large-scale functional brain networks.

Existing whole-brain models are generally tailored to the modelling of a particular data modality (e.g., fMRI or MEG/EEG). We propose that despite the differing aspects of neural activity each modality captures, they originate from shared network dynamics. Building on the universal principles of self-organising delay-coupled nonlinear systems, we aim to link distinct features of brain activity - captured across modalities - to the dynamics unfolding on a macroscopic structural connectome. To jointly predict connectivity, spatiotemporal and transient features of distinct signal modalities, we consider two large-scale models - the Stuart Landau and Wilson and Cowan models - which generate short-lived 40 Hz oscillations with varying levels of realism. To this end, we measure features of functional connectivity and metastable oscillatory modes (MOMs) in fMRI and MEG signals - and compare them against simulated data. We show that both models can represent MEG functional connectivity (FC), functional connectivity dynamics (FCD) and generate MOMs to a comparable degree. This is achieved by adjusting the global coupling and mean conduction time delay and, in the WC model, through the inclusion of balance between excitation and inhibition. For both models, the omission of delays dramatically decreased the performance. For fMRI, the SL model performed worse for FCD and MOMs, highlighting the importance of balanced dynamics for the emergence of spatiotemporal and transient patterns of ultra-slow dynamics. Notably, optimal working points varied across modalities and no model was able to achieve a correlation with empirical FC higher than 0.4 across modalities for the same set of parameters. Nonetheless, both displayed the emergence of FC patterns that extended beyond the constraints of the anatomical structure. Finally, we show that both models can generate MOMs with empirical-like properties such as size (number of brain regions engaging in a mode) and duration (continuous time interval during which a mode appears). Our results demonstrate the emergence of static and dynamic properties of neural activity at different timescales from networks of delay-coupled oscillators at 40 Hz. Given the higher dependence of simulated FC on the underlying structural connectivity, we suggest that mesoscale heterogeneities in neural circuitry may be critical for the emergence of parallel cross-modal functional networks and should be accounted for in future modelling endeavours.

Detection of pathological high-frequency oscillations in refractory epilepsy patients undergoing simultaneous stereo-electroencephalography and magnetoencephalography.

BACKGROUND: Stereo-electroencephalography (SEEG) and magnetoencephalography (MEG) have generally been used independently as part of the pre-surgical evaluation of drug-resistant epilepsy (DRE) patients. However, the possibility of simultaneously employing these recording techniques to determine whether MEG has the potential of offering the same information as SEEG less invasively, or whether it could offer a greater spatial indication of the epileptogenic zone (EZ) to aid surgical planning, has not been previously evaluated. METHODS: Data from 24 paediatric and adult DRE patients, undergoing simultaneous SEEG and MEG as part of their pre-surgical evaluation, was analysed employing manual and automated high-frequency oscillations (HFOs) detection, and spectral and source localisation analyses. RESULTS: Twelve patients (50%) were included in the analysis (4 males; mean age=25.08 years) and showed interictal SEEG and MEG HFOs. HFOs detection was concordant between the two recording modalities, but SEEG displayed higher ability of differentiating between deep and superficial epileptogenic sources. Automated HFO detector in MEG recordings was validated against the manual MEG detection method. Spectral analysis revealed that SEEG and MEG detect distinct epileptic events. The EZ was well correlated with the simultaneously recorded data in 50% patients, while 25% patients displayed poor correlation or discordance. CONCLUSION: MEG recordings can detect HFOs, and simultaneous use of SEEG and MEG HFO identification facilitates EZ localisation during the presurgical planning stage for DRE patients. Further studies are necessary to validate these findings and support the translation of automated HFO detectors into routine clinical practice.

Dynamic sensitivity analysis: Defining personalised strategies to drive brain state transitions via whole brain modelling.

Traditionally, in neuroimaging, model-free analyses are used to find significant differences between brain states via signal detection theory. Depending on the a priori assumptions about the underlying data, different spatio-temporal features can be analysed. Alternatively, model-based techniques infer features from the data and compare significance from model parameters. However, to assess transitions from one brain state to another remains a challenge in current paradigms. Here, we introduce a "Dynamic Sensitivity Analysis" framework that quantifies transitions between brain states in terms of stimulation ability to rebalance spatio-temporal brain activity towards a target state such as healthy brain dynamics. In practice, it means building a whole-brain model fitted to the spatio-temporal description of brain dynamics, and applying systematic stimulations in-silico to assess the optimal strategy to drive brain dynamics towards a target state. Further, we show how Dynamic Sensitivity Analysis extends to various brain stimulation paradigms, ultimately contributing to improving the efficacy of personalised clinical interventions.

Functional connectivity maps of theta/alpha and beta coherence within the subthalamic nucleus region.

The subthalamic nucleus (STN) is a primary target for deep brain stimulation in Parkinson's disease (PD). Although small in size, the STN is commonly partitioned into sensorimotor, cognitive/associative, and limbic subregions based on its structural connectivity profile to cortical areas. We investigated whether such a regional specialization is also supported by functional connectivity between local field potential recordings and simultaneous magnetoencephalography. Using a novel data set of 21 PD patients, we replicated previously reported cortico-STN coherence networks in the theta/alpha and beta frequency ranges, and looked for the spatial distribution of these networks within the STN region. Although theta/alpha and beta coherence peaks were both observed in on-medication recordings from electrode contacts at several locations within and around the STN, sites with theta/alpha coherence peaks were situated at significantly more inferior MNI coordinates than beta coherence peaks. Sites with only theta/alpha coherence peaks, i.e. without distinct beta coherence, were mostly located near the border of sensorimotor and cognitive/associative subregions as defined by a tractography-based atlas of the STN. Peak coherence values were largely unaltered by the medication state of the subject, however, theta/alpha peaks were more often identified in recordings obtained after administration of dopaminergic medication. Our findings suggest the existence of a frequency-specific topography of cortico-STN coherence within the STN, albeit with considerable spatial overlap between functional networks. Consequently, optimization of deep brain stimulation targeting might remain a trade-off between alleviating motor symptoms and avoiding adverse neuropsychiatric side effects.

Conflict Detection in a Sequential Decision Task Is Associated with Increased Cortico-Subthalamic Coherence and Prolonged Subthalamic Oscillatory Response in the ß Band.

Making accurate decisions often involves the integration of current and past evidence. Here, we examine the neural correlates of conflict and evidence integration during sequential decision-making. Female and male human patients implanted with deep-brain stimulation (DBS) electrodes and age-matched and gender-matched healthy controls performed an expanded judgment task, in which they were free to choose how many cues to sample. Behaviorally, we found that while patients sampled numerically more cues, they were less able to integrate evidence and showed suboptimal performance. Using recordings of magnetoencephalography (MEG) and local field potentials (LFPs; in patients) in the subthalamic nucleus (STN), we found that ß oscillations signaled conflict between cues within a sequence. Following cues that differed from previous cues, ß power in the STN and cortex first decreased and then increased. Importantly, the conflict signal in the STN outlasted the cortical one, carrying over to the next cue in the sequence. Furthermore, after a conflict, there was an increase in coherence between the dorsal premotor cortex and STN in the ß band. These results extend our understanding of cortico-subcortical dynamics of conflict processing, and do so in a context where evidence must be accumulated in discrete steps, much like in real life. Thus, the present work leads to a more nuanced picture of conflict monitoring systems in the brain and potential changes because of disease.SIGNIFICANCE STATEMENT Decision-making often involves the integration of multiple pieces of information over time to make accurate predictions. We simultaneously recorded whole-head magnetoencephalography (MEG) and local field potentials (LFPs) from the human subthalamic nucleus (STN) in a novel task which required integrating sequentially presented pieces of evidence. Our key finding is prolonged ß oscillations in the STN, with a concurrent increase in communication with frontal cortex, when presented with conflicting information. These neural effects reflect the behavioral profile of reduced tendency to respond after conflict, as well as relate to suboptimal cue integration in patients, which may be directly linked to clinically reported side-effects of deep-brain stimulation (DBS) such as impaired decision-making and impulsivity.

Separating Neural Oscillations from Aperiodic 1/f Activity: Challenges and Recommendations.

Electrophysiological power spectra typically consist of two components: An aperiodic part usually following an 1/f power law [Formula: see text] and periodic components appearing as spectral peaks. While the investigation of the periodic parts, commonly referred to as neural oscillations, has received considerable attention, the study of the aperiodic part has only recently gained more interest. The periodic part is usually quantified by center frequencies, powers, and bandwidths, while the aperiodic part is parameterized by the y-intercept and the 1/f exponent [Formula: see text]. For investigation of either part, however, it is essential to separate the two components. In this article, we scrutinize two frequently used methods, FOOOF (Fitting Oscillations & One-Over-F) and IRASA (Irregular Resampling Auto-Spectral Analysis), that are commonly used to separate the periodic from the aperiodic component. We evaluate these methods using diverse spectra obtained with electroencephalography (EEG), magnetoencephalography (MEG), and local field potential (LFP) recordings relating to three independent research datasets. Each method and each dataset poses distinct challenges for the extraction of both spectral parts. The specific spectral features hindering the periodic and aperiodic separation are highlighted by simulations of power spectra emphasizing these features. Through comparison with the simulation parameters defined a priori, the parameterization error of each method is quantified. Based on the real and simulated power spectra, we evaluate the advantages of both methods, discuss common challenges, note which spectral features impede the separation, assess the computational costs, and propose recommendations on how to use them.

Stimulating at the right time to recover network states in a model of the cortico-basal ganglia-thalamic circuit.

Synchronization of neural oscillations is thought to facilitate communication in the brain. Neurodegenerative pathologies such as Parkinson's disease (PD) can result in synaptic reorganization of the motor circuit, leading to altered neuronal dynamics and impaired neural communication. Treatments for PD aim to restore network function via pharmacological means such as dopamine replacement, or by suppressing pathological oscillations with deep brain stimulation. We tested the hypothesis that brain stimulation can operate beyond a simple "reversible lesion" effect to augment network communication. Specifically, we examined the modulation of beta band (14-30 Hz) activity, a known biomarker of motor deficits and potential control signal for stimulation in Parkinson's. To do this we setup a neural mass model of population activity within the cortico-basal ganglia-thalamic (CBGT) circuit with parameters that were constrained to yield spectral features comparable to those in experimental Parkinsonism. We modulated the connectivity of two major pathways known to be disrupted in PD and constructed statistical summaries of the spectra and functional connectivity of the resulting spontaneous activity. These were then used to assess the network-wide outcomes of closed-loop stimulation delivered to motor cortex and phase locked to subthalamic beta activity. Our results demonstrate that the spatial pattern of beta synchrony is dependent upon the strength of inputs to the STN. Precisely timed stimulation has the capacity to recover network states, with stimulation phase inducing activity with distinct spectral and spatial properties. These results provide a theoretical basis for the design of the next-generation brain stimulators that aim to restore neural communication in disease.

Spontaneous transient states of fronto-temporal and default-mode networks altered by suicide attempt in major depressive disorder.

Major depressive disorder (MDD) is associated with increased suicidality, and it's still challenging to identify suicide in clinical practice. Although suicide attempt (SA) is the most relevant precursor with multiple functional abnormalities reported from neuroimaging studies, little is known about how the spontaneous transient activated patterns organize and coordinate brain networks underlying SA. Thus, we obtained resting-state magnetoencephalography data for two MDD subgroups of 44 non-suicide patients and 34 suicide-attempted patients, together with 49 matched health-controls. For the source-space signals, Hidden Markov Model (HMM) helped to capture the sub-second dynamic activity via a hidden sequence of finite number of states. Temporal parameters and spectral activation were acquired for each state and then compared between groups. Here, HMM states characterized the spatiotemporal signatures of eight networks. The activity of suicide attempters switches more frequently into the fronto-temporal network, as the time spent occupancy of fronto-temporal state is increased and interval time is decreased compared with the non-suicide patients. Moreover, these changes are significantly correlated with Nurses' Global Assessment of Suicide Risk scores. Suicide attempters also exhibit increased state-wise activations in the theta band (4-8 Hz) in the posterior default mode network centered on posterior cingulate cortex, which can't be detected in the static spectral analysis. These alternations may disturb the time allocations of cognitive control regulations and cause inflexible decision making to SA. As the better sensitivity of dynamic study in reflecting SA diathesis than the static is validated, dynamic stability could serve as a potential neuronal marker for SA.

Balance between competing spectral states in subthalamic nucleus is linked to motor impairment in Parkinson's disease.

Exaggerated local field potential bursts of activity at frequencies in the low beta band are a well-established phenomenon in the subthalamic nucleus of patients with Parkinson's disease. However, such activity is only moderately correlated with motor impairment. Here we test the hypothesis that beta bursts are just one of several dynamic states in the subthalamic nucleus local field potential in Parkinson's disease, and that together these different states predict motor impairment with high fidelity. Local field potentials were recorded in 32 patients (64 hemispheres) undergoing deep brain stimulation surgery targeting the subthalamic nucleus. Recordings were performed following overnight withdrawal of anti-parkinsonian medication, and after administration of levodopa. Local field potentials were analysed using hidden Markov modelling to identify transient spectral states with frequencies under 40 Hz. Findings in the low beta frequency band were similar to those previously reported; levodopa reduced occurrence rate and duration of low beta states, and the greater the reductions, the greater the improvement in motor impairment. However, additional local field potential states were distinguished in the theta, alpha and high beta bands, and these behaved in an opposite manner. They were increased in occurrence rate and duration by levodopa, and the greater the increases, the greater the improvement in motor impairment. In addition, levodopa favoured the transition of low beta states to other spectral states. When all local field potential states and corresponding features were considered in a multivariate model it was possible to predict 50% of the variance in patients' hemibody impairment OFF medication, and in the change in hemibody impairment following levodopa. This only improved slightly if signal amplitude or gamma band features were also included in the multivariate model. In addition, it compares with a prediction of only 16% of the variance when using beta bursts alone. We conclude that multiple spectral states in the subthalamic nucleus local field potential have a bearing on motor impairment, and that levodopa-induced shifts in the balance between these states can predict clinical change with high fidelity. This is important in suggesting that some states might be upregulated to improve parkinsonism and in suggesting how local field potential feedback can be made more informative in closed-loop deep brain stimulation systems.

A unified view on beamformers for M/EEG source reconstruction.

Beamforming is a popular method for functional source reconstruction using magnetoencephalography (MEG) and electroencephalography (EEG) data. Beamformers, which were first proposed for MEG more than two decades ago, have since been applied in hundreds of studies, demonstrating that they are a versatile and robust tool for neuroscience. However, certain characteristics of beamformers remain somewhat elusive and there currently does not exist a unified documentation of the mathematical underpinnings and computational subtleties of beamformers as implemented in the most widely used academic open source software packages for MEG analysis (Brainstorm, FieldTrip, MNE, and SPM). Here, we provide such documentation that aims at providing the mathematical background of beamforming and unifying the terminology. Beamformer implementations are compared across toolboxes and pitfalls of beamforming analyses are discussed. Specifically, we provide details on handling rank deficient covariance matrices, prewhitening, the rank reduction of forward fields, and on the combination of heterogeneous sensor types, such as magnetometers and gradiometers. The overall aim of this paper is to contribute to contemporary efforts towards higher levels of computational transparency in functional neuroimaging.

Metastable oscillatory modes emerge from synchronization in the brain spacetime connectome.

A rich repertoire of oscillatory signals is detected from human brains with electro- and magnetoencephalography (EEG/MEG). However, the principles underwriting coherent oscillations and their link with neural activity remain under debate. Here, we revisit the mechanistic hypothesis that transient brain rhythms are a signature of metastable synchronization, occurring at reduced collective frequencies due to delays between brain areas. We consider a system of damped oscillators in the presence of background noise - approximating the short-lived gamma-frequency oscillations generated within neuronal circuits - coupled according to the diffusion weighted tractography between brain areas. Varying the global coupling strength and conduction speed, we identify a critical regime where spatially and spectrally resolved metastable oscillatory modes (MOMs) emerge at sub-gamma frequencies, approximating the MEG power spectra from 89 healthy individuals at rest. Further, we demonstrate that the frequency, duration, and scale of MOMs - as well as the frequency-specific envelope functional connectivity - can be controlled by global parameters, while the connectome structure remains unchanged. Grounded in the physics of delay-coupled oscillators, these numerical analyses demonstrate how interactions between locally generated fast oscillations in the connectome spacetime structure can lead to the emergence of collective brain rhythms organized in space and time.

Neural signatures of hyperdirect pathway activity in Parkinson's disease.

Parkinson's disease (PD) is characterised by the emergence of beta frequency oscillatory synchronisation across the cortico-basal-ganglia circuit. The relationship between the anatomy of this circuit and oscillatory synchronisation within it remains unclear. We address this by combining recordings from human subthalamic nucleus (STN) and internal globus pallidus (GPi) with magnetoencephalography, tractography and computational modelling. Coherence between supplementary motor area and STN within the high (21-30 Hz) but not low (13-21 Hz) beta frequency range correlated with 'hyperdirect pathway' fibre densities between these structures. Furthermore, supplementary motor area activity drove STN activity selectively at high beta frequencies suggesting that high beta frequencies propagate from the cortex to the basal ganglia via the hyperdirect pathway. Computational modelling revealed that exaggerated high beta hyperdirect pathway activity can provoke the generation of widespread pathological synchrony at lower beta frequencies. These findings suggest a spectral signature and a pathophysiological role for the hyperdirect pathway in PD.

The missing N1 or jittered P2: Electrophysiological correlates of pattern glare in the time and frequency domain.

Excessive sensitivity to certain visual stimuli (cortical hyperexcitability) is associated with a number of neurological disorders including migraine, epilepsy, multiple sclerosis, autism and possibly dyslexia. Others show disruptive sensitivity to visual stimuli with no other obvious pathology or symptom profile (visual stress) which can extend to discomfort and nausea. We used event-related potentials (ERPs) to explore the neural correlates of visual stress and headache proneness. We analysed ERPs in response to thick (0.37 cycles per degree [c/deg]), medium (3 c/deg) and thin (12 c/deg) gratings, using mass univariate analysis, considering three factors in the general population: headache proneness, visual stress and discomfort. We found relationships between ERP features and the headache and discomfort factors. Stimulus main effects were driven by the medium stimulus regardless of participant characteristics. Participants with high discomfort ratings had larger P1 components for the initial presentation of medium stimuli, suggesting initial cortical hyperexcitability that is later suppressed. The participants with high headache ratings showed atypical N1-P2 components for medium stripes relative to the other stimuli. This effect was present only after repeated stimulus presentation. These effects were also explored in the frequency domain, suggesting variations in intertrial theta band phase coherence. Our results suggest that discomfort and headache in response to striped stimuli are related to different neural processes; however, more exploration is needed to determine whether the results translate to a clinical migraine population.

Watching Movies Unfold, a Frame-by-Frame Analysis of the Associated Neural Dynamics.

Our lives unfold as sequences of events. We experience these events as seamless, although they are composed of individual images captured in between the interruptions imposed by eye blinks and saccades. Events typically involve visual imagery from the real world (scenes), and the hippocampus is frequently engaged in this context. It is unclear, however, whether the hippocampus would be similarly responsive to unfolding events that involve abstract imagery. Addressing this issue could provide insights into the nature of its contribution to event processing, with relevance for theories of hippocampal function. Consequently, during magnetoencephalography (MEG), we had female and male humans watch highly matched unfolding movie events composed of either scene image frames that reflected the real world, or frames depicting abstract patterns. We examined the evoked neuronal responses to each image frame along the time course of the movie events. Only one difference between the two conditions was evident, and that was during the viewing of the first image frame of events, detectable across frontotemporal sensors. Further probing of this difference using source reconstruction revealed greater engagement of a set of brain regions across parietal, frontal, premotor, and cerebellar cortices, with the largest change in broadband (1-30 Hz) power in the hippocampus during scene-based movie events. Hippocampal engagement during the first image frame of scene-based events could reflect its role in registering a recognizable context perhaps based on templates or schemas. The hippocampus, therefore, may help to set the scene for events very early on.

The use of simultaneous stereo-electroencephalography and magnetoencephalography in localizing the epileptogenic focus in refractory focal epilepsy.

Both magnetoencephalography and stereo-electroencephalography are used in presurgical epilepsy assessment, with contrasting advantages and limitations. It is not known whether simultaneous stereo-electroencephalography-magnetoencephalography recording confers an advantage over both individual modalities, in particular whether magnetoencephalography can provide spatial context to epileptiform activity seen on stereo-electroencephalography. Twenty-four adult and paediatric patients who underwent stereo-electroencephalography study for pre-surgical evaluation of drug-resistant focal epilepsy, were recorded using simultaneous stereo-electroencephalography-magnetoencephalography, of which 14 had abnormal interictal activity during recording. The 14 patients were divided into two groups; those with detected superficial (n = 7) and deep (n = 7) brain interictal activity. Interictal spikes were independently identified in stereo-electroencephalography and magnetoencephalography. Magnetoencephalography dipoles were derived using a distributed inverse method. There was no significant difference between stereo-electroencephalography and magnetoencephalography in detecting superficial spikes (P = 0.135) and stereo-electroencephalography was significantly better at detecting deep spikes (P = 0.002). Mean distance across patients between stereo-electroencephalography channel with highest average spike amplitude and magnetoencephalography dipole was 20.7 ± 4.4 mm. for superficial sources, and 17.8 ± 3.7 mm. for deep sources, even though for some of the latter (n = 4) no magnetoencephalography spikes were detected and magnetoencephalography dipole was fitted to a stereo-electroencephalography interictal activity triggered average. Removal of magnetoencephalography dipole was associated with 1 year seizure freedom in 6/7 patients with superficial source, and 5/6 patients with deep source. Although stereo-electroencephalography has greater sensitivity in identifying interictal activity from deeper sources, a magnetoencephalography source can be localized using stereo-electroencephalography information, thereby providing useful whole brain context to stereo-electroencephalography and potential role in epilepsy surgery planning.

Inference of brain networks with approximate Bayesian computation - assessing face validity with an example application in Parkinsonism.

This paper describes and validates a novel framework using the Approximate Bayesian Computation (ABC) algorithm for parameter estimation and model selection in models of mesoscale brain network activity. We provide a proof of principle, first pass validation of this framework using a set of neural mass models of the cortico-basal ganglia thalamic circuit inverted upon spectral features from experimental, in vivo recordings. This optimization scheme relaxes an assumption of fixed-form posteriors (i.e. the Laplace approximation) taken in previous approaches to inverse modelling of spectral features. This enables the exploration of model dynamics beyond that approximated from local linearity assumptions and so fit to explicit, numerical solutions of the underlying non-linear system of equations. In this first paper, we establish a face validation of the optimization procedures in terms of: (i) the ability to approximate posterior densities over parameters that are plausible given the known causes of the data; (ii) the ability of the model comparison procedures to yield posterior model probabilities that can identify the model structure known to generate the data; and (iii) the robustness of these procedures to local minima in the face of different starting conditions. Finally, as an illustrative application we show (iv) that model comparison can yield plausible conclusions given the known neurobiology of the cortico-basal ganglia-thalamic circuit in Parkinsonism. These results lay the groundwork for future studies utilizing highly nonlinear or brittle models that can explain time dependant dynamics, such as oscillatory bursts, in terms of the underlying neural circuits.

Dynamic analysis on simultaneous iEEG-MEG data via hidden Markov model.

BACKGROUND: Intracranial electroencephalography (iEEG) recordings are used for clinical evaluation prior to surgical resection of the focus of epileptic seizures and also provide a window into normal brain function. A major difficulty with interpreting iEEG results at the group level is inconsistent placement of electrodes between subjects making it difficult to select contacts that correspond to the same functional areas. Recent work using time delay embedded hidden Markov model (HMM) applied to magnetoencephalography (MEG) resting data revealed a distinct set of brain states with each state engaging a specific set of cortical regions. Here we use a rare group dataset with simultaneously acquired resting iEEG and MEG to test whether there is correspondence between HMM states and iEEG power changes that would allow classifying iEEG contacts into functional clusters. METHODS: Simultaneous MEG-iEEG recordings were performed at rest on 11 patients with epilepsy whose intracranial electrodes were implanted for pre-surgical evaluation. Pre-processed MEG sensor data was projected to source space. Time delay embedded HMM was then applied to MEG time series. At the same time, iEEG time series were analyzed with time-frequency decomposition to obtain spectral power changes with time. To relate MEG and iEEG results, correlations were computed between HMM probability time courses of state activation and iEEG power time course from the mid contact pair for each electrode in equally spaced frequency bins and presented as correlation spectra for the respective states and iEEG channels. Association of iEEG electrodes with HMM states based on significant correlations was compared to that based on the distance to peaks in subject-specific state topographies. RESULTS: Five HMM states were inferred from MEG. Two of them corresponded to the left and the right temporal activations and had a spectral signature primarily in the theta/alpha frequency band. All the electrodes had significant correlations with at least one of the states (p < 0.05 uncorrected) and for 27/50 electrodes these survived within-subject FDR correction (q < 0.05). These correlations peaked in the theta/alpha band. There was a highly significant dependence between the association of states and electrodes based on functional correlations and that based on spatial proximity (p = 5.6e-6,χ2 test for independence). Despite the potentially atypical functional anatomy and physiological abnormalities related to epilepsy, HMM model estimated from the patient group was very similar to that estimated from healthy subjects. CONCLUSION: Epilepsy does not preclude HMM analysis of interictal data. The resulting group functional states are highly similar to those reported for healthy controls. Power changes recorded with iEEG correlate with HMM state time courses in the alpha-theta band and the presence of this correlation can be related to the spatial location of electrode contacts close to the individual peaks of the corresponding state topographies. Thus, the hypothesized relation between iEEG contacts and HMM states exists and HMM could be further explored as a method for identifying comparable iEEG channels across subjects for the purposes of group analysis.