Tracking dynamic brain networks using high temporal resolution MEG measures of functional connectivity

CiftiStorm pipeline: facilitating reproducible EEG/MEG source connectomics

We present CiftiStorm, an electrophysiological source imaging (ESI) pipeline incorporating recently developed methods to improve forward and inverse solutions. The CiftiStorm pipeline produces Human Connectome Project (HCP) and megconnectome-compliant outputs from dataset inputs with varying degrees of spatial resolution. The input data can range from low-sensor-density electroencephalogram (EEG) or magnetoencephalogram (MEG) recordings without structural magnetic resonance imaging (sMRI) to high-density EEG/MEG recordings with an HCP multimodal sMRI compliant protocol. CiftiStorm introduces a numerical quality control of the lead field and geometrical corrections to the head and source models for forward modeling. For the inverse modeling, we present a Bayesian estimation of the cross-spectrum of sources based on multiple priors. We facilitate ESI in the T1w/FSAverage32k high-resolution space obtained from individual sMRI. We validate this feature by comparing CiftiStorm outputs for EEG and MRI data from the Cuban Human Brain Mapping Project (CHBMP) acquired with technologies a decade before the HCP MEG and MRI standardized dataset.

Local signal variability and functional connectivity: Sensitive measures of the excitation-inhibition ratio?

A novel network version of permutation entropy, the inverted joint permutation entropy (JPEinv), holds potential as non-invasive biomarker of abnormal excitation-inhibition (E-I) ratio in Alzheimer's disease (AD). In this computational modelling study, we test the hypotheses that this metric, and related measures of signal variability and functional connectivity, are sensitive to altered E-I ratios. The E-I ratio in each neural mass of a whole-brain computational network model was systematically varied. We evaluated whether JPEinv, local signal variability (by permutation entropy) and functional connectivity (by weighted symbolic mutual information (wsMI)) were related to E-I ratio, on whole-brain and regional level. The hub disruption index can identify regions primarily affected in terms of functional connectivity strength (or: degree) by the altered E-I ratios. Analyses were performed for a range of coupling strengths, filter and time-delay settings. On whole-brain level, higher E-I ratios were associated with higher functional connectivity (by JPEinv and wsMI) and lower local signal variability. These relationships were nonlinear and depended on the coupling strength, filter and time-delay settings. On regional level, hub-like regions showed a selective decrease in functional degree (by JPEinv and wsMI) upon a lower E-I ratio, and non-hub-like regions showed a selective increase in degree upon a higher E-I ratio. These results suggest that abnormal functional connectivity and signal variability, as previously reported in patients across the AD continuum, can inform us about altered E-I ratios.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11571-023-10003-x.

Multiplex dynamic networks in the newborn brain disclose latent links with neurobehavioral phenotypes

The higher brain functions arise from coordinated neural activity between distinct brain regions, but the spatial, temporal, and spectral complexity of these functional connectivity networks (FCNs) has challenged the identification of correlates with neurobehavioral phenotypes. Characterizing behavioral correlates of early life FCNs is important to understand the activity dependent emergence of neurodevelopmental performance and for improving health outcomes. Here, we develop an analysis pipeline for identifying multiplex dynamic FCNs that combine spectral and spatiotemporal characteristics of the newborn cortical activity. This data-driven approach automatically uncovers latent networks that show robust neurobehavioral correlations and consistent effects by in utero drug exposure. Altogether, the proposed pipeline provides a robust end-to-end solution for an objective assessment and quantitation of neurobehaviorally meaningful network constellations in the highly dynamic cortical functions.

Temporal and spatial variability of dynamic microstate brain network in disorders of consciousness

BACKGROUND

Accurately diagnosing patients with the vegetative state (VS) and the minimally conscious state (MCS) reached a misdiagnosis of approximately 40%.

METHODS

A method combined microstate and dynamic functional connectivity (dFC) to study the spatiotemporal variability of the brain in disorders of consciousness (DOC) patients was proposed. Resting-state EEG data were obtained from 16 patients with MCS and 16 patients with VS. Mutual information (MI) was used to assess the EEG connectivity in each microstate. MI-based features with statistical differences were selected as the total feature subset (TFS), then the TFS was utilized to feature selection and fed into the classifier, obtaining the optimal feature subsets (OFS) in each microstate. Subsequently, an OFS-based MI functional connectivity network (MIFCN) was constructed in the cortex.

RESULTS

The group-average MI connectivity matrix focused on all channels revealed that all five microstates exhibited stronger information interaction in the MCS when comparing with the VS. While OFS-based MIFCN, which only focused on a few channels, revealed greater MI flow in VS patients than in MCS patients under microstates A, B, C, and E, except for microstate D. Additionally, the average classification accuracy of OFS in the five microstates was 96.2%.

CONCLUSION

Constructing features based on microstates to distinguish between two categories of DOC patients had effectiveness.

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

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

Solving large-scale MEG/EEG source localisation and functional connectivity problems simultaneously using state-space models

State-space models are widely employed across various research disciplines to study unobserved dynamics. Conventional estimation techniques, such as Kalman filtering and expectation maximisation, offer valuable insights but incur high computational costs in large-scale analyses. Sparse inverse covariance estimators can mitigate these costs, but at the expense of a trade-off between enforced sparsity and increased estimation bias, necessitating careful assessment in low signal-to-noise ratio (SNR) situations. To address these challenges, we propose a three-fold solution: (1) Introducing multiple penalised state-space (MPSS) models that leverage data-driven regularisation; (2) Developing novel algorithms derived from backpropagation, gradient descent, and alternating least squares to solve MPSS models; (3) Presenting a K-fold cross-validation extension for evaluating regularisation parameters. We validate this MPSS regularisation framework through lower and more complex simulations under varying SNR conditions, including a large-scale synthetic magneto- and electro-encephalography (MEG/EEG) data analysis. In addition, we apply MPSS models to concurrently solve brain source localisation and functional connectivity problems for real event-related MEG/EEG data, encompassing thousands of sources on the cortical surface. The proposed methodology overcomes the limitations of existing approaches, such as constraints to small-scale and region-of-interest analyses. Thus, it may enable a more accurate and detailed exploration of cognitive brain functions.

Toward stability of dynamic FC estimates in neuroimaging and electrophysiology: Solutions and limits

Time-varying functional connectivity (FC) methods are used to map the spatiotemporal organization of brain activity. However, their estimation can be unstable, in the sense that different runs of the inference may yield different solutions. But to draw meaningful relations to behavior, estimates must be robust and reproducible. Here, we propose two solutions using the hidden Markov model (HMM) as a descriptive model of time-varying FC. The first, best ranked HMM, involves running the inference multiple times and selecting the best model based on a quantitative measure combining fitness and model complexity. The second, hierarchical-clustered HMM, generates stable cluster state time series by applying hierarchical clustering to the state time series obtained from multiple runs. Experimental results on fMRI and magnetoencephalography data demonstrate that these approaches substantially improve the stability of time-varying FC estimations. Overall, hierarchical-clustered HMM is preferred when the inference variability is high, while the best ranked HMM performs better otherwise.

Connectivity-based Meta-Bands: A new approach for automatic frequency band identification in connectivity analyses

