Diagnosis and exploration of massively univariate neuroimaging models

Grey-matter correlates of empathy in 4-Repeat Tauopathies

Loss of empathy is an early and central symptom of frontotemporal lobar degeneration spectrum diseases. We aimed to investigate the topographical distribution of morphometric brain changes associated with empathy in Progressive Supranuclear Palsy (PSP) and Corticobasal Syndrome (CBS) patients. Twenty-seven participants with CBS and 31 with PSP were evaluated using Interpersonal Reactivity Index scales in correlation with gray matter atrophy using a voxel-based morphometry approach. Lower levels of empathy were associated with an increased atrophy in fronto-temporal cortical structures. At subcortical level, empathy scores were positively correlated with gray matter volume in the amygdala, hippocampus and the cerebellum. These findings allow to extend the traditional cortico-centric view of cognitive empathy to the cerebellar regions in patients with neurodegenerative disorders and suggest that the cerebellum may play a more prominent role in social cognition than previously appreciated.

Identifying the neural correlates of anticipatory postural control: A novel fMRI paradigm

Altered postural control in the trunk/hip musculature is a characteristic of multiple neurological and musculoskeletal conditions. Previously it was not possible to determine if altered cortical and subcortical sensorimotor brain activation underlies impairments in postural control. This study used a novel fMRI-compatible paradigm to identify the brain activation associated with postural control in the trunk and hip musculature. BOLD fMRI imaging was conducted as participants performed two versions of a lower limb task involving lifting the left leg to touch the foot to a target. For the supported leg raise (SLR) the leg is raised from the knee while the thigh remains supported. For the unsupported leg raise (ULR) the leg is raised from the hip, requiring postural muscle activation in the abdominal/hip extensor musculature. Significant brain activation during the SLR task occurred predominantly in the right primary and secondary sensorimotor cortical regions. Brain activation during the ULR task occurred bilaterally in the primary and secondary sensorimotor cortical regions, as well as cerebellum and putamen. In comparison with the SLR, the ULR was associated with significantly greater activation in the right premotor/SMA, left primary motor and cingulate cortices, primary somatosensory cortex, supramarginal gyrus/parietal operculum, superior parietal lobule, cerebellar vermis, and cerebellar hemispheres. Cortical and subcortical regions activated during the ULR, but not during the SLR, were consistent with the planning, and execution of a task involving multisegmental, bilateral postural control. Future studies using this paradigm will determine mechanisms underlying impaired postural control in patients with neurological and musculoskeletal dysfunction.

Highlight results, don't hide them: Enhance interpretation, reduce biases and improve reproducibility

Most neuroimaging studies display results that represent only a tiny fraction of the collected data. While it is conventional to present "only the significant results" to the reader, here we suggest that this practice has several negative consequences for both reproducibility and understanding. This practice hides away most of the results of the dataset and leads to problems of selection bias and irreproducibility, both of which have been recognized as major issues in neuroimaging studies recently. Opaque, all-or-nothing thresholding, even if well-intentioned, places undue influence on arbitrary filter values, hinders clear communication of scientific results, wastes data, is antithetical to good scientific practice, and leads to conceptual inconsistencies. It is also inconsistent with the properties of the acquired data and the underlying biology being studied. Instead of presenting only a few statistically significant locations and hiding away the remaining results, studies should "highlight" the former while also showing as much as possible of the rest. This is distinct from but complementary to utilizing data sharing repositories: the initial presentation of results has an enormous impact on the interpretation of a study. We present practical examples and extensions of this approach for voxelwise, regionwise and cross-study analyses using publicly available data that was analyzed previously by 70 teams (NARPS; Botvinik-Nezer, et al., 2020), showing that it is possible to balance the goals of displaying a full set of results with providing the reader reasonably concise and "digestible" findings. In particular, the highlighting approach sheds useful light on the kind of variability present among the NARPS teams' results, which is primarily a varied strength of agreement rather than disagreement. Using a meta-analysis built on the informative "highlighting" approach shows this relative agreement, while one using the standard "hiding" approach does not. We describe how this simple but powerful change in practice-focusing on highlighting results, rather than hiding all but the strongest ones-can help address many large concerns within the field, or at least to provide more complete information about them. We include a list of practical suggestions for results reporting to improve reproducibility, cross-study comparisons and meta-analyses.

autohrf-an R package for generating data-informed event models for general linear modeling of task-based fMRI data

The analysis of task-related fMRI data at the level of individual participants is commonly based on general linear modeling (GLM), which allows us to estimate the extent to which the BOLD signal can be explained by the task response predictors specified in the event model. The predictors are constructed by convolving the hypothesized time course of neural activity with an assumed hemodynamic response function (HRF). However, our assumptions about the components of brain activity, including their onset and duration, may be incorrect. Their timing may also differ across brain regions or from person to person, leading to inappropriate or suboptimal models, poor fit of the model to actual data, and invalid estimates of brain activity. Here, we present an approach that uses theoretically driven models of task response to define constraints on which the final model is computationally derived using actual fMRI data. Specifically, we developed autohrf-an R package that enables the evaluation and data-driven estimation of event models for GLM analysis. The highlight of the package is the automated parameter search that uses genetic algorithms to find the onset and duration of task predictors that result in the highest fitness of GLM based on the fMRI signal under predefined constraints. We evaluated the usefulness of the autohrf package on two original datasets of task-related fMRI activity, a slow event-related spatial working memory study and a mixed state-item study using the flanker task, and on a simulated slow event-related working memory data. Our results suggest that autohrf can be used to efficiently construct and evaluate better task-related brain activity models to gain a deeper understanding of BOLD task response and improve the validity of model estimates. Our study also highlights the sensitivity of fMRI analysis with GLM to precise event model specification and the need for model evaluation, especially in complex and overlapping event designs.

Deep networks for behavioral variant frontotemporal dementia identification from multiple acquisition sources

Behavioral variant frontotemporal dementia (bvFTD) is a neurodegenerative syndrome whose clinical diagnosis remains a challenging task especially in the early stage of the disease. Currently, the presence of frontal and anterior temporal lobe atrophies on magnetic resonance imaging (MRI) is part of the diagnostic criteria for bvFTD. However, MRI data processing is usually dependent on the acquisition device and mostly require human-assisted crafting of feature extraction. Following the impressive improvements of deep architectures, in this study we report on bvFTD identification using various classes of artificial neural networks, and present the results we achieved on classification accuracy and obliviousness on acquisition devices using extensive hyperparameter search. In particular, we will demonstrate the stability and generalization of different deep networks based on the attention mechanism, where data intra-mixing confers models the ability to identify the disorder even on MRI data in inter-device settings, i.e., on data produced by different acquisition devices and without model fine tuning, as shown from the very encouraging performance evaluations that dramatically reach and overcome the 90% value on the AuROC and balanced accuracy metrics.

