A Time-Varying Connectivity Analysis from Distributed EEG Sources: A Simulation Study

The strength of anticipated distractors shapes EEG alpha and theta oscillations in a Working Memory task

Working Memory (WM) requires maintenance of task-relevant information and suppression of task-irrelevant/distracting information. Alpha and theta oscillations have been extensively investigated in relation to WM. However, studies that examine both theta and alpha bands in relation to distractors, encompassing not only power modulation but also connectivity modulation, remain scarce. Here, we depicted, at the EEG-source level, the increase in power and connectivity in theta and alpha bands induced by strong relative to weak distractors during a visual Sternberg-like WM task involving the encoding of verbal items. During retention, a strong or weak distractor was presented, predictable in time and nature. Analysis focused on the encoding and retention phases before distractor presentation. Theta and alpha power were computed in cortical regions of interest, and connectivity networks estimated via spectral Granger causality and synthetized using in/out degree indices. The following modulations were observed for strong vs. weak distractors. In theta band during encoding, the power in frontal regions increased, together with frontal-to-frontal and bottom-up occipital-to-temporal-to-frontal connectivity; even during retention, bottom-up theta connectivity increased. In alpha band during retention, but not during encoding, the power in temporal-occipital regions increased, together with top-down frontal-to-occipital and temporal-to-occipital connectivity. From our results, we postulate a proactive cooperation between theta and alpha mechanisms: the first would mediate enhancement of target representation both during encoding and retention, and the second would mediate increased inhibition of sensory areas during retention only, to suppress the processing of imminent distractor without interfering with the processing of ongoing target stimulus during encoding.

Effective Connectivity Network of Aberrant Prediction Error Processing in Auditory Phantom Perception

Introduction: Prediction error (PE) is key to perception in the predictive coding framework. However, previous studies indicated the varied neural activities evoked by PE in tinnitus patients. Here, we aimed to reconcile the conflict by (1) a more nuanced view of PE, which could be driven by changing stimulus (stimulus-driven PE [sPE]) and violation of current context (context-driven PE [cPE]) and (2) investigating the aberrant connectivity networks that are engaged in the processing of the two types of PEs in tinnitus patients. Methods: Ten tinnitus patients with normal hearing and healthy controls were recruited, and a local-global auditory oddball paradigm was applied to measure the electroencephalographic difference between the two groups during sPE and cPE conditions. Results: Overall, the sPE condition engaged bottom-up and top-down connections, whereas the cPE condition engaged mostly top-down connections. The tinnitus group showed decreased sensitivity to the sPE and increased sensitivity to the cPE condition. Particularly, the auditory cortex and posterior cingulate cortex were the hubs for processing cPE in the control and tinnitus groups, respectively, showing the orientation to an internal state in tinnitus. Furthermore, tinnitus patients showed stronger connectivity to the parahippocampus and pregenual anterior cingulate cortex for the establishment of the prediction during the cPE condition. Conclusion: These results begin to dissect the role of changes in stimulus characteristics versus changes in the context of processing the same stimulus in mechanisms of tinnitus generation. Impact Statement This study delves into the number dynamics of prediction error (PE) in tinnitus, proposing a dual framework distinguishing between stimulus-driven PE (sPE) and context-driven PE (cPE). Electroencephalographic data from tinnitus patients and controls revealed distinct connectivity patterns during sPE and cPE conditions. Tinnitus patients exhibited reduced sensitivity to sPE and increased sensitivity to cPE. The auditory cortex and posterior cingulate cortex emerged as pivotal regions for cPE processing in controls and tinnitus patients, indicative of an internal state orientation in tinnitus. Enhanced connectivity to the parahippocampus and pregenual anterior cingulate cortex underscores the role of context in tinnitus pathophysiology.

An evaluation of inter-brain EEG coupling methods in hyperscanning studies

