Analysis of activity in fMRI data using affinity propagation clustering

Clustering analysis is a promising data-driven method for the analysis of functional magnetic resonance imaging (fMRI) data. The huge computation load, however, makes it difficult for the practical use. We use affinity propagation clustering (APC), a new clustering algorithm especially for large data sets to detect brain functional activation from fMRI. It considers all data points as possible exemplars through the minimisation of an energy function and message-passing architecture, and obtains the optimal set of exemplars and their corresponding clusters. Four simulation studies and three in vivo fMRI studies reveal that brain functional activation can be effectively detected and that different response patterns can be distinguished using this method. Our results demonstrate that APC is superior to the k-centres clustering, as revealed by their performance measures in the weighted Jaccard coefficient and average squared error. These results suggest that the proposed APC will be useful in detecting brain functional activation from fMRI data.     

Granger mediation analysis of multiple time series with an application to functional magnetic resonance imaging

This paper presents Granger mediation analysis, a new framework for causal mediation analysis of multiple time series. This framework is motivated by a functional magnetic resonance imaging (fMRI) experiment where we are interested in estimating the mediation effects between a randomized stimulus time series and brain activity time series from two brain regions. The independent observation assumption is thus unrealistic for this type of time‐series data. To address this challenge, our framework integrates two types of models: causal mediation analysis across the mediation variables, and vector autoregressive (VAR) models across the temporal observations. We use “Granger” to refer to VAR correlations modeled in this paper. We further extend this framework to handle multilevel data, in order to model individual variability and correlated errors between the mediator and the outcome variables. Using Rubin's potential outcome framework, we show that the causal mediation effects are identifiable under our time‐series model. We further develop computationally efficient algorithms to maximize our likelihood‐based estimation criteria. Simulation studies show that our method reduces the estimation bias and improves statistical power, compared with existing approaches. On a real fMRI data set, our approach quantifies the causal effects through a brain pathway, while capturing the dynamic dependence between two brain regions.     

Assessment of brain interactivity in the motor cortex from the concept of functional connectivity and spectral analysis of fMRI data

Functional magnetic resonance imaging (fMRI) was used to assess the contributions of movement preparation and execution of a visuomotor task in a cerebral motor network. The functional connectivity of the voxel time series between brain regions in the frequency space was investigated by performing spectral analysis of fMRI time series. The regional interactivities between the two portions of the supplementary motor area (pre-SMA and SMA-proper) and the primary motor cortex (M1), defined as a seed region, were evaluated. The spectral parameter of coherence was used to describe a correlation structure in the frequency domain between two voxel-based time series and to infer the strength of the functional interaction within our presumed motor network of connections. The results showed meaningful differences of the functional interactions between the two portions of the SMA and the M1 area depending on the task conditions. This approach demonstrated the existence of a functional dissociation between the pre-SMA and SMA-proper subregions. We therefore conclude that spectral analysis is useful for identifying functional interactions of brain regions and might provide a powerful tool to quantify changes in connectivity profiles associated with various components of an experimental task.     

Multiplexed echo planar imaging for sub-second whole brain FMRI and fast diffusion imaging

