Accurate modeling of temporal correlations in rapidly sampled fMRI time series

Does higher sampling rate (multiband + SENSE) improve group statistics - An example from social neuroscience block design at 3T

Multiband (MB) or Simultaneous multi-slice (SMS) acquisition schemes allow the acquisition of MRI signals from more than one spatial coordinate at a time. Commercial availability has brought this technique within the reach of many neuroscientists and psychologists. Most early evaluation of the performance of MB acquisition employed resting state fMRI or the most basic tasks. In this study, we tested whether the advantages of using MB acquisition schemes generalize to group analyses using a cognitive task more representative of typical cognitive neuroscience applications. Twenty-three subjects were scanned on a Philips 3 ​T scanner using five sequences, up to eight-fold acceleration with MB-factors 1 to 4, SENSE factors up to 2 and corresponding TRs of 2.45s down to 0.63s, while they viewed (i) movie blocks showing complex actions with hand object interactions and (ii) control movie blocks without hand object interaction. Data were processed using a widely used analysis pipeline implemented in SPM12 including the unified segmentation and canonical HRF modelling. Using random effects group-level, voxel-wise analysis we found that all sequences were able to detect the basic action observation network known to be recruited by our task. The highest t-values were found for sequences with MB4 acceleration. For the MB1 sequence, a 50% bigger voxel volume was needed to reach comparable t-statistics. The group-level t-values for resting state networks (RSNs) were also highest for MB4 sequences. Here the MB1 sequence with larger voxel size did not perform comparable to the MB4 sequence. Altogether, we can thus recommend the use of MB4 (and SENSE 1.5 or 2) on a Philips scanner when aiming to perform group-level analyses using cognitive block design fMRI tasks and voxel sizes in the range of cortical thickness (e.g. 2.7 ​mm isotropic). While results will not be dramatically changed by the use of multiband, our results suggest that MB will bring a moderate but significant benefit.

Improved autoregressive model for correction of noise serial correlation in fast fMRI

PURPOSE

In rapidly acquired functional MRI (fast fMRI) data, the noise serial correlations (SC) can produce problematically overestimated T-statistics which lead to invalid statistical inferences. This study aims to evaluate and improve the accuracy of high-order autoregressive model (AR(p), where p is the model order) based prewhitening method in the SC correction.

METHODS

Fast fMRI images were acquired at rest (null data) using a multiband simultaneous multi-slice echo planar imaging pulse sequence with repetition time (TR) = 300 and 500 ms. The SC effect in the fast fMRI data was corrected using the prewhitening method based on two AR(p) models: (1) the conventional model (fixed AR(p)) which preselects a constant p for all the image voxels; (2) an improved model (AR_AICc ) that employs the corrected Akaike information criterion voxel-wise to automatically select the model orders for each voxel. To evaluate accuracy of SC correction, false positive characteristics were measured by assuming the presence of block and event-related tasks in the null data without image smoothing. The performance of prewhitening was also examined in smoothed images by adding pseudo task fMRI signals into the null data and comparing the detected to simulated activations (ground truth).

RESULTS

The measured false positive characteristics agreed well with the theoretical curve when using the AR_AICc , and the activation maps in the smoothed data matched the ground truth. The AR_AICc showed improved performance than the fixed AR(p) method.

CONCLUSION

The AR_AICc can effectively remove noise SC, and accurate statistical analysis results can be obtained with the AR_AICc correction in fast fMRI.

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.

Comparison of SMS-EPI and 3D-EPI at 7T in an fMRI localizer study with matched spatiotemporal resolution and homogenized excitation profiles

