EEG Signatures of Dynamic Functional Network Connectivity States

Spatial-temporal-spectral EEG patterns of BOLD functional network connectivity dynamics

OBJECTIVE

Growing interest in the examination of large-scale brain network functional connectivity dynamics is accompanied by an effort to find the electrophysiological correlates. The commonly used constraints applied to spatial and spectral domains during electroencephalogram (EEG) data analysis may leave part of the neural activity unrecognized. We propose an approach that blindly reveals multimodal EEG spectral patterns that are related to the dynamics of the BOLD functional network connectivity.

APPROACH

The blind decomposition of EEG spectrogram by parallel factor analysis has been shown to be a useful technique for uncovering patterns of neural activity. The simultaneously acquired BOLD fMRI data were decomposed by independent component analysis. Dynamic functional connectivity was computed on the component's time series using a sliding window correlation, and between-network connectivity states were then defined based on the values of the correlation coefficients. ANOVA tests were performed to assess the relationships between the dynamics of between-network connectivity states and the fluctuations of EEG spectral patterns.

MAIN RESULTS

We found three patterns related to the dynamics of between-network connectivity states. The first pattern has dominant peaks in the alpha, beta, and gamma bands and is related to the dynamics between the auditory, sensorimotor, and attentional networks. The second pattern, with dominant peaks in the theta and low alpha bands, is related to the visual and default mode network. The third pattern, also with peaks in the theta and low alpha bands, is related to the auditory and frontal network.

SIGNIFICANCE

Our previous findings revealed a relationship between EEG spectral pattern fluctuations and the hemodynamics of large-scale brain networks. In this study, we suggest that the relationship also exists at the level of functional connectivity dynamics among large-scale brain networks when no standard spatial and spectral constraints are applied on the EEG data.

EEG-Informed fMRI: A Review of Data Analysis Methods

The simultaneous acquisition of electroencephalography (EEG) with functional magnetic resonance imaging (fMRI) is a very promising non-invasive technique for the study of human brain function. Despite continuous improvements, it remains a challenging technique, and a standard methodology for data analysis is yet to be established. Here we review the methodologies that are currently available to address the challenges at each step of the data analysis pipeline. We start by surveying methods for pre-processing both EEG and fMRI data. On the EEG side, we focus on the correction for several MR-induced artifacts, particularly the gradient and pulse artifacts, as well as other sources of EEG artifacts. On the fMRI side, we consider image artifacts induced by the presence of EEG hardware inside the MR scanner, and the contamination of the fMRI signal by physiological noise of non-neuronal origin, including a review of several approaches to model and remove it. We then provide an overview of the approaches specifically employed for the integration of EEG and fMRI when using EEG to predict the blood oxygenation level dependent (BOLD) fMRI signal, the so-called EEG-informed fMRI integration strategy, the most commonly used strategy in EEG-fMRI research. Finally, we systematically review methods used for the extraction of EEG features reflecting neuronal phenomena of interest.

Temporal and spectral characteristics of dynamic functional connectivity between resting-state networks reveal information beyond static connectivity

Estimation of functional connectivity (FC) has become an increasingly powerful tool for investigating healthy and abnormal brain function. Static connectivity, in particular, has played a large part in guiding conclusions from the majority of resting-state functional MRI studies. However, accumulating evidence points to the presence of temporal fluctuations in FC, leading to increasing interest in estimating FC as a dynamic quantity. One central issue that has arisen in this new view of connectivity is the dramatic increase in complexity caused by dynamic functional connectivity (dFC) estimation. To computationally handle this increased complexity, a limited set of dFC properties, primarily the mean and variance, have generally been considered. Additionally, it remains unclear how to integrate the increased information from dFC into pattern recognition techniques for subject-level prediction. In this study, we propose an approach to address these two issues based on a large number of previously unexplored temporal and spectral features of dynamic functional connectivity. A Generalized Autoregressive Conditional Heteroskedasticity (GARCH) model is used to estimate time-varying patterns of functional connectivity between resting-state networks. Time-frequency analysis is then performed on dFC estimates, and a large number of previously unexplored temporal and spectral features drawn from signal processing literature are extracted for dFC estimates. We apply the investigated features to two neurologic populations of interest, healthy controls and patients with temporal lobe epilepsy, and show that the proposed approach leads to substantial increases in predictive performance compared to both traditional estimates of static connectivity as well as current approaches to dFC. Variable importance is assessed and shows that there are several quantities that can be extracted from dFC signal which are more informative than the traditional mean or variance of dFC. This work illuminates many previously unexplored facets of the dynamic properties of functional connectivity between resting-state networks, and provides a platform for dynamic functional connectivity analysis that facilitates its usage as an investigative measure for healthy as well as abnormal brain function.