Echo planar imaging (EPI) is an MRI technique of particular value to neuroscience, with its use for virtually all functional MRI (fMRI) and diffusion imaging of fiber connections in the human brain. EPI generates a single 2D image in a fraction of a second; however, it requires 2-3 seconds to acquire multi-slice whole brain coverage for fMRI and even longer for diffusion imaging. Here we report on a large reduction in EPI whole brain scan time at 3 and 7 Tesla, without significantly sacrificing spatial resolution, and while gaining functional sensitivity. The multiplexed-EPI (M-EPI) pulse sequence combines two forms of multiplexing: temporal multiplexing (m) utilizing simultaneous echo refocused (SIR) EPI and spatial multiplexing (n) with multibanded RF pulses (MB) to achieve m×n images in an EPI echo train instead of the normal single image. This resulted in an unprecedented reduction in EPI scan time for whole brain fMRI performed at 3 Tesla, permitting TRs of 400 ms and 800 ms compared to a more conventional 2.5 sec TR, and 2-4 times reductions in scan time for HARDI imaging of neuronal fibertracks. The simultaneous SE refocusing of SIR imaging at 7 Tesla advantageously reduced SAR by using fewer RF refocusing pulses and by shifting fat signal out of the image plane so that fat suppression pulses were not required. In preliminary studies of resting state functional networks identified through independent component analysis, the 6-fold higher sampling rate increased the peak functional sensitivity by 60%. The novel M-EPI pulse sequence resulted in a significantly increased temporal resolution for whole brain fMRI, and as such, this new methodology can be used for studying non-stationarity in networks and generally for expanding and enriching the functional information.     

Analysis of activity in fMRI data using affinity propagation clustering

Clustering analysis is a promising data-driven method for the analysis of functional magnetic resonance imaging (fMRI) data. The huge computation load, however, makes it difficult for the practical use. We use affinity propagation clustering (APC), a new clustering algorithm especially for large data sets to detect brain functional activation from fMRI. It considers all data points as possible exemplars through the minimisation of an energy function and message-passing architecture, and obtains the optimal set of exemplars and their corresponding clusters. Four simulation studies and three in vivo fMRI studies reveal that brain functional activation can be effectively detected and that different response patterns can be distinguished using this method. Our results demonstrate that APC is superior to the k-centres clustering, as revealed by their performance measures in the weighted Jaccard coefficient and average squared error. These results suggest that the proposed APC will be useful in detecting brain functional activation from fMRI data.     

Identifying Neural Drivers with Functional MRI: An Electrophysiological Validation

Whether functional magnetic resonance imaging (fMRI) allows the identification of neural drivers remains an open question of particular importance to refine physiological and neuropsychological models of the brain, and/or to understand neurophysiopathology. Here, in a rat model of absence epilepsy showing spontaneous spike-and-wave discharges originating from the first somatosensory cortex (S1BF), we performed simultaneous electroencephalographic (EEG) and fMRI measurements, and subsequent intracerebral EEG (iEEG) recordings in regions strongly activated in fMRI (S1BF, thalamus, and striatum). fMRI connectivity was determined from fMRI time series directly and from hidden state variables using a measure of Granger causality and Dynamic Causal Modelling that relates synaptic activity to fMRI. fMRI connectivity was compared to directed functional coupling estimated from iEEG using asymmetry in generalised synchronisation metrics. The neural driver of spike-and-wave discharges was estimated in S1BF from iEEG, and from fMRI only when hemodynamic effects were explicitly removed. Functional connectivity analysis applied directly on fMRI signals failed because hemodynamics varied between regions, rendering temporal precedence irrelevant. This paper provides the first experimental substantiation of the theoretical possibility to improve interregional coupling estimation from hidden neural states of fMRI. As such, it has important implications for future studies on brain connectivity using functional neuroimaging.     

Double‐wavelet transform for multisubject task‐induced functional magnetic resonance imaging data

The goal of this article is to model multisubject task‐induced functional magnetic resonance imaging (fMRI) response among predefined regions of interest (ROIs) of the human brain. Conventional approaches to fMRI analysis only take into account temporal correlations, but do not rigorously model the underlying spatial correlation due to the complexity of estimating and inverting the high dimensional spatio‐temporal covariance matrix. Other spatio‐temporal model approaches estimate the covariance matrix with the assumption of stationary time series, which is not always feasible. To address these limitations, we propose a double‐wavelet approach for modeling the spatio‐temporal brain process. Working with wavelet coefficients simplifies temporal and spatial covariance structure because under regularity conditions, wavelet coefficients are approximately uncorrelated. Different wavelet functions were used to capture different correlation structures in the spatio‐temporal model. The main advantages of the wavelet approach are that it is scalable and that it deals with nonstationarity in brain signals. Simulation studies showed that our method could reduce false‐positive and false‐negative rates by taking into account spatial and temporal correlations simultaneously. We also applied our method to fMRI data to study activation in prespecified ROIs in the prefontal cortex. Data analysis showed that the result using the double‐wavelet approach was more consistent than the conventional approach when sample size decreased.     

