Fine tuning breath-hold-based cerebrovascular reactivity analysis models

Long-term effects of multiple concussions on prefrontal cortex oxygenation during a hypercapnic challenge in retired contact sport athletes

This exploratory study aimed to investigate the long-term effects of multiple concussions on prefrontal cortex oxygenation during a five-minute hypercapnic challenge using Near Infrared Spectroscopy (NIRS). 55 physically active retired contact sport male athletes with three or more previous concussions (mTBI) were recruited along with 29 physically active males with no concussions history (CTRL). Participants completed five minutes of seated rest prior to the five-minute hypercapnic challenge (20-second breath-hold, 40-second recovery breathing; five times). NIRS measured right and left side oxygenated (O2Hb), deoxygenated (HHb), total (tHb) haemoglobin, and haemoglobin difference (HbDiff) with all parameters analysed through changes in average maximal and minimal values (ΔMAX), Z-scores, and standard deviations. Right prefrontal cortex HbDiff ΔMAX was significantly higher in the mTBI compared to CTRL (p = 0.045) group. Left prefrontal cortex O2Hb ΔMAX (p = 0.040), HHb Z-Scores (p = 0.008), and HbDiff ΔMAX(p = 0.014) were significantly higher in the mTBI group. Within-group right vs left analyses demonstrated significantly lower left HbDiff ΔMAX (p = 0.048) and HbDiff Z-scores (p = 0.002) in the mTBI group, while the CTRL group had significantly lower left HHb Z-scores (p = 0.003) and left tHb Z-scores (p = 0.042). This study provides preliminary evidence that athletes with a history of three or more concussions may have impaired prefrontal cortex oxygenation parameters during a hypercapnic challenge.

Evaluation of the cerebrovascular reactivity in patients with Moyamoya Angiopathy by use of breath-hold fMRI: investigation of voxel-wise hemodynamic delay correction in comparison to [15O]water PET

PURPOSE

Patients with Moyamoya Angiopathy (MMA) require hemodynamic assessment to evaluate the risk of stroke. Hemodynamic evaluation by use of breath-hold-triggered fMRI (bh-fMRI) was proposed as a readily available alternative to the diagnostic standard [¹⁵O]water PET. Recent studies suggest voxel-wise hemodynamic delay correction in hypercapnia-triggered fMRI. The aim of this study was to evaluate the effect of delay correction of bh-fMRI in patients with MMA and to compare the results with [¹⁵O]water PET.

METHODS

bh-fMRI data sets of 22 patients with MMA were evaluated without and with voxel-wise delay correction within different shift ranges and compared to the corresponding [¹⁵O]water PET data sets. The effects were evaluated combined and in subgroups of data sets with most severely impaired CVR (apparent steal phenomenon), data sets with territorial time delay, and data sets with neither steal phenomenon nor delay between vascular territories.

RESULTS

The study revealed a high mean cross-correlation (r = 0.79, p < 0.001) between bh-fMRI and [¹⁵O]water PET. The correlation was strongly dependent on the choice of the shift range. Overall, no shift range revealed a significantly improved correlation between bh-fMRI and [¹⁵O]water PET compared to the correlation without delay correction. Delay correction within shift ranges with positive high high cutoff revealed a lower agreement between bh-fMRI and PET overall and in all subgroups.

CONCLUSION

Voxel-wise delay correction, in particular with shift ranges with high cutoff, should be used critically as it can lead to false-negative results in regions with impaired CVR and a lower correlation to the diagnostic standard [¹⁵O]water PET.

Hemodynamic evaluation of patients with Moyamoya Angiopathy: comparison of resting-state fMRI to breath-hold fMRI and [15O]water PET

Purpose

Patients with Moyamoya Angiopathy (MMA) require hemodynamic evaluation to assess the risk of stroke. Assessment of cerebral blood flow with [¹⁵O]water PET and acetazolamide challenge is the diagnostic standard for the evaluation of the cerebral perfusion reserve (CPR). Estimation of the cerebrovascular reactivity (CVR) by use of breath-hold-triggered fMRI (bh-fMRI) as an index of CPR has been proposed as a reliable and more readily available approach. Recent findings suggest the use of resting-state fMRI (rs-fMRI) which requires minimum patient compliance. The aim of this study was to compare rs-fMRI to bh-fMRI and [¹⁵O]water PET in patients with MMA.

