Ordinal methods for a characterization of evolving functional brain networks

Ordinal time series analysis is based on the idea to map time series to ordinal patterns, i.e., order relations between the values of a time series and not the values themselves, as introduced in 2002 by C. Bandt and B. Pompe. Despite a resulting loss of information, this approach captures meaningful information about the temporal structure of the underlying system dynamics as well as about properties of interactions between coupled systems. This-together with its conceptual simplicity and robustness against measurement noise-makes ordinal time series analysis well suited to improve characterization of the still poorly understood spatiotemporal dynamics of the human brain. This minireview briefly summarizes the state-of-the-art of uni- and bivariate ordinal time-series-analysis techniques together with applications in the neurosciences. It will highlight current limitations to stimulate further developments, which would be necessary to advance characterization of evolving functional brain networks.

Dynamic Reconfiguration of Dominant Intrinsic Coupling Modes in Elderly at Prodromal Alzheimer's Disease Risk

Large-scale human brain networks interact across both spatial and temporal scales. Especially for electro- and magnetoencephalography (EEG/MEG), there are many evidences that there is a synergy of different subnetworks that oscillate on a dominant frequency within a quasi-stable brain temporal frame. Intrinsic cortical-level integration reflects the reorganization of functional brain networks that support a compensation mechanism for cognitive decline. Here, a computerized intervention integrating different functions of the medial temporal lobes, namely, object-level and scene-level representations, was conducted. One hundred fifty-eight patients with mild cognitive impairment underwent 90 min of training per day over 10 weeks. An active control (AC) group of 50 subjects was exposed to documentaries, and a passive control group of 55 subjects did not engage in any activity. Following a dynamic functional source connectivity analysis, the dynamic reconfiguration of intra- and cross-frequency coupling mechanisms before and after the intervention was revealed. After the neuropsychological and resting state electroencephalography evaluation, the ratio of inter versus intra-frequency coupling modes and also the contribution of β1 frequency was higher for the target group compared to its pre-intervention period. These frequency-dependent contributions were linked to neuropsychological estimates that were improved due to intervention. Additionally, the time-delays of the cortical interactions were improved in {δ, θ, α2, β1} compared to the pre-intervention period. Finally, dynamic networks of the target group further improved their efficiency over the total cost of the network. This is the first study that revealed a dynamic reconfiguration of intrinsic coupling modes and an improvement of time-delays due to a target intervention protocol.

Transcerebral information coordination in directional hippocampus-prefrontal cortex network during working memory based on bimodal neural electrical signals

Working memory (WM) is a kind of advanced cognitive function, which requires the participation and cooperation of multiple brain regions. Hippocampus and prefrontal cortex are the main responsible brain regions for WM. Exploring information coordination between hippocampus and prefrontal cortex during WM is a frontier problem in cognitive neuroscience. In this paper, an advanced information theory analysis based on bimodal neural electrical signals (local field potentials, LFPs and spikes) was employed to characterize the transcerebral information coordination across the two brain regions. Firstly, LFPs and spikes were recorded simultaneously from rat hippocampus and prefrontal cortex during the WM task by using multi-channel in vivo recording technique. Then, from the perspective of information theory, directional hippocampus-prefrontal cortex networks were constructed by using transfer entropy algorithm based on spectral coherence between LFPs and spikes. Finally, transcerebral coordination of bimodal information at the brain-network level was investigated during acquisition and performance of the WM task. The results show that the transfer entropy in directional hippocampus-prefrontal cortex networks is related to the acquisition and performance of WM. During the acquisition of WM, the information flow, local information transmission ability and information transmission efficiency of the directional hippocampus-prefrontal networks increase over learning days. During the performance of WM, the transfer entropy from the hippocampus to prefrontal cortex plays a leading role for bimodal information coordination across brain regions and hippocampus has a driving effect on prefrontal cortex. Furthermore, bimodal information coordination in the hippocampus → prefrontal cortex network could predict WM during the successful performance of WM.

Universal Lifespan Trajectories of Source-Space Information Flow Extracted from Resting-State MEG Data