The majority of electroencephalographic (EEG) and magnetoencephalographic (MEG) studies filter and analyse neural signals in specific frequency ranges, known as "canonical" frequency bands. However, this segmentation, is not exempt from limitations, mainly due to the lack of adaptation to the neural idiosyncrasies of each individual. In this study, we introduce a new data-driven method to automatically identify frequency ranges based on the topological similarity of the frequency-dependent functional neural network. The resting-state neural activity of 195 cognitively healthy subjects from three different databases (MEG: 123 subjects; EEG1: 27 subjects; EEG2: 45 subjects) was analysed. In a first step, MEG and EEG signals were filtered with a narrow-band filter bank (1 Hz bandwidth) from 1 to 70 Hz with a 0.5 Hz step. Next, the connectivity in each of these filtered signals was estimated using the orthogonalized version of the amplitude envelope correlation to obtain the frequency-dependent functional neural network. Finally, a community detection algorithm was used to identify communities in the frequency domain showing a similar network topology. We have called this approach the "Connectivity-based Meta-Bands" (CMB) algorithm. Additionally, two types of synthetic signals were used to configure the hyper-parameters of the CMB algorithm. We observed that the classical approaches to band segmentation are partially aligned with the underlying network topologies at group level for the MEG signals, but they are missing individual idiosyncrasies that may be biasing previous studies, as revealed by our methodology. On the other hand, the sensitivity of EEG signals to reflect this underlying frequency-dependent network structure is limited, revealing a simpler frequency parcellation, not aligned with that defined by the "canonical" frequency bands. To the best of our knowledge, this is the first study that proposes an unsupervised band segmentation method based on the topological similarity of functional neural network across frequencies. This methodology fully accounts for subject-specific patterns, providing more robust and personalized analyses, and paving the way for new studies focused on exploring the frequency-dependent structure of brain connectivity.

Comparing individual and group-level simulated neurophysiological brain connectivity using the Jansen and Rit neural mass model

Computational models are often used to assess how functional connectivity (FC) patterns emerge from neuronal population dynamics and anatomical brain connections. It remains unclear whether the commonly used group-averaged data can predict individual FC patterns. The Jansen and Rit neural mass model was employed, where masses were coupled using individual structural connectivity (SC). Simulated FC was correlated to individual magnetoencephalography-derived empirical FC. FC was estimated using phase-based (phase lag index (PLI), phase locking value (PLV)), and amplitude-based (amplitude envelope correlation (AEC)) metrics to analyze their goodness of fit for individual predictions. Individual FC predictions were compared against group-averaged FC predictions, and we tested whether SC of a different participant could equally well predict participants' FC patterns. The AEC provided a better match between individually simulated and empirical FC than phase-based metrics. Correlations between simulated and empirical FC were higher using individual SC compared to group-averaged SC. Using SC from other participants resulted in similar correlations between simulated and empirical FC compared to using participants' own SC. This work underlines the added value of FC simulations using individual instead of group-averaged SC for this particular computational model and could aid in a better understanding of mechanisms underlying individual functional network trajectories.

Identifying oscillatory brain networks with hidden Gaussian graphical spectral models of MEEG

Identifying the functional networks underpinning indirectly observed processes poses an inverse problem for neurosciences or other fields. A solution of such inverse problems estimates as a first step the activity emerging within functional networks from EEG or MEG data. These EEG or MEG estimates are a direct reflection of functional brain network activity with a temporal resolution that no other in vivo neuroimage may provide. A second step estimating functional connectivity from such activity pseudodata unveil the oscillatory brain networks that strongly correlate with all cognition and behavior. Simulations of such MEG or EEG inverse problem also reveal estimation errors of the functional connectivity determined by any of the state-of-the-art inverse solutions. We disclose a significant cause of estimation errors originating from misspecification of the functional network model incorporated into either inverse solution steps. We introduce the Bayesian identification of a Hidden Gaussian Graphical Spectral (HIGGS) model specifying such oscillatory brain networks model. In human EEG alpha rhythm simulations, the estimation errors measured as ROC performance do not surpass 2% in our HIGGS inverse solution and reach 20% in state-of-the-art methods. Macaque simultaneous EEG/ECoG recordings provide experimental confirmation for our results with 1/3 times larger congruence according to Riemannian distances than state-of-the-art methods.

Reward network dysfunction is associated with cognitive impairment after stroke

Stroke survivors not only suffer from severe motor, speech and neurocognitive deficits, but in many cases also from a "lack of pleasure" and a reduced motivational level. Especially apathy and anhedonic symptoms can be linked to a dysfunction of the reward system. Rewards are considered as important co-factor for learning, so the question arises as to why and how this affects the rehabilitation of stroke patients. We investigated reward behaviour, learning ability and brain network connectivity in acute (3-7d) mild to moderate stroke patients (n = 28) and age-matched healthy controls (n = 26). Reward system activity was assessed using the Monetary Incentive Delay task (MID) during magnetoencephalography (MEG). Coherence analyses were used to demonstrate reward effects on brain functional network connectivity. The MID-task showed that stroke survivors had lower reward sensitivity and required greater monetary incentives to improve performance and showed deficits in learning improvement. MEG-analyses showed a reduced network connectivity in frontal and temporoparietal regions. All three effects (reduced reward sensitivity, reduced learning ability and altered cerebral connectivity) were found to be closely related and differed strongly from the healthy group. Our results reinforce the notion that acute stroke induces reward network dysfunction, leading to functional impairment of behavioural systems. These findings are representative of a general pattern in mild strokes and are independent of the specific lesion localisation. For stroke rehabilitation, these results represent an important point to identify the reduced learning capacity after stroke and to implement individualised recovery exercises accordingly.

Breakdown of long-range spatial correlations of infraslow amplitude fluctuations of EEG oscillations in patients with current and past major depressive disorder

Introduction

Identifying biomarkers for depression from brain activity is important for the diagnosis and treatment of depression disorders. We investigated spatial correlations of the amplitude fluctuations of electroencephalography (EEG) oscillations as a potential biomarker of depression. The amplitude fluctuations of EEG oscillations intrinsically reveal both temporal and spatial correlations, indicating rapid and functional organization of the brain networks. Amid these correlations, long-range temporal correlations are reportedly impaired in patients with depression, exhibiting amplitude fluctuations closer to a random process. Based on this occurrence, we hypothesized that the spatial correlations of amplitude fluctuations would also be altered by depression.

Methods

In the present study, we extracted the amplitude fluctuations of EEG oscillations by filtering them through infraslow frequency band (0.05-0.1 Hz).

Results

We found that the amplitude fluctuations of theta oscillations during eye-closed rest depicted lower levels of spatial correlation in patients with major depressive disorder (MDD) compared to control individuals. This breakdown of spatial correlations was most prominent in the left fronto - temporal network, specifically in patients with current MDD rather than in those with past MDD. We also found that the amplitude fluctuations of alpha oscillations during eye-open rest exhibited lower levels of spatial correlation in patients with past MDD compared to control individuals or patients with current MDD.

Discussion

Our results suggest that breakdown of long-range spatial correlations may offer a biomarker for the diagnosis of depression (current MDD), as well as the tracking of the recovery from depression (past MDD).

Dissociation between phase and power correlation networks in the human brain is driven by co-occurrent bursts

Well-known haemodynamic resting-state networks are better mirrored in power correlation networks than phase coupling networks in electrophysiological data. However, what do these power correlation networks reflect? We address this long-outstanding question in neuroscience using rigorous mathematical analysis, biophysical simulations with ground truth and application of these mathematical concepts to empirical magnetoencephalography (MEG) data. Our mathematical derivations show that for two non-Gaussian electrophysiological signals, their power correlation depends on their coherence, cokurtosis and conjugate-coherence. Only coherence and cokurtosis contribute to power correlation networks in MEG data, but cokurtosis is less affected by artefactual signal leakage and better mirrors haemodynamic resting-state networks. Simulations and MEG data show that cokurtosis may reflect co-occurrent bursting events. Our findings shed light on the origin of the complementary nature of power correlation networks to phase coupling networks and suggests that the origin of resting-state networks is partly reflected in co-occurent bursts in neuronal activity.

Minimizing the distortions in electrophysiological source imaging of cortical oscillatory activity via Spectral Structured Sparse Bayesian Learning

