The effect of model order selection in group PICA

Merging clinical neuropsychology and functional neuroimaging to evaluate the construct validity and neural network engagement of the n-back task

The n-back task is a widely used neuroimaging paradigm for studying the neural basis of working memory (WM); however, its neuropsychometric properties have received little empirical investigation. The present study merged clinical neuropsychology and functional magnetic resonance imaging (fMRI) to explore the construct validity of the letter variant of the n-back task (LNB) and to further identify the task-evoked networks involved in WM. Construct validity of the LNB task was investigated using a bootstrapping approach to correlate LNB task performance across clinically validated neuropsychological measures of WM to establish convergent validity, as well as measures of related but distinct cognitive constructs (i.e., attention and short-term memory) to establish discriminant validity. Independent component analysis (ICA) identified brain networks active during the LNB task in 34 healthy control participants, and general linear modeling determined task-relatedness of these networks. Bootstrap correlation analyses revealed moderate to high correlations among measures expected to converge with LNB (|ρ|≥ 0.37) and weak correlations among measures expected to discriminate (|ρ|≤ 0.29), controlling for age and education. ICA identified 35 independent networks, 17 of which demonstrated engagement significantly related to task condition, controlling for reaction time variability. Of these, the bilateral frontoparietal networks, bilateral dorsolateral prefrontal cortices, bilateral superior parietal lobules including precuneus, and frontoinsular network were preferentially recruited by the 2-back condition compared to 0-back control condition, indicating WM involvement. These results support the use of the LNB as a measure of WM and confirm its use in probing the network-level neural correlates of WM processing.

Isoflurane induces dose-dependent alterations in the cortical connectivity profiles and dynamic properties of the brain's functional architecture

Despite their widespread use, the effect of anesthetic agents on the brain's functional architecture remains poorly understood. This is particularly true of alterations that occur beyond the point of induced unconsciousness. Here, we examined the distributed intrinsic connectivity of macaques across six isoflurane levels using resting-state functional MRI (fMRI) following the loss of consciousness. The results from multiple analysis strategies showed stable functional connectivity (FC) patterns between 1.00% and 1.50% suggesting this as a suitable range for anesthetized nonhuman primate resting-state investigations. Dose-dependent effects were evident at moderate to high dosages showing substantial alteration of the functional topology and a decrease or complete loss of interhemispheric cortical FC strength including that of contralateral homologues. The assessment of dynamic FC patterns revealed that the functional repertoire of brain states is related to anesthesia depth and most strikingly, that the number of state transitions linearly decreases with increased isoflurane dosage. Taken together, the results indicate dose-specific spatial and temporal alterations of FC that occur beyond the typically defined endpoint of consciousness. Future work will be necessary to determine how these findings generalize across anesthetic types and extend to the transition between consciousness and unconsciousness.

The spatial structure of resting state connectivity stability on the scale of minutes

Resting state functional MRI (rsfMRI) connectivity patterns are not temporally stable, but fluctuate in time at scales shorter than most common rest scan durations (5–10 min). Consequently, connectivity patterns for two different portions of the same scan can differ drastically. To better characterize this temporal variability and understand how it is spatially distributed across the brain, we scanned subjects continuously for 60 min, at a temporal resolution of 1 s, while they rested inside the scanner. We then computed connectivity matrices between functionally-defined regions of interest for non-overlapping 1 min windows, and classified connections according to their strength, polarity, and variability. We found that the most stable connections correspond primarily to inter-hemispheric connections between left/right homologous ROIs. However, only 32% of all within-network connections were classified as most stable. This shows that resting state networks have some long-term stability, but confirms the flexible configuration of these networks, particularly those related to higher order cognitive functions. The most variable connections correspond primarily to inter-hemispheric, across-network connections between non-homologous regions in occipital and frontal cortex. Finally we found a series of connections with negative average correlation, but further analyses revealed that such average negative correlations may be related to the removal of CSF signals during pre-processing. Using the same dataset, we also evaluated how similarity of within-subject whole-brain connectivity matrices changes as a function of window duration (used here as a proxy for scan duration). Our results suggest scanning for a minimum of 10 min to optimize within-subject reproducibility of connectivity patterns across the entire brain, rather than a few predefined networks.

Ongoing activity in temporally coherent networks predicts intra-subject fluctuation of response time to sporadic executive control demands

Can ongoing fMRI BOLD signals predict fluctuations in swiftness of a person's response to sporadic cognitive demands? This is an important issue because it clarifies whether intrinsic brain dynamics, for which spatio-temporal patterns are expressed as temporally coherent networks (TCNs), have effects not only on sensory or motor processes, but also on cognitive processes. Predictivity has been affirmed, although to a limited extent. Expecting a predictive effect on executive performance for a wider range of TCNs constituting the cingulo-opercular, fronto-parietal, and default mode networks, we conducted an fMRI study using a version of the color-word Stroop task that was specifically designed to put a higher load on executive control, with the aim of making its fluctuations more detectable. We explored the relationships between the fluctuations in ongoing pre-trial activity in TCNs and the task response time (RT). The results revealed the existence of TCNs in which fluctuations in activity several seconds before the onset of the trial predicted RT fluctuations for the subsequent trial. These TCNs were distributed in the cingulo-opercular and fronto-parietal networks, as well as in perceptual and motor networks. Our results suggest that intrinsic brain dynamics in these networks constitute "cognitive readiness," which plays an active role especially in situations where information for anticipatory attention control is unavailable. Fluctuations in these networks lead to fluctuations in executive control performance.

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.

Thalamus and posterior temporal lobe show greater inter-network connectivity at rest and across sensory paradigms in schizophrenia

