Robust brain parcellation using sparse representation on resting-state fMRI

Sensory and motor cortices parcellations estimated via distance-weighted sparse representation with application to autism spectrum disorder

BACKGROUND

Motor impairments and sensory processing abnormalities are prevalent in autism spectrum disorder (ASD), closely related to the core functions of the primary motor cortex (M1) and the primary somatosensory cortex (S1). Currently, there is limited knowledge about potential therapeutic targets in the subregions of M1 and S1 in ASD patients. This study aims to map clinically significant functional subregions of M1 and S1.

METHODS

Resting-state functional magnetic resonance imaging data (NTD = 266) from Autism Brain Imaging Data Exchange (ABIDE) were used for subregion modeling. We proposed a distance-weighted sparse representation algorithm to construct brain functional networks. Functional subregions of M1 and S1 were identified through consensus clustering at the group level. Differences in the characteristics of functional subregions were analyzed, along with their correlation with clinical scores.

RESULTS

We observed symmetrical and continuous subregion organization from dorsal to ventral aspects in M1 and S1, with M1 subregions conforming to the functional pattern of the motor homunculus. Significant intergroup differences and clinical correlations were found in the dorsal and ventral aspects of M1 (p < 0.05/3, Bonferroni correction) and the ventromedial BA3 of S1 (p < 0.05/5). These functional characteristics were positively correlated with autism severity. All subregions showed significant results in the ROI-to-ROI intergroup differential analysis (p < 0.05/80).

LIMITATIONS

The generalizability of the segmentation model requires further evaluation.

CONCLUSIONS

This study highlights the significance of M1 and S1 in ASD treatment and may provide new insights into brain parcellation and the identification of therapeutic targets for ASD.

Prior-guided individualized thalamic parcellation based on local diffusion characteristics

Comprising numerous subnuclei, the thalamus intricately interconnects the cortex and subcortex, orchestrating various facets of brain functions. Extracting personalized parcellation patterns for these subnuclei is crucial, as different thalamic nuclei play varying roles in cognition and serve as therapeutic targets for neuromodulation. However, accurately delineating the thalamic nuclei boundary at the individual level is challenging due to intersubject variability. In this study, we proposed a prior-guided parcellation (PG-par) method to achieve robust individualized thalamic parcellation based on a central-boundary prior. We first constructed probabilistic atlas of thalamic nuclei using high-quality diffusion MRI datasets based on the local diffusion characteristics. Subsequently, high-probability voxels in the probabilistic atlas were utilized as prior guidance to train unique multiple classification models for each subject based on a multilayer perceptron. Finally, we employed the trained model to predict the parcellation labels for thalamic voxels and construct individualized thalamic parcellation. Through a test-retest assessment, the proposed prior-guided individualized thalamic parcellation exhibited excellent reproducibility and the capacity to detect individual variability. Compared with group atlas registration and individual clustering parcellation, the proposed PG-par demonstrated superior parcellation performance under different scanning protocols and clinic settings. Furthermore, the prior-guided individualized parcellation exhibited better correspondence with the histological staining atlas. The proposed prior-guided individualized thalamic parcellation method contributes to the personalized modeling of brain parcellation.

Eagle-449: A volumetric, whole-brain compilation of brain atlases for vestibular functional MRI research

Human vestibular processing involves distributed networks of cortical and subcortical regions which perform sensory and multimodal integrative functions. These functional hubs are also interconnected with areas subserving cognitive, affective, and body-representative domains. Analysis of these diverse components of the vestibular and vestibular-associated networks, and synthesis of their holistic functioning, is therefore vital to our understanding of the genesis of vestibular dysfunctions and aid treatment development. Novel neuroimaging methodologies, including functional and structural connectivity analyses, have provided important contributions in this area, but often require the use of atlases which are comprised of well-defined a priori regions of interest. Investigating vestibular dysfunction requires a more detailed atlas that encompasses cortical, subcortical, cerebellar, and brainstem regions. The present paper represents an effort to establish a compilation of existing, peer-reviewed brain atlases which collectively afford comprehensive coverage of these regions while explicitly focusing on vestibular substrates. It is expected that this compilation will be iteratively improved with additional contributions from researchers in the field.