Source activity was extracted from resting-state magnetoencephalography data of 103 subjects aged 18-60 years. The directionality of information flow was computed from the regional time courses using delay symbolic transfer entropy and phase entropy. The analysis yielded a dynamic source connectivity profile, disentangling the direction, strength, and time delay of the underlying causal interactions, producing independent time delays for cross-frequency amplitude-to-amplitude and phase-to-phase coupling. The computation of the dominant intrinsic coupling mode (DoCM) allowed me to estimate the probability distribution of the DoCM independently of phase and amplitude. The results support earlier observations of a posterior-to-anterior information flow for phase dynamics in {α1, α2, β, γ} and an opposite flow (anterior to posterior) in θ. Amplitude dynamics reveal posterior-to-anterior information flow in {α1, α2, γ}, a sensory-motor β-oriented pattern, and an anterior-to-posterior pattern in {δ, θ}. The DoCM between intra- and cross-frequency couplings (CFC) are reported here for the first time and independently for amplitude and phase; in both domains {δ, θ, α1}, frequencies are the main contributors to DoCM. Finally, a novel brain age index (BAI) is introduced, defined as the ratio of the probability distribution of inter- over intra-frequency couplings. This ratio shows a universal age trajectory: a rapid rise from the end of adolescence, reaching a peak in adulthood, and declining slowly thereafter. The universal pattern is seen in the BAI of each frequency studied and for both amplitude and phase domains. No such universal age dependence was previously reported.

EEG emotion recognition based on cross-frequency granger causality feature extraction and fusion in the left and right hemispheres

EEG emotion recognition based on Granger causality (GC) brain networks mainly focus on the EEG signal from the same-frequency bands, however, there are still some causality relationships between EEG signals in the cross-frequency bands. Considering the functional asymmetric of the left and right hemispheres to emotional response, this paper proposes an EEG emotion recognition scheme based on cross-frequency GC feature extraction and fusion in the left and right hemispheres. Firstly, we calculate the GC relationship of EEG signals according to the frequencies and hemispheres, and mainly focus on the causality of the cross-frequency EEG signals in left and right hemispheres. Then, to remove the redundant connections of the GC brain network, an adaptive two-stage decorrelation feature extraction scheme is proposed under the condition of maintaining the best emotion recognition performance. Finally, a multi-GC feature fusion scheme is designed to balance the recognition accuracy and feature number of each GC feature, which comprehensively considers the influence of the recognition accuracy and computational complexity. Experimental results on the DEAP emotion dataset show that the proposed scheme can achieve an average accuracy of 84.91% for four classifications, which improved the classification accuracy by up to 8.43% compared with that of the traditional same-frequency band GC features.

A modified phase transfer entropy for cross-frequency directed coupling estimation in brain network

Cross-frequency coupling of neural oscillation is widespread during the complex cognitive process. Therefore, identifying cross-frequency information flow is essential for revealing neural dynamics mechanisms in the brain network. A current method based on the information theory, phase transfer entropy (PTE), has been proved its effectiveness in estimating directional coupling in several recent studies. However, there remains some limits in PTE: (1)lack of multivariable effect, (2) poor robustness, (3)curse of dimensionality in the high dimensional system. This study introduced a novel multivariate phase transfer entropy method named "MPTE_NUE" to solve the above issues. In MPTE_NUE, it considered the influence of remaining confounding variables, which guaranteed its applicability in a multivariable system. Meanwhile, a nonuniform embedding (NUE) approach for state reconstruction was adopted to eliminate the dimensional curse problem. We performed a series of numerical simulations based on the typical Hénon map model. The results proved that the MPTE_NUE achieved better noise robustness and effectively avoided the curse of dimension; meanwhile, the accuracy and sensitivity can reach 96.9% and 99.2%, respectively.

Altered directed connectivity during processing of predictive stimuli in psychiatric patient populations

OBJECTIVES

The study investigated the role of top-down versus bottom-up connectivity, during the processing of predictive information, in three different psychiatric disorders.

METHODS

