On the Quantization of Time-Varying Phase Synchrony Patterns into Distinct Functional Connectivity Microstates (FCμstates) in a Multi-trial Visual ERP Paradigm

Atypical instantaneous spatio-temporal patterns of neural dynamics in Alzheimer's disease

Cognitive functions produced by large-scale neural integrations are the most representative 'emergence phenomena' in complex systems. A novel approach focusing on the instantaneous phase difference of brain oscillations across brain regions has succeeded in detecting moment-to-moment dynamic functional connectivity. However, it is restricted to pairwise observations of two brain regions, contrary to large-scale spatial neural integration in the whole-brain. In this study, we introduce a microstate analysis to capture whole-brain instantaneous phase distributions instead of pairwise differences. Upon applying this method to electroencephalography signals of Alzheimer's disease (AD), which is characterised by progressive cognitive decline, the AD-specific state transition among the four states defined as the leading phase location due to the loss of brain regional interactions could be promptly characterised. In conclusion, our synthetic analysis approach, focusing on the microstate and instantaneous phase, enables the capture of the instantaneous spatiotemporal neural dynamics of brain activity and characterises its pathological conditions.

Temporal segmentation of EEG based on functional connectivity network structure

In the study of brain functional connectivity networks, it is assumed that a network is built from a data window in which activity is stationary. However, brain activity is non-stationary over sufficiently large time periods. Addressing the analysis electroencephalograph (EEG) data, we propose a data segmentation method based on functional connectivity network structure. The goal of segmentation is to ensure that within a window of analysis, there is similar network structure. We designed an intuitive and flexible graph distance measure to quantify the difference in network structure between two analysis windows. This measure is modular: a variety of node importance indices can be plugged into it. We use a reference window versus sliding window comparison approach to detect changes, as indicated by outliers in the distribution of graph distance values. Performance of our segmentation method was tested in simulated EEG data and real EEG data from a drone piloting experiment (using correlation or phase-locking value as the functional connectivity strength metric). We compared our method under various node importance measures and against matrix-based dissimilarity metrics that use singular value decomposition on the connectivity matrix. The results show the graph distance approach worked better than matrix-based approaches; graph distance based on partial node centrality was most sensitive to network structural changes, especially when connectivity matrix values change little. The proposed method provides EEG data segmentation tailored for detecting changes in terms of functional connectivity networks. Our study provides a new perspective on EEG segmentation, one that is based on functional connectivity network structure differences.

On the characterization of cognitive tasks using activity-specific short-lived synchronization between electroencephalography channels

Accurate characterization of brain activity during a cognitive task is challenging due to the dynamically changing and the complex nature of the brain. The majority of the proposed approaches assume stationarity in brain activity and disregard the systematic timing organization among brain regions during cognitive tasks. In this study, we propose a novel cognitive activity recognition method that captures the activity-specific timing parameters from training data that elicits maximal average short-lived pairwise synchronization between electroencephalography signals. We evaluated the characterization power of the activity-specific timing parameter triplets in a motor imagery activity recognition framework. The activity-specific timing parameter triplets consist of latency of the maximally synchronized signal segments from activity onset Δt, the time lag between maximally synchronized signal segments τ, and the duration of the maximally synchronized signal segments w. We used cosine-based similarity, wavelet bi-coherence, phase-locking value, phase coherence value, linearized mutual information, and cross-correntropy to calculate the channel synchronizations at the specific timing parameters. Recognition performances as well as statistical analyses on both BCI Competition-III dataset IVa and PhysioNet Motor Movement/Imagery dataset, indicate that the inter-channel short-lived synchronization calculated using activity-specific timing parameter triplets elicit significantly distinct synchronization profiles for different motor imagery tasks and can thus reliably be used for cognitive task recognition purposes.

Dyconnmap: Dynamic connectome mapping-A neuroimaging python module

Despite recent progress in the analysis of neuroimaging data sets, our comprehension of the main mechanisms and principles which govern human brain cognition and function remains incomplete. Network neuroscience makes substantial efforts to manipulate these challenges and provide real answers. For the last decade, researchers have been modelling brain structure and function via a graph or network that comprises brain regions that are either anatomically connected via tracts or functionally via a more extensive repertoire of functional associations. Network neuroscience is a relatively new multidisciplinary scientific avenue of the study of complex systems by pursuing novel ways to analyze, map, store and model the essential elements and their interactions in complex neurobiological systems, particularly the human brain, the most complex system in nature. Due to a rapid expansion of neuroimaging data sets' size and complexity, it is essential to propose and adopt new empirical tools to track dynamic patterns between neurons and brain areas and create comprehensive maps. In recent years, there is a rapid growth of scientific interest in moving functional neuroimaging analysis beyond simplified group or time‐averaged approaches and sophisticated algorithms that can capture the time‐varying properties of functional connectivity. We describe algorithms and network metrics that can capture the dynamic evolution of functional connectivity under this perspective. We adopt the word ‘chronnectome’ (integration of the Greek word ‘Chronos’, which means time, and connectome) to describe this specific branch of network neuroscience that explores how mutually informed brain activity correlates across time and brain space in a functional way. We also describe how good temporal mining of temporally evolved dynamic functional networks could give rise to the detection of specific brain states over which our brain evolved. This characteristic supports our complex human mind. The temporal evolution of these brain states and well‐known network metrics could give rise to new analytic trends. Functional brain networks could also increase the multi‐faced nature of the dynamic networks revealing complementary information. Finally, we describe a python module (https://github.com/makism/dyconnmap) which accompanies this article and contains a collection of dynamic complex network analytics and measures and demonstrates its great promise for the study of a healthy subject's repeated fMRI scans.

