BSMac: a MATLAB toolbox implementing a Bayesian spatial model for brain activation and connectivity

HINT: A hierarchical independent component analysis toolbox for investigating brain functional networks using neuroimaging data

BACKGROUND

Independent component analysis (ICA) is a popular tool for investigating brain organization in neuroscience research. In fMRI studies, an important goal is to study how brain networks are modulated by subjects' clinical and demographic variables. Existing ICA methods and toolboxes don't incorporate subjects' covariates effects in ICA estimation of brain networks, which potentially leads to loss in accuracy and statistical power in detecting brain network differences between subjects' groups.

NEW METHOD

We introduce a Matlab toolbox, HINT (Hierarchical INdependent component analysis Toolbox), that provides a hierarchical covariate-adjusted ICA (hc-ICA) for modeling and testing covariate effects and generates model-based estimates of brain networks on both the population- and individual-level. HINT provides a user-friendly Matlab GUI that allows users to easily load images, specify covariate effects, monitor model estimation via an EM algorithm, specify hypothesis tests, and visualize results. HINT also has a command line interface which allows users to conveniently run and reproduce the analysis with a script.

COMPARISON TO EXISTING METHODS

HINT implements a new multi-level probabilistic ICA model for group ICA. It provides a statistically principled ICA modeling framework for investigating covariate effects on brain networks. HINT can also generate and visualize model-based network estimates for user-specified subject groups, which greatly facilitates group comparisons.

RESULTS

We demonstrate the steps and functionality of HINT with an fMRI example data to estimate treatment effects on brain networks while controlling for other covariates. Results demonstrate estimated brain networks and model-based comparisons between the treatment and control groups. In comparisons using synthetic fMRI data, HINT shows desirable statistical power in detecting group differences in networks especially in small sample sizes, while maintaining a low false positive rate. HINT also demonstrates similar or increased accuracy in reconstructing both population- and individual-level source signal maps as compared to some state-of-the-art group ICA methods.

CONCLUSION

HINT can provide a useful tool for both statistical and neuroscience researchers to evaluate and test differences in brain networks between subject groups.

Cortical Thinning and Cognitive Impairment in Parkinson's Disease without Dementia

Parkinson's disease (PD) is a progressive neurodegenerative disorder characterized clinically by motor dysfunction (bradykinesia, rigidity, tremor, and postural instability), and pathologically by the loss of dopaminergic neurons in the substantia nigra of the basal ganglia. Growing literature supports that cognitive deficits may also be present in PD, even in non-demented patients. Gray matter (GM) atrophy has been reported in PD and may be related to cognitive decline. This study investigated cortical thickness in non-demented PD subjects and elucidated its relationship to cognitive impairment using high-resolution T1-weighted brain MRI and comprehensive cognitive function scores from 71 non-demented PD and 48 control subjects matched for age, gender, and education. Cortical thickness was compared between groups using a flexible hierarchical multivariate Bayesian model, which accounts for correlations between brain regions. Correlation analyses were performed among brain areas and cognitive domains as well, which showed significant group differences in the PD population. Compared to Controls, PD subjects demonstrated significant age-adjusted cortical thinning predominantly in inferior and superior parietal areas and extended to superior frontal, superior temporal, and precuneus areas (posterior probability >0.9). Cortical thinning was also found in the left precentral and lateral occipital, and right postcentral, middle frontal, and fusiform regions (posterior probability >0.9). PD patients showed significantly reduced cognitive performance in executive function, including set shifting (p = 0.005) and spontaneous flexibility (p = 0.02), which were associated with the above cortical thinning regions (p < 0.05).

Brain Imaging Analysis

The increasing availability of brain imaging technologies has led to intense neuroscientific inquiry into the human brain. Studies often investigate brain function related to emotion, cognition, language, memory, and numerous other externally induced stimuli as well as resting-state brain function. Studies also use brain imaging in an attempt to determine the functional or structural basis for psychiatric or neurological disorders and, with respect to brain function, to further examine the responses of these disorders to treatment. Neuroimaging is a highly interdisciplinary field, and statistics plays a critical role in establishing rigorous methods to extract information and to quantify evidence for formal inferences. Neuroimaging data present numerous challenges for statistical analysis, including the vast amounts of data collected from each individual and the complex temporal and spatial dependence present. We briefly provide background on various types of neuroimaging data and analysis objectives that are commonly targeted in the field. We present a survey of existing methods targeting these objectives and identify particular areas offering opportunities for future statistical contribution.

Efficient posterior probability mapping using Savage-Dickey ratios

Statistical Parametric Mapping (SPM) is the dominant paradigm for mass-univariate analysis of neuroimaging data. More recently, a Bayesian approach termed Posterior Probability Mapping (PPM) has been proposed as an alternative. PPM offers two advantages: (i) inferences can be made about effect size thus lending a precise physiological meaning to activated regions, (ii) regions can be declared inactive. This latter facility is most parsimoniously provided by PPMs based on Bayesian model comparisons. To date these comparisons have been implemented by an Independent Model Optimization (IMO) procedure which separately fits null and alternative models. This paper proposes a more computationally efficient procedure based on Savage-Dickey approximations to the Bayes factor, and Taylor-series approximations to the voxel-wise posterior covariance matrices. Simulations show the accuracy of this Savage-Dickey-Taylor (SDT) method to be comparable to that of IMO. Results on fMRI data show excellent agreement between SDT and IMO for second-level models, and reasonable agreement for first-level models. This Savage-Dickey test is a Bayesian analogue of the classical SPM-F and allows users to implement model comparison in a truly interactive manner.

Predicting brain activity using a Bayesian spatial model

Increasing the clinical applicability of functional neuroimaging technology is an emerging objective, e.g. for diagnostic and treatment purposes. We propose a novel Bayesian spatial hierarchical framework for predicting follow-up neural activity based on an individual's baseline functional neuroimaging data. Our approach attempts to overcome some shortcomings of the modeling methods used in other neuroimaging settings, by borrowing strength from the spatial correlations present in the data. Our proposed methodology is applicable to data from various imaging modalities including functional magnetic resonance imaging and positron emission tomography, and we provide an illustration here using positron emission tomography data from a study of Alzheimer's disease to predict disease progression.