The simultaneous multi-slice EPI (SMS-EPI, a.k.a. MB-EPI) sequence has met immense popularity recently in functional neuroimaging. A still less common alternative is the use of 3D-EPI, which offers similar acceleration capabilities. The aim of this work was to compare the SMS-EPI and the 3D-EPI sequences in terms of sampling strategies for the detection of task-evoked activations at 7T using detection theory. To this end, the spatial and temporal resolutions of the sequences were matched (1.6 mm isotropic resolution, TR = 1200 ms) and their excitation profiles were homogenized by means of calibration-free parallel-transmission (Universal Pulses). We used a fast-event “localizer” paradigm of 5:20 min in order to probe sensorimotor functions (visual, auditory and motor tasks) as well as higher level functions (language comprehension, mental calculation), where results from a previous large-scale study at 3T (N = 81) served as ground-truth reference for the brain areas implicated in each cognitive function. In the current study, ten subjects were scanned while their activation maps were generated for each cognitive function with the GLM analysis. The SMS-EPI and 3D-EPI sequences were compared in terms of raw tSNR, t-score testing for the mean signal, activation strength and accuracy of the robust sensorimotor functions. To this end, the sensitivity and specificity of these contrasts were computed by comparing their activation maps to the reference brain areas obtained in the 3T study. Estimated flip angle distributions in the brain reported a normalized root mean square deviation from the target value below 10% for both sequences. The analysis of the t-score testing for the mean signal revealed temporal noise correlations, suggesting the use of this metric instead of the traditional tSNR for testing fMRI sequences. The SMS-EPI and 3D-EPI thereby yielded similar performance from a detection theory perspective.

Model-based physiological noise removal in fast fMRI

Recent improvements in the speed and sensitivity of fMRI acquisition techniques suggest that fast fMRI can be used to detect and precisely localize sub-second neural dynamics. This enhanced temporal resolution has enormous potential for neuroscientists. However, physiological noise poses a major challenge for the analysis of fast fMRI data. Physiological noise scales with sensitivity, and its autocorrelation structure is altered in rapidly sampled data, suggesting that new approaches are needed for physiological noise removal in fast fMRI. Existing strategies either rely on external physiological recordings, which can be noisy or difficult to collect, or employ data-driven approaches which make assumptions that may not hold true in fast fMRI. We created a statistical model of harmonic regression with autoregressive noise (HRAN) to estimate and remove cardiac and respiratory noise from the fMRI signal directly. This technique exploits the fact that cardiac and respiratory noise signals are fully sampled (rather than aliasing) when imaging at fast rates, allowing us to track and model physiology over time without requiring external physiological measurements. We then created a joint model of neural hemodynamics, and physiological and autocorrelated noise to more accurately remove noise. We first verified that HRAN accurately estimates cardiac and respiratory dynamics and that our model demonstrates goodness-of-fit in fast fMRI data. In task-driven data, we then demonstrated that HRAN is able to remove physiological noise while leaving the neural signal intact, thereby increasing detection of task-driven voxels. Finally, we established that in both simulations and fast fMRI data HRAN is able to improve statistical inferences as compared with gold-standard physiological noise removal techniques. In conclusion, we created a tool that harnesses the novel information in fast fMRI to remove physiological noise, enabling broader use of the technology to study human brain function.

Impact of sampling rate on statistical significance for single subject fMRI connectivity analysis

A typical time series in functional magnetic resonance imaging (fMRI) exhibits autocorrelation, that is, the samples of the time series are dependent. In addition, temporal filtering, one of the crucial steps in preprocessing of functional magnetic resonance images, induces its own autocorrelation. While performing connectivity analysis in fMRI, the impact of the autocorrelation is largely ignored. Recently, autocorrelation has been addressed by variance correction approaches, which are sensitive to the sampling rate. In this article, we aim to investigate the impact of the sampling rate on the variance correction approaches. Toward this end, we first derived a generalized expression for the variance of the sample Pearson correlation coefficient (SPCC) in terms of the sampling rate and the filter cutoff frequency, in addition to the autocorrelation and cross-covariance functions of the time series. Through simulations, we illustrated the importance of the variance correction for a fixed sampling rate. Using the real resting state fMRI data sets, we demonstrated that the data sets with higher sampling rates were more prone to false positives, in agreement with the existing empirical reports. We further demonstrated with single subject results that for the data sets with higher sampling rates, the variance correction strategy restored the integrity of true connectivity.

