Efficient Computation of Functional Brain Networks: toward Real-Time Functional Connectivity

Gamma amplitude-envelope correlations are strongly elevated within hyperexcitable networks in focal epilepsy

Methods to quantify cortical hyperexcitability are of enormous interest for mapping epileptic networks in patients with focal epilepsy. We hypothesize that, in the resting state, cortical hyperexcitability increases firing-rate correlations between neuronal populations within seizure onset zones (SOZs). This hypothesis predicts that in the gamma frequency band (40-200 Hz), amplitude envelope correlations (AECs), a relatively straightforward measure of functional connectivity, should be elevated within SOZs compared to other areas. To test this prediction, we analyzed archived samples of interictal electrocorticographic (ECoG) signals recorded from patients who became seizure-free after surgery targeting SOZs identified by multiday intracranial recordings. We show that in the gamma band, AECs between nodes within SOZs are markedly elevated relative to those elsewhere. AEC-based node strength, eigencentrality, and clustering coefficient are also robustly increased within the SOZ with maxima in the low-gamma band (permutation test Z-scores > 8) and yield moderate discriminability of the SOZ using ROC analysis (maximal mean AUC ~ 0.73). By contrast to AECs, phase locking values (PLVs), a measure of narrow-band phase coupling across sites, and PLV-based graph metrics discriminate the seizure onset nodes weakly. Our results suggest that gamma band AECs may provide a clinically useful marker of cortical hyperexcitability in focal epilepsy.

Disentanglement of Resting State Brain Networks for Localizing Epileptogenic Zone in Focal Epilepsy

The objective of this study is to extract pathological brain networks from interictal period of E/MEG recordings to localize epileptic foci for presurgical evaluation. We proposed here a resting state E/MEG analysis framework, to disentangle brain functional networks represented by neural oscillations. By using an Embedded Hidden Markov Model, we constructed a state space for resting state recordings consisting of brain states with different spatiotemporal patterns. Functional connectivity analysis along with graph theory was applied on the extracted brain states to quantify the network features of the extracted brain states, based on which the source location of pathological states is determined. The method is evaluated by computer simulations and our simulation results revealed the proposed framework can extract brain states with high accuracy regarding both spatial and temporal profiles. We further evaluated the framework as compared with intracranial EEG defined seizure onset zone in 10 patients with drug-resistant focal epilepsy who underwent MEG recordings and were seizure free after surgical resection. The real patient data analysis showed very good localization results using the extracted pathological brain states in 6/10 patients, with localization error of about 15 mm as compared to the seizure onset zone. We show that the pathological brain networks can be disentangled from the resting-state electromagnetic recording and could be identified based on the connectivity features. The framework can serve as a useful tool in extracting brain functional networks from noninvasive resting state electromagnetic recordings, and promises to offer an alternative to aid presurgical evaluation guiding intracranial EEG electrodes implantation.

Exploring the Neural Processes behind Narrative Engagement: An EEG Study

Past cognitive neuroscience studies using naturalistic stimuli have considered narratives holistically and focused on cognitive processes. In this study, we incorporated the narrative structure, the dramatic arc, as an object of investigation, to examine how engagement levels fluctuate across a narrative-aligned dramatic arc. We explored the possibility of predicting self-reported engagement ratings from neural activity and investigated the idiosyncratic effects of each phase of the dramatic arc on brain responses as well as the relationship between engagement and brain responses. We presented a movie excerpt following the six-phase narrative arc structure to female and male participants while collecting EEG signals. We then asked this group of participants to recall the excerpt, another group to segment the video based on the dramatic arc model, and a third to rate their engagement levels while watching the movie. The results showed that the self-reported engagement ratings followed the pattern of the narrative dramatic arc. Moreover, while EEG amplitude could not predict group-averaged engagement ratings, other features comprising dynamic intersubject correlation (dISC), including certain frequency bands, dynamic functional connectivity patterns and graph features were able to achieve this. Furthermore, neural activity in the last two phases of the dramatic arc significantly predicted engagement patterns. This study is the first to explore the cognitive processes behind the dramatic arc and its phases. By demonstrating how neural activity predicts self-reported engagement, which itself aligns with the narrative structure, this study provides insights on the interrelationships between narrative structure, neural responses, and viewer engagement.

