Coupling of simultaneously acquired electrophysiological and haemodynamic responses during visual stimulation

Separating vascular and neuronal effects of age on fMRI BOLD signals

Accurate identification of brain function is necessary to understand the neurobiology of cognitive ageing, and thereby promote well-being across the lifespan. A common tool used to investigate neurocognitive ageing is functional magnetic resonance imaging (fMRI). However, although fMRI data are often interpreted in terms of neuronal activity, the blood oxygenation level-dependent (BOLD) signal measured by fMRI includes contributions of both vascular and neuronal factors, which change differentially with age. While some studies investigate vascular ageing factors, the results of these studies are not well known within the field of neurocognitive ageing and therefore vascular confounds in neurocognitive fMRI studies are common. Despite over 10 000 BOLD-fMRI papers on ageing, fewer than 20 have applied techniques to correct for vascular effects. However, neurovascular ageing is not only a confound in fMRI, but an important feature in its own right, to be assessed alongside measures of neuronal ageing. We review current approaches to dissociate neuronal and vascular components of BOLD-fMRI of regional activity and functional connectivity. We highlight emerging evidence that vascular mechanisms in the brain do not simply control blood flow to support the metabolic needs of neurons, but form complex neurovascular interactions that influence neuronal function in health and disease. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

Hemispheric differences in electrical and hemodynamic responses during hemifield visual stimulation with graded contrasts

A multimodal neuroimaging technique based on electroencephalography (EEG) and functional near-infrared spectroscopy (fNIRS) was used with horizontal hemifield visual stimuli with graded contrasts to investigate the retinotopic mapping more fully as well as to explore hemispheric differences in neuronal activity, the hemodynamic response, and the neurovascular coupling relationship in the visual cortex. The fNIRS results showed the expected activation over the contralateral hemisphere for both the left and right hemifield visual stimulations. However, the EEG results presented a paradoxical lateralization, with the maximal response located over the ipsilateral hemisphere but with the polarity inversed components located over the contralateral hemisphere. Our results suggest that the polarity inversion as well as the latency advantage over the contralateral hemisphere cause the amplitude of the VEP over the contralateral hemisphere to be smaller than that over the ipsilateral hemisphere. Both the neuronal and hemodynamic responses changed logarithmically with the level of contrast in the hemifield visual stimulations. Moreover, the amplitudes and latencies of the visual evoked potentials (VEPs) were linearly correlated with the hemodynamic responses despite differences in the slopes.

Correlation between electrical and hemodynamic responses during visual stimulation with graded contrasts

Brain functional activity involves complex cellular, metabolic, and vascular chain reactions, making it difficult to comprehend. Electroencephalography (EEG) and functional near infrared spectroscopy (fNIRS) have been combined into a multimodal neuroimaging method that captures both electrophysiological and hemodynamic information to explore the spatiotemporal characteristics of brain activity. Because of the significance of visually evoked functional activity in clinical applications, numerous studies have explored the amplitude of the visual evoked potential (VEP) to clarify its relationship with the hemodynamic response. However, relatively few studies have investigated the influence of latency, which has been frequently used to diagnose visual diseases, on the hemodynamic response. Moreover, because the latency and the amplitude of VEPs have different roles in coding visual information, investigating the relationship between latency and the hemodynamic response should be helpful. In this study, checkerboard reversal tasks with graded contrasts were used to evoke visual functional activity. Both EEG and fNIRS were employed to investigate the relationship between neuronal electrophysiological activities and the hemodynamic responses. The VEP amplitudes were linearly correlated with the hemodynamic response, but the VEP latency showed a negative linear correlation with the hemodynamic response.

The effect of ageing on fMRI: Correction for the confounding effects of vascular reactivity evaluated by joint fMRI and MEG in 335 adults

