Detecting overlapped functional clusters in resting state fMRI with Connected Iterative Scan: A graph theory based clustering algorithm

Large-scale functional network overlap is a general property of brain functional organization: Reconciling inconsistent fMRI findings from general-linear-model-based analyses

Functional magnetic resonance imaging (fMRI) studies regularly use univariate general-linear-model-based analyses (GLM). Their findings are often inconsistent across different studies, perhaps because of several fundamental brain properties including functional heterogeneity, balanced excitation and inhibition (E/I), and sparseness of neuronal activities. These properties stipulate heterogeneous neuronal activities in the same voxels and likely limit the sensitivity and specificity of GLM. This paper selectively reviews findings of histological and electrophysiological studies and fMRI spatial independent component analysis (sICA) and reports new findings by applying sICA to two existing datasets. The extant and new findings consistently demonstrate several novel features of brain functional organization not revealed by GLM. They include overlap of large-scale functional networks (FNs) and their concurrent opposite modulations, and no significant modulations in activity of most FNs across the whole brain during any task conditions. These novel features of brain functional organization are highly consistent with the brain's properties of functional heterogeneity, balanced E/I, and sparseness of neuronal activity, and may help reconcile inconsistent GLM findings.

A novel pattern mining approach for identifying cognitive activity in EEG based functional brain networks

The complex nature of neuronal interactions of the human brain has posed many challenges to the research community. To explore the underlying mechanisms of neuronal activity of cohesive brain regions during different cognitive activities, many innovative mathematical and computational models are required. This paper presents a novel Common Functional Pattern Mining approach to demonstrate the similar patterns of interactions due to common behavior of certain brain regions. The electrode sites of EEG-based functional brain network are modeled as a set of transactions and node-based complex network measures as itemsets. These itemsets are transformed into a graph data structure called Functional Pattern Graph. By mining this Functional Pattern Graph, the common functional patterns due to specific brain functioning can be identified. The empirical analyses show the efficiency of the proposed approach in identifying the extent to which the electrode sites (transactions) are similar during various cognitive load states.

An Evolutionary Computation Approach to Examine Functional Brain Plasticity

One common research goal in systems neurosciences is to understand how the functional relationship between a pair of regions of interest (ROIs) evolves over time. Examining neural connectivity in this way is well-suited for the study of developmental processes, learning, and even in recovery or treatment designs in response to injury. For most fMRI based studies, the strength of the functional relationship between two ROIs is defined as the correlation between the average signal representing each region. The drawback to this approach is that much information is lost due to averaging heterogeneous voxels, and therefore, the functional relationship between a ROI-pair that evolve at a spatial scale much finer than the ROIs remain undetected. To address this shortcoming, we introduce a novel evolutionary computation (EC) based voxel-level procedure to examine functional plasticity between an investigator defined ROI-pair by simultaneously using subject-specific BOLD-fMRI data collected from two sessions seperated by finite duration of time. This data-driven procedure detects a sub-region composed of spatially connected voxels from each ROI (a so-called sub-regional-pair) such that the pair shows a significant gain/loss of functional relationship strength across the two time points. The procedure is recursive and iteratively finds all statistically significant sub-regional-pairs within the ROIs. Using this approach, we examine functional plasticity between the default mode network (DMN) and the executive control network (ECN) during recovery from traumatic brain injury (TBI); the study includes 14 TBI and 12 healthy control subjects. We demonstrate that the EC based procedure is able to detect functional plasticity where a traditional averaging based approach fails. The subject-specific plasticity estimates obtained using the EC-procedure are highly consistent across multiple runs. Group-level analyses using these plasticity estimates showed an increase in the strength of functional relationship between DMN and ECN for TBI subjects, which is consistent with prior findings in the TBI-literature. The EC-approach also allowed us to separate sub-regional-pairs contributing to positive and negative plasticity; the detected sub-regional-pairs significantly overlap across runs thus highlighting the reliability of the EC-approach. These sub-regional-pairs may be useful in performing nuanced analyses of brain-behavior relationships during recovery from TBI.

Stable Overlapping Replicator Dynamics for Brain Community Detection

