Optimizing affinity measures for parcellating brain structures based on resting state fMRI data: a validation on medial superior frontal cortex

Concurrent brain parcellation and connectivity estimation via co-clustering of resting state fMRI data: A novel approach

Connectional topography mapping has been gaining widespread attention in human brain imaging studies. However, existing methods might not effectively utilize the information from neuroimaging data, thus hindering the understanding of the underlying connectional organization in the brain and uncovering the optimal clustering number from the data. In this study, we propose a novel method for the automated construction of inherent functional connectivity topography in a data‐driven manner by leveraging the power of co‐clustering‐based on resting state fMRI (rs‐fMRI) data. We propose the co‐clustering‐based method not only for concurrently parcellating two interconnected brain regions of interest (ROIs) under consideration into functionally homogenous subregions, but also for estimating the connectivity between these subregions from the two brain ROIs. In particular, we first model the connectional topography mapping as a co‐clustering‐based bipartite graph partitioning problem for constructing the inherent functional connectivity topography between the two interconnected brain ROIs. We also adopt an objective criterion, that is, silhouette width index measuring clustering quality, for determining the optimal number of clusters. The proposed method has been validated for mapping thalamocortical connectional topography based on rs‐fMRI data of 57 subjects. Validation results have demonstrated that our method identified the optimal solution with five pairs of mutually connected subregions of the thalamocortical system from the rs‐fMRI data, and could yield more meaningful, interpretable, and homogenous connectional topography than existing methods. The proposed method was further validated by the high symmetry of the mapped connectional topography between two hemispheres.

Functional parcellation of the hippocampus by semi-supervised clustering of resting state fMRI data

Many unsupervised methods are widely used for parcellating the brain. However, unsupervised methods aren’t able to integrate prior information, obtained from such as exiting functional neuroanatomy studies, to parcellate the brain, whereas the prior information guided semi-supervised method can generate more reliable brain parcellation. In this study, we propose a novel semi-supervised clustering method for parcellating the brain into spatially and functionally consistent parcels based on resting state functional magnetic resonance imaging (fMRI) data. Particularly, the prior supervised and spatial information is integrated into spectral clustering to achieve reliable brain parcellation. The proposed method has been validated in the hippocampus parcellation based on resting state fMRI data of 20 healthy adult subjects. The experimental results have demonstrated that the proposed method could successfully parcellate the hippocampus into head, body and tail parcels. The distinctive functional connectivity patterns of these parcels have further demonstrated the validity of the parcellation results. The effects of aging on the three hippocampus parcels’ functional connectivity were also explored across the healthy adult subjects. Compared with state-of-the-art methods, the proposed method had better performance on functional homogeneity. Furthermore, the proposed method had good test–retest reproducibility validated by parcellating the hippocampus based on three repeated resting state fMRI scans from 24 healthy adult subjects.

Regularized-Ncut: Robust and homogeneous functional parcellation of neonate and adult brain networks

Brain network parcellation based on resting-state functional MRI (rs-fMRI) is affected by noise, resulting in spurious small patches and decreased functional homogeneity within each network. Obtaining robust and homogeneous parcellation of neonate brain is more difficult, because neonate rs-fMRI is associated with relatively higher level of noise and no prior knowledge from a functional neonate atlas is available as spatial constraints. To meet these challenges, we developed a novel data-driven Regularized Normalized-cut (RNcut) method. RNcut is formulated by adding two regularization terms, a smoothing term using Markov random fields and a small-patch removal term, to conventional normalized-cut (Ncut) method. The RNcut and competing methods were tested with simulated datasets with known ground truth and then applied to both adult and neonate rs-fMRI datasets. Based on the parcellated networks generated by RNcut, intra-network connectivity was quantified. The test results from simulated datasets demonstrated that the RNcut method is more robust (p < 0.01) to noise and can delineate parcellated functional networks with significantly better (p < 0.01) spatial contiguity and significantly higher (p < 0.01) functional homogeneity than competing methods. Application of RNcut to neonate and adult rs-fMRI dataset revealed distinctive functional brain organization of neonate brains from that of adult brains. Collectively, we developed a novel data-driven RNcut method by integrating conventional Ncut with two regularization terms, generating robust and homogeneous functional parcellation without imposing spatial constraints. A broad range of brain network applications and analyses, especially neonate and infant brain parcellation with noisy and large sample of datasets, can potentially benefit from this RNcut method.