Although a number of recent studies have examined functional connectivity at rest, few have assessed differences between connectivity both during rest and across active task paradigms. Therefore, the question of whether cortical connectivity patterns remain stable or change with task engagement continues to be unaddressed. We collected multi-scan fMRI data on healthy controls (N=53) and schizophrenia patients (N=42) during rest and across paradigms arranged hierarchically by sensory load. We measured functional network connectivity among 45 non-artifactual distinct brain networks. Then, we applied a novel analysis to assess cross paradigm connectivity patterns applied to healthy controls and patients with schizophrenia. To detect these patterns, we fit a group by task full factorial ANOVA model to the group average functional network connectivity values. Our approach identified both stable (static effects) and state-based differences (dynamic effects) in brain connectivity providing a better understanding of how individuals' reactions to simple sensory stimuli are conditioned by the context within which they are presented. Our findings suggest that not all group differences observed during rest are detectable in other cognitive states. In addition, the stable differences of heightened connectivity between multiple brain areas with thalamus across tasks underscore the importance of the thalamus as a gateway to sensory input and provide new insight into schizophrenia.

Intrinsic brain activity of cognitively normal older persons resembles more that of patients both with and at risk for Alzheimer's disease than that of healthy younger persons

In Alzheimer's disease (AD), recent findings suggest that amyloid-β (Aβ)-pathology might start 20-30 years before first cognitive symptoms arise. To account for age as most relevant risk factor for sporadic AD, it has been hypothesized that lifespan intrinsic (i.e., ongoing) activity of hetero-modal brain areas with highest levels of functional connectivity triggers Aβ-pathology. This model induces the simple question whether in older persons without any cognitive symptoms intrinsic activity of hetero-modal areas is more similar to that of symptomatic patients with AD or to that of younger healthy persons. We hypothesize that due to advanced age and therefore potential impact of pre-clinical AD, intrinsic activity of older persons resembles more that of patients than that of younger controls. We tested this hypothesis in younger (ca. 25 years) and older healthy persons (ca. 70 years) and patients with mild cognitive impairment and AD-dementia (ca. 70 years) by the use of resting-state functional magnetic resonance imaging, distinct measures of intrinsic brain activity, and different hierarchical clustering approaches. Independently of applied methods and involved areas, healthy older persons' intrinsic brain activity was consistently more alike that of patients than that of younger controls. Our result provides evidence for larger similarity in intrinsic brain activity between healthy older persons and patients with or at-risk for AD than between older and younger ones, suggesting a significant proportion of pre-clinical AD cases in the group of cognitively normal older people. The observed link of aging and AD with intrinsic brain activity supports the view that lifespan intrinsic activity may contribute critically to the pathogenesis of AD.

Restricted Boltzmann machines for neuroimaging: an application in identifying intrinsic networks

Matrix factorization models are the current dominant approach for resolving meaningful data-driven features in neuroimaging data. Among them, independent component analysis (ICA) is arguably the most widely used for identifying functional networks, and its success has led to a number of versatile extensions to group and multimodal data. However there are indications that ICA may have reached a limit in flexibility and representational capacity, as the majority of such extensions are case-driven, custom-made solutions that are still contained within the class of mixture models. In this work, we seek out a principled and naturally extensible approach and consider a probabilistic model known as a restricted Boltzmann machine (RBM). An RBM separates linear factors from functional brain imaging data by fitting a probability distribution model to the data. Importantly, the solution can be used as a building block for more complex (deep) models, making it naturally suitable for hierarchical and multimodal extensions that are not easily captured when using linear factorizations alone. We investigate the capability of RBMs to identify intrinsic networks and compare its performance to that of well-known linear mixture models, in particular ICA. Using synthetic and real task fMRI data, we show that RBMs can be used to identify networks and their temporal activations with accuracy that is equal or greater than that of factorization models. The demonstrated effectiveness of RBMs supports its use as a building block for deeper models, a significant prospect for future neuroimaging research.

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.

Aberrant dependence of default mode/central executive network interactions on anterior insular salience network activity in schizophrenia

In schizophrenia, consistent structural and functional changes have been demonstrated for the insula including aberrant salience processing, which is critical for psychosis. Interactions within and across default mode and central executive network (DMN, CEN) are impaired in schizophrenia. The question arises whether these 2 types of changes are related. Recently, the anterior insula has been demonstrated to control DMN/CEN interactions. We hypothesized that aberrant insula and DMN/CEN activity in schizophrenia is associated with an impaired dependence of DMN/CEN interactions on anterior insular salience network (SN) activity. Eighteen patients with schizophrenia during psychosis and 20 healthy controls were studied by resting-state-fMRI and psychometric examination. High-model-order independent component analysis of fMRI data revealed spatiotemporal patterns of synchronized ongoing blood-oxygenation-level-dependent (BOLD) activity including SN, DMN, and CEN. Scores of functional and time-lagged connectivity across networks' time courses were calculated. Connectivity scores and spatial network maps were compared between groups and related with patients' hallucination and delusion severity. Spatial BOLD-synchronicity was altered in patients' SN, DMN, and CEN, including decreased activity in the right anterior insula (rAI). Patients' functional connectivity between DMN and CEN was increased and related with hallucinations severity. Importantly, patients' time-lagged connectivity between SN and DMN/CEN was reduced, and decreased rAI activity of the SN was associated with both hallucinations and increased functional connectivity between DMN and CEN. Data provide evidence for an aberrant dependence of DMN/CEN interactions on anterior insular SN activity, linking impaired insula, DMN, CEN activity, and psychosis in schizophrenia.

Insular dysfunction within the salience network is associated with severity of symptoms and aberrant inter-network connectivity in major depressive disorder