Associations between Functional Connectivity Dynamics and BOLD Dynamics Are Heterogeneous Across Brain Networks

Functional brain imaging has revealed two types of dynamic patterns of brain in the resting-state: the dynamics of spontaneous brain activities and the dynamics of functional interconnections between spontaneous brain activities. Although these two types of brain dynamics are usually investigated separately in the literature, recent functional magnetic resonance imaging (fMRI) studies have shown that they exhibit similar spatial patterns, suggesting the dynamics of spontaneous brain activities and the dynamics of their interconnections are associated with each other. In this study, we characterized the local blood oxygenation level dependent (BOLD) dynamics and the functional connectivity dynamics (FCD) in the resting-state, and then investigated their between-region associations. Time-varying FC was estimated as time-varying correlation coefficients using a sliding-window method, and the temporal variability of BOLD and time-varying FC were used to quantify the BOLD dynamics and the FCD, respectively. Our results showed that the BOLD dynamics and the FCD exhibit similar spatial patterns, and they are significantly associated across brain regions. Importantly, such associations are opposite for different types of FC (e.g., within-network FCD are negatively correlated with the BOLD dynamics but the between-network FCD are positively correlated with the BOLD dynamics), and they are spatially heterogeneous across brain networks. The identified heterogeneous associations between the BOLD dynamics and the FCD appear to convey related or even distinct information and might underscore the potential mechanism of brain coordination and co-evolution.

Schizophrenia Shows Disrupted Links between Brain Volume and Dynamic Functional Connectivity

Studies featuring multimodal neuroimaging data fusion for understanding brain function and structure, or disease characterization, leverage the partial information available in each of the modalities to reveal data variations not exhibited through the independent analyses. Similar to other complex syndromes, the characteristic brain abnormalities in schizophrenia may be better understood with the help of the additional information conveyed by leveraging an advanced modeling method involving multiple modalities. In this study, we propose a novel framework to fuse feature spaces corresponding to functional magnetic resonance imaging (functional) and gray matter (structural) data from 151 schizophrenia patients and 163 healthy controls. In particular, the features for the functional and structural modalities include dynamic (i.e., time-varying) functional network connectivity (dFNC) maps and the intensities of the gray matter (GM) maps, respectively. The dFNC maps are estimated from group independent component analysis (ICA) network time-courses by first computing windowed functional correlations using a sliding window approach, and then estimating subject specific states from this windowed data using temporal ICA followed by spatio-temporal regression. For each subject, the functional data features are horizontally concatenated with the corresponding GM features to form a combined feature space that is subsequently decomposed through a symmetric multimodal fusion approach involving a combination of multiset canonical correlation analysis (mCCA) and joint ICA (jICA). Our novel combined analyses successfully linked changes in the two modalities and revealed significantly disrupted links between GM volumes and time-varying functional connectivity in schizophrenia. Consistent with prior research, we found significant group differences in GM comprising regions in the superior parietal lobule, precuneus, postcentral gyrus, medial/superior frontal gyrus, superior/middle temporal gyrus, insula and fusiform gyrus, and several significant aberrations in the inter-regional functional connectivity strength as well. Importantly, structural and dFNC measures have independently shown changes associated with schizophrenia, and in this work we begin the process of evaluating the links between the two, which could shed light on the illness beyond what we can learn from a single imaging modality. In future work, we plan to evaluate replication of the inferred structure-function relationships in independent partitions of larger multi-modal schizophrenia datasets.

Abnormal intrinsic brain functional network dynamics in Parkinson's disease

See Nieuwhof and Helmich (doi:10.1093/brain/awx267 ) for a scientific commentary on this article .