EEG-based hyperscanning technology has been increasingly applied to analyze interpersonal interactions in social neuroscience in recent years. However, different methods are employed in various of studies without a complete investigation of the suitability of these methods. Our study aimed to systematically compare typical inter-brain EEG coupling methods, with simulated EEG data generated by real EEG data. In particular, two critical metrics of noise level and time delay were manipulated, and three different coupling models were tested. The results revealed that: (1) under certain conditions, various methods were leveraged by noise level and time delay, leading to different performances; (2) most algorithms achieved better experimental results and performance under high coupling degree; (3) with our simulation process, temporal and spectral models showed relatively good results, while data simulated with phase coupling model performed worse. This is the first systematic comparison of typical inter-brain EEG coupling methods, with simulated EEG data generated by real EEG data from different subjects. Existing methods mainly focused on intra-brain coupling. To our knowledge, there was only one previous study that compared five inter-brain EEG coupling methods (Burgess in Front Human Neurosci 7:881, 2013). However, the simulated data used in this study were generated time series with varied degrees of phase coupling without considering any EEG characteristics. For future research, appropriate methods need to be selected based on possible underlying mechanisms (temporal, spectral and phase coupling model hypothesis) of a specific study, as well as the expected coupling degree and conditions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11571-022-09911-1.

Identifying good practices for detecting inter-regional linear functional connectivity from EEG

Aggregating voxel-level statistical dependencies between multivariate time series is an important intermediate step when characterising functional connectivity (FC) between larger brain regions. However, there are numerous ways in which voxel-level data can be aggregated into inter-regional FC, and the advantages of each of these approaches are currently unclear. In this study we generate ground-truth data and compare the performances of various pipelines that estimate directed and undirected linear phase-to-phase FC between regions. We test the ability of several existing and novel FC analysis pipelines to identify the true regions within which connectivity was simulated. We test various inverse modelling algorithms, strategies to aggregate time series within regions, and connectivity metrics. Furthermore, we investigate the influence of the number of interactions, the signal-to-noise ratio, the noise mix, the interaction time delay, and the number of active sources per region on the ability of detecting phase-to-phase FC. Throughout all simulated scenarios, lowest performance is obtained with pipelines involving the absolute value of coherency. Further, the combination of dynamic imaging of coherent sources (DICS) beamforming with directed FC metrics that aggregate information across multiple frequencies leads to unsatisfactory results. Pipelines that show promising results with our simulated pseudo-EEG data involve the following steps: (1) Source projection using the linearly-constrained minimum variance (LCMV) beamformer. (2) Principal component analysis (PCA) using the same fixed number of components within every region. (3) Calculation of the multivariate interaction measure (MIM) for every region pair to assess undirected phase-to-phase FC, or calculation of time-reversed Granger Causality (TRGC) to assess directed phase-to-phase FC. We formulate recommendations based on these results that may increase the validity of future experimental connectivity studies. We further introduce the free ROIconnect plugin for the EEGLAB toolbox that includes the recommended methods and pipelines that are presented here. We show an exemplary application of the best performing pipeline to the analysis of EEG data recorded during motor imagery.

Modulations of Cortical Power and Connectivity in Alpha and Beta Bands during the Preparation of Reaching Movements

Planning goal-directed movements towards different targets is at the basis of common daily activities (e.g., reaching), involving visual, visuomotor, and sensorimotor brain areas. Alpha (8-13 Hz) and beta (13-30 Hz) oscillations are modulated during movement preparation and are implicated in correct motor functioning. However, how brain regions activate and interact during reaching tasks and how brain rhythms are functionally involved in these interactions is still limitedly explored. Here, alpha and beta brain activity and connectivity during reaching preparation are investigated at EEG-source level, considering a network of task-related cortical areas. Sixty-channel EEG was recorded from 20 healthy participants during a delayed center-out reaching task and projected to the cortex to extract the activity of 8 cortical regions per hemisphere (2 occipital, 2 parietal, 3 peri-central, 1 frontal). Then, we analyzed event-related spectral perturbations and directed connectivity, computed via spectral Granger causality and summarized using graph theory centrality indices (in degree, out degree). Results suggest that alpha and beta oscillations are functionally involved in the preparation of reaching in different ways, with the former mediating the inhibition of the ipsilateral sensorimotor areas and disinhibition of visual areas, and the latter coordinating disinhibition of the contralateral sensorimotor and visuomotor areas.

