Tracking whole-brain connectivity dynamics in the resting state

Estimates of segregation and overlap of functional connectivity networks in the human cerebral cortex

The organization of the human cerebral cortex has recently been explored using techniques for parcellating the cortex into distinct functionally coupled networks. The divergent and convergent nature of cortico-cortical anatomic connections suggests the need to consider the possibility of regions belonging to multiple networks and hierarchies among networks. Here we applied the Latent Dirichlet Allocation (LDA) model and spatial independent component analysis (ICA) to solve for functionally coupled cerebral networks without assuming that cortical regions belong to a single network. Data analyzed included 1000 subjects from the Brain Genomics Superstruct Project (GSP) and 12 high quality individual subjects from the Human Connectome Project (HCP). The organization of the cerebral cortex was similar regardless of whether a winner-take-all approach or the more relaxed constraints of LDA (or ICA) were imposed. This suggests that large-scale networks may function as partially isolated modules. Several notable interactions among networks were uncovered by the LDA analysis. Many association regions belong to at least two networks, while somatomotor and early visual cortices are especially isolated. As examples of interaction, the precuneus, lateral temporal cortex, medial prefrontal cortex and posterior parietal cortex participate in multiple paralimbic networks that together comprise subsystems of the default network. In addition, regions at or near the frontal eye field and human lateral intraparietal area homologue participate in multiple hierarchically organized networks. These observations were replicated in both datasets and could be detected (and replicated) in individual subjects from the HCP.

Detecting functional connectivity change points for single-subject fMRI data

Recently in functional magnetic resonance imaging (fMRI) studies there has been an increased interest in understanding the dynamic manner in which brain regions communicate with one another, as subjects perform a set of experimental tasks or as their psychological state changes. Dynamic Connectivity Regression (DCR) is a data-driven technique used for detecting temporal change points in functional connectivity between brain regions where the number and location of the change points are unknown a priori. After finding the change points, DCR estimates a graph or set of relationships between the brain regions for data that falls between pairs of change points. In previous work, the method was predominantly validated using multi-subject data. In this paper, we concentrate on single-subject data and introduce a new DCR algorithm. The new algorithm increases accuracy for individual subject data with a small number of observations and reduces the number of false positives in the estimated undirected graphs. We also introduce a new Likelihood Ratio test for comparing sparse graphs across (or within) subjects; thus allowing us to determine whether data should be combined across subjects. We perform an extensive simulation analysis on vector autoregression (VAR) data as well as to an fMRI data set from a study (n = 23) of a state anxiety induction using a socially evaluative threat challenge. The focus on single-subject data allows us to study the variation between individuals and may provide us with a deeper knowledge of the workings of the brain.

Dynamic functional connectivity: promise, issues, and interpretations

The brain must dynamically integrate, coordinate, and respond to internal and external stimuli across multiple time scales. Non-invasive measurements of brain activity with fMRI have greatly advanced our understanding of the large-scale functional organization supporting these fundamental features of brain function. Conclusions from previous resting-state fMRI investigations were based upon static descriptions of functional connectivity (FC), and only recently studies have begun to capitalize on the wealth of information contained within the temporal features of spontaneous BOLD FC. Emerging evidence suggests that dynamic FC metrics may index changes in macroscopic neural activity patterns underlying critical aspects of cognition and behavior, though limitations with regard to analysis and interpretation remain. Here, we review recent findings, methodological considerations, neural and behavioral correlates, and future directions in the emerging field of dynamic FC investigations.

Task-based core-periphery organization of human brain dynamics