Parkinson’s disease is a neurodegenerative disorder characterized by nigrostriatal dopamine depletion. Previous studies measuring spontaneous brain activity using resting state functional magnetic resonance imaging have reported abnormal changes in broadly distributed whole-brain networks. Although resting state functional connectivity, estimating temporal correlations between brain regions, is measured with the assumption that intrinsic fluctuations throughout the scan are stable, dynamic changes of functional connectivity have recently been suggested to reflect aspects of functional capacity of neural systems, and thus may serve as biomarkers of disease. The present work is the first study to investigate the dynamic functional connectivity in patients with Parkinson’s disease, with a focus on the temporal properties of functional connectivity states as well as the variability of network topological organization using resting state functional magnetic resonance imaging. Thirty-one Parkinson’s disease patients and 23 healthy controls were studied using group spatial independent component analysis, a sliding windows approach, and graph-theory methods. The dynamic functional connectivity analyses suggested two discrete connectivity configurations: a more frequent, sparsely connected within-network state (State I) and a less frequent, more strongly interconnected between-network state (State II). In patients with Parkinson’s disease, the occurrence of the sparsely connected State I dropped by 12.62%, while the expression of the more strongly interconnected State II increased by the same amount. This was consistent with the altered temporal properties of the dynamic functional connectivity characterized by a shortening of the dwell time of State I and by a proportional increase of the dwell time pattern in State II. These changes are suggestive of a reduction in functional segregation among networks and are correlated with the clinical severity of Parkinson’s disease symptoms. Additionally, there was a higher variability in the network global efficiency, suggesting an abnormal global integration of the brain networks. The altered functional segregation and abnormal global integration in brain networks confirmed the vulnerability of functional connectivity networks in Parkinson’s disease.

Ten Key Observations on the Analysis of Resting-state Functional MR Imaging Data Using Independent Component Analysis

For more than 20 years, the powerful, flexible family of independent component analysis (ICA) techniques has been used to examine spatial, temporal, and subject variation in functional magnetic resonance (fMR) imaging data. This article provides an overview of 10 key principles in the basic and advanced application of ICA to resting-state fMR imaging. ICA's core advantages include robustness to artifact; false-positives and autocorrelation; adaptability to variant study designs; agnosticism to the temporal evolution of fMR imaging signals; and ability to extract, identify, and analyze neural networks. ICA remains in the vanguard of fMRI methods development.

Characterizing dynamic amplitude of low-frequency fluctuation and its relationship with dynamic functional connectivity: An application to schizophrenia

The human brain is a highly dynamic system with non-stationary neural activity and rapidly-changing neural interaction. Resting-state dynamic functional connectivity (dFC) has been widely studied during recent years, and the emerging aberrant dFC patterns have been identified as important features of many mental disorders such as schizophrenia (SZ). However, only focusing on the time-varying patterns in FC is not enough, since the local neural activity itself (in contrast to the inter-connectivity) is also found to be highly fluctuating from research using high-temporal-resolution imaging techniques. Exploring the time-varying patterns in brain activity and their relationships with time-varying brain connectivity is important for advancing our understanding of the co-evolutionary property of brain network and the underlying mechanism of brain dynamics. In this study, we introduced a framework for characterizing time-varying brain activity and exploring its associations with time-varying brain connectivity, and applied this framework to a resting-state fMRI dataset including 151 SZ patients and 163 age- and gender matched healthy controls (HCs). In this framework, 48 brain regions were first identified as intrinsic connectivity networks (ICNs) using group independent component analysis (GICA). A sliding window approach was then adopted for the estimation of dynamic amplitude of low-frequency fluctuation (dALFF) and dFC, which were used to measure time-varying brain activity and time-varying brain connectivity respectively. The dALFF was further clustered into six reoccurring states by the k-means clustering method and the group difference in occurrences of dALFF states was explored. Lastly, correlation coefficients between dALFF and dFC were calculated and the group difference in these dALFF-dFC correlations was explored. Our results suggested that 1) ALFF of brain regions was highly fluctuating during the resting-state and such dynamic patterns are altered in SZ, 2) dALFF and dFC were correlated in time and their correlations are altered in SZ. The overall results support and expand prior work on abnormalities of brain activity, static FC (sFC) and dFC in SZ, and provide new evidence on aberrant time-varying brain activity and its associations with brain connectivity in SZ, which might underscore the disrupted brain cognitive functions in this mental disorder.