Inhibitory and excitatory responses in the dorso-medial prefrontal cortex during threat processing

Objective

To evaluate cortical excitability during instructed threat processing.

Methods

Single and paired transcranial magnetic stimulation (TMS) pulses were applied to the right dorsomedial prefrontal cortex (dmPFC) during high-density electroencephalography (EEG) recording in young healthy participants (n = 17) performing an instructed threat paradigm in which one of two conditioned stimuli (CS+ but not CS-) was paired with an electric shock (unconditioned stimulus [US]). We assessed TMS-induced EEG responses with spectral power (both at electrode and source level) and information flow (effective connectivity) using Time-resolved Partial Directed Coherence (TPDC). Support vector regression (SVR) was used to predict behavioral fear ratings for CS+ based on TMS impact on excitability.

Results

During intracortical facilitation (ICF), frontal lobe theta power was enhanced for CS+ compared to single pulse TMS for the time window 0-0.5 s after TMS pulse onset (t(16) = 3.9, p < 0.05). At source level, ICF led to an increase and short intracortical inhibition (SICI) to a decrease of theta power in the bilateral dmPFC, relative to single pulse TMS during 0-0.5 s. Compared to single pulse TMS, ICF increased information flows, whereas SICI reduced the information flows in theta band between dmPFC, amygdala, and hippocampus (all at p < 0.05). The magnitude of information flows between dmPFC to amygdala and dmPFC to hippocampus during ICF (0-0.5 s), predicted individual behavioral fear ratings (CS+; coefficient above 0.75).

Conclusion

Distinct excitatory and inhibitory mechanisms take place in the dmPFC. These findings may facilitate future research attempting to investigate inhibitory/facilitatory mechanisms alterations in psychiatric disorders and their behavioral correlates.

Age-related hearing loss is associated with alterations in temporal envelope processing in different neural generators along the auditory pathway

People with age-related hearing loss suffer from speech understanding difficulties, even after correcting for differences in hearing audibility. These problems are not only attributed to deficits in audibility but are also associated with changes in central temporal processing. The goal of this study is to obtain an understanding of potential alterations in temporal envelope processing for middle-aged and older persons with and without hearing impairment. The time series of activity of subcortical and cortical neural generators was reconstructed using a minimum-norm imaging technique. This novel technique allows for reconstructing a wide range of neural generators with minimal prior assumptions regarding the number and location of the generators. The results indicated that the response strength and phase coherence of middle-aged participants with hearing impairment (HI) were larger than for normal-hearing (NH) ones. In contrast, for the older participants, a significantly smaller response strength and phase coherence were observed in the participants with HI than the NH ones for most modulation frequencies. Hemispheric asymmetry in the response strength was also altered in middle-aged and older participants with hearing impairment and showed asymmetry toward the right hemisphere. Our brain source analyses show that age-related hearing loss is accompanied by changes in the temporal envelope processing, although the nature of these changes varies with age.

Dynamic Connectivity Analysis Using Adaptive Window Size

