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

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

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 Xn and Xn+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 Xn and Xn+k, removing the influence of the intermediate values Xn+k−1, …, Xn+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

On the Reliability of the EEG Microstate Approach.

EEG microstates represent functional brain networks observable in resting EEG recordings that remain stable for 40-120ms before rapidly switching into another network. It is assumed that microstate characteristics (i.e., durations, occurrences, percentage coverage, and transitions) may serve as neural markers of mental and neurological disorders and psychosocial traits. However, robust data on their retest-reliability are needed to provide the basis for this assumption. Furthermore, researchers currently use different methodological approaches that need to be compared regarding their consistency and suitability to produce reliable results. Based on an extensive dataset largely representative of western societies (2 days with two resting EEG measures each; day one: n = 583; day two: n = 542) we found good to excellent short-term retest-reliability of microstate durations, occurrences, and coverages (average ICCs = 0.874-0.920). There was good overall long-term retest-reliability of these microstate characteristics (average ICCs = 0.671-0.852), even when the interval between measures was longer than half a year, supporting the longstanding notion that microstate durations, occurrences, and coverages represent stable neural traits. Findings were robust across different EEG systems (64 vs. 30 electrodes), recording lengths (3 vs. 2 min), and cognitive states (before vs. after experiment). However, we found poor retest-reliability of transitions. There was good to excellent consistency of microstate characteristics across clustering procedures (except for transitions), and both procedures produced reliable results. Grand-mean fitting yielded more reliable results compared to individual fitting. Overall, these findings provide robust evidence for the reliability of the microstate approach

Degradation of EEG microstates patterns in subjective cognitive decline and mild cognitive impairment: Early biomarkers along the Alzheimer's Disease continuum?

Alzheimer's disease (AD) pathological changes may begin up to decades earlier than the appearance of the first symptoms of cognitive decline. Subjective cognitive decline (SCD) could be the first pre-clinical sign of possible AD, which might be followed by mild cognitive impairment (MCI), the initial stage of clinical cognitive decline. However, the neural correlates of these prodromic stages are not completely clear yet. Recent studies suggest that EEG analysis tools characterizing the cortical activity as a whole, such as microstates and cortical regions connectivity, might support a characterization of SCD and MCI conditions. Here we test this approach by performing a broad set of analyses to identify the prominent EEG markers differentiating SCD (n = 57), MCI (n = 46) and healthy control subjects (HC, n = 19). We found that the salient differences were in the temporal structure of the microstates patterns, with MCI being associated with less complex sequences due to the altered transition probability, frequency and duration of canonic microstate C. Spectral content of EEG, network connectivity, and spatial arrangement of microstates were instead largely similar in the three groups. Interestingly, comparing properties of EEG microstates in different cerebrospinal fluid (CSF) biomarkers profiles, we found that canonic microstate C displayed significant differences in topography in AD-like profile. These results show that the progression of dementia might be associated with a degradation of the cortical organization captured by microstates analysis, and that this leads to altered transitions between cortical states. Overall, our approach paves the way for the use of non-invasive EEG recordings in the identification of possible biomarkers of progression to AD from its prodromal states

Complexity of multi-dimensional spontaneous EEG decreases during propofol induced general anaesthesia

Emerging neural theories of consciousness suggest a correlation between a specific type of neural dynamical complexity and the level of consciousness: When awake and aware, causal interactions between brain regions are both integrated (all regions are to a certain extent connected) and differentiated (there is inhomogeneity and variety in the interactions). In support of this, recent work by Casali et al (2013) has shown that Lempel-Ziv complexity correlates strongly with conscious level, when computed on the EEG response to transcranial magnetic stimulation. Here we investigated complexity of spontaneous high-density EEG data during propofol-induced general anaesthesia. We consider three distinct measures: (i) Lempel-Ziv complexity, which is derived from how compressible the data are; (ii) amplitude coalition entropy, which measures the variability in the constitution of the set of active channels; and (iii) the novel synchrony coalition entropy (SCE), which measures the variability in the constitution of the set of synchronous channels. After some simulations on Kuramoto oscillator models which demonstrate that these measures capture distinct ‘flavours’ of complexity, we show that there is a robustly measurable decrease in the complexity of spontaneous EEG during general anaesthesia