Reconfiguration of αmplitude driven dominant coupling modes (DoCM) mediated by α-band in adolescents with schizophrenia spectrum disorders

Electroencephalography (EEG) based biomarkers have been shown to correlate with the presence of psychotic disorders. Increased delta and decreased alpha power in psychosis indicate an abnormal arousal state. We investigated brain activity across the basic EEG frequencies and also dynamic functional connectivity of both intra and cross-frequency coupling that could reveal a neurophysiological biomarker linked to an aberrant modulating role of alpha frequency in adolescents with schizophrenia spectrum disorders (SSDs). A dynamic functional connectivity graph (DFCG) has been estimated using the imaginary part of phase lag value (iPLV) and correlation of the envelope (corrEnv). We analyzed DFCG profiles of electroencephalographic resting state (eyes closed) recordings of healthy controls (HC) (n = 39) and SSDs subjects (n = 45) in basic frequency bands {δ,θ,α₁,α₂,β₁,β₂,γ}. In our analysis, we incorporated both intra and cross-frequency coupling modes. Adopting our recent Dominant Coupling Mode (DοCM) model leads to the construction of an integrated DFCG (iDFCG) that encapsulates the functional strength and the DοCM of every pair of brain areas. We revealed significantly higher ratios of delta/alpha₁,₂ power spectrum in SSDs subjects versus HC. The probability distribution (PD) of amplitude driven DoCM mediated by alpha frequency differentiated SSDs from HC with absolute accuracy (100%). The network Flexibility Index (FI) was significantly lower for subjects with SSDs compared to the HC group. Our analysis supports the central role of alpha frequency alterations in the neurophysiological mechanisms of SSDs. Currents findings open up new diagnostic pathways to clinical detection of SSDs and support the design of rational neurofeedback training.

Applications of Convolutional Neural Networks in neurodegeneration and physiological aging

The process of aging is linked with significant changes in a human's physiological organization and structure. This is more evident in the case of the brain whose functions generally vary between young and old individuals. Detecting such patterns can be of significant importance especially during the Mild Cognitive Impairment (MCI) stage which is a transition state before the clinical onset of dementia. Intervening in that stage may delay or eventually prevent dementia onset. In this paper we propose a new methodology based in electroencephalographic (EEG) recordings, aiming to classify individuals into healthy, pathological (patients diagnosed with MCI or Mild Dementia) and young, old groups (healthy individuals over and under 50 years of age) through functional connectivity and macro-architecture features. These features are calculated on the estimated brain region activations through the inverse problem solution, enabling us to transform the sensor level EEG recordings through an appropriate transformation matrix. Afterwards, Synchronization Likelihood and Relative Wavelet Entropy values were calculated along with the graph metrics corresponding to the functional connectivity values, as well as the relative energy contributions of five EEG bands (delta, theta, alpha, beta and gamma). These features were organized in Red, Green, Blue (RGB) image-like data structures. Therefore, it was possible to classify each individual into one of the two groups per experiment employing Convolutional Neural Networks. From the maximum classification accuracy achieved on the test set, 90.48% for the pathological aging group and 91.19% for the physiological aging, it is evident that the proposed approach is capable of providing adequate health and age group classification.

EEG source-space synchrostate transitions and Markov modeling in the math-gifted brain during a long-chain reasoning task

To reveal transition dynamics of global neuronal networks of math-gifted adolescents in handling long-chain reasoning, this study explores momentary phase-synchronized patterns, that is, electroencephalogram (EEG) synchrostates, of intracerebral sources sustained in successive 50 ms time windows during a reasoning task and non-task idle process. Through agglomerative hierarchical clustering for functional connectivity graphs and nested iterative cosine similarity tests, this study identifies seven general and one reasoning-specific prototypical functional connectivity patterns from all synchrostates. Markov modeling is performed for the time-sequential synchrostates of each trial to characterize the interstate transitions. The analysis reveals that default mode network, central executive network (CEN), dorsal attention network, cingulo-opercular network, left/right ventral frontoparietal network, and ventral visual network aperiodically recur over non-task or reasoning process, exhibiting high predictability in interactively reachable transitions. Compared to non-gifted subjects, math-gifted adolescents show higher fractional occupancy and mean duration in CEN and reasoning-triggered transient right frontotemporal network (rFTN) in the time course of the reasoning process. Statistical modeling of Markov chains reveals that there are more self-loops in CEN and rFTN of the math-gifted brain, suggesting robust state durability in temporally maintaining the topological structures. Besides, math-gifted subjects show higher probabilities in switching from the other types of synchrostates to CEN and rFTN, which represents more adaptive reconfiguration of connectivity pattern in the large-scale cortical network for focused task-related information processing, which underlies superior executive functions in controlling goal-directed persistence and high predictability of implementing imagination and creative thinking during long-chain reasoning.

The Musical Structure of Time in the Brain: Repetition, Rhythm, and Harmony in fMRI During Rest and Passive Movie Viewing