Topographic Factor Analysis: A Bayesian Model for Inferring Brain Networks from Neural Data

The neural patterns recorded during a neuroscientific experiment reflect complex interactions between many brain regions, each comprising millions of neurons. However, the measurements themselves are typically abstracted from that underlying structure. For example, functional magnetic resonance imaging (fMRI) datasets comprise a time series of three-dimensional images, where each voxel in an image (roughly) reflects the activity of the brain structure(s)–located at the corresponding point in space–at the time the image was collected. FMRI data often exhibit strong spatial correlations, whereby nearby voxels behave similarly over time as the underlying brain structure modulates its activity. Here we develop topographic factor analysis (TFA), a technique that exploits spatial correlations in fMRI data to recover the underlying structure that the images reflect. Specifically, TFA casts each brain image as a weighted sum of spatial functions. The parameters of those spatial functions, which may be learned by applying TFA to an fMRI dataset, reveal the locations and sizes of the brain structures activated while the data were collected, as well as the interactions between those structures.     

Self-regulation of local brain activity using real-time functional magnetic resonance imaging (fMRI)

Functional magnetic resonance imaging (fMRI) measures the blood oxygen level-dependent (BOLD) signal related to neuronal activity. So far, this technique has been limited by time-consuming data analysis impeding on-line analysis. In particular, no brain-computer interface (BCI) was available which provided on-line feedback to learn physiological self-regulation of the BOLD signal. Recently, studies have shown that fMRI feedback is feasible and facilitates voluntary control of brain activity. Here we review these studies to make the fMRI feedback methodology accessible to a broader scientific community such as researchers concerned with functional brain imaging and the neurobiology of learning. Methodological and conceptual limitations were substantially reduced by artefact control, sensitivity improvements, real-time algorithms, and adapted experimental designs. Physiological self-regulation of the local BOLD response is a new paradigm for cognitive neuroscience to study brain plasticity and the functional relevance of regulated brain areas by modification of behaviour. Voluntary control of abnormal activity in circumscribed brain areas may even be applied as psychophysiological treatment.     

Bayesian Spatiotemporal Inference in Functional Magnetic Resonance Imaging

Mapping of the human brain by means of functional magnetic resonance imaging (fMRI) is an emerging field in cognitive and clinical neuroscience. Current techniques to detect activated areas of the brain mostly proceed in two steps. First, conventional methods of correlation, regression, and time series analysis are used to assess activation by a separate, pixelwise comparison of the fMRI signal time courses to the reference function of a presented stimulus. Spatial aspects caused by correlations between neighboring pixels are considered in a separate second step, if at all. The aim of this article is to present hierarchical Bayesian approaches that allow one to simultaneously incorporate temporal and spatial dependencies between pixels directly in the model formulation. For reasons of computational feasibility, models have to be comparatively parsimonious, without oversimplifying. We introduce parametric and semiparametric spatial and spatiotemporal models that proved appropriate and illustrate their performance applied to visual fMRI data.     

FMRI Effective Connectivity and TMS Chronometry: Complementary Accounts of Causality in the Visuospatial Judgment Network

Background

While traditionally quite distinct, functional neuroimaging (e.g. functional magnetic resonance imaging: fMRI) and functional interference techniques (e.g. transcranial magnetic stimulation: TMS) increasingly address similar questions of functional brain organization, including connectivity, interactions, and causality in the brain. Time-resolved TMS over multiple brain network nodes can elucidate the relative timings of functional relevance for behavior (“TMS chronometry”), while fMRI functional or effective connectivity (fMRI EC) can map task-specific interactions between brain regions based on the interrelation of measured signals. The current study empirically assessed the relation between these different methods.

