Spatially Adaptive Mixture Modeling for Analysis of fMRI Time Series.

International audienceWithin-subject analysis in fMRI essentially addresses two problems, the detection of brain regions eliciting evoked activity and the estimation of the underlying dynamics. In [1, 2], a detection-estimation framework has been proposed to tackle these problems jointly, since they are connected to one another. In the Bayesian formalism, detection is achieved by modeling activating and non-activating voxels through independent mixture models (IMM) within each region while hemodynamic response estimation is performed at a regional scale in a nonparametric way. Instead of IMMs, in this paper we take advantage of spatial mixture models (SMM) for their non-linear spatial regularizing properties. The proposed method is unsupervised and spatially adaptive in the sense that the amount of spatial correlation is automatically tuned from the data and this setting automatically varies across brain regions. In addition, the level of regularization is specific to each experimental condition since both the signal-to-noise ratio and the activation pattern may vary across stimulus types in a given brain region. These aspects require the precise estimation of multiple partition functions of underlying Ising fields. This is addressed efficiently using first path sampling for a small subset of fields and then using a recently developed fast extrapolation technique for the large remaining set. Simulation results emphasize that detection relying on supervised SMM outperforms its IMM counterpart and that unsupervised spatial mixture models achieve similar results without any hand-tuning of the correlation parameter. On real datasets, the gain is illustrated in a localizer fMRI experiment: brain activations appear more spatially resolved using SMM in comparison with classical General Linear Model (GLM)-based approaches, while estimating a specific parcel-based HRF shape. Our approach therefore validates the treatment of unsmoothed fMRI data without fixed GLM definition at the subject level and makes also the classical strategy of spatial Gaussian filtering deprecated

Combining spatial priors and anatomical information for fMRI detection