Methods

Patients with MMA underwent rs-fMRI and bh-fMRI in the same MRI session. Maps of the CVR gained by both modalities were compared retrospectively by calculating the correlation between the mean CVR of 12 volumes of interest. Additionally, the rs-maps of a subgroup of patients were compared to CPR-maps gained by [¹⁵O]water PET.

Results

The comparison of the rs-maps and the bh-maps of 24 patients revealed a good correlation (Pearson’s r = 0.71 ± 0.13; preoperative patients: Pearson’s r = 0.71 ± 0.17; postoperative patients: Pearson’s r = 0.71 ± 0.11). The comparison of 7 rs-fMRI data sets to the corresponding [¹⁵O]water PET data sets also revealed a high level of agreement (Pearson’s r = 0.80 ± 0.19).

Conclusion

The present analysis indicates that rs-fMRI might be a promising non-invasive method with almost no patient cooperation needed to evaluate the CVR. Further prospective studies are required.

Effects of Thresholding on Voxel-Wise Correspondence of Breath-Hold and Resting-State Maps of Cerebrovascular Reactivity

Functional magnetic resonance imaging for presurgical brain mapping enables neurosurgeons to identify viable tissue near a site of operable pathology which might be at risk of surgery-induced damage. However, focal brain pathology (e.g., tumors) may selectively disrupt neurovascular coupling while leaving the underlying neurons functionally intact. Such neurovascular uncoupling can result in false negatives on brain activation maps thereby compromising their use for surgical planning. One way to detect potential neurovascular uncoupling is to map cerebrovascular reactivity using either an active breath-hold challenge or a passive resting-state scan. The equivalence of these two methods has yet to be fully established, especially at a voxel level of resolution. To quantitatively compare breath-hold and resting-state maps of cerebrovascular reactivity, we first identified threshold settings that optimized coverage of gray matter while minimizing false responses in white matter. When so optimized, the resting-state metric had moderately better gray matter coverage and specificity. We then assessed the spatial correspondence between the two metrics within cortical gray matter, again, across a wide range of thresholds. Optimal spatial correspondence was strongly dependent on threshold settings which if improperly set tended to produce statistically biased maps. When optimized, the two CVR maps did have moderately good correspondence with each other (mean accuracy of 73.6%). Our results show that while the breath-hold and resting-state maps may appear qualitatively similar they are not quantitatively identical at a voxel level of resolution.

Intraoperative BOLD-fMRI Cerebrovascular Reactivity Assessment

Blood oxygenation-level dependent cerebrovascular reactivity (BOLD-CVR) has gained attention in recent years as an effective way to investigate CVR, a measure of the hemodynamic state of the brain, with high spatial and temporal resolution. An association between impaired CVR and diverse pathologies has been observed, especially in ischemic cerebrovascular diseases and brain gliomas. The ability to obtain this information intraoperatively is novel and has not been widely tested. We report our first experience with this intraoperative technique in vascular and oncologic neurosurgical patients, discuss the results of its feasibility, and the possible developments of the intraoperative employment of BOLD-CVR.

Improving the Breath-Holding CVR Measurement Using the Multiband Multi-Echo EPI Sequence

