Groupwise whole-brain parcellation from resting-state fMRI data for network node identification

Topographic organization of the cerebral cortex and brain cartography

One of the most specific but also challenging properties of the brain is its topographic organization into distinct modules or cortical areas. In this paper, we first review the concept of topographic organization and its historical development. Next, we provide a critical discussion of the current definition of what constitutes a cortical area, why the concept has been so central to the field of neuroimaging and the challenges that arise from this view. A key aspect in this discussion is the issue of spatial scale and hierarchy in the brain. Focusing on in-vivo brain parcellation as a rapidly expanding field of research, we highlight potential limitations of the classical concept of cortical areas in the context of multi-modal parcellation and propose a revised interpretation of cortical areas building on the concept of neurobiological atoms that may be aggregated into larger units within and across modalities. We conclude by presenting an outlook on the implication of this revised concept for future mapping studies and raise some open questions in the context of brain parcellation.

Longitudinal test-retest neuroimaging data from healthy young adults in southwest China

Multimodal magnetic resonance imaging (mMRI) has been widely used to map the structure and function of the human brain, as well as its behavioral associations. However, to date, a large sample with a long-term longitudinal design and a narrow age-span has been lacking for the assessment of test-retest reliability and reproducibility of brain-behavior correlations, as well as the development of novel causal insights into these correlational findings. Here we describe the SLIM dataset, which includes brain and behavioral data across a long-term retest-duration within three and a half years, mMRI scans provided a set of structural, diffusion and resting-state functional MRI images, along with rich samples of behavioral assessments addressed-demographic, cognitive and emotional information. Together with the Consortium for Reliability and Reproducibility (CoRR), the SLIM is expected to accelerate the reproducible sciences of the human brain by providing an open resource for brain-behavior discovery sciences with big-data approaches.

Using connectome-based predictive modeling to predict individual behavior from brain connectivity

Neuroimaging is a fast-developing research area in which anatomical and functional images of human brains are collected using techniques such as functional magnetic resonance imaging (fMRI), diffusion tensor imaging (DTI), and electroencephalography (EEG). Technical advances and large-scale data sets have allowed for the development of models capable of predicting individual differences in traits and behavior using brain connectivity measures derived from neuroimaging data. Here, we present connectome-based predictive modeling (CPM), a data-driven protocol for developing predictive models of brain-behavior relationships from connectivity data using cross-validation. This protocol includes the following steps: (i) feature selection, (ii) feature summarization, (iii) model building, and (iv) assessment of prediction significance. We also include suggestions for visualizing the most predictive features (i.e., brain connections). The final result should be a generalizable model that takes brain connectivity data as input and generates predictions of behavioral measures in novel subjects, accounting for a considerable amount of the variance in these measures. It has been demonstrated that the CPM protocol performs as well as or better than many of the existing approaches in brain-behavior prediction. As CPM focuses on linear modeling and a purely data-driven approach, neuroscientists with limited or no experience in machine learning or optimization will find it easy to implement these protocols. Depending on the volume of data to be processed, the protocol can take 10-100 min for model building, 1-48 h for permutation testing, and 10-20 min for visualization of results.

Connectivity-based parcellation reveals distinct cortico-striatal connectivity fingerprints in Autism Spectrum Disorder

Autism Spectrum Disorder (ASD) has been associated with abnormal synaptic development causing a breakdown in functional connectivity. However, when measured at the macro scale using resting state fMRI, these alterations are subtle and often difficult to detect due to the large heterogeneity of the pathology. Recently, we outlined a novel approach for generating robust biomarkers of resting state functional magnetic resonance imaging (RS-fMRI) using connectivity based parcellation of gross morphological structures to improve single-subject reproducibility and generate more robust connectivity fingerprints. Here we apply this novel approach to investigating the organization and connectivity strength of the cortico-striatal system in a large sample of ASD individuals and typically developed (TD) controls (N=130 per group). Our results showed differences in the parcellation of the striatum in ASD. Specifically, the putamen was found to be one single structure in ASD, whereas this was split into anterior and posterior segments in an age, IQ, and head movement matched TD group. An analysis of the connectivity fingerprints revealed that the group differences in clustering were driven by differential connectivity between striatum and the supplementary motor area, posterior cingulate cortex, and posterior insula. Our approach for analysing RS-fMRI in clinical populations has provided clear evidence that cortico-striatal circuits are organized differently in ASD. Based on previous task-based segmentations of the striatum, we believe that the anterior putamen cluster present in TD, but not in ASD, likely contributes to social and language processes.

Mapping the functional connectome traits of levels of consciousness

Examining task-free functional connectivity (FC) in the human brain offers insights on how spontaneous integration and segregation of information relate to human cognition, and how this organization may be altered in different conditions, and neurological disorders. This is particularly relevant for patients in disorders of consciousness (DOC) following severe acquired brain damage and coma, one of the most devastating conditions in modern medical care. We present a novel data-driven methodology, connICA, which implements Independent Component Analysis (ICA) for the extraction of robust independent FC patterns (FC-traits) from a set of individual functional connectomes, without imposing any a priori data stratification into groups. We here apply connICA to investigate associations between network traits derived from task-free FC and cognitive/clinical features that define levels of consciousness. Three main independent FC-traits were identified and linked to consciousness-related clinical features. The first one represents the functional configuration of a "resting" human brain, and it is associated to a sedative (sevoflurane), the overall effect of the pathology and the level of arousal. The second FC-trait reflects the disconnection of the visual and sensory-motor connectivity patterns. It also relates to the time since the insult and to the ability of communicating with the external environment. The third FC-trait isolates the connectivity pattern encompassing the fronto-parietal and the default-mode network areas as well as the interaction between left and right hemispheres, which are also associated to the awareness of the self and its surroundings. Each FC-trait represents a distinct functional process with a role in the degradation of conscious states of functional brain networks, shedding further light on the functional sub-circuits that get disrupted in severe brain-damage.

Correspondent Functional Topography of the Human Left Inferior Parietal Lobule at Rest and Under Task Revealed Using Resting-State fMRI and Coactivation Based Parcellation

The human left inferior parietal lobule (LIPL) plays a pivotal role in many cognitive functions and is an important node in the default mode network (DMN). Although many previous studies have proposed different parcellation schemes for the LIPL, the detailed functional organization of the LIPL and the exact correspondence between the DMN and LIPL subregions remain unclear. Mounting evidence indicates that spontaneous fluctuations in the brain are strongly associated with cognitive performance at the behavioral level. However, whether a consistent functional topographic organization of the LIPL during rest and under task can be revealed remains unknown. Here, they used resting-state functional connectivity (RSFC) and task-related coactivation patterns separately to parcellate the LIPL and identified seven subregions. Four subregions were located in the supramarginal gyrus (SMG) and three subregions were located in the angular gyrus (AG). The subregion-specific networks and functional characterization revealed that the four anterior subregions were found to be primarily involved in sensorimotor processing, movement imagination and inhibitory control, audition perception and speech processing, and social cognition, whereas the three posterior subregions were mainly involved in episodic memory, semantic processing, and spatial cognition. The results revealed a detailed functional organization of the LIPL and suggested that the LIPL is a functionally heterogeneous area. In addition, the present study demonstrated that the functional architecture of the LIPL during rest corresponds with that found in task processing. Hum Brain Mapp 38:1659-1675, 2017. © 2017 Wiley Periodicals, Inc.