Electroencephalography (EEG) was recorded during the performance of a task, which evaluates the ability to use predictive information in order to facilitate predictable versus random target detection. We evaluated EEG event-related directed connectivity, in patients with schizophrenia (SZ), major depressive disorder (MDD), and autism spectrum disorder (ASD), compared with healthy age-matched controls. Directed connectivity was evaluated using phase transfer entropy.

RESULTS

We showed that top-down frontal-parietal connectivity was weaker in SZ (theta and beta bands) and ASD (alpha band) compared to control subjects, during the processing of stimuli consisting of the predictive sequence. In SZ patients, top-down connectivity was also attenuated, during the processing of predictive targets in the beta frequency band. In contrast, compared with controls, MDD patients displayed an increased top-down flow of information, during the processing of predicted targets (alpha band).

CONCLUSIONS

The findings suggest that top-down frontal-parietal connectivity is altered differentially across three major psychiatric disorders, specifically during the processing of predictive stimuli.

SIGNIFICANCE

Altered top-down connectivity may contribute to the specific prediction deficits observed in each of the patient populations.

Altered directed connectivity during processing of implicit versus explicit predictive stimuli in Parkinson's disease patients

The study investigated the role of top-down versus bottom-up connectivity, during the processing of implicit or explicit predictive information, in Parkinson's disease (PD). EEG was recorded during the performance of a task, which evaluated the ability to utilize either implicit or explicit predictive contextual information in order to facilitate the detection of predictable versus random targets. Thus, subjects performed an implicit and explicit session, where subjects were either unaware or made aware of a predictive sequence that signals the presentation of a subsequent target, respectively. We evaluated EEG event-related directed connectivity, in PD patients compared with healthy age-matched controls, using phase transfer entropy. PD patients showed increased top-down frontal-parietal connectivity, compared to control subjects, during the processing of the last (most informative) stimulus of the predictive sequence and of random standards, in the implicit and explicit session, respectively. These findings suggest that PD is associated with compensatory top-down connectivity, specifically during the processing of implicit predictive stimuli. During the explicit session, PD patients seem to allocate more attentional resources to non-informative standard stimuli, compared to controls. These connectivity changes shed further light on the cognitive deficits, associated with the processing of predictive contextual information, that are observed in PD patients.

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.

Covariation Informed Graph Slepians for Motor Imagery Decoding

Graph signal processing (GSP) provides signal analytic tools for data defined in irregular domains, as is the case of non-invasive electroencephalography (EEG). In this work, the recently introduced technique of Graph Slepian functions is exploited for the robust decoding of motor imagery (MI) brain activity. The particular technique builds over the concept of graph Fourier transform (GFT) and provides additional flexibility in the subsequent data analysis by incorporating domain knowledge. Based on contrastive learning, we introduce an algorithmic pipeline that attains a data driven and subject specific design of Graph Slepian functions. These functions, by incorporating both the topology of the sensor array and the empirical evidence about the differential functional covariation, act as spatial filters that enhance the information conveyed by the multichannel signal and specifically relates to the participant's intention. The proposed technique for crafting Graph Slepians is incorporated in a MI-decoding scheme, in which the informed projections are fed to a support vector machine (SVM) that casts a prediction regarding the type of intended movement. The employed MI-decoder is evaluated based on two publicly available datasets and its superiority against popular alternatives in the field is established. Computational efficiency is listed among its main advantages, since it involves only simple matrix operations, allowing to consider its use in real-time implementations.

Effective brain connectivity for fNIRS data analysis based on multi-delays symbolic phase transfer entropy

OBJECTIVE

Recently, effective connectivity (EC) calculation methods for functional near-infrared spectroscopy (fNIRS) data mainly face two problems: the first problem is that noise can seriously affect the EC calculation and even lead to false connectivity; the second problem is that it ignores the various real neurotransmission delays between the brain region, and instead uses a fixed delay coefficient for calculation.

APPROACH