Space generally overshadows time in the construction of theories in cognitive neuroscience. In this paper, we pivot from the spatial axes to the temporal, analyzing fMRI image series to reveal structures in time rather than space. To determine affinities among global brain patterns at different times, core concepts in network analysis (derived from graph theory) were applied temporally, as relations among brain images at every time point during an fMRI scanning epoch. To explore the temporal structures observed through this adaptation of network analysis, data from 180 subjects in the Human Connectome Project were examined, during two experimental conditions: passive movie viewing and rest. The temporal brain, like the spatial brain, exhibits a modular structure, where "modules" are intermittent (distributed in time). These temporal entities are here referred to as themes. Short sequences of themes - motifs - were studied in sequences from 4 to 11 s in length. Many motifs repeated at constant intervals, and are therefore rhythmic; rhythms, converted to frequencies, were often harmonic. We speculate that the structure and interaction of these global oscillations underwrites the capacity to experience and navigate a world which is both recognizably stable and noticeably changing at every moment - a temporal world. In its temporal structure, this brain-constituted world resembles music.

Aberrant Whole-Brain Transitions and Dynamics of Spontaneous Network Microstates in Mild Traumatic Brain Injury

Dynamic Functional Connectivity (DFC) analysis is a promising approach for the characterization of brain electrophysiological activity. In this study, we investigated abnormal alterations due to mild Traumatic Brain Injury (mTBI) using DFC of the source reconstructed magnetoencephalographic (MEG) resting-state recordings. Brain activity in several well-known frequency bands was first reconstructed using beamforming of the MEG data to determine ninety anatomical brain regions of interest. A DFC graph was formulated using the imaginary part of phase-locking values, which were obtained from 30 mTBI patients and 50 healthy controls (HC). Subsequently, we estimated normalized Laplacian transformations of individual, statistically and topologically filtered quasi-static graphs. The corresponding eigenvalues of each node synchronization were then computed and through the neural-gas algorithm, we quantized the evolution of the eigenvalues resulting in distinct network microstates (NMstates). The discrimination level between the two groups was assessed using an iterative cross-validation classification scheme with features either the NMstates in each frequency band, or the combination of the so-called chronnectomics (flexibility index, occupancy time of NMstate, and Dwell time) with the complexity index over the evolution of the NMstates across all frequency bands. Classification performance based on chronnectomics showed 80% accuracy, 99% sensitivity, and 49% specificity. However, performance was much higher (accuracy: 91-97%, sensitivity: 100%, and specificity: 77-93%) when focusing on the microstates. Exploring the mean node degree within and between brain anatomical networks (default mode network, frontoparietal, occipital, cingulo-opercular, and sensorimotor), a reduced pattern occurred from lower to higher frequency bands, with statistically significant stronger degrees for the HC than the mTBI group. A higher entropic profile on the temporal evolution of the modularity index was observed for both NMstates for the mTBI group across frequencies. A significant difference in the flexibility index was observed between the two groups for the β frequency band. The latter finding may support a central role of the thalamus impairment in mTBI. The current study considers a complete set of frequency-dependent connectomic markers of mTBI-caused alterations in brain connectivity that potentially could serve as markers to assess the return of an injured subject back to normality.

Complexity of brain activity and connectivity in functional neuroimaging

Understanding the complexity of human brain dynamics and brain connectivity across the repertoire of functional neuroimaging and various conditions, is of paramount importance. Novel measures should be designed tailored to the input focusing on multichannel activity and dynamic functional brain connectivity (DFBC). Here, we defined a novel complexity index (CI) from the field of symbolic dynamics that quantifies patterns of different words up to a length from a symbolic sequence. The CI characterizes the complexity of the brain activity. We analysed DFBC by adopting the sliding window approach using imaginary part of phase locking value (iPLV) for EEG/ECoG/MEG and wavelet coherence (WC) for fMRI. Both intra and cross-frequency couplings (CFC) namely phase-to-amplitude were estimated using iPLV/WC at every snapshot of the DFBC. Using proper surrogate analysis, we defined the dominant intrinsic coupling mode (DICM) per pair of regions-of-interest (ROI). The spatiotemporal probability distribution of DICM were reported to reveal the most prominent coupling modes per condition and modality. Finally, a novel flexibility index is defined that quantifies the transition of DICM per pair of ROIs between consecutive time windows. The whole methodology was demonstrated using four neuroimaging datasets (EEG/ECoG/MEG/fMRI). Finally, we succeeded to totally discriminate healthy controls from schizophrenic using FI and dynamic reconfiguration of DICM. Anaesthesia independently of the drug caused a global decreased of complexity in all frequency bands with the exception in δ and alters the dynamic reconfiguration of DICM. CI and DICM of MEG/fMRI resting-state recordings in two spatial scales were high reliable.

Modeling the Switching Behavior of Functional Connectivity Microstates (FCμstates) as a Novel Biomarker for Mild Cognitive Impairment

