Parietal cortical alpha/beta suppression during prospective memory retrieval.

Prospective memory (PM) impairment is among the most frequent memory complaints, yet little is known about the underlying neural mechanisms. PM for a planned intention may be achieved through strategic monitoring of the environment for cues, involving ongoing attentional processes, or through spontaneous retrieval. We hypothesized that parietal spectral power modulation accompanies prospectively encoded intention retrieval, irrespective of PM retrieval approach. A cognitively engaging arithmetic-based ongoing task (OGT) was employed to encourage spontaneous retrieval, with a focal, internally generated PM cue to eliminate OGT/PM trial differentiation based on perceptual or conceptual PM cue features. Two PM repetition frequencies were used to vary the extent of strategic monitoring. We observed a transient parietal alpha/beta spectral power reduction directly preceding the response, which was distinguishable on a single trial basis, as revealed by an OGT/PM trial classification rate exceeding 70% using linear discriminant analysis. The alpha/beta idling rhythm reflects cortical inhibition. A disengagement of task-relevant neural assemblies from this rhythm, reflected in alpha/beta power reduction, is deemed to increase information content, facilitate information integration, and enable engagement of neural assemblies in task-related cortical networks. The observed power reduction is consistent with the Dual Pathways model, where PM strategies converge at the PM retrieval stage.

Deep brain stimulation of the ventrointermediate nucleus of the thalamus to treat essential tremor improves motor sequence learning.

The network of brain structures engaged in motor sequence learning comprises the same structures as those involved in tremor, including basal ganglia, cerebellum, thalamus, and motor cortex. Deep brain stimulation (DBS) of the ventrointermediate nucleus of the thalamus (VIM) reduces tremor, but the effects on motor sequence learning are unknown. We investigated whether VIM stimulation has an impact on motor sequence learning and hypothesized that stimulation effects depend on the laterality of electrode location. Twenty patients (age: 38-81 years; 12 female) with VIM electrodes implanted to treat essential tremor (ET) successfully performed a serial reaction time task, varying whether the stimuli followed a repeating pattern or were selected at random, during which VIM-DBS was either on or off. Analyses of variance were applied to evaluate motor sequence learning performance according to reaction times (RTs) and accuracy. An interaction was observed between whether the sequence was repeated or random and whether VIM-DBS was on or off (F[1,18] = 7.89, p = .012). Motor sequence learning, reflected by reduced RTs for repeated sequences, was greater with DBS on than off (T[19] = 2.34, p = .031). Stimulation location correlated with the degree of motor learning, with greater motor learning when stimulation targeted the lateral VIM (n = 23, ρ = 0.46; p = .027). These results demonstrate the beneficial effects of VIM-DBS on motor sequence learning in ET patients, particularly with lateral VIM electrode location, and provide evidence for a role for the VIM in motor sequence learning.

Music-induced emotions can be predicted from a combination of brain activity and acoustic features.

It is widely acknowledged that music can communicate and induce a wide range of emotions in the listener. However, music is a highly-complex audio signal composed of a wide range of complex time- and frequency-varying components. Additionally, music-induced emotions are known to differ greatly between listeners. Therefore, it is not immediately clear what emotions will be induced in a given individual by a piece of music. We attempt to predict the music-induced emotional response in a listener by measuring the activity in the listeners electroencephalogram (EEG). We combine these measures with acoustic descriptors of the music, an approach that allows us to consider music as a complex set of time-varying acoustic features, independently of any specific music theory. Regression models are found which allow us to predict the music-induced emotions of our participants with a correlation between the actual and predicted responses of up to r=0.234,p<0.001. This regression fit suggests that over 20% of the variance of the participant's music induced emotions can be predicted by their neural activity and the properties of the music. Given the large amount of noise, non-stationarity, and non-linearity in both EEG and music, this is an encouraging result. Additionally, the combination of measures of brain activity and acoustic features describing the music played to our participants allows us to predict music-induced emotions with significantly higher accuracies than either feature type alone (p<0.01).

Ventrointermediate thalamic stimulation improves motor learning in humans.

Ventrointermediate thalamic stimulation (VIM-DBS) modulates oscillatory activity in a cortical network including primary motor cortex, premotor cortex, and parietal cortex. Here we show that, beyond the beneficial effects of VIM-DBS on motor execution, this form of invasive stimulation facilitates production of sequential finger movements that follow a repeated sequence. These results highlight the role of thalamo-cortical activity in motor learning.

