Elena A. Allen

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.

Tracking Whole-Brain Connectivity Dynamics in the Resting State

Spontaneous fluctuations are a hallmark of recordings of neural signals, emergent over time scales spanning milliseconds and tens of minutes. However, investigations of intrinsic brain organization based on resting-state functional magnetic resonance imaging have largely not taken into account the presence and potential of temporal variability, as most current approaches to examine functional connectivity (FC) implicitly assume that relationships are constant throughout the length of the recording. In this work, we describe an approach to assess whole-brain FC dynamics based on spatial independent component analysis, sliding time window correlation, and k-means clustering of windowed correlation matrices. The method is applied to resting-state data from a large sample (n = 405) of young adults. Our analysis of FC variability highlights particularly flexible connections between regions in lateral parietal and cingulate cortex, and argues against a labeling scheme where such regions are treated as separate and antagonistic entities. Additionally, clustering analysis reveals unanticipated FC states that in part diverge strongly from stationary connectivity patterns and challenge current descriptions of interactions between large-scale networks. Temporal trends in the occurrence of different FC states motivate theories regarding their functional roles and relationships with vigilance/arousal. Overall, we suggest that the study of time-varying aspects of FC can unveil flexibility in the functional coordination between different neural systems, and that the exploitation of these dynamics in further investigations may improve our understanding of behavioral shifts and adaptive processes.

A Baseline for the Multivariate Comparison of Resting-State Networks

As the size of functional and structural MRI datasets expands, it becomes increasingly important to establish a baseline from which diagnostic relevance may be determined, a processing strategy that efficiently prepares data for analysis, and a statistical approach that identifies important effects in a manner that is both robust and reproducible. In this paper, we introduce a multivariate analytic approach that optimizes sensitivity and reduces unnecessary testing. We demonstrate the utility of this mega-analytic approach by identifying the effects of age and gender on the resting-state networks (RSNs) of 603 healthy adolescents and adults (mean age: 23.4 years, range: 12–71 years). Data were collected on the same scanner, preprocessed using an automated analysis pipeline based in SPM, and studied using group independent component analysis. RSNs were identified and evaluated in terms of three primary outcome measures: time course spectral power, spatial map intensity, and functional network connectivity. Results revealed robust effects of age on all three outcome measures, largely indicating decreases in network coherence and connectivity with increasing age. Gender effects were of smaller magnitude but suggested stronger intra-network connectivity in females and more inter-network connectivity in males, particularly with regard to sensorimotor networks. These findings, along with the analysis approach and statistical framework described here, provide a useful baseline for future investigations of brain networks in health and disease.

Comparison of multi-subject ICA methods for analysis of fMRI data.

Spatial independent component analysis (ICA) applied to functional magnetic resonance imaging (fMRI) data identifies functionally connected networks by estimating spatially independent patterns from their linearly mixed fMRI signals. Several multi‐subject ICA approaches estimating subject‐specific time courses (TCs) and spatial maps (SMs) have been developed, however, there has not yet been a full comparison of the implications of their use. Here, we provide extensive comparisons of four multi‐subject ICA approaches in combination with data reduction methods for simulated and fMRI task data. For multi‐subject ICA, the data first undergo reduction at the subject and group levels using principal component analysis (PCA). Comparisons of subject‐specific, spatial concatenation, and group data mean subject‐level reduction strategies using PCA and probabilistic PCA (PPCA) show that computationally intensive PPCA is equivalent to PCA, and that subject‐specific and group data mean subject‐level PCA are preferred because of well‐estimated TCs and SMs. Second, aggregate independent components are estimated using either noise‐free ICA or probabilistic ICA (PICA). Third, subject‐specific SMs and TCs are estimated using back‐reconstruction. We compare several direct group ICA (GICA) back‐reconstruction approaches (GICA1‐GICA3) and an indirect back‐reconstruction approach, spatio‐temporal regression (STR, or dual regression). Results show the earlier group ICA (GICA1) approximates STR, however STR has contradictory assumptions and may show mixed‐component artifacts in estimated SMs. Our evidence‐based recommendation is to use GICA3, introduced here, with subject‐specific PCA and noise‐free ICA, providing the most robust and accurate estimated SMs and TCs in addition to offering an intuitive interpretation. Hum Brain Mapp, 2011.