The need for designing and validating novel biomarkers for the detection of mild cognitive impairment (MCI) is evident. MCI patients have a high risk of developing Alzheimer's disease (AD), and for that reason the introduction of novel and reliable biomarkers is of significant clinical importance. Motivated by recent findings on the rich information of dynamic functional connectivity graphs (DFCGs) about brain (dys) function, we introduced a novel approach of identifying MCI based on magnetoencephalographic (MEG) resting state recordings. The activity of different brain rhythms {δ, 𝜃, α1, α2, β1, β2, γ1, γ2} was first beamformed with linear constrained minimum norm variance in the MEG data to determine 90 anatomical regions of interest (ROIs). A DFCG was then estimated using the imaginary part of phase lag value (iPLV) for both intra-frequency coupling (8) and cross-frequency coupling pairs (28). We analyzed DFCG profiles of neuromagnetic resting state recordings of 18 MCI patients and 22 healthy controls. We followed our model of identifying the dominant intrinsic coupling mode (DICM) across MEG sources and temporal segments, which further leads to the construction of an integrated DFCG (iDFCG). We then filtered statistically and topologically every snapshot of the iDFCG with data-driven approaches. An estimation of the normalized Laplacian transformation for every temporal segment of the iDFCG and the related eigenvalues created a 2D map based on the network metric time series of the eigenvalues (NMTSeigs). The NMTSeigs preserves the non-stationarity of the fluctuated synchronizability of iDCFG for each subject. Employing the initial set of 20 healthy elders and 20 MCI patients, as training set, we built an overcomplete dictionary set of network microstates (n μstates). Afterward, we tested the whole procedure in an extra blind set of 20 subjects for external validation. We succeeded in gaining a high classification accuracy on the blind dataset (85%), which further supports the proposed Markovian modeling of the evolution of brain states. The adaptation of appropriate neuroinformatic tools that combine advanced signal processing and network neuroscience tools could properly manipulate the non-stationarity of time-resolved FC patterns revealing a robust biomarker for MCI.

Exploiting the heightened phase synchrony in patients with neuromuscular disease for the establishment of efficient motor imagery BCIs

Background

Phase synchrony has extensively been studied for understanding neural coordination in health and disease. There are a few studies concerning the implications in the context of BCIs, but its potential for establishing a communication channel in patients suffering from neuromuscular disorders remains totally unexplored. We investigate, here, this possibility by estimating the time-resolved phase connectivity patterns induced during a motor imagery (MI) task and adopting a supervised learning scheme to recover the subject’s intention from the streaming data.

Methods

Electroencephalographic activity from six patients suffering from neuromuscular disease (NMD) and six healthy individuals was recorded during two randomly alternating, externally cued, MI tasks (clenching either left or right fist) and a rest condition. The metric of Phase locking value (PLV) was used to describe the functional coupling between all recording sites. The functional connectivity patterns and the associate network organization was first compared between the two cohorts. Next, working at the level of individual patients, we trained support vector machines (SVMs) to discriminate between “left” and “right” based on different instantiations of connectivity patterns (depending on the encountered brain rhythm and the temporal interval). Finally, we designed and realized a novel brain decoding scheme that could interpret the intention from streaming connectivity patterns, based on an ensemble of SVMs.

Results

The group-level analysis revealed increased phase synchrony and richer network organization in patients. This trend was also seen in the performance of the employed classifiers. Time-resolved connectivity led to superior performance, with distinct SVMs acting as local experts, specialized in the patterning emerged within specific temporal windows (defined with respect to the external trigger). This empirical finding was further exploited in implementing a decoding scheme that can be activated without the need of the precise timing of a trigger.

Conclusion

The increased phase synchrony in NMD patients can turn to a valuable tool for MI decoding. Considering the fast implementation for the PLV pattern computation in multichannel signals, we can envision the development of efficient personalized BCI systems in assistance of these patients.

Electronic supplementary material

The online version of this article (10.1186/s12984-018-0431-6) contains supplementary material, which is available to authorized users.

Source-Informed Segmentation: A Data-Driven Approach for the Temporal Segmentation of EEG

GOAL

Understanding the dynamics of brain function through non-invasive monitoring techniques requires the development of computational methods that can deal with the non-stationary properties of recorded activities. As a solution to this problem, a new data-driven segmentation method for recordings obtained through electroencephalography (EEG) is presented.

METHODS

The proposed method utilizes singular value decomposition (SVD) to identify the time intervals in the EEG recordings during which the spatial distribution of clusters of active cortical neurons remains quasi-stationary. Theoretical analysis shows that the spatial locality features of these clusters can be, asymptotically, captured by the most significant left singular subspace of the EEG data. A reference/sliding window approach is employed to dynamically extract this feature subspace, and the running projection error is monitored for significant changes using Kolmogorov-Smirnov test.

RESULTS

Simulation results, for a wide range of possible scenarios regarding the spatial distribution of active cortical neurons, show that the algorithm is successful in accurately detecting the segmental structure of the simulated EEG data. The algorithm is also applied to experimental EEG recordings of a modified visual oddball task. Results identify a unique sequence of dynamic patterns in the event-related potential (ERP) response to each of the three involved stimuli.

CONCLUSION

The proposed method, without using source localization methods or scalp topographical maps, is able to identify intervals of quasi-stationarity in the EEG recordings.

SIGNIFICANCE

The proposed segmentation technique can offer new insights on the dynamics of functional organization of the brain in action.

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.

Reliability of Static and Dynamic Network Metrics in the Resting-State: A MEG-Beamformed Connectivity Analysis