A fundamental means for understanding the brain's organizational structure is to group its spatially disparate regions into functional subnetworks based on their interactions. Most community detection techniques are designed for generating partitions, but certain brain regions are known to interact with multiple subnetworks. Thus, the brain's underlying subnetworks necessarily overlap. In this paper, we propose a technique for identifying overlapping subnetworks from weighted graphs with statistical control over false node inclusion. Our technique improves upon the replicator dynamics formulation by incorporating a graph augmentation strategy to enable subnetwork overlaps, and a graph incrementation scheme for merging subnetworks that might be falsely split by replicator dynamics due to its stringent mutual similarity criterion in defining subnetworks. To statistically control for inclusion of false nodes into the detected subnetworks, we further present a procedure for integrating stability selection into our subnetwork identification technique. We refer to the resulting technique as stable overlapping replicator dynamics (SORD). Our experiments on synthetic data show significantly higher accuracy in subnetwork identification with SORD than several state-of-the-art techniques. We also demonstrate higher test-retest reliability in multiple network measures on the Human Connectome Project data. Further, we illustrate that SORD enables identification of neuroanatomically-meaningful subnetworks and network hubs.

Large-Scale Functional Networks Identified from Resting-State EEG Using Spatial ICA

Several methods have been applied to EEG or MEG signals to detect functional networks. In recent works using MEG/EEG and fMRI data, temporal ICA analysis has been used to extract spatial maps of resting-state networks with or without an atlas-based parcellation of the cortex. Since the links between the fMRI signal and the electromagnetic signals are not fully established, and to avoid any bias, we examined whether EEG alone was able to derive the spatial distribution and temporal characteristics of functional networks. To do so, we propose a two-step original method: 1) An individual multi-frequency data analysis including EEG-based source localisation and spatial independent component analysis, which allowed us to characterize the resting-state networks. 2) A group-level analysis involving a hierarchical clustering procedure to identify reproducible large-scale networks across the population. Compared with large-scale resting-state networks obtained with fMRI, the proposed EEG-based analysis revealed smaller independent networks thanks to the high temporal resolution of EEG, hence hierarchical organization of networks. The comparison showed a substantial overlap between EEG and fMRI networks in motor, premotor, sensory, frontal, and parietal areas. However, there were mismatches between EEG-based and fMRI-based networks in temporal areas, presumably resulting from a poor sensitivity of fMRI in these regions or artefacts in the EEG signals. The proposed method opens the way for studying the high temporal dynamics of networks at the source level thanks to the high temporal resolution of EEG. It would then become possible to study detailed measures of the dynamics of connectivity.

Implications of cortical balanced excitation and inhibition, functional heterogeneity, and sparseness of neuronal activity in fMRI

Blood-oxygenation-level-dependent (BOLD) functional magnetic resonance imaging (fMRI) studies often report inconsistent findings, probably due to brain properties such as balanced excitation and inhibition and functional heterogeneity. These properties indicate that different neurons in the same voxels may show variable activities including concurrent activation and deactivation, that the relationships between BOLD signal and neural activity (i.e., neurovascular coupling) are complex, and that increased BOLD signal may reflect reduced deactivation, increased activation, or both. The traditional general-linear-model-based-analysis (GLM-BA) is a univariate approach, cannot separate different components of BOLD signal mixtures from the same voxels, and may contribute to inconsistent findings of fMRI. Spatial independent component analysis (sICA) is a multivariate approach, can separate the BOLD signal mixture from each voxel into different source signals and measure each separately, and thus may reconcile previous conflicting findings generated by GLM-BA. We propose that methods capable of separating mixed signals such as sICA should be regularly used for more accurately and completely extracting information embedded in fMRI datasets.

Automated iterative reclustering framework for determining hierarchical functional networks in resting state fMRI

To spatially cluster resting state-functional magnetic resonance imaging (rs-fMRI) data into potential networks, there are only a few general approaches that determine the number of networks/clusters, despite a wide variety of techniques proposed for clustering. For individual subjects, extraction of a large number of spatially disjoint clusters results in multiple small networks that are spatio-temporally homogeneous but irreproducible across subjects. Alternatively, extraction of a small number of clusters creates spatially large networks that are temporally heterogeneous but spatially reproducible across subjects. We propose a fully automatic, iterative reclustering framework in which a small number of spatially large, heterogeneous networks are initially extracted to maximize spatial reproducibility. Subsequently, the large networks are iteratively subdivided to create spatially reproducible subnetworks until the overall within-network homogeneity does not increase substantially. The proposed approach discovers a rich network hierarchy in the brain while simultaneously optimizing spatial reproducibility of networks across subjects and individual network homogeneity. We also propose a novel metric to measure the connectivity of brain regions, and in a simulation study show that our connectivity metric and framework perform well in the face of low signal to noise and initial segmentation errors. Experimental results generated using real fMRI data show that the proposed metric improves stability of network clusters across subjects, and generates a meaningful pattern for spatially hierarchical structure of the brain.