Blood oxygen level-dependent (BOLD) functional MRI (fMRI) is commonly used to measure cerebrovascular reactivity (CVR), which can convey insightful information about neurovascular health. Breath-holding (BH) has been shown to be a practical vasodilatory stimulus for measuring CVR in clinical settings. The conventional BOLD fMRI approach has some limitations, however, such as susceptibility-induced signal dropout at air tissue interfaces and low BOLD sensitivity especially in areas of low T 2 * . These drawbacks can potentially be mitigated with multi-echo sequences, which acquire several images at different echo times in one shot. When combined with multiband techniques, high temporal resolution images can be acquired. This study compared an advanced multiband multi-echo (MBME) echo planar imaging (EPI) sequence with an existing multiband single-echo (MB) sequence to evaluate the repeatability and sensitivity of BH activation and CVR mapping. Images were acquired from 28 healthy volunteers, of which 18 returned for repeat imaging. Both MBME and MB data were pre-processed using both standard and advanced denoising techniques. The MBME data was further processed by combining echoes using a T 2 * -weighted approach and denoising using multi-echo independent component analysis. BH activation was calculated using a general linear model and the respiration response function. CVR was computed as the percent change related to the activation. To account for differences in CVR related to TE, relative CVR (rCVR) was computed and normalized to the mean gray matter CVR. Test-retest metrics were assessed with the Dice coefficient, rCVR difference, within subject coefficient of variation, and the intraclass correlation coefficient. Our findings demonstrate that rCVR for MBME scans were significantly higher than for MB scans across most of the gray matter. In areas of high susceptibility-induced signal dropout, however, MBME rCVR was significantly less than MB rCVR due to artifactually high rCVR for MB scans in these regions. MBME rCVR showed improved test-retest metrics compared with MB. Overall, the MBME sequence displayed superior BOLD sensitivity, improved specificity in areas of signal dropout on MBME scans, enhanced reliability, and reduced variability across subjects compared with MB acquisitions. Our results suggest that the MBME EPI sequence is a promising tool for imaging CVR.

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.

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.

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'.

Feasibility and safety of intraoperative BOLD functional MRI cerebrovascular reactivity to evaluate extracranial-to-intracranial bypass efficacy

Blood oxygenation level-dependent functional MRI cerebrovascular reactivity (BOLD-CVR) is a contemporary technique to assess brain tissue hemodynamic changes after extracranial- intracranial (EC-IC) bypass flow augmentation surgery. The authors conducted a preliminary study to investigate the feasibility and safety of intraoperative 3-T MRI BOLD-CVR after EC-IC bypass flow augmentation surgery. Five consecutive patients selected for EC-IC bypass revascularization underwent an intraoperative BOLD-CVR examination to assess early hemodynamic changes after revascularization and to confirm the safety of this technique. All patients had a normal postoperative course, and none of the patients exhibited complications or radiological alterations related to prolonged anesthesia time. In addition to intraoperative flow measurements of the bypass graft, BOLD-CVR maps added information on the hemodynamic status and changes at the brain tissue level. Intraoperative BOLD-CVR is feasible and safe in patients undergoing EC-IC bypass revascularization. This technique can offer immediate hemodynamic feedback on brain tissue revascularization after bypass flow augmentation surgery.

Co-localized impaired regional cerebrovascular reactivity in chronic concussion is associated with BOLD activation differences during a working memory task

The purpose of this study was to quantify differences in blood oxygen level dependent (BOLD) activation on a working memory task, baseline cerebral blood flow (CBF₀), and cerebrovascular reactivity (CVR) between participants with and without a history of concussion. A dual-echo pseudo-continuous arterial spin labelling (pCASL) sequence was performed on a group of 10 subjects with a previous concussion (126 ± 15 days prior) and on a control group (n = 10) during a visual working memory protocol. A separate dual-echo pCASL sequence was used to derive CVR and CBF₀ measurements from a boxcar hypercapnic breathing protocol. Brain areas with significant activation differences on the working memory task between groups were identified and combined as an aggregate region of interest for CBF and CVR analyses. Areas of reduced BOLD activation during the working memory task in the concussed group included the ventral anterior cingulate cortex (ACC), the medial temporal gyrus (MTG), and the lateral occipital cortex in two loci. A single area of increased activation was located in the parietal operculum. Further analyses of CBF₀ and CVR in these regions revealed reduced CVR in the concussed group in the MTG and ACC, while CBF₀ did not differ. The differences in CVR between the two groups in these regions suggest that concussive injury may result in microvascular dysfunction. In turn, the decreased BOLD response during the task could be due to altered neurovascular coupling, rather than an impairment in neural activation alone. However, in other regions associated with working memory, unchanged CBF₀ and CVR suggests that neural injury also persists after concussion. In the future, BOLD results should be normalized to CVR in order achieve a clearer understanding of the neural and vascular contributions to the differences in the signal.

Improving the Assessment of Breath-Holding Induced Cerebral Vascular Reactivity Using a Multiband Multi-echo ASL/BOLD Sequence