Cortico-muscular coherence and brain networks in familial adult myoclonic epilepsy and progressive myoclonic epilepsy

OBJECTIVE

Familial Adult Myoclonic Epilepsy (FAME) presents with action-activated myoclonus, often associated with epilepsy, sharing various features with Progressive Myoclonic Epilepsy (PMEs), but with slower course and limited motor disability. We aimed our study to identify measures suitable to explain the different severity of FAME2 compared to EPM1, the most common PME, and to detect the signature of the distinctive brain networks.

METHODS

We analyzed the EEG-EMG coherence (CMC) during segmental motor activity and indexes of connectivity in the two patient groups, and in healthy subjects (HS). We also investigated the regional and global properties of the network.

RESULTS

In FAME2, differently from EPM1, we found a well-localized distribution of beta-CMC and increased betweenness-centrality (BC) on the sensorimotor region contralateral to the activated hand. In both patient groups, compared to HS, there was a decline in the network connectivity indexes in the beta and gamma band, which was more obvious in FAME2.

CONCLUSIONS

In FAME2, better localized CMC and increased BC in comparison with EPM1 patients could counteract the severity and the spreading of the myoclonus. Decreased indexes of cortical integration were more severe in FAME2.

SIGNIFICANCE

Our measures correlated with different motor disabilities and identified distinctive brain network impairments.

Modeling Brain Functional Connectivity Patterns during an Isometric Arm Force Exertion Task at Different Levels of Perceived Exertion: A Graph Theoretical Approach

The perception of physical exertion is the cognitive sensation of work demands associated with voluntary muscular actions. Measurements of exerted force are crucial for avoiding the risk of overexertion and understanding human physical capability. For this purpose, various physiological measures have been used; however, the state-of-the-art in-force exertion evaluation lacks assessments of underlying neurophysiological signals. The current study applied a graph theoretical approach to investigate the topological changes in the functional brain network induced by predefined force exertion levels for twelve female participants during an isometric arm task and rated their perceived physical comfort levels. The functional connectivity under predefined force exertion levels was assessed using the coherence method for 84 anatomical brain regions of interest at the electroencephalogram (EEG) source level. Then, graph measures were calculated to quantify the network topology for two frequency bands. The results showed that high-level force exertions are associated with brain networks characterized by more significant clustering coefficients (6%), greater modularity (5%), higher global efficiency (9%), and less distance synchronization (25%) under alpha coherence. This study on the neurophysiological basis of physical exertions with various force levels suggests that brain regions communicate and cooperate higher when muscle force exertions increase to meet the demands of physically challenging tasks.

A New PLV-Spatial Filtering to Improve the Classification Performance in BCI Systems

OBJECTIVE

The performance of an EEG-based brain-computer interface (BCI) system is highly dependent on signal preprocessing. This manuscript presents a filtering method to improve the feature classification algorithms typically used in BCI.

METHODS

A graph Laplacian quadratic form using the Phase Locking Value (PLV) is applied to generate a new filtered signal in the preprocessing stage.

RESULTS

The accuracy of the classification algorithms improved significantly (up to 27.18% in the BCI Competition IV dataset, and up to 42.56% with records made with an Emotiv EPOC+). In addition, the proposed filtering algorithm has similar or better results when compared with the Filter Bank Common Spatial Pattern (FBCSP), which has disadvantages in a multiclass classification.

CONCLUSION

This paper shows how our PLV-based filtering between EEG channels could improve the performance of a BCI.

Real-Time Algorithm for Detrended Cross-Correlation Analysis of Long-Range Coupled Processes