Major depressive disorder (MDD) is characterized by altered intrinsic functional connectivity within (intra-iFC) intrinsic connectivity networks (ICNs), such as the Default Mode- (DMN), Salience- (SN) and Central Executive Network (CEN). It has been proposed that aberrant switching between DMN-mediated self-referential and CEN-mediated goal-directed cognitive processes might contribute to MDD, possibly explaining patients' difficulties to disengage the processing of self-focused, often negatively biased thoughts. Recently, it has been shown that the right anterior insula (rAI) within the SN is modulating DMN/CEN interactions. Since structural and functional alterations within the AI have been frequently reported in MDD, we hypothesized that aberrant intra-iFC in the SN's rAI is associated with both aberrant iFC between DMN and CEN (inter-iFC) and severity of symptoms in MDD. Twenty-five patients with MDD and 25 healthy controls were assessed using resting-state fMRI (rs-fMRI) and psychometric examination. High-model-order independent component analysis (ICA) of rs-fMRI data was performed to identify ICNs including DMN, SN, and CEN. Intra-iFC within and inter-iFC between distinct subsystems of the DMN, SN, and CEN were calculated, compared between groups and correlated with the severity of symptoms. Patients with MDD showed (1) decreased intra-iFC within the SN's rAI, (2) decreased inter-iFC between the DMN and CEN, and (3) increased inter-iFC between the SN and DMN. Moreover, decreased intra-iFC in the SN's rAI was associated with severity of symptoms and aberrant DMN/CEN interactions, with the latter losing significance after correction for multiple comparisons. Our results provide evidence for a relationship between aberrant intra-iFC in the salience network's rAI, aberrant DMN/CEN interactions and severity of symptoms, suggesting a link between aberrant salience mapping, abnormal coordination of DMN/CEN based cognitive processes and psychopathology in MDD.

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.

Altered resting-state activity in seasonal affective disorder

At present, our knowledge about seasonal affective disorder (SAD) is based mainly up on clinical symptoms, epidemiology, behavioral characteristics and light therapy. Recently developed measures of resting-state functional brain activity might provide neurobiological markers of brain disorders. Studying functional brain activity in SAD could enhance our understanding of its nature and possible treatment strategies. Functional network connectivity (measured using ICA-dual regression), and amplitude of low-frequency fluctuations (ALFF) were measured in 45 antidepressant-free patients (39.78 ± 10.64, 30 ♀, 15 ♂) diagnosed with SAD and compared with age-, gender- and ethnicity-matched healthy controls (HCs) using resting-state functional magnetic resonance imaging. After correcting for Type 1 error at high model orders (inter-RSN correction), SAD patients showed significantly increased functional connectivity in 11 of the 47 identified RSNs. Increased functional connectivity involved RSNs such as visual, sensorimotor, and attentional networks. Moreover, our results revealed that SAD patients compared with HCs showed significant higher ALFF in the visual and right sensorimotor cortex. Abnormally altered functional activity detected in SAD supports previously reported attentional and psychomotor symptoms in patients suffering from SAD. Further studies, particularly under task conditions, are needed in order to specifically investigate cognitive deficits in SAD.

ICA model order selection of task co-activation networks

Independent component analysis (ICA) has become a widely used method for extracting functional networks in the brain during rest and task. Historically, preferred ICA dimensionality has widely varied within the neuroimaging community, but typically varies between 20 and 100 components. This can be problematic when comparing results across multiple studies because of the impact ICA dimensionality has on the topology of its resultant components. Recent studies have demonstrated that ICA can be applied to peak activation coordinates archived in a large neuroimaging database (i.e., BrainMap Database) to yield whole-brain task-based co-activation networks. A strength of applying ICA to BrainMap data is that the vast amount of metadata in BrainMap can be used to quantitatively assess tasks and cognitive processes contributing to each component. In this study, we investigated the effect of model order on the distribution of functional properties across networks as a method for identifying the most informative decompositions of BrainMap-based ICA components. Our findings suggest dimensionality of 20 for low model order ICA to examine large-scale brain networks, and dimensionality of 70 to provide insight into how large-scale networks fractionate into sub-networks. We also provide a functional and organizational assessment of visual, motor, emotion, and interoceptive task co-activation networks as they fractionate from low to high model-orders.

Resting state fMRI reveals a default mode dissociation between retrosplenial and medial prefrontal subnetworks in ASD despite motion scrubbing

In resting state functional magnetic resonance imaging (fMRI) studies of autism spectrum disorders (ASDs) decreased frontal-posterior functional connectivity is a persistent finding. However, the picture of the default mode network (DMN) hypoconnectivity remains incomplete. In addition, the functional connectivity analyses have been shown to be susceptible even to subtle motion. DMN hypoconnectivity in ASD has been specifically called for re-evaluation with stringent motion correction, which we aimed to conduct by so-called scrubbing. A rich set of default mode subnetworks can be obtained with high dimensional group independent component analysis (ICA) which can potentially provide more detailed view of the connectivity alterations. We compared the DMN connectivity in high-functioning adolescents with ASDs to typically developing controls using ICA dual-regression with decompositions from typical to high dimensionality. Dual-regression analysis within DMN subnetworks did not reveal alterations but connectivity between anterior and posterior DMN subnetworks was decreased in ASD. The results were very similar with and without motion scrubbing thus indicating the efficacy of the conventional motion correction methods combined with ICA dual-regression. Specific dissociation between DMN subnetworks was revealed on high ICA dimensionality, where networks centered at the medial prefrontal cortex and retrosplenial cortex showed weakened coupling in adolescents with ASDs compared to typically developing control participants. Generally the results speak for disruption in the anterior-posterior DMN interplay on the network level whereas local functional connectivity in DMN seems relatively unaltered.

Resting-state networks and the functional connectome of the human brain in agenesis of the corpus callosum