Biological Relevance of Network Architecture

Mathematical representations of brain networks in neuroscience through the use of graph theory may be very useful for the understanding of neurological diseases and disorders and such an explanatory power is currently under intense investigation. Graph metrics are expected to vary across subjects and are likely to reflect behavioural and cognitive performances. The challenge is to set up a framework that can explain how behaviour, cognition, memory, and other brain properties can emerge through the combined interactions of neurons, ensembles of neurons, and larger-scale brain regions that make information transfer possible. "Hidden" graph theoretic properties in the construction of brain networks may limit or enhance brain functionality and may be representative of aspects of human psychology. As theorems emerge from simple mathematical properties of graphs, similarly, cognition and behaviour may emerge from the molecular, cellular and brain region substrate interactions. In this review report, we identify some studies in the current literature that have used graph theoretical metrics to extract neurobiological conclusions, we briefly discuss the link with the human connectome project as an effort to integrate human data that may aid the study of emergent patterns and we suggest a way to start categorizing diseases according to their brain network pathologies as these are measured by graph theory.

Resting-State Functional Connectivity in the Human Connectome Project: Current Status and Relevance to Understanding Psychopathology

A key tenet of modern psychiatry is that psychiatric disorders arise from abnormalities in brain circuits that support human behavior. Our ability to examine hypotheses around circuit-level abnormalities in psychiatric disorders has been made possible by advances in human neuroimaging technologies. These advances have provided the basis for recent efforts to develop a more complex understanding of the function of brain circuits in health and of their relationship to behavior-providing, in turn, a foundation for our understanding of how disruptions in such circuits contribute to the development of psychiatric disorders. This review focuses on the use of resting-state functional connectivity MRI to assess brain circuits, on the advances generated by the Human Connectome Project, and on how these advances potentially contribute to understanding neural circuit dysfunction in psychopathology. The review gives particular attention to the methods developed by the Human Connectome Project that may be especially relevant to studies of psychopathology; it outlines some of the key findings about what constitutes a brain region; and it highlights new information about the nature and stability of brain circuits. Some of the Human Connectome Project's new findings particularly relevant to psychopathology-about neural circuits and their relationships to behavior-are also presented. The review ends by discussing the extension of Human Connectome Project methods across the lifespan and into manifest illness. Potential treatment implications are also considered.

A Supervoxel-Based Method for Groupwise Whole Brain Parcellation with Resting-State fMRI Data

Node definition is a very important issue in human brain network analysis and functional connectivity studies. Typically, the atlases generated from meta-analysis, random criteria, and structural criteria are utilized as nodes in related applications. However, these atlases are not originally designed for such purposes and may not be suitable. In this study, we combined normalized cut (Ncut) and a supervoxel method called simple linear iterative clustering (SLIC) to parcellate whole brain resting-state fMRI data in order to generate appropriate brain atlases. Specifically, Ncut was employed to extract features from connectivity matrices, and then SLIC was applied on the extracted features to generate parcellations. To obtain group level parcellations, two approaches named mean SLIC and two-level SLIC were proposed. The cluster number varied in a wide range in order to generate parcellations with multiple granularities. The two SLIC approaches were compared with three state-of-the-art approaches under different evaluation metrics, which include spatial contiguity, functional homogeneity, and reproducibility. Both the group-to-group reproducibility and the group-to-subject reproducibility were evaluated in our study. The experimental results showed that the proposed approaches obtained relatively good overall clustering performances in different conditions that included different weighting functions, different sparsifying schemes, and several confounding factors. Therefore, the generated atlases are appropriate to be utilized as nodes for network analysis. The generated atlases and major source codes of this study have been made publicly available at http://www.nitrc.org/projects/slic/.

Cognitive complaints in older adults at risk for Alzheimer's disease are associated with altered resting-state networks

INTRODUCTION

Pathophysiological changes that accompany early clinical symptoms in prodromal Alzheimer's disease (AD) may have a disruptive influence on brain networks. We investigated resting-state functional magnetic resonance imaging (rsfMRI), combined with brain connectomics, to assess changes in whole-brain functional connectivity (FC) in relation to neurocognitive variables.

METHODS

Participants included 58 older adults who underwent rsfMRI. Individual FC matrices were computed based on a 278-region parcellation. FastICA decomposition was performed on a matrix combining all subjects' FC. Each FC pattern was then used as a response in a multilinear regression model including neurocognitive variables associated with AD (cognitive complaint index [CCI] scores from self and informant, an episodic memory score, and an executive function score).

RESULTS

Three connectivity independent component analysis (connICA) components (RSN, VIS, and FP-DMN FC patterns) associated with neurocognitive variables were identified based on prespecified criteria. RSN-pattern, characterized by increased FC within all resting-state networks, was negatively associated with self CCI. VIS-pattern, characterized by an increase in visual resting-state network, was negatively associated with CCI self or informant scores. FP-DMN-pattern, characterized by an increased interaction of frontoparietal and default mode networks (DMN), was positively associated with verbal episodic memory.

DISCUSSION

Specific patterns of FC were differently associated with neurocognitive variables thought to change early in the course of AD. An integrative connectomics approach relating cognition to changes in FC may help identify preclinical and early prodromal stages of AD and help elucidate the complex relationship between subjective and objective indices of cognitive decline and differences in brain functional organization.

Let's Not Waste Time: Using Temporal Information in Clustered Activity Estimation with Spatial Adjacency Restrictions (CAESAR) for Parcellating FMRI Data

We have proposed a Bayesian approach for functional parcellation of whole-brain FMRI measurements which we call Clustered Activity Estimation with Spatial Adjacency Restrictions (CAESAR). We use distance-dependent Chinese restaurant processes (dd-CRPs) to define a flexible prior which partitions the voxel measurements into clusters whose number and shapes are unknown a priori. With dd-CRPs we can conveniently implement spatial constraints to ensure that our parcellations remain spatially contiguous and thereby physiologically meaningful. In the present work, we extend CAESAR by using Gaussian process (GP) priors to model the temporally smooth haemodynamic signals that give rise to the measured FMRI data. A challenge for GP inference in our setting is the cubic scaling with respect to the number of time points, which can become computationally prohibitive with FMRI measurements, potentially consisting of long time series. As a solution we describe an efficient implementation that is practically as fast as the corresponding time-independent non-GP model with typically-sized FMRI data sets. We also employ a population Monte-Carlo algorithm that can significantly speed up convergence compared to traditional single-chain methods. First we illustrate the benefits of CAESAR and the GP priors with simulated experiments. Next, we demonstrate our approach by parcellating resting state FMRI data measured from twenty participants as taken from the Human Connectome Project data repository. Results show that CAESAR affords highly robust and scalable whole-brain clustering of FMRI timecourses.