The cluster depth tests: Toward point-wise strong control of the family-wise error rate in massively univariate tests with application to M/EEG

The cluster mass test has been widely used for massively univariate tests in M/EEG, fMRI and, recently, pupillometry analysis. It is a powerful method for detecting effects while controlling weakly the family-wise error rate (FWER), although its correct interpretation can only be performed at the cluster level without any point-wise conclusion. It implies that the discoveries of a cluster mass test cannot be precisely localized in time or in space. We propose a new multiple comparisons procedure, the cluster depth tests, that both controls the FWER while allowing an interpretation at the time point level. We show the conditions for a strong control of the FWER, and a simulation study shows that the cluster depth tests achieve large power and guarantee the FWER even in the presence of physiologically plausible effects. By having an interpretation at the time point/voxel level, the cluster depth tests make it possible to take full advantage of the high temporal resolution of EEG recording and give a precise timing of the start and end of the significant effects.

Functional dynamics of dopamine synthesis during monetary reward and punishment processing

The assessment of dopamine release with the PET competition model is thoroughly validated but entails disadvantages for the investigation of cognitive processes. We introduce a novel approach incorporating 6-[¹⁸F]FDOPA uptake as index of the dynamic regulation of dopamine synthesis enzymes by neuronal firing. The feasibility of this approach is demonstrated by assessing widely described sex differences in dopamine neurotransmission. Reward processing was behaviorally investigated in 36 healthy participants, of whom 16 completed fPET and fMRI during the monetary incentive delay task. A single 50 min fPET acquisition with 6-[¹⁸F]FDOPA served to quantify task-specific changes in dopamine synthesis. In men monetary gain induced stronger increases in ventral striatum dopamine synthesis than loss. Interestingly, the opposite effect was discovered in women. These changes were further associated with reward (men) and punishment sensitivity (women). As expected, fMRI showed robust task-specific neuronal activation but no sex difference. Our findings provide a neurobiological basis for known behavioral sex differences in reward and punishment processing, with important implications in psychiatric disorders showing sex-specific prevalence, altered reward processing and dopamine signaling. The high temporal resolution and magnitude of task-specific changes make fPET a promising tool to investigate functional neurotransmitter dynamics during cognitive processing and in brain disorders.

Contribution of the multi-echo approach in accelerated functional magnetic resonance imaging multiband acquisition

We wanted to verify the effect of combining multi‐echo (ME) functional magnetic resonance imaging (fMRI) with slice acceleration in simultaneous multi‐slice acquisition. The aim was to shed light on the benefits of multiple echoes for various acquisition settings, especially for levels of slice acceleration and flip angle. Whole‐brain ME fMRI data were obtained from 26 healthy volunteers (using three echoes; seven runs with slice acceleration 1, 4, 6, and 8; and two different flip angles for each of the first three acceleration factors) and processed as single‐echo (SE) data and ME data based on optimal combinations weighted by the contrast‐to‐noise ratio. Global metrics (temporal signal‐to‐noise ratio, signal‐to‐noise separation, number of active voxels, etc.) and local characteristics in regions of interest were used to evaluate SE and ME data. ME results outperformed SE results in all runs; the differences became more apparent for higher acceleration, where a significant decrease in data quality is observed. ME fMRI can improve the observed data quality metrics over SE fMRI for a wide range of accelerated fMRI acquisitions.

A Bayesian group lasso classification for ADNI volumetrics data

The primary objective of this paper is to develop a statistically valid classification procedure for analyzing brain image volumetrics data obtained from the Alzheimer's Disease Neuroimaging Initiative (ADNI) in elderly subjects with cognitive impairments. The Bayesian group lasso method thereby proposed for logistic regression efficiently selects an optimal model with the use of a spike and slab type prior. This method selects groups of attributes of a brain subregion encouraged by the group lasso penalty. We conduct simulation studies for high- and low-dimensional scenarios where our method is always able to select the true parameters that are truly predictive among a large number of parameters. The method is then applied on dichotomous response ADNI data which selects predictive atrophied brain regions and classifies Alzheimer's disease patients from healthy controls. Our analysis is able to give an accuracy rate of 80% for classifying Alzheimer's disease. The suggested method selects 29 brain subregions. The medical literature indicates that all these regions are associated with Alzheimer's patients. The Bayesian method of model selection further helps selecting only the subregions that are statistically significant, thus obtaining an optimal model.

Neuroplastic effects of a selective serotonin reuptake inhibitor in relearning and retrieval

Animal studies using selective serotonin reuptake inhibitors (SSRIs) and learning paradigms have demonstrated that serotonin is important for flexibility in executive functions and learning. SSRIs might facilitate relearning through neuroplastic processes and thus exert their clinical effects in psychiatric diseases where cognitive functioning is affected. However, translation of these mechanisms to humans is missing. In this randomized placebo-controlled trial, we assessed functional brain activation during learning and memory retrieval in healthy volunteers performing associative learning tasks aiming to translate facilitated relearning by SSRIs.

To this extent, seventy-six participants underwent three MRI scanning sessions: (1) at baseline, (2) after three weeks of daily associative learning and subsequent retrieval (face-matching or Chinese character–noun matching) and (3) after three weeks of relearning under escitalopram (10 mg/day) or placebo. Associative learning and retrieval tasks were performed during each functional MRI (fMRI) session. Statistical modeling was done using a repeated-measures ANOVA, to test for content-by-treatment-by-time interaction effects.

During the learning task, a significant substance-by-time interaction was found in the right insula showing a greater deactivation in the SSRI cohort after 21 days of relearning compared to the learning phase. In the retrieval task, there was a significant content-by-time interaction in the left angular gyrus (AG) with an increased activation in face-matching compared to Chinese-character matching for both learning and relearning phases. A further substance-by-time interaction was found in task performance after 21 days of relearning, indicating a greater decrease of performance in the placebo group.

Our findings that escitalopram modulate insula activation demonstrates successful translation of relearning as a mechanism of SSRIs in human. Furthermore, we show that the left AG is an active component of correct memory retrieval, which coincides with previous literature. We extend the function of this region by demonstrating its activation is not only stimulus dependent but also time constrained. Finally, we were able to show that escitalopram aids in relearning, irrespective of content.

Sailing in rough waters: Examining volatility of fMRI noise

Background

The assumption that functional magnetic resonance imaging (fMRI) noise has constant volatility has recently been challenged by studies examining heteroscedasticity arising from head motion and physiological noise. The present study builds on this work using latest methods from the field of financial mathematics to model fMRI noise volatility.

Methods

Multi-echo phantom and human fMRI scans were used and realised volatility was estimated. The Hurst parameter H ∈ (0, 1), which governs the roughness/irregularity of realised volatility time series, was estimated. Calibration of H was performed pathwise, using well-established neural network calibration tools.

Results