The corpus callosum is the largest white matter fiber bundle connecting the two cerebral hemispheres. In this work, we investigate the effect of callosal dysgenesis on functional magnetic resonance imaging (fMRI) resting-state networks and the functional connectome. Since alternate commissural routes between the cerebral hemispheres exist, we hypothesize that bilateral cortical networks can still be maintained in partial or even complete agenesis of the corpus callosum (AgCC). However, since these commissural routes are frequently indirect, requiring polysynaptic pathways, we hypothesize that quantitative measurements of interhemispheric functional connectivity in bilateral networks will be reduced in AgCC compared with matched controls, especially in the most highly interconnected cortical regions that are the hubs of the connectome. Seventeen resting-state networks were extracted from fMRI of 11 subjects with partial or complete AgCC and 11 matched controls. The results show that the qualitative organization of resting-state networks is very similar between controls and AgCC. However, interhemispheric functional connectivity of precuneus, posterior cingulate cortex, and insular-opercular regions was significantly reduced in AgCC. The preserved network organization was confirmed with a connectomic analysis of the resting-state fMRI data, showing five functional modules that are largely consistent across the control and AgCC groups. Hence, the reduction or even complete absence of callosal connectivity does not affect the qualitative organization of bilateral resting-state networks or the modular organization of the functional connectome, although quantitatively reduced functional connectivity can be demonstrated by measurements within bilateral cortical hubs, supporting the hypothesis that indirect polysynaptic pathways are utilized to preserve interhemispheric temporal synchrony.

State-related functional integration and functional segregation brain networks in schizophrenia

Altered topological properties of brain connectivity networks have emerged as important features of schizophrenia. The aim of this study was to investigate how the state-related modulations to graph measures of functional integration and functional segregation brain networks are disrupted in schizophrenia. Firstly, resting state and auditory oddball discrimination (AOD) fMRI data of healthy controls (HCs) and schizophrenia patients (SZs) were decomposed into spatially independent components (ICs) by group independent component analysis (ICA). Then, weighted positive and negative functional integration (inter-component networks) and functional segregation (intra-component networks) brain networks were built in each subject. Subsequently, connectivity strength, clustering coefficient, and global efficiency of all brain networks were statistically compared between groups (HCs and SZs) in each state and between states (rest and AOD) within group. We found that graph measures of negative functional integration brain network and several positive functional segregation brain networks were altered in schizophrenia during AOD task. The metrics of positive functional integration brain network and one positive functional segregation brain network were higher during the resting state than during the AOD task only in HCs. These findings imply that state-related characteristics of both functional integration and functional segregation brain networks are impaired in schizophrenia which provides new insight into the altered brain performance in this brain disorder.

Shifted intrinsic connectivity of central executive and salience network in borderline personality disorder

Borderline personality disorder (BPD) is characterized by “stable instability” of emotions and behavior and their regulation. This emotional and behavioral instability corresponds with a neurocognitive triple network model of psychopathology, which suggests that aberrant emotional saliency and cognitive control is associated with aberrant interaction across three intrinsic connectivity networks [i.e., the salience network (SN), default mode network (DMN), and central executive network (CEN)]. The objective of the current study was to investigate whether and how such triple network intrinsic functional connectivity (iFC) is changed in patients with BPD. We acquired resting-state functional magnetic resonance imaging (rs-fMRI) data from 14 patients with BPD and 16 healthy controls. High-model order independent component analysis was used to extract spatiotemporal patterns of ongoing, coherent blood-oxygen-level-dependent signal fluctuations from rs-fMRI data. Main outcome measures were iFC within networks (intra-iFC) and between networks (i.e., network time course correlation inter-iFC). Aberrant intra-iFC was found in patients’ DMN, SN, and CEN, consistent with previous findings. While patients’ inter-iFC of the CEN was decreased, inter-iFC of the SN was increased. In particular, a balance index reflecting the relationship of CEN- and SN-inter-iFC across networks was strongly shifted from CEN to SN connectivity in patients. Results provide first preliminary evidence for aberrant triple network iFC in BPD. Our data suggest a shift of inter-network iFC from networks involved in cognitive control to those of emotion-related activity in BPD, potentially reflecting the persistent instability of emotion regulation in patients.

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.

Applying independent component analysis to clinical FMRI at 7 t

Increased BOLD sensitivity at 7 T offers the possibility to increase the reliability of fMRI, but ultra-high field is also associated with an increase in artifacts related to head motion, Nyquist ghosting, and parallel imaging reconstruction errors. In this study, the ability of independent component analysis (ICA) to separate activation from these artifacts was assessed in a 7 T study of neurological patients performing chin and hand motor tasks. ICA was able to isolate primary motor activation with negligible contamination by motion effects. The results of General Linear Model (GLM) analysis of these data were, in contrast, heavily contaminated by motion. Secondary motor areas, basal ganglia, and thalamus involvement were apparent in ICA results, but there was low capability to isolate activation in the same brain regions in the GLM analysis, indicating that ICA was more sensitive as well as more specific. A method was developed to simplify the assessment of the large number of independent components. Task-related activation components could be automatically identified via these intuitive and effective features. These findings demonstrate that ICA is a practical and sensitive analysis approach in high field fMRI studies, particularly where motion is evoked. Promising applications of ICA in clinical fMRI include presurgical planning and the study of pathologies affecting subcortical brain areas.

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.

Impact of analysis methods on the reproducibility and reliability of resting-state networks

Though previous examinations of intrinsic resting-state networks (RSNs) in healthy populations have consistently identified several RSNs that represent connectivity patterns evoked by cognitive and sensory tasks, the effects of different analytic approaches on the reliability and reproducibility of these RSNs have yet to be fully explored. Thus, the primary aim of the current study was to investigate the effect of method (independent component analyses [ICA] vs. seed-based analyses) on RSN reproducibility (independent datasets) for ICA and reliability (independent time points) in both methods using functional magnetic resonance imaging. Good to excellent reproducibility was observed in 9 out of 10 commonly identified RSNs, indicating the robustness of these intrinsic fluctuations at the group level. Reliability analyses showed that results were dependent on three main methodological factors: (1) group versus subject-level analyses (group>subject); (2) whether data from different visits were analyzed separately or jointly with ICA (combined>separate ICA); and (3) whether ICA output was used to directly assess reliability or to inform seed-based analyses (seed-based>ICA). These results suggest that variations in the analytic technique have a significant impact on individual reliability measurements, but do not significantly affect the reproducibility or reliability of RSNs at the group level. Further investigation into the effect of the analytic technique on RSN quantification is warranted to increase the utility of RSN analyses in clinical studies.