In this paper, we propose a new method to study and evaluate the time-varying brain network dynamics. The proposed RICI-imCPCC method (relative intersection of confidence intervals for the imaginary component of the complex Pearson correlation coefficient) is based on an adaptive window size and the imaginary part of the complex Pearson correlation coefficient. It reduces the weaknesses of the existing method of constant sliding window analysis with narrow and wide windows. These are the low temporal precision and low reliability for short connectivity periods for wide windows, and high susceptibility to noise for narrow windows, all resulting in low estimation accuracy. The proposed method overcomes these shortcomings by dynamically adjusting the window width using the RICI rule, which is based on the statistical properties of the area around the observed sample. In this paper, we compare the RICI-imCPCC with the existing constant sliding window analysis method and describe its advantages. First, the mathematical principles are established. Then, the comparison between the existing and the proposed method using synthetic and real electroencephalography (EEG) data is presented. The results show that the proposed RICI-imCPCC method has improved temporal resolution and estimation accuracy compared to the existing method and is less affected by the noise. The estimation error energy calculated for the RICI-imCPCC method on synthetic signals was lower by a factor of 1.22 compared to the error of the constant sliding window analysis using narrow window size imCPCC, by a factor of 2.87 compared to using wide window size imCPCC, by a factor of 6.69 compared to using narrow window size wPLI, and by a factor of 4.72 compared to using wide window size wPLI. Analysis of the real signals shows the ability of the proposed method to detect a P300 response and to detect a decrease in dynamic connectivity due to desynchronization and blockage of mu-rhythms.

Bottom-Up and Top-Down Attention Impairment Induced by Long-Term Exposure to Noise in the Absence of Threshold Shifts

Objective

We aimed to assess the effect of noise exposure on bottom-up and top-down attention functions in industrial workers based on behavioral and brain responses recorded by the multichannel electroencephalogram (EEG).

Method

In this cross-sectional study, 563 shipyard noise-exposed workers with clinical normal hearing were recruited for cognitive testing. Personal cumulative noise exposure (CNE) was calculated with the long-term equivalent noise level and employment duration. The performance of cognitive tests was compared between the high CNE group (H-CNE, >92.2) and the low CNE group; additionally, brain responses were recorded with a 256-channel EEG from a subgroup of 20 noise-exposed (NG) workers, who were selected from the cohort with a pure tone threshold <25 dB HL from 0.25 to 16 kHz and 20 healthy controls matched for age, sex, and education. P300 and mismatch negativity (MMN) evoked by auditory stimuli were obtained to evaluate the top-down and bottom-up attention functions. The sources of P300 and MMN were investigated using GeoSource.

Results

The total score of the cognitive test (24.55 ± 3.71 vs. 25.32 ± 2.62, p < 0.01) and the subscale of attention score (5.43 ± 1.02 vs. 5.62 ± 0.67, p < 0.001) were significantly lower in the H-CNE group than in the L-CNE group. The attention score has the fastest decline of all the cognitive domain dimensions (slope = -0.03 in individuals under 40 years old, p < 0.001; slope = -0.06 in individuals older than 40 years old, p < 0.001). When NG was compared with controls, the P300 amplitude was significantly decreased in NG at Cz (3.9 ± 2.1 vs. 6.7 ± 2.3 μV, p < 0.001). In addition, the latency of P300 (390.7 ± 12.1 vs. 369.4 ± 7.5 ms, p < 0.001) and MMN (172.8 ± 15.5 vs. 157.8 ± 10.5 ms, p < 0.01) was significantly prolonged in NG compared with controls. The source for MMN for controls was in the left BA11, whereas the noise exposure group's source was lateralized to the BA20.

Conclusion

Long-term exposure to noise deteriorated the bottom-up and top-down attention functions even in the absence of threshold shifts, as evidenced by behavioral and brain responses.

Brain functional and effective connectivity based on electroencephalography recordings: A review

