Resting-State fMRI: Emerging Concepts for Future Clinical Application.

Resting-state functional magnetic resonance imaging (rsfMRI) has been developed as a method of investigating spontaneous neural activity. Based on its low-frequency signal synchronization, rsfMRI has made it possible to identify multiple macroscopic structures termed resting-state networks (RSNs) on a single scan of less than 10 minutes. It is easy to implement even in clinical practice, in which assigning tasks to patients can be challenging. These advantages have accelerated the adoption and growth of rsfMRI. Recently, studies on the global rsfMRI signal have attracted increasing attention. Because it primarily arises from physiological events, less attention has hitherto been paid to the global signal than to the local network (i.e., RSN) component. However, the global signal is not a mere nuisance or a subsidiary component. On the contrary, it is quantitatively the dominant component that accounts for most of the variance in the rsfMRI signal throughout the brain and provides rich information on local hemodynamics that can serve as an individual-level diagnostic biomarker. Moreover, spatiotemporal analyses of the global signal have revealed that it is closely and fundamentally associated with the organization of RSNs, thus challenging the basic assumptions made in conventional rsfMRI analyses and views on RSNs. This review introduces new concepts emerging from rsfMRI spatiotemporal analyses focusing on the global signal and discusses how they may contribute to future clinical medicine. EVIDENCE LEVEL: 5 TECHNICAL EFFICACY: Stage 1.

Neural efficiency and proficiency adaptation of effective connectivity corresponding to early and advanced skill levels in athletes of racket sports.

This study explored how the neural efficiency and proficiency worked in athletes with different skill levels from the perspective of effective connectivity brain network in resting state. The deconvolved conditioned Granger causality (GC) analysis was applied to functional magnetic resonance imaging (fMRI) data of 35 elite athletes (EAs) and 42 student-athletes (SAs) of racket sports as well as 39 normal controls (NCs), to obtain the voxel-wised hemodynamic response function (HRF) parameters representing the functional segregation and effective connectivity representing the functional integration. The results showed decreased time-to-peak of HRF in the visual attention brain regions in the two athlete groups compared with NC and decreased response height in the advanced motor control brain regions in EA comparing to the nonelite groups, suggesting the neural efficiency represented by the regional HRF was different in early and advanced skill levels. GC analysis demonstrated that the GC values within the middle occipital gyrus had a linear trend from negative to positive, suggesting a stepwise "neural proficiency" of the effective connectivity from NC to SA then to EA. The GC values of the inter-lobe circuits in EA had the trend to regress to NC levels, in agreement with the neural efficiency of these circuits in EA. Further feature selection approach suggested the important role of the cerebral-brainstem GC circuit for discriminating EA. Our findings gave new insight into the complementary neural mechanisms in brain functional segregation and integration, which was associated with early and advanced skill levels in athletes of racket sports.

Hemodynamic response function description in patients with glioma.

INTRODUCTION: Functional magnetic resonance imaging is a powerful tool that has provided many insights into cognitive sciences. Yet, as its analysis is mostly based on the knowledge of an a priori canonical hemodynamic response function (HRF), its reliability in patients' applications has been questioned. There have been reports of neurovascular uncoupling in patients with glioma, but no specific description of the Hemodynamic Response Function (HRF) in glioma has been reported so far. The aim of this work is to describe the HRF in patients with glioma. METHODS: Forty patients were included. MR images were acquired on a 1.5T scanner. Activated clusters were identified using a fuzzy general linear model; HRFs were adjusted with a double-gamma function. Analyses were undertaken considering the tumor grade, age, sex, tumor location, and activated location. RESULTS: Differences are found in the occipital, limbic, insular, and sub-lobar areas, but not in the frontal, temporal, and parietal lobes. The presence of a glioma slows the time-to-peak and onset times by 5.2 and 3.8 % respectively; high-grade gliomas present 8.1 % smaller HRF widths than low-grade gliomas. DISCUSSION AND CONCLUSION: There is significant HRF variation due to the presence of glioma, but the magnitudes of the observed differences are small. Most processing pipelines should be robust enough for this magnitude of variation and little if any impact should be visible on functional maps. The differences that have been observed in the literature between functional mapping obtained with magnetic resonance vs. that obtained with direct electrostimulation during awake surgery are more probably due to the intrinsic difference in the mapping process: fMRI mapping detects all recruited areas while intra-surgical mapping indicates only the areas indispensable for the realization of a certain task. Surgical mapping might not be the gold standard to use when trying to validate the fMRI mapping process.

