Juha Pajula

A versatile software package for inter-subject correlation based analyses of fMRI

In the inter-subject correlation (ISC) based analysis of the functional magnetic resonance imaging (fMRI) data, the extent of shared processing across subjects during the experiment is determined by calculating correlation coefficients between the fMRI time series of the subjects in the corresponding brain locations. This implies that ISC can be used to analyze fMRI data without explicitly modeling the stimulus and thus ISC is a potential method to analyze fMRI data acquired under complex naturalistic stimuli. Despite of the suitability of ISC based approach to analyze complex fMRI data, no generic software tools have been made available for this purpose, limiting a widespread use of ISC based analysis techniques among neuroimaging community. In this paper, we present a graphical user interface (GUI) based software package, ISC Toolbox, implemented in Matlab for computing various ISC based analyses. Many advanced computations such as comparison of ISCs between different stimuli, time window ISC, and inter-subject phase synchronization are supported by the toolbox. The analyses are coupled with re-sampling based statistical inference. The ISC based analyses are data and computation intensive and the ISC toolbox is equipped with mechanisms to execute the parallel computations in a cluster environment automatically and with an automatic detection of the cluster environment in use. Currently, SGE-based (Oracle Grid Engine, Son of a Grid Engine, or Open Grid Scheduler) and Slurm environments are supported. In this paper, we present a detailed account on the methods behind the ISC Toolbox, the implementation of the toolbox and demonstrate the possible use of the toolbox by summarizing selected example applications. We also report the computation time experiments both using a single desktop computer and two grid environments demonstrating that parallelization effectively reduces the computing time. The ISC Toolbox is available in https://code.google.com/p/isc-toolbox/

Inter-Subject Correlation in fMRI: Method Validation against Stimulus-Model Based Analysis

Within functional magnetic resonance imaging (fMRI), the use of the traditional general linear model (GLM) based analysis methods is often restricted to strictly controlled research setups requiring a parametric activation model. Instead, Inter-Subject Correlation (ISC) method is based on voxel-wise correlation between the time series of the subjects, which makes it completely non-parametric and thus suitable for naturalistic stimulus paradigms such as movie watching. In this study, we compared an ISC based analysis results with those of a GLM based in five distinct controlled research setups. We used International Consortium for Brain Mapping functional reference battery (FRB) fMRI data available from the Laboratory of Neuro Imaging image data archive. The selected data included measurements from 37 right-handed subjects, who all had performed the same five tasks from FRB. The GLM was expected to locate activations accurately in FRB data and thus provide good grounds for investigating relationship between ISC and stimulus induced fMRI activation. The statistical maps of ISC and GLM were compared with two measures. The first measure was the Pearsons correlation between the non-thresholded ISC test-statistics and absolute values of the GLM Z-statistics. The average correlation value over five tasks was 0.74. The second was the Dice index between the activation regions of the methods. The average Dice value over the tasks and three threshold levels was 0.73. The results of this study indicated how the data driven ISC analysis found the same foci as the model-based GLM analysis. The agreement of the results is highly interesting, because ISC is applicable in situations where GLM is not suitable, for example, when analyzing data from a naturalistic stimuli experiment.

Effects of spatial smoothing on inter-subject correlation based analysis of FMRI.

This study evaluates the effects of spatial smoothing on inter-subject correlation (ISC) analysis for FMRI data using the traditional model based analysis as a reference. So far within ISC analysis the effects of smoothing have not been studied systematically and linear Gaussian filters with varying kernel widths have been used without better knowledge about the effects of filtering. Instead, with the traditional general linear model (GLM) based analysis, the effects of smoothing have been studied extensively. In this study, ISC and GLM analyses were computed with two experimental and one simulated block-design datasets. The test statistics and the detected activation areas were compared numerically with correlation and Dice similarity measures, respectively. The study verified that (1) the choice of the filter substantially affected the activations detected by ISC analysis, (2) the detected activations according to ISC and GLM methods were highly similar regardless of the smoothing kernel and (3) the effect of spatial smoothing was mildly smaller on ISC than GLM analysis. Our results indicated that a good selection of the full width at half maximum of the Gaussian smoothing kernel for ISC was slightly larger than double the original voxel size.

Functional Brain Segmentation Using Inter-Subject Correlation in fMRI

The human brain continuously processes massive amounts of rich sensory information. To better understand such highly complex brain processes, modern neuroimaging studies are increasingly utilizing experimental setups that better mimic daily‐life situations. A new exploratory data‐analysis approach, functional segmentation inter‐subject correlation analysis (FuSeISC), was proposed to facilitate the analysis of functional magnetic resonance (fMRI) data sets collected in these experiments. The method provides a new type of functional segmentation of brain areas, not only characterizing areas that display similar processing across subjects but also areas in which processing across subjects is highly variable. FuSeISC was tested using fMRI data sets collected during traditional block‐design stimuli (37 subjects) as well as naturalistic auditory narratives (19 subjects). The method identified spatially local and/or bilaterally symmetric clusters in several cortical areas, many of which are known to be processing the types of stimuli used in the experiments. The method is not only useful for spatial exploration of large fMRI data sets obtained using naturalistic stimuli, but also has other potential applications, such as generation of a functional brain atlases including both lower‐ and higher‐order processing areas. Finally, as a part of FuSeISC, a criterion‐based sparsification of the shared nearest‐neighbor graph was proposed for detecting clusters in noisy data. In the tests with synthetic data, this technique was superior to well‐known clustering methods, such as Wards method, affinity propagation, and K‐means ++ . Hum Brain Mapp 38:2643–2665, 2017.