Methodology/Principal Findings

One group of 15 participants took part in two experiments: one fMRI EC study, and one TMS chronometry study, both of which used an established cognitive paradigm involving one visuospatial judgment task and one color judgment control task. Granger causality mapping (GCM), a data-driven variant of fMRI EC analysis, revealed a frontal-to-parietal flow of information, from inferior/middle frontal gyrus (MFG) to posterior parietal cortex (PPC). FMRI EC-guided Neuronavigated TMS had behavioral effects when applied to both PPC and to MFG, but the temporal pattern of these effects was similar for both stimulation sites. At first glance, this would seem in contradiction to the fMRI EC results. However, we discuss how TMS chronometry and fMRI EC are conceptually different and show how they can be complementary and mutually constraining, rather than contradictory, on the basis of our data.

Conclusions/Significance

The findings that fMRI EC could successfully localize functionally relevant TMS target regions on the single subject level, and conversely, that TMS confirmed an fMRI EC identified functional network to be behaviorally relevant, have important methodological and theoretical implications. Our results, in combination with data from earlier studies by our group (Sack et al., 2007, Cerebral Cortex), lead to informed speculations on complex brain mechanisms, and TMS disruption thereof, underlying visuospatial judgment. This first in-depth empirical and conceptual comparison of fMRI EC and TMS chronometry thereby shows the complementary insights offered by the two methods.     

Using Informational Connectivity to Measure the Synchronous Emergence of fMRI Multi-voxel Information Across Time

It is now appreciated that condition-relevant information can be present within distributed patterns of functional magnetic resonance imaging (fMRI) brain activity, even for conditions with similar levels of univariate activation. Multi-voxel pattern (MVP) analysis has been used to decode this information with great success. FMRI investigators also often seek to understand how brain regions interact in interconnected networks, and use functional connectivity (FC) to identify regions that have correlated responses over time. Just as univariate analyses can be insensitive to information in MVPs, FC may not fully characterize the brain networks that process conditions with characteristic MVP signatures. The method described here, informational connectivity (IC), can identify regions with correlated changes in MVP-discriminability across time, revealing connectivity that is not accessible to FC. The method can be exploratory, using searchlights to identify seed-connected areas, or planned, between pre-selected regions-of-interest. The results can elucidate networks of regions that process MVP-related conditions, can breakdown MVPA searchlight maps into separate networks, or can be compared across tasks and patient groups.     

Asynchronous functional linear regression models for longitudinal data in reproducing kernel Hilbert space

Motivated by the analysis of longitudinal neuroimaging studies, we study the longitudinal functional linear regression model under asynchronous data setting for modeling the association between clinical outcomes and functional (or imaging) covariates. In the asynchronous data setting, both covariates and responses may be measured at irregular and mismatched time points, posing methodological challenges to existing statistical methods. We develop a kernel weighted loss function with roughness penalty to obtain the functional estimator and derive its representer theorem. The rate of convergence, a Bahadur representation, and the asymptotic pointwise distribution of the functional estimator are obtained under the reproducing kernel Hilbert space framework. We propose a penalized likelihood ratio test to test the nullity of the functional coefficient, derive its asymptotic distribution under the null hypothesis, and investigate the separation rate under the alternative hypotheses. Simulation studies are conducted to examine the finite-sample performance of the proposed procedure. We apply the proposed methods to the analysis of multitype data obtained from the Alzheimer's Disease Neuroimaging Initiative (ADNI) study, which reveals significant association between 21 regional brain volume density curves and the cognitive function. Data used in preparation of this paper were obtained from the ADNI database (adni.loni.usc.edu).     