Dynamic variations of resting-state BOLD signal spectra in white matter.

Recent studies have demonstrated that the mathematical model used for analyzing and interpreting fMRI data in gray matter (GM) is inappropriate for detecting or describing blood-oxygenation-level-dependent (BOLD) signals in white matter (WM). In particular the hemodynamic response function (HRF) which serves as the regressor in general linear models is different in WM compared to GM. We recently reported measurements of the frequency contents of resting-state signal time courses in WM that showed distinct power spectra which depended on local structural-vascular-functional associations. In addition, multiple studies of GM have revealed how functional connectivity between regions, as measured by the correlation between BOLD time series, varies dynamically over time. We therefore investigated whether and how BOLD signals from WM in a resting state varied over time. We measured voxel-wise spectrograms, which reflect the time-varying spectral patterns of WM time courses. The results suggest that the spectral patterns are non-stationary but could be categorized into five modes that recurred over time. These modes showed distinct spatial distributions of their occurrences and durations, and the distributions were highly consistent across individuals. In addition, one of the modes exhibited a strong coupling of its occurrence between GM and WM across individuals, and two communities of WM voxels were identified according to the hierarchical structures of transitions among modes. Moreover, these modes are coupled to the shape of instantaneous HRFs. Our findings extend previous studies and reveal the non-stationary nature of spectral patterns of BOLD signals over time, providing a spatial-temporal-frequency characterization of resting-state signals in WM.

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

The task-evoked positive BOLD response (PBR) to a unilateral visual hemi-field stimulation is often accompanied by robust and sustained contralateral as well as ipsilateral negative BOLD responses (NBRs) in the visual cortex. The signal characteristics and the neural and/or vascular mechanisms that underlie these two types of NBRs are not completely understood. In this paper, we investigated the properties of these two types of NBRs. We first demonstrated the linearity of both NBRs with respect to stimulus duration. Next, we showed that the hemodynamic response functions (HRFs) of the two NBRs were similar to each other, but significantly different from that of the PBR. Moreover, the subject-wise expressions of the two NBRs were tightly coupled to the degree that the correlation between the two NBRs was significantly higher than the correlation between each NBR and the PBR. However, the activation patterns of the two NBRs did not show a high level of interhemispheric spatial similarity, and the functional connectivity between them was not different than the interhemispheric functional connectivity between the NBRs and PBR. Finally, while attention did modulate both NBRs, the attention-related changes in their HRFs were similar. Our findings suggest that the two NBRs might be generated through common neural and/or vascular mechanisms involving distal/deep brain regions that project to the two hemispheres.

Modeling the Hemodynamic Response Function Using EEG-fMRI Data During Eyes-Open Resting-State Conditions and Motor Task Execution.

Being able to accurately quantify the hemodynamic response function (HRF) that links the blood oxygen level dependent functional magnetic resonance imaging (BOLD-fMRI) signal to the underlying neural activity is important both for elucidating neurovascular coupling mechanisms and improving the accuracy of fMRI-based functional connectivity analyses. In particular, HRF estimation using BOLD-fMRI is challenging particularly in the case of resting-state data, due to the absence of information about the underlying neuronal dynamics. To this end, using simultaneously recorded electroencephalography (EEG) and fMRI data is a promising approach, as EEG provides a more direct measure of neural activations. In the present work, we employ simultaneous EEG-fMRI to investigate the regional characteristics of the HRF using measurements acquired during resting conditions. We propose a novel methodological approach based on combining distributed EEG source space reconstruction, which improves the spatial resolution of HRF estimation and using block-structured linear and nonlinear models, which enables us to simultaneously obtain HRF estimates and the contribution of different EEG frequency bands. Our results suggest that the dynamics of the resting-state BOLD signal can be sufficiently described using linear models and that the contribution of each band is region specific. Specifically, it was found that sensory-motor cortices exhibit positive HRF shapes, whereas the lateral occipital cortex and areas in the parietal cortex, such as the inferior and superior parietal lobule exhibit negative HRF shapes. To validate the proposed method, we repeated the analysis using simultaneous EEG-fMRI measurements acquired during execution of a unimanual hand-grip task. Our results reveal significant associations between BOLD signal variations and electrophysiological power fluctuations in the ipsilateral primary motor cortex, particularly for the EEG beta band, in agreement with previous studies in the literature.