Opposite modulation of brain functional networks implicated at low vs. high demand of attention and working memory

BACKGROUND

Functional magnetic resonance imaging (fMRI) studies indicate that the brain organizes its activity into multiple functional networks (FNs) during either resting condition or task-performance. However, the functions of these FNs are not fully understood yet.

METHODOLOGY/PRINCIPAL FINDINGS

To investigate the operation of these FNs, spatial independent component analysis (sICA) was used to extract FNs from fMRI data acquired from healthy participants performing a visual task with two levels of attention and working memory load. The task-related modulations of extracted FNs were assessed. A group of FNs showed increased activity at low-load conditions and reduced activity at high-load conditions. These FNs together involve the left lateral frontoparietal cortex, insula, and ventromedial prefrontal cortex. A second group of FNs showed increased activity at high-load conditions and reduced activity at low-load conditions. These FNs together involve the intraparietal sulcus, frontal eye field, lateral frontoparietal cortex, insula, and dorsal anterior cingulate, bilaterally. Though the two groups of FNs showed opposite task-related modulations, they overlapped extensively at both the lateral and medial frontoparietal cortex and insula. Such an overlap of FNs would not likely be revealed using standard general-linear-model-based analyses.

CONCLUSIONS

By assessing task-related modulations, this study differentiated the functional roles of overlapping FNs. Several FNs including the left frontoparietal network are implicated in task conditions of low attentional load, while another set of FNs including the dorsal attentional network is implicated in task conditions involving high attentional demands.

Graph independent component analysis reveals repertoires of intrinsic network components in the human brain

Does each cognitive task elicit a new cognitive network each time in the brain? Recent data suggest that pre-existing repertoires of a much smaller number of canonical network components are selectively and dynamically used to compute new cognitive tasks. To this end, we propose a novel method (graph-ICA) that seeks to extract these canonical network components from a limited number of resting state spontaneous networks. Graph-ICA decomposes a weighted mixture of source edge-sharing subnetworks with different weighted edges by applying an independent component analysis on cross-sectional brain networks represented as graphs. We evaluated the plausibility in our simulation study and identified 49 intrinsic subnetworks by applying it in the resting state fMRI data. Using the derived subnetwork repertories, we decomposed brain networks during specific tasks including motor activity, working memory exercises, and verb generation, and identified subnetworks associated with performance on these tasks. We also analyzed sex differences in utilization of subnetworks, which was useful in characterizing group networks. These results suggest that this method can effectively be utilized to identify task-specific as well as sex-specific functional subnetworks. Moreover, graph-ICA can provide more direct information on the edge weights among brain regions working together as a network, which cannot be directly obtained through voxel-level spatial ICA.

Task-related concurrent but opposite modulations of overlapping functional networks as revealed by spatial ICA

Animal studies indicate that different functional networks (FNs), each with a unique timecourse, may overlap at common brain regions. For understanding how different FNs overlap in the human brain and how the timecourses of overlapping FNs are modulated by cognitive tasks, we applied spatial independent component analysis (sICA) to functional magnetic resonance imaging (fMRI) data. These data were acquired from healthy participants while they performed a visual task with parametric loads of attention and working memory. sICA identified a total of 14 FNs, and they showed different extents of overlap at a majority of brain regions exhibiting any functional activity. More FNs overlapped at the higher-order association cortex including the anterior and posterior cingulate, precuneus, insula, and lateral and medial frontoparietal cortices (FPCs) than at the primary sensorimotor cortex. Furthermore, overlapping FNs exhibited concurrent but different task-related modulations of timecourses. FNs showing task-related up- vs. down-modulation of timecourses overlapped at both the lateral and medial FPCs and subcortical structures including the thalamus, striatum, and midbrain ventral tegmental area (VTA). Such task-related, concurrent, but opposite changes in timecourses in the same brain regions may not be detected by current analyses based on General-Linear-Model (GLM). The present findings indicate that multiple cognitive processes may associate with common brain regions and exhibit simultaneous but different modulations in timecourses during cognitive tasks.

