Beyond the “Pain Matrix,” inter-run synchronization during mechanical nociceptive stimulation

PurposeRecent studies have applied the new magnetic resonance encephalography (MREG) sequence to the study of interictal epileptic discharges (IEDs) in the electroencephalogram (EEG) of epileptic patients. However, there are no criteria to quantitatively evaluate different processing methods, to properly use the new sequence.MethodsWe evaluated different processing steps of this new sequence under the common generalized linear model (GLM) framework by assessing the reliability of results. A bootstrap sampling technique was first used to generate multiple replicated data sets; a GLM with different processing steps was then applied to obtain activation maps, and the reliability of these maps was assessed.ResultsWe applied our analysis in an event‐related GLM related to IEDs. A higher reliability was achieved by using a GLM with head motion confound regressor with 24 components rather than the usual 6, with an autoregressive model of order 5 and with a canonical hemodynamic response function (HRF) rather than variable latency or patient‐specific HRFs. Comparison of activation with IED field also favored the canonical HRF, consistent with the reliability analysis.ConclusionThe reliability analysis helps to optimize the processing methods for this fast fMRI sequence, in a context in which we do not know the ground truth of activation areas. Magn Reson Med 78:1–1, 2017. © 2016 International Society for Magnetic Resonance in Medicine

BOLD fMRI has become a prevalent method to study cerebral sensory processing in rodent disease models, including pain and mechanical hypersensitivity. fMRI data analysis is frequently combined with a general-linear-model (GLM) -based analysis, which uses the convolution of a hemodynamic response function (HRF) with the stimulus paradigm. However, several studies indicated that the HRF differs across species, sexes, brain structures, and experimental factors, including stimulation modalities or anesthesia, and hence might strongly affect the outcome of BOLD analyzes. While considerable work has been done in humans and rats to understand the HRF, much less is known in mice. As a prerequisite to investigate mechano-sensory processing and BOLD fMRI data in male and female mice, we (1) designed a rotating stimulator that allows application of two different mechanical modalities, including innocuous von Frey and noxious pinprick stimuli and (2) determined and statistically compared HRFs across 30 brain structures and experimental conditions, including sex and, stimulus modalities. We found that mechanical stimulation lead to brain-wide BOLD signal changes thereby allowing extraction of HRFs from multiple brain structures. However, we did not find differences in HRFs across all brain structures and experimental conditions. Hence, we computed a whole-brain mouse HRF, which is based on 88 functional scans from 30 mice. A comparison of this mouse-specific HRF with our previously reported rat-derived HRF showed significantly slower kinetics in mice. Finally, we detected pronounced differences in cerebral BOLD activation between male and female mice with mechanical stimulation, thereby exposing divergent processing of noxious and innocuous stimuli in both sexes.

In an event-related functional MRI data analysis, an accurate and robust extraction of the hemodynamic response function (HRF) and its associated statistics (e.g., magnitude, width, and time to peak) is critical to infer quantitative information about the relative timing of the neuronal events in different brain regions. The aim of this paper is to develop a multiscale adaptive smoothing model (MASM) to accurately estimate HRFs pertaining to each stimulus sequence across all voxels. MASM explicitly accounts for both spatial and temporal smoothness information, while incorporating such information to adaptively estimate HRFs in the frequency domain. One simulation study and a real data set are used to demonstrate the methodology and examine its finite sample performance in HRF estimation, which confirms that MASM significantly outperforms the existing methods including the smooth finite impulse response model, the inverse logit model and the canonical HRF.

The hemodynamic response function (HRF) links neuronal activity to functional magnetic resonance imaging (fMRI) signals. While most fMRI studies use a canonical HRF, increasing evidence from studies of healthy subjects suggests that the HRF depends on anatomical location and disease states. Here, we investigate how HRF variability relates to anatomical location and pathology in the epileptic brain, using a large simultaneous electroencephalogram and fMRI dataset. Applying HRF deconvolution and temporal decomposition, we built the first whole-brain HRF library specific to epilepsy, identifying four distinct shape groups. We mapped HRF features across parcellations of two atlases using novel Bayesian hierarchical models. In non-epileptogenic regions, HRF shape and spatial distributions align with findings from healthy subjects. Within pathological regions, they vary significantly according to pathology. Our results indicate that HRF variability is associated with pathology, in addition to its dependence on anatomical location, motivating region- and pathology-based HRF modulation in epilepsy studies.

Using patient-specific hemodynamic response function in epileptic spike analysis of human epilepsy: a study based on EEG–fNIRS