Spontaneous brain activity predicts learning ability of foreign sounds

Can learning capacity of the human brain be predicted from initial spontaneous functional connectivity (FC) between brain areas involved in a task? We combined task-related functional magnetic resonance imaging (fMRI) and resting-state fMRI (rs-fMRI) before and after training with a Hindi dental-retroflex nonnative contrast. Previous fMRI results were replicated, demonstrating that this learning recruited the left insula/frontal operculum and the left superior parietal lobe, among other areas of the brain. Crucially, resting-state FC (rs-FC) between these two areas at pretraining predicted individual differences in learning outcomes after distributed (Experiment 1) and intensive training (Experiment 2). Furthermore, this rs-FC was reduced at posttraining, a change that may also account for learning. Finally, resting-state network analyses showed that the mechanism underlying this reduction of rs-FC was mainly a transfer in intrinsic activity of the left frontal operculum/anterior insula from the left frontoparietal network to the salience network. Thus, rs-FC may contribute to predict learning ability and to understand how learning modifies the functioning of the brain. The discovery of this correspondence between initial spontaneous brain activity in task-related areas and posttraining performance opens new avenues to find predictors of learning capacities in the brain using task-related fMRI and rs-fMRI combined.

Abnormal functional network connectivity among resting-state networks in children with frontal lobe epilepsy

BACKGROUND AND PURPOSE

Epilepsy is considered a disorder of neural networks. The aims of this study were to assess functional connectivity within resting-state networks and functional network connectivity across resting-state networks by use of resting-state fMRI in children with frontal lobe epilepsy and to relate changes in resting-state networks with neuropsychological function.

MATERIALS AND METHODS

Fifteen patients with frontal lobe epilepsy and normal MR imaging and 14 healthy control subjects were recruited. Spatial independent component analysis was used to identify the resting-state networks, including frontal, attention, default mode network, sensorimotor, visual, and auditory networks. The Z-maps of resting-state networks were compared between patients and control subjects. The relation between abnormal connectivity and neuropsychological function was assessed. Correlations from all pair-wise combinations of independent components were performed for each group and compared between groups.

RESULTS

The frontal network was the only network that showed reduced connectivity in patients relative to control subjects. The remaining 5 networks demonstrated both reduced and increased functional connectivity within resting-state networks in patients. There was a weak association between connectivity in frontal network and executive function (P = .029) and a significant association between sensorimotor network and fine motor function (P = .004). Control subjects had 79 pair-wise independent components that showed significant temporal coherence across all resting-state networks except for default mode network-auditory network. Patients had 66 pairs of independent components that showed significant temporal coherence across all resting-state networks. Group comparison showed reduced functional network connectivity between default mode network-attention, frontal-sensorimotor, and frontal-visual networks and increased functional network connectivity between frontal-attention, default mode network-sensorimotor, and frontal-visual networks in patients relative to control subjects.

CONCLUSIONS

We found abnormal functional connectivity within and across resting-state networks in children with frontal lobe epilepsy. Impairment in functional connectivity was associated with impaired neuropsychological function.

Functional parcellation of the inferior frontal and midcingulate cortices in a flanker-stop-change paradigm

Conflict monitoring and motor inhibition are engaged in the performance of complex tasks. The midcingulate cortex (MCC) has been suggested to detect conflicts, whereas the right inferior frontal cortex (IFC) seems to be of relevance for the inhibition process. The current experiment investigates the neural underpinnings of their interplay via a modified flanker paradigm. Conflict was manipulated by the congruency of flanking stimuli relative to a target (congruent vs. incongruent) and motor inhibition by a within-trial response change of the initiated response (keep response vs. stop-change). We used event-related functional magnetic resonance imaging, decomposition with high model order ICA, and single trial analysis to derive a functional parcellation of the whole-brain data. Results demonstrate the segmentation of the MCC into anterior and posterior subregions, and of the IFC into the pars opercularis, pars triangularis, and pars orbitalis. The pars opercularis and pars triangularis of the right IFC constituted the foundation of inhibition-related networks. With high conflict on incongruent trials, activity in the posterior MCC network, as well as in one right IFC network was observed. Stop-change trials modulated both the MCC as well as networks covering extended parts of the IFC. Whereas conflict processing and inhibition most often are studied separately, this study provides a synopsis of functionally coupled brain regions acting in concert to enable an optimal performance in situations involving interference and inhibition.

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.

Insular Dysfunction Reflects Altered Between-Network Connectivity and Severity of Negative Symptoms in Schizophrenia during Psychotic Remission

Schizophrenia is characterized by aberrant intrinsic functional connectivity (iFC) within and between intrinsic connectivity networks (ICNs), including the Default Mode- (DMN), Salience- (SN), and Central Executive Network (CEN). The anterior insula (AI) of the SN has been demonstrated to modulate DMN/CEN interactions. Recently, we found that the dependence of DMN/CEN interactions on SN's right AI activity is altered in patients with schizophrenia in acute psychosis and related to psychotic symptoms, indicating a link between aberrant AI, DMN, CEN, and psychosis. However, since structural alterations of the insula are also present during psychotic remission and associated with negative symptoms, impaired AI interaction might be relevant even for psychotic remission and corresponding symptoms. Twelve patients with schizophrenia during psychotic remission (SR) and 12 healthy controls were assessed using resting-state fMRI and psychometric examination. High-model-order independent component analysis of fMRI data revealed ICNs including DMN, SN, and CEN. Scores of iFC within (intra-iFC) and between (inter-iFC) distinct subsystems of the DMN, SN, and CEN were calculated, compared between groups and correlated with the severity of symptoms. Intra-iFC was altered in patients' SN, DMN, and CEN, including decreased intra-iFC in the left AI within the SN. Patients' inter-iFC between SN and CEN was increased and correlated with the severity of negative symptoms. Furthermore, decreased intra-iFC of the left AI correlated with both severity of negative symptoms and increased inter-iFC between SN and CEN. Our result provides first evidence for a relationship between AI dysfunction and altered between-network interactions in schizophrenia during psychotic remission, which is related to the severity of negative symptoms. Together with our previous results, data suggest specific SN/DMN/CEN reorganization in schizophrenia with distinct insular pathways for distinct symptom dimensions.

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.