The potential of MR-Encephalography for BCI/Neurofeedback applications with high temporal resolution

Real-time functional magnetic resonance imaging (rt-fMRI) enables the update of various brain-activity measures during an ongoing experiment as soon as a new brain volume is acquired. However, the recorded Blood-oxygen-level dependent (BOLD) signal also contains physiological artifacts such as breathing and heartbeat, which potentially cause misleading false positive effects especially problematic in brain-computer interface (BCI) and neurofeedback (NF) setups. The low temporal resolution of echo planar imaging (EPI) sequences (which is in the range of seconds) prevents a proper separation of these artifacts from the BOLD signal. MR-Encephalography (MREG) has been shown to provide the high temporal resolution required to unalias and correct for physiological fluctuations and leads to increased specificity and sensitivity for mapping task-based activation and functional connectivity as well as for detecting dynamic changes in connectivity over time. By comparing a simultaneous multislice echo planar imaging (SMS-EPI) sequence and an MREG sequence using the same nominal spatial resolution in an offline analysis for three different experimental fMRI paradigms (perception of house and face stimuli, motor imagery, Stroop task), the potential of this novel technique for future BCI and NF applications was investigated. First, adapted general linear model pre-whitening which accounts for the high temporal resolution in MREG was implemented to calculate proper statistical results and be able to compare these with the SMS-EPI sequence. Furthermore, the respiration- and cardiac pulsation-related signals were successfully separated from the MREG signal using independent component analysis which were then included as regressors for a GLM analysis. Only the MREG sequence allowed to clearly separate cardiac pulsation and respiration components from the signal time course. It could be shown that these components highly correlate with the recorded respiration and cardiac pulsation signals using a respiratory belt and fingertip pulse plethysmograph. Temporal signal-to-noise ratios of SMS-EPI and MREG were comparable. Functional connectivity analysis using partial correlation showed a reduced standard error in MREG compared to SMS-EPI. Also, direct time course comparisons by down-sampling the MREG signal to the SMS-EPI temporal resolution showed lower variance in MREG. In general, we show that the higher temporal resolution is beneficial for fMRI time course modeling and this aspect can be exploited in offline application but also, is especially attractive, for real-time BCI and NF applications.

Brain networks underlying tactile softness perception: A functional magnetic resonance imaging study

Humans are adept at perceiving physical properties of an object through touch. Tangible object properties can be categorized into two types: macro-spatial properties, including shape and orientation; and material properties, such as roughness, softness, and temperature. Previous neuroimaging studies have shown that roughness and temperature are extracted at nodes of a network, such as that involving the parietal operculum and insula, which is different from the network engaged in processing macro-spatial properties. However, it is unclear whether other perceptual dimensions pertaining to material properties engage the same regions. Here, we conducted a functional magnetic resonance imaging study to test whether the parietal operculum and insula were involved in extracting tactually-perceived softness magnitude. Fifty-six healthy right-handed participants estimated perceived softness magnitude using their right middle finger. We presented three stimuli that had the same shape but different compliances. The force applied to the finger was manipulated at two levels. Classical mass-univariate analysis showed that activity in the parietal operculum, insula, and medial prefrontal cortex was positively associated with perceived softness magnitude, regardless of the applied force. Softness-related activity was stronger in the ventral striatum in the high-force condition than in the low-force condition. The multivariate voxel pattern analysis showed higher accuracy than chance levels and control regions in the parietal operculum/insula, postcentral gyrus, posterior parietal lobule, and middle occipital gyrus. These results indicate that a distributed set of the brain regions, including the parietal operculum and insula, is involved in representing perceived softness.

On the analysis of rapidly sampled fMRI data