To overcome these two issues, a delay symbolic phase transfer entropy (dSPTE) is proposed by developing traditional transfer entropy (TE) to estimate EC for fNIRS. Firstly, the phase time sequence was obtained from the original sequence by the Hilbert transform and state-space reconstruction was realized using a uniform embedding scheme. Then, a symbolization technique was applied based on a neural-gas algorithm to improve its noise robustness. Finally, the EC was calculated on multiple time delay scales to match different inter-region neurotransmission delays.

MAIN RESULTS

A linear AR model, a nonlinear model and a multivariate hybrid model were introduced to simulate the performance of dSPTE, and the results showed that the accuracy of dSPTE was the highest, up to 74.27%, and specificity was 100% which means no false connectivity. The results confirmed that the dSPTE method realized better noise robustness, higher accuracy, and correct identification even if there was a long delay between series. Finally, we applied dSPTE to fNIRS dataset to analyse the EC during the finger-tapping task, the results showed that EC strength of task state significantly increased compared with the resting state.

SIGNIFICANCE

The proposed dSPTE method is a promising way to measure the EC for fNIRS. It incorporates the phase information TE with a symbolic process for fNIRS analysis for the first time. It has been confirmed to be noise robust and suitable for the complex network with different coupling delays.

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.

A Data-Driven Measure of Effective Connectivity Based on Renyi's α-Entropy

Transfer entropy (TE) is a model-free effective connectivity measure based on information theory. It has been increasingly used in neuroscience because of its ability to detect unknown non-linear interactions, which makes it well suited for exploratory brain effective connectivity analyses. Like all information theoretic quantities, TE is defined regarding the probability distributions of the system under study, which in practice are unknown and must be estimated from data. Commonly used methods for TE estimation rely on a local approximation of the probability distributions from nearest neighbor distances, or on symbolization schemes that then allow the probabilities to be estimated from the symbols' relative frequencies. However, probability estimation is a challenging problem, and avoiding this intermediate step in TE computation is desirable. In this work, we propose a novel TE estimator using functionals defined on positive definite and infinitely divisible kernels matrices that approximate Renyi's entropy measures of order α. Our data-driven approach estimates TE directly from data, sidestepping the need for probability distribution estimation. Also, the proposed estimator encompasses the well-known definition of TE as a sum of Shannon entropies in the limiting case when α → 1. We tested our proposal on a simulation framework consisting of two linear models, based on autoregressive approaches and a linear coupling function, respectively, and on the public electroencephalogram (EEG) database BCI Competition IV, obtained under a motor imagery paradigm. For the synthetic data, the proposed kernel-based TE estimation method satisfactorily identifies the causal interactions present in the data. Also, it displays robustness to varying noise levels and data sizes, and to the presence of multiple interaction delays in the same connected network. Obtained results for the motor imagery task show that our approach codes discriminant spatiotemporal patterns for the left and right-hand motor imagination tasks, with classification performances that compare favorably to the state-of-the-art.

Connectivity steered graph Fourier transform for motor imagery BCI decoding

OBJECTIVE

Graph signal processing (GSP) concepts are exploited for brain activity decoding and particularly the detection and recognition of a motor imagery (MI) movement. A novel signal analytic technique that combines graph Fourier transform (GFT) with estimates of cross-frequency coupling (CFC) and discriminative learning is introduced as a means to recover the subject's intention from the multichannel signal.

APPROACH

Adopting a multi-view perspective, based on the popular concept of co-existing and interacting brain rhythms, a multilayer network model is first built from empirical data and its connectivity graph is used to derive the GFT-basis. A personalized decoding scheme supporting a binary decision, either 'left versus right' or 'rest versus MI', is crafted from a small set of training trials. Electroencephalographic (EEG) activity from 12 volunteers recorded during two randomly alternating, externally cued, MI tasks (clenching either left or right fist) and a rest condition is used to introduce and validate our methodology. In addition, the introduced methodology was further validated based on dataset IVa of BCI III competition.

MAIN RESULTS

Our GFT-domain decoding scheme achieves nearly optimal performance and proves superior to alternative techniques that are very popular in the field.

SIGNIFICANCE

At a conceptual level, our work suggests a fruitful way to introduce network neuroscience in BCI research. At a more practical level, it is characterized by efficiency. Training is realized using a small number of exemplar trials and decoding requires very simple operations that leaves room for real-time implementation.

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.