As a person learns a new skill, distinct synapses, brain regions, and circuits are engaged and change over time. In this paper, we develop methods to examine patterns of correlated activity across a large set of brain regions. Our goal is to identify properties that enable robust learning of a motor skill. We measure brain activity during motor sequencing and characterize network properties based on coherent activity between brain regions. Using recently developed algorithms to detect time-evolving communities, we find that the complex reconfiguration patterns of the brain's putative functional modules that control learning can be described parsimoniously by the combined presence of a relatively stiff temporal core that is composed primarily of sensorimotor and visual regions whose connectivity changes little in time and a flexible temporal periphery that is composed primarily of multimodal association regions whose connectivity changes frequently. The separation between temporal core and periphery changes over the course of training and, importantly, is a good predictor of individual differences in learning success. The core of dynamically stiff regions exhibits dense connectivity, which is consistent with notions of core-periphery organization established previously in social networks. Our results demonstrate that core-periphery organization provides an insightful way to understand how putative functional modules are linked. This, in turn, enables the prediction of fundamental human capacities, including the production of complex goal-directed behavior.

When someone learns a new skill, his/her brain dynamically alters individual synapses, regional activity, and larger-scale circuits. In this paper, we capture some of these dynamics by measuring and characterizing patterns of coherent brain activity during the learning of a motor skill. We extract time-evolving communities from these patterns and find that a temporal core that is composed primarily of primary sensorimotor and visual regions reconfigures little over time, whereas a periphery that is composed primarily of multimodal association regions reconfigures frequently. The core consists of densely connected nodes, and the periphery consists of sparsely connected nodes. Individual participants with a larger separation between core and periphery learn better in subsequent training sessions than individuals with a smaller separation. Conceptually, core-periphery organization provides a framework in which to understand how putative functional modules are linked. This, in turn, enables the prediction of fundamental human capacities, including the production of complex goal-directed behavior.

Quasi-periodic patterns (QPP): large-scale dynamics in resting state fMRI that correlate with local infraslow electrical activity

Functional connectivity measurements from resting state blood-oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) are proving a powerful tool to probe both normal brain function and neuropsychiatric disorders. However, the neural mechanisms that coordinate these large networks are poorly understood, particularly in the context of the growing interest in network dynamics. Recent work in anesthetized rats has shown that the spontaneous BOLD fluctuations are tightly linked to infraslow local field potentials (LFPs) that are seldom recorded but comparable in frequency to the slow BOLD fluctuations. These findings support the hypothesis that long-range coordination involves low frequency neural oscillations and establishes infraslow LFPs as an excellent candidate for probing the neural underpinnings of the BOLD spatiotemporal patterns observed in both rats and humans. To further examine the link between large-scale network dynamics and infraslow LFPs, simultaneous fMRI and microelectrode recording were performed in anesthetized rats. Using an optimized filter to isolate shared components of the signals, we found that time-lagged correlation between infraslow LFPs and BOLD is comparable in spatial extent and timing to a quasi-periodic pattern (QPP) found from BOLD alone, suggesting that fMRI-measured QPPs and the infraslow LFPs share a common mechanism. As fMRI allows spatial resolution and whole brain coverage not available with electroencephalography, QPPs can be used to better understand the role of infraslow oscillations in normal brain function and neurological or psychiatric disorders.

Dynamic cooperation and competition between brain systems during cognitive control

The human brain is characterized by a remarkable ability to adapt its information processing based on current goals. This ability, which is encompassed by the psychological construct of cognitive control, involves activity throughout large-scale, specialized brain systems that support segregated functions at rest and during active task performance. Based on recent research, we propose an account in which control functions rely on transitory changes in patterns of cooperation and competition between neural systems. This account challenges current conceptualizations of control as relying on segregated or antagonistic activity of specialized brain systems. Accordingly, we argue that the study of transitory task-based interactions between brain systems is critical to understanding the flexibility of normal cognitive control and its disruption in pathological conditions.

Structure and function of complex brain networks

An increasing number of theoretical and empirical studies approach the function of the human brain from a network perspective. The analysis of brain networks is made feasible by the development of new imaging acquisition methods as well as new tools from graph theory and dynamical systems. This review surveys some of these methodological advances and summarizes recent findings on the architecture of structural and functional brain networks. Studies of the structural connectome reveal several modules or network communities that are interlinked by hub regions mediating communication processes between modules. Recent network analyses have shown that network hubs form a densely linked collective called a "rich club," centrally positioned for attracting and dispersing signal traffic. In parallel, recordings of resting and task-evoked neural activity have revealed distinct resting-state networks that contribute to functions in distinct cognitive domains. Network methods are increasingly applied in a clinical context, and their promise for elucidating neural substrates of brain and mental disorders is discussed.