Recent advances in parallel imaging and simultaneous multi-slice techniques have permitted whole-brain fMRI acquisitions at sub-second sampling intervals, without significantly sacrificing the spatial coverage and resolution. Apart from probing brain function at finer temporal scales, faster sampling rates may potentially lead to enhanced functional sensitivity, owing possibly to both cleaner neural representations (due to less aliased physiological noise) and additional statistical benefits (due to more degrees of freedom for a fixed scan duration). Accompanying these intriguing aspects of fast acquisitions, however, confusion has also arisen regarding (1) how to preprocess/analyze these fast fMRI data, and (2) what exactly is the extent of benefits with fast acquisitions, i.e., how fast is fast enough for a specific research aim? The first question is motivated by the altered spectral distribution and noise characteristics at short sampling intervals, while the second question seeks to reconcile the complicated trade-offs between the functional contrast-to-noise ratio and the effective degrees of freedom. Although there have been recent efforts to empirically approach different aspects of these two questions, in this work we discuss, from a theoretical perspective accompanied by some illustrative, proof-of-concept experimental in vivo human fMRI data, a few considerations that are rarely mentioned, yet are important for both preprocessing and optimizing statistical inferences for studies that employ acquisitions with sub-second sampling intervals. Several summary recommendations include concerns regarding advisability of relying on low-pass filtering to de-noise physiological contributions, employment of statistical models with sufficient complexity to account for the substantially increased serial correlation, and cautions regarding using rapid sampling to enhance functional sensitivity given that different analysis models may associate with distinct trade-offs between contrast-to-noise ratios and the effective degrees of freedom. As an example, we demonstrate that as TR shortens, the intrinsic differences in how noise is accommodated in general linear models and Pearson correlation analyses (assuming Gaussian distributed stochastic signals and noise) can result in quite different outcomes, either gaining or losing statistical power.

Intracortical smoothing of small-voxel fMRI data can provide increased detection power without spatial resolution losses compared to conventional large-voxel fMRI data

Continued improvement in MRI acquisition technology has made functional MRI (fMRI) with small isotropic voxel sizes down to 1 mm and below more commonly available. Although many conventional fMRI studies seek to investigate regional patterns of cortical activation for which conventional voxel sizes of 3 mm and larger provide sufficient spatial resolution, smaller voxels can help avoid contamination from adjacent white matter (WM) and cerebrospinal fluid (CSF), and thereby increase the specificity of fMRI to signal changes within the gray matter. Unfortunately, temporal signal-to-noise ratio (tSNR), a metric of fMRI sensitivity, is reduced in high-resolution acquisitions, which offsets the benefits of small voxels. Here we introduce a framework that combines small, isotropic fMRI voxels acquired at 7 T field strength with a novel anatomically-informed, surface mesh-navigated spatial smoothing that can provide both higher detection power and higher resolution than conventional voxel sizes. Our smoothing approach uses a family of intracortical surface meshes and allows for kernels of various shapes and sizes, including curved 3D kernels that adapt to and track the cortical folding pattern. Our goal is to restrict smoothing to the cortical gray matter ribbon and avoid noise contamination from CSF and signal dilution from WM via partial volume effects. We found that the intracortical kernel that maximizes tSNR does not maximize percent signal change (ΔS/S), and therefore the kernel configuration that optimizes detection power cannot be determined from tSNR considerations alone. However, several kernel configurations provided a favorable balance between boosting tSNR and ΔS/S, and allowed a 1.1-mm isotropic fMRI acquisition to have higher performance after smoothing (in terms of both detection power and spatial resolution) compared to an unsmoothed 3.0-mm isotropic fMRI acquisition. Overall, the results of this study support the strategy of acquiring voxels smaller than the cortical thickness, even for studies not requiring high spatial resolution, and smoothing them down within the cortical ribbon with a kernel of an appropriate shape to achieve the best performance-thus decoupling the choice of fMRI voxel size from the spatial resolution requirements of the particular study. The improvement of this new intracortical smoothing approach over conventional surface-based smoothing is expected to be modest for conventional resolutions, however the improvement is expected to increase with higher resolutions. This framework can also be applied to anatomically-informed intracortical smoothing of higher-resolution data (e.g. along columns and layers) in studies with prior information about the spatial structure of activation.