Investigating mechanisms of fast BOLD responses: The effects of stimulus intensity and of spatial heterogeneity of hemodynamics.

Recent studies have demonstrated that fast fMRI can track neural activity well above the temporal limit predicted by the canonical hemodynamic response model. While these findings are promising, the biophysical mechanisms underlying these fast fMRI phenomena remain underexplored. In this study, we discuss two aspects of the hemodynamic response, complementary to several existing hypotheses, that can accommodate faster fMRI dynamics beyond those predicted by the canonical model. First, we demonstrate, using both visual and somatosensory paradigms, that the timing and shape of hemodynamic response functions (HRFs) vary across graded levels of stimulus intensity-with lower-intensity stimulation eliciting faster and narrower HRFs. Second, we show that as the spatial resolution of fMRI increases, voxel-wise HRFs begin to deviate from the canonical model, with a considerable portion of voxels exhibiting faster temporal dynamics than predicted by the canonical HRF. Collectively, both stimulus/task intensity and image resolution can affect the sensitivity of fMRI to fast brain activity, which may partly explain recent observations of fast fMRI signals. It is further noteworthy that, while the present investigations focus on fast neural responses, our findings suggest that a revised hemodynamic model may benefit the many fMRI studies using paradigms with wide ranges of contrast levels (e.g., resting or naturalistic conditions) or with modern, high-resolution MR acquisitions.

Augmenting interictal mapping with neurovascular coupling biomarkers by structured factorization of epileptic EEG and fMRI data.

EEG-correlated fMRI analysis is widely used to detect regional BOLD fluctuations that are synchronized to interictal epileptic discharges, which can provide evidence for localizing the ictal onset zone. However, the typical, asymmetrical and mass-univariate approach cannot capture the inherent, higher order structure in the EEG data, nor multivariate relations in the fMRI data, and it is nontrivial to accurately handle varying neurovascular coupling over patients and brain regions. We aim to overcome these drawbacks in a data-driven manner by means of a novel structured matrix-tensor factorization: the single-subject EEG data (represented as a third-order spectrogram tensor) and fMRI data (represented as a spatiotemporal BOLD signal matrix) are jointly decomposed into a superposition of several sources, characterized by space-time-frequency profiles. In the shared temporal mode, Toeplitz-structured factors account for a spatially specific, neurovascular 'bridge' between the EEG and fMRI temporal fluctuations, capturing the hemodynamic response's variability over brain regions. By analyzing interictal data from twelve patients, we show that the extracted source signatures provide a sensitive localization of the ictal onset zone (10/12). Moreover, complementary parts of the IOZ can be uncovered by inspecting those regions with the most deviant neurovascular coupling, as quantified by two entropy-like metrics of the hemodynamic response function waveforms (9/12). Hence, this multivariate, multimodal factorization provides two useful sets of EEG-fMRI biomarkers, which can assist the presurgical evaluation of epilepsy. We make all code required to perform the computations available at https://github.com/svaneynd/structured-cmtf.

Simultaneous fMRI and fast-scan cyclic voltammetry bridges evoked oxygen and neurotransmitter dynamics across spatiotemporal scales.