Estimating the Memory Order of Electrocorticography Recordings

OBJECTIVE

This paper presents a data-driven method for estimating the memory order (the average length of the statistical dependence of a given sample on previous samples) of a recorded electrocorticography (ECoG) sequence.

METHODS

The proposed inference method is based on the relationship between the loss in predicting the next sample in a time-series and the dependence of this sample on the previous samples. Specifically, the memory order is estimated to be the number of past samples that minimize the least squares error (LSE) in predicting the next sample. To deal with the lack of an analytical model for ECoG recordings, the proposed method combines a collection of different predictors, thereby achieving LSE at least as low as the LSE achieved by each of the different predictors.

RESULTS

ECoG recordings from six patients with epilepsy were analyzed, and the empirical cumulative density functions (ECDFs) of the memory orders estimated from these recordings were generated, for rest as well as pre-ictal time intervals. For pre-ictal time intervals, the electrodes corresponding to the seizure-onset-zone were separately analyzed. The estimated ECDFs were different between patients and between different types of blocks. For all the analyzed patients, the estimated memory orders were on the order of tens of milliseconds (up to 100 ms).

SIGNIFICANCE

The proposed method facilitates the estimation of the causal associations between ECoG recordings, as these associations strongly depend on the recordings' memory. An improved estimation of causal associations can improve the performance of algorithms that use ECoG recordings to localize the epileptogenic zone. Such algorithms can aid doctors in their pre-surgical planning for the surgery of patients with epilepsy.

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.

Data-Driven Topological Filtering Based on Orthogonal Minimal Spanning Trees: Application to Multigroup Magnetoencephalography Resting-State Connectivity

In the present study, a novel data-driven topological filtering technique is introduced to derive the backbone of functional brain networks relying on orthogonal minimal spanning trees (OMSTs). The method aims to identify the essential functional connections to ensure optimal information flow via the objective criterion of global efficiency minus the cost of surviving connections. The OMST technique was applied to multichannel, resting-state neuromagnetic recordings from four groups of participants: healthy adults (n = 50), adults who have suffered mild traumatic brain injury (n = 30), typically developing children (n = 27), and reading-disabled children (n = 25). Weighted interactions between network nodes (sensors) were computed using an integrated approach of dominant intrinsic coupling modes based on two alternative metrics (symbolic mutual information and phase lag index), resulting in excellent discrimination of individual cases according to their group membership. Classification results using OMST-derived functional networks were clearly superior to results using either relative power spectrum features or functional networks derived through the conventional minimal spanning tree algorithm.

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.

fMRI classification method with multiple feature fusion based on minimum spanning tree analysis

Resting state functional brain networks have been widely studied in brain disease research. Conventional network analysis methods are hampered by differences in network size, density and normalization. Minimum spanning tree (MST) analysis has been recently suggested to ameliorate these limitations. Moreover, common MST analysis methods involve calculating quantifiable attributes and selecting these attributes as features in the classification. However, a disadvantage of these methods is that information about the topology of the network is not fully considered, limiting further improvement of classification performance. To address this issue, we propose a novel method combining brain region and subgraph features for classification, utilizing two feature types to quantify two properties of the network. We experimentally validated our proposed method using a major depressive disorder (MDD) patient dataset. The results indicated that MSTs of MDD patients were more similar to random networks and exhibited significant differences in certain regions involved in the limbic-cortical-striatal-pallidal-thalamic (LCSPT) circuit, which is considered to be a major pathological circuit of depression. Moreover, we demonstrated that this novel classification method could effectively improve classification accuracy and provide better interpretability. Overall, the current study demonstrated that different forms of feature representation provide complementary information.

Propofol-induced sedation diminishes the strength of frontal-parietal-occipital EEG network