Hebbian plasticity induced by temporally coincident BCI enhances post-stroke motor recovery.

Functional electrical stimulation (FES) can support functional restoration of a paretic limb post-stroke. Hebbian plasticity depends on temporally coinciding pre- and post-synaptic activity. A tight temporal relationship between motor cortical (MC) activity associated with attempted movement and FES-generated visuo-proprioceptive feedback is hypothesized to enhance motor recovery. Using a brain-computer interface (BCI) to classify MC spectral power in electroencephalographic (EEG) signals to trigger FES-delivery with detection of movement attempts improved motor outcomes in chronic stroke patients. We hypothesized that heightened neural plasticity earlier post-stroke would further enhance corticomuscular functional connectivity and motor recovery. We compared subcortical non-dominant hemisphere stroke patients in BCI-FES and Random-FES (FES temporally independent of MC movement attempt detection) groups. The primary outcome measure was the Fugl-Meyer Assessment, Upper Extremity (FMA-UE). We recorded high-density EEG and transcranial magnetic stimulation-induced motor evoked potentials before and after treatment. The BCI group showed greater: FMA-UE improvement; motor evoked potential amplitude; beta oscillatory power and long-range temporal correlation reduction over contralateral MC; and corticomuscular coherence with contralateral MC. These changes are consistent with enhanced post-stroke motor improvement when movement is synchronized with MC activity reflecting attempted movement.

Endogenous cholinergic tone modulates spontaneous network level neuronal activity in primary cortical cultures grown on multi-electrode arrays.

BACKGROUND: Cortical cultures grown long-term on multi-electrode arrays (MEAs) are frequently and extensively used as models of cortical networks in studies of neuronal firing activity, neuropharmacology, toxicology and mechanisms underlying synaptic plasticity. However, in contrast to the predominantly asynchronous neuronal firing activity exhibited by intact cortex, electrophysiological activity of mature cortical cultures is dominated by spontaneous epileptiform-like global burst events which hinders their effective use in network-level studies, particularly for neurally-controlled animat ('artificial animal') applications. Thus, the identification of culture features that can be exploited to produce neuronal activity more representative of that seen in vivo could increase the utility and relevance of studies that employ these preparations. Acetylcholine has a recognised neuromodulatory role affecting excitability, rhythmicity, plasticity and information flow in vivo although its endogenous production by cortical cultures and subsequent functional influence upon neuronal excitability remains unknown. RESULTS: Consequently, using MEA electrophysiological recording supported by immunohistochemical and RT-qPCR methods, we demonstrate for the first time, the presence of intrinsic cholinergic neurons and significant, endogenous cholinergic tone in cortical cultures with a characterisation of the muscarinic and nicotinic components that underlie modulation of spontaneous neuronal activity. We found that tonic muscarinic ACh receptor (mAChR) activation affects global excitability and burst event regularity in a culture age-dependent manner whilst, in contrast, tonic nicotinic ACh receptor (nAChR) activation can modulate burst duration and the proportion of spikes occurring within bursts in a spatio-temporal fashion. CONCLUSIONS: We suggest that the presence of significant endogenous cholinergic tone in cortical cultures and the comparability of its modulatory effects to those seen in intact brain tissues support emerging, exploitable commonalities between in vivo and in vitro preparations. We conclude that experimental manipulation of endogenous cholinergic tone could offer a novel opportunity to improve the use of cortical cultures for studies of network-level mechanisms in a manner that remains largely consistent with its functional role.

Testing for significance of phase synchronisation dynamics in the EEG.

A number of tests exist to check for statistical significance of phase synchronisation within the Electroencephalogram (EEG); however, the majority suffer from a lack of generality and applicability. They may also fail to account for temporal dynamics in the phase synchronisation, regarding synchronisation as a constant state instead of a dynamical process. Therefore, a novel test is developed for identifying the statistical significance of phase synchronisation based upon a combination of work characterising temporal dynamics of multivariate time-series and Markov modelling. We show how this method is better able to assess the significance of phase synchronisation than a range of commonly used significance tests. We also show how the method may be applied to identify and classify significantly different phase synchronisation dynamics in both univariate and multivariate datasets.

Emergence of a small-world functional network in cultured neurons.