Breath holding (BH) is a viable vasodilatory stimulus for calculating functional MRI-derived cerebral vascular reactivity (CVR). The BH technique suffers from reduced repeatability compared with gas inhalation techniques; however, extra equipment is needed to perform gas inhalation techniques, and this equipment is not available at all institutions. This study aimed to determine the sensitivity and repeatability of BH activation and CVR using a multiband multi-echo simultaneous arterial spin labelling/blood oxygenation level dependent (ASL/BOLD) sequence. Whole-brain images were acquired in 14 volunteers. Ten subjects returned for repeat imaging. Each subject performed four cycles of 16 s BH on expiration interleaved with paced breathing. Following standard preprocessing, the echoes were combined using a T2*-weighted approach. BOLD and ASL BH activation was computed, and CVR was then determined as the percent signal change related to the activation. The "M" parameter from the Davis Model was also computed by incorporating the ASL signal. Our results showed higher BH activation strength, volume, and repeatability for the combined multi-echo (MEC) data compared with the single-echo data. MEC CVR also had higher repeatability, sensitivity, specificity, and reliability compared with the single-echo BOLD data. These data support the usefulness of an MBME ASL/BOLD acquisition for BH CVR and M measurements.

Hypercapnic BOLD MRI compared to H215O PET/CT for the hemodynamic evaluation of patients with Moyamoya Disease

Background and purpose

Patients with Moyamoya Disease (MMD) need hemodynamic evaluation of vascular territories at risk of stroke. Today's investigative standards include H₂¹⁵O PET/CT with pharmacological challenges with acetazolamide (ACZ). Recent developments suggest that CO₂-triggered blood‑oxygen-level-dependent (BOLD) functional MRI might provide comparable results to current standard methods for evaluation of territorial hemodynamics, while being a more widely available and easily implementable method. This study examines results of a newly developed quantifiable analysis algorithm for CO₂-triggered BOLD MRI in Moyamoya patients and correlates the results with H₂¹⁵O PET/CT with ACZ challenge to assess comparability between both modalities.

Methods

CO₂-triggered BOLD MRI was performed and compared to H₂¹⁵O PET/CT with ACZ challenge in patients with angiographically proven MMD. Images of both modalities were analyzed retrospectively in a blinded, standardized fashion by visual inspection, as well as with a semi-quantitative analysis using stimuli-induced approximated regional perfusion-weighted data and BOLD-signal changes with reference to cerebellum.

Results

20 consecutive patients fulfilled the inclusion criteria, a total of 160 vascular territories were analyzed retrospectively. Visual analysis (4-step visual rating system) of standardized, color-coded cerebrovascular reserve/reactivity maps showed a very strong correlation (Spearman's rho = 0.9, P < 0.001) between both modalities. Likewise, comparison of approximated regional perfusion changes across vascular territories (normalized to cerebellar change) reveal a highly significant correlation between both methods (Pearson's r = 0.71, P < 0.001).

Conclusions

The present analysis indicates that CO₂-triggered BOLD MRI is a very promising tool for the hemodynamic evaluation of MMD patients with results comparable to those seen in H₂¹⁵O PET/CT with ACZ challenge. It therefore holds future potential in becoming a routine examination in the pre- and postoperative evaluation of MMD patients after further prospective evaluation.

•

Non-invasive cerebrovascular reactivity measurement with BOLD MRI.

•

CO₂-triggered BOLD MRI correlates strongly with H₂¹⁵O PET/CT with ACZ challenge

•

Widely-available tool for the hemodynamic evaluation of Moyamoya patients.

BOLD cerebrovascular reactivity as a novel marker for crossed cerebellar diaschisis

Objective

To study blood oxygen level–dependent cerebrovascular reactivity (BOLD-CVR) as a surrogate imaging marker for crossed cerebellar diaschisis (CCD).

Methods

Twenty-five participants with symptomatic unilateral cerebrovascular steno-occlusive disease underwent a BOLD-CVR and an acetazolamide challenged (15O)-H2O-PET study. CCD and cerebellar asymmetry index were determined from PET and compared to BOLD-CVR quantitative values. Neurologic status at admission and outcome after 3 months were determined with NIH Stroke Scale (NIHSS) and modified Rankin Scale (mRS) scores.

Results

