Reproducibility of hypocapnic cerebrovascular reactivity measurements using BOLD fMRI in combination with a paced deep breathing task

Hemodynamic timing in resting-state and breathing-task BOLD fMRI

The blood flow response to a vasoactive stimulus demonstrates regional heterogeneity across both the healthy brain and in cerebrovascular pathology. The timing of a regional hemodynamic response is emerging as an important biomarker of cerebrovascular dysfunction, as well as a confound within fMRI analyses. Previous research demonstrated that hemodynamic timing is more robustly characterized when a larger systemic vascular response is evoked by a breathing challenge, compared to when only spontaneous fluctuations in vascular physiology are present (i.e., in resting-state data). However, it is not clear whether hemodynamic delays in these two conditions are physiologically interchangeable, and how methodological signal-to-noise factors may limit their agreement. To address this, we generated whole-brain maps of hemodynamic delays in nine healthy adults. We assessed the agreement of voxel-wise gray matter (GM) hemodynamic delays between two conditions: resting-state and breath-holding. We found that delay values demonstrated poor agreement when considering all GM voxels, but increasingly greater agreement when limiting analyses to voxels showing strong correlation with the GM mean time-series. Voxels showing the strongest agreement with the GM mean time-series were primarily located near large venous vessels, however these voxels explain some, but not all, of the observed agreement in timing. Increasing the degree of spatial smoothing of the fMRI data enhanced the correlation between individual voxel time-series and the GM mean time-series. These results suggest that signal-to-noise factors may be limiting the accuracy of voxel-wise timing estimates and hence their agreement between the two data segments. In conclusion, caution must be taken when using voxel-wise delay estimates from resting-state and breathing-task data interchangeably, and additional work is needed to evaluate their relative sensitivity and specificity to aspects of vascular physiology and pathology.

Comparing end-tidal CO2, respiration volume per time (RVT), and average gray matter signal for mapping cerebrovascular reactivity amplitude and delay with breath-hold task BOLD fMRI

Cerebrovascular reactivity (CVR), defined as the cerebral blood flow response to a vasoactive stimulus, is an imaging biomarker with demonstrated utility in a range of diseases and in typical development and aging processes. A robust and widely implemented method to map CVR involves using a breath-hold task during a BOLD fMRI scan. Recording end-tidal CO₂ (P_ETCO₂) changes during the breath-hold task is recommended to be used as a reference signal for modeling CVR amplitude in standard units (%BOLD/mmHg) and CVR delay in seconds. However, obtaining reliable P_ETCO₂ recordings requires equipment and task compliance that may not be achievable in all settings. To address this challenge, we investigated two alternative reference signals to map CVR amplitude and delay in a lagged general linear model (lagged-GLM) framework: respiration volume per time (RVT) and average gray matter BOLD response (GM-BOLD). In 8 healthy adults with multiple scan sessions, we compare spatial agreement of CVR maps from RVT and GM-BOLD to those generated with P_ETCO₂. We define a threshold to determine whether a P_ETCO₂ recording has "sufficient" quality for CVR mapping and perform these comparisons in 16 datasets with sufficient P_ETCO₂ and 6 datasets with insufficient P_ETCO₂. When P_ETCO₂ quality is sufficient, both RVT and GM-BOLD produce CVR amplitude maps that are nearly identical to those from P_ETCO₂ (after accounting for differences in scale), with the caveat they are not in standard units to facilitate between-group comparisons. CVR delays are comparable to P_ETCO₂ with an RVT regressor but may be underestimated with the average GM-BOLD regressor. Importantly, when P_ETCO₂ quality is insufficient, RVT and GM-BOLD CVR recover reasonable CVR amplitude and delay maps, provided the participant attempted the breath-hold task. Therefore, our framework offers a solution for achieving high quality CVR maps in both retrospective and prospective studies where sufficient P_ETCO₂ recordings are not available and especially in populations where obtaining reliable measurements is a known challenge (e.g., children). Our results have the potential to improve the accessibility of CVR mapping and to increase the prevalence of this promising metric of vascular health.

Lag-Optimized Blood Oxygenation Level Dependent Cerebrovascular Reactivity Estimates Derived From Breathing Task Data Have a Stronger Relationship With Baseline Cerebral Blood Flow

