Medullary vein architecture modulates the white matter BOLD cerebrovascular reactivity signal response to CO2: Observations from high-resolution T2* weighted imaging at 7T

Cerebrovascular Reactivity Following Spinal Cord Injury

Background

Spinal cord injuries (SCI) often result in cardiovascular issues, increasing the risk of stroke and cognitive deficits.

Objectives

This study assessed cerebrovascular reactivity (CVR) using functional magnetic resonance imaging (fMRI) during a hypercapnic challenge in SCI participants compared to noninjured controls.

Methods

Fourteen participants were analyzed (n = 8 with SCI [unless otherwise noted], median age = 44 years; n = 6 controls, median age = 33 years). CVR was calculated through fMRI signal changes.

Results

The results showed a longer CVR component (tau) in the grey matter of SCI participants (n = 7) compared to controls (median difference = 3.0 s; p < .05). Time since injury (TSI) correlated negatively with steady-state CVR in the grey matter and brainstem of SCI participants (RS = -0.81, p = .014; RS = -0.84, p = .009, respectively). Lower steady-state CVR in the brainstem of the SCI group (n = 7) correlated with lower diastolic blood pressure (RS = 0.76, p = .046). Higher frequency of hypotensive episodes (n = 7) was linked to lower CVR outcomes in the grey matter (RS = -0.86, p = .014) and brainstem (RS = -0.89, p = .007).

Conclusion

Preliminary findings suggest a difference in the dynamic CVR component, tau, between the SCI and noninjured control groups, potentially explaining the higher cerebrovascular health burden in SCI individuals. Exploratory associations indicate that longer TSI, lower diastolic blood pressure, and more hypotensive episodes may lead to poorer CVR outcomes. However, further research is necessary to establish causality and support these observations.

Contrast-free dynamic susceptibility contrast using sinusoidal and bolus oxygenation challenges

Deoxygenation-based dynamic susceptibility contrast (dDSC) MRI uses respiratory challenges as a source of endogenous contrast as an alternative to gadolinium injection. These gas challenges induce T2*-weighted MRI signal losses, after which tracer kinetics modeling was applied to calculate cerebral perfusion. This work compares three gas challenges, desaturation (transient hypoxia), resaturation (transient normoxia), and SineO2 (sinusoidal modulation of end-tidal oxygen pressures) in a cohort of 10 healthy volunteers (age 37 ± 11 years; 60% female). Perfusion estimates consisted of cerebral blood flow (CBF), cerebral blood volume (CBV), and mean transit time (MTT). Calculations were computed using a traditional tracer kinetics model in the time domain for desaturation and resaturation and in the frequency domain for SineO2. High correlations and limits of agreement were observed among the three deoxygenation-based paradigms for CBV, although MTT and CBF estimates varied with the hypoxic stimulus. Cross-modality correlation with gadolinium DSC was lower, particularly for MTT, but on a par with agreement between the other perfusion references. Overall, this work demonstrated the feasibility and reliability of oxygen respiratory challenges to measure brain perfusion. Additional work is needed to assess the utility of dDSC in the diagnostic evaluation of various pathologies such as ischemic strokes, brain tumors, and neurodegenerative diseases.

Hemodynamic imaging parameters in brain metastases patients - Agreement between multi-delay ASL and hypercapnic BOLD

Arterial spin labeling (ASL) MRI is a routine clinical imaging technique that provides quantitative cerebral blood flow (CBF) information. A related technique is blood oxygenation level-dependent (BOLD) MRI during hypercapnia, which can assess cerebrovascular reactivity (CVR). ASL is weighted towards arteries, whereas BOLD is weighted towards veins. Their associated parameters in heterogeneous tissue types or under different hemodynamic conditions remains unclear. Baseline multi-delay ASL MRI and BOLD MRI during hypercapnia were performed in fourteen patients with brain metastases. In the ROI analysis, the CBF and CVR values were positively correlated in regions showing sufficient reserve capacity (i.e. non-steal regions, rrm = 0.792). Additionally, longer hemodynamic lag times were related to lower baseline CBF (rrm = -0.822) and longer arterial arrival time (AAT; rrm = 0.712). In contrast, in regions exhibiting vascular steal an inverse relationship was found with higher baseline CBF related to more negative CVR (rrm = -0.273). These associations were confirmed in voxelwise analyses. The relationship between CBF, AAT and CVR measures seems to be dependent on the vascular status of the underlying tissue. Healthy tissue relationships do not hold in tissues experiencing impaired or exhausted autoregulation. CVR metrics can possibly identify at-risk areas before perfusion deficiencies become visible on ASL MRI, specifically within vascular steal regions.