Group-guided individual functional parcellation of the hippocampus and application to normal aging

Aging is closely associated with cognitive decline affecting attention, memory and executive functions. The hippocampus is the core brain area for human memory, learning, and cognition processing. To delineate the individual functional patterns of hippocampus is pivotal to reveal the neural basis of aging. In this study, we developed a group‐guided individual parcellation approach based on semisupervised affinity propagation clustering using the resting‐state functional magnetic resonance imaging to identify individual functional subregions of hippocampus and to identify the functional patterns of each subregion during aging. A three‐way group parcellation was yielded and was taken as prior information to guide individual parcellation of hippocampus into head, body, and tail in each subject. The superiority of individual parcellation of hippocampus is validated by higher intraregional functional similarities by compared to group‐level parcellation results. The individual variations of hippocampus were associated with coactivation patterns of three typical functions of hippocampus. Moreover, the individual functional connectivities of hippocampus subregions with predefined target regions could better predict age than group‐level functional connectivities. Our study provides a novel framework for individual brain functional parcellations, which may facilitate the future individual researches for brain cognitions and brain disorders and directing accurate neuromodulation.

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.

CBPtools: a Python package for regional connectivity-based parcellation

Regional connectivity-based parcellation (rCBP) is a widely used procedure for investigating the structural and functional differentiation within a region of interest (ROI) based on its long-range connectivity. No standardized software or guidelines currently exist for applying rCBP, making the method only accessible to those who develop their own tools. As such, there exists a discrepancy between the laboratories applying the procedure each with their own software solutions, making it difficult to compare and interpret the results. Here, we outline an rCBP procedure accompanied by an open source software package called CBPtools. CBPtools is a Python (version 3.5+) package that allows users to run an extensively evaluated rCBP analysis workflow on a given ROI. It currently supports two modalities: resting-state functional connectivity and structural connectivity based on diffusion-weighted imaging, along with support for custom connectivity matrices. Analysis parameters are customizable and the workflow can be scaled to a large number of subjects using a parallel processing environment. Parcellation results with corresponding validity metrics are provided as textual and graphical output. Thus, CBPtools provides a simple plug-and-play, yet customizable way to conduct rCBP analyses. By providing an open-source software we hope to promote reproducible and comparable rCBP analyses and, importantly, make the rCBP procedure readily available. Here, we demonstrate the utility of CBPtools using a voluminous data set on an average compute-cluster infrastructure by performing rCBP on three ROIs prominently featured in parcellation literature.

Parcellation of the human hippocampus based on gray matter volume covariance: Replicable results on healthy young adults

The hippocampus is a key brain region that participates in a range of cognitive and affective functions, and is involved in the etiopathogenesis of numerous neuropsychiatric disorders. The structural complexity and functional diversity of the hippocampus suggest the existence of structural and functional subdivisions within this structure. For the first time, we parcellated the human hippocampus with two independent data sets, each of which consisted of 198 T1-weighted structural magnetic resonance imaging (sMRI) images of healthy young subjects. The method was based on gray matter volume (GMV) covariance, which was quantified by a bivariate voxel-to-voxel linear correlation approach, as well as a multivariate masked independent component analysis approach. We subsequently interrogated the relationship between the GMV covariance patterns and the functional connectivity patterns of the hippocampal subregions using sMRI and resting-state functional MRI (fMRI) data from the same participants. Seven distinct GMV covariance-based subregions were identified for bilateral hippocampi, with robust reproducibility across the two data sets. We further demonstrated that the structural covariance patterns of the hippocampal subregions had a correspondence with the intrinsic functional connectivity patterns of these subregions. Together, our results provide a topographical configuration of the hippocampus with converging structural and functional support. The resulting subregions may improve our understanding of the hippocampal connectivity and functions at a subregional level, which provides useful parcellations and masks for future neuroscience and clinical research on the structural and/or functional connectivity of the hippocampus.