Functional imaging of memory processes in humans: positron emission tomography and functional magnetic resonance imaging

Human memory is not a unitary function; it consists of multiple memory systems, with different characteristics and specialisations that are implemented in the brain. The cognitive neuroscience of human memory tries to comprehend how we encode, store, and retrieve memory items within and across those systems. The emergence of functional neuroimaging techniques offered the unprecedented opportunity to directly observe the brain regions engaged in memory functions. Brain imaging techniques can roughly be divided into those measuring the electric or magnetic fields generated by neuronal activity (EEG, magnetencephalography [MEG]) and those measuring the haemodynamic or metabolic sequelae of neuronal activity (positron emission tomography [PET], functional magnetic resonance imaging [fMRI]). Out of these techniques, the following two will be discussed in detail: fMRI and PET. Although functional neuroimaging is able to acquire images of the brain engaged in consolidating or retrieving memories, these processes are not clearly visible in the data. Statistical techniques are needed to reduce the complexity of the data and to extract the processes of interest. This article outlines the experimental and analytical procedures of neuroimaging studies with PET and fMRI. We will use a PET-study on episodic memory in human volunteers to illustrate design, analysis, and interpretation of functional imaging studies on memory.     

fMRI Validation of fNIRS Measurements During a Naturalistic Task

We present a method to compare brain activity recorded with near-infrared spectroscopy (fNIRS) in a dance video game task to that recorded in a reduced version of the task using fMRI (functional magnetic resonance imaging). Recently, it has been shown that fNIRS can accurately record functional brain activities equivalent to those concurrently recorded with functional magnetic resonance imaging for classic psychophysical tasks and simple finger tapping paradigms. However, an often quoted benefit of fNIRS is that the technique allows for studying neural mechanisms of complex, naturalistic behaviors that are not possible using the constrained environment of fMRI. Our goal was to extend the findings of previous studies that have shown high correlation between concurrently recorded fNIRS and fMRI signals to compare neural recordings obtained in fMRI procedures to those separately obtained in naturalistic fNIRS experiments. Specifically, we developed a modified version of the dance video game Dance Dance Revolution (DDR) to be compatible with both fMRI and fNIRS imaging procedures. In this methodology we explain the modifications to the software and hardware for compatibility with each technique as well as the scanning and calibration procedures used to obtain representative results. The results of the study show a task-related increase in oxyhemoglobin in both modalities and demonstrate that it is possible to replicate the findings of fMRI using fNIRS in a naturalistic task. This technique represents a methodology to compare fMRI imaging paradigms which utilize a reduced-world environment to fNIRS in closer approximation to naturalistic, full-body activities and behaviors. Further development of this technique may apply to neurodegenerative diseases, such as Parkinson’s disease, late states of dementia, or those with magnetic susceptibility which are contraindicated for fMRI scanning.     

Interleaved EPI Based fMRI Improved by Multiplexed Sensitivity Encoding (MUSE) and Simultaneous Multi-Band Imaging

Functional magnetic resonance imaging (fMRI) is a non-invasive and powerful imaging tool for detecting brain activities. The majority of fMRI studies are performed with single-shot echo-planar imaging (EPI) due to its high temporal resolution. Recent studies have demonstrated that, by increasing the spatial-resolution of fMRI, previously unidentified neuronal networks can be measured. However, it is challenging to improve the spatial resolution of conventional single-shot EPI based fMRI. Although multi-shot interleaved EPI is superior to single-shot EPI in terms of the improved spatial-resolution, reduced geometric distortions, and sharper point spread function (PSF), interleaved EPI based fMRI has two main limitations: 1) the imaging throughput is lower in interleaved EPI; 2) the magnitude and phase signal variations among EPI segments (due to physiological noise, subject motion, and B₀ drift) are translated to significant in-plane aliasing artifact across the field of view (FOV). Here we report a method that integrates multiple approaches to address the technical limitations of interleaved EPI-based fMRI. Firstly, the multiplexed sensitivity-encoding (MUSE) post-processing algorithm is used to suppress in-plane aliasing artifacts resulting from time-domain signal instabilities during dynamic scans. Secondly, a simultaneous multi-band interleaved EPI pulse sequence, with a controlled aliasing scheme incorporated, is implemented to increase the imaging throughput. Thirdly, the MUSE algorithm is then generalized to accommodate fMRI data obtained with our multi-band interleaved EPI pulse sequence, suppressing both in-plane and through-plane aliasing artifacts. The blood-oxygenation-level-dependent (BOLD) signal detectability and the scan throughput can be significantly improved for interleaved EPI-based fMRI. Our human fMRI data obtained from 3 Tesla systems demonstrate the effectiveness of the developed methods. It is expected that future fMRI studies requiring high spatial-resolvability and fidelity will largely benefit from the reported techniques.     