Functional connectivity and effective connectivity of the human brain, representing statistical dependence and directed information flow between cortical regions, significantly contribute to the study of the intrinsic brain network and its functional mechanism. Many recent studies on electroencephalography (EEG) have been focusing on modeling and estimating brain connectivity due to increasing evidence that it can help better understand various brain neurological conditions. However, there is a lack of a comprehensive updated review on studies of EEG-based brain connectivity, particularly on visualization options and associated machine learning applications, aiming to translate those techniques into useful clinical tools. This article reviews EEG-based functional and effective connectivity studies undertaken over the last few years, in terms of estimation, visualization, and applications associated with machine learning classifiers. Methods are explored and discussed from various dimensions, such as either linear or nonlinear, parametric or nonparametric, time-based, and frequency-based or time-frequency-based. Then it is followed by a novel review of brain connectivity visualization methods, grouped by Heat Map, data statistics, and Head Map, aiming to explore the variation of connectivity across different brain regions. Finally, the current challenges of related research and a roadmap for future related research are presented.

A survey of brain network analysis by electroencephalographic signals

Brain network analysis is one efficient tool in exploring human brain diseases and can differentiate the alterations from comparative networks. The alterations account for time, mental states, tasks, individuals, and so forth. Furthermore, the changes determine the segregation and integration of functional networks that lead to network reorganization (or reconfiguration) to extend the neuroplasticity of the brain. Exploring related brain networks should be of interest that may provide roadmaps for brain research and clinical diagnosis. Recent electroencephalogram (EEG) studies have revealed the secrets of the brain networks and diseases (or disorders) within and between subjects and have provided instructive and promising suggestions and methods. This review summarized the corresponding algorithms that had been used to construct functional or effective networks on the scalp and cerebral cortex. We reviewed EEG network analysis that unveils more cognitive functions and neural disorders of the human and then explored the relationship between brain science and artificial intelligence which may fuel each other to accelerate their advances, and also discussed some innovations and future challenges in the end.

Stimulus evoked causality estimation in stereo-EEG

Objective.Stereo-electroencephalography (SEEG) has recently gained importance in analyzing brain functions. Its high temporal resolution and spatial specificity make it a powerful tool to investigate the strength, direction, and spectral content of brain networks interactions, especially when these connections are stimulus-evoked. However, choosing the best approach to evaluate the flow of information is not trivial, due to the lack of validated methods explicitly designed for SEEG.Approach.We propose a novel non-parametric statistical test for event-related causality (ERC) assessment on SEEG recordings. Here, we refer to the ERC as the causality evoked by a particular part of the stimulus (a response window (RW)). We also present a data surrogation method to evaluate the performance of a causality estimation algorithm. We finally validated our pipeline using surrogate SEEG data derived from an experimentally collected dataset, and compared the most used and successful measures to estimate effective connectivity, belonging to the Geweke-Granger causality framework.Main results.Here we show that our workflow correctly identified all the directed connections in the RW of the surrogate data and prove the robustness of the procedure against synthetic noise with amplitude exceeding physiological-plausible values. Among the causality measures tested, partial directed coherence performed best.Significance.This is the first non-parametric statistical test for ERC estimation explicitly designed for SEEG datasets. The pipeline, in principle, can also be applied to the analysis of any type of time-varying estimator, if there exists a clearly defined RW.

Changes in Brain Electrical Activity According to Post-traumatic Stress Symptoms in Survivors of the Sewol Ferry Disaster: A 1-year Longitudinal Study

Objective

The pathology of post-traumatic stress disorder (PTSD) is associated with changes in brain structure and function, especially in the amygdala, medial prefrontal cortex, hippocampus, and insula. Survivors of tragic accidents often experience psychological stress and develop post-traumatic stress symptoms (PTSS), regardless of the diagnosis of PTSD. This study aimed to evaluate electroencephalographic changes according to PTSS in victims of a single traumatic event.

Methods

This study enrolled 60 survivors of the Sewol ferry disaster that occurred in 2014 from Danwon High School and collected electroencephalographic data through 19 channels twice for each person in 2014 and 2015 (mean 451.88 [standard deviation 25.77] days of follow-up). PTSS was assessed using the PTSD Checklist-Civilian Version (PCL-C) and the participants were divided into two groups according to the differences in PCL-C scores between 2014 and 2015. Electroencephalographic data were converted to three-dimensional data to perform low-resolution electrical tomographic analysis.