A Robust Reduced Rank Graph Regression Method for Neuroimaging Genetic Analysis

To characterize associations between genetic and neuroimaging data, a variety of analytic methods have been proposed in neuroimaging genetic studies. These methods have achieved promising performance by taking into account inherent correlation in either the neuroimaging data or the genetic data alone. In this study, we propose a novel robust reduced rank graph regression based method in a linear regression framework by considering correlations inherent in neuroimaging data and genetic data jointly. Particularly, we model the association analysis problem in a reduced rank regression framework with the genetic data as a feature matrix and the neuroimaging data as a response matrix by jointly considering correlations among the neuroimaging data as well as correlations between the genetic data and the neuroimaging data. A new graph representation of genetic data is adopted to exploit their inherent correlations, in addition to robust loss functions for both the regression and the data representation tasks, and a square-root-operator applied to the robust loss functions for achieving adaptive sample weighting. The resulting optimization problem is solved using an iterative optimization method whose convergence has been theoretically proved. Experimental results on the Alzheimer's Disease Neuroimaging Initiative (ADNI) dataset have demonstrated that our method could achieve competitive performance in terms of regression performance between brain structural measures and the Single Nucleotide Polymorphisms (SNPs), compared with state-of-the-art alternative methods.

Automatic Segmentation of Corpus Callosum in Midsagittal Based on Bayesian Inference Consisting of Sparse Representation Error and Multi-Atlas Voting

In this paper, we introduce a novel automatic method for Corpus Callosum (CC) in midsagittal plane segmentation. The robust segmentation of CC in midsagittal plane is key role for quantitative study of structural features of CC associated with various neurological disorder such as epilepsy, autism, Alzheimer's disease, and so on. Our approach is based on Bayesian inference using sparse representation and multi-atlas voting which both methods are used in various medical imaging, and show outstanding performance. Prior information in the proposed Bayesian inference is obtained from probability map generated from multi-atlas voting. The probability map contains the information of shape and location of CC of target image. Likelihood in the proposed Bayesian inference is obtained from gamma distribution function, generated from reconstruction errors (or sparse representation error), which are calculated in sparse representation of target patch using foreground dictionary and background dictionary each. Unlike the usual sparse representation method, we added gradient magnitude and gradient direction information to the patches of dictionaries and target, which had better segmentation performance than when not added. We compared three main segmentation results as follow: (1) the joint label fusion (JLF) method which is state-of-art method in multi-atlas voting based segmentation for evaluation of our method; (2) prior information estimated from multi-atlas voting only; (3) likelihood estimated from comparison of the reconstruction errors from sparse representation error only; (4) the proposed Bayesian inference. The methods were evaluated using two data sets of T1-weighted images, which one data set consists of 100 normal young subjects and the other data set consist of 25 normal old subjects and 22 old subjects with heavy drinker. In both data sets, the proposed Bayesian inference method has significantly the best segmentation performance than using each method separately.

Functional mapping and effective connectivity of the human operculum