Resting-State Magnetoencephalography Reveals Different Patterns of Aberrant Functional Connectivity in Combat-Related Mild Traumatic Brain Injury

Blast mild traumatic brain injury (mTBI) is a leading cause of sustained impairment in military service members and veterans. However, the mechanism of persistent disability is not fully understood. The present study investigated disturbances in brain functioning in mTBI participants using a source-imaging-based approach to analyze functional connectivity (FC) from resting-state magnetoencephalography (rs-MEG). Study participants included 26 active-duty service members or veterans who had blast mTBI with persistent post-concussive symptoms, and 22 healthy control active-duty service members or veterans. The source time courses from regions of interest (ROIs) were used to compute ROI to whole-brain (ROI-global) FC for different frequency bands using two different measures: 1) time-lagged cross-correlation and 2) phase-lock synchrony. Compared with the controls, blast mTBI participants showed increased ROI-global FC in beta, gamma, and low-frequency bands, but not in the alpha band. Sources of abnormally increased FC included the: 1) prefrontal cortex (right ventromedial prefrontal cortex [vmPFC], right rostral anterior cingulate cortex [rACC]), and left ventrolateral and dorsolateral prefrontal cortex; 2) medial temporal lobe (bilateral parahippocampus, hippocampus, and amygdala); and 3) right putamen and cerebellum. In contrast, the blast mTBI group also showed decreased FC of the right frontal pole. Group differences were highly consistent across the two different FC measures. FC of the left ventrolateral prefrontal cortex correlated with executive functioning and processing speed in mTBI participants. Altogether, our findings of increased and decreased regionalpatterns of FC suggest that disturbances in intrinsic brain connectivity may be the result of multiple mechanisms, and are associated with cognitive sequelae of the injury.

sGraSP: A graph-based method for the derivation of subject-specific functional parcellations of the brain

BACKGROUND

Resting-state fMRI (rs-fMRI) has emerged as a prominent tool for the study of functional connectivity. The identification of the regions associated with the different brain functions has received significant interest. However, most of the studies conducted so far have focused on the definition of a common set of regions, valid for an entire population. The variation of the functional regions within a population has rarely been accounted for.

NEW METHOD

In this paper, we propose sGraSP, a graph-based approach for the derivation of subject-specific functional parcellations. Our method generates first a common parcellation for an entire population, which is then adapted to each subject individually.

RESULTS

Several cortical parcellations were generated for 859 children being part of the Philadelphia Neurodevelopmental Cohort. The stability of the parcellations generated by sGraSP was tested by mixing population and subject rs-fMRI signals, to generate subject-specific parcels increasingly closer to the population parcellation. We also checked if the parcels generated by our method were better capturing a development trend underlying our data than the original parcels, defined for the entire population.

COMPARISON WITH EXISTING METHODS

We compared sGraSP with a simpler and faster approach based on a Voronoi tessellation, by measuring their ability to produce functionally coherent parcels adapted to the subject data.

CONCLUSIONS

Our parcellations outperformed the Voronoi tessellations. The parcels generated by sGraSP vary consistently with respect to signal mixing, the results are highly reproducible and the neurodevelopmental trend is better captured with the subject-specific parcellation, under all the signal mixing conditions.

Inferring Individual-Level Variations in the Functional Parcellation of the Cerebral Cortex

OBJECTIVE

Functional parcellation of the cerebral cortex is variable across different subjects or between cognitive states. Ignoring individual-or state-dependent variations in the functional parcellation may lead to inaccurate representations of individual functional connectivity, limiting the precision of interpretations of differences in individual connectivity profiles. However, it is difficult to infer the individual-level variations due to the relatively low robustness of methods for parcellation of individual subjects.

METHODS

We propose a method called "joint K-means" to robustly parcellate the cerebral cortex using functional magnetic resonance imaging (fMRI) data for contrasts between two states or subjects that intended to characterize variance in individual functional parcellations. The key idea of the proposed method is to jointly infer parcellations in contrasted datasets by iterative descent, while constraining the similarity of the two pathways in searches for local minima to reduce spurious variations.

RESULTS

Parcellations of resting-state fMRI datasets from the Human Connectome Project show that the similarity of parcellations for an individual subject studied on two sessions is greater than that between different subjects. Differences in parcellations between subjects are nonuniformly distributed across the cerebral cortex, with clusters of higher variance in the prefrontal, lateral temporal, and occipito-parietal cortices. This pattern is reproducible across sessions, between groups, and using different numbers of parcels.

CONCLUSION

The individual-level variations inferred by the proposed method are plausible and consistent with the previously reported functional connectivity variability.

SIGNIFICANCE

The proposed method is a promising tool for investigating relationships between the cerebral functional organization and behavioral differences.

Multisite reliability of MR-based functional connectivity

Recent years have witnessed an increasing number of multisite MRI functional connectivity (fcMRI) studies. While multisite studies provide an efficient way to accelerate data collection and increase sample sizes, especially for rare clinical populations, any effects of site or MRI scanner could ultimately limit power and weaken results. Little data exists on the stability of functional connectivity measurements across sites and sessions. In this study, we assess the influence of site and session on resting state functional connectivity measurements in a healthy cohort of traveling subjects (8 subjects scanned twice at each of 8 sites) scanned as part of the North American Prodrome Longitudinal Study (NAPLS). Reliability was investigated in three types of connectivity analyses: (1) seed-based connectivity with posterior cingulate cortex (PCC), right motor cortex (RMC), and left thalamus (LT) as seeds; (2) the intrinsic connectivity distribution (ICD), a voxel-wise connectivity measure; and (3) matrix connectivity, a whole-brain, atlas-based approach to assessing connectivity between nodes. Contributions to variability in connectivity due to subject, site, and day-of-scan were quantified and used to assess between-session (test-retest) reliability in accordance with Generalizability Theory. Overall, no major site, scanner manufacturer, or day-of-scan effects were found for the univariate connectivity analyses; instead, subject effects dominated relative to the other measured factors. However, summaries of voxel-wise connectivity were found to be sensitive to site and scanner manufacturer effects. For all connectivity measures, although subject variance was three times the site variance, the residual represented 60-80% of the variance, indicating that connectivity differed greatly from scan to scan independent of any of the measured factors (i.e., subject, site, and day-of-scan). Thus, for a single 5min scan, reliability across connectivity measures was poor (ICC=0.07-0.17), but increased with increasing scan duration (ICC=0.21-0.36 at 25min). The limited effects of site and scanner manufacturer support the use of multisite studies, such as NAPLS, as a viable means of collecting data on rare populations and increasing power in univariate functional connectivity studies. However, the results indicate that aggregation of fcMRI data across longer scan durations is necessary to increase the reliability of connectivity estimates at the single-subject level.