Cerebrovascular reactivity measurements using 3T BOLD MRI and a fixed inhaled CO2 gas challenge: Repeatability and impact of processing strategy

Introduction: Cerebrovascular reactivity (CVR) measurements using blood oxygen level dependent (BOLD) magnetic resonance imaging (MRI) are commonly used to assess the health of cerebral blood vessels, including in patients with cerebrovascular diseases; however, evidence and consensus regarding reliability and optimal processing are lacking. We aimed to assess the repeatability, accuracy and precision of voxel- and region-based CVR measurements at 3 T using a fixed inhaled (FI) CO₂ stimulus in a healthy cohort. Methods: We simulated the effect of noise, delay constraints and voxel- versus region-based analysis on CVR parameters. Results were verified in 15 healthy volunteers (28.1±5.5 years, female: 53%) with a test-retest MRI experiment consisting of two CVR scans. CVR magnitude and delay in grey matter (GM) and white matter were computed for both analyses assuming a linear relationship between the BOLD signal and time-shifted end-tidal CO₂ (EtCO₂) profile. Results: Test-retest repeatability was high [mean (95% CI) inter-scan difference: -0.01 (-0.03, -0.00) %/mmHg for GM CVR magnitude; -0.3 (-1.2,0.6) s for GM CVR delay], but we detected a small systematic reduction in CVR magnitude at scan 2 versus scan 1, accompanied by a greater EtCO2 change [±1.0 (0.4,1.5) mmHg] and lower heart rate [-5.5 (-8.6,-2.4] bpm]. CVR magnitude estimates were higher for voxel- versus region-based analysis [difference in GM: ±0.02 (0.01,0.03) %/mmHg]. Findings were supported by simulation results, predicting a positive bias for voxel-based CVR estimates dependent on temporal contrast-to-noise ratio and delay fitting constraints and an underestimation for region-based CVR estimates. Discussion: BOLD CVR measurements using FI stimulus have good within-day repeatability in healthy volunteers. However, measurements may be influenced by physiological effects and the analysis protocol. Voxel-based analyses should be undertaken with care due to potential for systematic bias; region-based analyses are more reliable in such cases.

Using carpet plots to analyze blood transit times in the brain during hypercapnic challenge magnetic resonance imaging

Blood arrival time and blood transit time are useful metrics in characterizing hemodynamic behaviors in the brain. Functional magnetic resonance imaging in combination with a hypercapnic challenge has been proposed as a non-invasive imaging tool to determine blood arrival time and replace dynamic susceptibility contrast (DSC) magnetic resonance imaging, a current gold-standard imaging tool with the downsides of invasiveness and limited repeatability. Using a hypercapnic challenge, blood arrival times can be computed by cross-correlating the administered CO₂ signal with the fMRI signal, which increases during elevated CO₂ due to vasodilation. However, whole-brain transit times derived from this method can be significantly longer than the known cerebral transit time for healthy subjects (nearing 20 s vs. the expected 5-6 s). To address this unrealistic measurement, we here propose a novel carpet plot-based method to compute improved blood transit times derived from hypercapnic blood oxygen level dependent fMRI, demonstrating that the method reduces estimated blood transit times to an average of 5.32 s. We also investigate the use of hypercapnic fMRI with cross-correlation to compute the venous blood arrival times in healthy subjects and compare the computed delay maps with DSC-MRI time to peak maps using the structural similarity index measure (SSIM). The strongest delay differences between the two methods, indicated by low structural similarity index measure, were found in areas of deep white matter and the periventricular region. SSIM measures throughout the remainder of the brain reflected a similar arrival sequence derived from the two methods despite the exaggerated spread of voxel delays computed using CO₂ fMRI.