Assessing power-law cross-correlations between a pair - or among a set - of processes is of great significance in diverse fields of analyses ranging from neuroscience to financial markets. In most cases such analyses are computationally expensive and thus carried out offline once the entire signal is obtained. However, many applications - such as mental state monitoring or financial forecasting - call for fast algorithms capable of estimating scale-free coupling in real time. Detrended cross-correlation analysis (DCCA), a generalization of the detrended fluctuation analysis (DFA) to the bivariate domain, has been introduced as a method designed to quantify power-law cross-correlations between a pair of non-stationary signals. Later, in analogy with the Pearson cross-correlation coefficient, DCCA was adapted to the detrended cross-correlation coefficient (DCCC), however as of now no online algorithms were provided for either of these analysis techniques. Here we introduce a new formula for obtaining the scaling functions in real time for DCCA. Moreover, the formula can be generalized via matrix notation to obtain the scaling relationship between not only a pair of signals, but also all possible pairs among a set of signals at the same time. This includes parallel estimation of the DFA scaling function of each individual process as well, thus allowing also for real-time acquisition of DCCC. The proposed algorithm matches its offline variants in precision, while being substantially more efficient in terms of execution time. We demonstrate that the method can be utilized for mental state monitoring on multi-channel electroencephalographic recordings obtained in eyes-closed and eyes-open resting conditions.

Cortico-muscular and cortico-cortical coherence changes resulting from Perampanel treatment in patients with cortical myoclonus

OBJECTIVE

To investigate the mechanisms by which Perampanel (PER) reduces the severity of action myoclonus, we studied on MEG signals the changes occurring in cortico-muscular coherence (CMC) and cortico-cortical connectivity in patients with progressive myoclonus epilepsies.

METHODS

The subjects performed an isometric extension of the hand; CMC and cortico-cortical connectivity were assessed using autoregressive models and generalized partial-directed coherence. The contralateral (Co) sensors showing average CMC values >0.7 of the maximum (set to 1) were grouped as central (C) regions of interest (ROI), while adjacent sensors showing CMC values >0.3 were grouped as Surrounding (Sr) ROIs.

RESULTS

Under PER treatment, CMC decreased on Co C and Sr ROIs, but also on homologous ipsilateral (Ip) ROIs; out-degrees and betweenness centrality increased in Co ROIs and decreased in Ip ROIs. The flow from Ip to Co ROIs and from activated muscles to Ip C ROI decreased.

CONCLUSION

The improvement of myoclonus corresponded to decreased CMC and recovered leadership of the cortical regions directly involved in the motor task, with a reduced interference of ipsilateral areas.

SIGNIFICANCE

Our study highlights on mechanisms suitable to treating myoclonus and suggests the role of a reduced local synchronization together a better control of distant synaptic effects.

Deep Neural Regression Prediction of Motor Imagery Skills Using EEG Functional Connectivity Indicators

Motor imaging (MI) induces recovery and neuroplasticity in neurophysical regulation. However, a non-negligible portion of users presents insufficient coordination skills of sensorimotor cortex control. Assessments of the relationship between wakefulness and tasks states are conducted to foster neurophysiological and mechanistic interpretation in MI-related applications. Thus, to understand the organization of information processing, measures of functional connectivity are used. Also, models of neural network regression prediction are becoming popular, These intend to reduce the need for extracting features manually. However, predicting MI practicing's neurophysiological inefficiency raises several problems, like enhancing network regression performance because of the overfitting risk. Here, to increase the prediction performance, we develop a deep network regression model that includes three procedures: leave-one-out cross-validation combined with Monte Carlo dropout layers, subject clustering of MI inefficiency, and transfer learning between neighboring runs. Validation is performed using functional connectivity predictors extracted from two electroencephalographic databases acquired in conditions close to real MI applications (150 users), resulting in a high prediction of pretraining desynchronization and initial training synchronization with adequate physiological interpretability.

A novel index of functional connectivity: phase lag based on Wilcoxon signed rank test