Cerebrovascular reactivity (CVR), an important indicator of cerebrovascular health, is commonly studied with the Blood Oxygenation Level Dependent functional MRI (BOLD-fMRI) response to a vasoactive stimulus. Theoretical and empirical evidence suggests that baseline cerebral blood flow (CBF) modulates BOLD signal amplitude and may influence BOLD-CVR estimates. We address how acquisition and modeling choices affect the relationship between baseline cerebral blood flow (bCBF) and BOLD-CVR: whether BOLD-CVR is modeled with the inclusion of a breathing task, and whether BOLD-CVR amplitudes are optimized for hemodynamic lag effects. We assessed between-subject correlations of average GM values and within-subject spatial correlations across cortical regions. Our results suggest that a breathing task addition to a resting-state acquisition, alongside lag-optimization within BOLD-CVR modeling, can improve BOLD-CVR correlations with bCBF, both between- and within-subjects, likely because these CVR estimates are more physiologically accurate. We report positive correlations between bCBF and BOLD-CVR, both between- and within-subjects. The physiological explanation of this positive correlation is unclear; research with larger samples and tightly controlled vasoactive stimuli is needed. Insights into what drives variability in BOLD-CVR measurements and related measurements of cerebrovascular function are particularly relevant when interpreting results in populations with altered vascular and/or metabolic baselines or impaired cerebrovascular reserve.

Reproducibility of cerebrovascular reactivity measurements: A systematic review of neuroimaging techniques. J Cereb Blood Flow Metab 2022;42:700-717. [PMID: 34806918 PMCID: PMC9254040 DOI: 10.1177/0271678x211056702]

Cerebrovascular reactivity (CVR), the capacity of the brain to increase cerebral blood flow (CBF) to meet changes in physiological demand, is an important biomarker to evaluate brain health. Typically, this brain "stress test" is performed by using a medical imaging modality to measure the CBF change between two states: at baseline and after vasodilation. However, since there are many imaging modalities and many ways to augment CBF, a wide range of CVR values have been reported. An understanding of CVR reproducibility is critical to determine the most reliable methods to measure CVR as a clinical biomarker. This review focuses on CVR reproducibility studies using neuroimaging techniques in 32 articles comprising 427 total subjects. The literature search was performed in PubMed, Embase, and Scopus. The review was conducted using the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA). We identified 5 factors of the experimental subjects (such as sex, blood characteristics, and smoking) and 9 factors of the measuring technique (such as the imaging modality, the type of the vasodilator, and the quantification method) that have strong effects on CVR reproducibility. Based on this review, we recommend several best practices to improve the reproducibility of CVR quantification in neuroimaging studies.

A practical modification to a resting state fMRI protocol for improved characterization of cerebrovascular function

Cerebrovascular reactivity (CVR), defined here as the Blood Oxygenation Level Dependent (BOLD) response to a CO₂ pressure change, is a useful metric of cerebrovascular function. Both the amplitude and the timing (hemodynamic lag) of the CVR response can bring insight into the nature of a cerebrovascular pathology and aid in understanding noise confounds when using functional Magnetic Resonance Imaging (fMRI) to study neural activity. This research assessed a practical modification to a typical resting-state fMRI protocol, to improve the characterization of cerebrovascular function. In 9 healthy subjects, we modelled CVR and lag in three resting-state data segments, and in data segments which added a 2–3 minute breathing task to the start of a resting-state segment. Two different breathing tasks were used to induce fluctuations in arterial CO₂ pressure: a breath-hold task to induce hypercapnia (CO₂ increase) and a cued deep breathing task to induce hypocapnia (CO₂ decrease). Our analysis produced voxel-wise estimates of the amplitude (CVR) and timing (lag) of the BOLD-fMRI response to CO₂ by systematically shifting the CO₂ regressor in time to optimize the model fit. This optimization inherently increases gray matter CVR values and fit statistics. The inclusion of a simple breathing task, compared to a resting-state scan only, increases the number of voxels in the brain that have a significant relationship between CO₂ and BOLD-fMRI signals, and improves our confidence in the plausibility of voxel-wise CVR and hemodynamic lag estimates. We demonstrate the clinical utility and feasibility of this protocol in an incidental finding of Moyamoya disease, and explore the possibilities and challenges of using this protocol in younger populations. This hybrid protocol has direct applications for CVR mapping in both research and clinical settings and wider applications for fMRI denoising and interpretation.

ICA-based denoising strategies in breath-hold induced cerebrovascular reactivity mapping with multi echo BOLD fMRI

