Exploring the post-stimulus undershoot with spin-echo fMRI: implications for models of neurovascular response

BOLD Response is more than just magnitude: Improving detection sensitivity through capturing hemodynamic profiles

Typical fMRI analyses often assume a canonical hemodynamic response function (HRF) that primarily focuses on the peak height of the overshoot, neglecting other morphological aspects. Consequently, reported analyses often reduce the overall response curve to a single scalar value. In this study, we take a data-driven approach to HRF estimation at the whole-brain voxel level, without assuming a response profile at the individual level. We then employ a roughness penalty at the population level to estimate the response curve, aiming to enhance predictive accuracy, inferential efficiency, and cross-study reproducibility. By examining a fast event-related FMRI dataset, we demonstrate the shortcomings and information loss associated with adopting the canonical approach. Furthermore, we address the following key questions: 1) To what extent does the HRF shape vary across different regions, conditions, and participant groups? 2) Does the data-driven approach improve detection sensitivity compared to the canonical approach? 3) Can analyzing the HRF shape help validate the presence of an effect in conjunction with statistical evidence? 4) Does analyzing the HRF shape offer evidence for whole-brain response during a simple task?

Determining Excitatory and Inhibitory Neuronal Activity from Multimodal fMRI Data Using a Generative Hemodynamic Model

Hemodynamic responses, in general, and the blood oxygenation level-dependent (BOLD) fMRI signal, in particular, provide an indirect measure of neuronal activity. There is strong evidence that the BOLD response correlates well with post-synaptic changes, induced by changes in the excitatory and inhibitory (E-I) balance between active neuronal populations. Typical BOLD responses exhibit transients, such as the early-overshoot and post-stimulus undershoot, that can be linked to transients in neuronal activity, but they can also result from vascular uncoupling between cerebral blood flow (CBF) and venous cerebral blood volume (venous CBV). Recently, we have proposed a novel generative hemodynamic model of the BOLD signal within the dynamic causal modeling framework, inspired by physiological observations, called P-DCM (Havlicek et al., 2015). We demonstrated the generative model's ability to more accurately model commonly observed neuronal and vascular transients in single regions but also effective connectivity between multiple brain areas (Havlicek et al., 2017b). In this paper, we additionally demonstrate the versatility of the generative model to jointly explain dynamic relationships between neuronal and hemodynamic physiological variables underlying the BOLD signal using multi-modal data. For this purpose, we utilized three distinct data-sets of experimentally induced responses in the primary visual areas measured in human, cat, and monkey brain, respectively: (1) CBF and BOLD responses; (2) CBF, total CBV, and BOLD responses (Jin and Kim, 2008); and (3) positive and negative neuronal and BOLD responses (Shmuel et al., 2006). By fitting the generative model to the three multi-modal experimental data-sets, we showed that the presence or absence of dynamic features in the BOLD signal is not an unambiguous indication of presence or absence of those features on the neuronal level. Nevertheless, the generative model that takes into account the dynamics of the physiological mechanisms underlying the BOLD response allowed dissociating neuronal from vascular transients and deducing excitatory and inhibitory neuronal activity time-courses from BOLD data alone and from multi-modal data.

Metamemory function in chemotherapy-treated patients with breast cancer: an explanation for the dissociation between subjective and objective memory measures?

BACKGROUND

The purpose of this study was to determine if a deficit in metamemory could account for the disparity between subjective and objective measures of memory function commonly observed in patients with breast cancer (BC). Metamemory refers to the awareness and management of one's own memory function. It is considered an aspect of executive functioning, one of the most common areas of cognitive compromise associated with BC and its treatment.

METHODS

Fifty-four women with early stage BC who had recently completed chemotherapy were compared with 54 healthy women matched on age and education. Cognitive function was objectively assessed with a neuropsychological test battery and subjectively assessed with the Functional Assessment of Cancer Therapy Cognitive Scale. Metamemory was assessed with a Feeling of Knowing (FOK) paradigm.

RESULTS

The patients with BC scored significantly lower than the controls on both the objective and subjective cognitive measures, as well as on free recall and recognition conditions of the FOK, suggesting some decline in primary memory functions such as working memory, encoding, and retrieval. The discrepancy between the objective and subjective measures was larger in the patients with BC than in the controls, but there was no difference between the groups on the FOK metamemory index.

CONCLUSIONS

