A simple and objective method for reproducible resting state network (RSN) detection in fMRI

Spatial Independent Component Analysis (ICA) decomposes the time by space functional MRI (fMRI) matrix into a set of 1-D basis time courses and their associated 3-D spatial maps that are optimized for mutual independence. When applied to resting state fMRI (rsfMRI), ICA produces several spatial independent components (ICs) that seem to have biological relevance - the so-called resting state networks (RSNs). The ICA problem is well posed when the true data generating process follows a linear mixture of ICs model in terms of the identifiability of the mixing matrix. However, the contrast function used for promoting mutual independence in ICA is dependent on the finite amount of observed data and is potentially non-convex with multiple local minima. Hence, each run of ICA could produce potentially different IC estimates even for the same data. One technique to deal with this run-to-run variability of ICA was proposed by [1] in their algorithm RAICAR which allows for the selection of only those ICs that have a high run-to-run reproducibility. We propose an enhancement to the original RAICAR algorithm that enables us to assign reproducibility p-values to each IC and allows for an objective assessment of both within subject and across subjects reproducibility. We call the resulting algorithm RAICAR-N (N stands for null hypothesis test), and we have applied it to publicly available human rsfMRI data (http://www.nitrc.org). Our reproducibility analyses indicated that many of the published RSNs in rsfMRI literature are highly reproducible. However, we found several other RSNs that are highly reproducible but not frequently listed in the literature.     

Altered topological properties of functional network connectivity in schizophrenia during resting state: a small-world brain network study

Aberrant topological properties of small-world human brain networks in patients with schizophrenia (SZ) have been documented in previous neuroimaging studies. Aberrant functional network connectivity (FNC, temporal relationships among independent component time courses) has also been found in SZ by a previous resting state functional magnetic resonance imaging (fMRI) study. However, no study has yet determined if topological properties of FNC are also altered in SZ. In this study, small-world network metrics of FNC during the resting state were examined in both healthy controls (HCs) and SZ subjects. FMRI data were obtained from 19 HCs and 19 SZ. Brain images were decomposed into independent components (ICs) by group independent component analysis (ICA). FNC maps were constructed via a partial correlation analysis of ICA time courses. A set of undirected graphs were built by thresholding the FNC maps and the small-world network metrics of these maps were evaluated. Our results demonstrated significantly altered topological properties of FNC in SZ relative to controls. In addition, topological measures of many ICs involving frontal, parietal, occipital and cerebellar areas were altered in SZ relative to controls. Specifically, topological measures of whole network and specific components in SZ were correlated with scores on the negative symptom scale of the Positive and Negative Symptom Scale (PANSS). These findings suggest that aberrant architecture of small-world brain topology in SZ consists of ICA temporally coherent brain networks.     

Default Mode Network,Motor Network,Dorsal and Ventral Basal Ganglia Networks in the Rat Brain: Comparison to Human Networks Using Resting State-fMRI

Rodent models are developed to enhance understanding of the underlying biology of different brain disorders. However, before interpreting findings from animal models in a translational aspect to understand human disease, a fundamental step is to first have knowledge of similarities and differences of the biological systems studied. In this study, we analyzed and verified four known networks termed: default mode network, motor network, dorsal basal ganglia network, and ventral basal ganglia network using resting state functional MRI (rsfMRI) in humans and rats. Our work supports the notion that humans and rats have common robust resting state brain networks and that rsfMRI can be used as a translational tool when validating animal models of brain disorders. In the future, rsfMRI may be used, in addition to short-term interventions, to characterize longitudinal effects on functional brain networks after long-term intervention in humans and rats.     

Resting State fMRI Reveals Diminished Functional Connectivity in a Mouse Model of Amyloidosis

Introduction

Functional connectivity (FC) studies have gained immense popularity in the evaluation of several neurological disorders, such as Alzheimer’s disease (AD). AD is a complex disorder, characterised by several pathological features. The problem with FC studies in patients is that it is not straightforward to focus on a specific aspect of pathology. In the current study, resting state functional magnetic resonance imaging (rsfMRI) is applied in a mouse model of amyloidosis to assess the effects of amyloid pathology on FC in the mouse brain.

Methods

Nine APP/PS1 transgenic and nine wild-type mice (average age 18.9 months) were imaged on a 7T MRI system. The mice were anesthetized with medetomidine and rsfMRI data were acquired using a gradient echo EPI sequence. The data were analysed using a whole brain seed correlation analysis and interhemispheric FC was evaluated using a pairwise seed analysis. Qualitative histological analyses were performed to assess amyloid pathology, inflammation and synaptic deficits.

Results

The whole brain seed analysis revealed an overall decrease in FC in the brains of transgenic mice compared to wild-type mice. The results showed that interhemispheric FC was relatively preserved in the motor cortex of the transgenic mice, but decreased in the somatosensory cortex and the hippocampus when compared to the wild-type mice. The pairwise seed analysis confirmed these results. Histological analyses confirmed the presence of amyloid pathology, inflammation and synaptic deficits in the transgenic mice.

Conclusions

In the current study, rsfMRI demonstrated decreased FC in APP/PS1 transgenic mice compared to wild-type mice in several brain regions. The APP/PS1 transgenic mice had advanced amyloid pathology across the brain, as well as inflammation and synaptic deficits surrounding the amyloid plaques. Future studies should longitudinally evaluate APP/PS1 transgenic mice and correlate the rsfMRI findings to specific stages of amyloid pathology.     

Robust Data Driven Model Order Estimation for Independent Component Analysis of fMRI Data with Low Contrast to Noise

Independent component analysis (ICA) has been successfully utilized for analysis of functional MRI (fMRI) data for task related as well as resting state studies. Although it holds the promise of becoming an unbiased data-driven analysis technique, a few choices have to be made prior to performing ICA, selection of a method for determining the number of independent components (nIC) being one of them. Choice of nIC has been shown to influence the ICA maps, and various approaches (mostly relying on information theoretic criteria) have been proposed and implemented in commonly used ICA analysis packages, such as MELODIC and GIFT. However, there has been no consensus on the optimal method for nIC selection, and many studies utilize arbitrarily chosen values for nIC. Accurate and reliable determination of true nIC is especially important in the setting where the signals of interest contribute only a small fraction of the total variance, i.e. very low contrast-to-noise ratio (CNR), and/or very focal response. In this study, we evaluate the performance of different model order selection criteria and demonstrate that the model order selected based upon bootstrap stability of principal components yields more reliable and accurate estimates of model order. We then demonstrate the utility of this fully data-driven approach to detect weak and focal stimulus-driven responses in real data. Finally, we compare the performance of different multi-run ICA approaches using pseudo-real data.     

Resting-state quantitative electroencephalography reveals increased neurophysiologic connectivity in depression

Symptoms of Major Depressive Disorder (MDD) are hypothesized to arise from dysfunction in brain networks linking the limbic system and cortical regions. Alterations in brain functional cortical connectivity in resting-state networks have been detected with functional imaging techniques, but neurophysiologic connectivity measures have not been systematically examined. We used weighted network analysis to examine resting state functional connectivity as measured by quantitative electroencephalographic (qEEG) coherence in 121 unmedicated subjects with MDD and 37 healthy controls. Subjects with MDD had significantly higher overall coherence as compared to controls in the delta (0.5-4 Hz), theta (4-8 Hz), alpha (8-12 Hz), and beta (12-20 Hz) frequency bands. The frontopolar region contained the greatest number of "hub nodes" (surface recording locations) with high connectivity. MDD subjects expressed higher theta and alpha coherence primarily in longer distance connections between frontopolar and temporal or parietooccipital regions, and higher beta coherence primarily in connections within and between electrodes overlying the dorsolateral prefrontal cortical (DLPFC) or temporal regions. Nearest centroid analysis indicated that MDD subjects were best characterized by six alpha band connections primarily involving the prefrontal region. The present findings indicate a loss of selectivity in resting functional connectivity in MDD. The overall greater coherence observed in depressed subjects establishes a new context for the interpretation of previous studies showing differences in frontal alpha power and synchrony between subjects with MDD and normal controls. These results can inform the development of qEEG state and trait biomarkers for MDD.     

Independent EEG sources are dipolar

Independent component analysis (ICA) and blind source separation (BSS) methods are increasingly used to separate individual brain and non-brain source signals mixed by volume conduction in electroencephalographic (EEG) and other electrophysiological recordings. We compared results of decomposing thirteen 71-channel human scalp EEG datasets by 22 ICA and BSS algorithms, assessing the pairwise mutual information (PMI) in scalp channel pairs, the remaining PMI in component pairs, the overall mutual information reduction (MIR) effected by each decomposition, and decomposition 'dipolarity' defined as the number of component scalp maps matching the projection of a single equivalent dipole with less than a given residual variance. The least well-performing algorithm was principal component analysis (PCA); best performing were AMICA and other likelihood/mutual information based ICA methods. Though these and other commonly-used decomposition methods returned many similar components, across 18 ICA/BSS algorithms mean dipolarity varied linearly with both MIR and with PMI remaining between the resulting component time courses, a result compatible with an interpretation of many maximally independent EEG components as being volume-conducted projections of partially-synchronous local cortical field activity within single compact cortical domains. To encourage further method comparisons, the data and software used to prepare the results have been made available (http://sccn.ucsd.edu/wiki/BSSComparison).     

Evaluating the effective connectivity of resting state networks using conditional Granger causality

The human brain has been documented to be spatially organized in a finite set of specific coherent patterns, namely resting state networks (RSNs). The interactions among RSNs, being potentially dynamic and directional, may not be adequately captured by simple correlation or anticorrelation. In order to evaluate the possible effective connectivity within those RSNs, we applied a conditional Granger causality analysis (CGCA) to the RSNs retrieved by independent component analysis (ICA) from resting state functional magnetic resonance imaging (fMRI) data. Our analysis provided evidence for specific causal influences among the detected RSNs: default-mode, dorsal attention, core, central-executive, self-referential, somatosensory, visual, and auditory networks. In particular, we identified that self-referential and default-mode networks (DMNs) play distinct and crucial roles in the human brain functional architecture. Specifically, the former RSN exerted the strongest causal influence over the other RSNs, revealing a top-down modulation of self-referential mental activity (SRN) over sensory and cognitive processing. In quite contrast, the latter RSN was profoundly affected by the other RSNs, which may underlie an integration of information from primary function and higher level cognition networks, consistent with previous task-related studies. Overall, our results revealed the causal influences among these RSNs at different processing levels, and supplied information for a deeper understanding of the brain network dynamics.     

A Multi-Methodological MR Resting State Network Analysis to Assess the Changes in Brain Physiology of Children with ADHD

The purpose of this work was to highlight the neurological differences between the MR resting state networks of a group of children with ADHD (pre-treatment) and an age-matched healthy group. Results were obtained using different image analysis techniques. A sample of n = 46 children with ages between 6 and 12 years were included in this study (23 per cohort). Resting state image analysis was performed using ReHo, ALFF and ICA techniques. ReHo and ICA represent connectivity analyses calculated with different mathematical approaches. ALFF represents an indirect measurement of brain activity. The ReHo and ICA analyses suggested differences between the two groups, while the ALFF analysis did not. The ReHo and ALFF analyses presented differences with respect to the results previously reported in the literature. ICA analysis showed that the same resting state networks that appear in healthy volunteers of adult age were obtained for both groups. In contrast, these networks were not identical when comparing the healthy and ADHD groups. These differences affected areas for all the networks except the Right Memory Function network. All techniques employed in this study were used to monitor different cerebral regions which participate in the phenomenological characterization of ADHD patients when compared to healthy controls. Results from our three analyses indicated that the cerebellum and mid-frontal lobe bilaterally for ReHo, the executive function regions in ICA, and the precuneus, cuneus and the clacarine fissure for ALFF, were the “hubs” in which the main inter-group differences were found. These results do not just help to explain the physiology underlying the disorder but open the door to future uses of these methodologies to monitor and evaluate patients with ADHD.     

Default Mode,Dorsal Attention and Auditory Resting State Networks Exhibit Differential Functional Connectivity in Tinnitus and Hearing Loss

We investigated auditory, dorsal attention, and default mode networks in adults with tinnitus and hearing loss in a resting state functional connectivity study. Data were obtained using continuous functional magnetic resonance imaging (fMRI) while the participants were at “rest” and were not performing any task. Participants belonged to one of three groups: middle-aged adults with tinnitus and mild-to-moderate high frequency hearing loss (TIN), age-matched controls with normal hearing and no tinnitus (NH), and a second control group with mild-to-moderate high frequency hearing loss without tinnitus (HL). After standard preprocessing, (a) a group independent component analysis (ICA) using 30 components and (b) a seeding-based connectivity analysis were conducted. In the group ICA, the default mode network was the only network to display visual differences between subject groups. In the seeding analysis, we found increased connectivity between the left parahippocampus and the auditory resting state network in the TIN group when compared to NH controls. Similarly, there was also an increased correlation between the right parahippocampus and the dorsal attention network when compared to HL controls. Other group differences in this attention network included decreased correlations between the seed regions and the right supramarginal gyrus in TIN patients when compared to HL controls. In the default mode network, there was a strong decrease in correlation between the seed regions and the precuneus when compared to both control groups. The findings of this study identify specific alterations in the connectivity of the default mode, dorsal attention, and auditory resting state networks due to tinnitus. The results suggest that therapies for tinnitus that mitigate the increased connectivity of limbic regions with auditory and attention resting state networks and the decreased coherence of the default mode network could be effective at reducing tinnitus-related distress.     

Robust reproducible resting state networks in the awake rodent brain

Resting state networks (RSNs) have been studied extensively with functional MRI in humans in health and disease to reflect brain function in the un-stimulated state as well as reveal how the brain is altered with disease. Rodent models of disease have been used comprehensively to understand the biology of the disease as well as in the development of new therapies. RSN reported studies in rodents, however, are few, and most studies are performed with anesthetized rodents that might alter networks and differ from their non-anesthetized state. Acquiring RSN data in the awake rodent avoids the issues of anesthesia effects on brain function. Using high field fMRI we determined RSNs in awake rats using an independent component analysis (ICA) approach, however, ICA analysis can produce a large number of components, some with biological relevance (networks). We further have applied a novel method to determine networks that are robust and reproducible among all the components found with ICA. This analysis indicates that 7 networks are robust and reproducible in the rat and their putative role is discussed.     

Changes in community structure of resting state functional connectivity in unipolar depression

Major depression is a prevalent disorder that imposes a significant burden on society, yet objective laboratory-style tests to assist in diagnosis are lacking. We employed network-based analyses of "resting state" functional neuroimaging data to ascertain group differences in the endogenous cortical activity between healthy and depressed subjects.We additionally sought to use machine learning techniques to explore the ability of these network-based measures of resting state activity to provide diagnostic information for depression. Resting state fMRI data were acquired from twenty two depressed outpatients and twenty two healthy subjects matched for age and gender. These data were anatomically parcellated and functional connectivity matrices were then derived using the linear correlations between the BOLD signal fluctuations of all pairs of cortical and subcortical regions.We characterised the hierarchical organization of these matrices using network-based matrics, with an emphasis on their mid-scale "modularity" arrangement. Whilst whole brain measures of organization did not differ between groups, a significant rearrangement of their community structure was observed. Furthermore we were able to classify individuals with a high level of accuracy using a support vector machine, primarily through the use of a modularity-based metric known as the participation index.In conclusion, the application of machine learning techniques to features of resting state fMRI network activity shows promising potential to assist in the diagnosis of major depression, now suggesting the need for validation in independent data sets.     

Iterative Cross-Correlation Analysis of Resting State Functional Magnetic Resonance Imaging Data

Seed-based cross-correlation analysis (sCCA) and independent component analysis have been widely employed to extract functional networks from the resting state functional magnetic resonance imaging data. However, the results of sCCA, in terms of both connectivity strength and network topology, can be sensitive to seed selection variations. ICA avoids the potential problems due to seed selection, but choosing which component(s) to represent the network of interest could be subjective and problematic. In this study, we proposed a seed-based iterative cross-correlation analysis (siCCA) method for resting state brain network analysis. The method was applied to extract default mode network (DMN) and stable task control network (STCN) in two independent datasets acquired from normal adults. Compared with the networks obtained by traditional sCCA and ICA, the resting state networks produced by siCCA were found to be highly stable and independent on seed selection. siCCA was used to analyze DMN in first-episode major depressive disorder (MDD) patients. It was found that, in the MDD patients, the volume of DMN negatively correlated with the patients'' social disability screening schedule scores.     

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.     

Resting‐state networks of the neonate brain identified using independent component analysis

Resting‐state functional magnetic resonance imaging (rs‐fMRI) has been successfully used to probe the intrinsic functional organization of the brain and to study brain development. Here, we implemented a combination of individual and group independent component analysis (ICA) of FSL on a 6‐min resting‐state data set acquired from 21 naturally sleeping term‐born (age 26 ± 6.7 d), healthy neonates to investigate the emerging functional resting‐state networks (RSNs). In line with the previous literature, we found evidence of sensorimotor, auditory/language, visual, cerebellar, thalmic, parietal, prefrontal, anterior cingulate as well as dorsal and ventral aspects of the default‐mode‐network. Additionally, we identified RSNs in frontal, parietal, and temporal regions that have not been previously described in this age group and correspond to the canonical RSNs established in adults. Importantly, we found that careful ICA‐based denoising of fMRI data increased the number of networks identified with group‐ICA, whereas the degree of spatial smoothing did not change the number of identified networks. Our results show that the infant brain has an established set of RSNs soon after birth.     

Analysis of EEG characteristic at early stages of depression using method of independent components

Independent Component Analysis (ICA) was used for 19-channel resting EEG analysis 111 patients at early stages of depressive disorder and 526 age-matched healthy subjects. Comparison of independent components power spectra in depressed patients and healthy subjects in two states: Eyes closed and Eyes open, has revealed significant differences between groups for three frequency bands: Theta (4-7.5 Hz), Alpha (7.5-14 Hz), and Beta (14-20 Hz). Increased power of alpha and theta activity in depressed patients at parietal and occipital sites may be caused by decreased cortical activation of these regions. Diffuse enhancement of beta activity level can correlate with anxiety symptoms which take an important place in clinical picture of depressive disorder at early stages. Using of ICA method for comparison of spectral characteristics of EEG in groups of patients with different brain pathology and healthy subjects gives a possibility to localize more precisely the discovered differences as compare to traditional analysis of EEG spectra.     

Spatiotemporal Segregation of Neural Response to Auditory Stimulation: An fMRI Study Using Independent Component Analysis and Frequency-Domain Analysis

Although auditory processing has been widely studied with conventional parametric methods, there have been a limited number of independent component analysis (ICA) applications in this area. The purpose of this study was to examine spatiotemporal behavior of brain networks in response to passive auditory stimulation using ICA. Continuous broadband noise was presented binaurally to 19 subjects with normal hearing. ICA was performed to segregate spatial networks, which were subsequently classified according to their temporal relation to the stimulus using power spectrum analysis. Classification of separated networks resulted in 3 stimulus-activated, 9 stimulus-deactivated, 2 stimulus-neutral (stimulus-dependent but not correlated with the stimulation timing), and 2 stimulus-unrelated (fluctuations that did not follow the stimulus cycles) components. As a result of such classification, spatiotemporal subdivisions were observed in a number of cortical structures, namely auditory, cingulate, and sensorimotor cortices, where parts of the same cortical network responded to the stimulus with different temporal patterns. The majority of the classified networks seemed to comprise subparts of the known resting-state networks (RSNs); however, they displayed different temporal behavior in response to the auditory stimulus, indicating stimulus-dependent temporal segregation of RSNs. Only one of nine deactivated networks coincided with the “classic” default-mode network, suggesting the existence of a stimulus-dependent default-mode network, different from that commonly accepted.     

Cortical network dynamics with time delays reveals functional connectivity in the resting brain

In absence of all goal-directed behavior, a characteristic network of cortical regions involving prefrontal and cingulate cortices consistently shows temporally coherent fluctuations. The origin of these fluctuations is unknown, but has been hypothesized to be of stochastic nature. In the present paper we test the hypothesis that time delays in the network dynamics play a crucial role in the generation of these fluctuations. By tuning the propagation velocity in a network based on primate connectivity, we scale the time delays and demonstrate the emergence of the resting state networks for biophysically realistic parameters.     

Anatomical connectivity defines the organization of clusters of cortical areas in the macaque monkey and the cat

The number of different cortical structures in mammalian brains and the number of extrinsic fibres linking these regions are both large. As with any complex system, systematic analysis is required to draw reliable conclusions about the organization of the complex neural networks comprising these numerous elements. One aspect of organization that has long been suspected is that cortical networks are organized into 'streams' or 'systems'. Here we report computational analyses capable of showing whether clusters of strongly interconnected areas are aspects of the global organization of cortical systems in macaque and cat. We used two different approaches to analyse compilations of corticocortical connection data from the macaque and the cat. The first approach, optimal set analysis, employed an explicit definition of a neural 'system' or 'stream', which was based on differential connectivity. We defined a two-component cost function that described the cost of the global cluster arrangement of areas in terms of the areas' connectivity within and between candidate clusters. Optimal cluster arrangements of cortical areas were then selected computationally from the very many possible arrangements, using an evolutionary optimization algorithm. The second approach, non-parametric cluster analysis (NPCA), grouped cortical areas on the basis of their proximity in multidimensional scaling representations. We used non-metric multidimensional scaling to represent the cortical connectivity structures metrically in two and five dimensions. NPCA then analysed these representations to determine the nature of the clusters for a wide range of different cluster shape parameters. The results from both approaches largely agreed. They showed that macaque and cat cortices are organized into densely intra-connected clusters of areas, and identified the constituent members of the clusters. These clusters reflected functionally specialized sets of cortical areas, suggesting that structure and function are closely linked at this gross, systems level.