Individual variation in functional brain connectivity: implications for personalized approaches to psychiatric disease

Functional brain connectivity measured with functional magnetic resonance imaging (fMRI) is a popular technique for investigating neural organization in both healthy subjects and patients with mental illness. Despite a rapidly growing body of literature, however, functional connectivity research has yet to deliver biomarkers that can aid psychiatric diagnosis or prognosis at the single-subject level. One impediment to developing such practical tools has been uncertainty regarding the ratio of intra- to interindividual variability in functional connectivity; in other words, how much variance is state- versus trait-related. Here, we review recent evidence that functional connectivity profiles are both reliable within subjects and unique across subjects, and that features of these profiles relate to behavioral phenotypes. Together, these results suggest the potential to discover reliable correlates of present and future illness and/or response to treatment in the strength of an individual's functional brain connections. Ultimately, this work could help develop personalized approaches to psychiatric illness.

Fluctuations in Global Brain Activity Are Associated With Changes in Whole-Brain Connectivity of Functional Networks

OBJECTIVE

The aim of this study was to explore the relationship between global brain activity, changes in whole-brain connectivity, and changes in brain states across subjects using resting-state functional magnetic resonance imaging.

METHODS

We extended current methods that use a sparse set of coactivation patterns to extract critical time points in global brain activity. Critical activity time points were defined as points where the global signal is greater than one standard deviation above or below the average global signal. Four categories of critical points were defined along dimensions of global signal intensity and trajectory. Voxel-based methods were used to interrogate differences in connectivity between these critical points.

RESULTS

Several differences in connectivity were found in functional resting-state networks (RSNs) as a function of global activity. RSNs associated with cognitive functions in frontal, parietal, and subcortical regions exhibited greater whole-brain connectivity during lower global activity states. Meanwhile, RSNs associated with sensory functions exhibited greater whole-brain connectivity during the higher global activity states. Moreover, we present evidence that these results depend in part upon the standard deviation threshold used to define the critical points, suggesting critical points at different thresholds represent unique brain states.

CONCLUSION

Overall, the findings support the hypothesis that the brain oscillates through different states over the course of a resting-state study reflecting differences in RSN connectivity associated with global brain activity.

SIGNIFICANCE

Increased understanding of brain dynamics may help to elucidate individual differences in behavior and dysfunction.

A multi-modal parcellation of human cerebral cortex

Understanding the amazingly complex human cerebral cortex requires a map (or parcellation) of its major subdivisions, known as cortical areas. Making an accurate areal map has been a century-old objective in neuroscience. Using multi-modal magnetic resonance images from the Human Connectome Project (HCP) and an objective semi-automated neuroanatomical approach, we delineated 180 areas per hemisphere bounded by sharp changes in cortical architecture, function, connectivity, and/or topography in a precisely aligned group average of 210 healthy young adults. We characterized 97 new areas and 83 areas previously reported using post-mortem microscopy or other specialized study-specific approaches. To enable automated delineation and identification of these areas in new HCP subjects and in future studies, we trained a machine-learning classifier to recognize the multi-modal 'fingerprint' of each cortical area. This classifier detected the presence of 96.6% of the cortical areas in new subjects, replicated the group parcellation, and could correctly locate areas in individuals with atypical parcellations. The freely available parcellation and classifier will enable substantially improved neuroanatomical precision for studies of the structural and functional organization of human cerebral cortex and its variation across individuals and in development, aging, and disease.

Functional atlas of the awake rat brain: A neuroimaging study of rat brain specialization and integration

Connectivity-based parcellation approaches present an innovative method to segregate the brain into functionally specialized regions. These approaches have significantly advanced our understanding of the human brain organization. However, parallel progress in animal research is sparse. Using resting-state fMRI data and a novel, data-driven parcellation method, we have obtained robust functional parcellations of the rat brain. These functional parcellations reveal the regional specialization of the rat brain, which exhibited high within-parcel homogeneity and high reproducibility across animals. Graph analysis of the whole-brain network constructed based on these functional parcels indicates that the rat brain has a topological organization similar to humans, characterized by both segregation and integration. Our study also provides compelling evidence that the cingulate cortex is a functional hub region conserved from rodents to humans. Together, this study has characterized the rat brain specialization and integration, and has significantly advanced our understanding of the rat brain organization. In addition, it is valuable for studies of comparative functional neuroanatomy in mammalian brains.

Lateral-Medial Dissociation in Orbitofrontal Cortex-Hypothalamus Connectivity

The orbitofrontal cortex (OFC) is involved in cognitive functions, and is also closely related to autonomic functions. The OFC is densely connected with the hypothalamus, a heterogeneous structure controlling autonomic functions that can be divided into two major parts: the lateral and the medial. Resting-state functional connectivity has allowed us to parcellate the cerebral cortex into putative functional areas based on the changes in the spatial pattern of connectivity in the cerebral cortex when a seed point is moved from one voxel to another. In the present high spatial-resolution fMRI study, we investigate the connectivity-based organization of the OFC with reference to the hypothalamus. The OFC was parcellated using resting-state functional connectivity in an individual subject approach, and then the functional connectivity was examined between the parcellated areas in the OFC and the lateral/medial hypothalamus. We found a functional double dissociation in the OFC: the lateral OFC (the lateral orbital gyrus) was more likely connected with the lateral hypothalamus, whereas the medial OFC (the medial orbital and rectal gyri) was more likely connected with the medial hypothalamus. These results demonstrate the fundamental heterogeneity of the OFC, and suggest a potential neural basis of the OFC–hypothalamic functional interaction.

Group-wise parcellation of the cortex through multi-scale spectral clustering

The delineation of functionally and structurally distinct regions as well as their connectivity can provide key knowledge towards understanding the brain's behaviour and function. Cytoarchitecture has long been the gold standard for such parcellation tasks, but has poor scalability and cannot be mapped in vivo. Functional and diffusion magnetic resonance imaging allow in vivo mapping of brain's connectivity and the parcellation of the brain based on local connectivity information. Several methods have been developed for single subject connectivity driven parcellation, but very few have tackled the task of group-wise parcellation, which is essential for uncovering group specific behaviours. In this paper, we propose a group-wise connectivity-driven parcellation method based on spectral clustering that captures local connectivity information at multiple scales and directly enforces correspondences between subjects. The method is applied to diffusion Magnetic Resonance Imaging driven parcellation on two independent groups of 50 subjects from the Human Connectome Project. Promising quantitative and qualitative results in terms of information loss, modality comparisons, group consistency and inter-group similarities demonstrate the potential of the method.

Cerebellar Functional Parcellation Using Sparse Dictionary Learning Clustering