Scale-Free and Multifractal Time Dynamics of fMRI Signals during Rest and Task

Scaling temporal dynamics in functional MRI (fMRI) signals have been evidenced for a decade as intrinsic characteristics of ongoing brain activity (Zarahn et al., 1997). Recently, scaling properties were shown to fluctuate across brain networks and to be modulated between rest and task (He, 2011): notably, Hurst exponent, quantifying long memory, decreases under task in activating and deactivating brain regions. In most cases, such results were obtained: First, from univariate (voxelwise or regionwise) analysis, hence focusing on specific cognitive systems such as Resting-State Networks (RSNs) and raising the issue of the specificity of this scale-free dynamics modulation in RSNs. Second, using analysis tools designed to measure a single scaling exponent related to the second order statistics of the data, thus relying on models that either implicitly or explicitly assume Gaussianity and (asymptotic) self-similarity, while fMRI signals may significantly depart from those either of those two assumptions (Ciuciu et al., 2008; Wink et al., 2008). To address these issues, the present contribution elaborates on the analysis of the scaling properties of fMRI temporal dynamics by proposing two significant variations. First, scaling properties are technically investigated using the recently introduced Wavelet Leader-based Multifractal formalism (WLMF; Wendt et al., 2007). This measures a collection of scaling exponents, thus enables a richer and more versatile description of scale invariance (beyond correlation and Gaussianity), referred to as multifractality. Also, it benefits from improved estimation performance compared to tools previously used in the literature. Second, scaling properties are investigated in both RSN and non-RSN structures (e.g., artifacts), at a broader spatial scale than the voxel one, using a multivariate approach, namely the Multi-Subject Dictionary Learning (MSDL) algorithm (Varoquaux et al., 2011) that produces a set of spatial components that appear more sparse than their Independent Component Analysis (ICA) counterpart. These tools are combined and applied to a fMRI dataset comprising 12 subjects with resting-state and activation runs (Sadaghiani et al., 2009). Results stemming from those analysis confirm the already reported task-related decrease of long memory in functional networks, but also show that it occurs in artifacts, thus making this feature not specific to functional networks. Further, results indicate that most fMRI signals appear multifractal at rest except in non-cortical regions. Task-related modulation of multifractality appears only significant in functional networks and thus can be considered as the key property disentangling functional networks from artifacts. These finding are discussed in the light of the recent literature reporting scaling dynamics of EEG microstate sequences at rest and addressing non-stationarity issues in temporally independent fMRI modes

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

Multivariate Hurst Exponent Estimation in FMRI. Application to Brain Decoding of Perceptual Learning

International audienceSo far considered as noise in neuroscience, irregular arrhyth-mic field potential activity accounts for the majority of the signal power recorded in EEG or MEG [1, 2]. This brain activity follows a power law spectrum P (f) ∼ 1/f β in the limit of low frequencies, which is a hallmark of scale invariance. Recently, several studies [1, 3–6] have shown that the slope β (or equivalently Hurst exponent H) tends to be modulated by task performance or cognitive state (eg, sleep vs awake). These observations were confirmed in fMRI [7–9] although the short length of fMRI time series makes these findings less reliable. In this paper, to compensate for the slower sampling rate in fMRI, we extend univariate wavelet-based Hurst exponent estimator to a multivariate setting using spatial regular-ization. Next, we demonstrate the relevance of the proposed tools on resting-state fMRI data recorded in three groups of individuals once they were specifically trained to a visual discrimination task during a MEG experiment [10]. In a supervised classification framework, our multivariate approach permits to better predict the type of training the participants received as compared to their univariate counterpart