Performing a BOLD functional MRI (fMRI) acquisition during breath-hold (BH) tasks is a non-invasive, robust method to estimate cerebrovascular reactivity (CVR). However, movement and breathing-related artefacts caused by the BH can substantially hinder CVR estimates due to their high temporal collinearity with the effect of interest, and attention has to be paid when choosing which analysis model should be applied to the data. In this study, we evaluate the performance of multiple analysis strategies based on lagged general linear models applied on multi-echo BOLD fMRI data, acquired in ten subjects performing a BH task during ten sessions, to obtain subject-specific CVR and haemodynamic lag estimates. The evaluated approaches range from conventional regression models, i.e. including drifts and motion timecourses as nuisance regressors, applied on single-echo or optimally-combined data, to more complex models including regressors obtained from multi-echo independent component analysis with different grades of orthogonalization in order to preserve the effect of interest, i.e. the CVR. We compare these models in terms of their ability to make signal intensity changes independent from motion, as well as the reliability as measured by voxelwise intraclass correlation coefficients of both CVR and lag maps over time. Our results reveal that a conservative independent component analysis model applied on the optimally-combined multi-echo fMRI signal offers the largest reduction of motion-related effects in the signal, while yielding reliable CVR amplitude and lag estimates, although a conventional regression model applied on the optimally-combined data results in similar estimates. This work demonstrates the usefulness of multi-echo based fMRI acquisitions and independent component analysis denoising for precision mapping of CVR in single subjects based on BH paradigms, fostering its potential as a clinically-viable neuroimaging tool for individual patients. It also proves that the way in which data-driven regressors should be incorporated in the analysis model is not straight-forward due to their complex interaction with the BH-induced BOLD response.

Cerebrovascular Reactivity Measurement Using Magnetic Resonance Imaging: A Systematic Review

Cerebrovascular reactivity (CVR) magnetic resonance imaging (MRI) probes cerebral haemodynamic changes in response to a vasodilatory stimulus. CVR closely relates to the health of the vasculature and is therefore a key parameter for studying cerebrovascular diseases such as stroke, small vessel disease and dementias. MRI allows in vivo measurement of CVR but several different methods have been presented in the literature, differing in pulse sequence, hardware requirements, stimulus and image processing technique. We systematically reviewed publications measuring CVR using MRI up to June 2020, identifying 235 relevant papers. We summarised the acquisition methods, experimental parameters, hardware and CVR quantification approaches used, clinical populations investigated, and corresponding summary CVR measures. CVR was investigated in many pathologies such as steno-occlusive diseases, dementia and small vessel disease and is generally lower in patients than in healthy controls. Blood oxygen level dependent (BOLD) acquisitions with fixed inspired CO₂ gas or end-tidal CO₂ forcing stimulus are the most commonly used methods. General linear modelling of the MRI signal with end-tidal CO₂ as the regressor is the most frequently used method to compute CVR. Our survey of CVR measurement approaches and applications will help researchers to identify good practice and provide objective information to inform the development of future consensus recommendations.

Altered cerebrovascular reactivity due to respiratory rate and breath holding: a BOLD-fMRI study on healthy adults

Cerebrovascular reactivity (CVR) is of great significance for the treatment and prevention of cerebrovascular diseases. CVR can be mapped using the blood oxygenation level-dependent (BOLD) signal of fMRI. Breath holding (BH) is a reliable method to produce the desired increase in arterial CO₂, while its application in clinical research is limited due to subject's compliance and variability. BH task with variable respiratory rates could allow more flexibility in clinical populations. In this study, 50 healthy volunteers were scanned for end-inspiration BH tasks with three different respiration rates. For the three respiratory rates BH tasks, the CVR was estimated based on the BOLD signal and general linear model (GLM) separately. Specifically, the extra time delay was considered for the hemodynamic response function, and the optimal delay was estimated for each voxel. To measure CVR in grey matter, BOLD signals of end-inspiration BH were used as regressors in general linear models to quantify their impact on CVR. This was performed for regions and voxels. Systematic differences were observed between the three end-inspiratory breathing rates. The greatest increase in activation intensity was found in fast breathing followed by self-paced and slow breathing. We conclude that the BH task of variable respiratory rates allows for CVR measurement, making breath-holding challenges more flexible and appropriate for routine practice.

Cerebrovascular Reactivity Mapping Without Gas Challenges: A Methodological Guide