Discrepancy in objective and subjective measures of cognition in patients with BC cannot be accounted for in terms of a deficit in meta-cognition. Objective and subjective measures are complementary, and a comprehensive cognitive assessment in patients with BC requires both. Copyright © 2015 John Wiley & Sons, Ltd.

The neural basis of humour comprehension and humour appreciation: The roles of the temporoparietal junction and superior frontal gyrus

Psychological well-being and social acumen benefit from the recognition of humourous intent and its enjoyment. The enjoyment of humour requires recognition, but humour recognition is not necessarily accompanied by humour enjoyment. Humour recognition is crucial during social interactions, while the associated enjoyment is less critical. Few neuroimaging studies have explicitly differentiated between the neural foundations of humour comprehension and humour appreciation. Among such studies, design limitations have obscured the specification of neural correlates to humour comprehension or appreciation. We implemented a trichotomous response option to address these design limitations. Twenty-four participants rated 120 comics (90 unaltered with humourous intent and 30 caption-altered without humourous intent) as either funny jokes (FJ), not funny jokes but intended to be funny (NFJ), or not intended to be funny or non-jokes (NJ). We defined humour comprehension by NFJ minus NJ and humour appreciation by FJ minus NFJ. We measured localized blood oxygen level dependent (BOLD) neural responses with a 3T MRI scanner. We tested for BOLD responses in humour comprehension brain regions of interest (ROIs), humour appreciation ROIs, and across the whole-brain. We found significant NFJ-NJ BOLD responses in our humour comprehension ROIs and significant FJ-NFJ BOLD responses in select humour appreciation ROIs. One key finding is that comprehension accuracy levels correlated with humour-comprehension responses in the left temporo-parietal junction (TPJ). This finding represents a novel and precise neural linkage to humour comprehension. A second key finding is that the superior frontal gyrus (SFG) was uniquely associated with humour-appreciation. The SFG response suggests that complex cognitive processing underlies humour appreciation and that current models of humour appreciation be revised. Finally, our research design provides an operational distinction between humour comprehension and appreciation and a sensitive measure of individual differences in humour comprehension accuracy.

Neuronal or hemodynamic? Grappling with the functional MRI signal

Magnetic resonance imaging (MRI) and functional MRI (fMRI) continue to advance because creative physicists, engineers, neuroscientists, clinicians, and physiologists find new ways for extracting more information from the signal. Innovations in pulse sequence design, paradigm design, and processing methods have advanced the field and firmly established fMRI as a cornerstone for understanding the human brain. In this article, the field of fMRI is described through consideration of the central problem of separating hemodynamic from neuronal information. Discussed here are examples of how pulse sequences, activation paradigms, and processing methods are integrated such that novel, high-quality information can be obtained. Examples include the extraction of information such as activation onset latency, metabolic rate, neuronal adaptation, vascular patency, vessel diameter, vigilance, and subvoxel activation. Experimental measures include time series latency, hemodynamic shape, MR phase, multivoxel patterns, ratios of activation-related R2* to R2, metabolic rate changes, fluctuation correlations and frequencies, changes in fluctuation correlations and frequencies over time, resting correlation states, echo time dependence, and more.

Investigating the properties of the hemodynamic response function after mild traumatic brain injury

Abstract Although several functional magnetic resonance imaging (fMRI) studies have been conducted in human models of mild traumatic brain injury (mTBI), to date no studies have explicitly examined how injury may differentially affect both the positive phase of the hemodynamic response function (HRF) as well as the post-stimulus undershoot (PSU). Animal models suggest that the acute and semi-acute stages of mTBI are associated with significant disruptions in metabolism and to the microvasculature, both of which could impact on the HRF. Therefore, fMRI data were collected on a cohort of 30 semi-acute patients with mTBI (16 males; 27.83±9.97 years old; 13.00±2.18 years of education) and 30 carefully matched healthy controls (HC; 16 males; 27.17±10.08 years old; 13.37±2.31 years of education) during a simple sensory-motor task. Patients reported increased cognitive, somatic, and emotional symptoms relative to controls, although no group differences were detected on traditional neuropsychological examination. There were also no differences between patients with mTBI and controls on fMRI data using standard analytic techniques, although mTBI exhibited a greater volume of activation during the task qualitatively. A significant Group×Time interaction was observed in the right supramarginal gyrus, bilateral primary and secondary visual cortex, and the right parahippocampal gyrus. The interaction was the result of an earlier time-to-peak and positive magnitude shift throughout the estimated HRF in patients with mTBI relative to HC. This difference in HRF shape combined with the greater volume of activated tissue may be indicative of a potential compensatory mechanism to injury. The current study demonstrates that direct examination and modeling of HRF characteristics beyond magnitude may provide additional information about underlying neuropathology that is not available with more standard fMRI analyses.

