Information-Theoretical Analysis of EEG Microstate Sequences in Python

MICROSTATELAB: The EEGLAB Toolbox for Resting-State Microstate Analysis

Microstate analysis is a multivariate method that enables investigations of the temporal dynamics of large-scale neural networks in EEG recordings of human brain activity. To meet the enormously increasing interest in this approach, we provide a thoroughly updated version of the first open source EEGLAB toolbox for the standardized identification, visualization, and quantification of microstates in resting-state EEG data. The toolbox allows scientists to (i) identify individual, mean, and grand mean microstate maps using topographical clustering approaches, (ii) check data quality and detect outlier maps, (iii) visualize, sort, and label individual, mean, and grand mean microstate maps according to published maps, (iv) compare topographical similarities of group and grand mean microstate maps and quantify shared variances, (v) obtain the temporal dynamics of the microstate classes in individual EEGs, (vi) export quantifications of these temporal dynamics of the microstates for statistical tests, and finally, (vii) test for topographical differences between groups and conditions using topographic analysis of variance (TANOVA). Here, we introduce the toolbox in a step-by-step tutorial, using a sample dataset of 34 resting-state EEG recordings that are publicly available to follow along with this tutorial. The goals of this manuscript are (a) to provide a standardized, freely available toolbox for resting-state microstate analysis to the scientific community, (b) to allow researchers to use best practices for microstate analysis by following a step-by-step tutorial, and (c) to improve the methodological standards of microstate research by providing previously unavailable functions and recommendations on critical decisions required in microstate analyses.

Establishment and validation of a bad outcomes prediction model based on EEG and clinical parameters in prolonged disorder of consciousness

Objective

This study aimed to explore the electroencephalogram (EEG) indicators and clinical factors that may lead to poor prognosis in patients with prolonged disorder of consciousness (pDOC), and establish and verify a clinical predictive model based on these factors.

Methods

This study included 134 patients suffering from prolonged disorder of consciousness enrolled in our department of neurosurgery. We collected the data of sex, age, etiology, coma recovery scales (CRS-R) score, complications, blood routine, liver function, coagulation and other laboratory tests, resting EEG data and follow-up after discharge. These patients were divided into two groups: training set (n = 107) and verification set (n = 27). These patients were divided into a training set of 107 and a validation set of 27 for this study. Univariate and multivariate regression analysis were used to determine the factors affecting the poor prognosis of pDOC and to establish nomogram model. We use the receiver operating characteristic (ROC) and calibration curves to quantitatively test the effectiveness of the training set and the verification set. In order to further verify the clinical practical value of the model, we use decision curve analysis (DCA) to evaluate the model.

Result

The results from univariate and multivariate logistic regression analyses suggested that an increased frequency of occurrence microstate A, reduced CRS-R scores at the time of admission, the presence of episodes associated with paroxysmal sympathetic hyperactivity (PSH), and decreased fibrinogen levels all function as independent prognostic factors. These factors were used to construct the nomogram. The training and verification sets had areas under the curve of 0.854 and 0.920, respectively. Calibration curves and DCA demonstrated good model performance and significant clinical benefits in both sets.

Conclusion

This study is based on the use of clinically available and low-cost clinical indicators combined with EEG to construct a highly applicable and accurate model for predicting the adverse prognosis of patients with prolonged disorder of consciousness. It provides an objective and reliable tool for clinicians to evaluate the prognosis of prolonged disorder of consciousness, and helps clinicians to provide personalized clinical care and decision-making for patients with prolonged disorder of consciousness and their families.

Complexity Measures for EEG Microstate Sequences: Concepts and Algorithms