Quantifying cerebral blood arrival times using hypoxia-mediated arterial BOLD contrast

Cerebral blood arrival and tissue transit times are sensitive measures of the efficiency of tissue perfusion and can provide clinically meaningful information on collateral blood flow status. We exploit the arterial blood oxygen level dependent (BOLD) signal contrast established by precisely decreasing, and then increasing, arterial hemoglobin saturation using respiratory re-oxygenation challenges to quantify arterial blood arrival times throughout the brain. We term this approach the Step Hemoglobin re-Oxygenation Contrast Stimulus (SHOCS). Carpet plot analysis yielded measures of signal onset (blood arrival), global transit time (gTT) and calculations of relative total blood volume. Onset times averaged across 12 healthy subjects were 1.1 ± 0.4 and 1.9 ± 0.6 for cortical gray and deep white matter, respectively. The average whole brain gTT was 4.5 ± 0.9 seconds. The SHOCS response was 1.7 fold higher in grey versus white matter; in line with known differences in tissue-specific blood volume fraction. SHOCS was also applied in a patient with unilateral carotid artery occlusion revealing ipsilateral prolonged signal onset with normal perfusion in the unaffected hemisphere. We anticipate that SHOCS will further inform on the extent of collateral blood flow in patients with upstream steno-occlusive vascular disease, including those already known to manifest reductions in vasodilatory reserve capacity or vascular steal.

The Spatiotemporal Dynamics of Cerebral Autoregulation in Functional Magnetic Resonance Imaging

The thigh-cuff release (TCR) maneuver is a physiological challenge that is widely used to assess dynamic cerebral autoregulation (dCA). It is often applied in conjunction with Transcranial Doppler ultrasound (TCD), which provides temporal information of the global flow response in the brain. This established method can only yield very limited insights into the regional variability of dCA, whereas functional MRI (fMRI) has the ability to reveal the spatial distribution of flow responses in the brain with high spatial resolution. The aim of this study was to use whole-brain blood-oxygenation-level-dependent (BOLD) fMRI to characterize the spatiotemporal dynamics of the flow response to the TCR challenge, and thus pave the way toward mapping dCA in the brain. We used a data driven approach to derive a novel basis set that was then used to provide a voxel-wise estimate of the TCR associated haemodynamic response function (HRF TCR ). We found that the HRF TCR evolves with a specific spatiotemporal pattern, with gray and white matter showing an asynchronous response, which likely reflects the anatomical structure of cerebral blood supply. Thus, we propose that TCR challenge fMRI is a promising method for mapping spatial variability in dCA, which will likely prove to be clinically advantageous.

Coupling between cerebrovascular oscillations and CSF flow fluctuations during wakefulness: An fMRI study

It is commonly believed that cerebrospinal fluid (CSF) movement is facilitated by blood vessel wall movements (i.e., hemodynamic oscillations) in the brain. A coherent pattern of low frequency hemodynamic oscillations and CSF movement was recently found during non-rapid eye movement (NREM) sleep via functional MRI. This finding raises other fundamental questions: 1) the explanation of coupling between hemodynamic oscillations and CSF movement from fMRI signals; 2) the existence of the coupling during wakefulness; 3) the direction of CSF movement. In this resting state fMRI study, we proposed a mechanical model to explain the coupling between hemodynamics and CSF movement through the lens of fMRI. Time delays between CSF movement and global hemodynamics were calculated. The observed delays between hemodynamics and CSF movement match those predicted by the model. Moreover, by conducting separate fMRI scans of the brain and neck, we confirmed the low frequency CSF movement at the fourth ventricle is bidirectional. Our finding also demonstrates that CSF movement is facilitated by changes in cerebral blood volume mainly in the low frequency range, even when the individual is awake.