In all experiments the volatility calibrated to values within the rough case, H < 0.5, and on average fMRI noise was very rough with 0.03 < H < 0.05. Some edge effects were also observed, whereby H was larger near the edges of the phantoms.

Discussion

The findings suggest that fMRI volatility is not only non-constant, but also substantially more irregular than a standard Brownian motion. Thus, further research is needed to examine the impact such pronounced oscillations in the volatility of fMRI noise have on data analyses.

Multivariate pattern classification on BOLD activation pattern induced by deep brain stimulation in motor, associative, and limbic brain networks

Deep brain stimulation (DBS) has been shown to be an effective treatment for movement disorders and it is now being extended to the treatment of psychiatric disorders. Functional magnetic resonance imaging (fMRI) studies indicate that DBS stimulation targets dependent brain network effects, in networks that respond to stimulation. Characterizing these patterns is crucial for linking DBS-induced therapeutic and adverse effects. Conventional DBS-fMRI, however, lacks the sensitivity needed for decoding multidimensional information such as spatially diffuse patterns. We report here on the use of a multivariate pattern analysis (MVPA) to demonstrate that stimulation of three DBS targets (STN, subthalamic nucleus; GPi, globus pallidus internus; NAc, nucleus accumbens) evoked a sufficiently distinctive blood-oxygen-level-dependent (BOLD) activation in swine brain. The findings indicate that STN and GPi evoke a similar motor network pattern, while NAc shows a districted associative and limbic pattern. The findings show that MVPA could be effectively applied to overlapping or sparse BOLD patterns which are often found in DBS. Future applications are expected employ MVPA fMRI to identify the proper stimulation target dependent brain circuitry for a DBS outcome.

Accurate autocorrelation modeling substantially improves fMRI reliability

Given the recent controversies in some neuroimaging statistical methods, we compare the most frequently used functional Magnetic Resonance Imaging (fMRI) analysis packages: AFNI, FSL and SPM, with regard to temporal autocorrelation modeling. This process, sometimes known as pre-whitening, is conducted in virtually all task fMRI studies. Here, we employ eleven datasets containing 980 scans corresponding to different fMRI protocols and subject populations. We found that autocorrelation modeling in AFNI, although imperfect, performed much better than the autocorrelation modeling of FSL and SPM. The presence of residual autocorrelated noise in FSL and SPM leads to heavily confounded first level results, particularly for low-frequency experimental designs. SPM’s alternative pre-whitening method, FAST, performed better than SPM’s default. The reliability of task fMRI studies could be improved with more accurate autocorrelation modeling. We recommend that fMRI analysis packages provide diagnostic plots to make users aware of any pre-whitening problems.

There has been recent controversy over the validity of commonly-used software packages for functional MRI (fMRI) data analysis. Here, the authors compare the performance of three leading packages (AFNI, FSL, SPM) in terms of temporal autocorrelation modeling, a key statistical step in fMRI analysis.

Fast Adaptive Smoothing and Thresholding for Improved Activation Detection in Low-Signal fMRI

Functional magnetic resonance imaging is a noninvasive tool for studying cerebral function. Many factors challenge activation detection, especially in low-signal scenarios that arise in the performance of high-level cognitive tasks. We provide a fully automated fast adaptive smoothing and thresholding (FAST) algorithm that uses smoothing and extreme value theory on correlated statistical parametric maps for thresholding. Performance on experiments spanning a range of low-signal settings is very encouraging. The methodology also performs well in a study to identify the cerebral regions that perceive only-auditory-reliable or only-visual-reliable speech stimuli.

Functional Connectivity of the Striatum in Schizophrenia and Psychotic Bipolar Disorder

BACKGROUND

The striatum is abnormal in schizophrenia and possibly represents a common neurobiological mechanism underlying psychotic disorders. Resting-state functional magnetic resonance imaging studies have not reached a consensus regarding striatal dysconnectivity in schizophrenia, although these studies generally find impaired frontoparietal and salience network connectivity. The goal of the current study was to clarify the pattern of corticostriatal connectivity, including whether corticostriatal dysconnectivity is transdiagnostic and extends into psychotic bipolar disorder.

METHODS

We examined corticostriatal functional connectivity in 60 healthy subjects and 117 individuals with psychosis, including 77 with a schizophrenia spectrum illness and 40 with psychotic bipolar disorder. We conducted a cortical seed-based region-of-interest analysis with follow-up voxelwise analysis for any significant results. Further, a striatum seed-based analysis was conducted to examine group differences in connectivity between the striatum and the whole cortex.

RESULTS

Cortical region-of-interest analysis indicated that overall connectivity of the salience network with the striatum was reduced in psychotic disorders, which follow-up voxelwise analysis localized to the left putamen. Striatum seed-based analyses showed reduced ventral rostral putamen connectivity with the salience network portion of the medial prefrontal cortex in both schizophrenia and psychotic bipolar disorder.

CONCLUSIONS

The current study found evidence of transdiagnostic corticostriatal dysconnectivity in both schizophrenia and psychotic bipolar disorder, including reduced salience network connectivity, as well as reduced connectivity between the putamen and the medial prefrontal cortex. Overall, the current study points to the relative importance of salience network hypoconnectivity in psychotic disorders.

Track density imaging in progressive supranuclear palsy: A pilot study

Progressive supranuclear palsy (PSP) is a neurodegenerative disorder characterized by white matter (WM) changes in different supra- and infratentorial brain structures. We used track density imaging (TDI) to characterize WM microstructural alterations in patients with PSP-Richardson's Syndrome (PSP-RS). Moreover, we investigated the diagnostic utility of TDI in distinguishing patients with PSP-RS from those with Parkinson's disease and healthy controls (HC). Twenty PSP-RS patients, 21 PD patients, and 23 HC underwent a 3 T MRI diffusion-weighted (DW) imaging. Then, we combined constrained spherical deconvolution and WM probabilistic tractography to reconstruct track density maps by calculating the number of WM streamlines traversing each voxel. Voxel-wise analysis was performed to assess group differences in track density maps. A support vector machine (SVM) approach was also used to evaluate the performance of TDI for discriminating between groups. Relative to PD patients, decreases in track density in PSP-RS patients were found in brainstem, cerebellum, thalamus, corpus callosum, and corticospinal tract. Similar findings were obtained between PSP-RS patients and HC. No differences in TDI were observed between PD and HC. SVM approach based on whole-brain analysis differentiated PD patients from PSP-RS with an area under the curve (AUC) of 0.82. The AUC reached a value of 0.98 considering only the voxels belonging to the superior cerebellar peduncle. This study shows that TDI may represent a useful approach for characterizing WM alterations in PSP-RS patients. Moreover, track density decrease in PSP could be considered a new feature for the differentiation of patients with PSP-RS from those with PD.

MACS - a new SPM toolbox for model assessment, comparison and selection