The resting activity of the brain can be described by so-called intrinsic connectivity networks (ICNs), which consist of spatially and temporally distributed, but functionally connected, nodes. The coordinated activity of the resting state can be explored via magnetoencephalography (MEG) by studying frequency-dependent functional brain networks at the source level. Although many algorithms for the analysis of brain connectivity have been proposed, the reliability of network metrics derived from both static and dynamic functional connectivity is still unknown. This is a particular problem for studies of associations between ICN metrics and personality variables or other traits, and for studies of differences between patient and control groups, which both depend critically on the reliability of the metrics used. A detailed investigation of the reliability of metrics derived from resting-state MEG repeat scans is therefore a prerequisite for the development of connectomic biomarkers. Here, we first estimated both static (SFC) and dynamic functional connectivity (DFC) after beamforming source reconstruction using the imaginary part of the phase locking index (iPLV) and the correlation of the amplitude envelope (CorEnv). Using our approach, functional network microstates (FCμstates) were derived from the DFC and chronnectomics were computed from the evolution of FCμstates across experimental time. In both temporal scales, the reliability of network metrics (SFC), the FCμstates and the related chronnectomics were evaluated for every frequency band. Chronnectomic statistics and FCμstates were generally more reliable than node-wise static network metrics. CorEnv-based network metrics were more reproducible at the static approach. The reliability of chronnectomics have been evaluated also in a second dataset. This study encourages the analysis of MEG resting-state via DFC.

Enhancing Performance and Bit Rates in a Brain-Computer Interface System With Phase-to-Amplitude Cross-Frequency Coupling: Evidences From Traditional c-VEP, Fast c-VEP, and SSVEP Designs

A brain–computer interface (BCI) is a channel of communication that transforms brain activity into specific commands for manipulating a personal computer or other home or electrical devices. In other words, a BCI is an alternative way of interacting with the environment by using brain activity instead of muscles and nerves. For that reason, BCI systems are of high clinical value for targeted populations suffering from neurological disorders. In this paper, we present a new processing approach in three publicly available BCI data sets: (a) a well-known multi-class (N = 6) coded-modulated Visual Evoked potential (c-VEP)-based BCI system for able-bodied and disabled subjects; (b) a multi-class (N = 32) c-VEP with slow and fast stimulus representation; and (c) a steady-state Visual Evoked potential (SSVEP) multi-class (N = 5) flickering BCI system. Estimating cross-frequency coupling (CFC) and namely δ-θ [δ: (0.5–4 Hz), θ: (4–8 Hz)] phase-to-amplitude coupling (PAC) within sensor and across experimental time, we succeeded in achieving high classification accuracy and Information Transfer Rates (ITR) in the three data sets. Our approach outperformed the originally presented ITR on the three data sets. The bit rates obtained for both the disabled and able-bodied subjects reached the fastest reported level of 324 bits/min with the PAC estimator. Additionally, our approach outperformed alternative signal features such as the relative power (29.73 bits/min) and raw time series analysis (24.93 bits/min) and also the original reported bit rates of 10–25 bits/min. In the second data set, we succeeded in achieving an average ITR of 124.40 ± 11.68 for the slow 60 Hz and an average ITR of 233.99 ± 15.75 for the fast 120 Hz. In the third data set, we succeeded in achieving an average ITR of 106.44 ± 8.94. Current methodology outperforms any previous methodologies applied to each of the three free available BCI datasets.

Mining Time-Resolved Functional Brain Graphs to an EEG-Based Chronnectomic Brain Aged Index (CBAI)

The brain at rest consists of spatially and temporal distributed but functionally connected regions that called intrinsic connectivity networks (ICNs). Resting state electroencephalography (rs-EEG) is a way to characterize brain networks without confounds associated with task EEG such as task difficulty and performance. A novel framework of how to study dynamic functional connectivity under the notion of functional connectivity microstates (FCμstates) and symbolic dynamics is further discussed. Furthermore, we introduced a way to construct a single integrated dynamic functional connectivity graph (IDFCG) that preserves both the strength of the connections between every pair of sensors but also the type of dominant intrinsic coupling modes (DICM). The whole methodology is demonstrated in a significant and unexplored task for EEG which is the definition of an objective Chronnectomic Brain Aged index (CBAI) extracted from resting-state data (N = 94 subjects) with both eyes-open and eyes-closed conditions. Novel features have been defined based on symbolic dynamics and the notion of DICM and FCμstates. The transition rate of FCμstates, the symbolic dynamics based on the evolution of FCμstates (the Markovian Entropy, the complexity index), the probability distribution of DICM, the novel Flexibility Index that captures the dynamic reconfiguration of DICM per pair of EEG sensors and the relative signal power constitute a valuable pool of features that can build the proposed CBAI. Here we applied a feature selection technique and Extreme Learning Machine (ELM) classifier to discriminate young adults from middle-aged and a Support Vector Regressor to build a linear model of the actual age based on EEG-based spatio-temporal features. The most significant type of features for both prediction of age and discrimination of young vs. adults age groups was the dynamic reconfiguration of dominant coupling modes derived from a subset of EEG sensor pairs. Specifically, our results revealed a very high prediction of age for eyes-open (R² = 0.60; y = 0.79x + 8.03) and lower for eyes-closed (R² = 0.48; y = 0.71x + 10.91) while we succeeded to correctly classify young vs. middle-age group with 97.8% accuracy in eyes-open and 87.2% for eyes-closed. Our results were reproduced also in a second dataset for further external validation of the whole analysis. The proposed methodology proved valuable for the characterization of the intrinsic properties of dynamic functional connectivity through the age untangling developmental differences using EEG resting-state recordings.