The vascular contributions of neurotransmitters to the hemodynamic response are gaining more attention in neuroimaging studies, as many neurotransmitters are vasomodulatory. To date, well-established electrochemical techniques that detect neurotransmission in high magnetic field environments are limited. Here, we propose an experimental setting enabling simultaneous fast-scan cyclic voltammetry (FSCV) and blood oxygenation level-dependent functional magnetic imaging (BOLD fMRI) to measure both local tissue oxygen and dopamine responses, and global BOLD changes, respectively. By using MR-compatible materials and the proposed data acquisition schemes, FSCV detected physiological analyte concentrations with high temporal resolution and spatial specificity inside of a 9.4 T MRI bore. We found that tissue oxygen and BOLD correlate strongly, and brain regions that encode dopamine amplitude differences can be identified via modeling simultaneously acquired dopamine FSCV and BOLD fMRI time-courses. This technique provides complementary neurochemical and hemodynamic information and expands the scope of studying the influence of local neurotransmitter release over the entire brain.

Retinotopic variations of the negative blood-oxygen-level dependent hemodynamic response function in human primary visual cortex.

Functional magnetic resonance imaging (fMRI) measures blood-oxygen-level-dependent (BOLD) contrast that is generally assumed to be linearly related to excitatory neural activity. The positive hemodynamic response function (pHRF) is the positive BOLD response (PBR) evoked by a brief neural stimulation; the pHRF is often used as the impulse response for linear analysis of neural excitation. Many fMRI studies have observed a negative BOLD response (NBR) that is often associated with neural suppression. However, the temporal dynamics of the NBR evoked by a brief stimulus, the negative HRF (nHRF), remains unclear. Here, a unilateral visual stimulus was presented in a slow event-related design to elicit both pHRFs in the stimulus representation (SR), and nHRFs elsewhere. The observed nHRFs were not inverted versions of the pHRF previously reported. They were characterized by a stronger initial negative response followed by a significantly later positive peak. In contralateral primary visual cortex (V1), these differences varied with eccentricity from the SR. Similar nHRFs were observed in ipsilateral V1 with less eccentricity variation. Experiments with the blocked version of the same stimulus confirmed that brain regions presenting the unexpected nHRF dynamics correspond to those presenting a strong NBR. These data demonstrated that shift-invariant temporal linearity did not hold for the NBR while confirming that the PBR maintained rough linearity. Modeling indicated that the observed nHRFs can be created by suppression of both blood flow and oxygen metabolism. Critically, the nHRF can be misinterpreted as a pHRF due to their similarity, which could confound linear analysis for event-related fMRI experiments.NEW & NOTEWORTHY We investigate dynamics of the negative hemodynamic response function (nHRF), the negative blood-oxygen-level-dependent (BOLD) response (NBR) evoked by a brief stimulus, in human early visual cortex. Here, we show that the nHRFs are not inverted versions of the corresponding pHRFs. The nHRF has complex dynamics that varied significantly with eccentricity. The results also show shift-invariant temporal linearity does not hold for the NBR.

A cortical rat hemodynamic response function for improved detection of BOLD activation under common experimental conditions.

For a reliable estimation of neuronal activation based on BOLD fMRI measurements an accurate model of the hemodynamic response is essential. Since a large part of basic neuroscience research is based on small animal data, it is necessary to characterize a hemodynamic response function (HRF) which is optimized for small animals. Therefore, we have determined and investigated the HRFs of rats obtained under a variety of experimental conditions in the primary somatosensory cortex. Measurements were performed on animals of different sex and strain, under different anesthetics, with and without ventilation and using different stimulation modalities. All modalities of stimulation used in this study induced neuronal activity in the primary somatosensory cortex or in subcortical regions. Since the HRFs of the BOLD responses in the primary somatosensory cortex showed a close concordance for the different conditions, we were able to determine a cortical rat HRF. This HRF is based on 143 BOLD measurements of 76 rats and can be used for statistical parametric mapping. It showed substantially faster progression than the human HRF, with a maximum after 2.8 ± 0.8 s, and a following undershoot after 6.1 ± 3.7 s. If the rat HRF was used statistical analysis of rat data showed a significantly improved detection performance in the somatosensory cortex in comparison to the commonly used HRF based on measurements in humans.

Changes in Hemodynamic Response Function Resulting From Chronic Alcohol Consumption.