The human cerebellum has recently been discovered to contribute to cognition and emotion beyond the planning and execution of movement, suggesting its functional heterogeneity. We aimed to identify the functional parcellation of the cerebellum using information from resting-state functional magnetic resonance imaging (rs-fMRI). For this, we introduced a new data-driven decomposition-based functional parcellation algorithm, called Sparse Dictionary Learning Clustering (SDLC). SDLC integrates dictionary learning, sparse representation of rs-fMRI, and k-means clustering into one optimization problem. The dictionary is comprised of an over-complete set of time course signals, with which a sparse representation of rs-fMRI signals can be constructed. Cerebellar functional regions were then identified using k-means clustering based on the sparse representation of rs-fMRI signals. We solved SDLC using a multi-block hybrid proximal alternating method that guarantees strong convergence. We evaluated the reliability of SDLC and benchmarked its classification accuracy against other clustering techniques using simulated data. We then demonstrated that SDLC can identify biologically reasonable functional regions of the cerebellum as estimated by their cerebello-cortical functional connectivity. We further provided new insights into the cerebello-cortical functional organization in children.

Change in brain network topology as a function of treatment response in schizophrenia: a longitudinal resting-state fMRI study using graph theory

A number of neuroimaging studies have provided evidence in support of the hypothesis that faulty interactions between spatially disparate brain regions underlie the pathophysiology of schizophrenia, but it remains unclear to what degree antipsychotic medications affect these. We hypothesized that the balance between functional integration and segregation of brain networks is impaired in unmedicated patients with schizophrenia, but that it can be partially restored by antipsychotic medications. We included 32 unmedicated patients with schizophrenia (SZ) and 32 matched healthy controls (HC) in this study. We obtained resting-state scans while unmedicated, and again after 6 weeks of treatment with risperidone to assess functional integration and functional segregation of brain networks using graph theoretical measures. Compared with HC, unmedicated SZ showed reduced global efficiency and increased clustering coefficients. This pattern of aberrant functional network integration and segregation was modulated with antipsychotic medications, but only in those who responded to treatment. Our work lends support to the concept of schizophrenia as a dysconnectivity syndrome, and suggests that faulty brain network topology in schizophrenia is modulated by antipsychotic medication as a function of treatment response.

Group-wise consistent cortical parcellation based on connectional profiles

For decades, seeking common, consistent and corresponding anatomical/functional regions across individual brains via cortical parcellation has been a longstanding challenging problem. In our opinion, two major barriers to solve this problem are determining meaningful cortical boundaries that segregate homogeneous regions and establishing correspondences among parcellated regions of multiple brains. To establish a corresponding system across subjects, we recently developed the Dense Individualized and Common Connectivity-based Cortical Landmarks (DICCCOL) system which possesses group-wise consistent white matter fiber connection patterns across individuals and thus provides a dense map of corresponding cortical landmarks. Despite this useful property, however, the DICCCOL landmarks are still far from covering the whole cerebral cortex and do not provide clear structural/functional cortical boundaries. To address the above limitation while leveraging the advantage of DICCCOL, in this paper, we present a novel approach for group-wise consistent parcellation of the cerebral cortex via a hierarchical scheme. In each hierarchical level, DICCCOLs are used as corresponding samples to automatically determine the cluster number so that other cortical surface vertices are iteratively classified into corresponding clusters across subjects within a group-wise classification framework. Experimental results showed that this approach can achieve consistent fine-granularity cortical parcellation with intrinsically-established structural correspondences across individual brains. Besides, comparisons with resting-state and task-based fMRI datasets demonstrated that the group-wise parcellation boundaries segregate functionally homogeneous areas.

Spin-Echo Resting-State Functional Connectivity in High-Susceptibility Regions: Accuracy, Reliability, and the Impact of Physiological Noise

Gradient-echo (GE) echo-planar imaging (EPI) is the method of choice in blood-oxygenation level-dependent (BOLD) functional MRI (fMRI) studies, as it demonstrates substantially higher BOLD sensitivity than its spin-echo (SE) counterpart. However, it is also well known that the GE-EPI signal is prone to signal dropouts and shifts due to susceptibility effects near air-tissue interfaces. SE-EPI, in contrast, is minimally affected by these artifacts. In this study, we quantify, for the first time, the sensitivity and specificity of SE and GE EPI for resting-state fMRI functional connectivity (fcMRI) mapping, using the 1000-brain fcMRI atlas (Yeo et al., 2011 ) as the pseudoground truth. Moreover, we assess the influence of physiological processes on resting-state BOLD measured using both regular and ultrafast GE and SE acquisitions. Our work demonstrates that SE-EPI and GE-EPI are associated with similar sensitivities, specificities, and intersubject reproducibility in fcMRI for most brain networks, generated using both seed-based analysis and independent component analysis. More importantly, SE-based fcMRI measurements demonstrated significantly higher sensitivity, specificity, and intersubject reproducibility in high-susceptibility regions, spanning the limbic and frontal networks in the 1000-brain atlas. In addition, SE-EPI is significantly less sensitive to prominent sources of physiological noise, including low-frequency respiratory volume and heart rate variations. Our work suggests that SE-EPI should be increasingly adopted in the study of networks spanning susceptibility-affected brain regions, including those that are important to memory, language, and emotion.

Groupwise connectivity-based parcellation of the whole human cortical surface using watershed-driven dimension reduction

Segregating the human cortex into distinct areas based on structural connectivity criteria is of widespread interest in neuroscience. This paper presents a groupwise connectivity-based parcellation framework for the whole cortical surface using a new high quality diffusion dataset of 79 healthy subjects. Our approach performs gyrus by gyrus to parcellate the whole human cortex. The main originality of the method is to compress for each gyrus the connectivity profiles used for the clustering without any anatomical prior information. This step takes into account the interindividual cortical and connectivity variability. To this end, we consider intersubject high density connectivity areas extracted using a surface-based watershed algorithm. A wide validation study has led to a fully automatic pipeline which is robust to variations in data preprocessing (tracking type, cortical mesh characteristics and boundaries of initial gyri), data characteristics (including number of subjects), and the main algorithmic parameters. A remarkable reproducibility is achieved in parcellation results for the whole cortex, leading to clear and stable cortical patterns. This reproducibility has been tested across non-overlapping subgroups and the validation is presented mainly on the pre- and postcentral gyri.

The modular and integrative functional architecture of the human brain

Network-based analyses of brain imaging data consistently reveal distinct modules and connector nodes with diverse global connectivity across the modules. How discrete the functions of modules are, how dependent the computational load of each module is to the other modules' processing, and what the precise role of connector nodes is for between-module communication remains underspecified. Here, we use a network model of the brain derived from resting-state functional MRI (rs-fMRI) data and investigate the modular functional architecture of the human brain by analyzing activity at different types of nodes in the network across 9,208 experiments of 77 cognitive tasks in the BrainMap database. Using an author-topic model of cognitive functions, we find a strong spatial correspondence between the cognitive functions and the network's modules, suggesting that each module performs a discrete cognitive function. Crucially, activity at local nodes within the modules does not increase in tasks that require more cognitive functions, demonstrating the autonomy of modules' functions. However, connector nodes do exhibit increased activity when more cognitive functions are engaged in a task. Moreover, connector nodes are located where brain activity is associated with many different cognitive functions. Connector nodes potentially play a role in between-module communication that maintains the modular function of the brain. Together, these findings provide a network account of the brain's modular yet integrated implementation of cognitive functions.