Networks of task co-activations

Recent progress in neuroimaging informatics and meta-analytic techniques has enabled a novel domain of human brain connectomics research that focuses on task-dependent co-activation patterns across behavioral tasks and cognitive domains. Here, we review studies utilizing the BrainMap database to investigate data trends in the activation literature using methods such as meta-analytic connectivity modeling (MACM), connectivity-based parcellation (CPB), and independent component analysis (ICA). We give examples of how these methods are being applied to learn more about the functional connectivity of areas such as the amygdala, the default mode network, and visual area V5. Methods for analyzing the behavioral metadata corresponding to regions of interest and to their intrinsically connected networks are described as a tool for local functional decoding. We finally discuss the relation of observed co-activation connectivity results to resting state connectivity patterns, and provide implications for future work in this domain.

Extracting intrinsic functional networks with feature-based group independent component analysis

There is increasing use of functional imaging data to understand the macro-connectome of the human brain. Of particular interest is the structure and function of intrinsic networks (regions exhibiting temporally coherent activity both at rest and while a task is being performed), which account for a significant portion of the variance in functional MRI data. While networks are typically estimated based on the temporal similarity between regions (based on temporal correlation, clustering methods, or independent component analysis [ICA]), some recent work has suggested that these intrinsic networks can be extracted from the inter-subject covariation among highly distilled features, such as amplitude maps reflecting regions modulated by a task or even coordinates extracted from large meta analytic studies. In this paper our goal was to explicitly compare the networks obtained from a first-level ICA (ICA on the spatio-temporal functional magnetic resonance imaging (fMRI) data) to those from a second-level ICA (i.e., ICA on computed features rather than on the first-level fMRI data). Convergent results from simulations, task-fMRI data, and rest-fMRI data show that the second-level analysis is slightly noisier than the first-level analysis but yields strikingly similar patterns of intrinsic networks (spatial correlations as high as 0.85 for task data and 0.65 for rest data, well above the empirical null) and also preserves the relationship of these networks with other variables such as age (for example, default mode network regions tended to show decreased low frequency power for first-level analyses and decreased loading parameters for second-level analyses). In addition, the best-estimated second-level results are those which are the most strongly reflected in the input feature. In summary, the use of feature-based ICA appears to be a valid tool for extracting intrinsic networks. We believe it will become a useful and important approach in the study of the macro-connectome, particularly in the context of data fusion.

Testing independent component patterns by inter-subject or inter-session consistency

Independent component analysis (ICA) is increasingly used to analyze patterns of spontaneous activity in brain imaging. However, there are hardly any methods for answering the fundamental question: are the obtained components statistically significant? Most methods considering the significance of components either consider group-differences or use arbitrary thresholds with weak statistical justification. In previous work, we proposed a statistically principled method for testing if the coefficients in the mixing matrix are similar in different subjects or sessions. In many applications of ICA, however, we would like to test the reliability of the independent components themselves and not the mixing coefficients. Here, we develop a test for such an inter-subject consistency by extending our previous theory. The test is applicable, for example, to the spatial activity patterns obtained by spatial ICA in resting-state fMRI. We further improve both this and the previously proposed testing method by introducing a new way of correcting for multiple testing, new variants of the clustering method, and a computational approximation which greatly reduces the memory and computation required.

Disrupted correlation between low frequency power and connectivity strength of resting state brain networks in schizophrenia

Altered brain connectivity has emerged as a central feature of schizophrenia. Low frequency oscillations and connectivity strength (CS) of resting state brain networks are altered in patients with schizophrenia (SZs). However, the relationship between these two measures has not yet been studied. Such work may be helpful in understanding the so-called "rich club" organization (i.e. high-CS nodes are more densely connected among themselves than are nodes of a lower CS in the human brain) in healthy controls (HCs) and SZs. Here we present a study of HCs and SZs examining low frequency oscillations and CS by first decomposing resting state fMRI (R-fMRI) data into independent components (ICs) using group independent component analysis (ICA) and computing the low frequency power ratio (LFPR) of each ICA time course. Weighted brain graphs consisting of ICs were built based on correlations between ICA time courses. Positive CS and negative CS of each node in the brain graphs were then examined. The correlations between LFPR and CSs as well as "rich club" coefficients of group mean brain graphs were assessed. Results demonstrate that the LFPR of some ICs were lower in SZs compared to HCs. In addition, LFPR was correlated with positive CS in HCs, but to a lesser extent in SZs. HCs showed higher normalized rich club parameter than SZs. The findings provide new insight into disordered intrinsic brain graphs in schizophrenia.

The spatiospectral characterization of brain networks: fusing concurrent EEG spectra and fMRI maps

Different imaging modalities capture different aspects of brain activity. Functional magnetic resonance imaging (fMRI) reveals intrinsic networks whose BOLD signals have periods from 100 s (0.01 Hz) to about 10s (0.1 Hz). Electroencephalographic (EEG) recordings, in contrast, commonly reflect cortical electrical fluctuations with periods up to 20 ms (50 Hz) or above. We examined the correspondence between intrinsic fMRI and EEG network activity at rest in order to characterize brain networks both spatially (with fMRI) and spectrally (with EEG). Brain networks were separately identified within the concurrently recorded fMRI and EEG at the aggregate group level with group independent component analysis and the association between spatial fMRI and frequency by spatial EEG sources was examined by deconvolving their component time courses. The two modalities are considered linked if the estimated impulse response function (IRF) is significantly non-zero at biologically plausible delays. We found that negative associations were primarily present within two of five alpha components, which highlights the importance of considering multiple alpha sources in EEG-fMRI. Positive associations were primarily present within the lower (e.g. delta and theta) and higher (e.g. upper beta and lower gamma) spectral regions, sometimes within the same fMRI components. Collectively, the results demonstrate a promising approach to characterize brain networks spatially and spectrally, and reveal that positive and negative associations appear within partially distinct regions of the EEG spectrum.

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.