For both the BOLD-CVR and PET examination, a significant cerebellar asymmetry index was found for participants exhibiting CCD (CCD+ vs CCD−: for BOLD-CVR 13.11 ± 9.46 vs 1.52 ± 4.97, p < 0.001; and for PET 7.31 ± 2.75 vs 1.68 ± 2.98, p < 0.001). The area under the curve for BOLD-CVR was 0.89 (95% confidence interval: 0.75–1.0) with 0.91 sensitivity and 0.81 specificity to detect CCD. Participants exhibiting CCD were in poorer clinical condition at baseline (CCD+ vs CCD−: NIHSS 7 vs 1, p = 0.003; mRS 3 vs 1, p = 0.001) and after 3-month follow-up (NIHSS 2 vs 0, p = 0.02; mRS 1 vs 0, p = 0.04). Worse performance on both scores showed an agreement with a larger BOLD-CVR cerebellar asymmetry index. This was not found for PET.

Conclusions

BOLD-CVR demonstrates similar sensitivity to detect CCD as compared to (15O)-H2O-PET in patients with symptomatic unilateral cerebrovascular steno-occlusive disease. Furthermore, participants exhibiting CCD had a poorer baseline neurologic performance and neurologic outcome at 3 months.

Classification of evidence

This study provides Class II evidence that BOLD-CVR identifies CCD in patients with symptomatic unilateral cerebrovascular steno-occlusive disease.

A Resilient, Non-neuronal Source of the Spatiotemporal Lag Structure Detected by BOLD Signal-Based Blood Flow Tracking

Recent evidence has suggested that blood oxygenation level-dependent (BOLD) signals convey information about brain circulation via low frequency oscillation of systemic origin (sLFO) that travels through the vascular structure ("lag mapping"). Prompted by its promising application in both physiology and pathology, we examined this signal component using multiple approaches. A total of 30 healthy volunteers were recruited to perform two reproducibility experiments at 3 Tesla using multiband echo planar imaging. The first experiment investigated the effect of denoising and the second was designed to study the effect of subject behavior on lag mapping. The lag map's intersession test-retest reproducibility and image contrast were both diminished by removal of either the neuronal or the non-neuronal (e.g., cardiac, respiratory) components by independent component analysis-based denoising, suggesting that the neurovascular coupling also comprises a part of the BOLD lag structure. The lag maps were, at the same time, robust against local perfusion increases due to visuomotor task and global changes in perfusion induced by breath-holding at the same level as the intrasession reliability. The lag structure was preserved after time-locked averaging to the visuomotor task and breath-holding events, while any preceding signal changes were canceled out for the visuomotor task, consistent with the passive effect of neurovascular coupling in the venous side of the vasculature. These findings support the current assumption that lag mapping primarily reflects vascular structure despite the presence of sLFO perturbation of neuronal or non-neuronal origin and, thus, emphasize the vascular origin of the lag map, encouraging application of BOLD-based blood flow tracking.

Fine tuning breath-hold-based cerebrovascular reactivity analysis models

INTRODUCTION

We elaborate on existing analysis methods for breath-hold (BH)-derived cerebrovascular reactivity (CVR) measurements and describe novel insights and models toward more exact CVR interpretation.

METHODS

Five blood-oxygen-level-dependent (BOLD) fMRI datasets of neurovascular patients with unilateral hemispheric hemodynamic impairment were used to test various BH CVR analysis methods. Temporal lag (phase), percent BOLD signal change (CVR), and explained variance (coherence) maps were calculated using three different sine models and two novel "Optimal Signal" model-free methods based on the unaffected hemisphere and the sagittal sinus fMRI signal time series, respectively.

RESULTS

All models showed significant differences in CVR and coherence between the affected-hemodynamic impaired-and unaffected hemisphere. Voxel-wise phase determination significantly increases CVR (0.60 ± 0.18 vs. 0.82 ± 0.27; P < 0.05). Incorporating different durations of breath hold and resting period in one sine model (two-task) did increase coherence in the unaffected hemisphere, as well as eliminating negative phase commonly obtained by one-task frequency models. The novel model-free "optimal signal" methods both explained the BOLD MR data similar to the two task sine model.

CONCLUSIONS

Our CVR analysis demonstrates an improved CVR and coherence after implementation of voxel-wise phase and frequency adjustment. The novel "optimal signal" methods provide a robust and feasible alternative to the sine models, as both are model-free and independent of compliance. Here, the sagittal sinus model may be advantageous, as it is independent of hemispheric CVR impairment.