Sensitivity of BOLD response to increasing visual contrast: spin echo versus gradient echo EPI

Previous evidence showed that spin-echo (SE) BOLD signals offer an increased linearity and promptness with respect to gradient-echo (GE) acquisition, possibly providing a more accurate estimate of the amplitude of neuronal activity. However there is no evidence that the two sequences differ in representing different activation levels due to changes in stimulus intensity. In this study at 3T we compared GE and SE BOLD responses to visual stimuli at increasing contrast levels (5%, 20%, 60%, and 100%). Both sequences showed a monotonic increase of the BOLD response with stimulus contrast. While the larger sensitivity of GE yielded overall larger signal changes, step-wise increase in activation for GE was significant only when comparing 20% with 5% contrast, whereas for SE a significant increase was observed also when comparing 60% with 20% contrast. Moreover, BOLD responses normalized to the lowest contrast showed that relative increases of SE fMRI signal with increasing stimulus strength are larger with respect to the corresponding values of GE signal. This difference was observed also when excluding voxels attributed to large vessels, suggesting a non negligible role of the extravascular contribution to the modulation of SE fMRI signal with stimulus intensity. These results are shown to be in agreement with theoretical predictions that we derived from a recently proposed model of GE and SE functional signals. The present findings suggest that, despite the limited increase in functional localization accuracy at low field, SE fMRI might offer a potential advantage in distinguishing different levels of stimulus-evoked brain activity.

The BOLD post-stimulus undershoot, one of the most debated issues in fMRI

This paper provides a brief overview of how we got involved in fMRI work and of our efforts to elucidate the mechanisms underlying BOLD signal changes. The phenomenon discussed here in particular is the post-stimulus undershoot (PSU), the interpretation of which has captivated many fMRI scientists and is still under debate to date. This controversy is caused both by the convoluted physiological origin of the BOLD effect, which allows many possible explanations, and the lack of comprehensive data in the early years. BOLD effects reflect changes in cerebral blood flow (CBF), volume (CBV), metabolic rate of oxygen (CMRO(2)), and hematocrit fraction (Hct). However, the size of such effects is modulated by vascular origin such as intravascular, extravascular, macro and microvascular, venular and capillary, the relative contributions of which depend not only on the spatial resolution of the measurements, but also on stimulus duration, on magnetic field strength and on whether spin echo (SE) or gradient echo (GRE) detection is used. The two most dominant explanations of the PSU have been delayed vascular compliance (first venular, later arteriolar, and recently capillary) and sustained increases in CMRO(2), while post-activation reduction in CBF is a distant third. MRI has the capability to independently measure CBF and arteriolar, venous, and total CBV contributions in humans and animals, which has been of great assistance in improving the understanding of BOLD phenomena. Using currently available MRI and optical data, we conclude that the predominant PSU origin is a sustained increase in CMRO(2). However, some contributions from delayed vascular compliance are likely, and small CBF undershoot contributions that are difficult to detect with current arterial spin labeling technology can also not be excluded. The relative contribution of these different processes, which are not mutually exclusive and can act together, is likely to vary with stimulus duration and type.

Dynamic models of BOLD contrast

This personal recollection looks at the evolution of ideas about the dynamics of the blood oxygenation level dependent (BOLD) signal, with an emphasis on the balloon model. From the first detection of the BOLD response it has been clear that the signal exhibits interesting dynamics, such as a pronounced and long-lasting post-stimulus undershoot. The BOLD response, reflecting a change in local deoxyhemoglobin, is a combination of a hemodynamic response, related to changes in blood flow and venous blood volume, and a metabolic response related to oxygen metabolism. Modeling is potentially a way to understand the complex path from changes in neural activity to the BOLD signal. In the early days of fMRI it was hoped that the hemodynamic/metabolic response could be modeled in a unitary way, with blood flow, oxygen metabolism, and venous blood volume-the physiological factors that affect local deoxyhemoglobin-all tightly linked. The balloon model was an attempt to do this, based on the physiological ideas of limited oxygen delivery at baseline and a slow recovery of venous blood volume after the stimulus (the balloon effect), and this simple model of the physiology worked well to simulate the BOLD response. However, subsequent experiments suggest a more complicated picture of the underlying physiology, with blood flow and oxygen metabolism driven in parallel, possibly by different aspects of neural activity. In addition, it is still not clear whether the post-stimulus undershoot is a hemodynamic or a metabolic phenomenon, although the original venous balloon effect is unlikely to be the full explanation, and a flow undershoot is likely to be important. Although our understanding of the physics of the BOLD response is now reasonably solid, our understanding of the underlying physiological relationships is still relatively poor, and this is the primary hurdle for future models of BOLD dynamics.