BACKGROUND

In cognitive neuroscience, functional magnetic resonance imaging (fMRI) data are widely analyzed using general linear models (GLMs). However, model quality of GLMs for fMRI is rarely assessed, in part due to the lack of formal measures for statistical model inference.

NEW METHOD

We introduce a new SPM toolbox for model assessment, comparison and selection (MACS) of GLMs applied to fMRI data. MACS includes classical, information-theoretic and Bayesian methods of model assessment previously applied to GLMs for fMRI as well as recent methodological developments of model selection and model averaging in fMRI data analysis.

RESULTS

The toolbox - which is freely available from GitHub - directly builds on the Statistical Parametric Mapping (SPM) software package and is easy-to-use, general-purpose, modular, readable and extendable. We validate the toolbox by reproducing model selection and model averaging results from earlier publications.

COMPARISON WITH EXISTING METHODS

A previous toolbox for model diagnosis in fMRI has been discontinued and other approaches to model comparison between GLMs have not been translated into reusable computational resources in the past.

CONCLUSIONS

Increased attention on model quality will lead to lower false-positive rates in cognitive neuroscience and increased application of the MACS toolbox will increase the reproducibility of GLM analyses and is likely to increase the replicability of fMRI studies.

Assessing assumptions of multivariate linear regression framework implemented for directionality analysis of fMRI

Directionality analysis of time-series, recorded from task-activated regions-of-interest (ROIs) during functional Magnetic Resonance Imaging (fMRI), has helped in gaining insights of complex human behavior and human brain functioning. The most widely used standard method of Granger Causality for evaluating directionality employ linear regression modeling of temporal processes. Such a parameter-driven approach rests on various underlying assumptions about the data. The short-comings can arise when misleading conclusions are reached after exploration of data for which the assumptions are getting violated. In this study, we assess assumptions of Multivariate Autoregressive (MAR) framework which is employed for evaluating directionality among fMRI time-series recorded during a Sensory-Motor (SM) task. The fMRI time-series here is an averaged time-series from a user-defined ROI of multiple voxels. The "aim" is to establish a step-by-step procedure using statistical methods in conjunction with graphical methods to seek the validity of MAR models, specifically in the context of directionality analysis of fMRI data which has not been done previously to the best of our knowledge. Here, in our case of SM task (block design paradigm) there is violation of assumptions, indicating the inadequacy of MAR models to find directional interactions among different task-activated regions of brain.

Improving mass-univariate analysis of neuroimaging data by modelling important unknown covariates: Application to Epigenome-Wide Association Studies

Statistical inference on neuroimaging data is often conducted using a mass-univariate model, equivalent to fitting a linear model at every voxel with a known set of covariates. Due to the large number of linear models, it is challenging to check if the selection of covariates is appropriate and to modify this selection adequately. The use of standard diagnostics, such as residual plotting, is clearly not practical for neuroimaging data. However, the selection of covariates is crucial for linear regression to ensure valid statistical inference. In particular, the mean model of regression needs to be reasonably well specified. Unfortunately, this issue is often overlooked in the field of neuroimaging. This study aims to adopt the existing Confounder Adjusted Testing and Estimation (CATE) approach and to extend it for use with neuroimaging data. We propose a modification of CATE that can yield valid statistical inferences using Principal Component Analysis (PCA) estimators instead of Maximum Likelihood (ML) estimators. We then propose a non-parametric hypothesis testing procedure that can improve upon parametric testing. Monte Carlo simulations show that the modification of CATE allows for more accurate modelling of neuroimaging data and can in turn yield a better control of False Positive Rate (FPR) and Family-Wise Error Rate (FWER). We demonstrate its application to an Epigenome-Wide Association Study (EWAS) on neonatal brain imaging and umbilical cord DNA methylation data obtained as part of a longitudinal cohort study. Software for this CATE study is freely available at http://www.bioeng.nus.edu.sg/cfa/Imaging_Genetics2.html.

A Bayesian heteroscedastic GLM with application to fMRI data with motion spikes

We propose a voxel-wise general linear model with autoregressive noise and heteroscedastic noise innovations (GLMH) for analyzing functional magnetic resonance imaging (fMRI) data. The model is analyzed from a Bayesian perspective and has the benefit of automatically down-weighting time points close to motion spikes in a data-driven manner. We develop a highly efficient Markov Chain Monte Carlo (MCMC) algorithm that allows for Bayesian variable selection among the regressors to model both the mean (i.e., the design matrix) and variance. This makes it possible to include a broad range of explanatory variables in both the mean and variance (e.g., time trends, activation stimuli, head motion parameters and their temporal derivatives), and to compute the posterior probability of inclusion from the MCMC output. Variable selection is also applied to the lags in the autoregressive noise process, making it possible to infer the lag order from the data simultaneously with all other model parameters. We use both simulated data and real fMRI data from OpenfMRI to illustrate the importance of proper modeling of heteroscedasticity in fMRI data analysis. Our results show that the GLMH tends to detect more brain activity, compared to its homoscedastic counterpart, by allowing the variance to change over time depending on the degree of head motion.

The Basolateral Amygdalae and Frontotemporal Network Functions for Threat Perception

Although the amygdalae play a central role in threat perception and reactions, the direct contributions of the amygdalae to specific aspects of threat perception, from ambiguity resolution to reflexive or deliberate action, remain ill understood in humans. Animal studies show that a detailed understanding requires a focus on the different subnuclei, which is not yet achieved in human research. Given the limits of human imaging methods, the crucial contribution needs to come from individuals with exclusive and selective amygdalae lesions. The current study investigated the role of the basolateral amygdalae and their connection with associated frontal and temporal networks in the automatic perception of threat. Functional activation and connectivity of five individuals with Urbach–Wiethe disease with focal basolateral amygdalae damage and 12 matched controls were measured with functional MRI while they attended to the facial expression of a threatening face–body compound stimuli. Basolateral amygdalae damage was associated with decreased activation in the temporal pole but increased activity in the ventral and dorsal medial prefrontal and medial orbitofrontal cortex. This dissociation between the prefrontal and temporal networks was also present in the connectivity maps. Our results contribute to a dynamic, multirole, subnuclei-based perspective on the involvement of the amygdalae in fear perception. Damage to the basolateral amygdalae decreases activity in the temporal network while increasing activity in the frontal network, thereby potentially triggering a switch from resolving ambiguity to dysfunctional threat signaling and regulation, resulting in hypersensitivity to threat.

Resolving the neural dynamics of visual and auditory scene processing in the human brain: a methodological approach