Oscillatory processes at all spatial scales and on all frequencies underpin brain function. Electrophysiological Source Imaging (ESI) is the data-driven brain imaging modality that provides the inverse solutions to the source processes of the EEG, MEG, or ECoG data. This study aimed to carry out an ESI of the source cross-spectrum while controlling common distortions of the estimates. As with all ESI-related problems under realistic settings, the main obstacle we faced is a severely ill-conditioned and high-dimensional inverse problem. Therefore, we opted for Bayesian inverse solutions that posited a priori probabilities on the source process. Indeed, rigorously specifying both the likelihoods and a priori probabilities of the problem leads to the proper Bayesian inverse problem of cross-spectral matrices. These inverse solutions are our formal definition for cross-spectral ESI (cESI), which requires a priori of the source cross-spectrum to counter the severe ill-condition and high-dimensionality of matrices. However, inverse solutions for this problem were NP-hard to tackle or approximated within iterations with bad-conditioned matrices in the standard ESI setup. We introduce cESI with a joint a priori probability upon the source cross-spectrum to avoid these problems. cESI inverse solutions are low-dimensional ones for the set of random vector instances and not random matrices. We achieved cESI inverse solutions through the variational approximations via our Spectral Structured Sparse Bayesian Learning (ssSBL) algorithm https://github.com/CCC-members/Spectral-Structured-Sparse-Bayesian-Learning. We compared low-density EEG (10-20 system) ssSBL inverse solutions with reference cESIs for two experiments: (a) high-density MEG that were used to simulate EEG and (b) high-density macaque ECoG that were recorded simultaneously with EEG. The ssSBL resulted in two orders of magnitude with less distortion than the state-of-the-art ESI methods. Our cESI toolbox, including the ssSBL method, is available at https://github.com/CCC-members/BC-VARETA_Toolbox.

Analysis of task-related MEG functional brain networks using dynamic mode decomposition

Objective.Functional connectivity networks explain the different brain states during the diverse motor, cognitive, and sensory functions. Extracting connectivity network configurations and their temporal evolution is crucial for understanding brain function during diverse behavioral tasks.Approach.In this study, we introduce the use of dynamic mode decomposition (DMD) to extract the dynamics of brain networks. We compared DMD with principal component analysis (PCA) using real magnetoencephalography data during motor and memory tasks.Main results.The framework generates dominant connectivity brain networks and their time dynamics during simple tasks, such as button press and left-hand movement, as well as more complex tasks, such as picture naming and memory tasks. Our findings show that the proposed methodology with both the PCA-based and DMD-based approaches extracts similar dominant connectivity networks and their corresponding temporal dynamics.Significance.We believe that the proposed methodology with both the PCA and the DMD approaches has a very high potential for deciphering the spatiotemporal dynamics of electrophysiological brain network states during tasks.

The trend of disruption in the functional brain network topology of Alzheimer's disease

Alzheimer's disease (AD) is a progressive disorder associated with cognitive dysfunction that alters the brain's functional connectivity. Assessing these alterations has become a topic of increasing interest. However, a few studies have examined different stages of AD from a complex network perspective that cover different topological scales. This study used resting state fMRI data to analyze the trend of functional connectivity alterations from a cognitively normal (CN) state through early and late mild cognitive impairment (EMCI and LMCI) and to Alzheimer's disease. The analyses had been done at the local (hubs and activated links and areas), meso (clustering, assortativity, and rich-club), and global (small-world, small-worldness, and efficiency) topological scales. The results showed that the trends of changes in the topological architecture of the functional brain network were not entirely proportional to the AD progression. There were network characteristics that have changed non-linearly regarding the disease progression, especially at the earliest stage of the disease, i.e., EMCI. Further, it has been indicated that the diseased groups engaged somatomotor, frontoparietal, and default mode modules compared to the CN group. The diseased groups also shifted the functional network towards more random architecture. In the end, the methods introduced in this paper enable us to gain an extensive understanding of the pathological changes of the AD process.

Mixtures of large-scale dynamic functional brain network modes

Accurate temporal modelling of functional brain networks is essential in the quest for understanding how such networks facilitate cognition. Researchers are beginning to adopt time-varying analyses for electrophysiological data that capture highly dynamic processes on the order of milliseconds. Typically, these approaches, such as clustering of functional connectivity profiles and Hidden Markov Modelling (HMM), assume mutual exclusivity of networks over time. Whilst a powerful constraint, this assumption may be compromising the ability of these approaches to describe the data effectively. Here, we propose a new generative model for functional connectivity as a time-varying linear mixture of spatially distributed statistical "modes". The temporal evolution of this mixture is governed by a recurrent neural network, which enables the model to generate data with a rich temporal structure. We use a Bayesian framework known as amortised variational inference to learn model parameters from observed data. We call the approach DyNeMo (for Dynamic Network Modes), and show using simulations it outperforms the HMM when the assumption of mutual exclusivity is violated. In resting-state MEG, DyNeMo reveals a mixture of modes that activate on fast time scales of 100-150 ms, which is similar to state lifetimes found using an HMM. In task MEG data, DyNeMo finds modes with plausible, task-dependent evoked responses without any knowledge of the task timings. Overall, DyNeMo provides decompositions that are an approximate remapping of the HMM's while showing improvements in overall explanatory power. However, the magnitude of the improvements suggests that the HMM's assumption of mutual exclusivity can be reasonable in practice. Nonetheless, DyNeMo provides a flexible framework for implementing and assessing future modelling developments.

A Fast and Efficient Ensemble Transfer Entropy and Applications in Neural Signals

The ensemble transfer entropy (TEensemble) refers to the transfer entropy estimated from an ensemble of realizations. Due to its time-resolved analysis, it is adapted to analyze the dynamic interaction between brain regions. However, in the traditional TEensemble, multiple sets of surrogate data should be used to construct the null hypothesis distribution, which dramatically increases the computational complexity. To reduce the computational cost, a fast, efficient TEensemble with a simple statistical test method is proposed here, in which just one set of surrogate data is involved. To validate the improved efficiency, the simulated neural signals are used to compare the characteristics of the novel TEensemble with those of the traditional TEensemble. The results show that the time consumption is reduced by two or three magnitudes in the novel TEensemble. Importantly, the proposed TEensemble could accurately track the dynamic interaction process and detect the strength and the direction of interaction robustly even in the presence of moderate noises. The novel TEensemble reaches its steady state with the increased samples, which is slower than the traditional method. Furthermore, the effectiveness of the novel TEensemble was verified in the actual neural signals. Accordingly, the TEensemble proposed in this work may provide a suitable way to investigate the dynamic interactions between brain regions.

A review of MEG dynamic brain network research

The dynamic description of neural networks has attracted the attention of researchers for dynamic networks may carry more information compared with resting-state networks. As a non-invasive electrophysiological data with high temporal and spatial resolution, magnetoencephalogram (MEG) can provide rich information for the analysis of dynamic functional brain networks. In this review, the development of MEG brain network was summarized. Several analysis methods such as sliding window, Hidden Markov model, and time-frequency based methods used in MEG dynamic brain network studies were discussed. Finally, the current research about multi-modal brain network analysis and their applications with MEG neurophysiology, which are prospected to be one of the research directions in the future, were concluded.

Closed-Loop Construction and Analysis of Cortico-Muscular-Cortical Functional Network After Stroke