Phase synchronization has been an effective measurement of functional connectivity, detecting similar dynamics over time among distinct brain regions. However, traditional phase synchronization-based functional connectivity indices have been proved to have some drawbacks. For example, the phase locking value (PLV) index is sensitive to volume conduction, while the phase lag index (PLI) and the weighted phase lag index (wPLI) are easily affected by noise perturbations. In addition, thresholds need to be applied to these indices to obtain the binary adjacency matrix that determines the connections. However, the selection of the thresholds is generally arbitrary. To address these issues, in this paper we propose a novel index of functional connectivity, named the phase lag based on the Wilcoxon signed-rank test (PLWT). Specifically, it characterizes the functional connectivity based on the phase lag with a weighting procedure to reduce the influence of volume conduction and noise. Besides, it automatically identifies the important connections without relying on thresholds, by taking advantage of the framework of the Wilcoxon signed-rank test. The performance of the proposed PLWT index is evaluated on simulated electroencephalograph (EEG) datasets, as well as on two resting-state EEG datasets. The experimental results on the simulated EEG data show that the PLWT index is robust to volume conduction and noise. Furthermore, the brain functional networks derived by PLWT on the real EEG data exhibit a reasonable scale-free characteristic and high test-retest (TRT) reliability of graph measures. We believe that the proposed PLWT index provides a useful and reliable tool to identify the underlying neural interactions, while effectively diminishing the influence of volume conduction and noise.

Investigation of Brain Functional Networks in Children Suffering from Attention Deficit Hyperactivity Disorder

ADHD defects the recognition of facial emotions. This study assesses the neurophysiological differences between children with ADHD and matched healthy controls during a face emotional recognition task. The study also explores how brain connectivity is affected by ADHD. Electroencephalogram (EEG) signals were recorded from 64 scalp electrodes. Event-related phase coherence (ERPCOH) method was applied to pre-processed signals, and functional connectivity between any pair of electrodes was computed in different frequency bands. A logistic regression (LR) classifier with elastic net regularization (ENR) was trained to classify ADHD and HC participants using the functional connectivity of frequency bands as a potential biomarker. Subsequently, the brain network is constructed using graph-theoretic techniques, and graph indices such as clustering coefficient (C) and shortest path length (L) were calculated. Significant intra-hemispheric and the inter-hemispheric discrepancy between ADHD and healthy control (HC) groups in the beta band was observed. The graph features indicate that the clustering coefficient is significantly higher in the ADHD group than that in the HC group. At the same time, the shortest path length is significantly lower in the beta band. ADHD's brain networks have a problem in transferring information among various neural regions, which can cause a deficiency in the processing of facial emotions. The beta band seems better to reflect the differences between ADHD and HC. The observed functional connectivity and graph differences could also be helpful in ADHD investigations.

Warped phase coherence: An empirical synchronization measure combining phase and amplitude information

The entrainment between weakly coupled nonlinear oscillators, as well as between complex signals such as those representing physiological activity, is frequently assessed in terms of whether a stable relationship is detectable between the instantaneous phases extracted from the measured or simulated time-series via the analytic signal. Here, we demonstrate that adding a possibly complex constant value to this normally null-mean signal has a non-trivial warping effect. Among other consequences, this introduces a level of sensitivity to the amplitude fluctuations and average relative phase. By means of simulations of Rössler systems and experiments on single-transistor oscillator networks, it is shown that the resulting coherence measure may have an empirical value in improving the inference of the structural couplings from the dynamics. When tentatively applied to the electroencephalogram recorded while performing imaginary and real movements, this straightforward modification of the phase locking value substantially improved the classification accuracy. Hence, its possible practical relevance in brain-computer and brain-machine interfaces deserves consideration.

The blessing of Dimensionality: Feature Selection outperforms functional connectivity-based feature transformation to classify ADHD subjects from EEG patterns of phase synchronisation

Functional connectivity (FC) characterizes brain activity from a multivariate set of N brain signals by means of an NxN matrix A, whose elements estimate the dependence within each possible pair of signals. Such matrix can be used as a feature vector for (un)supervised subject classification. Yet if N is large, A is highly dimensional. Little is known on the effect that different strategies to reduce its dimensionality may have on its classification ability. Here, we apply different machine learning algorithms to classify 33 children (age [6-14 years]) into two groups (healthy controls and Attention Deficit Hyperactivity Disorder patients) using EEG FC patterns obtained from two phase synchronisation indices. We found that the classification is highly successful (around 95%) if the whole matrix A is taken into account, and the relevant features are selected using machine learning methods. However, if FC algorithms are applied instead to transform A into a lower dimensionality matrix, the classification rate drops to less than 80%. We conclude that, for the purpose of pattern classification, the relevant features should be selected among the elements of A by using appropriate machine learning algorithms.