Parcellating cortical functional networks in individuals

The capacity to identify the unique functional architecture of an individual's brain is a crucial step toward personalized medicine and understanding the neural basis of variation in human cognition and behavior. Here we developed a cortical parcellation approach to accurately map functional organization at the individual level using resting-state functional magnetic resonance imaging (fMRI). A population-based functional atlas and a map of inter-individual variability were employed to guide the iterative search for functional networks in individual subjects. Functional networks mapped by this approach were highly reproducible within subjects and effectively captured the variability across subjects, including individual differences in brain lateralization. The algorithm performed well across different subject populations and data types, including task fMRI data. The approach was then validated by invasive cortical stimulation mapping in surgical patients, suggesting potential for use in clinical applications.

A neuromarker of sustained attention from whole-brain functional connectivity

Although attention plays a ubiquitous role in perception and cognition, researchers lack a simple way to measure a person's overall attentional abilities. Because behavioral measures are diverse and difficult to standardize, we pursued a neuromarker of an important aspect of attention, sustained attention, using functional magnetic resonance imaging. To this end, we identified functional brain networks whose strength during a sustained attention task predicted individual differences in performance. Models based on these networks generalized to previously unseen individuals, even predicting performance from resting-state connectivity alone. Furthermore, these same models predicted a clinical measure of attention--symptoms of attention deficit hyperactivity disorder--from resting-state connectivity in an independent sample of children and adolescents. These results demonstrate that whole-brain functional network strength provides a broadly applicable neuromarker of sustained attention.

Functional connectome fingerprinting: identifying individuals using patterns of brain connectivity

Functional magnetic resonance imaging (fMRI) studies typically collapse data from many subjects, but brain functional organization varies between individuals. Here we establish that this individual variability is both robust and reliable, using data from the Human Connectome Project to demonstrate that functional connectivity profiles act as a 'fingerprint' that can accurately identify subjects from a large group. Identification was successful across scan sessions and even between task and rest conditions, indicating that an individual's connectivity profile is intrinsic, and can be used to distinguish that individual regardless of how the brain is engaged during imaging. Characteristic connectivity patterns were distributed throughout the brain, but the frontoparietal network emerged as most distinctive. Furthermore, we show that connectivity profiles predict levels of fluid intelligence: the same networks that were most discriminating of individuals were also most predictive of cognitive behavior. Results indicate the potential to draw inferences about single subjects on the basis of functional connectivity fMRI.

A human brain atlas derived via n-cut parcellation of resting-state and task-based fMRI data

The growth of functional MRI has led to development of human brain atlases derived by parcellating resting-state connectivity patterns into functionally independent regions of interest (ROIs). All functional atlases to date have been derived from resting-state fMRI data. But given that functional connectivity between regions varies with task, we hypothesized that an atlas incorporating both resting-state and task-based fMRI data would produce an atlas with finer characterization of task-relevant regions than an atlas derived from resting-state alone. To test this hypothesis, we derived parcellation atlases from twenty-nine healthy adult participants enrolled in the Cognitive Connectome project, an initiative to improve functional MRI's translation into clinical decision-making by mapping normative variance in brain-behavior relationships. Participants underwent resting-state and task-based fMRI spanning nine cognitive domains: motor, visuospatial, attention, language, memory, affective processing, decision-making, working memory, and executive function. Spatially constrained n-cut parcellation derived brain atlases using (1) all participants' functional data (Task) or (2) a single resting-state scan (Rest). An atlas was also derived from random parcellation for comparison purposes (Random). Two methods were compared: (1) a parcellation applied to the group's mean edge weights (mean), and (2) a two-stage approach with parcellation of individual edge weights followed by parcellation of mean binarized edges (two-stage). The resulting Task and Rest atlases had significantly greater similarity with each other (mean Jaccard indices JI=0.72-0.85) than with the Random atlases (JI=0.59-0.63; all p<0.001 after Bonferroni correction). Task and Rest atlas similarity was greatest for the two-stage method (JI=0.85), which has been shown as more robust than the mean method; these atlases also better reproduced voxelwise seed maps of the left dorsolateral prefrontal cortex during rest and performing the n-back working memory task (r=0.75-0.80) than the Random atlases (r=0.64-0.72), further validating their utility. We expected regions governing higher-order cognition (such as frontal and anterior temporal lobes) to show greatest difference between Task and Rest atlases; contrary to expectations, these areas had greatest similarity between atlases. Our findings indicate that atlases derived from parcellation of task-based and resting-state fMRI data are highly comparable, and existing resting-state atlases are suitable for task-based analyses. We introduce an anatomically labeled fMRI-derived whole-brain human atlas for future Cognitive Connectome analyses.

Multimodal description of whole brain connectivity: A comparison of resting state MEG, fMRI, and DWI

Structural and functional connectivity (SC and FC) have received much attention over the last decade, as they offer unique insight into the coordination of brain functioning. They are often assessed independently with three imaging modalities: SC using diffusion‐weighted imaging (DWI), FC using functional magnetic resonance imaging (fMRI), and magnetoencephalography/electroencephalography (MEG/EEG). DWI provides information about white matter organization, allowing the reconstruction of fiber bundles. fMRI uses blood‐oxygenation level‐dependent (BOLD) contrast to indirectly map neuronal activation. MEG and EEG are direct measures of neuronal activity, as they are sensitive to the synchronous inputs in pyramidal neurons. Seminal studies have targeted either the electrophysiological substrate of BOLD or the anatomical basis of FC. However, multimodal comparisons have been scarcely performed, and the relation between SC, fMRI‐FC, and MEG‐FC is still unclear. Here we present a systematic comparison of SC, resting state fMRI‐FC, and MEG‐FC between cortical regions, by evaluating their similarities at three different scales: global network, node, and hub distribution. We obtained strong similarities between the three modalities, especially for the following pairwise combinations: SC and fMRI‐FC; SC and MEG‐FC at theta, alpha, beta and gamma bands; and fMRI‐FC and MEG‐FC in alpha and beta. Furthermore, highest node similarity was found for regions of the default mode network and primary motor cortex, which also presented the highest hubness score. Distance was partially responsible for these similarities since it biased all three connectivity estimates, but not the unique contributor, since similarities remained after controlling for distance. Hum Brain Mapp 37:20–34, 2016. © 2015 The Authors Human Brain Mapping Published by Wiley Periodicals, Inc.

Nodal centrality of functional network in the differentiation of schizophrenia