Brain networks allow a topological understanding into the pathophysiology of stroke-induced motor deficits, and have been an influential tool for investigating brain functions. Unfortunately, currently applied methods generally lack in the recognition of the dynamic changes in the cortical networks related to muscle activity, which is crucial to clarify the alterations of the cooperative working patterns in the motor control system after stroke. In this study, we integrate corticomuscular and intermuscular interactions to cortico-cortical network and propose a novel closed-loop construction of cortico-muscular-cortical functional network, named closed-loop network (CLN). Directional characteristic in terms of differentiating causal interactions is endowed on basis of the CLN framework, further expanding the definition of functional connectivity (FC) and effective connectivity (EC) dedicated to CLN. Next, CLN is applied to stroke patients to reveal the underlying after-effects mechanism of low frequency repetitive transcranial magnetic stimulation (rTMS) induced alterations of cortical physiologic functions during movement. Results show that the short-term modulation of rTMS is reflected in the enhancement of information interaction within the interhemispheric primary motor regions and inhibition of the coupling between motor cortex and effector muscles. CLN provides a new perspective for the study of motor-related cortical networks with muscle activities involvement instead of being restricted to brain network analysis of behaviors.

Predicting time-resolved electrophysiological brain networks from structural eigenmodes

How temporal modulations in functional interactions are shaped by the underlying anatomical connections remains an open question. Here, we analyse the role of structural eigenmodes, in the formation and dissolution of temporally evolving functional brain networks using resting-state magnetoencephalography and diffusion magnetic resonance imaging data at the individual subject level. Our results show that even at short timescales, phase and amplitude connectivity can partly be expressed by structural eigenmodes, but hardly by direct structural connections. Albeit a stronger relationship was found between structural eigenmodes and time-resolved amplitude connectivity. Time-resolved connectivity for both phase and amplitude was mostly characterised by a stationary process, superimposed with very brief periods that showed deviations from this stationary process. For these brief periods, dynamic network states were extracted that showed different expressions of eigenmodes. Furthermore, the eigenmode expression was related to overall cognitive performance and co-occurred with fluctuations in community structure of functional networks. These results implicate that ongoing time-resolved resting-state networks, even at short timescales, can to some extent be understood in terms of activation and deactivation of structural eigenmodes and that these eigenmodes play a role in the dynamic integration and segregation of information across the cortex, subserving cognitive functions.

Spatio-temporal dynamics of large-scale electrophysiological networks during cognitive action control in healthy controls and Parkinson's disease patients

Among the cognitive symptoms that are associated with Parkinson's disease (PD), alterations in cognitive action control (CAC) are commonly reported in patients. CAC enables the suppression of an automatic action, in favor of a goal-directed one. The implementation of CAC is time-resolved and arguably associated with dynamic changes in functional brain networks. However, the electrophysiological functional networks involved, their dynamic changes, and how these changes are affected by PD, still remain unknown. In this study, to address this gap of knowledge, 10 PD patients and 10 healthy controls (HC) underwent a Simon task while high-density electroencephalography (HD-EEG) was recorded. Source-level dynamic connectivity matrices were estimated using the phase-locking value in the beta (12-25 Hz) and gamma (30-45 Hz) frequency bands. Temporal independent component analyses were used as a dimension reduction tool to isolate the task-related brain network states. Typical microstate metrics were quantified to investigate the presence of these states at the subject-level. Our results first confirmed that PD patients experienced difficulties in inhibiting automatic responses during the task. At the group-level, we found three functional network states in the beta band that involved fronto-temporal, temporo-cingulate and fronto-frontal connections with typical CAC-related prefrontal and cingulate nodes (e.g., inferior frontal cortex). The presence of these networks did not differ between PD patients and HC when analyzing microstates metrics, and no robust correlations with behavior were found. In the gamma band, five networks were found, including one fronto-temporal network that was identical to the one found in the beta band. These networks also included CAC-related nodes previously identified in different neuroimaging modalities. Similarly to the beta networks, no subject-level differences were found between PD patients and HC. Interestingly, in both frequency bands, the dominant network at the subject-level was never the one that was the most durably modulated by the task. Altogether, this study identified the dynamic functional brain networks observed during CAC, but did not highlight PD-related changes in these networks that might explain behavioral changes. Although other new methods might be needed to investigate the presence of task-related networks at the subject-level, this study still highlights that task-based dynamic functional connectivity is a promising approach in understanding the cognitive dysfunctions observed in PD and beyond.

In sync with your child: The potential of parent-child electroencephalography in developmental research

Healthy interaction between parent and child is foundational for the child's socioemotional development. Recently, an innovative paradigm shift in electroencephalography (EEG) research has enabled the simultaneous measurement of neural activity in caregiver and child. This dual-EEG or hyperscanning approach, termed parent-child dual-EEG, combines the strength of both behavioral observations and EEG methods. In this review, we aim to inform on the potential of dual-EEG in parents and children (0-6 years) for developmental researchers. We first provide a general overview of the dual-EEG technique and continue by reviewing the first empirical work on the emerging field of parent-child dual-EEG, discussing the limited but fascinating findings on parent-child brain-to-behavior and brain-to-brain synchrony. We then continue by providing an overview of dual-EEG analysis techniques, including the technical challenges and solutions one may encounter. We finish by discussing the potential of parent-child dual-EEG for the future of developmental research. The analysis of multiple EEG data is technical and challenging, but when performed well, parent-child EEG may transform the way we understand how caregiver and child connect on a neurobiological level. Importantly, studying objective physiological measures of parent-child interactions could lead to the identification of novel brain-to-brain synchrony markers of interaction quality.

A mixed-modeling framework for whole-brain dynamic network analysis

The emerging area of dynamic brain network analysis has gained considerable attention in recent years. However, development of multivariate statistical frameworks that allow for examining the associations between phenotypic traits and dynamic patterns of system-level properties of the brain, and drawing statistical inference about such associations, has largely lagged behind. To address this need we developed a mixed-modeling framework that allows for assessing the relationship between any desired phenotype and dynamic patterns of whole-brain connectivity and topology. This novel framework also allows for simulating dynamic brain networks with respect to desired covariates. Unlike current tools, which largely use data-driven methods, our model-based method enables aligning neuroscientific hypotheses with the analytic approach. We demonstrate the utility of this model in identifying the relationship between fluid intelligence and dynamic brain networks by using resting-state fMRI (rfMRI) data from 200 participants in the Human Connectome Project (HCP) study. We also demonstrate the utility of this model to simulate dynamic brain networks at both group and individual levels. To our knowledge, this approach provides the first model-based statistical method for examining dynamic patterns of system-level properties of the brain and their relationships to phenotypic traits as well as simulating dynamic brain networks.

In recent years, a growing body of studies have aimed at analyzing the brain as a complex dynamic system by using various neuroimaging data. This has opened new avenues to answer compelling questions about the brain function in health and disease. However, methods that allow for providing statistical inference about how the complex interactions of the brain are associated with desired phenotypes are to be developed for a more profound insight. This study introduces a promising regression-based model to relate dynamic brain networks to desired phenotypes and provide statistical inference. Moreover, it can be used for simulating dynamic brain networks with respect to desired phenotypes at the group and individual levels.

MEG cortical microstates: spatiotemporal characteristics, dynamic functional connectivity and stimulus-evoked responses

•

EEG microstate analysis is a method to study brain states in health and disease.

•

We present a source-space microstate pipeline, giving additional anatomical insight.

•

Simulations validate source-space microstates.

•

We identify 10 resting microstates and their patterns of functional connectivity.

•

Microstate activation likelihoods change in response to auditory stimulus.

EEG microstate analysis is an approach to study brain states and their fast transitions in healthy cognition and disease. A key limitation of conventional microstate analysis is that it must be performed at the sensor level, and therefore gives limited anatomical insight. Here, we generalise the microstate methodology to be applicable to source-reconstructed electrophysiological data. Using simulations of a neural-mass network model, we first established the validity and robustness of the proposed method. Using MEG resting-state data, we uncovered ten microstates with distinct spatial distributions of cortical activation. Multivariate pattern analysis demonstrated that source-level microstates were associated with distinct functional connectivity patterns. We further demonstrated that the occurrence probability of MEG microstates were altered by auditory stimuli, exhibiting a hyperactivity of the microstate including the auditory cortex. Our results support the use of source-level microstates as a method for investigating brain dynamic activity and connectivity at the millisecond scale.