Phase locking value revisited: teaching new tricks to an old dog

OBJECTIVE

Despite the increase in calculation power over the last few decades, the estimation of brain connectivity is still a tedious task. The high computational cost of the algorithms escalates with the square of the number of signals evaluated, usually in the range of thousands. In this work we propose a re-formulation of a widely used algorithm that allows the estimation of whole brain connectivity in much smaller times.

APPROACH

We start from the original implementation of phase locking value (PLV) and re-formulated it in a computationally very efficient way. What is more, this formulation stresses its strong similarity with coherence, which we used to introduce two new metrics insensitive to zero lag synchronization: the imaginary part of PLV (iPLV) and its corrected counterpart (ciPLV).

MAIN RESULTS

The new implementation of PLV avoids some highly CPU-expensive operations and achieves a 100-fold speedup over the original algorithm. The new derived metrics were highly robust in the presence of volume conduction. Moreover, ciPLV proved capable of ignoring zero-lag connectivity, while correctly estimating nonzero-lag connectivity.

SIGNIFICANCE

Our implementation of PLV makes it possible to calculate whole-brain connectivity in much shorter times. The results of the simulations using ciPLV suggest that this metric is ideal to measure synchronization in the presence of volume conduction or source leakage effects.

Differential Path-Length Factor's Effect on the Characterization of Brain's Hemodynamic Response Function: A Functional Near-Infrared Study

Functional near-infrared spectroscopy (fNIRS) has evolved as a neuro-imaging modality over the course of the past two decades. The removal of superfluous information accompanying the optical signal, however, remains a challenge. A comprehensive analysis of each step is necessary to ensure the extraction of actual information from measured fNIRS waveforms. A slight change in shape could alter the features required for fNIRS-BCI applications. In the present study, the effect of the differential path-length factor (DPF) values on the characteristics of the hemodynamic response function (HRF) was investigated. Results were compiled for both simulated data sets and healthy human subjects over a range of DPF values from three to eight. Different sets of activation durations and stimuli were used to generate the simulated signals for further analysis. These signals were split into optical densities under a constrained environment utilizing known values of DPF. Later, different values of DPF were used to analyze the variations of actual HRF. The results, as summarized into four categories, suggest that the DPF can change the main and post-stimuli responses in addition to other interferences. Six healthy subjects participated in this study. Their observed optical brain time-series were fed into an iterative optimization problem in order to estimate the best possible fit of HRF and physiological noises present in the measured signals with free parameters. A series of solutions was derived for different values of DPF in order to analyze the variations of HRF. It was observed that DPF change is responsible for HRF creep from actual values as well as changes in HRF characteristics.

Functional Connectivity Disruption in Subjective Cognitive Decline and Mild Cognitive Impairment: A Common Pattern of Alterations

Functional connectivity (FC) alterations represent a key feature in Alzheimer's Disease (AD) and provide a useful tool to characterize and predict the course of the disease. Those alterations have been also described in Mild Cognitive Impairment (MCI), a prodromal stage of AD. There is a growing interest in detecting AD pathology in the brain in the very early stages of the disorder. Subjective Cognitive Decline (SCD) could represent a preclinical asymptomatic stage of AD but very little is known about this population. In the present work we assessed whether FC disruptions are already present in this stage, and if they share any spatial distribution properties with MCI alterations (a condition known to be highly related to AD). To this end, we measured electromagnetic spontaneous activity with MEG in 39 healthy control elders, 41 elders with SCD and 51 MCI patients. The results showed FC alterations in both SCD and MCI compared to the healthy control group. Interestingly, both groups exhibited a very similar spatial pattern of altered links: a hyper-synchronized anterior network and a posterior network characterized by a decrease in FC. This decrease was more pronounced in the MCI group. These results highlight that elders with SCD present FC alterations. More importantly, those disruptions affected AD typically related areas and showed great overlap with the alterations exhibited by MCI patients. These results support the consideration of SCD as a preclinical stage of AD and may indicate that FC alterations appear very early in the course of the disease.