Independent component analysis of broadband near-infrared spectroscopy data acquired on adult human head

The goal of this study was to investigate the ability of independent component analysis in the time-spectral domain to isolate physiological sources of functional near infrared spectroscopy signals. We apply independent component analysis to the broadband fNIRS data acquired on the human forehead at 650 different wavelengths between 700 nm and 950 nm. To induce cerebral oxygenation changes we use the breath holding paradigm. We found one major independent component during baseline and two major components during exercise. Each independent component corresponds to one oxy-hemoglobin and one deoxy-hemoglobin time courses. The corresponding characteristic spectra of changes in optical absorption suggested that one component represented vasodilation of cerebral arterioles while the delayed component represented the washout of deoxyhemoglobin either in cerebral capillaries and venules or in extra cerebral tissue. We found that both broadband and isolated wavelength data can produce similar independent components.

Physiological origin for the BOLD poststimulus undershoot in human brain: vascular compliance versus oxygen metabolism

The poststimulus blood oxygenation level-dependent (BOLD) undershoot has been attributed to two main plausible origins: delayed vascular compliance based on delayed cerebral blood volume (CBV) recovery and a sustained increased oxygen metabolism after stimulus cessation. To investigate these contributions, multimodal functional magnetic resonance imaging was employed to monitor responses of BOLD, cerebral blood flow (CBF), total CBV, and arterial CBV (CBV(a)) in human visual cortex after brief breath hold and visual stimulation. In visual experiments, after stimulus cessation, CBV(a) was restored to baseline in 7.9±3.4 seconds, and CBF and CBV in 14.8±5.0 seconds and 16.1±5.8 seconds, respectively, all significantly faster than BOLD signal recovery after undershoot (28.1±5.5 seconds). During the BOLD undershoot, postarterial CBV (CBV(pa), capillaries and venules) was slightly elevated (2.4±1.8%), and cerebral metabolic rate of oxygen (CMRO(2)) was above baseline (10.6±7.4%). Following breath hold, however, CBF, CBV, CBV(a) and BOLD signals all returned to baseline in ∼20 seconds. No significant BOLD undershoot, and residual CBV(pa) dilation were observed, and CMRO(2) did not substantially differ from baseline. These data suggest that both delayed CBV(pa) recovery and enduring increased oxidative metabolism impact the BOLD undershoot. Using a biophysical model, their relative contributions were estimated to be 19.7±15.9% and 78.7±18.6%, respectively.

Linear coupling of undershoot with BOLD response in ER-fMRI and nonlinear BOLD response in rapid-presentation ER-fMRI

In event-related (ER) BOLD-fMRI brain activation studies, understanding the relationship between the elicited BOLD signal and its underlying neuronal activity is essential for any quantitative interpretation of the neural events from the BOLD measurements. This requires a better understanding of the dynamic BOLD response. Besides the neuronal activity-induced positive BOLD response, the dynamic response is also characterized by a profound post-stimulus undershoot. The relationship between the positive response and the post-stimulus undershoot, however, remains poorly understood. Earlier studies using block-design paradigms with long stimulation durations (>10s) do not suggest a quantitative relationship. Using an ER paradigm, this study revealed a linear coupling between the positive BOLD response and the post-stimulus undershoot across the human visual cortex. The voxelwise linear coupling across the visual cortex strongly supports a homogeneous hemodynamic response in ER paradigms, though the BOLD response magnitude varies substantially over a wide range across the visual cortex. Although underlying neuronal activity is responsible for a BOLD response, the blood volume fraction affects the magnitude of the BOLD response; the larger the blood volume fraction, the larger the magnitude. This effect needs to be accounted for in any quantitative interpretation of the BOLD measurements. In the absence of nonlinear neuronal activities, the nonlinear vascular response renders the estimated BOLD responses smaller in rapid presentation (RP) ER paradigms compared to that in ER paradigms, and this reduction effect also needs to be considered when interpreting the estimated BOLD responses in RP-ER paradigms. Interestingly, this nonlinear effect might be simply accounted for by a scaling factor across the visual cortex.