How many is enough? effect of sample size in inter-subject correlation analysis of fMRI

Inter-subject correlation (ISC) is a widely used method for analyzing functional magnetic resonance imaging (fMRI) data acquired during naturalistic stimuli. A challenge in ISC analysis is to define the required sample size in the way that the results are reliable. We studied the effect of the sample size on the reliability of ISC analysis and additionally addressed the following question: How many subjects are needed for the ISC statistics to converge to the ISC statistics obtained using a large sample? The study was realized using a large block design data set of 130 subjects. We performed a split-half resampling based analysis repeatedly sampling two nonoverlapping subsets of 10–65 subjects and comparing the ISC maps between the independent subject sets. Our findings suggested that with 20 subjects, on average, the ISC statistics had converged close to a large sample ISC statistic with 130 subjects. However, the split-half reliability of unthresholded and thresholded ISC maps improved notably when the number of subjects was increased from 20 to 30 or more.

Brain hemodynamic activity during viewing and re-viewing of comedy movies explained by experienced humor.

Humor is crucial in human social interactions. To study the underlying neural processes, three comedy clips were shown twice to 20 volunteers during functional magnetic resonance imaging (fMRI). Inter-subject similarities in humor ratings, obtained immediately after fMRI, explained inter-subject correlation of hemodynamic activity in right frontal pole and in a number of other brain regions. General linear model analysis also indicated activity in right frontal pole, as well as in additional cortical areas and subcortically in striatum, explained by humorousness. The association of the right frontal pole with experienced humorousness is a novel finding, which might be related to humor unfolding over longer time scales in the movie clips. Specifically, frontal pole has been shown to exhibit longer temporal receptive windows than, e.g., sensory areas, which might have enabled processing of humor in the clips based on holding information and reinterpreting that in light of new information several (even tens of) seconds later. As another novel finding, medial and lateral prefrontal areas, frontal pole, posterior-inferior temporal areas, posterior parietal areas, posterior cingulate, striatal structures and amygdala showed reduced activity upon re-viewing of the clips, suggesting involvement in processing of humor related to novelty of the comedic events.

Platform for systems medicine research and diagnostic applications in psychotic disorders—The METSY project

Psychotic disorders are associated with metabolic abnormalities including alterations in glucose and lipid metabolism. A major challenge in the treatment of psychosis is to identify patients with vulnerable metabolic profiles who may be at risk of developing cardiometabolic co-morbidities. It is established that both central and peripheral metabolic organs use lipids to control energy balance and regulate peripheral insulin sensitivity. The endocannabinoid system, implicated in the regulation of glucose and lipid metabolism, has been shown to be dysregulated in psychosis. It is currently unclear how these endocannabinoid abnormalities relate to metabolic changes in psychosis. Here we review recent research in the field of metabolic co-morbidities in psychotic disorders as well as the methods to study them and potential links to the endocannabinoid system. We also describe the bioinformatics platforms developed in the EU project METSY for the investigations of the biological etiology in patients at risk of psychosis and in first episode psychosis patients. The METSY project was established with the aim to identify and evaluate multi-modal peripheral and neuroimaging markers that may be able to predict the onset and prognosis of psychiatric and metabolic symptoms in patients at risk of developing psychosis and first episode psychosis patients. Given the intrinsic complexity and widespread role of lipid metabolism, a systems biology approach which combines molecular, structural and functional neuroimaging methods with detailed metabolic characterisation and multi-variate network analysis is essential in order to identify how lipid dysregulation may contribute to psychotic disorders. A decision support system, integrating clinical, neuropsychological and neuroimaging data, was also developed in order to aid clinical decision making in psychosis. Knowledge of common and specific mechanisms may aid the etiopathogenic understanding of psychotic and metabolic disorders, facilitate early disease detection, aid treatment selection and elucidate new targets for pharmacological treatments.

Comparing fMRI inter-subject correlations between groups

Inter-subject correlation (ISC) based analysis is a conceptually simple approach to analyze functional magnetic resonance imaging (fMRI) data acquired under naturalistic stimuli such as a movie. We describe and validate the statistical approaches for comparing ISCs between two groups of subjects implemented in the ISC toolbox, which is an open source software package for ISC-based analysis of fMRI data. The approaches are based on permutation tests. We validated the approaches using five different data sets from the ICBM functional reference battery tasks. In these experiments, we created two matched groups of subjects and assumed that no group difference exists. Based on the experiments, we recommend the usage of subject-wise permutations, instead of element-wise permutations following Chen et al. (2016). However, we observed that the null-distributions should be voxel-specific and not based on pooling all voxels across the brain as is typical in fMRI. This was the case even if studentized permutation tests were used. Additionally, we experimented with an fMRI dataset acquired using a dance movie stimulus for comparison of the group of adult males in autism spectrum to the matched control group. The experiment confirmed the differences between voxel-based permutation tests and global model based permutation tests.