Looking through the windows: a study about the dependency of phase-coupling estimates on the data length

OBJECTIVE

Being able to characterize functional connectivity (FC) state dynamics in a real-time setting, such as in brain-computer interface, neurofeedback or closed-loop neurostimulation frameworks, requires the rapid detection of the statistical dependencies that quantify FC in short windows of data. The aim of this study is to characterize, through extensive realistic simulations, the reliability of FC estimation as a function of the data length. In particular, we focused on FC as measured by phase-coupling (PC) of neuronal oscillations, one of the most functionally relevant neural coupling modes.

APPROACH

We generated synthetic data corresponding to different scenarios by varying the data length, the signal-to-noise ratio, the phase difference value, the spectral analysis approach (Hilbert or Fourier) and the fractional bandwidth. We compared seven PC metrics, i.e. imaginary part of phase locking value (PLV), PLV of orthogonalized signals, phase lag index (PLI), debiased weighted PLI, imaginary part of coherency, coherence of orthogonalized signals and lagged coherence.

MAIN RESULTS

Our findings show that, for a signal-to-noise-ratio of at least 10 dB, a data window that contains 5 to 8 cycles of the oscillation of interest (e.g. a 500-800ms window at 10Hz) is generally required to achieve reliable PC estimates. In general, Hilbert-based approaches were associated with higher performance than Fourier-based approaches. Furthermore, the results suggest that, when the analysis is performed in a narrow frequency range, a larger window is required.

SIGNIFICANCE

The achieved results pave the way to the introduction of best-practice guidelines to be followed when a real-time frequency-specific PC assessment is at target.

Schizophrenia induces abnormal frequency-dependent patterns of dynamic brain network reconfiguration during an auditory oddball task

OBJECTIVE

Schizophrenia is a psychiatric disorder that has been shown to disturb the dynamic top-down processing of sensory information. Various imaging techniques have revealed abnormalities in brain activity associated with this disorder, both locally and between cerebral regions. However, there is increasing interest in investigating dynamic network response to novel and relevant events at the network level during an attention-demanding task with high-temporal-resolution techniques. The aim of the work was: (i) to test the capacity of a novel algorithm to detect recurrent brain meta-states from auditory oddball task recordings; and (ii) to evaluate how the dynamic activation and behavior of the aforementioned meta-states were altered in schizophrenia, since it has been shown to impair top-down processing of sensory information.

APPROACH

A novel unsupervised method for the detection of brain meta-states based on recurrence plots and community detection algorithms, previously tested on resting-state data, was used on auditory oddball task recordings. Brain meta-states and several properties related to their activation during target trials in the task were extracted from electroencephalography (EEG) data from patients with schizophrenia and cognitively healthy controls.

MAIN RESULTS

The methodology successfully detected meta-states during an auditory oddball task, and they appeared to show both frequency-dependent time-locked and non-time-locked activity with respect to the stimulus onset. Moreover, patients with schizophrenia displayed higher network diversity, and showed more sluggish meta-state transitions, reflected in increased dwell times, less complex meta-state sequences, decreased meta-state space speed, and abnormal ratio of negative meta-state correlations.

SIGNIFICANCE

Abnormal cognition in schizophrenia is also reflected in decreased brain flexibility at the dynamic network level, which may hamper top-down processing, possibly indicating impaired decision-making linked to dysfunctional predictive coding. Moreover, the results showed the ability of the methodology to find meaningful and task-relevant changes in dynamic connectivity and pathology-related group differences.

Effect of segment length, sampling frequency, and imaging modality on the estimation of measures of brain meta-state activation: an MEG/EEG study

The main objective of this study was to examine the influence that recording length, sampling frequency, and imaging modality have on the estimation and characterization of spontaneous brain meta-states during rest. To this end, a recently developed method of meta-state extraction and characterization was applied to a subset of 16 healthy elderly subjects from two independent electroencephalographic and magnetoencephalographic (EEG/MEG) databases. The recordings were segmented into the first 5, 10, 15, 20, 25, 30, 60 and 90-s of artifact-free activity and meta-states were extracted. Temporal activation sequence (TAS) complexity, which characterizes the complexity of the metastateactivation sequences during rest, was calculated. Then, its stability as a function of segment length, sampling frequency, and imaging modality was assessed. The results showed that, in general, the minimum segment length needed to fully characterize resting-state meta-state activation complexity ranged from 15 to 25 seconds. Moreover, it was found that the sampling frequency has a slight influence on the complexity measure, and that results were similar across EEG and MEG. The findings indicate that the proposed methodology can be applied to both EEG and MEG recordings and displays stable behavior with relatively short segments. However, methodological choices, such as sampling frequency, should be carefully considered.

Impaired saccadic eye movements in multiple sclerosis are related to altered functional connectivity of the oculomotor brain network

•

Impaired eye movements in multiple sclerosis (MS) and functional connectivity (FC)

•

Eye movements related to altered FC of the oculomotor brain network.

•

Lower (beta band) and higher (theta/delta band) FC related to abnormal eye movements.

•

Regional changes were more informative than whole-network measures.

•

Eye movement parameters also related to disability and cognitive dysfunction.

Background

Impaired eye movements in multiple sclerosis (MS) are common and could represent a non-invasive and accurate measure of (dys)functioning of interconnected areas within the complex brain network. The aim of this study was to test whether altered saccadic eye movements are related to changes in functional connectivity (FC) in patients with MS.

Methods

Cross-sectional eye movement (pro-saccades and anti-saccades) and magnetoencephalography (MEG) data from the Amsterdam MS cohort were included from 176 MS patients and 33 healthy controls. FC was calculated between all regions of the Brainnetome atlas in six conventional frequency bands. Cognitive function and disability were evaluated by previously validated measures. The relationships between saccadic parameters and both FC and clinical scores in MS patients were analysed using multivariate linear regression models.

Results

In MS pro- and anti-saccades were abnormal compared to healthy controls A relationship of saccadic eye movements was found with FC of the oculomotor network, which was stronger for regional than global FC. In general, abnormal eye movements were related to higher delta and theta FC but lower beta FC. Strongest associations were found for pro-saccadic latency and FC of the precuneus (beta band β = -0.23, p = .006), peak velocity and FC of the parietal eye field (theta band β = -0.25, p = .005) and gain and FC of the inferior frontal eye field (theta band β = -0.25, p = .003). Pro-saccadic latency was also strongly associated with disability scores and cognitive dysfunction.

Conclusions

Impaired saccadic eye movements were related to functional connectivity of the oculomotor network and clinical performance in MS. This study also showed that, in addition to global network connectivity, studying regional changes in MEG studies could yield stronger correlations.

Sensitivity of amplitude and phase based MEG measures of interhemispheric connectivity during unilateral finger movements