The functional networks of cultured neurons exhibit complex network properties similar to those found in vivo. Starting from random seeding, cultures undergo significant reorganization during the initial period in vitro, yet despite providing an ideal platform for observing developmental changes in neuronal connectivity, little is known about how a complex functional network evolves from isolated neurons. In the present study, evolution of functional connectivity was estimated from correlations of spontaneous activity. Network properties were quantified using complex measures from graph theory and used to compare cultures at different stages of development during the first 5 weeks in vitro. Networks obtained from young cultures (14 days in vitro) exhibited a random topology, which evolved to a small-world topology during maturation. The topology change was accompanied by an increased presence of highly connected areas (hubs) and network efficiency increased with age. The small-world topology balances integration of network areas with segregation of specialized processing units. The emergence of such network structure in cultured neurons, despite a lack of external input, points to complex intrinsic biological mechanisms. Moreover, the functional network of cultures at mature ages is efficient and highly suited to complex processing tasks.

Assessment of inter-examiner agreement and variability in the manual classification of auditory brainstem response.

BACKGROUND: The analysis of the Auditory Brainstem Response (ABR) is of fundamental importance to the investigation of the auditory system behaviour, though its interpretation has a subjective nature because of the manual process employed in its study and the clinical experience required for its analysis. When analysing the ABR, clinicians are often interested in the identification of ABR signal components referred to as Jewett waves. In particular, the detection and study of the time when these waves occur (i.e., the wave latency) is a practical tool for the diagnosis of disorders affecting the auditory system. Significant differences in inter-examiner results may lead to completely distinct clinical interpretations of the state of the auditory system. In this context, the aim of this research was to evaluate the inter-examiner agreement and variability in the manual classification of ABR. METHODS: A total of 160 ABR data samples were collected, for four different stimulus intensity (80dBHL, 60dBHL, 40dBHL and 20dBHL), from 10 normal-hearing subjects (5 men and 5 women, from 20 to 52 years). Four examiners with expertise in the manual classification of ABR components participated in the study. The Bland-Altman statistical method was employed for the assessment of inter-examiner agreement and variability. The mean, standard deviation and error for the bias, which is the difference between examiners' annotations, were estimated for each pair of examiners. Scatter plots and histograms were employed for data visualization and analysis. RESULTS: In most comparisons the differences between examiner's annotations were below 0.1 ms, which is clinically acceptable. In four cases, it was found a large error and standard deviation (>0.1 ms) that indicate the presence of outliers and thus, discrepancies between examiners. CONCLUSIONS: Our results quantify the inter-examiner agreement and variability of the manual analysis of ABR data, and they also allows for the determination of different patterns of manual ABR analysis.

Brain and brain-heart Granger causality during wakefulness and sleep.

In this exploratory study we apply Granger Causality (GC) to investigate the brain-brain and brain-heart interactions during wakefulness and sleep. Our analysis includes electroencephalogram (EEG) and electrocardiogram (ECG) data during all-night polysomnographic recordings from volunteers with apnea, available from the Massachusetts General Hospital's Computational Clinical Neurophysiology Laboratory and the Clinical Data Animation Laboratory. The data is manually annotated by clinical staff at the MGH in 30 second contiguous intervals (wakefulness and sleep stages 1, 2, 3, and rapid eye movement (REM). We applied GC to 4-s non-overlapping segments of available EEG and ECG across all-night recordings of 50 randomly chosen patients. To identify differences in GC between the different sleep stages, the GC for each sleep stage was subtracted from the GC during wakefulness. Positive (negative) differences indicated that GC was greater (lower) during wakefulness compared to the specific sleep stage. The application of GC to study brain-brain and brain-heart bidirectional connections during wakefulness and sleep confirmed the importance of fronto-posterior connectivity during these two states, but has also revealed differences in ipsilateral and contralateral mechanisms of these connections. It has also confirmed the existence of bidirectional brain-heart connections that are more prominent in the direction from brain to heart. Our exploratory study has shown that GC can be successfully applied to sleep data analysis and captures the varying physiological mechanisms that are related to wakefulness and different sleep stages.

Subthalamic low beta bursts differ in Parkinson's disease phenotypes.