In functional magnetic resonance imaging (fMRI) research one is typically interested in neural activity. However, the blood‐oxygenation level‐dependent (BOLD) signal is a composite of both neural and vascular activity. As factors such as age or medication may alter vascular function, it is essential to account for changes in neurovascular coupling when investigating neurocognitive functioning with fMRI. The resting‐state fluctuation amplitude (RSFA) in the fMRI signal (rsfMRI) has been proposed as an index of vascular reactivity. The RSFA compares favourably with other techniques such as breath‐hold and hypercapnia, but the latter are more difficult to perform in some populations, such as older adults. The RSFA is therefore a candidate for use in adjusting for age‐related changes in vascular reactivity in fMRI studies. The use of RSFA is predicated on its sensitivity to vascular rather than neural factors; however, the extent to which each of these factors contributes to RSFA remains to be characterized. The present work addressed these issues by comparing RSFA (i.e., rsfMRI variability) to proxy measures of (i) cardiovascular function in terms of heart rate (HR) and heart rate variability (HRV) and (ii) neural activity in terms of resting state magnetoencephalography (rsMEG). We derived summary scores of RSFA, a sensorimotor task BOLD activation, cardiovascular function and rsMEG variability for 335 healthy older adults in the population‐based Cambridge Centre for Ageing and Neuroscience cohort (Cam‐CAN; http://www.cam-can.com). Mediation analysis revealed that the effects of ageing on RSFA were significantly mediated by vascular factors, but importantly not by the variability in neuronal activity. Furthermore, the converse effects of ageing on the rsMEG variability were not mediated by vascular factors. We then examined the effect of RSFA scaling of task‐based BOLD in the sensorimotor task. The scaling analysis revealed that much of the effects of age on task‐based activation studies with fMRI do not survive correction for changes in vascular reactivity, and are likely to have been overestimated in previous fMRI studies of ageing. The results from the mediation analysis demonstrate that RSFA is modulated by measures of vascular function and is not driven solely by changes in the variance of neural activity. Based on these findings we propose that the RSFA scaling method is articularly useful in large scale and longitudinal neuroimaging studies of ageing, or with frail participants, where alternative measures of vascular reactivity are impractical. Hum Brain Mapp 36:2248–2269, 2015. © 2015 The Authors Human Brain Mapping Published by Wiley Periodicals, Inc.

Comparison of BCG artifact removal methods for evoked responses in simultaneous EEG-fMRI

Simultaneous recording of electroencephalogram (EEG) and functional magnetic resonance imaging (fMRI) has gained attention due to the complimentary properties of the two imaging modalities. Their combined recording enables the study of brain function while taking advantage of the high temporal resolution of EEG and high spatial resolution of fMRI. However EEG data recorded inside the MR scanner is significantly contaminated by two main sources of artifacts: MR gradient artifacts and ballistocardiogram (BCG) artifacts. Most existing removal approaches for these artifacts fall into two main categories: average artifact subtraction (AAS) and optimal basis selection (OBS). While these techniques can improve the data quality significantly, highly effective removal of artifacts - particularly the BCG artifact - from the data is still lacking. Here, we compared two of the most commonly used algorithms for BCG artifact removal (OBS and AAS) based on the estimated signal-to-noise ratio (SNR) of auditory and visual evoked responses recorded during fMRI acquisition. We also further compared optimization of OBS for groups, and at the individual subject and run level. The results suggest that performance of the OBS algorithm can be significantly improved by choosing the optimum number of principal components. Furthermore, optimizing the number of principal components at the individual participant and run level results in significant improvements in the SNR of evoked responses compared to group optimization.

Oxygen Level and LFP in Task-Positive and Task-Negative Areas: Bridging BOLD fMRI and Electrophysiology

The human default mode network (DMN) shows decreased blood oxygen level dependent (BOLD) signals in response to a wide range of attention-demanding tasks. Our understanding of the specifics regarding the neural activity underlying these "task-negative" BOLD responses remains incomplete. We paired oxygen polarography, an electrode-based oxygen measurement technique, with standard electrophysiological recording to assess the relationship of oxygen and neural activity in task-negative posterior cingulate cortex (PCC), a hub of the DMN, and visually responsive task-positive area V3 in the awake macaque. In response to engaging visual stimulation, oxygen, LFP power, and multi-unit activity in PCC showed transient activation followed by sustained suppression. In V3, oxygen, LFP power, and multi-unit activity showed an initial phasic response to the stimulus followed by sustained activation. Oxygen responses were correlated with LFP power in both areas, although the apparent hemodynamic coupling between oxygen level and electrophysiology differed across areas. Our results suggest that oxygen responses reflect changes in LFP power and multi-unit activity and that either the coupling of neural activity to blood flow and metabolism differs between PCC and V3 or computing a linear transformation from a single LFP band to oxygen level does not capture the true physiological process.

Evidence that the negative BOLD response is neuronal in origin: a simultaneous EEG-BOLD-CBF study in humans

Unambiguous interpretation of changes in the BOLD signal is challenging because of the complex neurovascular coupling that translates changes in neuronal activity into the subsequent haemodynamic response. In particular, the neurophysiological origin of the negative BOLD response (NBR) remains incompletely understood. Here, we simultaneously recorded BOLD, EEG and cerebral blood flow (CBF) responses to 10 s blocks of unilateral median nerve stimulation (MNS) in order to interrogate the NBR. Both negative BOLD and negative CBF responses to MNS were observed in the same region of the ipsilateral primary sensorimotor cortex (S1/M1) and calculations showed that MNS induced a decrease in the cerebral metabolic rate of oxygen consumption (CMRO2) in this NBR region. The ∆CMRO2/∆CBF coupling ratio (n) was found to be significantly larger in this ipsilateral S1/M1 region (n=0.91±0.04, M=10.45%) than in the contralateral S1/M1 (n=0.65±0.03, M=10.45%) region that exhibited a positive BOLD response (PBR) and positive CBF response, and a consequent increase in CMRO2 during MNS. The fMRI response amplitude in ipsilateral S1/M1 was negatively correlated with both the power of the 8-13 Hz EEG mu oscillation and somatosensory evoked potential amplitude. Blocks in which the largest magnitude of negative BOLD and CBF responses occurred therefore showed greatest mu power, an electrophysiological index of cortical inhibition, and largest somatosensory evoked potentials. Taken together, our results suggest that a neuronal mechanism underlies the NBR, but that the NBR may originate from a different neurovascular coupling mechanism to the PBR, suggesting that caution should be taken in assuming the NBR simply represents the neurophysiological inverse of the PBR.

Coupling between resting cerebral perfusion and EEG

While several studies have investigated interactions between the electroencephalography (EEG) and functional magnetic resonance imaging BOLD signal fluctuations, less is known about the associations between EEG oscillations and baseline brain haemodynamics, and few studies have examined the link between EEG power outside the alpha band and baseline perfusion. Here we compare whole-brain arterial spin labelling perfusion MRI and EEG in a group of healthy adults (n = 16, ten females, median age: 27 years, range 21-48) during an eyes closed rest condition. Correlations emerged between perfusion and global average EEG power in low (delta: 2-4 Hz and theta: 4-7 Hz), middle (alpha: 8-13 Hz), and high (beta: 13-30 Hz and gamma: 30-45 Hz) frequency bands in both cortical and sub-cortical regions. The correlations were predominately positive in middle and high-frequency bands, and negative in delta. In addition, central alpha frequency positively correlated with perfusion in a network of brain regions associated with the modulation of attention and preparedness for external input, and central theta frequency correlated negatively with a widespread network of cortical regions. These results indicate that the coupling between average EEG power/frequency and local cerebral blood flow varies in a frequency specific manner. Our results are consistent with longstanding concepts that decreasing EEG frequencies which in general map onto decreasing levels of activation.

Correlation between single-trial visual evoked potentials and the blood oxygenation level dependent response in simultaneously recorded electroencephalography-functional magnetic resonance imaging

To compare different electroencephalography (EEG)-based regressors and their ability to predict the simultaneously recorded blood oxygenation level dependent response during blocked visual stimulation, simultaneous EEG-functional magnetic resonance imaging in 10 healthy volunteers was performed. The performance of different single-trial EEG regressors was compared in terms of predicting the measured blood oxygenation level dependent response. The EEG-based regressors were the amplitude and latency of the primary positive (P1) and negative (N2) peaks of the visual evoked potential, the combined P1-N2 amplitude, and the alpha power. Apart from peak latencies, all regressors showed significant positive or negative correlation with the blood oxygenation level dependent response in visual cortex. In addition, several EEG-based regressors were found to predict blood oxygenation level dependent variations in different occipital and extraoccipital cortical areas not explained by the boxcar regressor. The results suggest that the P1-N2 regressor is the best EEG-based regressor to model the visual paradigm, but when looking for additional effects like habituation or attention modulation that cannot be modeled by the boxcar regressor, it is better to include regressors based on individual peaks or alpha power.