In natural environments, visual and auditory stimulation elicit responses across a large set of brain regions in a fraction of a second, yielding representations of the multimodal scene and its properties. The rapid and complex neural dynamics underlying visual and auditory information processing pose major challenges to human cognitive neuroscience. Brain signals measured non-invasively are inherently noisy, the format of neural representations is unknown, and transformations between representations are complex and often nonlinear. Further, no single non-invasive brain measurement technique provides a spatio-temporally integrated view. In this opinion piece, we argue that progress can be made by a concerted effort based on three pillars of recent methodological development: (i) sensitive analysis techniques such as decoding and cross-classification, (ii) complex computational modelling using models such as deep neural networks, and (iii) integration across imaging methods (magnetoencephalography/electroencephalography, functional magnetic resonance imaging) and models, e.g. using representational similarity analysis. We showcase two recent efforts that have been undertaken in this spirit and provide novel results about visual and auditory scene analysis. Finally, we discuss the limits of this perspective and sketch a concrete roadmap for future research.

This article is part of the themed issue ‘Auditory and visual scene analysis’.

The role of the basolateral amygdala in the perception of faces in natural contexts

The amygdala is a complex structure that plays its role in perception and threat-related behaviour by activity of its specific nuclei and their separate networks. In the present functional magnetic resonance imaging study, we investigated the role of the basolateral amygdala in face and context processing. Five individuals with focal basolateral amygdala damage and 12 matched controls viewed fearful or neutral faces in a threatening or neutral context. We tested the hypothesis that basolateral amygdala damage modifies the relation between face and threatening context, triggering threat-related activation in the dorsal stream. The findings supported this hypothesis. First, activation was increased in the right precentral gyrus for threatening versus neutral scenes in the basolateral amygdala damage group compared with the control group. Second, activity in the bilateral middle frontal gyrus, and left anterior inferior parietal lobule was enhanced for neutral faces presented in a threatening versus neutral scene in the group with basolateral amygdala damage compared with controls. These findings provide the first evidence for the neural consequences of basolateral amygdala damage during the processing of complex emotional situations.

Support vector machine classification of arterial volume-weighted arterial spin tagging images

INTRODUCTION

In recent years, machine-learning techniques have gained growing popularity in medical image analysis. Temporal brain-state classification is one of the major applications of machine-learning techniques in functional magnetic resonance imaging (fMRI) brain data. This article explores the use of support vector machine (SVM) classification technique with motor-visual activation paradigm to perform brain-state classification into activation and rest with an emphasis on different acquisition techniques.

METHODS

Images were acquired using a recently developed variant of traditional pseudocontinuous arterial spin labeling technique called arterial volume-weighted arterial spin tagging (AVAST). The classification scheme is also performed on images acquired using blood oxygenation-level dependent (BOLD) and traditional perfusion-weighted arterial spin labeling (ASL) techniques for comparison.

RESULTS

The AVAST technique outperforms traditional pseudocontinuous ASL, achieving classification accuracy comparable to that of BOLD contrast images.

CONCLUSION

This study demonstrates that AVAST has superior signal-to-noise ratio and improved temporal resolution as compared with traditional perfusion-weighted ASL and reduced sensitivity to scanner drift as compared with BOLD. Owing to these characteristics, AVAST lends itself as an ideal choice for dynamic fMRI and real-time neurofeedback experiments with sustained activation periods.

Increased fMRI Sensitivity at Equal Data Burden Using Averaged Shifted Echo Acquisition

There is growing evidence as to the benefits of collecting BOLD fMRI data with increased sampling rates. However, many of the newly developed acquisition techniques developed to collect BOLD data with ultra-short TRs require hardware, software, and non-standard analytic pipelines that may not be accessible to all researchers. We propose to incorporate the method of shifted echo into a standard multi-slice, gradient echo EPI sequence to achieve a higher sampling rate with a TR of <1 s with acceptable spatial resolution. We further propose to incorporate temporal averaging of consecutively acquired EPI volumes to both ameliorate the reduced temporal signal-to-noise inherent in ultra-fast EPI sequences and reduce the data burden. BOLD data were collected from 11 healthy subjects performing a simple, event-related visual-motor task with four different EPI sequences: (1) reference EPI sequence with TR = 1440 ms, (2) shifted echo EPI sequence with TR = 700 ms, (3) shifted echo EPI sequence with every two consecutively acquired EPI volumes averaged and effective TR = 1400 ms, and (4) shifted echo EPI sequence with every four consecutively acquired EPI volumes averaged and effective TR = 2800 ms. Both the temporally averaged sequences exhibited increased temporal signal-to-noise over the shifted echo EPI sequence. The shifted echo sequence with every two EPI volumes averaged also had significantly increased BOLD signal change compared with the other three sequences, while the shifted echo sequence with every four EPI volumes averaged had significantly decreased BOLD signal change compared with the other three sequences. The results indicated that incorporating the method of shifted echo into a standard multi-slice EPI sequence is a viable method for achieving increased sampling rate for collecting event-related BOLD data. Further, consecutively averaging every two consecutively acquired EPI volumes significantly increased the measured BOLD signal change and the subsequently calculated activation map statistics.

Bootstrapping fMRI Data: Dealing with Misspecification

The validity of inference based on the General Linear Model (GLM) for the analysis of functional magnetic resonance imaging (fMRI) time series has recently been questioned. Bootstrap procedures that partially avoid modeling assumptions may offer a welcome solution. We empirically compare two voxelwise GLM-based bootstrap approaches: a semi-parametric approach, relying solely on a model for the expected signal; and a fully parametric bootstrap approach, requiring an additional parameterization of the temporal structure. While the fully parametric approach assumes independent whitened residuals, the semi-parametric approach relies on independent blocks of residuals. The evaluation is based on inferential properties and the potential to reproduce important data characteristics. Different noise structures and data-generating mechanisms for the signal are simulated. When the model for the noise and expected signal is correct, we find that the fully parametric approach works well, with respect to both inference and reproduction of data characteristics. However, in the presence of misspecification, the fully parametric approach can be improved with additional blocking. The semi-parametric approach performs worse than the (fully) parametric approach with respect to inference but achieves comparable results as the parametric approach with additional blocking with respect to image reproducibility. We demonstrate that when the expected signal is incorrect GLM-based bootstrapping can overcome the poor performance of classical (non-bootstrap) parametric inference. We illustrate both approaches on a study exploring the neural representation of object representation in the visual pathway.

Evaluation of Second-Level Inference in fMRI Analysis

We investigate the impact of decisions in the second-level (i.e., over subjects) inferential process in functional magnetic resonance imaging on (1) the balance between false positives and false negatives and on (2) the data-analytical stability, both proxies for the reproducibility of results. Second-level analysis based on a mass univariate approach typically consists of 3 phases. First, one proceeds via a general linear model for a test image that consists of pooled information from different subjects. We evaluate models that take into account first-level (within-subjects) variability and models that do not take into account this variability. Second, one proceeds via inference based on parametrical assumptions or via permutation-based inference. Third, we evaluate 3 commonly used procedures to address the multiple testing problem: familywise error rate correction, False Discovery Rate (FDR) correction, and a two-step procedure with minimal cluster size. Based on a simulation study and real data we find that the two-step procedure with minimal cluster size results in most stable results, followed by the familywise error rate correction. The FDR results in most variable results, for both permutation-based inference and parametrical inference. Modeling the subject-specific variability yields a better balance between false positives and false negatives when using parametric inference.