Altered Rich-Club and Frequency-Dependent Subnetwork Organization in Mild Traumatic Brain Injury: A MEG Resting-State Study

Functional brain connectivity networks exhibit "small-world" characteristics and some of these networks follow a "rich-club" organization, whereby a few nodes of high connectivity (hubs) tend to connect more densely among themselves than to nodes of lower connectivity. The Current study followed an "attack strategy" to compare the rich-club and small-world network organization models using Magnetoencephalographic (MEG) recordings from mild traumatic brain injury (mTBI) patients and neurologically healthy controls to identify the topology that describes the underlying intrinsic brain network organization. We hypothesized that the reduction in global efficiency caused by an attack targeting a model's hubs would reveal the "true" underlying topological organization. Connectivity networks were estimated using mutual information as the basis for cross-frequency coupling. Our results revealed a prominent rich-club network organization for both groups. In particular, mTBI patients demonstrated hyper-synchronization among rich-club hubs compared to controls in the δ band and the δ-γ1, θ-γ1, and β-γ2 frequency pairs. Moreover, rich-club hubs in mTBI patients were overrepresented in right frontal brain areas, from θ to γ1 frequencies, and underrepresented in left occipital regions in the δ-β, δ-γ1, θ-β, and β-γ2 frequency pairs. These findings indicate that the rich-club organization of resting-state MEG, considering its role in information integration and its vulnerability to various disorders like mTBI, may have a significant predictive value in the development of reliable biomarkers to help the validation of the recovery from mTBI. Furthermore, the proposed approach might be used as a validation tool to assess patient recovery.

Topological Filtering of Dynamic Functional Brain Networks Unfolds Informative Chronnectomics: A Novel Data-Driven Thresholding Scheme Based on Orthogonal Minimal Spanning Trees (OMSTs)

The human brain is a large-scale system of functionally connected brain regions. This system can be modeled as a network, or graph, by dividing the brain into a set of regions, or “nodes,” and quantifying the strength of the connections between nodes, or “edges,” as the temporal correlation in their patterns of activity. Network analysis, a part of graph theory, provides a set of summary statistics that can be used to describe complex brain networks in a meaningful way. The large-scale organization of the brain has features of complex networks that can be quantified using network measures from graph theory. The adaptation of both bivariate (mutual information) and multivariate (Granger causality) connectivity estimators to quantify the synchronization between multichannel recordings yields a fully connected, weighted, (a)symmetric functional connectivity graph (FCG), representing the associations among all brain areas. The aforementioned procedure leads to an extremely dense network of tens up to a few hundreds of weights. Therefore, this FCG must be filtered out so that the “true” connectivity pattern can emerge. Here, we compared a large number of well-known topological thresholding techniques with the novel proposed data-driven scheme based on orthogonal minimal spanning trees (OMSTs). OMSTs filter brain connectivity networks based on the optimization between the global efficiency of the network and the cost preserving its wiring. We demonstrated the proposed method in a large EEG database (N = 101 subjects) with eyes-open (EO) and eyes-closed (EC) tasks by adopting a time-varying approach with the main goal to extract features that can totally distinguish each subject from the rest of the set. Additionally, the reliability of the proposed scheme was estimated in a second case study of fMRI resting-state activity with multiple scans. Our results demonstrated clearly that the proposed thresholding scheme outperformed a large list of thresholding schemes based on the recognition accuracy of each subject compared to the rest of the cohort (EEG). Additionally, the reliability of the network metrics based on the fMRI static networks was improved based on the proposed topological filtering scheme. Overall, the proposed algorithm could be used across neuroimaging and multimodal studies as a common computationally efficient standardized tool for a great number of neuroscientists and physicists working on numerous of projects.

Brain Network Activation Analysis Utilizing Spatiotemporal Features for Event Related Potentials Classification

The purpose of this study was to introduce an improved tool for automated classification of event-related potentials (ERPs) using spatiotemporally parcellated events incorporated into a functional brain network activation (BNA) analysis. The auditory oddball ERP paradigm was selected to demonstrate and evaluate the improved tool. Methods: The ERPs of each subject were decomposed into major dynamic spatiotemporal events. Then, a set of spatiotemporal events representing the group was generated by aligning and clustering the spatiotemporal events of all individual subjects. The temporal relationship between the common group events generated a network, which is the spatiotemporal reference BNA model. Scores were derived by comparing each subject's spatiotemporal events to the reference BNA model and were then entered into a support vector machine classifier to classify subjects into relevant subgroups. The reliability of the BNA scores (test-retest repeatability using intraclass correlation) and their utility as a classification tool were examined in the context of Target-Novel classification. Results: BNA intraclass correlation values of repeatability ranged between 0.51 and 0.82 for the known ERP components N100, P200, and P300. Classification accuracy was high when the trained data were validated on the same subjects for different visits (AUCs 0.93 and 0.95). The classification accuracy remained high for a test group recorded at a different clinical center with a different recording system (AUCs 0.81, 0.85 for 2 visits). Conclusion: The improved spatiotemporal BNA analysis demonstrates high classification accuracy. The BNA analysis method holds promise as a tool for diagnosis, follow-up and drug development associated with different neurological conditions.