The relation between electrophysiology and BOLD‐fMRI requires further elucidation. One approach for studying this relation is to find time‐frequency features from electrophysiology that explain the variance of BOLD time‐series. Convolution of these features with a canonical hemodynamic response function (HRF) is often required to model neurovascular coupling mechanisms and thus account for time shifts between electrophysiological and BOLD‐fMRI data. We propose a framework for extracting the spatial distribution of these time‐frequency features while also estimating more flexible, region‐specific HRFs. The core component of this method is the decomposition of a tensor containing impulse response functions using the Canonical Polyadic Decomposition. The outputs of this decomposition provide insight into the relation between electrophysiology and BOLD‐fMRI and can be used to construct estimates of BOLD time‐series. We demonstrated the performance of this method on simulated data while also examining the effects of simulated measurement noise and physiological confounds. Afterwards, we validated our method on publicly available task‐based and resting‐state EEG‐fMRI data. We adjusted our method to accommodate the multisubject nature of these datasets, enabling the investigation of inter‐subject variability with regards to EEG‐to‐BOLD neurovascular coupling mechanisms. We thus also demonstrate how EEG features for modelling the BOLD signal differ across subjects.

.SignificanceAdvances in electronics have allowed the recent development of compact, high channel count time domain functional near-infrared spectroscopy (TD-fNIRS) systems. Temporal moment analysis has been proposed for increased brain sensitivity due to the depth selectivity of higher order temporal moments. We propose a general linear model (GLM) incorporating TD moment data and auxiliary physiological measurements, such as short separation channels, to improve the recovery of the HRF.AimsWe compare the performance of previously reported multi-distance TD moment techniques to commonly used techniques for continuous wave (CW) fNIRS hemodynamic response function (HRF) recovery, namely block averaging and CW GLM. Additionally, we compare the multi-distance TD moment technique to TD moment GLM.ApproachWe augmented resting TD-fNIRS moment data (six subjects) with known synthetic HRFs. We then employed block averaging and GLM techniques with “short-separation regression” designed both for CW and TD to recover the HRFs. We calculated the root mean square error (RMSE) and the correlation of the recovered HRF to the ground truth. We compared the performance of equivalent CW and TD techniques with paired t-tests.ResultsWe found that, on average, TD moment HRF recovery improves correlations by 98% and 48% for HbO and HbR respectively, over CW GLM. The improvement on the correlation for TD GLM over TD moment is 12% (HbO) and 27% (HbR). RMSE decreases 56% and 52% (HbO and HbR) for TD moment compared to CW GLM. We found no statistically significant improvement in the RMSE for TD GLM compared to TD moment.ConclusionsProperly covariance-scaled TD moment techniques outperform their CW equivalents in both RMSE and correlation in the recovery of the synthetic HRFs. Furthermore, our proposed TD GLM based on moments outperforms regular TD moment analysis, while allowing the incorporation of auxiliary measurements of the confounding physiological signals from the scalp.

In this paper, we used time-domain functional near infrared spectroscopy (TD-fNIRS) to evaluate the haemodynamic response function (HRF) in the occipital cortex following visual stimulation in glaucomatous eyes as compared to healthy eyes. A total of 98 subjects were enrolled in the study and clinically classified as healthy subjects, glaucoma patients (primary open-angle glaucoma) and mixed subjects (i.e. with a different classification for the two eyes). After quality check data were used from HRF of 73 healthy and 62 glaucomatous eyes. The amplitudes of the oxygenated and deoxygenated haemoglobin concentrations, together with their latencies with respect to the stimulus onset, were estimated by fitting their time course with a canonical HRF. Statistical analysis showed that the amplitudes of both haemodynamic parameters show a significant association with the pathology and a significant discriminating ability, while no significant result was found for latencies. Overall, our findings together with the ease of use and noninvasiveness of TD-NIRS, make this technique a promising candidate as a supporting tool for a better evaluation of the glaucoma pathology.

The hemodynamic response function (HRF) describes the dynamic process of the blood oxygen level dependent (BOLD) signal evoked by a very short stimulus of unit intensity. Understanding HRF is a key issue for exploring the underlying dynamics of nervous system activation. Some researches suggested that the shape of HRF varied with different brain regions and different subjects, even at different time point in the same brain region of a subject. In this work, epilepsy discharges were adopted as markers of short stimuli, and then we use EEG-fMRI method to investigate HRF in temporal lobe epilepsy with interictal regional slow wave. Five patients with idiopathic temporal lobe epilepsy were enrolled in this study, the peak delay and amplitude of the estimated HRF was compared with Gamma, Glover and the SPM canonical HRF, and the results show that the delay of the estimated HRF is different from either one of the three theoretical model, thus further study on both theoretical model and real data estimation is needed.

The impact of temporal regularization on estimates of the BOLD hemodynamic response function: A comparative analysis

Evidence suggesting common mechanisms underlie contralateral and ipsilateral negative BOLD responses in the human visual cortex