Cerebrovascular reactivity (CVR) is defined as the ability of vessels to alter their caliber in response to vasoactive factors, by means of dilating or constricting, in order to increase or decrease regional cerebral blood flow (CBF). Importantly, CVR may provide a sensitive biomarker for pathologies where vasculature is compromised. Furthermore, the spatiotemporal dynamics of CVR observed in healthy subjects, reflecting regional differences in cerebral vascular tone and response, may also be important in functional MRI studies based on neurovascular coupling mechanisms. Assessment of CVR is usually based on the use of a vasoactive stimulus combined with a CBF measurement technique. Although transcranial Doppler ultrasound has been frequently used to obtain global flow velocity measurements, MRI techniques are being increasingly employed for obtaining CBF maps. For the vasoactive stimulus, vasodilatory hypercapnia is usually induced through the manipulation of respiratory gases, including the inhalation of increased concentrations of carbon dioxide. However, most of these methods require an additional apparatus and complex setups, which not only may not be well-tolerated by some populations but are also not widely available. For these reasons, strategies based on voluntary breathing fluctuations without the need for external gas challenges have been proposed. These include the task-based methodologies of breath holding and paced deep breathing, as well as a new generation of methods based on spontaneous breathing fluctuations during resting-state. Despite the multitude of alternatives to gas challenges, existing literature lacks definitive conclusions regarding the best practices for the vasoactive modulation and associated analysis protocols. In this work, we perform an extensive review of CVR mapping techniques based on MRI and CO2 variations without gas challenges, focusing on the methodological aspects of the breathing protocols and corresponding data analysis. Finally, we outline a set of practical guidelines based on generally accepted practices and available data, extending previous reports and encouraging the wider application of CVR mapping methodologies in both clinical and academic MRI settings.

Voxelwise optimization of hemodynamic lags to improve regional CVR estimates in breath-hold fMRI

Cerebrovascular Reactivity (CVR), the responsiveness of blood vessels to a vasodilatory stimulus, is an important indicator of cerebrovascular health. Assessing CVR with fMRI, we can measure the change in the Blood Oxygen Level Dependent (BOLD) response induced by a change in CO2 pressure (%BOLD/mmHg). However, there exists a temporal offset between the recorded CO2 pressure and the local BOLD response, due to both measurement and physiological delays. If this offset is not corrected for, voxel-wise CVR values will not be accurate. In this paper, we propose a framework for mapping hemodynamic lag in breath-hold fMRI data. As breath-hold tasks drive task-correlated head motion artifacts in BOLD fMRI data, our framework for lag estimation fits a model that includes polynomial terms and head motion parameters, as well as a shifted variant of the CO2 regressor (±9 s in 0.3 s increments), and the hemodynamic lag at each voxel is the shift producing the maximum total model R2 within physiological constraints. This approach is evaluated in 8 subjects with multi-echo fMRI data, resulting in robust maps of hemodynamic delay that show consistent regional variation across subjects, and improved contrast-to-noise compared to methods where motion regression is ignored or performed earlier in preprocessing.Clinical Relevance- We map hemodynamic lag using breathhold fMRI, providing insight into vascular transit times and improving the regional accuracy of cerebrovascular reactivity measurements.

Measuring vascular reactivity with resting-state blood oxygenation level-dependent (BOLD) signal fluctuations: A potential alternative to the breath-holding challenge?

Measurement of the ability of blood vessels to dilate and constrict, known as vascular reactivity, is often performed with breath-holding tasks that transiently raise arterial blood carbon dioxide (P_aCO₂) levels. However, following the proper commands for a breath-holding experiment may be difficult or impossible for many patients. In this study, we evaluated two approaches for obtaining vascular reactivity information using blood oxygenation level-dependent signal fluctuations obtained from resting-state functional magnetic resonance imaging data: physiological fluctuation regression and coefficient of variation of the resting-state functional magnetic resonance imaging signal. We studied a cohort of 28 older adults (69 ± 7 years) and found that six of them (21%) could not perform the breath-holding protocol, based on an objective comparison with an idealized respiratory waveform. In the subjects that could comply, we found a strong linear correlation between data extracted from spontaneous resting-state functional magnetic resonance imaging signal fluctuations and the blood oxygenation level-dependent percentage signal change during breath-holding challenge ( R²= 0.57 and 0.61 for resting-state physiological fluctuation regression and resting-state coefficient of variation methods, respectively). This technique may eliminate the need for subject cooperation, thus allowing the evaluation of vascular reactivity in a wider range of clinical and research conditions in which it may otherwise be impractical.

Fourier modeling of the BOLD response to a breath-hold task: Optimization and reproducibility