Mnemonic strategy training of the elderly at risk for dementia enhances integration of information processing via cross-frequency coupling

Introduction

We sought to identify whether intensive 10-week mobile health mnemonic strategy training (MST) could shift the resting-state brain network more toward cortical-level integration, which has recently been proven to reflect the reorganization of the brain networks compensating the cognitive decline.

Methods

One hundred fifty-eight patients with mild cognitive impairment (MCI) were selected and participated in 10-week training lasting 90 min/d of memory training. They benefited from an initial and a follow-up neuropsychological evaluation and resting-state electroencephalography (EEG).

Results

At follow-up, MST revealed an extensive significant training effect that changed the network with an increase of synchronization between parietotemporal and frontal areas; frontal^θ-parietal^α2 causal strengthening as part of top-down inhibitory control; enhancement of sensorimotor connections in β band; and a general increase of cortical-level integration. More precisely, MST induced gain as an increase of the global cost efficiency (GCE) of the whole cortical network and a neuropsychological performance improvement, which was correlated with it (r = 0.32, P = .0001). The present study unfolded intervention changes based on EEG source activity via novel neuroinformatic tools for revealing intrinsic coupling modes in both amplitude-phase representations and in the mixed spectrospatiotemporal domain.

Discussion

Further work should identify whether the GCE enhancement of the functional cortical brain networks is a compensation mechanism to the brain network dysfunction or a more permanent neuroplasticity effect.

Mining cross-frequency coupling microstates (CFCμstates) from EEG recordings during resting state and mental arithmetic tasks

The functional brain connectivity has been studied by analyzing synchronization between dynamic oscillations of identical frequency or between different frequencies of distinct brain areas. It has been hypothesized that cross-frequency coupling (CFC) between different frequency bands is the carrier mechanism for the coordination of global and local neural processes and hence supports the distributed information processing in the brain. In the present study, we attempt to study the dynamic evolution of CFC at resting-state and during a mental task. The concept of CFC microstates (CFCμstates) is introduced as emerged short-lived patterns of CFC. We analyzed dynamic CFC (dCFC) at resting-state and during a comparison task by adopting a phase-amplitude coupling (PAC) estimator for [δ phase-γ-amplitude] coupling at every sensor. Modifying a well-established framework for mining brain dynamics, we show that a small sized repertoire of CFCμstates can be derived so as to encapsulate connectivity variations and further provide novel insights into network's functional reorganization. By analyzing the transition dynamics among CFCμstates, in both tasks, we provided a clear evidence about intrinsic networks that may play a crucial role in information integration.

Mining cross-frequency coupling microstates from resting state MEG: An application to mild traumatic brain injury

Recent studies have investigated the possible role of dynamic functional connectivity and the role of cross-frequency coupling (CFC) to provide the substrate for reliable biomarkers of brain disorders. In this study, we analyzed time-varying CFC profiles from resting state Magnetoencephal-ographic recordings of 30 mild Traumatic Brain Injury (mTBI) patients and 50 normal controls. Interactions among sensors at specific pairs of frequency bands were computed via estimation of phase-to-amplitude couplings. We then computed time-varying functional connectivity graphs that were described in terms of segregation (local efficiency, LE) and integration (global efficiency, GE) and mapped those graphs to time series of GE/LE estimates. The resulting dynamic network revealed transitions between a limited number of microstates for mTBI subjects compared to controls. The significant differences in transition probability between the two groups, along with the limited repertoire of possible states, can form the basis for a robust dynamic connectomic biomarker for the diagnosis of mTBI.

Phase to amplitude coupling as a potential biomarker for creative ideation: an EEG study

The most consistent finding of creative ideation in the neuroscientific study of creativity is the increment of EEG α power. However, the majority of existing studies focused only on ERP experimental paradigms while only a few analyzed time-related changes of EEG α power patterns during the time unlocked creation of ideas. Here, we designed an experimental paradigm where the participants were asked to generate alternative uses of everyday objects (AU task). For the control task, we adopted an Object Characteristics (OC) task, for which participants were asked to list typical characteristics or properties of an object. We estimated relative power spectrum, global efficiency from brain networks constructed with the imaginary part of coherence and phase-to-amplitude coupling (PAC) as potential biomarkers of creativity. Both relative power spectrum and nodal global efficiency failed to reach significant level by comparing AU with OC. In contrast, statistically significant differences between AU and OC were detected with PAC estimated within sensors in frequency pairs θ-γ and α2-γ. Our results can be the ground for both detecting and designing a connectomic biomarker of creativity.

Greater Repertoire and Temporal Variability of Cross-Frequency Coupling (CFC) Modes in Resting-State Neuromagnetic Recordings among Children with Reading Difficulties