OBJECTIVE: Parkinson's disease (PD) patients may be categorized into tremor-dominant (TD) and postural-instability and gait disorder (PIGD) motor phenotypes, but the dynamical aspects of subthalamic nucleus local field potentials (STN-LFP) and the neural correlates of this phenotypical classification remain unclear. METHODS: 35 STN-LFP (20 PIGD and 15 TD) were investigated through continuous wavelet transform and machine-learning-based methods. The beta oscillation - the main band associated with motor impairment in PD - dynamics was characterized through beta burst parameters across phenotypes and burst intervals under specific proposed criteria for optimal burst threshold definition. RESULTS: Low-frequency (13-22 Hz) beta burst probability was the best predictor for PD phenotypes (75% accuracy). PIGD patients presented higher average burst duration (p = 0.018), while TD patients exhibited higher burst probability (p = 0.014). Categorization into shorter and longer than 400 ms bursts led to significant interaction between burst length categories and the phenotypes (p < 0.050) as revealed by mixed-effects models. Long burst durations and short bursts probability positively correlated, respectively, with rigidity-bradykinesia (p = 0.029) and tremor (p = 0.038) scores. CONCLUSIONS: Subthalamic low-frequency beta bursts differed between TD and PIGD phenotypes and correlated with motor symptoms. SIGNIFICANCE: These findings improve the PD phenotypes' electrophysiological characterization and may define new criteria for adaptive deep brain stimulation.

Dynamics of Long-Range Temporal Correlations in Broadband EEG During Different Motor Execution and Imagery Tasks.

Brain activity is composed of oscillatory and broadband arrhythmic components; however, there is more focus on oscillatory sensorimotor rhythms to study movement, but temporal dynamics of broadband arrhythmic electroencephalography (EEG) remain unexplored. We have previously demonstrated that broadband arrhythmic EEG contains both short- and long-range temporal correlations that change significantly during movement. In this study, we build upon our previous work to gain a deeper understanding of these changes in the long-range temporal correlation (LRTC) in broadband EEG and contrast them with the well-known LRTC in alpha oscillation amplitude typically found in the literature. We investigate and validate changes in LRTCs during five different types of movements and motor imagery tasks using two independent EEG datasets recorded with two different paradigms-our finger tapping dataset with single self-initiated asynchronous finger taps and publicly available EEG dataset containing cued continuous movement and motor imagery of fists and feet. We quantified instantaneous changes in broadband LRTCs by detrended fluctuation analysis on single trial 2 s EEG sliding windows. The broadband LRTC increased significantly (p < 0.05) during all motor tasks as compared to the resting state. In contrast, the alpha oscillation LRTC, which had to be computed on longer stitched EEG segments, decreased significantly (p < 0.05) consistently with the literature. This suggests the complementarity of underlying fast and slow neuronal scale-free dynamics during movement and motor imagery. The single trial broadband LRTC gave high average binary classification accuracy in the range of 70.54±10.03% to 76.07±6.40% for all motor execution and imagery tasks and hence can be used in brain-computer interface (BCI). Thus, we demonstrate generalizability, robustness, and reproducibility of novel motor neural correlate, the single trial broadband LRTC, during different motor execution and imagery tasks in single asynchronous and cued continuous motor-BCI paradigms and its contrasting behavior with LRTC in alpha oscillation amplitude.

Synchronization-based control for a collaborative robot.

This article introduces a new control scheme for controlling a robotic manipulator in a collaborative task, allowing it to respond proactively to its partner's movements. Unlike conventional robotic systems, humans can operate in an unstructured, dynamic environment due to their ability to anticipate changes before they occur and react accordingly. Recreating this artificially by using a forward model would lead to the huge computational task of simulating a world full of complex nonlinear dynamics and autonomous human agents. In this study, a controller based on anticipating synchronization, where a 'leader' dynamical system is predicted by a coupled 'follower' with delayed self-feedback, is used to modify a robot's dynamical behaviour to follow that of a series of leaky integrators and harmonic oscillators. This allows the robot (follower) to be coupled with a collaborative partner (leader) to anticipate its movements, without a complete model of its behaviour. This is tested by tasking a simulated Baxter robot with performing a collaborative manual coordination task with an autonomous partner under a range of feedback delay conditions, confirming its ability to anticipate using oscillators instead of a detailed forward model.

Neural and physiological data from participants listening to affective music.