EEG microstate sequence analysis quantifies properties of ongoing brain electrical activity which is known to exhibit complex dynamics across many time scales. In this report we review recent developments in quantifying microstate sequence complexity, we classify these approaches with regard to different complexity concepts, and we evaluate excess entropy as a yet unexplored quantity in microstate research. We determined the quantities entropy rate, excess entropy, Lempel-Ziv complexity (LZC), and Hurst exponents on Potts model data, a discrete statistical mechanics model with a temperature-controlled phase transition. We then applied the same techniques to EEG microstate sequences from wakefulness and non-REM sleep stages and used first-order Markov surrogate data to determine which time scales contributed to the different complexity measures. We demonstrate that entropy rate and LZC measure the Kolmogorov complexity (randomness) of microstate sequences, whereas excess entropy and Hurst exponents describe statistical complexity which attains its maximum at intermediate levels of randomness. We confirmed the equivalence of entropy rate and LZC when the LZ-76 algorithm is used, a result previously reported for neural spike train analysis (Amigó et al., Neural Comput 16:717-736, https://doi.org/10.1162/089976604322860677 , 2004). Surrogate data analyses prove that entropy-based quantities and LZC focus on short-range temporal correlations, whereas Hurst exponents include short and long time scales. Sleep data analysis reveals that deeper sleep stages are accompanied by a decrease in Kolmogorov complexity and an increase in statistical complexity. Microstate jump sequences, where duplicate states have been removed, show higher randomness, lower statistical complexity, and no long-range correlations. Regarding the practical use of these methods, we suggest that LZC can be used as an efficient entropy rate estimator that avoids the estimation of joint entropies, whereas entropy rate estimation via joint entropies has the advantage of providing excess entropy as the second parameter of the same linear fit. We conclude that metrics of statistical complexity are a useful addition to microstate analysis and address a complexity concept that is not yet covered by existing microstate algorithms while being actively explored in other areas of brain research.

Propofol Reversibly Attenuates Short-Range Microstate Ordering and 20 Hz Microstate Oscillations

Microstate sequences summarize the changing voltage patterns measured by electroencephalography, using a clustering approach to reduce the high dimensionality of the underlying data. A common approach is to restrict the pattern matching step to local maxima of the global field power (GFP) and to interpolate the microstate fit in between. In this study, we investigate how the anesthetic propofol affects microstate sequence periodicity and predictability, and how these metrics are changed by interpolation. We performed two frequency analyses on microstate sequences, one based on time-lagged mutual information, the other based on Fourier transform methodology, and quantified the effects of interpolation. Resting-state microstate sequences had a 20 Hz frequency peak related to dominant 10 Hz (alpha) rhythms, and the Fourier approach demonstrated that all five microstate classes followed this frequency. The 20 Hz periodicity was reversibly attenuated under moderate propofol sedation, as shown by mutual information and Fourier analysis. Characteristic microstate frequencies could only be observed in non-interpolated microstate sequences and were masked by smoothing effects of interpolation. Information-theoretic analysis revealed faster microstate dynamics and larger entropy rates under propofol, whereas Shannon entropy did not change significantly. In moderate sedation, active information storage decreased for non-interpolated sequences. Signatures of non-equilibrium dynamics were observed in non-interpolated sequences, but no changes were observed between sedation levels. All changes occurred while subjects were able to perform an auditory perception task. In summary, we show that low dose propofol reversibly increases the randomness of microstate sequences and attenuates microstate oscillations without correlation to cognitive task performance. Microstate dynamics between GFP peaks reflect physiological processes that are not accessible in interpolated sequences.

Frequency Analysis of EEG Microstate Sequences in Wakefulness and NREM Sleep

The majority of EEG microstate analyses concern wakefulness, and the existing sleep studies have focused on changes in spatial microstate properties and on microstate transitions between adjacent time points, the shortest available time scale. We present a more extensive time series analysis of unsmoothed EEG microstate sequences in wakefulness and non-REM sleep stages across many time scales. Very short time scales are assessed with Markov tests, intermediate time scales by the entropy rate and long time scales by a spectral analysis which identifies characteristic microstate frequencies. During the descent from wakefulness to sleep stage N3, we find that the increasing mean microstate duration is a gradual phenomenon explained by a continuous slowing of microstate dynamics as described by the relaxation time of the transition probability matrix. The finite entropy rate, which considers longer microstate histories, shows that microstate sequences become more predictable (less random) with decreasing vigilance level. Accordingly, the Markov property is absent in wakefulness but in sleep stage N3, 10/19 subjects have microstate sequences compatible with a second-order Markov process. A spectral microstate analysis is performed by comparing the time-lagged mutual information coefficients of microstate sequences with the autocorrelation function of the underlying EEG. We find periodic microstate behavior in all vigilance states, linked to alpha frequencies in wakefulness, theta activity in N1, sleep spindle frequencies in N2, and in the delta frequency band in N3. In summary, we show that EEG microstates are a dynamic phenomenon with oscillatory properties that slow down in sleep and are coupled to specific EEG frequencies across several sleep stages.

Mindfulness meditation styles differently modulate source-level MEG microstate dynamics and complexity

Background

The investigation of mindfulness meditation practice, classically divided into focused attention meditation (FAM), and open monitoring meditation (OMM) styles, has seen a long tradition of theoretical, affective, neurophysiological and clinical studies. In particular, the high temporal resolution of magnetoencephalography (MEG) or electroencephalography (EEG) has been exploited to fill the gap between the personal experience of meditation practice and its neural correlates. Mounting evidence, in fact, shows that human brain activity is highly dynamic, transiting between different brain states (microstates). In this study, we aimed at exploring MEG microstates at source-level during FAM, OMM and in the resting state, as well as the complexity and criticality of dynamic transitions between microstates.

Methods

Ten right-handed Theravada Buddhist monks with a meditative expertise of minimum 2,265 h participated in the experiment. MEG data were acquired during a randomized block design task (6 min FAM, 6 min OMM, with each meditative block preceded and followed by 3 min resting state). Source reconstruction was performed using eLORETA on individual cortical space, and then parcellated according to the Human Connect Project atlas. Microstate analysis was then applied to parcel level signals in order to derive microstate topographies and indices. In addition, from microstate sequences, the Hurst exponent and the Lempel-Ziv complexity (LZC) were computed.

Results

Our results show that the coverage and occurrence of specific microstates are modulated either by being in a meditative state or by performing a specific meditation style. Hurst exponent values in both meditation conditions are reduced with respect to the value observed during rest, LZC shows significant differences between OMM, FAM, and REST, with a progressive increase from REST to FAM to OMM.

Discussion

Importantly, we report changes in brain criticality indices during meditation and between meditation styles, in line with a state-like effect of meditation on cognitive performance. In line with previous reports, we suggest that the change in cognitive state experienced in meditation is paralleled by a shift with respect to critical points in brain dynamics.

Machine learning using multimodal clinical, electroencephalographic, and magnetic resonance imaging data can predict incident depression in adults with epilepsy: A pilot study

OBJECTIVE

This study was undertaken to develop a multimodal machine learning (ML) approach for predicting incident depression in adults with epilepsy.

METHODS

We randomly selected 200 patients from the Calgary Comprehensive Epilepsy Program registry and linked their registry-based clinical data to their first-available clinical electroencephalogram (EEG) and magnetic resonance imaging (MRI) study. We excluded patients with a clinical or Neurological Disorders Depression Inventory for Epilepsy (NDDI-E)-based diagnosis of major depression at baseline. The NDDI-E was used to detect incident depression over a median of 2.4 years of follow-up (interquartile range [IQR] = 1.5-3.3 years). A ReliefF algorithm was applied to clinical as well as quantitative EEG and MRI parameters for feature selection. Six ML algorithms were trained and tested using stratified threefold cross-validation. Multiple metrics were used to assess model performances.

RESULTS

Of 200 patients, 150 had EEG and MRI data of sufficient quality for ML, of whom 59 were excluded due to prevalent depression. Therefore, 91 patients (41 women) were included, with a median age of 29 (IQR = 22-44) years. A total of 42 features were selected by ReliefF, none of which was a quantitative MRI or EEG variable. All models had a sensitivity > 80%, and five of six had an F1 score ≥ .72. A multilayer perceptron model had the highest F1 score (median = .74, IQR = .71-.78) and sensitivity (84.3%). Median area under the receiver operating characteristic curve and normalized Matthews correlation coefficient were .70 (IQR = .64-.78) and .57 (IQR = .50-.65), respectively.

SIGNIFICANCE

Multimodal ML using baseline features can predict incident depression in this population. Our pilot models demonstrated high accuracy for depression prediction. However, overall performance and calibration can be improved. This model has promise for identifying those at risk for incident depression during follow-up, although efforts to refine it in larger populations along with external validation are required.

EEG microstate changes during hyperbaric oxygen therapy in patients with chronic disorders of consciousness

Hyperbaric oxygen (HBO) therapy is an effective treatment for patients with disorders of consciousness (DOC). In this study, real-time electroencephalogram (EEG) recordings were obtained from patients with DOC during HBO therapy. EEG microstate indicators including mean microstate duration (MMD), ratio of total time covered (RTT), global explained variance (GEV), transition probability, mean occurrence, and mean global field power (GFP) were compared before and during HBO therapy. The results showed that the duration of microstate C in all patients with DOC increased after 20 min of HBO therapy (p < 0.05). Further statistical analysis found that the duration of microstate C was longer in the higher CRS-R group (≥8, 17 cases) than in the lower group (<8, 24 cases) during HBO treatment. In the higher CRS-R group, the transition probabilities from microstate A to microstate C and from microstate C to microstate A also increased significantly compared with the probability before treatment (p < 0.05). Microstate C is generally considered to be related to a salience network; an increase in the transition probability between microstate A and microstate C indicates increased information exchange between the auditory network and the salience network. The results of this study show that HBO therapy has a specific activating effect on attention and cognitive control in patients and causes increased activity in the primary sensory cortex (temporal lobe and occipital lobe). This study demonstrates that real-time EEG detection and analysis during HBO is a clinically feasible method for assessing brain function in patients with DOC. During HBO therapy, some EEG microstate indicators show significant changes related to the state of consciousness in patients with chronic DOC. This will be complementary to important electrophysiological indicators for assessing consciousness and may also provide an objective foundation for the precise treatment of patients with DOC.

Resting-state EEG functional connectivity predicts post-traumatic stress disorder subtypes in veterans

Objective. Post-traumatic stress disorder (PTSD) is highly heterogeneous, and identification of quantifiable biomarkers that could pave the way for targeted treatment remains a challenge. Most previous electroencephalography (EEG) studies on PTSD have been limited to specific handpicked features, and their findings have been highly variable and inconsistent. Therefore, to disentangle the role of promising EEG biomarkers, we developed a machine learning framework to investigate a wide range of commonly used EEG biomarkers in order to identify which features or combinations of features are capable of characterizing PTSD and potential subtypes.Approach. We recorded 5 min of eyes-closed and 5 min of eyes-open resting-state EEG from 202 combat-exposed veterans (53% with probable PTSD and 47% combat-exposed controls). Multiple spectral, temporal, and connectivity features were computed and logistic regression, random forest, and support vector machines with feature selection methods were employed to classify PTSD. To obtain robust results, we performed repeated two-layer cross-validation to test on an entirely unseen test set.Main results. Our classifiers obtained a balanced test accuracy of up to 62.9% for predicting PTSD patients. In addition, we identified two subtypes within PTSD: one where EEG patterns were similar to those of the combat-exposed controls, and another that were characterized by increased global functional connectivity. Our classifier obtained a balanced test accuracy of 79.4% when classifying this PTSD subtype from controls, a clear improvement compared to predicting the whole PTSD group. Interestingly, alpha connectivity in the dorsal and ventral attention network was particularly important for the prediction, and these connections were positively correlated with arousal symptom scores, a central symptom cluster of PTSD.Significance. Taken together, the novel framework presented here demonstrates how unsupervised subtyping can delineate heterogeneity and improve machine learning prediction of PTSD, and may pave the way for better identification of quantifiable biomarkers.

EEG microstate features according to performance on a mental arithmetic task

In this study, we hypothesized that task performance could be evaluated applying EEG microstate to mental arithmetic task. This pilot study also aimed at evaluating the efficacy of microstates as novel features to discriminate task performance. Thirty-six subjects were divided into good and poor performers, depending on how well they performed the task. Microstate features were derived from EEG recordings during resting and task states. In the good performers, there was a decrease in type C and an increase in type D features during the task compared to the resting state. Mean duration and occurrence decreased and increased, respectively. In the poor performers, occurrence of type D feature, mean duration and occurrence showed greater changes. We investigated whether microstate features were suitable for task performance classification and eleven features including four archetypes were selected by recursive feature elimination (RFE). The model that implemented them showed the highest classification performance for differentiating between groups. Our pilot findings showed that the highest mean Area Under Curve (AUC) was 0.831. This study is the first to apply EEG microstate features to specific cognitive tasks in healthy subjects, suggesting that EEG microstate features can reflect task achievement.

EEG microstate periodicity explained by rotating phase patterns of resting-state alpha oscillations

Spatio-temporal patterns in electroencephalography (EEG) can be described by microstate analysis, a discrete approximation of the continuous electric field patterns produced by the cerebral cortex. Resting-state EEG microstates are largely determined by alpha frequencies (8-12 Hz) and we recently demonstrated that microstates occur periodically with twice the alpha frequency. To understand the origin of microstate periodicity, we analyzed the analytic amplitude and the analytic phase of resting-state alpha oscillations independently. In continuous EEG data we found rotating phase patterns organized around a small number of phase singularities which varied in number and location. The spatial rotation of phase patterns occurred with the underlying alpha frequency. Phase rotors coincided with periodic microstate motifs involving the four canonical microstate maps. The analytic amplitude showed no oscillatory behaviour and was almost static across time intervals of 1-2 alpha cycles, resulting in the global pattern of a standing wave. In n=23 healthy adults, time-lagged mutual information analysis of microstate sequences derived from amplitude and phase signals of awake eyes-closed EEG records showed that only the phase component contributed to the periodicity of microstate sequences. Phase sequences showed mutual information peaks at multiples of 50 ms and the group average had a main peak at 100 ms (10 Hz), whereas amplitude sequences had a slow and monotonous information decay. This result was confirmed by an independent approach combining temporal principal component analysis (tPCA) and autocorrelation analysis. We reproduced our observations in a generic model of EEG oscillations composed of coupled non-linear oscillators (Stuart-Landau model). Phase-amplitude dynamics similar to experimental EEG occurred when the oscillators underwent a supercritical Hopf bifurcation, a common feature of many computational models of the alpha rhythm. These findings explain our previous description of periodic microstate recurrence and its relation to the time scale of alpha oscillations. Moreover, our results corroborate the predictions of computational models and connect experimentally observed EEG patterns to properties of critical oscillator networks.

Partial Autoinformation to Characterize Symbolic Sequences

An information-theoretic approach to numerically determine the Markov order of discrete stochastic processes defined over a finite state space is introduced. To measure statistical dependencies between different time points of symbolic time series, two information-theoretic measures are proposed. The first measure is time-lagged mutual information between the random variables X_n and X_n+k, representing the values of the process at time points n and n + k, respectively. The measure will be termed autoinformation, in analogy to the autocorrelation function for metric time series, but using Shannon entropy rather than linear correlation. This measure is complemented by the conditional mutual information between X_n and X_n+k, removing the influence of the intermediate values X_n+k−1, …, X_n+1. The second measure is termed partial autoinformation, in analogy to the partial autocorrelation function (PACF) in metric time series analysis. Mathematical relations with known quantities such as the entropy rate and active information storage are established. Both measures are applied to a number of examples, ranging from theoretical Markov and non-Markov processes with known stochastic properties, to models from statistical physics, and finally, to a discrete transform of an EEG data set. The combination of autoinformation and partial autoinformation yields important insights into the temporal structure of the data in all test cases. For first- and higher-order Markov processes, partial autoinformation correctly identifies the order parameter, but also suggests extended, non-Markovian effects in the examples that lack the Markov property. For three hidden Markov models (HMMs), the underlying Markov order is found. The combination of both quantities may be used as an early step in the analysis of experimental, non-metric time series and can be employed to discover higher-order Markov dependencies, non-Markovianity and periodicities in symbolic time series.

EEG Microstate Sequences From Different Clustering Algorithms Are Information-Theoretically Invariant

We analyse statistical and information-theoretical properties of EEG microstate sequences, as seen through the lens of five different clustering algorithms. Microstate sequences are computed for n = 20 resting state EEG recordings during wakeful rest. The input for all clustering algorithms is the set of EEG topographic maps obtained at local maxima of the spatial variance. This data set is processed by two classical microstate clustering algorithms (1) atomize and agglomerate hierarchical clustering (AAHC) and (2) a modified K-means algorithm, as well as by (3) K-medoids, (4) principal component analysis (PCA) and (5) fast independent component analysis (Fast-ICA). Using this technique, EEG topographies can be substituted with microstate labels by competitive fitting based on spatial correlation, resulting in a symbolic, non-metric time series, the microstate sequence. Microstate topographies and symbolic time series are further analyzed statistically, including static and dynamic properties. Static properties, which do not contain information about temporal dependencies of the microstate sequence include the maximum similarity of microstate maps within and between the tested clustering algorithms, the global explained variance and the Shannon entropy of the microstate sequences. Dynamic properties are sensitive to temporal correlations between the symbols and include the mixing time of the microstate transition matrix, the entropy rate of the microstate sequences and the location of the first local maximum of the autoinformation function. We also test the Markov property of microstate sequences, the time stationarity of the transition matrix and detect periodicities by means of time-lagged mutual information. Finally, possible long-range correlations of microstate sequences are assessed via Hurst exponent estimation. We find that while static properties partially reflect properties of the clustering algorithms, information-theoretical quantities are largely invariant with respect to the clustering method used. As each clustering algorithm has its own profile of computational speed, ease of implementation, determinism vs. stochasticity and theoretical underpinnings, our results convey a positive message concerning the free choice of method and the comparability of results obtained from different algorithms. The invariance of these quantities implies that the tested properties are algorithm-independent, inherent features of resting state EEG derived microstate sequences.