Modulatory interactions of resting-state brain functional connectivity

The functional brain connectivity studies are generally based on the synchronization of the resting-state functional magnetic resonance imaging (fMRI) signals. Functional connectivity measures usually assume a stable relationship over time; however, accumulating studies have reported time-varying properties of strength and spatial distribution of functional connectivity. The present study explored the modulation of functional connectivity between two regions by a third region using the physiophysiological interaction (PPI) technique. We first identified eight brain networks and two regions of interest (ROIs) representing each of the networks using a spatial independent component analysis. A voxel-wise analysis was conducted to identify regions that showed modulatory interactions (PPI) with the two ROIs of each network. Mostly, positive modulatory interactions were observed within regions involved in the same system. For example, the two regions of the dorsal attention network revealed modulatory interactions with the regions related to attention, while the two regions of the extrastriate network revealed modulatory interactions with the regions in the visual cortex. In contrast, the two regions of the default mode network (DMN) revealed negative modulatory interactions with the regions in the executive network, and vice versa, suggesting that the activities of one network may be associated with smaller within network connectivity of the competing network. These results validate the use of PPI analysis to study modulation of resting-state functional connectivity by a third region. The modulatory effects may provide a better understanding of complex brain functions.

Functional connectivity in the developing brain: a longitudinal study from 4 to 9months of age

We characterize the development of intrinsic connectivity networks (ICNs) from 4 to 9months of age with resting state magnetic resonance imaging performed on sleeping infants without sedative medication. Data is analyzed with independent component analysis (ICA). Using both low (30 components) and high (100 components) ICA model order decompositions, we find that the functional network connectivity (FNC) map is largely similar at both 4 and 9months. However at 9months the connectivity strength decreases within local networks and increases between more distant networks. The connectivity within the default-mode network, which contains both local and more distant nodes, also increases in strength with age. The low frequency power spectrum increases with age only in the posterior cingulate cortex and posterior default mode network. These findings are consistent with a general developmental pattern of increasing longer distance functional connectivity over the first year of life and raise questions regarding the developmental importance of the posterior cingulate at this age.

GroupICA dual regression analysis of resting state networks in a behavioral variant of frontotemporal dementia

Functional MRI studies have revealed changes in default-mode and salience networks in neurodegenerative dementias, especially in Alzheimer's disease (AD). The purpose of this study was to analyze the whole brain cortex resting state networks (RSNs) in patients with behavioral variant frontotemporal dementia (bvFTD) by using resting state functional MRI (rfMRI). The group specific RSNs were identified by high model order independent component analysis (ICA) and a dual regression technique was used to detect between-group differences in the RSNs with p < 0.05 threshold corrected for multiple comparisons. A y-concatenation method was used to correct for multiple comparisons for multiple independent components, gray matter differences as well as the voxel level. We found increased connectivity in several networks within patients with bvFTD compared to the control group. The most prominent enhancement was seen in the right frontotemporal area and insula. A significant increase in functional connectivity was also detected in the left dorsal attention network (DAN), in anterior paracingulate—a default mode sub-network as well as in the anterior parts of the frontal pole. Notably the increased patterns of connectivity were seen in areas around atrophic regions. The present results demonstrate abnormal increased connectivity in several important brain networks including the DAN and default-mode network (DMN) in patients with bvFTD. These changes may be associated with decline in executive functions and attention as well as apathy, which are the major cognitive and neuropsychiatric defects in patients with frontotemporal dementia.

High-speed real-time resting-state FMRI using multi-slab echo-volumar imaging