Cross-frequency, phase-to-amplitude coupling (PAC) between neuronal oscillations at rest may serve as the substrate that supports information exchange between functionally specialized neuronal populations both within and between cortical regions. The study utilizes novel algorithms to identify prominent instantaneous modes of cross-frequency coupling and their temporal stability in resting state magnetoencephalography (MEG) data from 25 students experiencing severe reading difficulties (RD) and 27 age-matched non-impaired readers (NI). Phase coherence estimates were computed in order to identify the prominent mode of PAC interaction for each sensor, sensor pair, and pair of frequency bands (from δ to γ) at successive time windows of the continuous MEG record. The degree of variability in the characteristic frequency-pair PAC^f1−f2 modes over time was also estimated. Results revealed a wider repertoire of prominent PAC interactions in RD as compared to NI students, suggesting an altered functional substrate for information exchange between neuronal assemblies in the former group. Moreover, RD students showed significant variability in PAC modes over time. This temporal instability of PAC values was particularly prominent: (a) within and between right hemisphere temporo-parietal and occipito-temporal sensors and, (b) between left hemisphere frontal, temporal, and occipito-temporal sensors and corresponding right hemisphere sites. Altered modes of neuronal population coupling may help account for extant data revealing reduced, task-related neurophysiological and hemodynamic activation in left hemisphere regions involved in the reading network in RD. Moreover, the spatial distribution of pronounced instability of cross-frequency coupling modes in this group may provide an explanation for previous reports suggesting the presence of inefficient compensatory mechanisms to support reading.

A novel biomarker of amnestic MCI based on dynamic cross-frequency coupling patterns during cognitive brain responses

The detection of mild cognitive impairment (MCI), the transitional stage between normal cognitive changes of aging and the cognitive decline caused by AD, is of paramount clinical importance, since MCI patients are at increased risk of progressing into AD. Electroencephalographic (EEG) alterations in the spectral content of brainwaves and connectivity at resting state have been associated with early-stage AD. Recently, cognitive event-related potentials (ERPs) have entered into the picture as an easy to perform screening test. Motivated by the recent findings about the role of cross-frequency coupling (CFC) in cognition, we introduce a relevant methodological approach for detecting MCI based on cognitive responses from a standard auditory oddball paradigm. By using the single trial signals recorded at Pz sensor and comparing the responses to target and non-target stimuli, we first demonstrate that increased CFC is associated with the cognitive task. Then, considering the dynamic character of CFC, we identify instances during which the coupling between particular pairs of brainwave frequencies carries sufficient information for discriminating between normal subjects and patients with MCI. In this way, we form a multiparametric signature of impaired cognition. The new composite biomarker was tested using data from a cohort that consists of 25 amnestic MCI patients and 15 age-matched controls. Standard machine-learning algorithms were employed so as to implement the binary classification task. Based on leave-one-out cross-validation, the measured classification rate was found reaching very high levels (95%). Our approach compares favorably with the traditional alternative of using the morphology of averaged ERP response to make the diagnosis and the usage of features from spectro-temporal analysis of single-trial responses. This further indicates that task-related CFC measurements can provide invaluable analytics in AD diagnosis and prognosis.

Transition dynamics of EEG-based network microstates during mental arithmetic and resting wakefulness reflects task-related modulations and developmental changes

We studied how maturation influences the organization of functional brain networks engaged during mental calculations and in resting state. Surface EEG measurements from 20 children (8-12 years) and 25 students (21-26 years) were analyzed. Interregional synchronization of brain activity was quantified by means of Phase Lag Index and for various frequency bands. Based on these pairwise estimates of functional connectivity, we formed graphs which were then characterized in terms of local structure [local efficiency (LE)] and overall integration (global efficiency). The overall data analytic scheme was applied twice, in a static and time-varying mode. Our results showed a characteristic trend: functional segregation dominates the network organization of younger brains. Moreover, in childhood, the overall functional network possesses more prominent small-world network characteristics than in early acorrect in xmldulthood in accordance with the Neural Efficiency Hypothesis. The above trends were intensified by the time-varying approach and identified for the whole set of tested frequency bands (from δ to low γ). By mapping the time-indexed connectivity patterns to multivariate timeseries of nodal LE measurements, we carried out an elaborate study of the functional segregation dynamics and demonstrated that the underlying network undergoes transitions between a restricted number of stable states, that can be thought of as "network-level microstates". The rate of these transitions provided a robust marker of developmental and task-induced alterations, that was found to be insensitive to reference montage and independent component analysis denoising.

A graph theoretic approach to dynamic functional connectivity tracking and network state identification

With the advances in neuroimaging technology, it is now possible to measure human brain activity with increasing temporal and spatial resolution. This vast amount of spatio-temporal data requires the development of computational methods capable of building an integrated picture of the functional networks for a better understanding of the healthy and diseased brain [1]. Although the construction of these networks from neuroimaging data is well-established [2], current approaches are limited to the characterization of the global topology of static networks where the links between different brain regions represent average connectivity over a long time period [3], [2]. Recent studies suggest that human cognition arises from the rapid formation and dissociation of synchronized neural activity on short time scales in the order of milliseconds [4]. There is a strong need for new electroencephalogram (EEG)-based analytic frameworks for monitoring dynamic functional network activity. In this paper, we propose a graph theoretic approach for tracking the changing topology of functional connectivity networks across time. First, we introduce an event detection algorithm based on node level feature extraction and principal components analysis of time-dependent node correlation matrices. Then, we propose a k-means based clustering approach to characterize each time interval with the most common connectivity states. Finally, the proposed methodology is applied to the study of the dynamics of functional connectivity networks during error-related negativity (ERN).