Cerebrovascular reactivity (CVR) reflects the capacity of blood vessels to adjust their caliber in order to maintain a steady supply of brain perfusion, and it may provide a sensitive disease biomarker. Measurement of the blood oxygen level dependent (BOLD) response to a hypercapnia-inducing breath-hold (BH) task has been frequently used to map CVR noninvasively using functional magnetic resonance imaging (fMRI). However, the best modeling approach for the accurate quantification of CVR maps remains an open issue. Here, we compare and optimize Fourier models of the BOLD response to a BH task with a preparatory inspiration, and assess the test-retest reproducibility of the associated CVR measurements, in a group of 10 healthy volunteers studied over two fMRI sessions. Linear combinations of sine-cosine pairs at the BH task frequency and its successive harmonics were added sequentially in a nested models approach, and were compared in terms of the adjusted coefficient of determination and corresponding variance explained (VE) of the BOLD signal, as well as the number of voxels exhibiting significant BOLD responses, the estimated CVR values, and their test-retest reproducibility. The brain average VE increased significantly with the Fourier model order, up to the 3rd order. However, the number of responsive voxels increased significantly only up to the 2nd order, and started to decrease from the 3rd order onwards. Moreover, no significant relative underestimation of CVR values was observed beyond the 2nd order. Hence, the 2nd order model was concluded to be the optimal choice for the studied paradigm. This model also yielded the best test-retest reproducibility results, with intra-subject coefficients of variation of 12 and 16% and an intra-class correlation coefficient of 0.74. In conclusion, our results indicate that a Fourier series set consisting of a sine-cosine pair at the BH task frequency and its two harmonics is a suitable model for BOLD-fMRI CVR measurements based on a BH task with preparatory inspiration, yielding robust estimates of this important physiological parameter.

Reproducibility of cerebrovascular reactivity measures in children using BOLD MRI

PURPOSE

To evaluate the reproducibility of cerebrovascular reactivity (CVR) measurements acquired in children using magnetic resonance imaging (MRI) in combination with a computer-controlled carbon dioxide (CO2 ) stimulus.

MATERIALS AND METHODS

Ten healthy children (age 16.1 ± 1.6 years) underwent CVR imaging on a 3T scanner using a blood-oxygen level-dependent (BOLD) MRI sequence. Targeted hypercapnia was induced during imaging with a CO2 gas challenge delivered using a specialized gas sequencer (RespirAct). A total of four BOLD scans were performed over 2 separate days to test within-day and between-day consistency of the data. CVR values were computed by correlating the relative change in BOLD signal in response to the CO2 stimulus delivered to the each subject.

RESULTS

Intraclass correlation coefficients (ICCs) of within-day values show highly reproducible measures in both the gray matter (ICC = 0.857, P < 0.001) and white matter (ICC = 0.895, P < 0.001). Relatively lower between-day reproducibility was observed in both the gray matter (ICC = 0.776, P = 0.001) and white matter (ICC = 0.719, P = 0.004).

CONCLUSION

Using a computer-controlled CO2 stimulus, we have demonstrated the reliability of BOLD-CVR measurements in pediatric subjects. Within-day and between-day metrics of reproducibility were comparable to adult data.

Comparing cerebrovascular reactivity measured using BOLD and cerebral blood flow MRI: The effect of basal vascular tension on vasodilatory and vasoconstrictive reactivity

Cerebrovascular reactivity (CVR) is an important metric of cerebrovascular health. While the BOLD fMRI method in conjunction with carbon-dioxide (CO2) based vascular manipulation has been the most commonly used, the BOLD signal is not a direct measure of vascular changes, and the use of arterial-spin labeling (ASL) cerebral blood flow (CBF) imaging is increasingly advocated. Nonetheless, given the differing dependencies of BOLD and CBF on vascular baseline conditions and the diverse CO2 manipulation types currently used in the literature, knowledge of potential biases introduced by each technique is critical for the interpretation of CVR measurements. In this work, we use simultaneous BOLD-CBF acquisitions during both vasodilatory (hypercapnic) and vasoconstrictive (hypocapnic) stimuli to measure CVR. We further imposed different levels of baseline vascular tension by inducing hypercapnic and hypocapnic baselines, separately from normocapnia by 4mmHg. We saw significant and diverse dependencies on vascular stimulus and baseline condition in both BOLD and CBF CVR measurements: (i) BOLD-based CVR is more sensitive to basal vascular tension than CBF-based CVR; (ii) the use of a combination of vasodilatory and vasoconstrictive stimuli maximizes the sensitivity of CBF-based CVR to vascular tension changes; (iii) the BOLD and CBF vascular response delays are both significantly lengthened at predilated baseline. As vascular tension can often be altered by potential pathology, our findings are important considerations when interpreting CVR measurements in health and disease.