We recently demonstrated that ultra-high-speed real-time fMRI using multi-slab echo-volumar imaging (MEVI) significantly increases sensitivity for mapping task-related activation and resting-state networks (RSNs) compared to echo-planar imaging (Posse et al., 2012). In the present study we characterize the sensitivity of MEVI for mapping RSN connectivity dynamics, comparing independent component analysis (ICA) and a novel seed-based connectivity analysis (SBCA) that combines sliding-window correlation analysis with meta-statistics. This SBCA approach is shown to minimize the effects of confounds, such as movement, and CSF and white matter signal changes, and enables real-time monitoring of RSN dynamics at time scales of tens of seconds. We demonstrate highly sensitive mapping of eloquent cortex in the vicinity of brain tumors and arterio-venous malformations, and detection of abnormal resting-state connectivity in epilepsy. In patients with motor impairment, resting-state fMRI provided focal localization of sensorimotor cortex compared with more diffuse activation in task-based fMRI. The fast acquisition speed of MEVI enabled segregation of cardiac-related signal pulsation using ICA, which revealed distinct regional differences in pulsation amplitude and waveform, elevated signal pulsation in patients with arterio-venous malformations and a trend toward reduced pulsatility in gray matter of patients compared with healthy controls. Mapping cardiac pulsation in cortical gray matter may carry important functional information that distinguishes healthy from diseased tissue vasculature. This novel fMRI methodology is particularly promising for mapping eloquent cortex in patients with neurological disease, having variable degree of cooperation in task-based fMRI. In conclusion, ultra-high-real-time speed fMRI enhances the sensitivity of mapping the dynamics of resting-state connectivity and cerebro-vascular pulsatility for clinical and neuroscience research applications.

Dynamical intrinsic functional architecture of the brain during absence seizures

Epilepsy is characterized by recurrent and temporary brain dysfunction due to discharges of interconnected groups of neurons. The brain of epilepsy patients has a dynamic bifurcation that switches between epileptic and normal states. The dysfunctional state involves large-scale brain networks. It is very important to understand the network mechanisms of seizure initiation, maintenance, and termination in epilepsy. Absence epilepsy provides a unique model for neuroimaging investigation on dynamic evolutions of brain networks over seizure repertoire. By using a dynamic functional connectivity and graph theoretical analyses to study absence seizures (AS), we aimed to obtain transition of network properties that account for seizure onset and offset. We measured resting-state functional magnetic resonance imaging and simultaneous electroencephalography (EEG) from children with AS. We used simultaneous EEG to define the preictal, ictal and postictal intervals of seizures. We measured dynamic connectivity maps of the thalamus network and the default mode network (DMN), as well as functional connectome topologies, during the three different seizure intervals. The analysis of dynamic changes of anti-correlation between the thalamus and the DMN is consistent with an inhibitory effect of seizures on the default mode of brain function, which gradually fades out after seizure onset. Also, we observed complex transitions of functional network topology, implicating adaptive reconfiguration of functional brain networks. In conclusion, our work revealed novel insights into modifications in large-scale functional connectome during AS, which may contribute to a better understanding the network mechanisms of state bifurcations in epileptogenesis.

A comprehensive assessment of regional variation in the impact of head micromovements on functional connectomics

Functional connectomics is one of the most rapidly expanding areas of neuroimaging research. Yet, concerns remain regarding the use of resting-state fMRI (R-fMRI) to characterize inter-individual variation in the functional connectome. In particular, recent findings that "micro" head movements can introduce artifactual inter-individual and group-related differences in R-fMRI metrics have raised concerns. Here, we first build on prior demonstrations of regional variation in the magnitude of framewise displacements associated with a given head movement, by providing a comprehensive voxel-based examination of the impact of motion on the BOLD signal (i.e., motion-BOLD relationships). Positive motion-BOLD relationships were detected in primary and supplementary motor areas, particularly in low motion datasets. Negative motion-BOLD relationships were most prominent in prefrontal regions, and expanded throughout the brain in high motion datasets (e.g., children). Scrubbing of volumes with FD>0.2 effectively removed negative but not positive correlations; these findings suggest that positive relationships may reflect neural origins of motion while negative relationships are likely to originate from motion artifact. We also examined the ability of motion correction strategies to eliminate artifactual differences related to motion among individuals and between groups for a broad array of voxel-wise R-fMRI metrics. Residual relationships between motion and the examined R-fMRI metrics remained for all correction approaches, underscoring the need to covary motion effects at the group-level. Notably, global signal regression reduced relationships between motion and inter-individual differences in correlation-based R-fMRI metrics; Z-standardization (mean-centering and variance normalization) of subject-level maps for R-fMRI metrics prior to group-level analyses demonstrated similar advantages. Finally, our test-retest (TRT) analyses revealed significant motion effects on TRT reliability for R-fMRI metrics. Generally, motion compromised reliability of R-fMRI metrics, with the exception of those based on frequency characteristics - particularly, amplitude of low frequency fluctuations (ALFF). The implications of our findings for decision-making regarding the assessment and correction of motion are discussed, as are insights into potential differences among volume-based metrics of motion.