A disturbance in the integration of information during mental processing has been implicated in schizophrenia, possibly due to faulty communication within and between brain regions. Graph theoretic measures allow quantification of functional brain networks. Functional networks are derived from correlations between time courses of brain regions. Group differences between SZ and control groups have been reported for functional network properties, but the potential of such measures to classify individual cases has been little explored. We tested whether the network measure of betweenness centrality could classify persons with schizophrenia and normal controls. Functional networks were constructed for 19 schizophrenic patients and 29 non-psychiatric controls based on resting state functional MRI scans. The betweenness centrality of each node, or fraction of shortest-paths that pass through it, was calculated in order to characterize the centrality of the different regions. The nodes with high betweenness centrality agreed well with hub nodes reported in previous studies of structural and functional networks. Using a linear support vector machine algorithm, the schizophrenia group was differentiated from non-psychiatric controls using the ten nodes with the highest betweenness centrality. The classification accuracy was around 80%, and stable against connectivity thresholding. Better performance was achieved when using the ranks as feature space as opposed to the actual values of betweenness centrality. Overall, our findings suggest that changes in functional hubs are associated with schizophrenia, reflecting a variation of the underlying functional network and neuronal communications. In addition, a specific network property, betweenness centrality, can classify persons with SZ with a high level of accuracy.

Connectivity-based parcellation: Critique and implications

Regional specialization and functional integration are often viewed as two fundamental principles of human brain organization. They are closely intertwined because each functionally specialized brain region is probably characterized by a distinct set of long-range connections. This notion has prompted the quickly developing family of connectivity-based parcellation (CBP) methods in neuroimaging research. CBP assumes that there is a latent structure of parcels in a region of interest (ROI). First, connectivity strengths are computed to other parts of the brain for each voxel/vertex within the ROI. These features are then used to identify functionally distinct groups of ROI voxels/vertices. CBP enjoys increasing popularity for the in-vivo mapping of regional specialization in the human brain. Due to the requirements of different applications and datasets, CBP has diverged into a heterogeneous family of methods. This broad overview critically discusses the current state as well as the commonalities and idiosyncrasies of the main CBP methods. We target frequent concerns faced by novices and veterans to provide a reference for the investigation and review of CBP studies.

Preterm birth alters neonatal, functional rich club organization

Alterations in neural networks are associated with the cognitive difficulties of the prematurely born. Using functional magnetic resonance imaging, we analyzed functional connectivity for preterm (PT) and term neonates at term equivalent age. Specifically, we constructed whole-brain networks and examined rich club (RC) organization, a common construct among complex systems where important (or "rich") nodes connect preferentially to other important nodes. Both PT and term neonates showed RC organization with PT neonates exhibiting significantly reduced connections between these RC nodes. Additionally, PT neonates showed evidence of weaker functional segregation. Our results suggest that PT birth is associated with fundamental changes of functional organization in the developing brain.

Integrated Analysis and Visualization of Group Differences in Structural and Functional Brain Connectivity: Applications in Typical Ageing and Schizophrenia

Structural and functional brain connectivity are increasingly used to identify and analyze group differences in studies of brain disease. This study presents methods to analyze uni- and bi-modal brain connectivity and evaluate their ability to identify differences. Novel visualizations of significantly different connections comparing multiple metrics are presented. On the global level, “bi-modal comparison plots” show the distribution of uni- and bi-modal group differences and the relationship between structure and function. Differences between brain lobes are visualized using “worm plots”. Group differences in connections are examined with an existing visualization, the “connectogram”. These visualizations were evaluated in two proof-of-concept studies: (1) middle-aged versus elderly subjects; and (2) patients with schizophrenia versus controls. Each included two measures derived from diffusion weighted images and two from functional magnetic resonance images. The structural measures were minimum cost path between two anatomical regions according to the “Statistical Analysis of Minimum cost path based Structural Connectivity” method and the average fractional anisotropy along the fiber. The functional measures were Pearson’s correlation and partial correlation of mean regional time series. The relationship between structure and function was similar in both studies. Uni-modal group differences varied greatly between connectivity types. Group differences were identified in both studies globally, within brain lobes and between regions. In the aging study, minimum cost path was highly effective in identifying group differences on all levels; fractional anisotropy and mean correlation showed smaller differences on the brain lobe and regional levels. In the schizophrenia study, minimum cost path and fractional anisotropy showed differences on the global level and within brain lobes; mean correlation showed small differences on the lobe level. Only fractional anisotropy and mean correlation showed regional differences. The presented visualizations were helpful in comparing and evaluating connectivity measures on multiple levels in both studies.

Influence of ROI selection on resting state functional connectivity: an individualized approach for resting state fMRI analysis

The differences in how our brain is connected are often thought to reflect the differences in our individual personalities and cognitive abilities. Individual differences in brain connectivity has long been recognized in the neuroscience community however it has yet to manifest itself in the methodology of resting state analysis. This is evident as previous studies use the same region of interest (ROIs) for all subjects. In this paper we demonstrate that the use of ROIs which are standardized across individuals leads to inaccurate calculations of functional connectivity. We also show that this problem can be addressed by taking an individualized approach by using subject-specific ROIs. Finally we show that ROI selection can affect the way we interpret our data by showing different changes in functional connectivity with aging.

The (in)stability of functional brain network measures across thresholds

The large-scale organization of the brain has features of complex networks that can be quantified using network measures from graph theory. However, many network measures were designed to be calculated on binary graphs, whereas functional brain organization is typically inferred from a continuous measure of correlations in temporal signal between brain regions. Thresholding is a necessary step to use binary graphs derived from functional connectivity data. However, there is no current consensus on what threshold to use, and network measures and group contrasts may be unstable across thresholds. Nevertheless, whole-brain network analyses are being applied widely with findings typically reported at an arbitrary threshold or range of thresholds. This study sought to evaluate the stability of network measures across thresholds in a large resting state functional connectivity dataset. Network measures were evaluated across absolute (correlation-based) and proportional (sparsity-based) thresholds, and compared between sex and age groups. Overall, network measures were found to be unstable across absolute thresholds. For example, the direction of group differences in a given network measure may change depending on the threshold. Network measures were found to be more stable across proportional thresholds. These results demonstrate that caution should be used when applying thresholds to functional connectivity data and when interpreting results from binary graph models.

Increased segregation of brain networks in focal epilepsy: An fMRI graph theory finding

Focal epilepsy is conceived of as activating local areas of the brain as well as engaging regional brain networks. Graph theory represents a powerful quantitative framework for investigation of brain networks. Here we investigate whether functional network changes are present in extratemporal focal epilepsy.

Task-free functional magnetic resonance imaging data from 15 subjects with extratemporal epilepsy and 26 age and gender matched healthy controls were used for analysis. Local network properties were calculated using local efficiency, clustering coefficient and modularity metrics. Global network properties were assessed with global efficiency and betweenness centrality metrics. Cost-efficiency of the networks at both local and global levels was evaluated by estimating the physical distance between functionally connected nodes, in addition to the overall numbers of connections in the network.

Clustering coefficient, local efficiency and modularity were significantly higher in individuals with focal epilepsy than healthy control subjects, while global efficiency and betweenness centrality were not significantly different between the two groups. Local network properties were also highly efficient, at low cost, in focal epilepsy subjects compared to healthy controls.