Dynamic functional connectivity analysis reveals transient states of dysconnectivity in schizophrenia.

Schizophrenia is a psychotic disorder characterized by functional dysconnectivity or abnormal integration between distant brain regions. Recent functional imaging studies have implicated large-scale thalamo-cortical connectivity as being disrupted in patients. However, observed connectivity differences in schizophrenia have been inconsistent between studies, with reports of hyperconnectivity and hypoconnectivity between the same brain regions. Using resting state eyes-closed functional imaging and independent component analysis on a multi-site data that included 151 schizophrenia patients and 163 age- and gender matched healthy controls, we decomposed the functional brain data into 100 components and identified 47 as functionally relevant intrinsic connectivity networks. We subsequently evaluated group differences in functional network connectivity, both in a static sense, computed as the pairwise Pearson correlations between the full network time courses (5.4 minutes in length), and a dynamic sense, computed using sliding windows (44 s in length) and k-means clustering to characterize five discrete functional connectivity states. Static connectivity analysis revealed that compared to healthy controls, patients show significantly stronger connectivity, i.e., hyperconnectivity, between the thalamus and sensory networks (auditory, motor and visual), as well as reduced connectivity (hypoconnectivity) between sensory networks from all modalities. Dynamic analysis suggests that (1), on average, schizophrenia patients spend much less time than healthy controls in states typified by strong, large-scale connectivity, and (2), that abnormal connectivity patterns are more pronounced during these connectivity states. In particular, states exhibiting cortical–subcortical antagonism (anti-correlations) and strong positive connectivity between sensory networks are those that show the group differences of thalamic hyperconnectivity and sensory hypoconnectivity. Group differences are weak or absent during other connectivity states. Dynamic analysis also revealed hypoconnectivity between the putamen and sensory networks during the same states of thalamic hyperconnectivity; notably, this finding cannot be observed in the static connectivity analysis. Finally, in post-hoc analyses we observed that the relationships between sub-cortical low frequency power and connectivity with sensory networks is altered in patients, suggesting different functional interactions between sub-cortical nuclei and sensorimotor cortex during specific connectivity states. While important differences between patients with schizophrenia and healthy controls have been identified, one should interpret the results with caution given the history of medication in patients. Taken together, our results support and expand current knowledge regarding dysconnectivity in schizophrenia, and strongly advocate the use of dynamic analyses to better account for and understand functional connectivity differences.

Capturing inter-subject variability with group independent component analysis of fMRI data: a simulation study.

A key challenge in functional neuroimaging is the meaningful combination of results across subjects. Even in a sample of healthy participants, brain morphology and functional organization exhibit considerable variability, such that no two individuals have the same neural activation at the same location in response to the same stimulus. This inter-subject variability limits inferences at the group-level as average activation patterns may fail to represent the patterns seen in individuals. A promising approach to multi-subject analysis is group independent component analysis (GICA), which identifies group components and reconstructs activations at the individual level. GICA has gained considerable popularity, particularly in studies where temporal response models cannot be specified. However, a comprehensive understanding of the performance of GICA under realistic conditions of inter-subject variability is lacking. In this study we use simulated functional magnetic resonance imaging (fMRI) data to determine the capabilities and limitations of GICA under conditions of spatial, temporal, and amplitude variability. Simulations, generated with the SimTB toolbox, address questions that commonly arise in GICA studies, such as: (1) How well can individual subject activations be estimated and when will spatial variability preclude estimation? (2) Why does component splitting occur and how is it affected by model order? (3) How should we analyze component features to maximize sensitivity to intersubject differences? Overall, our results indicate an excellent capability of GICA to capture between-subject differences and we make a number of recommendations regarding analytic choices for application to functional imaging data.

Deep learning for neuroimaging: a validation study

Deep learning methods have recently made notable advances in the tasks of classification and representation learning. These tasks are important for brain imaging and neuroscience discovery, making the methods attractive for porting to a neuroimagers toolbox. Success of these methods is, in part, explained by the flexibility of deep learning models. However, this flexibility makes the process of porting to new areas a difficult parameter optimization problem. In this work we demonstrate our results (and feasible parameter ranges) in application of deep learning methods to structural and functional brain imaging data. These methods include deep belief networks and their building block the restricted Boltzmann machine. We also describe a novel constraint-based approach to visualizing high dimensional data. We use it to analyze the effect of parameter choices on data transformations. Our results show that deep learning methods are able to learn physiologically important representations and detect latent relations in neuroimaging data.