The level of conscious experience can be effectively and reversibly altered by the administration of sedative agents. Several studies attempted to explore the variations in frontal-parietal network during propofol-induced sedation. However, contradictory outcomes warrant further investigations. In this study, we implemented the Neural Gas algorithm-based delay symbolic transfer entropy (NG-dSTE) for investigation of frontal-parietal-occipital (F-P-O) network using scalp EEG signals recorded during altered levels of consciousness. Our results show significant disruption of the F-P-O network during mild and moderate levels of propofol sedation. In particular, the interaction between frontal and parietal-occipital region is highly disturbed. Moreover, we found measurable effect of sedation on local interactions in the frontal network whereas parietal-occipital network experienced least variations. The results support the conclusion that the connectivity based features can be utilized as reliable biomarker for assessment of sedation levels effectively.

How to Build a Functional Connectomic Biomarker for Mild Cognitive Impairment From Source Reconstructed MEG Resting-State Activity: The Combination of ROI Representation and Connectivity Estimator Matters

Our work aimed to demonstrate the combination of machine learning and graph theory for the designing of a connectomic biomarker for mild cognitive impairment (MCI) subjects using eyes-closed neuromagnetic recordings. The whole analysis based on source-reconstructed neuromagnetic activity. As ROI representation, we employed the principal component analysis (PCA) and centroid approaches. As representative bi-variate connectivity estimators for the estimation of intra and cross-frequency interactions, we adopted the phase locking value (PLV), the imaginary part (iPLV) and the correlation of the envelope (CorrEnv). Both intra and cross-frequency interactions (CFC) have been estimated with the three connectivity estimators within the seven frequency bands (intra-frequency) and in pairs (CFC), correspondingly. We demonstrated how different versions of functional connectivity graphs single-layer (SL-FCG) and multi-layer (ML-FCG) can give us a different view of the functional interactions across the brain areas. Finally, we applied machine learning techniques with main scope to build a reliable connectomic biomarker by analyzing both SL-FCG and ML-FCG in two different options: as a whole unit using a tensorial extraction algorithm and as single pair-wise coupling estimations. We concluded that edge-weighed feature selection strategy outperformed the tensorial treatment of SL-FCG and ML-FCG. The highest classification performance was obtained with the centroid ROI representation and edge-weighted analysis of the SL-FCG reaching the 98% for the CorrEnv in α1:α2 and 94% for the iPLV in α2. Classification performance based on the multi-layer participation coefficient, a multiplexity index reached 52% for iPLV and 52% for CorrEnv. Selected functional connections that build the multivariate connectomic biomarker in the edge-weighted scenario are located in default-mode, fronto-parietal, and cingulo-opercular network. Our analysis supports the notion of analyzing FCG simultaneously in intra and cross-frequency whole brain interactions with various connectivity estimators in beamformed recordings.

Improving the Reliability of Network Metrics in Structural Brain Networks by Integrating Different Network Weighting Strategies into a Single Graph

Structural brain networks estimated from diffusion MRI (dMRI) via tractography have been widely studied in healthy controls and patients with neurological and psychiatric diseases. However, few studies have addressed the reliability of derived network metrics both node-specific and network-wide. Different network weighting strategies (NWS) can be adopted to weight the strength of connection between two nodes yielding structural brain networks that are almost fully-weighted. Here, we scanned five healthy participants five times each, using a diffusion-weighted MRI protocol and computed edges between 90 regions of interest (ROI) from the Automated Anatomical Labeling (AAL) template. The edges were weighted according to nine different methods. We propose a linear combination of these nine NWS into a single graph using an appropriate diffusion distance metric. We refer to the resulting weighted graph as an Integrated Weighted Structural Brain Network (ISWBN). Additionally, we consider a topological filtering scheme that maximizes the information flow in the brain network under the constraint of the overall cost of the surviving connections. We compared each of the nine NWS and the ISWBN based on the improvement of: (a) intra-class correlation coefficient (ICC) of well-known network metrics, both node-wise and per network level; and (b) the recognition accuracy of each subject compared to the remainder of the cohort, as an attempt to access the uniqueness of the structural brain network for each subject, after first applying our proposed topological filtering scheme. Based on a threshold where the network level ICC should be >0.90, our findings revealed that six out of nine NWS lead to unreliable results at the network level, while all nine NWS were unreliable at the node level. In comparison, our proposed ISWBN performed as well as the best performing individual NWS at the network level, and the ICC was higher compared to all individual NWS at the node level. Importantly, both network and node-wise ICCs of network metrics derived from the topologically filtered ISBWN (ISWBN^TF), were further improved compared to the non-filtered ISWBN. Finally, in the recognition accuracy tests, we assigned each single ISWBN^TF to the correct subject. We also applied our methodology to a second dataset of diffusion-weighted MRI in healthy controls and individuals with psychotic experience. Following a binary classification scheme, the classification performance based on ISWBN^TF outperformed the nine different weighting strategies and the ISWBN. Overall, these findings suggest that the proposed methodology results in improved characterization of genuine between-subject differences in connectivity leading to the possibility of network-based structural phenotyping.

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 δ-γ₁, θ-γ₁, and β-γ₂ frequency pairs. Moreover, rich-club hubs in mTBI patients were overrepresented in right frontal brain areas, from θ to γ₁ frequencies, and underrepresented in left occipital regions in the δ-β, δ-γ₁, θ-β, and β-γ₂ 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.

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.