Permutation and parametric tests for effect sizes in voxel-based morphometry of gray matter volume in brain structural MRI

Permutation testing has been widely implemented in voxel-based morphometry (VBM) tools. However, this type of non-parametric inference has yet to be thoroughly compared with traditional parametric inference in VBM studies of brain structure. Here we compare both types of inference and investigate what influence the number of permutations in permutation testing has on results in an exemplar study of how gray matter proportion changes with age in a group of working age adults. High resolution T1-weighted volume scans were acquired from 80 healthy adults aged 25-64years. Using a validated VBM procedure and voxel-based permutation testing for Pearson product-moment coefficient, the effect sizes of changes in gray matter proportion with age were assessed using traditional parametric and permutation testing inference with 100, 500, 1000, 5000, 10000 and 20000 permutations. The statistical significance was set at P<0.05 and false discovery rate (FDR) was used to correct for multiple comparisons. Clusters of voxels with statistically significant (PFDR<0.05) declines in gray matter proportion with age identified with permutation testing inference (N≈6000) were approximately twice the size of those identified with parametric inference (N=3221voxels). Permutation testing with 10000 (N=6251voxels) and 20000 (N=6233voxels) permutations produced clusters that were generally consistent with each other. However, with 1000 permutations there were approximately 20% more statistically significant voxels (N=7117voxels) than with ≥10000 permutations. Permutation testing inference may provide a more sensitive method than traditional parametric inference for identifying age-related differences in gray matter proportion. Based on the results reported here, at least 10000 permutations should be used in future univariate VBM studies investigating age related changes in gray matter to avoid potential false findings. Additional studies using permutation testing in large imaging databanks are required to address the impact of model complexity, multivariate analysis, number of observations, sampling bias and data quality on the accuracy with which subtle differences in brain structure associated with normal aging can be identified.

Regional covariance of muscarinic acetylcholine receptors in Alzheimer's disease using (R, R) [(123)I]-QNB SPECT

Alzheimer's disease (AD) is characterised by deficits in cholinergic neurotransmission and subsequent receptor changes. We investigated (123)I-iodo-quinuclidinyl-benzilate (QNB) SPECT images using spatial covariance analysis (SCA), to detect an M1/M4 receptor spatial covariance pattern (SCP) that distinguished AD from controls. Furthermore, a corresponding regional cerebral blood flow (rCBF) SCP was also derived. Thirty-nine subjects (15 AD and 24 healthy elderly controls) underwent (123)I-QNB and (99m)Tc-exametazime SPECT. Voxel SCA was simultaneously applied to the set of smoothed/registered scans, generating a series of eigenimages representing common intercorrelated voxels across subjects. Linear regression identified individual M1/M4 and rCBF SCPs that discriminated AD from controls. The M1/M4 SCP showed concomitant decreased uptake in medial temporal, inferior frontal, basal forebrain and cingulate relative to concomitant increased uptake in frontal poles, occipital, pre-post central and precuneus/superior parietal regions (F1,37 = 85.7, p < 0.001). A largely different perfusion SCP was obtained showing concomitant decreased rCBF in medial and superior temporal, precuneus, inferior parietal and cingulate relative to concomitant increased rCBF in cerebellum, pre-post central, putamen, fusiform and brain stem/midbrain regions (F1,37 = 77.5, p < 0.001). The M1/M4 SCP expression correlated with the duration of cognitive symptoms (r = 0.90, p < 0.001), whereas the rCBF SCP expression negatively correlated with MMSE, CAMCOG and CAMCOGmemory (r ≥ |0.63|, p ≤ 0.006). (123)I-QNB SPECT revealed an M1/M4 basocortical covariance pattern, distinct from rCBF, reflecting the duration of disease rather than current clinical symptoms. This approach could be more sensitive than univariate methods in characterising the cholinergic/rCBF changes that underpin the clinical phenotype of AD.

Use of brain MRI atlases to determine boundaries of age-related pathology: the importance of statistical method

Introduction

Neurodegenerative disease diagnoses may be supported by the comparison of an individual patient’s brain magnetic resonance image (MRI) with a voxel-based atlas of normal brain MRI. Most current brain MRI atlases are of young to middle-aged adults and parametric, e.g., mean ±standard deviation (SD); these atlases require data to be Gaussian. Brain MRI data, e.g., grey matter (GM) proportion images, from normal older subjects are apparently not Gaussian. We created a nonparametric and a parametric atlas of the normal limits of GM proportions in older subjects and compared their classifications of GM proportions in Alzheimer’s disease (AD) patients.

Methods

Using publicly available brain MRI from 138 normal subjects and 138 subjects diagnosed with AD (all 55–90 years), we created: a mean ±SD atlas to estimate parametrically the percentile ranks and limits of normal ageing GM; and, separately, a nonparametric, rank order-based GM atlas from the same normal ageing subjects. GM images from AD patients were then classified with respect to each atlas to determine the effect statistical distributions had on classifications of proportions of GM in AD patients.

Results

The parametric atlas often defined the lower normal limit of the proportion of GM to be negative (which does not make sense physiologically as the lowest possible proportion is zero). Because of this, for approximately half of the AD subjects, 25–45% of voxels were classified as normal when compared to the parametric atlas; but were classified as abnormal when compared to the nonparametric atlas. These voxels were mainly concentrated in the frontal and occipital lobes.

Discussion

To our knowledge, we have presented the first nonparametric brain MRI atlas. In conditions where there is increasing variability in brain structure, such as in old age, nonparametric brain MRI atlases may represent the limits of normal brain structure more accurately than parametric approaches. Therefore, we conclude that the statistical method used for construction of brain MRI atlases should be selected taking into account the population and aim under study. Parametric methods are generally robust for defining central tendencies, e.g., means, of brain structure. Nonparametric methods are advisable when studying the limits of brain structure in ageing and neurodegenerative disease.

Robust regression for large-scale neuroimaging studies