Exploring the psychosis functional connectome: aberrant intrinsic networks in schizophrenia and bipolar disorder

Intrinsic functional brain networks (INs) are regions showing temporal coherence with one another. These INs are present in the context of a task (as opposed to an undirected task such as rest), albeit modulated to a degree both spatially and temporally. Prominent networks include the default mode, attentional fronto-parietal, executive control, bilateral temporal lobe, and motor networks. The characterization of INs has recently gained considerable momentum, however; most previous studies evaluate only a small subset of the INs (e.g., default mode). In this paper we use independent component analysis to study INs decomposed from functional magnetic resonance imaging data collected in a large group of schizophrenia patients, healthy controls, and individuals with bipolar disorder, while performing an auditory oddball task. Schizophrenia and bipolar disorder share significant overlap in clinical symptoms, brain characteristics, and risk genes which motivates our goal of identifying whether functional imaging data can differentiate the two disorders. We tested for group differences in properties of all identified INs including spatial maps, spectra, and functional network connectivity. A small set of default mode, temporal lobe, and frontal networks with default mode regions appearing to play a key role in all comparisons. Bipolar subjects showed more prominent changes in ventromedial and prefrontal default mode regions whereas schizophrenia patients showed changes in posterior default mode regions. Anti-correlations between left parietal areas and dorsolateral prefrontal cortical areas were different in bipolar and schizophrenia patients and amplitude was significantly different from healthy controls in both patient groups. Patients exhibited similar frequency behavior across multiple networks with decreased low frequency power. In summary, a comprehensive analysis of INs reveals a key role for the default mode in both schizophrenia and bipolar disorder.

SimTB, a simulation toolbox for fMRI data under a model of spatiotemporal separability

We introduce SimTB, a MATLAB toolbox designed to simulate functional magnetic resonance imaging (fMRI) datasets under a model of spatiotemporal separability. The toolbox meets the increasing need of the fMRI community to more comprehensively understand the effects of complex processing strategies by providing a ground truth that estimation methods may be compared against. SimTB captures the fundamental structure of real data, but data generation is fully parameterized and fully controlled by the user, allowing for accurate and precise comparisons. The toolbox offers a wealth of options regarding the number and configuration of spatial sources, implementation of experimental paradigms, inclusion of tissue-specific properties, addition of noise and head movement, and much more. A straightforward data generation method and short computation time (3-10 seconds for each dataset) allow a practitioner to simulate and analyze many datasets to potentially understand a problem from many angles. Beginning MATLAB users can use the SimTB graphical user interface (GUI) to design and execute simulations while experienced users can write batch scripts to automate and customize this process. The toolbox is freely available at http://mialab.mrn.org/software together with sample scripts and tutorials.

Correspondence between structure and function in the human brain at rest.

To further understanding of basic and complex cognitive functions, previous connectome research has identified functional and structural connections of the human brain. Functional connectivity is often measured by using resting-state functional magnetic resonance imaging (rs-fMRI) and is generally interpreted as an indirect measure of neuronal activity. Gray matter (GM) primarily consists of neuronal and glia cell bodies; therefore, it is surprising that the majority of connectome research has excluded GM measures. Therefore, we propose that by exploring where GM corresponds to function would aid in the understanding of both structural and functional connectivity and in turn the human connectome. A cohort of 603 healthy participants underwent structural and functional scanning on the same 3 T scanner at the Mind Research Network. To investigate the spatial correspondence between structure and function, spatial independent component analysis (ICA) was applied separately to both GM density (GMD) maps and to rs-fMRI data. ICA of GM delineates structural components based on the covariation of GMD regions among subjects. For the rs-fMRI data, ICA identified spatial patterns with common temporal features. These decomposed structural and functional components were then compared by spatial correlation. Basal ganglia components exhibited the highest structural to resting-state functional spatial correlation (r = 0.59). Cortical components generally show correspondence between a single structural component and several resting-state functional components. We also studied relationships between the weights of different structural components and identified the precuneus as a hub in GMD structural network correlations. In addition, we analyzed relationships between component weights, age, and gender; concluding that age has a significant effect on structural components.