Interactions between different brain regions can be revealed by dependencies between their neuronal oscillations. We examined the sensitivity of different oscillatory connectivity measures in revealing interhemispheric interactions between primary motor cortices (M1s) during unilateral finger movements. Based on frequency, amplitude, and phase of the oscillations, a number of metrics have been developed to measure connectivity between brain regions, and each metric has its own strengths, weaknesses, and pitfalls. Taking advantage of the well-known movement-related modulations of oscillatory amplitude in M1s, this study compared and contrasted a number of leading connectivity metrics during distinct phases of oscillatory power changes. Between M1s during unilateral movements, we found that phase-based metrics were effective at revealing connectivity during the beta (15-35 Hz) rebound period linked to movement termination, but not during the early period of beta desynchronization occurring during the movement itself. Amplitude correlation metrics revealed robust connectivity during both periods. Techniques for estimating the direction of connectivity had limited success. Granger Causality was not well suited to studying these connections because it was strongly confounded by differences in signal-to-noise ratio linked to modulation of beta amplitude occurring during the task. Phase slope index was suggestive but not conclusive of a unidirectional influence between motor cortices during the beta rebound. Our findings suggest that a combination of amplitude and phase-based metrics is likely required to fully characterize connectivity during task protocols that involve modulation of oscillatory power, and that amplitude-based metrics appear to be more sensitive despite the lack of directional information.

Magnetoencephalography Signal Processing, Forward Modeling, Magnetoencephalography Inverse Source Imaging, and Coherence Analysis

Magnetoencephalography (MEG) is a noninvasive functional imaging technique for the brain. MEG directly measures the magnetic signal due to neuronal activation in gray matter with high spatial localization accuracy. The first part of this article covers the overall concepts of MEG and the forward and inverse modeling techniques. It is followed by examples of analyzing evoked and resting-state MEG signals using a high-resolution MEG source imaging technique. Next, different techniques for connectivity and network analysis are reviewed with examples showing connectivity estimates from resting-state and epileptic activity.

Abnormal meta-state activation of dynamic brain networks across the Alzheimer spectrum

The characterization of the distinct dynamic functional connectivity (dFC) patterns that activate in the brain during rest can help to understand the underlying time-varying network organization. The presence and behavior of these patterns (known as meta-states) have been widely studied by means of functional magnetic resonance imaging (fMRI). However, modalities with high-temporal resolution, such as electroencephalography (EEG), enable the characterization of fast temporally evolving meta-state sequences. Mild cognitive impairment (MCI) and dementia due to Alzheimer's disease (AD) have been shown to disrupt spatially localized activation and dFC between different brain regions, but not much is known about how they affect meta-state network topologies and their network dynamics. The main hypothesis of the study was that MCI and dementia due to AD alter normal meta-state sequences by inducing a loss of structure in their patterns and a reduction of their dynamics. Moreover, we expected that patients with MCI would display more flexible behavior compared to patients with dementia due to AD. Thus, the aim of the current study was twofold: (i) to find repeating, distinctly organized network patterns (meta-states) in neural activity; and (ii) to extract information about meta-state fluctuations and how they are influenced by MCI and dementia due to AD. To accomplish these goals, we present a novel methodology to characterize dynamic meta-states and their temporal fluctuations by capturing aspects based on both their discrete activation and the continuous evolution of their individual strength. These properties were extracted from 60-s resting-state EEG recordings from 67 patients with MCI due to AD, 50 patients with dementia due to AD, and 43 cognitively healthy controls. First, the instantaneous amplitude correlation (IAC) was used to estimate instantaneous functional connectivity with a high temporal resolution. We then extracted meta-states by means of graph community detection based on recurrence plots (RPs), both at the individual- and group-level. Subsequently, a diverse set of properties of the continuous and discrete fluctuation patterns of the meta-states was extracted and analyzed. The main novelty of the methodology lies in the usage of Louvain GJA community detection to extract meta-states from IAC-derived RPs and the extended analysis of their discrete and continuous activation. Our findings showed that distinct dynamic functional connectivity meta-states can be found on the EEG time-scale, and that these were not affected by the oscillatory slowing induced by MCI or dementia due to AD. However, both conditions displayed a loss of meta-state modularity, coupled with shorter dwell times and higher complexity of the meta-state sequences. Furthermore, we found evidence that meta-state sequencing is not entirely random; it shows an underlying structure that is partially lost in MCI and dementia due to AD. These results show evidence that AD progression is associated with alterations in meta-state switching, and a degradation of dynamic brain flexibility.

Dynamics of task-related electrophysiological networks: a benchmarking study

Motor, sensory and cognitive functions rely on dynamic reshaping of functional brain networks. Tracking these rapid changes is crucial to understand information processing in the brain, but challenging due to the great variety of dimensionality reduction methods used at the network-level and the limited evaluation studies. Using Magnetoencephalography (MEG) combined with Source Separation (SS) methods, we present an integrated framework to track fast dynamics of electrophysiological brain networks. We evaluate nine SS methods applied to three independent MEG databases (N=95) during motor and memory tasks. We report differences between these methods at the group and subject level. We seek to help researchers in choosing objectively the appropriate SS method when tracking fast reconfiguration of functional brain networks, due to its enormous benefits in cognitive and clinical neuroscience.

Brain dysconnectivity relates to disability and cognitive impairment in multiple sclerosis

The pathophysiology of cognitive dysfunction in multiple sclerosis (MS) is still unclear. This magnetoencephalography (MEG) study investigates the impact of MS on brain resting-state functional connectivity (rsFC) and its relationship to disability and cognitive impairment. We investigated rsFC based on power envelope correlation within and between different frequency bands, in a large cohort of participants consisting of 99 MS patients and 47 healthy subjects. Correlations were investigated between rsFC and outcomes on disability, disease duration and 7 neuropsychological scores within each group, while stringently correcting for multiple comparisons and possible confounding factors. Specific dysconnections correlating with MS-induced physical disability and disease duration were found within the sensorimotor and language networks, respectively. Global network-level reductions in within- and cross-network rsFC were observed in the default-mode network. Healthy subjects and patients significantly differed in their scores on cognitive fatigue and verbal fluency. Healthy subjects and patients showed different correlation patterns between rsFC and cognitive fatigue or verbal fluency, both of which involved a shift in patients from the posterior default-mode network to the language network. Introducing electrophysiological rsFC in a regression model of verbal fluency and cognitive fatigue in MS patients significantly increased the explained variance compared to a regression limited to structural MRI markers (relative thalamic volume and lesion load). This MEG study demonstrates that MS induces distinct changes in the resting-state functional brain architecture that relate to disability, disease duration and specific cognitive functioning alterations. It highlights the potential value of electrophysiological intrinsic rsFC for monitoring the cognitive impairment in patients with MS.

High-dimensional brain-wide functional connectivity mapping in magnetoencephalography

BACKGROUND

Brain functional connectivity (FC) analyses based on magneto/electroencephalography (M/EEG) signals have yet to exploit the intrinsic high-dimensional information. Typically, these analyses are constrained to regions of interest to avoid the curse of dimensionality, with the latter leading to conservative hypothesis testing.

NEW METHOD

We removed such constraint by estimating high-dimensional source-based M/EEG-FC using cluster-permutation statistic (CPS) and demonstrated the feasibility of this approach by identifying resting-state changes in mild cognitive impairment (MCI), a prodromal stage of Alzheimer's disease. Particularly, we proposed a unified framework for CPS analysis together with a novel neighbourhood measure to estimate more compact and neurophysiological plausible neural communication. As clusters could more confidently reveal interregional communication, we proposed and tested a cluster-strength index to demonstrate other advantages of CPS analysis.

RESULTS

We found clusters of increased communication or hypersynchronization in MCI compared to healthy controls in delta (1-4 Hz) and higher-theta (6-8 Hz) bands oscillations. These mainly consisted of interactions between occipitofrontal and occipitotemporal regions in the left hemisphere, which may be critically affected in the early stages of Alzheimer's disease.

CONCLUSIONS

Our approach could be important to create high-resolution FC maps from neuroimaging studies in general, allowing the multimodal analysis of neural communication across multiple spatial scales. Particularly, FC clusters more robustly represent the interregional communication by identifying dense bundles of connections that are less sensitive to inter-individual anatomical and functional variability. Overall, this approach could help to better understand neural information processing in healthy and disease conditions as needed for developing biomarker research.