Multi-task connectivity reveals flexible hubs for adaptive task control

Extensive evidence suggests the human ability to adaptively implement a wide variety of tasks is preferentially due to the operation of a fronto-parietal brain network. We hypothesized that this network’s adaptability is made possible by ‘flexible hubs’ – brain regions that rapidly update their pattern of global functional connectivity according to task demands. We utilized recent advances in characterizing brain network organization and dynamics to identify mechanisms consistent with the flexible hub theory. We found that the fronto-parietal network’s brain-wide functional connectivity pattern shifted more than other networks’ across a variety of task states, and that these connectivity patterns could be used to identify the current task. Further, these patterns were consistent across practiced and novel tasks, suggesting reuse of flexible hub connectivity patterns facilitates adaptive (novel) task performance. Together, these findings support a central role for fronto-parietal flexible hubs in cognitive control and adaptive implementation of task demands generally.

Neural correlates of time-varying functional connectivity in the rat

Functional connectivity between brain regions, measured with resting state functional magnetic resonance imaging, holds great potential for understanding the basis of behavior and neuropsychiatric diseases. Recently it has become clear that correlations between the blood oxygenation level dependent (BOLD) signals from different areas vary over the course of a typical scan (6-10 min in length), though the changes are obscured by standard methods of analysis that assume the relationships are stationary. Unfortunately, because similar variability is observed in signals that share no temporal information, it is unclear which dynamic changes are related to underlying neural events. To examine this question, BOLD data were recorded simultaneously with local field potentials (LFP) from interhemispheric primary somatosensory cortex (SI) in anesthetized rats. LFP signals were converted into band-limited power (BLP) signals including delta, theta, alpha, beta and gamma. Correlation between signals from interhemispheric SI was performed in sliding windows to produce signals of correlation over time for BOLD and each BLP band. Both BOLD and BLP signals showed large changes in correlation over time and the changes in BOLD were significantly correlated to the changes in BLP. The strongest relationship was seen when using the theta, beta and gamma bands. Interestingly, while steady-state BOLD and BLP correlate with the global fMRI signal, dynamic BOLD becomes more like dynamic BLP after the global signal is regressed. As BOLD sliding window connectivity is partially reflecting underlying LFP changes, the present study suggests it may be a valuable method of studying dynamic changes in brain states.

BrainNet Viewer: a network visualization tool for human brain connectomics

The human brain is a complex system whose topological organization can be represented using connectomics. Recent studies have shown that human connectomes can be constructed using various neuroimaging technologies and further characterized using sophisticated analytic strategies, such as graph theory. These methods reveal the intriguing topological architectures of human brain networks in healthy populations and explore the changes throughout normal development and aging and under various pathological conditions. However, given the huge complexity of this methodology, toolboxes for graph-based network visualization are still lacking. Here, using MATLAB with a graphical user interface (GUI), we developed a graph-theoretical network visualization toolbox, called BrainNet Viewer, to illustrate human connectomes as ball-and-stick models. Within this toolbox, several combinations of defined files with connectome information can be loaded to display different combinations of brain surface, nodes and edges. In addition, display properties, such as the color and size of network elements or the layout of the figure, can be adjusted within a comprehensive but easy-to-use settings panel. Moreover, BrainNet Viewer draws the brain surface, nodes and edges in sequence and displays brain networks in multiple views, as required by the user. The figure can be manipulated with certain interaction functions to display more detailed information. Furthermore, the figures can be exported as commonly used image file formats or demonstration video for further use. BrainNet Viewer helps researchers to visualize brain networks in an easy, flexible and quick manner, and this software is freely available on the NITRC website (www.nitrc.org/projects/bnv/).