The operculum, defined as the cortex adjacent to the insula, is a large structure encompassing three lobes, with a recognized role in a variety of neurologic and psychiatric conditions. Its complex functions include sensory, motor, autonomic and cognitive processing. In humans, these are extended with the addition of language. These functions are implemented by highly specialized neuronal populations and their widespread connections, which our study aims at mapping in detail. We studied a group of 31 patients that were explored with intracranial electrodes during the pre-surgical workup for drug-resistant epilepsy. We have selected the subset of contacts implanted in non-epileptogenic opercular cortex and we analyzed the neurophysiological and behavioral responses to direct electrical stimulation. The functional mapping was performed by applying 1 Hz and 50 Hz electrical stimulation on 252 contact pairs and recording the threshold for evoking clinical effects. The effective connectivity was assessed using cortico-cortical evoked potentials elicited by single-pulse electrical stimulation in a subset of 19 patients. The locations of the effects grouped in twelve distinct semiological classes were analyzed. The most frequent effects evoked by stimulation of the frontal operculum were language related (29%). The Rolandic area produced most often oropharyngeal symptoms (47%), the parietal operculum produced somatosensory effects (67%), while the temporal evoked auditory (58%) semiology. The connectivity pattern was complex, with these structures having widespread ipsilateral and contralateral projections. The local connections between the opercular subregions and with the insula, as well as with more distant areas like the cingulate gyrus, were distinguished by strength and between-subjects consistency. In conclusion, we demonstrate specific opercular functionality, distinct from the one of the insular cortex. The study is complemented by a literature review on the opercular functional connectome in human and non-human primates.

Structural impacts on the timing and amplitude of the negative BOLD response

The positive (PBR) and negative BOLD responses (NBR) arising in task-fMRI display varying magnitudes and dynamics across voxels. While the effects of structure, particularly of veins, on the PBR have been studied, little is known of NBR-structure relationships. Like the PBR, the NBR is often used as a surrogate marker of neuronal activation in both basic and clinical research and assessing its relationship with cortical structure may help interpret group differences. We therefore investigated how local structure affects BOLD amplitude and timing in PBR and NBR areas using multi-band fMRI during visual stimulation to obtain high temporal resolution (TR=0.45s) data combined with T1 imaging and susceptibility-weighted imaging (SWI) to quantify the local densities of gray/white matter and veins, respectively. In both PBRs and NBRs, larger venous density was consistently associated with larger BOLD amplitude and delay, up to 1-2s larger relative to areas devoid of large veins. Both binary and sinusoidal visual stimulus modulation yielded similar activation maps and results, suggesting that underlying vasculature affects PBR and NBR temporal dynamics in the same manner. Accounting for structural impacts on PBR and NBR magnitude and timing could help enhance activation map accuracy, better assess functional connectivity, and better characterize neurovascular coupling.

Making group inferences using sparse representation of resting-state functional mRI data with application to sleep deprivation

Past studies on drawing group inferences for functional magnetic resonance imaging (fMRI) data usually assume that a brain region is involved in only one functional brain network. However, recent evidence has demonstrated that some brain regions might simultaneously participate in multiple functional networks. Here, we presented a novel approach for making group inferences using sparse representation of resting-state fMRI data and its application to the identification of changes in functional networks in the brains of 37 healthy young adult participants after 36 h of sleep deprivation (SD) in contrast to the rested wakefulness (RW) stage. Our analysis based on group-level sparse representation revealed that multiple functional networks involved in memory, emotion, attention, and vigilance processing were impaired by SD. Of particular interest, the thalamus was observed to contribute to multiple functional networks in which differentiated response patterns were exhibited. These results not only further elucidate the impact of SD on brain function but also demonstrate the ability of the proposed approach to provide new insights into the functional organization of the resting-state brain by permitting spatial overlap between networks and facilitating the description of the varied relationships of the overlapping regions with other regions of the brain in the context of different functional systems. Hum Brain Mapp 38:4671-4689, 2017. © 2017 Wiley Periodicals, Inc.

Anatomical and functional organization of the human substantia nigra and its connections

We investigated the anatomical and functional organization of the human substantia nigra (SN) using diffusion and functional MRI data from the Human Connectome Project. We identified a tripartite connectivity-based parcellation of SN with a limbic, cognitive, motor arrangement. The medial SN connects with limbic striatal and cortical regions and encodes value (greater response to monetary wins than losses during fMRI), while the ventral SN connects with associative regions of cortex and striatum and encodes salience (equal response to wins and losses). The lateral SN connects with somatomotor regions of striatum and cortex and also encodes salience. Behavioral measures from delay discounting and flanker tasks supported a role for the value-coding medial SN network in decisional impulsivity, while the salience-coding ventral SN network was associated with motor impulsivity. In sum, there is anatomical and functional heterogeneity of human SN, which underpins value versus salience coding, and impulsive choice versus impulsive action.