Results

Significant electroencephalographic changes over time were observed. The group of participants with worsened PCL-C score showed an increased change of delta slow waves in Brodmann areas 13 and 44, with the largest difference in the insula region, compared to those with improved PCL-C scores.

Conclusion

Our findings suggests that the electrophysiological changes in the insula are associated with PTSS changes.

SEED-G: Simulated EEG Data Generator for Testing Connectivity Algorithms

EEG signals are widely used to estimate brain circuits associated with specific tasks and cognitive processes. The testing of connectivity estimators is still an open issue because of the lack of a ground-truth in real data. Existing solutions such as the generation of simulated data based on a manually imposed connectivity pattern or mass oscillators can model only a few real cases with limited number of signals and spectral properties that do not reflect those of real brain activity. Furthermore, the generation of time series reproducing non-ideal and non-stationary ground-truth models is still missing. In this work, we present the SEED-G toolbox for the generation of pseudo-EEG data with imposed connectivity patterns, overcoming the existing limitations and enabling control of several parameters for data simulation according to the user's needs. We first described the toolbox including guidelines for its correct use and then we tested its performances showing how, in a wide range of conditions, datasets composed by up to 60 time series were successfully generated in less than 5 s and with spectral features similar to real data. Then, SEED-G is employed for studying the effect of inter-trial variability Partial Directed Coherence (PDC) estimates, confirming its robustness.

Brain mapping of auditory steady-state responses: A broad view of cortical and subcortical sources

Auditory steady-state responses (ASSRs) are evoked brain responses to modulated or repetitive acoustic stimuli. Investigating the underlying neural generators of ASSRs is important to gain in-depth insight into the mechanisms of auditory temporal processing. The aim of this study is to reconstruct an extensive range of neural generators, that is, cortical and subcortical, as well as primary and non-primary ones. This extensive overview of neural generators provides an appropriate basis for studying functional connectivity. To this end, a minimum-norm imaging (MNI) technique is employed. We also present a novel extension to MNI which facilitates source analysis by quantifying the ASSR for each dipole. Results demonstrate that the proposed MNI approach is successful in reconstructing sources located both within (primary) and outside (non-primary) of the auditory cortex (AC). Primary sources are detected in different stimulation conditions (four modulation frequencies and two sides of stimulation), thereby demonstrating the robustness of the approach. This study is one of the first investigations to identify non-primary sources. Moreover, we show that the MNI approach is also capable of reconstructing the subcortical activities of ASSRs. Finally, the results obtained using the MNI approach outperform the group-independent component analysis method on the same data, in terms of detection of sources in the AC, reconstructing the subcortical activities and reducing computational load.

SimMEEG software for simulating event-related MEG and EEG data with underlying functional connectivity

BACKGROUND

Simulating electromagnetic brain signals is becoming more and more prevalent because of the need to verify and validate electrophysiological analysis methods, especially with regards to measurements of functional connectivity and phase-amplitude coupling (PAC).

NEW METHOD

A Matlab software package with a graphical user interface, called SimMEEG, was developed so that novice and expert users could input basic signal parameters (e.g., multiple frequencies, amplitudes, durations, phase interactions, and phase-amplitude couplings) in order to simulate signals that approximate event-related brain activity. A major objective for developing this software was to create software to simulate interacting brain signals that could be used to verify functional connectivity analyses, specifically phase-locking value (PLV), phase-lag index (PLI) algorithms, and PAC. SimMEEG simulates three sources with various parameters defined by users. Users can define, multiple frequencies, event-related changes in amplitude, synchronizations (i.e., PLV and PLI), and PAC.

RESULTS

