Abstract. Inference of brain activation through the analysis of functional magnetic resonance imaging (fMRI) data is seriously confounded by the high level of noise in the observations. To mitigate the effects of noise, we propose incorporating anatomical connectivity into brain activation detection as motivated by how the functional integration of distinct brain areas is facilitated via neural fiber pathways. In this work, we formulate activation detection as a probabilistic graph-based segmentation problem with fiber networks estimated from diffusion MRI (dMRI) data serving as a prior. Our approach is reinforced with a data-driven scheme for refining the connectivity prior to reflect the fact that not all fibers are necessarily deployed during a given cognitive task as well as to account for false fiber tracts arising from limitations of dMRI tractography. Validating on real clinical data collected from 7 schizophrenia patients and 13 matched healthy controls, we show that incorporating anatomical connectivity significantly increases sensitivity in detecting task activation in controls compared to existing univariate techniques. Further, we illustrate how our model enables the detection of significant group activation differences between controls and patients that are missed with standard methods.Keywords: activation detection, connectivity, dMRI, fMRI, random walker
IntroductionFunctional magnetic resonance imaging (fMRI) has become the primary modality for studying human brain activity. To map brain regions to function, standard analysis models the fMRI observations at each voxel as a linear combination of expected temporal responses using the general linear model (GLM) [1]. This univariate approach does not model the integrative property of the brain, which is known to facilitate brain function [2]. To ameliorate this serious limitation, the use of local neighbourhood information has been proposed to regularize activation detection [3,4]. Although such methods help suppress false spatially-isolated activations by encouraging spatial continuity, they completely ignore long-range functional interactions. The incorporation of functional connectivity information into task activation detection has also been put forth [5], but the approach taken estimates both activation effects and functional connectivity from the same dataset, hence the information gain might be limited. Other works investigated the use of resting-state (RS) functional connectivity information to