Temporal Dynamics on Decoding Target Stimuli in Rapid Serial Visual Presentation using Magnetoencephalography

Rapid serial visual presentation (RSVP) is a high efficient paradigm in brain-computer interface (BCI). Target detection accuracy is the first consideration of RSVP-BCI. But the influence of different frequency bands and time ranges on decoding accuracy are still an open questions. Moreover, the underlying neural dynamic of the rapid target detecting process is still unclear. Methods: This work focused the temporal dynamic of the responses triggered by target stimuli in a static RSVP paradigm using paired structural Magnetic Resonance Imaging (MRI) and magnetoencephalography (MEG) signals with different frequency bands. Multivariate pattern analysis (MVPA) was applied on the MEG signal with different frequency bands and time points after stimuli onset. Cortical neuronal activation estimation technology was also applied to present the temporal-spatial dynamic on cortex surface. Results: The MVPA results showed that the low frequency signals (0.1 - 7 Hz) yield highest decoding accuracy, and the decoding power reached its peak at 0.4 second after target stimuli onset. The cortical neuronal activation method identified the target stimuli triggered regions, like bilateral parahippocampal cortex, precentral gyrus and insula cortex, and the averaged time series were presented.

Profound regional spectral, connectivity, and network changes reflect visual deficits in posterior cortical atrophy: an EEG study

Patients with posterior cortical atrophy (PCA-AD) show more severe visuospatial and perceptual deficits than those with typical AD (tAD). The aim of this study was to investigate whether functional alterations measured by electroencephalography can help understand the mechanisms that explain this clinical heterogeneity. 21-channel electroencephalography recordings of 29 patients with PCA-AD were compared with 29 patients with tAD and 29 controls matched for age, gender, and disease severity. Patients with PCA-AD and tAD both showed a global decrease in fast and increase in slow oscillatory activity compared with controls. This pattern was, however, more profound in patients with PCA-AD which was driven by more extensive slowing of the posterior regions. Alpha band functional connectivity showed a similar decrease in PCA-AD and tAD. Compared with controls, a less integrated network topology was observed in PCA-AD, with a decrease of posterior and an increase of frontal hubness. In PCA-AD, decreased right parietal peak frequency correlated with worse performance on visual tasks. Regional vulnerability of the posterior network might explain the atypical pattern of neurodegeneration in PCA-AD.

Connectome spectral analysis to track EEG task dynamics on a subsecond scale

We present an approach for tracking fast spatiotemporal cortical dynamics in which we combine white matter connectivity data with source-projected electroencephalographic (EEG) data. We employ the mathematical framework of graph signal processing in order to derive the Fourier modes of the brain structural connectivity graph, or "network harmonics". These network harmonics are naturally ordered by smoothness. Smoothness in this context can be understood as the amount of variation along the cortex, leading to a multi-scale representation of brain connectivity. We demonstrate that network harmonics provide a sparse representation of the EEG signal, where, at certain times, the smoothest 15 network harmonics capture 90% of the signal power. This suggests that network harmonics are functionally meaningful, which we demonstrate by using them as a basis for the functional EEG data recorded from a face detection task. There, only 13 network harmonics are sufficient to track the large-scale cortical activity during the processing of the stimuli with a 50 ​ms resolution, reproducing well-known activity in the fusiform face area as well as revealing co-activation patterns in somatosensory/motor and frontal cortices that an unconstrained ROI-by-ROI analysis fails to capture. The proposed approach is simple and fast, provides a means of integration of multimodal datasets, and is tied to a theoretical framework in mathematics and physics. Thus, network harmonics point towards promising research directions both theoretically - for example in exploring the relationship between structure and function in the brain - and practically - for example for network tracking in different tasks and groups of individuals, such as patients.

Altered Connectivity During a False-Belief Task in Adults With Autism Spectrum Disorder

BACKGROUND

Deficits in social communication are one of the main features of autism spectrum disorder (ASD). Adults with ASD show atypical brain activity during false-belief understanding, an aspect of social communication involving the ability to infer that an individual can have an incorrect belief about a situation. Our study is the first to investigate whether adults with ASD exhibit differences in frequency-specific functional connectivity patterns during false-belief reasoning.

METHODS

We used magnetoencephalography to contrast functional connectivity underlying false-belief understanding between 40 adults with ASD and 39 control adults. We examined whole-brain phase synchrony measures during a false-belief task in 3 frequency bands: theta (4-7 Hz), alpha (8-14 Hz), and beta (15-30 Hz).

RESULTS

Adults with ASD demonstrated reduced theta-band connectivity compared with control adults between several right-lateralized and midline regions such as the medial prefrontal cortex, right temporoparietal junction, right inferior frontal gyrus, and right superior temporal gyrus. During false-belief trials, they also recruited a network in the beta band that included primary visual regions such as the bilateral inferior occipital gyri and the left anterior temporoparietal junction.

CONCLUSIONS

Reduced theta-band synchrony between areas associated with mentalizing, inhibition, and visual processing implies some difficulty in communication among these functions in ASD. This impairment in top-down control in the theta band may be counterbalanced by their engagement of a beta-band network because both the left anterior temporoparietal junction and beta-band oscillations are associated with attentional processes. Thus, adults with ASD demonstrate alternative neural mechanisms for successful false-belief reasoning.

Discovering dynamic task-modulated functional networks with specific spectral modes using MEG

Efficient neuronal communication between brain regions through oscillatory synchronization at certain frequencies is necessary for cognition. Such synchronized networks are transient and dynamic, established on the timescale of milliseconds in order to support ongoing cognitive operations. However, few studies characterizing dynamic electrophysiological brain networks have simultaneously accounted for temporal non-stationarity, spectral structure, and spatial properties. Here, we propose an analysis framework for characterizing the large-scale phase-coupling network dynamics during task performance using magnetoencephalography (MEG). We exploit the high spatiotemporal resolution of MEG to measure time-frequency dynamics of connectivity between parcellated brain regions, yielding data in tensor format. We then use a tensor component analysis (TCA)-based procedure to identify the spatio-temporal-spectral modes of covariation among separate regions in the human brain. We validate our pipeline using MEG data recorded during a hand movement task, extracting a transient motor network with beta-dominant spectral mode, which is significantly modulated by the movement task. Next, we apply the proposed pipeline to explore brain networks that support cognitive operations during a working memory task. The derived results demonstrate the temporal formation and dissolution of multiple phase-coupled networks with specific spectral modes, which are associated with face recognition, vision, and movement. The proposed pipeline can characterize the spectro-temporal dynamics of functional connectivity in the brain on the subsecond timescale, commensurate with that of cognitive performance.

Whole-brain electrophysiological functional connectivity dynamics in resting-state EEG

OBJECTIVE

Functional connectivity (FC) dynamics have been studied in functional magnetic resonance imaging (fMRI) data, while it is largely unknown in electrophysiological data, e.g. EEG.

APPROACH

The present study proposed a novel analytic framework to study spatiotemporal dynamics of FC (dFC) in resting-state human EEG data, including independent component analysis, cortical source imaging, sliding-window correlation analysis, and k-means clustering.

MAIN RESULTS