Multi-subject datasets used in neuroimaging group studies have a complex structure, as they exhibit non-stationary statistical properties across regions and display various artifacts. While studies with small sample sizes can rarely be shown to deviate from standard hypotheses (such as the normality of the residuals) due to the poor sensitivity of normality tests with low degrees of freedom, large-scale studies (e.g. >100 subjects) exhibit more obvious deviations from these hypotheses and call for more refined models for statistical inference. Here, we demonstrate the benefits of robust regression as a tool for analyzing large neuroimaging cohorts. First, we use an analytic test based on robust parameter estimates; based on simulations, this procedure is shown to provide an accurate statistical control without resorting to permutations. Second, we show that robust regression yields more detections than standard algorithms using as an example an imaging genetics study with 392 subjects. Third, we show that robust regression can avoid false positives in a large-scale analysis of brain-behavior relationships with over 1500 subjects. Finally we embed robust regression in the Randomized Parcellation Based Inference (RPBI) method and demonstrate that this combination further improves the sensitivity of tests carried out across the whole brain. Altogether, our results show that robust procedures provide important advantages in large-scale neuroimaging group studies.

Evaluation of statistical inference on empirical resting state fMRI

Modern statistical inference techniques may be able to improve the sensitivity and specificity of resting state functional magnetic resonance imaging (rs-fMRI) connectivity analysis through more realistic assumptions. In simulation, the advantages of such methods are readily demonstrable. However, quantitative empirical validation remains elusive in vivo as the true connectivity patterns are unknown and noise distributions are challenging to characterize, especially in ultra-high field (e.g., 7T fMRI). Though the physiological characteristics of the fMRI signal are difficult to replicate in controlled phantom studies, it is critical that the performance of statistical techniques be evaluated. The SIMulation EXtrapolation (SIMEX) method has enabled estimation of bias with asymptotically consistent estimators on empirical finite sample data by adding simulated noise . To avoid the requirement of accurate estimation of noise structure, the proposed quantitative evaluation approach leverages the theoretical core of SIMEX to study the properties of inference methods in the face of diminishing data (in contrast to increasing noise). The performance of ordinary and robust inference methods in simulation and empirical rs-fMRI are compared using the proposed quantitative evaluation approach. This study provides a simple, but powerful method for comparing a proxy for inference accuracy using empirical data.

Model specification and the reliability of fMRI results: implications for longitudinal neuroimaging studies in psychiatry

Functional Magnetic Resonance Imagine (fMRI) is an important assessment tool in longitudinal studies of mental illness and its treatment. Understanding the psychometric properties of fMRI-based metrics, and the factors that influence them, will be critical for properly interpreting the results of these efforts. The current study examined whether the choice among alternative model specifications affects estimates of test-retest reliability in key emotion processing regions across a 6-month interval. Subjects (N = 46) performed an emotional-faces paradigm during fMRI in which neutral faces dynamically morphed into one of four emotional faces. Median voxelwise intraclass correlation coefficients (mvICCs) were calculated to examine stability over time in regions showing task-related activity as well as in bilateral amygdala. Four modeling choices were evaluated: a default model that used the canonical hemodynamic response function (HRF), a flexible HRF model that included additional basis functions, a modified CompCor (mCompCor) model that added corrections for physiological noise in the global signal, and a final model that combined the flexible HRF and mCompCor models. Model residuals were examined to determine the degree to which each pipeline met modeling assumptions. Results indicated that the choice of modeling approaches impacts both the degree to which model assumptions are met and estimates of test-retest reliability. ICC estimates in the visual cortex increased from poor (mvICC = 0.31) in the default pipeline to fair (mvICC = 0.45) in the full alternative pipeline – an increase of 45%. In nearly all tests, the models with the fewest assumption violations generated the highest ICC estimates. Implications for longitudinal treatment studies that utilize fMRI are discussed.

Mapping track density changes in nigrostriatal and extranigral pathways in Parkinson's disease

In Parkinson's disease (PD) the demonstration of neuropathological disturbances in nigrostriatal and extranigral brain pathways using magnetic resonance imaging remains a challenge. Here, we applied a novel diffusion-weighted imaging approach-track density imaging (TDI). Twenty-seven non-demented Parkinson's patients (mean disease duration: 5 years, mean score on the Hoehn & Yahr scale=1.5) were compared with 26 elderly controls matched for age, sex, and education level. Track density images were created by sampling each subject's spatially normalized fiber tracks in 1mm isotropic intervals and counting the fibers that passed through each voxel. Whole-brain voxel-based analysis was performed and significance was assessed with permutation testing. Statistically significant increases in track density were found in the Parkinson's patients, relative to controls. Clusters were distributed in disease-relevant areas including motor, cognitive, and limbic networks. From the lower medulla to the diencephalon and striatum, clusters encompassed the known location of the locus coeruleus and pedunculopontine nucleus in the pons, and from the substantia nigra up to medial aspects of the posterior putamen, bilaterally. The results identified in brainstem and nigrostriatal pathways show a large overlap with the known distribution of neuropathological changes in non-demented PD patients. Our results also support an early involvement of limbic and cognitive networks in Parkinson's disease.

REGULARIZED 3D FUNCTIONAL REGRESSION FOR BRAIN IMAGE DATA VIA HAAR WAVELETS

The primary motivation and application in this article come from brain imaging studies on cognitive impairment in elderly subjects with brain disorders. We propose a regularized Haar wavelet-based approach for the analysis of three-dimensional brain image data in the framework of functional data analysis, which automatically takes into account the spatial information among neighboring voxels. We conduct extensive simulation studies to evaluate the prediction performance of the proposed approach and its ability to identify related regions to the outcome of interest, with the underlying assumption that only few relatively small subregions are truly predictive of the outcome of interest. We then apply the proposed approach to searching for brain subregions that are associated with cognition using PET images of patients with Alzheimer's disease, patients with mild cognitive impairment, and normal controls.

Masking level differences--a diffusion tensor imaging and functional MRI study

In our previous study we investigated Masking Level Differences (MLD) using functional Magnetic Resonance Imaging (fMRI), but were unable to confirm neural correlations for the MLD within the auditory cortex and inferior colliculus. Here we have duplicated conditions from our previous study, but have included more participants and changed the study site to a new location with a newer scanner and presentation system. Additionally, Diffusion Tensor Imaging (DTI) is included to allow investigation of fiber tracts that may be involved with MLDs. Twenty participants were included and underwent audiometric testing and MRI scanning. The current study revealed regions of increased and decreased activity within the auditory cortex when comparing the combined noise and signal of the dichotic MLD stimuli (N0Sπ and NπS0) with N0S0. Furthermore, we found evidence of inferior colliculus involvement. Our DTI findings show strong correlations between DTI measures within the brainstem and signal detection threshold levels. Patterns of correlation when the signal was presented only to the right ear showed an extensive network in the left hemisphere; however, the opposite was not true for the signal presented only to the left ear. Our current study was able to confirm what we had previously hypothesized using fMRI, while extending our investigation of MLDs to include the characteristics of connecting neural pathways.

Bayesian spatial transformation models with applications in neuroimaging data