Fully exploratory network independent component analysis of the 1000 functional connectomes database

The 1000 Functional Connectomes Project is a collection of resting-state fMRI datasets from more than 1000 subjects acquired in more than 30 independent studies from around the globe. This large, heterogeneous sample of resting-state data offers the unique opportunity to study the consistencies of resting-state networks at both subject and study level. In extension to the seminal paper by Biswal et al. (2010), where a repeated temporal concatenation group independent component analysis (ICA) approach on reduced subsets (using 20 as a pre-specified number of components) was used due to computational resource limitations, we herein apply Fully Exploratory Network ICA (FENICA) to 1000 single-subject independent component analyses. This, along with the possibility of using datasets of different lengths without truncation, enabled us to benefit from the full dataset available, thereby obtaining 16 networks consistent over the whole group of 1000 subjects. Furthermore, we demonstrated that the most consistent among these networks at both subject and study level matched networks most often reported in the literature, and found additional components emerging in prefrontal and parietal areas. Finally, we identified the influence of scan duration on the number of components as a source of heterogeneity between studies.

Consistency of network modules in resting-state FMRI connectome data

At rest, spontaneous brain activity measured by fMRI is summarized by a number of distinct resting state networks (RSNs) following similar temporal time courses. Such networks have been consistently identified across subjects using spatial ICA (independent component analysis). Moreover, graph theory-based network analyses have also been applied to resting-state fMRI data, identifying similar RSNs, although typically at a coarser spatial resolution. In this work, we examined resting-state fMRI networks from 194 subjects at a voxel-level resolution, and examined the consistency of RSNs across subjects using a metric called scaled inclusivity (SI), which summarizes consistency of modular partitions across networks. Our SI analyses indicated that some RSNs are robust across subjects, comparable to the corresponding RSNs identified by ICA. We also found that some commonly reported RSNs are less consistent across subjects. This is the first direct comparison of RSNs between ICAs and graph-based network analyses at a comparable resolution.

Functional subdivision of group-ICA results of fMRI data collected during cinema viewing

Independent component analysis (ICA) can unravel functional brain networks from functional magnetic resonance imaging (fMRI) data. The number of the estimated components affects both the spatial pattern of the identified networks and their time-course estimates. Here group-ICA was applied at four dimensionalities (10, 20, 40, and 58 components) to fMRI data collected from 15 subjects who viewed a 15-min silent film ("At land" by Maya Deren). We focused on the dorsal attention network, the default-mode network, and the sensorimotor network. The lowest dimensionalities demonstrated most prominent activity within the dorsal attention network, combined with the visual areas, and in the default-mode network; the sensorimotor network only appeared with ICA comprising at least 20 components. The results suggest that even very low-dimensional ICA can unravel the most prominent functionally-connected brain networks. However, increasing the number of components gives a more detailed picture and functionally feasible subdivision of the major networks. These results improve our understanding of the hierarchical subdivision of brain networks during viewing of a movie that provides continuous stimulation embedded in an attention-directing narrative.

Monkey in the middle: why non-human primates are needed to bridge the gap in resting-state investigations

Resting-state investigations based on the evaluation of intrinsic low-frequency fluctuations of the BOLD fMRI signal have been extensively utilized to map the structure and dynamics of large-scale functional network organization in humans. In addition to increasing our knowledge of normal brain connectivity, disruptions of the spontaneous hemodynamic fluctuations have been suggested as possible diagnostic indicators of neurological and psychiatric disease states. Though the non-invasive technique has been received with much acclamation, open questions remain regarding the origin, organization, phylogenesis, as well as the basis of disease-related alterations underlying the signal patterns. Experimental work utilizing animal models, including the use of neurophysiological recordings and pharmacological manipulations, therefore, represents a critical component in the understanding and successful application of resting-state analysis, as it affords a range of experimental manipulations not possible in human subjects. In this article, we review recent rodent and non-human primate studies and based on the examination of the homologous brain architecture propose the latter to be the best-suited model for exploring these unresolved resting-state concerns. Ongoing work examining the correspondence of functional and structural connectivity, state-dependency and the neuronal correlates of the hemodynamic oscillations are discussed. We then consider the potential experiments that will allow insight into different brain states and disease-related network disruptions that can extend the clinical applications of resting-state fMRI (RS-fMRI).

Inferring task-related networks using independent component analysis in magnetoencephalography

A novel framework for analysing task-positive data in magnetoencephalography (MEG) is presented that can identify task-related networks. Techniques that combine beamforming, the Hilbert transform and temporal independent component analysis (ICA) have recently been applied to resting-state MEG data and have been shown to extract resting-state networks similar to those found in fMRI. Here we extend this approach in two ways. First, we systematically investigate optimisation of time-frequency windows for connectivity measurement. This is achieved by estimating the distribution of functional connectivity scores between nodes of known resting-state networks and contrasting it with a distribution of artefactual scores that are entirely due to spatial leakage caused by the inverse problem. We find that functional connectivity, both in the resting-state and during a cognitive task, is best estimated via correlations in the oscillatory envelope in the 8–20 Hz frequency range, temporally down-sampled with windows of 1–4 s. Second, we combine ICA with the general linear model (GLM) to incorporate knowledge of task structure into our connectivity analysis. The combination of ICA with the GLM helps overcome problems of these techniques when used independently: namely, the interpretation and separation of interesting independent components from those that represent noise in ICA and the correction for multiple comparisons when applying the GLM. We demonstrate the approach on a 2-back working memory task and show that this novel analysis framework is able to elucidate the functional networks involved in the task beyond that which is achieved using the GLM alone. We find evidence of localised task-related activity in the area of the hippocampus, which is difficult to detect reliably using standard methods. Task-positive ICA, coupled with the GLM, has the potential to be a powerful tool in the analysis of MEG data.