Cross-cultural consistency and diversity in intrinsic functional organization of Broca's Region

As a core language area, Broca's region was consistently activated in a variety of language studies even across different language systems. Moreover, a high degree of structural and functional heterogeneity in Broca's region has been reported in many studies. This raised the issue of how the intrinsic organization of Broca's region effects by different language experiences in light of its subdivisions. To address this question, we used multi-center resting-state fMRI data to explore the cross-cultural consistency and diversity of Broca's region in terms of its subdivisions, connectivity patterns and modularity organization in Chinese and German speakers. A consistent topological organization of the 13 subdivisions within the extended Broca's region was revealed on the basis of a new in-vivo parcellation map, which corresponded well to the previously reported receptorarchitectonic map. Based on this parcellation map, consistent functional connectivity patterns and modularity organization of these subdivisions were found. Some cultural difference in the functional connectivity patterns was also found, for instance stronger connectivity in Chinese subjects between area 6v2 and the motor hand area, as well as higher correlations between area 45p and middle frontal gyrus. Our study suggests that a generally invariant organization of Broca's region, together with certain regulations of different language experiences on functional connectivity, might exists to support language processing in human brain.

The Human Brainnetome Atlas: A New Brain Atlas Based on Connectional Architecture

The human brain atlases that allow correlating brain anatomy with psychological and cognitive functions are in transition from ex vivo histology-based printed atlases to digital brain maps providing multimodal in vivo information. Many current human brain atlases cover only specific structures, lack fine-grained parcellations, and fail to provide functionally important connectivity information. Using noninvasive multimodal neuroimaging techniques, we designed a connectivity-based parcellation framework that identifies the subdivisions of the entire human brain, revealing the in vivo connectivity architecture. The resulting human Brainnetome Atlas, with 210 cortical and 36 subcortical subregions, provides a fine-grained, cross-validated atlas and contains information on both anatomical and functional connections. Additionally, we further mapped the delineated structures to mental processes by reference to the BrainMap database. It thus provides an objective and stable starting point from which to explore the complex relationships between structure, connectivity, and function, and eventually improves understanding of how the human brain works. The human Brainnetome Atlas will be made freely available for download at http://atlas.brainnetome.org, so that whole brain parcellations, connections, and functional data will be readily available for researchers to use in their investigations into healthy and pathological states.

Identifying functional subdivisions in the human brain using meta-analytic activation modeling-based parcellation

Parcellation of the human brain into fine-grained units by grouping voxels into distinct clusters has been an effective approach for delineating specific brain regions and their subregions. Published neuroimaging studies employing coordinate-based meta-analyses have shown that the activation foci and their corresponding behavioral categories may contain useful information about the anatomical-functional organization of brain regions. Inspired by these developments, we proposed a new parcellation scheme called meta-analytic activation modeling-based parcellation (MAMP) that uses meta-analytically obtained information. The raw meta data, including the experiments and the reported activation coordinates related to a brain region of interest, were acquired from the Brainmap database. Using this data, we first obtained the "modeled activation" pattern by modeling the voxel-wise activation probability given spatial uncertainty for each experiment that featured at least one focus within the region of interest. Then, we processed these "modeled activation" patterns across the experiments with a K-means clustering algorithm to group the voxels into different subregions. In order to verify the reliability of the method, we employed our method to parcellate the amygdala and the left Brodmann area 44 (BA44). The parcellation results were quite consistent with previous cytoarchitectonic and in vivo neuroimaging findings. Therefore, the MAMP proposed in the current study could be a useful complement to other methods for uncovering the functional organization of the human brain.