The procedure to estimate brain activity based on fMR signals is a process based on many assumptions. Some of the methods such as GLM (General Linear Model) and ICA(Independent Component Analysis) used for this purpose contain several restrictions. In GLM, it is assumed that each active voxel responds similarly and linearly towards a given sitmulus. In ICA, an unsurmountable number of independent time series are produced, one of which is assumed to reflect the activity pattern. In this study, we used minimal number of assumptions to estimate an underlying HRF (hemodynamic response function) from a given fMR time series, and then used the estimated HRFs to classify voxels as active or passive. We have investigated results from simulations and real fMR experiments.

The objective of fMRI data analysis is to detect the region of the brain that gets activated in response to a specific stimulus presented to the subject. We develop a new algorithm for activation detection in event-related fMRI data. We utilize a forward model for fMRI data acquisition which explicitly incorporates physiological noise, scanner noise and the spatial blurring introduced by the scanner. After slice-by-slice image restoration procedure that independently restores each data slice corresponding to each time index, we estimate the parameters of the hemodynamic response function (HRF) model for each pixel of the restored data. In order to enforce spatial regularity in our estimates, we model the prior distribution of the HRF parameters as a generalized Gaussian Markov random field (GGMRF) model. We develop an algorithm to compute the maximum a posteriori (MAP) estimates of the parameters. We then threshold the amplitude parameters to obtain the final activation map. We illustrate our algorithm by comparing it with the widely used general linear model (GLM) method. In synthetic data experiments, under the same probability of false alarm, the probability of correct detection for our method is up to 15% higher than GLM. In real data experiments, through anatomical analysis and benchmark testing using block paradigm results, we demonstrate that our algorithm produces fewer false alarms than GLM.

The simultaneous acquisition and subsequent analysis of EEG and fMRI data is challenging owing to increased noise levels in the EEG data. A common method to integrate data from these two modalities is to use aspects of the EEG data, such as the amplitudes of event-related potentials (ERP) or oscillatory EEG activity, to predict fluctuations in the fMRI data. However, this relies on the acquisition of high quality datasets to ensure that only the correlates of neuronal activity are being studied. In this study, we investigate the effects of head-motion-related artefacts in the EEG signal on the predicted T2*-weighted signal variation. We apply our analyses to two independent datasets: 1) four participants were asked to move their feet in the scanner to generate small head movements, and 2) four participants performed an episodic memory task. We created T2*-weighted signal predictors from indicators of abrupt head motion using derivatives of the realignment parameters, from visually detected artefacts in the EEG as well as from three EEG frequency bands (theta, alpha and beta). In both datasets, we found little correlation between the T2*-weighted signal and EEG predictors that were not convolved with the canonical haemodynamic response function (cHRF). However, all convolved EEG predictors strongly correlated with the T2*-weighted signal variation in various regions including the bilateral superior temporal cortex, supplementary motor area, medial parietal cortex and cerebellum. The finding that movement onset spikes in the EEG predict T2*-weighted signal intensity only when the time course of movements is convolved with the cHRF, suggests that the correlated signal might reflect a BOLD response to neural activity associated with head movement. Furthermore, the observation that broad-spectral EEG spikes tend to occur at the same time as abrupt head movements, together with the finding that abrupt movements and EEG spikes show similar correlations with the T2*-weighted signal, indicates that the EEG spikes are produced by abrupt movement and that continuous regressors of EEG oscillations contain motion-related noise even after stringent correction of the EEG data. If not properly removed, these artefacts complicate the use of EEG data as a predictor of T2*-weighted signal variation.

Inter-subject parcellation of functional Magnetic Resonance Imaging (fMRI) data based on a standard General Linear Model (GLM) and spectral clustering was recently proposed as a means to alleviate the issues associated with spatial normalization in fMRI. However, for all its appeal, a GLM-based parcellation approach introduces its own biases, in the form of a priori knowledge about the shape of Hemodynamic Response Function (HRF) and task-related signal changes, or about the subject behaviour during the task. In this paper, we introduce a data-driven version of the spectral clustering parcellation, based on Independent Component Analysis (ICA) and Partial Least Squares (PLS) instead of the GLM. First, a number of independent components are automatically selected. Seed voxels are then obtained from the associated ICA maps and we compute the PLS latent variables between the fMRI signal of the seed voxels (which covers regional variations of the HRF) and the principal components of the signal across all voxels. Finally, we parcellate all subjects data with a spectral clustering of the PLS latent variables. We present results of the application of the proposed method on both single-subject and multi-subject fMRI datasets. Preliminary experimental results, evaluated with intra-parcel variance of GLM t-values and PLS derived t-values, indicate that this data-driven approach offers improvement in terms of parcellation accuracy over GLM based techniques.