BACKGROUND: Functional MRI (fMRI) task-related analyses rely on an estimate of the brain's hemodynamic response function (HRF) to model the brain's response to events. Although changes in the HRF have been found after acute alcohol administration, the effects of heavy chronic alcohol consumption on the HRF have not been explored, and the potential benefits or pitfalls of estimating each individual's HRF on fMRI analyses of chronic alcohol use disorder (AUD) are not known. METHODS: Participants with AUD and controls (CTL) received structural, functional, and vascular scans. During fMRI, participants were cued to tap their fingers, and averaged responses were extracted from the motor cortex. Curve fitting on these HRFs modeled them as a difference between 2 gamma distributions, and the temporal occurrence of the main peak and undershoot of the HRF was computed from the mean of the first and second gamma distributions, respectively. RESULTS: ANOVA and regression analyses found that the timing of the HRF undershoot increased significantly as a function of total lifetime drinking. Although gray matter volume in the motor cortex decreased with lifetime drinking, this was not sufficient to explain undershoot timing shifts, and vascular factors measured in the motor cortex did not differ among groups. Comparison of random-effects analyses using custom-fitted and canonical HRFs for CTL and AUD groups showed better results throughout the brain for custom-fitted versus canonical HRFs for CTL subjects. For AUD subjects, the same was true except for the basal ganglia. CONCLUSIONS: These findings suggest that excessive alcohol consumption is associated with changes in the HRF undershoot. HRF changes could provide a possible biomarker for the effects of lifetime drinking on brain function. Changes in HRF topography affect fMRI activation measures, and subject-specific HRFs generally improve fMRI activation results.

Whole brain measurements of the positive BOLD response variability during a finger tapping task at 7 T show regional differences in its profiles.

PURPOSE: The positive BOLD response can vary across brain regions. Here, the positive BOLD responses of motor regions, including the cerebellum, were investigated by fast fMRI acquisition. METHODS: The participants were asked to perform an event-related finger-tapping task in a 7T MRI scanner during a fast 3D-EPI controlled aliasing in parallel imaging acquisition protocol (CAIPI; TR = 399 ms). The positive BOLD responses of 6 motor regions were extracted and their timings and shapes measured. RESULTS: Compared with other brain regions, the positive BOLD responses in the cerebellum and secondary somatosensory cortex showed delayed onsets, but no differences were observed for the time to-peak. Additionally, variations of the undershoot and main peak amplitudes were also observed, and undershoot was quasi-absent in the cerebellum. CONCLUSION: This study confirms that care should be taken when drawing conclusions about neuronal activity from the BOLD signal, particularly for the cerebellum.

A low-rank multivariate general linear model for multi-subject fMRI data and a non-convex optimization algorithm for brain response comparison.

The focus of this paper is on evaluating brain responses to different stimuli and identifying brain regions with different responses using multi-subject, stimulus-evoked functional magnetic resonance imaging (fMRI) data. To jointly model many brain voxels' responses to designed stimuli, we present a new low-rank multivariate general linear model (LRMGLM) for stimulus-evoked fMRI data. The new model not only is flexible to characterize variation in hemodynamic response functions (HRFs) across different regions and stimulus types, but also enables information "borrowing" across voxels and uses much fewer parameters than typical nonparametric models for HRFs. To estimate the proposed LRMGLM, we introduce a new penalized optimization function, which leads to temporally and spatially smooth HRF estimates. We develop an efficient optimization algorithm to minimize the optimization function and identify the voxels with different responses to stimuli. We show that the proposed method can outperform several existing voxel-wise methods by achieving both high sensitivity and specificity. We apply the proposed method to the fMRI data collected in an emotion study, and identify anterior dACC to have different responses to a designed threat and control stimuli.

Macroscale variation in resting-state neuronal activity and connectivity assessed by simultaneous calcium imaging, hemodynamic imaging and electrophysiology.