Our results confirm that major fMRI intrinsic connectivity networks (ICNs) can be successfully reconstructed from EEG using our analytic framework. Prominent spatial and temporal variability were revealed in these ICNs. The mean dFC spatial patterns of individual ICNs resemble their corresponding static FC (sFC) patterns but show fewer cross-talks among distinct ICNs. Our investigation unveils evidences of time-domain variations in individual ICNs comparable to their mean FC level in terms of magnitude. The major contributors to these variations are from the frequency below 0.0156 Hz, in the similar range of FC dynamics from fMRI data. Among different ICNs, larger temporal variabilities are observed in the frontal attention and auditory/visual ICNs, while sensorimotor, salience, and default model networks showed less. Our analytic framework for the first time revealed quasi-stable states within individual EEG ICNs, with various strengths or spatial patterns that were reliably detected at both group and individual levels. These states all together reveal a more complete picture of EEG ICNs: (1) quasi-stable state spatial patterns as a whole for each EEG ICN are more consistent with the corresponding fMRI ICN in terms of the bilateral distribution and multi-nodes structure; (2) EEG ICNs reveal more transient patterns about within-ICN between-node communications than fMRI ICNs.

SIGNIFICANCE

The present findings highlight the fact that rich temporal and spatial dynamics exist in ICN that can be detected from EEG data. Future studies might extend investigations towards spectral dynamics of EEG ICNs.

Characterization of the dynamic behavior of neural activity in Alzheimer's disease: exploring the non-stationarity and recurrence structure of EEG resting-state activity

OBJECTIVE

Mild cognitive impairment (MCI) and dementia due to Alzheimer's disease (AD) have been shown to induce perturbations to normal neuronal behavior and disrupt neuronal networks. Recent work suggests that the dynamic properties of resting-state neuronal activity could be affected by MCI and AD-induced neurodegeneration. The aim of the study was to characterize these properties from different perspectives: (i) using the Kullback-Leibler divergence (KLD), a measure of non-stationarity derived from the continuous wavelet transform; and (ii) using the entropy of the recurrence point density ([Formula: see text]) and the median of the recurrence point density ([Formula: see text]), two novel metrics based on recurrence quantification analysis.

APPROACH

KLD, [Formula: see text] and [Formula: see text] were computed for 49 patients with dementia due to AD, 66 patients with MCI due to AD and 43 cognitively healthy controls from 60 s electroencephalographic (EEG) recordings with a 10 s sliding window with no overlap. Afterwards, we tested whether the measures reflected alterations to normal neuronal activity induced by MCI and AD.

MAIN RESULTS

Our results showed that frequency-dependent alterations to normal dynamic behavior can be found in patients with MCI and AD, both in non-stationarity and recurrence structure. Patients with MCI showed signs of patterns of abnormal state recurrence in the theta (4-8 Hz) and beta (13-30 Hz) frequency bands that became more marked in AD. Moreover, abnormal non-stationarity patterns were found in MCI patients, but not in patients with AD in delta (1-4 Hz), alpha (8-13 Hz), and gamma (30-70 Hz).

SIGNIFICANCE

The alterations in normal levels of non-stationarity in patients with MCI suggest an initial increase in cortical activity during the development of AD. This increase could possibly be due to an impairment in neuronal inhibition that is not present during later stages. MCI and AD induce alterations to the recurrence structure of cortical activity, suggesting that normal state switching during rest may be affected by these pathologies.

Neural circuits of idiopathic Normal Pressure Hydrocephalus: A perspective review of brain connectivity and symptoms meta-analysis

Idiopathic normal pressure hydrocephalus (iNPH) is a prevalent reversible neurological disorder characterized by impaired locomotion, cognition and urinary control with ventriculomegaly. Symptoms can be relieved with cerebrospinal fluid drainage, which makes iNPH the leading cause of reversible dementia. Because of a limited understanding of pathophysiological mechanisms, unspecific symptoms and the high prevalence of comorbidity (i.e. Alzheimer's disease), iNPH is largely underdiagnosed. For these reasons, there is an urgent need for developing noninvasive quantitative biomarkers for iNPH diagnosis and prognosis. Structural and functional changes of brain circuits in relation to symptoms and treatment response are expected to deliver major advances in this direction. We review structural and functional brain connectivity findings in iNPH and complement those findings with iNPH symptom meta-analyses in healthy populations. Our goal is to reinforce our conceptualization of iNPH as to brain network mechanisms and foster the development of new hypotheses for future research and treatment options.

Frontoparietal Network Connectivity During an N-Back Task in Adults With Autism Spectrum Disorder

BACKGROUND

Short-term and working memory (STM and WM) deficits have been demonstrated in individuals with autism spectrum disorder (ASD) and may emerge through atypical functional activity and connectivity of the frontoparietal network, which exerts top-down control necessary for successful STM and WM processes. Little is known regarding the spectral properties of the frontoparietal network during STM or WM processes in ASD, although certain neural frequencies have been linked to specific neural mechanisms.

METHODS

We analysed magnetoencephalographic data from 39 control adults (26 males; 27.15 ± 5.91 years old) and 40 adults with ASD (26 males; 27.17 ± 6.27 years old) during a 1-back condition (STM) of an n-back task, and from a subset of this sample during a 2-back condition (WM). We performed seed-based connectivity analyses using regions of the frontoparietal network. Interregional synchrony in theta, alpha, and beta bands was assessed with the phase difference derivative and compared between groups during periods of maintenance and recognition.

RESULTS

During maintenance of newly presented vs. repeated stimuli, the two groups did not differ significantly in theta, alpha, or beta phase synchrony for either condition. Adults with ASD showed alpha-band synchrony in a network containing the right dorsolateral prefrontal cortex, bilateral inferior parietal lobules (IPL), and precuneus in both 1- and 2-back tasks, whereas controls demonstrated alpha-band synchrony in a sparser set of regions, including the left insula and IPL, in only the 1-back task. During recognition of repeated vs. newly presented stimuli, adults with ASD exhibited decreased theta-band connectivity compared to controls in a network with hubs in the right inferior frontal gyrus and left IPL in the 1-back condition. Whilst there were no group differences in connectivity in the 2-back condition, adults with ASD showed no frontoparietal network recruitment during recognition, whilst controls activated networks in the theta and beta bands.

CONCLUSIONS

Our findings suggest that since adults with ASD performed well on the n-back task, their appropriate, but effortful recruitment of alpha-band mechanisms in the frontoparietal network to maintain items in STM and WM may compensate for atypical modulation of this network in the theta band to recognise previously presented items in STM.

How Sensitive Are Conventional MEG Functional Connectivity Metrics With Sliding Windows to Detect Genuine Fluctuations in Dynamic Functional Connectivity?

Despite advances in the field of dynamic connectivity, fixed sliding window approaches for the detection of fluctuations in functional connectivity are still widely used. The use of conventional connectivity metrics in conjunction with a fixed sliding window comes with the arbitrariness of the chosen window lengths. In this paper we use multivariate autoregressive and neural mass models with a priori defined ground truths to systematically analyze the sensitivity of conventional metrics in combination with different window lengths to detect genuine fluctuations in connectivity for various underlying state durations. Metrics of interest are the coherence, imaginary coherence, phase lag index, phase locking value and the amplitude envelope correlation. We performed analysis for two nodes and at the network level. We demonstrate that these metrics show indeed higher variability for genuine temporal fluctuations in connectivity compared to a static connectivity state superimposed by noise. Overall, the error of the connectivity estimates themselves decreases for longer state durations (order of seconds), while correlations of the connectivity fluctuations with the ground truth was higher for longer state durations. In general, metrics, in combination with a sliding window, perform poorly for very short state durations. Increasing the SNR of the system only leads to a moderate improvement. In addition, at the network level, only longer window widths were sufficient to detect plausible resting state networks that matched the underlying ground truth, especially for the phase locking value, amplitude envelope correlation and coherence. The length of these longer window widths did not necessarily correspond to the underlying state durations. For short window widths resting state network connectivity patterns could not be retrieved. We conclude that fixed sliding window approaches for connectivity can detect modulations of connectivity, but mostly if the underlying dynamics operate on moderate to slow timescales. In practice, this can be a drawback, as state durations can vary significantly in empirical data.