Music provides a means of communicating affective meaning. However, the neurological mechanisms by which music induces affect are not fully understood. Our project sought to investigate this through a series of experiments into how humans react to affective musical stimuli and how physiological and neurological signals recorded from those participants change in accordance with self-reported changes in affect. In this paper, the datasets recorded over the course of this project are presented, including details of the musical stimuli, participant reports of their felt changes in affective states as they listened to the music, and concomitant recordings of physiological and neurological activity. We also include non-identifying meta data on our participant populations for purposes of further exploratory analysis. This data provides a large and valuable novel resource for researchers investigating emotion, music, and how they affect our neural and physiological activity.

Modeling the Ongoing Dynamics of Short and Long-Range Temporal Correlations in Broadband EEG During Movement.

Electroencephalogram (EEG) undergoes complex temporal and spectral changes during voluntary movement intention. Characterization of such changes has focused mostly on narrowband spectral processes such as Event-Related Desynchronization (ERD) in the sensorimotor rhythms because EEG is mostly considered as emerging from oscillations of the neuronal populations. However, the changes in the temporal dynamics, especially in the broadband arrhythmic EEG have not been investigated for movement intention detection. The Long-Range Temporal Correlations (LRTC) are ubiquitously present in several neuronal processes, typically requiring longer timescales to detect. In this paper, we study the ongoing changes in the dynamics of long- as well as short-range temporal dependencies in the single trial broadband EEG during movement intention. We obtained LRTC in 2 s windows of broadband EEG and modeled it using the Autoregressive Fractionally Integrated Moving Average (ARFIMA) model which allowed simultaneous modeling of short- and long-range temporal correlations. There were significant (p < 0.05) changes in both broadband long- and short-range temporal correlations during movement intention and execution. We discovered that the broadband LRTC and narrowband ERD are complementary processes providing distinct information about movement because eliminating LRTC from the signal did not affect the ERD and conversely, eliminating ERD from the signal did not affect LRTC. Exploring the possibility of applications in Brain Computer Interfaces (BCI), we used hybrid features with combinations of LRTC, ARFIMA, and ERD to detect movement intention. A significantly higher (p < 0.05) classification accuracy of 88.3 ± 4.2% was obtained using the combination of ARFIMA and ERD features together, which also predicted the earliest movement at 1 s before its onset. The ongoing changes in the long- and short-range temporal correlations in broadband EEG contribute to effectively capturing the motor command generation and can be used to detect movement successfully. These temporal dependencies provide different and additional information about the movement.

Electroencephalography reflects the activity of sub-cortical brain regions during approach-withdrawal behaviour while listening to music.

The ability of music to evoke activity changes in the core brain structures that underlie the experience of emotion suggests that it has the potential to be used in therapies for emotion disorders. A large volume of research has identified a network of sub-cortical brain regions underlying music-induced emotions. Additionally, separate evidence from electroencephalography (EEG) studies suggests that prefrontal asymmetry in the EEG reflects the approach-withdrawal response to music-induced emotion. However, fMRI and EEG measure quite different brain processes and we do not have a detailed understanding of the functional relationships between them in relation to music-induced emotion. We employ a joint EEG - fMRI paradigm to explore how EEG-based neural correlates of the approach-withdrawal response to music reflect activity changes in the sub-cortical emotional response network. The neural correlates examined are asymmetry in the prefrontal EEG, and the degree of disorder in that asymmetry over time, as measured by entropy. Participants' EEG and fMRI were recorded simultaneously while the participants listened to music that had been specifically generated to target the elicitation of a wide range of affective states. While listening to this music, participants also continuously reported their felt affective states. Here we report on co-variations in the dynamics of these self-reports, the EEG, and the sub-cortical brain activity. We find that a set of sub-cortical brain regions in the emotional response network exhibits activity that significantly relates to prefrontal EEG asymmetry. Specifically, EEG in the pre-frontal cortex reflects not only cortical activity, but also changes in activity in the amygdala, posterior temporal cortex, and cerebellum. We also find that, while the magnitude of the asymmetry reflects activity in parts of the limbic and paralimbic systems, the entropy of that asymmetry reflects activity in parts of the autonomic response network such as the auditory cortex. This suggests that asymmetry magnitude reflects affective responses to music, while asymmetry entropy reflects autonomic responses to music. Thus, we demonstrate that it is possible to infer activity in the limbic and paralimbic systems from pre-frontal EEG asymmetry. These results show how EEG can be used to measure and monitor changes in the limbic and paralimbic systems. Specifically, they suggest that EEG asymmetry acts as an indicator of sub-cortical changes in activity induced by music. This shows that EEG may be used as a measure of the effectiveness of music therapy to evoke changes in activity in the sub-cortical emotion response network. This is also the first time that the activity of sub-cortical regions, normally considered "invisible" to EEG, has been shown to be characterisable directly from EEG dynamics measured during music listening.