Our results show that functional networks in focal epilepsy are altered in a way that the nodes of the network are more isolated. We postulate that network regularity, or segregation of the nodes of the networks, may be an adaptation that inhibits the conversion of the interictal state to seizures. It remains possible that this may be part of the epileptogenic process or an effect of medications.

•

Focal epilepsy is a neurological disorder known to affect neural networks.

•

We examine fMRI network changes during rest in extratemporal epilepsy patients.

•

Extratemporal epilepsy subjects show increased (isolated) local network regularity.

•

Network isolation may represent a ‘seizure-protection’ mechanism in focal epilepsy.

Direct imaging of functional networks

In blood-oxygenation-level-dependent functional magnetic resonance imaging (fMRI), current methods typically acquire ∼ 500,000 imaging voxels at each time point, and then use computer algorithms to reduce this data to the coefficients of a few hundred parcels or networks. This suggests that the amount of relevant information present in the fMRI signal is relatively small, and presents an opportunity to greatly improve the speed and signal to noise ratio (SNR) of the fMRI process. In this work, a theoretical framework is presented for calculating the coefficients of functional networks directly from highly undersampled fMRI data. Using predefined functional parcellations or networks and a compact k-space trajectory that samples data at optimal spatial scales, the problem of estimating network coefficients is reformulated to allow for direct least squares estimation, without Fourier encoding. By simulation, this approach is shown to allow for acceleration of the imaging process under ideal circumstances by nearly three orders of magnitude.

Development and validation of consensus clustering-based framework for brain segmentation using resting fMRI

BACKGROUND

Clustering methods are increasingly employed to segment brain regions into functional subdivisions using resting-state functional magnetic resonance imaging (rs-fMRI). However, these methods are highly sensitive to the (i) precise algorithms employed, (ii) their initializations, and (iii) metrics used for uncovering the optimal number of clusters from the data.

NEW METHOD

To address these issues, we develop a novel consensus clustering evidence accumulation (CC-EAC) framework, which effectively combines multiple clustering methods for segmenting brain regions using rs-fMRI data. Using extensive computer simulations, we examine the performance of widely used clustering algorithms including K-means, hierarchical, and spectral clustering as well as their combinations. We also examine the accuracy and validity of five objective criteria for determining the optimal number of clusters: mutual information, variation of information, modified silhouette, Rand index, and probabilistic Rand index.

RESULTS

A CC-EAC framework with a combination of base K-means clustering (KC) and hierarchical clustering (HC) with probabilistic Rand index as the criterion for choosing the optimal number of clusters, accurately uncovered the correct number of clusters from simulated datasets. In experimental rs-fMRI data, these methods reliably detected functional subdivisions of the supplementary motor area, insula, intraparietal sulcus, angular gyrus, and striatum.

COMPARISON WITH EXISTING METHODS

Unlike conventional approaches, CC-EAC can accurately determine the optimal number of stable clusters in rs-fMRI data, and is robust to initialization and choice of free parameters.

CONCLUSIONS

A novel CC-EAC framework is proposed for segmenting brain regions, by effectively combining multiple clustering methods and identifying optimal stable functional clusters in rs-fMRI data.

Joint Spectral Decomposition for the Parcellation of the Human Cerebral Cortex Using Resting-State fMRI

Identification of functional connections within the human brain has gained a lot of attention due to its potential to reveal neural mechanisms. In a whole-brain connectivity analysis, a critical stage is the computation of a set of network nodes that can effectively represent cortical regions. To address this problem, we present a robust cerebral cortex parcellation method based on spectral graph theory and resting-state fMRI correlations that generates reliable parcellations at the single-subject level and across multiple subjects. Our method models the cortical surface in each hemisphere as a mesh graph represented in the spectral domain with its eigenvectors. We connect cortices of different subjects with each other based on the similarity of their connectivity profiles and construct a multi-layer graph, which effectively captures the fundamental properties of the whole group as well as preserves individual subject characteristics. Spectral decomposition of this joint graph is used to cluster each cortical vertex into a subregion in order to obtain whole-brain parcellations. Using rs-fMRI data collected from 40 healthy subjects, we show that our proposed algorithm computes highly reproducible parcellations across different groups of subjects and at varying levels of detail with an average Dice score of 0.78, achieving up to 9% better reproducibility compared to existing approaches. We also report that our group-wise parcellations are functionally more consistent, thus, can be reliably used to represent the population in network analyses.

GraSP: geodesic Graph-based Segmentation with Shape Priors for the functional parcellation of the cortex

Resting-state functional MRI is a powerful technique for mapping the functional organization of the human brain. However, for many types of connectivity analysis, high-resolution voxelwise analyses are computationally infeasible and dimensionality reduction is typically used to limit the number of network nodes. Most commonly, network nodes are defined using standard anatomic atlases that do not align well with functional neuroanatomy or regions of interest covering a small portion of the cortex. Data-driven parcellation methods seek to overcome such limitations, but existing approaches are highly dependent on initialization procedures and produce spatially fragmented parcels or overly isotropic parcels that are unlikely to be biologically grounded. In this paper, we propose a novel graph-based parcellation method that relies on a discrete Markov Random Field framework. The spatial connectedness of the parcels is explicitly enforced by shape priors. The shape of the parcels is adapted to underlying data through the use of functional geodesic distances. Our method is initialization-free and rapidly segments the cortex in a single optimization. The performance of the method was assessed using a large developmental cohort of more than 850 subjects. Compared to two prevalent parcellation methods, our approach provides superior reproducibility for a similar data fit. Furthermore, compared to other methods, it avoids incoherent parcels. Finally, the method's utility is demonstrated through its ability to detect strong brain developmental effects that are only weakly observed using other methods.

Generation and Evaluation of a Cortical Area Parcellation from Resting-State Correlations

The cortical surface is organized into a large number of cortical areas; however, these areas have not been comprehensively mapped in the human. Abrupt transitions in resting-state functional connectivity (RSFC) patterns can noninvasively identify locations of putative borders between cortical areas (RSFC-boundary mapping; Cohen et al. 2008). Here we describe a technique for using RSFC-boundary maps to define parcels that represent putative cortical areas. These parcels had highly homogenous RSFC patterns, indicating that they contained one unique RSFC signal; furthermore, the parcels were much more homogenous than a null model matched for parcel size when tested in two separate datasets. Several alternative parcellation schemes were tested this way, and no other parcellation was as homogenous as or had as large a difference compared with its null model. The boundary map-derived parcellation contained parcels that overlapped with architectonic mapping of areas 17, 2, 3, and 4. These parcels had a network structure similar to the known network structure of the brain, and their connectivity patterns were reliable across individual subjects. These observations suggest that RSFC-boundary map-derived parcels provide information about the location and extent of human cortical areas. A parcellation generated using this method is available at http://www.nil.wustl.edu/labs/petersen/Resources.html.