Replicability of time-varying connectivity patterns in large resting state fMRI samples

The past few years have seen an emergence of approaches that leverage temporal changes in whole-brain patterns of functional connectivity (the chronnectome). In this chronnectome study, we investigate the replicability of the human brain’s inter-regional coupling dynamics during rest by evaluating two different dynamic functional network connectivity (dFNC) analysis frameworks using 7 500 functional magnetic resonance imaging (fMRI) datasets. To quantify the extent to which the emergent functional connectivity (FC) patterns are reproducible, we characterize the temporal dynamics by deriving several summary measures across multiple large, independent age-matched samples. Reproducibility was demonstrated through the existence of basic connectivity patterns (FC states) amidst an ensemble of inter-regional connections. Furthermore, application of the methods to conservatively configured (statistically stationary, linear and Gaussian) surrogate datasets revealed that some of the studied state summary measures were indeed statistically significant and also suggested that this class of null model did not explain the fMRI data fully. This extensive testing of reproducibility of similarity statistics also suggests that the estimated FC states are robust against variation in data quality, analysis, grouping, and decomposition methods. We conclude that future investigations probing the functional and neurophysiological relevance of time-varying connectivity assume critical importance.

Neural and metabolic basis of dynamic resting state fMRI

Resting state fMRI (rsfMRI) as a technique showed much initial promise for use in psychiatric and neurological diseases where diagnosis and treatment were difficult. To realize this promise, many groups have moved towards examining "dynamic rsfMRI," which relies on the assumption that rsfMRI measurements on short time scales remain relevant to the underlying neural and metabolic activity. Many dynamic rsfMRI studies have demonstrated differences between clinical or behavioral groups beyond what static rsfMRI measured, suggesting a neurometabolic basis. Correlative studies combining dynamic rsfMRI and other physiological measurements have supported this. However, they also indicate multiple mechanisms and, if using correlation alone, it is difficult to separate cause and effect. Hypothesis-driven studies are needed, a few of which have begun to illuminate the underlying neurometabolic mechanisms that shape observed differences in dynamic rsfMRI. While the number of potential noise sources, potential actual neurometabolic sources, and methodological considerations can seem overwhelming, dynamic rsfMRI provides a rich opportunity in systems neuroscience. Even an incrementally better understanding of the neurometabolic basis of dynamic rsfMRI would expand rsfMRI's research and clinical utility, and the studies described herein take the first steps on that path forward.

Ten Key Observations on the Analysis of Resting-state Functional MR Imaging Data Using Independent Component Analysis

For more than 20 years, the powerful, flexible family of independent component analysis (ICA) techniques has been used to examine spatial, temporal, and subject variation in functional magnetic resonance (fMR) imaging data. This article provides an overview of 10 key principles in the basic and advanced application of ICA to resting-state fMR imaging. ICA's core advantages include robustness to artifact; false-positives and autocorrelation; adaptability to variant study designs; agnosticism to the temporal evolution of fMR imaging signals; and ability to extract, identify, and analyze neural networks. ICA remains in the vanguard of fMRI methods development.

Chronnectomic patterns and neural flexibility underlie executive function

Despite extensive research into executive function (EF), the precise relationship between brain dynamics and flexible cognition remains unknown. Using a large, publicly available dataset (189 participants), we find that functional connections measured throughout 56min of resting state fMRI data comprise five distinct connectivity states. Elevated EF performance as measured outside of the scanner was associated with greater episodes of more frequently occurring connectivity states, and fewer episodes of less frequently occurring connectivity states. Frequently occurring states displayed metastable properties, where cognitive flexibility may be facilitated by attenuated correlations and greater functional connection variability. Less frequently occurring states displayed properties consistent with low arousal and low vigilance. These findings suggest that elevated EF performance may be associated with the propensity to occupy more frequently occurring brain configurations that enable cognitive flexibility, while avoiding less frequently occurring brain configurations related to low arousal/vigilance states. The current findings offer a novel framework for identifying neural processes related to individual differences in executive function.