CanICA: Model-based extraction of reproducible group-level ICA patterns from fMRI time series

Spatial Independent Component Analysis (ICA) is an increasingly used data-driven method to analyze functional Magnetic Resonance Imaging (fMRI) data. To date, it has been used to extract meaningful patterns without prior information. However, ICA is not robust to mild data variation and remains a parameter-sensitive algorithm. The validity of the extracted patterns is hard to establish, as well as the significance of differences between patterns extracted from different groups of subjects. We start from a generative model of the fMRI group data to introduce a probabilistic ICA pattern-extraction algorithm, called CanICA (Canonical ICA). Thanks to an explicit noise model and canonical correlation analysis, our method is auto-calibrated and identifies the group-reproducible data subspace before performing ICA. We compare our method to state-of-the-art multi-subject fMRI ICA methods and show that the features extracted are more reproducible

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 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 non-uniformly 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

Role of Alpha Oscillations During Short Time Memory Task Investigated by Graph Based Partitioning

In this study, we investigate the clustering pattern of alpha band (8 Hz - 12 Hz) electroencephalogram (EEG) oscillations obtained from healthy individuals during a short time memory task with 3 different memory loads. The retention period during which subjects were asked to memorize a pattern in a square matrix is analyzed with a graph theoretical approach. The functional coupling among EEG electrodes are quantified via mutual information in the time-frequency plane. A spectral clustering algorithm followed by bootstrapping is used to parcellate memory related circuits and for identifying significant clusters in the brain. The main outcome of the study is that the size of the significant clusters formed by alpha oscillations decreases as the memory load increases. This finding corroborates the active inhibition hypothesis about alpha oscillations

Which fMRI clustering gives good brain parcellations?

International audienceAnalysis and interpretation of neuroimaging data often require one to divide the brain into a number of regions, or parcels, with homogeneous characteristics, be these regions defined in the brain volume or on on the cortical surface. While predefined brain atlases do not adapt to the signal in the individual subjects images, parcellation approaches use brain activity (e.g. found in some functional contrasts of interest) and clustering techniques to define regions with some degree of signal homogeneity. In this work, we address the question of which clustering technique is appropriate and how to optimize the corresponding model. We use two principled criteria: goodness of fit (accuracy), and reproducibility of the parcellation across bootstrap samples. We study these criteria on both simulated and two task-based functional Magnetic Resonance Imaging datasets for the Ward, spectral and K-means clustering algorithms. We show that in general Ward's clustering performs better than alternative methods with regard to reproducibility and accuracy and that the two criteria diverge regarding the preferred models (reproducibility leading to more conservative solutions), thus deferring the practical decision to a higher level alternative, namely the choice of a trade-off between accuracy and stability

Improving Reliability of Subject-Level Resting-State fMRI Parcellation with Shrinkage Estimators

A recent interest in resting state functional magnetic resonance imaging (rsfMRI) lies in subdividing the human brain into anatomically and functionally distinct regions of interest. For example, brain parcellation is often used for defining the network nodes in connectivity studies. While inference has traditionally been performed on group-level data, there is a growing interest in parcellating single subject data. However, this is difficult due to the low signal-to-noise ratio of rsfMRI data, combined with typically short scan lengths. A large number of brain parcellation approaches employ clustering, which begins with a measure of similarity or distance between voxels. The goal of this work is to improve the reproducibility of single-subject parcellation using shrinkage estimators of such measures, allowing the noisy subject-specific estimator to "borrow strength" in a principled manner from a larger population of subjects. We present several empirical Bayes shrinkage estimators and outline methods for shrinkage when multiple scans are not available for each subject. We perform shrinkage on raw intervoxel correlation estimates and use both raw and shrinkage estimates to produce parcellations by performing clustering on the voxels. Our proposed method is agnostic to the choice of clustering method and can be used as a pre-processing step for any clustering algorithm. Using two datasets---a simulated dataset where the true parcellation is known and is subject-specific and a test-retest dataset consisting of two 7-minute rsfMRI scans from 20 subjects---we show that parcellations produced from shrinkage correlation estimates have higher reliability and validity than those produced from raw estimates. Application to test-retest data shows that using shrinkage estimators increases the reproducibility of subject-specific parcellations of the motor cortex by up to 30%.Comment: body 21 pages, 11 figure

Learning and comparing functional connectomes across subjects

Functional connectomes capture brain interactions via synchronized fluctuations in the functional magnetic resonance imaging signal. If measured during rest, they map the intrinsic functional architecture of the brain. With task-driven experiments they represent integration mechanisms between specialized brain areas. Analyzing their variability across subjects and conditions can reveal markers of brain pathologies and mechanisms underlying cognition. Methods of estimating functional connectomes from the imaging signal have undergone rapid developments and the literature is full of diverse strategies for comparing them. This review aims to clarify links across functional-connectivity methods as well as to expose different steps to perform a group study of functional connectomes