SimMEEG simulated signals that estimated the signal input parameters within specific target limits for PLV, PLI, and PAC. In addition, SimMEEG's synthetic signals closely approximated real human visual event-related responses. The Monte-Carlo module was used to simulate ten MEG and EEG datasets, conduct function connectivity analyses for comparison among several inverse solutions, and generate errors from ground truth. SimMEEG identified inverse solutions that had good and poor performances for functional connectivity analyses.

COMPARISON WITH EXISTING METHODS

Comparison between other simulation methods and SimMEEG are discussed.

CONCLUSION

SimMEEG software could be used to simulate real human EEG and MEG to aid in verification of analysis pipelines by comparing analyzed results to the simulated ground truth of the user-defined event-related brain signals. A core program within SimMEEG, called SimSignals, might also be used in other fields that would benefit from simulations of three interacting signals.

Neural Generators Underlying Temporal Envelope Processing Show Altered Responses and Hemispheric Asymmetry Across Age

Speech understanding problems are highly prevalent in the aging population, even when hearing sensitivity is clinically normal. These difficulties are attributed to changes in central temporal processing with age and can potentially be captured by age-related changes in neural generators. The aim of this study is to investigate age-related changes in a wide range of neural generators during temporal processing in middle-aged and older persons with normal audiometric thresholds. A minimum-norm imaging technique is employed to reconstruct cortical and subcortical neural generators of temporal processing for different acoustic modulations. The results indicate that for relatively slow modulations (<50 Hz), the response strength of neural sources is higher in older adults than in younger ones, while the phase-locking does not change. For faster modulations (80 Hz), both the response strength and the phase-locking of neural sources are reduced in older adults compared to younger ones. These age-related changes in temporal envelope processing of slow and fast acoustic modulations are possibly due to loss of functional inhibition, which is accompanied by aging. Both cortical (primary and non-primary) and subcortical neural generators demonstrate similar age-related changes in response strength and phase-locking. Hemispheric asymmetry is also altered in older adults compared to younger ones. Alterations depend on the modulation frequency and side of stimulation. The current findings at source level could have important implications for the understanding of age-related changes in auditory temporal processing and for developing advanced rehabilitation strategies to address speech understanding difficulties in the aging population.

Modeling time-varying brain networks with a self-tuning optimized Kalman filter

Brain networks are complex dynamical systems in which directed interactions between different areas evolve at the sub-second scale of sensory, cognitive and motor processes. Due to the highly non-stationary nature of neural signals and their unknown noise components, however, modeling dynamic brain networks has remained one of the major challenges in contemporary neuroscience. Here, we present a new algorithm based on an innovative formulation of the Kalman filter that is optimized for tracking rapidly evolving patterns of directed functional connectivity under unknown noise conditions. The Self-Tuning Optimized Kalman filter (STOK) is a novel adaptive filter that embeds a self-tuning memory decay and a recursive regularization to guarantee high network tracking accuracy, temporal precision and robustness to noise. To validate the proposed algorithm, we performed an extensive comparison against the classical Kalman filter, in both realistic surrogate networks and real electroencephalography (EEG) data. In both simulations and real data, we show that the STOK filter estimates time-frequency patterns of directed connectivity with significantly superior performance. The advantages of the STOK filter were even clearer in real EEG data, where the algorithm recovered latent structures of dynamic connectivity from epicranial EEG recordings in rats and human visual evoked potentials, in excellent agreement with known physiology. These results establish the STOK filter as a powerful tool for modeling dynamic network structures in biological systems, with the potential to yield new insights into the rapid evolution of network states from which brain functions emerge.

During normal behavior, brains transition between functional network states several times per second. This allows humans to quickly read a sentence, and a frog to catch a fly. Understanding these fast network dynamics is fundamental to understanding how brains work, but up to now it has proven very difficult to model fast brain dynamics for various methodological reasons. To overcome these difficulties, we designed a new Kalman filter (STOK) by innovating on previous solutions from control theory and state-space modelling. We show that STOK accurately models fast network changes in simulations and real neural data, making it an essential new tool for modelling fast brain networks in the time and frequency domain.