In this paper, we analyze Markov Random Field (MRF) as a spatial regularizer in fMRI detection. The low signal-to-noise ratio (SNR) in fMRI images presents a serious challenge for detection algorithms, making regularization necessary to achieve good detection accuracy. Gaussian smoothing, traditionally employed to boost SNR, often produces over-smoothed activation maps. Recently, the use of MRF priors has been suggested as an alternative regularization approach. However, solving for an optimal configuration of the MRF is NP-hard in general. In this work, we investigate fast inference algorithms based on the Mean Field approximation in application to MRF priors for fMRI detection. Furthermore, we propose a novel way to incorporate anatomical information into the MRF-based detection framework and into the traditional smoothing methods. Intuitively speaking, the anatomical evidence increases the likelihood of activation in the gray matter and improves spatial coherency of the resulting activation maps within each tissue type. Validation using the receiver operating characteristic (ROC) analysis and the confusion matrix analysis on simulated data illustrates substantial improvement in detection accuracy using the anatomically guided MRF spatial regularizer. We further demonstrate the potential benefits of the proposed method in real fMRI signals of reduced length. The anatomically guided MRF regularizer enables significant reduction of the scan length while maintaining the quality of the resulting activation maps.National Institutes of Health (U.S.) (National Institute for Biomedical Imaging and Bioengineering (U.S.)/National Alliance for Medical Image Computing (U.S.) Grant U54-EB005149)National Science Foundation (U.S.) (Grant IIS 9610249)National Institutes of Health (U.S.) (National Center for Research Resources (U.S.)/Biomedical Informatics Research Network Grant U24-RR021382)National Institutes of Health (U.S.) (National Center for Research Resources (U.S.)/Neuroimaging Analysis Center (U.S.) Grant P41-RR13218)National Institutes of Health (U.S.) (National Institute of Neurological Disorders and Stroke (U.S.) Grant R01-NS051826)National Science Foundation (U.S.) (CAREER Grant 0642971)National Science Foundation (U.S.). Graduate Research FellowshipNational Center for Research Resources (U.S.) (FIRST-BIRN Grant)Neuroimaging Analysis Center (U.S.

Bayesian analysis of functional magnetic resonance imaging data with spatially varying auto‐regressive orders

Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/148250/1/rssc12320.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/148250/2/rssc12320_am.pd

Application and validation of spatial mixture modelling for the joint detection-estimation of brain activity in fMRI.

International audienceWithin-subject analysis in event-related functional Magnetic Resonance Imaging (fMRI) first relies on (i) a detection step to localize which parts of the brain are activated by a given stimulus type, and second on (ii) an estimation step to recover the temporal dynamics of the brain response. Recently, a Bayesian detection-estimation approach that jointly addresses (i)-(ii) has been proposed in [1]. This work is based on an independent mixture model (IMM) and provides both a spatial activity map and an estimate of brain dynamics. In [2], we accounted for spatial correlation using a spatial mixture model (SMM) based on a binary Markov random field. Here, we assess the SMM robustness and flexibility on simulations which diverge from the priors and the generative BOLD model and further extend comparison between SMM and IMM on real fMRI data, focusing on a region of interest in the auditory cortex

Doctor of Philosophy

dissertationFunctional magnetic resonance imaging (fMRI) measures the change of oxygen consumption level in the blood vessels of the human brain, hence indirectly detecting the neuronal activity. Resting-state fMRI (rs-fMRI) is used to identify the intrinsic functional patterns of the brain when there is no external stimulus. Accurate estimation of intrinsic activity is important for understanding the functional organization and dynamics of the brain, as well as differences in the functional networks of patients with mental disorders. This dissertation aims to robustly estimate the functional connectivities and networks of the human brain using rs-fMRI data of multiple subjects. We use Markov random field (MRF), an undirected graphical model to represent the statistical dependency among the functional network variables. Graphical models describe multivariate probability distributions that can be factorized and represented by a graph. By defining the nodes and the edges along with their weights according to our assumptions, we build soft constraints into the graph structure as prior information. We explore various approximate optimization methods including variational Bayesian, graph cuts, and Markov chain Monte Carlo sampling (MCMC). We develop the random field models to solve three related problems. In the first problem, the goal is to detect the pairwise connectivity between gray matter voxels in a rs-fMRI dataset of the single subject. We define a six-dimensional graph to represent our prior information that two voxels are more likely to be connected if their spatial neighbors are connected. The posterior mean of the connectivity variables are estimated by variational inference, also known as mean field theory in statistical physics. The proposed method proves to outperform the standard spatial smoothing and is able to detect finer patterns of brain activity. Our second work aims to identify multiple functional systems. We define a Potts model, a special case of MRF, on the network label variables, and define von Mises-Fisher distribution on the normalized fMRI signal. The inference is significantly more difficult than the binary classification in the previous problem. We use MCMC to draw samples from the posterior distribution of network labels. In the third application, we extend the graphical model to the multiple subject scenario. By building a graph including the network labels of both a group map and the subject label maps, we define a hierarchical model that has richer structure than the flat single-subject model, and captures the shared patterns as well as the variation among the subjects. All three solutions are data-driven Bayesian methods, which estimate model parameters from the data. The experiments show that by the regularization of MRF, the functional network maps we estimate are more accurate and more consistent across multiple sessions

Studying Resting State and Stimulus Induced BOLD Activity using the Generalized Ising Model and Independent Component Graph Analysis

Although many technical advancements have been made, neuroscientists still struggle to explain the underlying behaviour of how brain regions communicate with each other to form large-scale functional networks. functional Magnetic Resonance Imaging (fMRI) has been commonly used to investigate changes between brain regions over time using the Blood Oxygen Level Dependent (BOLD) signal. The goal of this thesis is to show the applicability of novel techniques and tools, such as the generalized Ising model (GIM) and the independent component graph analysis (GraphICA), to obtain information on the functional connectivity of populations with altered perception of consciousness. The GIM was used to model brain activity in healthy brains during various stages of consciousness, as induced by an anesthetic agent, propofol, in the auditory paradigm. GraphICA, a tool that combines ICA and graph theory was used to investigate the functional connectivity of resting state networks (RSNs) in patients with altered perception caused by tinnitus and in patients with altered states of consciousness caused by severe brain injury. For the tinnitus patients, we examined whether a correlation exists between tinnitus behavioural scores and functional brain connectivity of RSNs. Moreover, for the severely brain injured patients, a multimodal neuroimaging approach with hybrid FDG-PET/MRI was implemented to study the functional connectivity changes of the RSNs. The GIM simulated with an external field was able to model the brain activity at different levels of consciousness under naturalistic stimulation, at a temperature in the super critical regime. Further, a strong correlation was observed between tinnitus distress and the connectivity pattern within and between the right executive control network and the other RSNs. This suggests that tinnitus distress is the consequence of a hyperactive attention condition. A variability was observed in the appearance and temporal/spatial patterns of RSNs among the two resting state fMRI acquisitions acquired within the same scanning session of the severely brain injured patients. This suggests the need for new strategies to be developed in order to pick the best RSN from each acquisition. Overall, this work demonstrated that the GIM and GraphICA were promising tools to investigate brain activity of populations with altered perception of consciousness and in future can be extended to investigate different neurological populations

Spatio-temporal wavelet regularization for parallel MRI reconstruction: application to functional MRI

Parallel MRI is a fast imaging technique that enables the acquisition of highly resolved images in space or/and in time. The performance of parallel imaging strongly depends on the reconstruction algorithm, which can proceed either in the original k-space (GRAPPA, SMASH) or in the image domain (SENSE-like methods). To improve the performance of the widely used SENSE algorithm, 2D- or slice-specific regularization in the wavelet domain has been deeply investigated. In this paper, we extend this approach using 3D-wavelet representations in order to handle all slices together and address reconstruction artifacts which propagate across adjacent slices. The gain induced by such extension (3D-Unconstrained Wavelet Regularized -SENSE: 3D-UWR-SENSE) is validated on anatomical image reconstruction where no temporal acquisition is considered. Another important extension accounts for temporal correlations that exist between successive scans in functional MRI (fMRI). In addition to the case of 2D+t acquisition schemes addressed by some other methods like kt-FOCUSS, our approach allows us to deal with 3D+t acquisition schemes which are widely used in neuroimaging. The resulting 3D-UWR-SENSE and 4D-UWR-SENSE reconstruction schemes are fully unsupervised in the sense that all regularization parameters are estimated in the maximum likelihood sense on a reference scan. The gain induced by such extensions is illustrated on both anatomical and functional image reconstruction, and also measured in terms of statistical sensitivity for the 4D-UWR-SENSE approach during a fast event-related fMRI protocol. Our 4D-UWR-SENSE algorithm outperforms the SENSE reconstruction at the subject and group levels (15 subjects) for different contrasts of interest (eg, motor or computation tasks) and using different parallel acceleration factors (R=2 and R=4) on 2x2x3mm3 EPI images.Comment: arXiv admin note: substantial text overlap with arXiv:1103.353

Physiologically Informed Bayesian Analysis of ASL fMRI Data

Arterial Spin Labelling (ASL) functional Magnetic Resonance Imaging (fMRI) data provides a quantitative measure of blood perfusion, that can be correlated to neuronal activation. In contrast to BOLD measure, it is a direct measure of cerebral blood flow. However, ASL data has a lower SNR and resolution so that the recovery of the perfusion response of interest suffers from the contamination by a stronger hemodynamic component in the ASL signal. In this work we consider a model of both hemodynamic and perfusion components within the ASL signal. A physiological link between these two components is analyzed and used for a more accurate estimation of the perfusion response function in particular in the usual ASL low SNR conditions