Functional imaging of spontaneous activity continues to play an important role in the field of connectomics. The most common imaging signal used for these experiments is the blood-oxygen-level-dependent (BOLD) functional MRI (fMRI) signal, but how this signal relates to spontaneous neuronal activity remains incompletely understood. Genetically encoded calcium indicators represent a promising tool to study this problem, as they can provide brain-wide measurements of neuronal activity compared to point measurements afforded by electrophysiological recordings. However, the relationship between the calcium signal and neurophysiological parameters at the mesoscopic scale requires further systematic characterization. Therefore, we collected simultaneous resting-state measurements of electrophysiology, along with calcium and hemodynamic imaging, in lightly anesthetized mice to investigate two aims. First, we examined the relationship between each imaging signal and the simultaneously recorded electrophysiological signal in a single brain region, finding that both signals are better correlated with multi-unit activity compared to local field potentials, with the calcium signal possessing greater signal-to-noise ratio and regional specificity. Second, we used the resting-state imaging data to model the relationship between the calcium and hemodynamic signals across the brain. We found that this relationship varied across brain regions in a way that is consistent across animals, with delays increasing by600 ms towards posterior cortical regions. Furthermore, while overall functional connectivity (FC) measured by the hemodynamic signal is significantly correlated with FC measured by calcium, the two estimates were found to be significantly different. We hypothesize that these differences arise at least in part from the observed regional variation in the hemodynamic response. In total, this work highlights some of the caveats needed in interpreting hemodynamic-based measurements of FC, as well as the need for improved modeling methods to reduce this potential source of bias.

Hemodynamic response function (HRF) variability confounds resting-state fMRI functional connectivity.

PURPOSE: fMRI is the convolution of the hemodynamic response function (HRF) and unmeasured neural activity. HRF variability (HRFv) across the brain could, in principle, alter functional connectivity (FC) estimates from resting-state fMRI (rs-fMRI). Given that HRFv is driven by both neural and non-neural factors, it is problematic when it confounds FC. However, this aspect has remained largely unexplored even though FC studies have grown exponentially. We hypothesized that HRFv confounds FC estimates in the brain's default-mode-network. METHODS: We tested this hypothesis using both simulations (where the ground truth is known and modulated) as well as rs-fMRI data obtained in a 7T MRI scanner (N = 47, healthy). FC was obtained using 2 pipelines: data with hemodynamic deconvolution (DC) to estimate the HRF and minimize HRFv, and data with no deconvolution (NDC, HRFv-ignored). DC and NDC FC networks were compared, along with regional HRF differences, revealing potential false connectivities that resulted from HRFv. RESULTS: We found evidence supporting our hypothesis using both simulations and experimental data. With simulations, we found that HRFv could cause a change of up to 50% in FC. With rs-fMRI, several potential false connectivities attributable to HRFv, with majority connections being between different lobes, were identified. We found a double exponential relationship between the magnitude of HRFv and its impact on FC, with a mean/median error of 30.5/11.5% caused in FC by HRF confounds. CONCLUSION: HRFv, if ignored, could cause identification of false FC. FC findings from HRFv-ignored data should be interpreted cautiously. We suggest deconvolution to minimize HRFv.

Effect of confounding variables on hemodynamic response function estimation using averaging and deconvolution analysis: An event-related NIRS study.

Slow and rapid event-related designs are used in fMRI and functional near-infrared spectroscopy (fNIRS) experiments to temporally characterize the brain hemodynamic response to discrete events. Conventional averaging (CA) and the deconvolution method (DM) are the two techniques commonly used to estimate the Hemodynamic Response Function (HRF) profile in event-related designs. In this study, we conducted a series of simulations using synthetic and real NIRS data to examine the effect of the main confounding factors, including event sequence timing parameters, different types of noise, signal-to-noise ratio (SNR), temporal autocorrelation and temporal filtering on the performance of these techniques in slow and rapid event-related designs. We also compared systematic errors in the estimates of the fitted HRF amplitude, latency and duration for both techniques. We further compared the performance of deconvolution methods based on Finite Impulse Response (FIR) basis functions and gamma basis sets. Our results demonstrate that DM was much less sensitive to confounding factors than CA. Event timing was the main parameter largely affecting the accuracy of CA. In slow event-related designs, deconvolution methods provided similar results to those obtained by CA. In rapid event-related designs, our results showed that DM outperformed CA for all SNR, especially above -5 dB regardless of the event sequence timing and the dynamics of background NIRS activity. Our results also show that periodic low-frequency systemic hemodynamic fluctuations as well as phase-locked noise can markedly obscure hemodynamic evoked responses. Temporal autocorrelation also affected the performance of both techniques by inducing distortions in the time profile of the estimated hemodynamic response with inflated t-statistics, especially at low SNRs. We also found that high-pass temporal filtering could substantially affect the performance of both techniques by removing the low-frequency components of HRF profiles. Our results emphasize the importance of characterization of event timing, background noise and SNR when estimating HRF profiles using CA and DM in event-related designs.