Causal Interactions between Frontal(θ) - Parieto-Occipital(α2) Predict Performance on a Mental Arithmetic Task

Many neuroimaging studies have demonstrated the different functional contributions of spatially distinct brain areas to working memory (WM) subsystems in cognitive tasks that demand both local information processing and interregional coordination. In WM cognitive task paradigms employing electroencephalography (EEG), brain rhythms such as θ and α have been linked to specific functional roles over given brain areas, but their functional coupling has not been extensively studied. Here we analyzed an arithmetic task with five cognitive workload levels (CWLs) and demonstrated functional/effective coupling between the two WM subsystems: the central executive located over frontal (F) brain areas that oscillates on the dominant θ rhythm (Frontal^θ/F^θ) and the storage buffer located over parieto-occipital (PO) brain areas that operates on the α₂ dominant brain rhythm (Parieto-Occipital^α2/PO^α2). We focused on important differences between and within WM subsystems in relation to behavioral performance. A repertoire of brain connectivity estimators was employed to elucidate the distinct roles of amplitude, phase within and between frequencies, and the hierarchical role of functionally specialized brain areas related to the task. Specifically, for each CWL, we conducted a) a conventional signal power analysis within both frequency bands at F^θ and PO^α2, b) the intra- and inter-frequency phase interactions between F^θ and PO^α2, and c) their causal phase and amplitude relationship. We found no significant statistical difference of signal power or phase interactions between correct and wrong answers. Interestingly, the study of causal interactions between F^θ and PO^α2 revealed frontal brain region(s) as the leader, while the strength differentiated between correct and wrong responses in every CWL with absolute accuracy. Additionally, zero time-lag between bilateral F^θ and right PO^a2 could serve as an indicator of mental calculation failure. Overall, our study highlights the significant role of coordinated activity between F^θ and PO^α2 via their causal interactions and the timing for arithmetic performance.

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.

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.

Altered cross-frequency coupling in resting-state MEG after mild traumatic brain injury

Cross-frequency coupling (CFC) is thought to represent a basic mechanism of functional integration of neural networks across distant brain regions. In this study, we analyzed CFC profiles from resting state Magnetoencephalographic (MEG) recordings obtained from 30 mild traumatic brain injury (mTBI) patients and 50 controls. We used mutual information (MI) to quantify the phase-to-amplitude coupling (PAC) of activity among the recording sensors in six nonoverlapping frequency bands. After forming the CFC-based functional connectivity graphs, we employed a tensor representation and tensor subspace analysis to identify the optimal set of features for subject classification as mTBI or control. Our results showed that controls formed a dense network of stronger local and global connections indicating higher functional integration compared to mTBI patients. Furthermore, mTBI patients could be separated from controls with more than 90% classification accuracy. These findings indicate that analysis of brain networks computed from resting-state MEG with PAC and tensorial representation of connectivity profiles may provide a valuable biomarker for the diagnosis of mTBI.