BOLD Frequency Power Indexes Working Memory Performance

Electrophysiology studies routinely investigate the relationship between neural oscillations and task performance. However, the sluggish nature of the BOLD response means that few researchers have investigated the spectral properties of the BOLD signal in a similar manner. For the first time we have applied group ICA to fMRI data collected during a standard working memory task (delayed match-to-sample) and using a multivariate analysis, we investigate the relationship between working memory performance (accuracy and reaction time) and BOLD spectral power within functional networks. Our results indicate that BOLD spectral power within specific networks (visual, temporal-parietal, posterior default-mode network, salience network, basal ganglia) correlated with task accuracy. Multivariate analyses show that the relationship between task accuracy and BOLD spectral power is stronger than the relationship between BOLD spectral power and other variables (age, gender, head movement, and neuropsychological measures). A traditional General Linear Model (GLM) analysis found no significant group differences, or regions that covaried in signal intensity with task accuracy, suggesting that BOLD spectral power holds unique information that is lost in a standard GLM approach. We suggest that the combination of ICA and BOLD spectral power is a useful novel index of cognitive performance that may be more sensitive to brain-behavior relationships than traditional approaches.

Multisubject independent component analysis of fMRI: a decade of intrinsic networks, default mode, and neurodiagnostic discovery

Since the discovery of functional connectivity in fMRI data (i.e., temporal correlations between spatially distinct regions of the brain) there has been a considerable amount of work in this field. One important focus has been on the analysis of brain connectivity using the concept of networks instead of regions. Approximately ten years ago, two important research areas grew out of this concept. First, a network proposed to be "a default mode of brain function" since dubbed the default mode network was proposed by Raichle. Secondly, multisubject or group independent component analysis (ICA) provided a data-driven approach to study properties of brain networks, including the default mode network. In this paper, we provide a focused review of how ICA has contributed to the study of intrinsic networks. We discuss some methodological considerations for group ICA and highlight multiple analytic approaches for studying brain networks. We also show examples of some of the differences observed in the default mode and resting networks in the diseased brain. In summary, we are in exciting times and still just beginning to reap the benefits of the richness of functional brain networks as well as available analytic approaches.

Electroconvulsive therapy response in major depressive disorder: a pilot functional network connectivity resting state FMRI investigation

Major depressive disorder (MDD) is associated with increased functional connectivity in specific neural networks. Electroconvulsive therapy (ECT), the gold-standard treatment for acute, treatment-resistant MDD, but temporal dependencies between networks associated with ECT response have yet to be investigated. In the present longitudinal, case-control investigation, we used independent component analysis to identify distinct networks of brain regions with temporally coherent hemodynamic signal change and functional network connectivity (FNC) to assess component time course correlations across these networks. MDD subjects completed imaging and clinical assessments immediately prior to the ECT series and a minimum of 5 days after the last ECT treatment. We focused our analysis on four networks affected in MDD: the subcallosal cingulate gyrus, default mode, dorsal lateral prefrontal cortex, and dorsal medial prefrontal cortex (DMPFC). In an older sample of ECT subjects (n = 12) with MDD, remission associated with the ECT series reverses the relationship from negative to positive between the posterior default mode (p_DM) and two other networks: the DMPFC and left dorsal lateral prefrontal cortex (l_DLPFC). Relative to demographically healthy subjects (n = 12), the FNC between the p_DM areas and the DMPFC normalizes with ECT response. The FNC changes following treatment did not correlate with symptom improvement; however, a direct comparison between ECT remitters and non-remitters showed the pattern of increased FNC between the p_DM and l_DLPFC following ECT to be specific to those who responded to the treatment. The differences between ECT remitters and non-remitters suggest that this increased FNC between p_DM areas and the left dorsolateral prefrontal cortex is a neural correlate and potential biomarker of recovery from a depressed episode.