Magnetic resonance imaging and graph theoretical analysis of complex brain networks in neuropsychiatric disorders

Neurological and psychiatric disorders disturb higher cognitive functions and are accompanied by aberrant cortico-cortical axonal pathways or synchronizations of neural activity. A large proportion of neuroimaging studies have focused on examining the focal morphological abnormalities of various gray and white matter structures or the functional activities of brain areas during goal-directed tasks or the resting state, which provides vast quantities of information on both the structural and functional alterations in the patients' brain. However, these studies often ignore the interactions among multiple brain regions that constitute complex brain networks underlying higher cognitive function. Information derived from recent advances of noninvasive magnetic resonance imaging (MRI) techniques and computational methodologies such as graph theory have allowed researchers to explore the patterns of structural and functional connectivity of healthy and diseased brains in vivo. In this article, we summarize the recent advances made in the studies of both structural (gray matter morphology and white matter fibers) and functional (synchronized neural activity) brain networks based on human MRI data pertaining to neuropsychiatric disorders. These studies bring a systems-level perspective to the alterations of the topological organization of complex brain networks and the underlying pathophysiological mechanisms. Specifically, noninvasive imaging of structural and functional brain networks and follow-up graph-theoretical analyses demonstrate the potential to establish systems-level biomarkers for clinical diagnosis, progression monitoring, and treatment effects evaluation for neuropsychiatric disorders.

A sliding time-window ICA reveals spatial variability of the default mode network in time

Recent evidence on resting-state networks in functional (connectivity) magnetic resonance imaging (fcMRI) suggests that there may be significant spatial variability of activity foci over time. This study used a sliding time window approach with the spatial domain-independent component analysis (SliTICA) to detect spatial maps of resting-state networks over time. The study hypothesis was that the spatial distribution of a functionally connected network would present marked variability over time. The spatial stability of successive sliding-window maps of the default mode network (DMN) from fcMRI data of 12 participants imaged in the resting state was analyzed. Control measures support previous findings on the stability of independent component analysis in measuring sliding-window sources accurately. The spatial similarity of successive DMN maps varied over time at low frequencies and presented a 1/f power spectral pattern. SliTICA maps show marked temporal variation within the DMN; a single voxel was detected inside a group DMN map in maximally 82% of time windows. Mapping of incidental connectivity reveals centrifugally increasing connectivity to the brain cortex outside the DMN core areas. In conclusion, SliTICA shows marked spatial variance of DMN activity in time, which may offer a more comprehensive measurement of the overall functional activity of a network.

Overlapping replicator dynamics for functional subnetwork identification

Functional magnetic resonance imaging (fMRI) has been widely used for inferring brain regions that tend to work in tandem and grouping them into subnetworks. Despite that certain brain regions are known to interact with multiple subnetworks, few existing techniques support identification of subnetworks with overlaps. To address this limitation, we propose a novel approach based on replicator dynamics that facilitates detection of sparse overlapping subnetworks. We refer to our approach as overlapping replicator dynamics (RDOL). On synthetic data, we show that RDOL achieves higher accuracy in subnetwork identification than state-of-the-art methods. On real data, we demonstrate that RDOL is able to identify major functional hubs that are known to serve as communication channels between brain regions, in addition to detecting commonly observed functional subnetworks. Moreover, we illustrate that knowing the subnetwork overlaps enables inference of functional pathways, e.g. from primary sensory areas to the integration hubs.

Protein-protein interaction network-based detection of functionally similar proteins within species

Although functionally similar proteins across species have been widely studied, functionally similar proteins within species showing low sequence similarity have not been examined in detail. Identification of these proteins is of significant importance for understanding biological functions, evolution of protein families, progression of co-evolution, and convergent evolution and others which cannot be obtained by detection of functionally similar proteins across species. Here, we explored a method of detecting functionally similar proteins within species based on graph theory. After denoting protein-protein interaction networks using graphs, we split the graphs into subgraphs using the 1-hop method. Proteins with functional similarities in a species were detected using a method of modified shortest path to compare these subgraphs and to find the eligible optimal results. Using seven protein-protein interaction networks and this method, some functionally similar proteins with low sequence similarity that cannot detected by sequence alignment were identified. By analyzing the results, we found that, sometimes, it is difficult to separate homologous from convergent evolution. Evaluation of the performance of our method by gene ontology term overlap showed that the precision of our method was excellent.