Time-varying effective EEG source connectivity: the optimization of model parameters

Adaptive estimation methods based on general Kalman filter are powerful tools to investigate brain networks dynamics given the non-stationary nature of neural signals. These methods rely on two parameters, the model order p and adaptation constant c, which determine the resolution and smoothness of the time-varying multivariate autoregressive estimates. A sub-optimal filtering may present consistent biases in the frequency domain and temporal distortions, leading to fallacious interpretations. Thus, the performance of these methods heavily depends on the accurate choice of these two parameters in the filter design. In this work, we sought to define an objective criterion for the optimal choice of these parameters. Since residual- and information-based criteria are not guaranteed to reach an absolute minimum, we propose to study the partial derivatives of these functions to guide the choice of p and c. To validate the performance of our method, we used a dataset of human visual evoked potentials during face perception where the generation and propagation of information in the brain is well understood and a set of simulated data where the ground truth is available.

A regularized and smoothed General Linear Kalman Filter for more accurate estimation of time-varying directed connectivity

Adaptive algorithms based on the Kalman filter are valuable tools to model the dynamic and directed Granger causal interactions between neurophysiological signals simultaneously recorded from multiple cortical regions. Among these algorithms, the General Linear Kalman Filter (GLKF) has proven to be the most accurate and reliable. Here we propose a regularized and smoothed GLKF (spsm-GLKF) with ℓ1 norm penalties based on lasso or group lasso and a fixedinterval smoother. We show that the group lasso penalty promotes sparse solutions by shrinking spurious connections to zero, while the smoothing increases the robustness of the estimates. Overall, our results demonstrate that spsm-GLKF outperforms the original GLKF, and represents a more accurate tool for the characterization of dynamical and sparse functional brain networks.

Overlapping connectivity patterns during semantic processing of abstract and concrete words revealed with multivariate Granger Causality analysis

. Abstract, unlike concrete, nouns refer to notions beyond our perception. Even though there is no consensus among linguists as to what exactly constitutes a concrete or abstract word, neuroscientists found clear evidence of a "concreteness" effect. This can, for instance, be seen in patients with language impairments due to brain injury or developmental disorder who are capable of perceiving one category better than another. Even though the results are inconclusive, neuroimaging studies on healthy subjects also provide a spatial and temporal account of differences in the processing of abstract versus concrete words. A description of the neural pathways during abstract word reading, the manner in which the connectivity patterns develop over the different stages of lexical and semantic processing compared to that of concrete word processing are still debated. We conducted a high-density EEG study on 24 healthy young volunteers using an implicit categorization task. From this, we obtained high spatio-temporal resolution data and, by means of source reconstruction, reduced the effect of signal mixing observed on scalp level. A multivariate, time-varying and directional method of analyzing connectivity based on the concept of Granger Causality (Partial Directed Coherence) revealed a dynamic network that transfers information from the right superior occipital lobe along the ventral and dorsal streams towards the anterior temporal and orbitofrontal lobes of both hemispheres. Some regions along these pathways appear to be primarily involved in either receiving or sending information. A clear difference in information transfer of abstract and concrete words was observed during the time window of semantic processing, specifically for information transferred towards the left anterior temporal lobe. Further exploratory analysis confirmed a generally stronger connectivity pattern for processing concrete words. We believe our study could guide future research towards a more refined theory of abstract word processing in the brain.

Electrophysiological Brain Connectivity: Theory and Implementation

We review the theory and algorithms of electrophysiological brain connectivity analysis. This tutorial is aimed at providing an introduction to brain functional connectivity from electrophysiological signals, including electroencephalography (EEG), magnetoencephalography (MEG), electrocorticography (ECoG), stereoelectroencephalography (SEEG). Various connectivity estimators are discussed, and algorithms introduced. Important issues for estimating and mapping brain functional connectivity with electrophysiology are discussed.