Generation of Individual Whole-Brain Atlases With Resting-State fMRI Data Using Simultaneous Graph Computation and Parcellation

The human brain can be characterized as functional networks. Therefore, it is important to subdivide the brain appropriately in order to construct reliable networks. Resting-state functional connectivity-based parcellation is a commonly used technique to fulfill this goal. Here we propose a novel individual subject-level parcellation approach based on whole-brain resting-state functional magnetic resonance imaging (fMRI) data. We first used a supervoxel method known as simple linear iterative clustering directly on resting-state fMRI time series to generate supervoxels, and then combined similar supervoxels to generate clusters using a clustering method known as graph-without-cut (GWC). The GWC approach incorporates spatial information and multiple features of the supervoxels by energy minimization, simultaneously yielding an optimal graph and brain parcellation. Meanwhile, it theoretically guarantees that the actual cluster number is exactly equal to the initialized cluster number. By comparing the results of the GWC approach and those of the random GWC approach, we demonstrated that GWC does not rely heavily on spatial structures, thus avoiding the challenges encountered in some previous whole-brain parcellation approaches. In addition, by comparing the GWC approach to two competing approaches, we showed that GWC achieved better parcellation performances in terms of different evaluation metrics. The proposed approach can be used to generate individualized brain atlases for applications related to cognition, development, aging, disease, personalized medicine, etc. The major source codes of this study have been made publicly available at https://github.com/yuzhounh/GWC.

Individualized Functional Parcellation of the Human Amygdala Using a Semi-supervised Clustering Method: A 7T Resting State fMRI Study

The amygdala plays an important role in emotional functions and its dysfunction is considered to be associated with multiple psychiatric disorders in humans. Cytoarchitectonic mapping has demonstrated that the human amygdala complex comprises several subregions. However, it's difficult to delineate boundaries of these subregions in vivo even if using state of the art high resolution structural MRI. Previous attempts to parcellate this small structure using unsupervised clustering methods based on resting state fMRI data suffered from the low spatial resolution of typical fMRI data, and it remains challenging for the unsupervised methods to define subregions of the amygdala in vivo. In this study, we developed a novel brain parcellation method to segment the human amygdala into spatially contiguous subregions based on 7T high resolution fMRI data. The parcellation was implemented using a semi-supervised spectral clustering (SSC) algorithm at an individual subject level. Under guidance of prior information derived from the Julich cytoarchitectonic atlas, our method clustered voxels of the amygdala into subregions according to similarity measures of their functional signals. As a result, three distinct amygdala subregions can be obtained in each hemisphere for every individual subject. Compared with the cytoarchitectonic atlas, our method achieved better performance in terms of subregional functional homogeneity. Validation experiments have also demonstrated that the amygdala subregions obtained by our method have distinctive, lateralized functional connectivity (FC) patterns. Our study has demonstrated that the semi-supervised brain parcellation method is a powerful tool for exploring amygdala subregional functions.

Topographical Information-Based High-Order Functional Connectivity and Its Application in Abnormality Detection for Mild Cognitive Impairment

Temporal synchronization-based functional connectivity (FC) has long been used by the neuroscience community. However, topographical FC information may provide additional information to characterize the advanced relationship between two brain regions. Accordingly, we proposed a novel method, namely high-order functional connectivity (HOFC), to capture this second-level relationship using inter-regional resemblance of the FC topographical profiles. Specifically, HOFC first calculates an FC profile for each brain region, notably between the given brain region and other brain regions. Based on these FC profiles, a second layer of correlations is computed between all pairs of brain regions (i.e., correlation's correlation). On this basis, we generated an HOFC network, where "high-order" network properties were computed. We found that HOFC was discordant with the traditional FC in several links, indicating additional information being revealed by the new metrics. We applied HOFC to identify biomarkers for early detection of Alzheimer's disease by comparing 77 mild cognitive impairment patients with 89 healthy individuals (control group). Sensitivity in detection of group difference was consistently improved by ∼25% using HOFC compared to using FC. An HOFC network analysis also provided complementary information to an FC network analysis. For example, HOFC between olfactory and orbitofrontal cortices was found significantly reduced in patients, besides extensive alterations in HOFC network properties. In conclusion, our results showed promise in using HOFC to comprehensively map the human brain connectome.

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/.