The neural basis of the blood-oxygen-level-dependent functional magnetic resonance imaging signal

Magnetic resonance imaging (MRI) has rapidly become an important tool in clinical medicine and biological research. Its functional variant (functional magnetic resonance imaging; fMRI) is currently the most widely used method for brain mapping and studying the neural basis of human cognition. While the method is widespread, there is insufficient knowledge of the physiological basis of the fMRI signal to interpret the data confidently with respect to neural activity. This paper reviews the basic principles of MRI and fMRI, and subsequently discusses in some detail the relationship between the blood-oxygen-level-dependent (BOLD) fMRI signal and the neural activity elicited during sensory stimulation. To examine this relationship, we conducted the first simultaneous intracortical recordings of neural signals and BOLD responses. Depending on the temporal characteristics of the stimulus, a moderate to strong correlation was found between the neural activity measured with microelectrodes and the BOLD signal averaged over a small area around the microelectrode tips. However, the BOLD signal had significantly higher variability than the neural activity, indicating that human fMRI combined with traditional statistical methods underestimates the reliability of the neuronal activity. To understand the relative contribution of several types of neuronal signals to the haemodynamic response, we compared local field potentials (LFPs), single- and multi-unit activity (MUA) with high spatio-temporal fMRI responses recorded simultaneously in monkey visual cortex. At recording sites characterized by transient responses, only the LFP signal was significantly correlated with the haemodynamic response. Furthermore, the LFPs had the largest magnitude signal and linear systems analysis showed that the LFPs were better than the MUAs at predicting the fMRI responses. These findings, together with an analysis of the neural signals, indicate that the BOLD signal primarily measures the input and processing of neuronal information within a region and not the output signal transmitted to other brain regions.     

Modulatory Interactions of Resting-State Brain Functional Connectivity

The functional brain connectivity studies are generally based on the synchronization of the resting-state functional magnetic resonance imaging (fMRI) signals. Functional connectivity measures usually assume a stable relationship over time; however, accumulating studies have reported time-varying properties of strength and spatial distribution of functional connectivity. The present study explored the modulation of functional connectivity between two regions by a third region using the physiophysiological interaction (PPI) technique. We first identified eight brain networks and two regions of interest (ROIs) representing each of the networks using a spatial independent component analysis. A voxel-wise analysis was conducted to identify regions that showed modulatory interactions (PPI) with the two ROIs of each network. Mostly, positive modulatory interactions were observed within regions involved in the same system. For example, the two regions of the dorsal attention network revealed modulatory interactions with the regions related to attention, while the two regions of the extrastriate network revealed modulatory interactions with the regions in the visual cortex. In contrast, the two regions of the default mode network (DMN) revealed negative modulatory interactions with the regions in the executive network, and vice versa, suggesting that the activities of one network may be associated with smaller within network connectivity of the competing network. These results validate the use of PPI analysis to study modulation of resting-state functional connectivity by a third region. The modulatory effects may provide a better understanding of complex brain functions.