Functional connectivity tracks clinical deterioration in Alzheimer's disease

While resting state functional connectivity has been shown to decrease in patients with mild and/or moderate Alzheimer's disease, it is not yet known how functional connectivity changes in patients as the disease progresses. Furthermore, it has been noted that the default mode network is not as homogenous as previously assumed and several fractionations of the network have been proposed. Here, we separately investigated the modulation of 3 default mode subnetworks, as identified with group independent component analysis, by comparing Alzheimer's disease patients to healthy controls and by assessing connectivity changes over time. Our results showed decreased connectivity at baseline in patients versus controls in the posterior default mode network, and increased connectivity in the anterior and ventral default mode networks. At follow-up, functional connectivity decreased across all default mode systems in patients. Our results suggest that earlier in the disease, regions of the posterior default mode network start to disengage whereas regions within the anterior and ventral networks enhance their connectivity. However, as the disease progresses, connectivity within all systems eventually deteriorates.

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.

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.

Functional connectivity-based identification of subdivisions of the basal ganglia and thalamus using multilevel independent component analysis of resting state fMRI

This study aimed to identify subunits of the basal ganglia and thalamus and to investigate the functional connectivity among these anatomically segregated subdivisions and the cerebral cortex in healthy subjects. For this purpose, we introduced multilevel independent component analysis (ICA) of the resting-state functional magnetic resonance imaging (fMRI). After applying ICA to the whole brain gray matter, we applied second-level ICA restrictively to the basal ganglia and the thalamus area to identify discrete functional subunits of those regions. As a result, the basal ganglia and the thalamus were parcelled into 31 functional subdivisions according to their temporal activity patterns. The extracted parcels showed functional network connectivity between hemispheres, between subdivisions of the basal ganglia and thalamus, and between the extracted subdivisions and cerebral functional components. Grossly, these findings correspond to cortico-striato-thalamo-cortical circuits in the brain. This study also showed the utility of multilevel ICA of resting state fMRI in brain network research.

Modular Organization of Functional Network Connectivity in Healthy Controls and Patients with Schizophrenia during the Resting State

Neuroimaging studies have shown that functional brain networks composed from select regions of interest have a modular community structure. However, the organization of functional network connectivity (FNC), comprising a purely data-driven network built from spatially independent brain components, is not yet clear. The aim of this study is to explore the modular organization of FNC in both healthy controls (HCs) and patients with schizophrenia (SZs). Resting state functional magnetic resonance imaging data of HCs and SZs were decomposed into independent components (ICs) by group independent component analysis (ICA). Then weighted brain networks (in which nodes are brain components) were built based on correlations between ICA time courses. Clustering coefficients and connectivity strength of the networks were computed. A dynamic branch cutting algorithm was used to identify modules of the FNC in HCs and SZs. Results show stronger connectivity strength and higher clustering coefficient in HCs with more and smaller modules in SZs. In addition, HCs and SZs had some different hubs. Our findings demonstrate altered modular architecture of the FNC in schizophrenia and provide insights into abnormal topological organization of intrinsic brain networks in this mental illness.

ICA-fNORM: Spatial Normalization of fMRI Data Using Intrinsic Group-ICA Networks

A common pre-processing challenge associated with group level fMRI analysis is spatial registration of multiple subjects to a standard space. Spatial normalization, using a reference image such as the Montreal Neurological Institute brain template, is the most common technique currently in use to achieve spatial congruence across multiple subjects. This method corrects for global shape differences preserving regional asymmetries, but does not account for functional differences. We propose a novel approach to co-register task-based fMRI data using resting state group-ICA networks. We posit that these intrinsic networks (INs) can provide to the spatial normalization process with important information about how each individual’s brain is organized functionally. The algorithm is initiated by the extraction of single subject representations of INs using group level independent component analysis (ICA) on resting state fMRI data. In this proof-of-concept work two of the robust, commonly identified, networks are chosen as functional templates. As an estimation step, the relevant INs are utilized to derive a set of normalization parameters for each subject. Finally, the normalization parameters are applied individually to a different set of fMRI data acquired while the subjects performed an auditory oddball task. These normalization parameters, although derived using rest data, generalize successfully to data obtained with a cognitive paradigm for each subject. The improvement in results is verified using two widely applied fMRI analysis methods: the general linear model and ICA. Resulting activation patterns from each analysis method show significant improvements in terms of detection sensitivity and statistical significance at the group level. The results presented in this article provide initial evidence to show that common functional domains from the resting state brain may be used to improve the group statistics of task-fMRI data.

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.

Functional connectivity during resting-state functional MR imaging: study of the correspondence between independent component analysis and region-of-interest-based methods

BACKGROUND AND PURPOSE

The connectivity across brain regions can be evaluated through fMRI either by using ICA or by means of correlation analysis of time courses measured in predefined ROIs. The purpose of this study was to investigate quantitatively the correspondence between the connectivity information provided by the 2 techniques.

MATERIALS AND METHODS

In this study, resting-state fMRI data from 40 healthy participants were independently analyzed by using spatial ICA and ROI-based analysis. To assess the correspondence between the results provided by the 2 methods, for all combinations of ROIs, we compared the time course correlation coefficient with the corresponding "ICA coactivation index."

RESULTS

A strongly significant correspondence of moderate intensity was found for 20 ICA components (r = 0.44, P < .001). Repeating the analysis with 10, 15, 25, 30, 35, and 40 components, we found that the correlation remained but was weaker (r = 0.35-0.41).

CONCLUSIONS

There is a significant but not complete correspondence between the results provided by ICA and ROI-based analysis of resting-state data.