The aim of this article is to develop a class of spatial transformation models (STM) to spatially model the varying association between imaging measures in a three-dimensional (3D) volume (or 2D surface) and a set of covariates. The proposed STM include a varying Box-Cox transformation model for dealing with the issue of non-Gaussian distributed imaging data and a Gaussian Markov random field model for incorporating spatial smoothness of the imaging data. Posterior computation proceeds via an efficient Markov chain Monte Carlo algorithm. Simulations and real data analysis demonstrate that the STM significantly outperforms the voxel-wise linear model with Gaussian noise in recovering meaningful geometric patterns. Our STM is able to reveal important brain regions with morphological changes in children with attention deficit hyperactivity disorder.

A quality control method for detecting and suppressing uncorrected residual motion in fMRI studies

Motion correction is an important step in the functional magnetic resonance imaging (fMRI) analysis pipeline. While many studies simply exclude subjects who are estimated to have moved beyond an arbitrary threshold, there exists no objective method for determining an appropriate threshold. Furthermore, any criterion based only upon motion estimation ignores the potential for proper realignment. The method proposed here uses unsupervised learning (specifically k-means clustering) on features derived from the mean square derivative (MSD) of the signal before and after realignment to identify problem data. These classifications are refined through analysis of correlation between subject activation maps and the mean activation map, as well as the relationship between tasking and motion as measured through regression of the canonical hemodynamic response functions to fit both estimated motion parameters and MSD. The MSD is further used to identify specific scans containing residual motion, data which is suppressed by adding nuisance regressors to the general linear model; this statistical suppression is performed for identified problem subjects, but has potential for use over all subjects. For problem subjects, our results show increased hemodynamic activity more consistent with group results; that is, the addition of nuisance regressors resulted in a doubling of the correlation between the activation map for the problem subjects and the activation map for all subjects. The proposed method should be useful in helping fMRI researchers make more efficient use of their data by reducing the need to exclude entire subjects from studies and thus collect new data to replace excluded subjects.

Functional anatomy of the masking level difference, an fMRI study

INTRODUCTION

Masking level differences (MLDs) are differences in the hearing threshold for the detection of a signal presented in a noise background, where either the phase of the signal or noise is reversed between ears. We use N0/Nπ to denote noise presented in-phase/out-of-phase between ears and S0/Sπ to denote a 500 Hz sine wave signal as in/out-of-phase. Signal detection level for the noise/signal combinations N0Sπ and NπS0 is typically 10-20 dB better than for N0S0. All combinations have the same spectrum, level, and duration of both the signal and the noise.

METHODS

Ten participants (5 female), age: 22-43, with N0Sπ-N0S0 MLDs greater than 10 dB, were imaged using a sparse BOLD fMRI sequence, with a 9 second gap (1 second quiet preceding stimuli). Band-pass (400-600 Hz) noise and an enveloped signal (.25 second tone burst, 50% duty-cycle) were used to create the stimuli. Brain maps of statistically significant regions were formed from a second-level analysis using SPM5.

RESULTS

The contrast NπS0- N0Sπ had significant regions of activation in the right pulvinar, corpus callosum, and insula bilaterally. The left inferior frontal gyrus had significant activation for contrasts N0Sπ-N0S0 and NπS0-N0S0. The contrast N0S0-N0Sπ revealed a region in the right insula, and the contrast N0S0-NπS0 had a region of significance in the left insula.

CONCLUSION

Our results extend the view that the thalamus acts as a gating mechanism to enable dichotic listening, and suggest that MLD processing is accomplished through thalamic communication with the insula, which communicate across the corpus callosum to either enhance or diminish the binaural signal (depending on the MLD condition). The audibility improvement of the signal with both MLD conditions is likely reflected by activation in the left inferior frontal gyrus, a late stage in the what/where model of auditory processing.

Identifying FMRI model violations with Lagrange multiplier tests

The standard modeling framework in functional magnetic resonance imaging (fMRI) is predicated on assumptions of linearity, time invariance and stationarity. These assumptions are rarely checked because doing so requires specialized software, although failure to do so can lead to bias and mistaken inference. Identifying model violations is an essential but largely neglected step in standard fMRI data analysis. Using Lagrange multiplier testing methods we have developed simple and efficient procedures for detecting model violations such as nonlinearity, nonstationarity and validity of the common double gamma specification for hemodynamic response. These procedures are computationally cheap and can easily be added to a conventional analysis. The test statistic is calculated at each voxel and displayed as a spatial anomaly map which shows regions where a model is violated. The methodology is illustrated with a large number of real data examples.

Accounting for movement increases sensitivity in detecting brain activity in Parkinson's disease

Parkinson's disease (PD) is manifested by motor impairment, which may impede the ability to accurately perform motor tasks during functional magnetic resonance imaging (fMRI). Both temporal and amplitude deviations of movement performance affect the blood oxygenation level-dependent (BOLD) response. We present a general approach for assessing PD patients' movement control employing simultaneously recorded fMRI time series and behavioral data of the patients' kinematics using MR-compatible gloves. Twelve male patients with advanced PD were examined with fMRI at 1.5T during epoch-based visually paced finger tapping. MR-compatible gloves were utilized online to quantify motor outcome in two conditions with or without dopaminergic medication. Modeling of individual-level brain activity included (i) a predictor consisting of a condition-specific, constant-amplitude boxcar function convolved with the canonical hemodynamic response function (HRF) as commonly used in fMRI statistics (standard model), or (ii) a custom-made predictor computed from glove time series convolved with the HRF (kinematic model). Factorial statistics yielded a parametric map for each modeling technique, showing the medication effect on the group level. Patients showed bilateral response to levodopa in putamen and globus pallidus during the motor experiment. Interestingly, kinematic modeling produced significantly higher activation in terms of both the extent and amplitude of activity. Our results appear to account for movement performance in fMRI motor experiments with PD and increase sensitivity in detecting brain response to levodopa. We strongly advocate quantitatively controlling for motor performance to reach more reliable and robust analyses in fMRI with PD patients.

An iterative jackknife approach for assessing reliability and power of FMRI group analyses

For functional magnetic resonance imaging (fMRI) group activation maps, so-called second-level random effect approaches are commonly used, which are intended to be generalizable to the population as a whole. However, reliability of a certain activation focus as a function of group composition or group size cannot directly be deduced from such maps. This question is of particular relevance when examining smaller groups (<20–27 subjects). The approach presented here tries to address this issue by iteratively excluding each subject from a group study and presenting the overlap of the resulting (reduced) second-level maps in a group percent overlap map. This allows to judge where activation is reliable even upon excluding one, two, or three (or more) subjects, thereby also demonstrating the inherent variability that is still present in second-level analyses. Moreover, when progressively decreasing group size, foci of activation will become smaller and/or disappear; hence, the group size at which a given activation disappears can be considered to reflect the power necessary to detect this particular activation. Systematically exploiting this effect allows to rank clusters according to their observable effect size. The approach is tested using different scenarios from a recent fMRI study (children performing a “dual-use” fMRI task, n = 39), and the implications of this approach are discussed.