Optimizing Data for Modeling Neuronal Responses

In this technical note, we address an unresolved challenge in neuroimaging statistics: how to determine which of several datasets is the best for inferring neuronal responses. Comparisons of this kind are important for experimenters when choosing an imaging protocol-and for developers of new acquisition methods. However, the hypothesis that one dataset is better than another cannot be tested using conventional statistics (based on likelihood ratios), as these require the data to be the same under each hypothesis. Here we present Bayesian data comparison (BDC), a principled framework for evaluating the quality of functional imaging data, in terms of the precision with which neuronal connectivity parameters can be estimated and competing models can be disambiguated. For each of several candidate datasets, neuronal responses are modeled using Bayesian (probabilistic) forward models, such as General Linear Models (GLMs) or Dynamic Casual Models (DCMs). Next, the parameters from subject-specific models are summarized at the group level using a Bayesian GLM. A series of measures, which we introduce here, are then used to evaluate each dataset in terms of the precision of (group-level) parameter estimates and the ability of the data to distinguish similar models. To exemplify the approach, we compared four datasets that were acquired in a study evaluating multiband fMRI acquisition schemes, and we used simulations to establish the face validity of the comparison measures. To enable people to reproduce these analyses using their own data and experimental paradigms, we provide general-purpose Matlab code via the SPM software.

A circular echo planar sequence for fast volumetric fMRI

PURPOSE

To demonstrate a circular EPI (CEPI) sequence as well as a generalized EPI reconstruction for fast fMRI with parallel imaging acceleration.

METHODS

The CEPI acquisition was constructed using variable readout lengths and maximum ramp sampling as well as blipped-CAIPI z-gradient encoding for simultaneous multislice (SMS) and 3D volumetric imaging. A signal equation model with constant and linear phase terms was used to iteratively reconstruct images with low ghosting. Simulation, phantom, and human imaging experiments including audio/visual fMRI were performed at 3T using a 52-channel coil.

RESULTS

Application of CEPI gradients with duration of 27 ms covering a 22-cm FOV at a 64 × 64 pixel resolution in SMS and 3D acquisitions resulted in images with comparable quality to those of standard Cartesian EPI. With parallel imaging techniques robust detection of BOLD fMRI activation with temporal sampling down to 275 ms was possible. The high temporal resolution enabled higher activation statistics at a penalty in increased noise and residual aliasing. The un-accelerated 3D acquisition showed large temporal instability compared with a standard 2D acquisition.

CONCLUSION

Nonuniform sampling and generalized image reconstructions can be applied to EPI acquisitions including those with blipped-CAIPI z gradients. The same gradients can be used for either SMS or 3D acquisitions providing identical coverage.

Accurate modeling of temporal correlations in rapidly sampled fMRI time series

Rapid imaging techniques are increasingly used in functional MRI studies because they allow a greater number of samples to be acquired per unit time, thereby increasing statistical power. However, temporal correlations limit the increase in functional sensitivity and must be accurately accounted for to control the false‐positive rate. A common approach to accounting for temporal correlations is to whiten the data prior to estimating fMRI model parameters. Models of white noise plus a first‐order autoregressive process have proven sufficient for conventional imaging studies, but more elaborate models are required for rapidly sampled data. Here we show that when the “FAST” model implemented in SPM is used with a well‐controlled number of parameters, it can successfully prewhiten 80% of grey matter voxels even with volume repetition times as short as 0.35 s. We further show that the temporal signal‐to‐noise ratio (tSNR), which has conventionally been used to assess the relative functional sensitivity of competing imaging approaches, can be augmented to account for the temporal correlations in the time series. This amounts to computing the t‐score testing for the mean signal. We show in a visual perception task that unlike the tSNR weighted by the number of samples, the t‐score measure is directly related to the t‐score testing for activation when the temporal correlations are correctly modeled. This score affords a more accurate means of evaluating the functional sensitivity of different data acquisition options.