Anticipation from sensation: using anticipating synchronization to stabilize a system with inherent sensory delay.

We present a novel way of using a dynamical model for predictive tracking control that can adapt to a wide range of delays without parameter update. This is achieved by incorporating the paradigm of anticipating synchronization (AS), where a 'slave' system predicts a 'master' via delayed self-feedback. By treating the delayed output of the plant as one half of a 'sensory' AS coupling, the plant and an internal dynamical model can be synchronized such that the plant consistently leads the target's motion. We use two simulated robotic systems with differing arrangements of the plant and internal model ('parallel' and 'serial') to demonstrate that this form of control adapts to a wide range of delays without requiring the parameters of the controller to be changed.

Publisher Correction: Inferring structural connectivity using Ising couplings in models of neuronal networks.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has not been fixed in the paper.

Exploration of neural correlates of movement intention based on characterisation of temporal dependencies in electroencephalography.

Brain computer interfaces (BCIs) provide a direct communication channel by using brain signals, enabling patients with motor impairments to interact with external devices. Motion intention detection is useful for intuitive movement-based BCI as movement is the fundamental mode of interaction with the environment. The aim of this paper is to investigate the temporal dynamics of brain processes using electroencephalography (EEG) to explore novel neural correlates of motion intention. We investigate the changes in temporal dependencies of the EEG by characterising the decay of autocorrelation during asynchronous voluntary finger tapping movement. The evolution of the autocorrelation function is characterised by its relaxation time, which is used as a robust marker for motion intention. We observed that there was reorganisation of temporal dependencies in EEG during motion intention. The autocorrelation decayed slower during movement intention and faster during the resting state. There was an increase in temporal dependence during movement intention. The relaxation time of the autocorrelation function showed significant (p < 0.05) discrimination between movement and resting state with the mean sensitivity of 78.37 ± 8.83%. The relaxation time provides movement related information that is complementary to the well-known event-related desynchronisation (ERD) by characterising the broad band EEG dynamics which is frequency independent in contrast to ERD. It can also detect motion intention on average 0.51s before the actual movement onset. We have thoroughly compared autocorrelation relaxation time features with ERD in four frequency bands. The relaxation time may therefore, complement the well-known features used in motion-based BCI leading to more robust and intuitive BCI solutions. The results obtained suggest that changes in autocorrelation decay may involve reorganisation of temporal dependencies of brain activity over longer duration during motion intention. This opens the possibilities of investigating further the temporal dynamics of fundamental neural processes underpinning motion intention.

Markov Model-Based Method to Analyse Time-Varying Networks in EEG Task-Related Data.

The dynamic nature of functional brain networks is being increasingly recognized in cognitive neuroscience, and methods to analyse such time-varying networks in EEG/MEG data are required. In this work, we propose a pipeline to characterize time-varying networks in single-subject EEG task-related data and further, evaluate its validity on both simulated and experimental datasets. Pre-processing is done to remove channel-wise and trial-wise differences in activity. Functional networks are estimated from short non-overlapping time windows within each "trial," using a sparse-MVAR (Multi-Variate Auto-Regressive) model. Functional "states" are then identified by partitioning the entire space of functional networks into a small number of groups/symbols via k-means clustering.The multi-trial sequence of symbols is then described by a Markov Model (MM). We show validity of this pipeline on realistic electrode-level simulated EEG data, by demonstrating its ability to discriminate "trials" from two experimental conditions in a range of scenarios. We then apply it to experimental data from two individuals using a Brain-Computer Interface (BCI) via a P300 oddball task. Using just the Markov Model parameters, we obtain statistically significant discrimination between target and non-target trials. The functional networks characterizing each 'state' were also highly similar between the two individuals. This work marks the first application of the Markov Model framework to infer time-varying networks from EEG/MEG data. Due to the pre-processing, results from the pipeline are orthogonal to those from conventional ERP averaging or a typical EEG microstate analysis. The results provide powerful proof-of-concept for a Markov model-based approach to analyzing the data, paving the way for its use to track rapid changes in interaction patterns as a task is being performed. MATLAB code for the entire pipeline has been made available.