BOLD neurovascular coupling does not change significantly with normal aging.

Studies of cognitive function that compare the blood oxygenation level dependent (BOLD) signal across age groups often require the assumption that neurovascular coupling does not change with age. Tests of this assumption have produced mixed results regarding the strength of the coupling and its relative time course. Using deconvolution, we found that age does not have a significant effect on the time course of the hemodynamic impulse response function or on the slope of the BOLD versus stimulus duration relationship. These results suggest that in cognitive studies of healthy aging, group differences in BOLD activation are likely due to age-related changes in cognitive-neural interactions and information processing rather than to impairments in neurovascular coupling. Hum Brain Mapp 38:3538-3551, 2017. © 2017 Wiley Periodicals, Inc.

Impact of task-related changes in heart rate on estimation of hemodynamic response and model fit.

The blood oxygen level dependent (BOLD) signal, as measured using functional magnetic resonance imaging (fMRI), is widely used as a proxy for changes in neural activity in the brain. Physiological variables such as heart rate (HR) and respiratory variation (RV) affect the BOLD signal in a way that may interfere with the estimation and detection of true task-related neural activity. This interference is of particular concern when these variables themselves show task-related modulations. We first establish that a simple movement task reliably induces a change in HR but not RV. In group data, the effect of HR on the BOLD response was larger and more widespread throughout the brain than were the effects of RV or phase regressors. The inclusion of HR regressors, but not RV or phase regressors, had a small but reliable effect on the estimated hemodynamic response function (HRF) in M1 and the cerebellum. We next asked whether the inclusion of a nested set of physiological regressors combining phase, RV, and HR significantly improved the model fit in individual participants' data sets. There was a significant improvement from HR correction in M1 for the greatest number of participants, followed by RV and phase correction. These improvements were more modest in the cerebellum. These results indicate that accounting for task-related modulation of physiological variables can improve the detection and estimation of true neural effects of interest.

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

Functional near-infrared spectroscopy (fNIRS) can be combined with electroencephalography (EEG) to continuously monitor the hemodynamic signal evoked by epileptic events such as seizures or interictal epileptiform discharges (IEDs, aka spikes). As estimation methods assuming a canonical shape of the hemodynamic response function (HRF) might not be optimal, we sought to model patient-specific HRF (sHRF) with a simple deconvolution approach for IED-related analysis with EEG-fNIRS data. Furthermore, a quadratic term was added to the model to account for the nonlinearity in the response when IEDs are frequent. Prior to analyzing clinical data, simulations were carried out to show that the HRF was estimable by the proposed deconvolution methods under proper conditions. EEG-fNIRS data of five patients with refractory focal epilepsy were selected due to the presence of frequent clear IEDs and their unambiguous focus localization. For each patient, both the linear sHRF and the nonlinear sHRF were estimated at each channel. Variability of the estimated sHRFs was seen across brain regions and different patients. Compared with the SPM8 canonical HRF (cHRF), including these sHRFs in the general linear model (GLM) analysis led to hemoglobin activations with higher statistical scores as well as larger spatial extents on all five patients. In particular, for patients with frequent IEDs, nonlinear sHRFs were seen to provide higher sensitivity in activation detection than linear sHRFs. These observations support using sHRFs in the analysis of IEDs with EEG-fNIRS data.