Cerebrovascular reactivity mapping using intermittent breath modulation

Non-invasive in-utero quantification of vascular reactivity in human placenta

OBJECTIVE

Placental vascular reactivity (PlVR) indicates the ability of the placental vasculature to match blood supply to fetal demand. Many pregnancy disorders alter the characteristics of PlVR, resulting in suboptimal oxygen delivery, although current understanding is limited by the lack of non-invasive, repeatable methods to measure PlVR in utero. Our objective was to quantify PlVR by measuring the placental response to transient changes in maternal carbon dioxide (CO2) using blood-oxygen-level-dependent (BOLD) magnetic resonance imaging (MRI). We hypothesized that PlVR will increase with gestational age to meet the changing demands of a growing fetus, and that PlVR will be driven by a maternal response to changes in CO2 concentration.

METHODS

This was a cross-sectional study of 35 women with a healthy singleton pregnancy, of whom 31 were included in the analysis. The median gestational age was 32.6 (range, 22.6-38.4) weeks. Pregnant women were instructed to follow audiovisual breathing cues during a MRI scan. Maternal end-tidal CO2 (EtCO2) was measured concurrently with resting placental BOLD MRI for a total of 7-8 min. Preprocessing of magnetic resonance images consisted of manual delineation of placental anatomy and motion correction. In each placental voxel, vascular reactivity was computed using a coherence-weighted general linear model between MRI signal and EtCO2 stimulus. Global PlVR was computed as the mean of voxel-wise PlVR values across the placenta.

RESULTS

PlVR, quantified by the placental response to induced, transient changes in maternal CO2, was consistently measured in utero using BOLD MRI. PlVR increased non-linearly with advancing gestational age (P < 0.001) and was higher on the fetal side of the placenta. PlVR was associated positively with fetal brain volume after accounting for gestational age. PlVR did not show any significant associations with maternal characteristics.

CONCLUSIONS

We present, for the first time, a non-invasive paradigm to quantify PlVR in ongoing human pregnancies without the use of exogenous gases or contrast agents. Our findings suggest that PlVR is driven by a fetal response to changes in maternal CO2. Ease of translation to the clinical setting makes PlVR a promising biomarker for the identification and management of high-risk pregnancies. © 2023 The Authors. Ultrasound in Obstetrics & Gynecology published by John Wiley & Sons Ltd on behalf of International Society of Ultrasound in Obstetrics and Gynecology.

Improved cerebrovascular reactivity mapping using coherence weighted general linear model in the frequency domain

Cerebrovascular reactivity (CVR) is a prognostic indicator of cerebrovascular health. Estimating CVR from endogenous end-tidal carbon dioxide (CO2) fluctuation and MRI signal recorded under resting state can be difficult due to the poor signal-to-noise ratio (SNR) of signals. Thus, we aimed to improve the method of estimating CVR from end-tidal CO2 and MRI signals. We proposed a coherence weighted general linear model (CW-GLM) to estimate CVR from the Fourier coefficients weighted by the signal coherence in frequency domain, which confers two advantages. First, it requires no signal alignment in time domain, which simplifies experimental methods. Second, it limits the GLM analysis within the frequency band where CO2 and MRI signals are highly correlated, which automatically suppresses noise and nuisance signals. We compared the performance of our method with time-domain GLM (TD-GLM) and frequency-domain GLM (FD-GLM) in both synthetic and in-vivo data; wherein we calculated CVR from signals recorded under both resting state and sinusoidal stimulus. In synthetic data, CW-GLM has a remarkable performance on CVR estimation from narrow band signals with a mean-absolute error of 0.7 % (gray matter) and 1.2 % (white matter), which was lower than all the other methods. Meanwhile, CW-GLM maintains a comparable performance on CVR estimation from resting signals, with a mean-absolute error of 4.1 % (gray matter) and 8 % (white matter). The superior performance was maintained across the 36 in-vivo measurements, with CW-GLM exhibiting limits of agreement of -16.7 % - 9.5 % between CVR calculated from the resting and sinusoidal CO2 paradigms which was 12 % - 209 % better than current time-domain methods. Evaluating of the cross-coherence spectrum revealed highest signal coherence within the frequency band from 0.01 Hz to 0.05 Hz, which overlaps with previously recommended frequency band (0.02 Hz to 0.04 Hz) for CVR analysis. Our data demonstrates that CW-GLM can work as a self-adaptive band-pass filter to improve CVR robustness, while also avoiding the need for signal temporal alignment.

Assessing depth sensitivity in laser interferometry speckle visibility spectroscopy (iSVS) through source-to-detector distance variation and cerebral blood flow monitoring in humans and rabbits

Recently, speckle visibility spectroscopy (SVS) was non-invasively applied on the head to monitor cerebral blood flow. The technique, using a multi-pixel detecting device (e.g., camera), allows the detection of a larger number of speckles, increasing the proportion of light that is detected. Due to this increase, it is possible to collect light that has propagated deeper through the brain. As a direct consequence, cerebral blood flow can be monitored. However, isolating the cerebral blood flow from the other layers, such as the scalp or skull components, remains challenging. In this paper, we report our investigations on the depth-sensitivity of laser interferometry speckle visibility spectroscopy (iSVS). Specifically, we varied the depth of penetration of the laser light into the head by tuning the source-to-detector distance, and identified the transition point at which cerebral blood flow in humans and rabbits starts to be detected.

Agreement in cerebrovascular reactivity assessed with diffuse correlation spectroscopy across experimental paradigms improves with short separation regression

Significance

Cerebrovascular reactivity (CVR), i.e., the ability of cerebral vasculature to dilate or constrict in response to vasoactive stimuli, is a biomarker of vascular health. Exogenous administration of inhaled carbon dioxide, i.e., hypercapnia (HC), remains the "gold-standard" intervention to assess CVR. More tolerable paradigms that enable CVR quantification when HC is difficult/contraindicated have been proposed. However, because these paradigms feature mechanistic differences in action, an assessment of agreement of these more tolerable paradigms to HC is needed.

Aim

We aim to determine the agreement of CVR assessed during HC, breath-hold (BH), and resting state (RS) paradigms.

Approach

Healthy adults were subject to HC, BH, and RS paradigms. End tidal carbon dioxide (EtCO2) and cerebral blood flow (CBF, assessed with diffuse correlation spectroscopy) were monitored continuously. CVR (%/mmHg) was quantified via linear regression of CBF versus EtCO2 or via a general linear model (GLM) that was used to minimize the influence of systemic and extracerebral signal contributions.

Results

Strong agreement ( CCC ≥ 0.69 ; R ≥ 0.76 ) among CVR paradigms was demonstrated when utilizing a GLM to regress out systemic/extracerebral signal contributions. Linear regression alone showed poor agreement across paradigms ( CCC ≤ 0.35 ; R ≤ 0.45 ).

Conclusions

More tolerable experimental paradigms coupled with regression of systemic/extracerebral signal contributions may offer a viable alternative to HC for assessing CVR.

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.

Quantitative cerebrovascular reactivity MRI in mice using acetazolamide challenge

PURPOSE

To develop a quantitative MRI method to estimate cerebrovascular reactivity (CVR) in mice.

METHODS

We described an MRI procedure to measure cerebral vasodilatory response to acetazolamide (ACZ), a vasoactive agent previously used in human clinical imaging. Vascular response was determined by cerebral blood flow (CBF) measured with phase-contrast or pseudo-continuous arterial spin labeling MRI. Vasodilatory input intensity was determined by plasma ACZ level using high-performance liquid chromatography. We verified the source of the CVR MRI signal by comparing ACZ injection to phosphate-buffered saline injection and noninjection experiments. Dose dependence and feasibility of regional CVR measurement were also investigated.

RESULTS

Cerebral blood flow revealed an exponential increase following intravenous ACZ injection, with a time constant of 1.62 min. In contrast, phosphate-buffered saline or noninjection exhibited a slow linear CBF increase, consistent with a gradual accumulation of anesthetic agent, isoflurane, used in this study. When comparing different ACZ doses, injections of 30, 60, 120, and 180 mg/kg yielded a linear increase in plasma ACZ concentration (p < 0.0001). On the other hand, CBF changes under these doses were not different from each other (p = 0.50). The pseudo-continuous arterial spin labeling MRI with multiple postlabeling delays revealed similar vascular responses at different postlabeling delay values. There was a regional difference in CVR (p = 0.005), with isocortex (0.81 ± 0.17%/[μg/ml]) showing higher CVR than deep-brain regions. Mice receiving multiple ACZ injections lived for a minimum of 6 months after the study without noticeable aberrant behavior or appearance.

CONCLUSIONS

We demonstrated the proof-of-principle of a new quantitative CVR mapping technique in mice.

fMRI studies evaluating central respiratory control in humans

A plethora of neural centers in the central nervous system control the fundamental respiratory pattern. This control is ensured by neurons that act as pacemakers, modulating activity through chemical control driven by changes in the O₂/CO₂ balance. Most of the respiratory neural centers are located in the brainstem, but difficult to localize on magnetic resonance imaging (MRI) due to their small size, lack of visually-detectable borders with neighboring areas, and significant physiological noise hampering detection of its activity with functional MRI (fMRI). Yet, several approaches make it possible to study the normal response to different abnormal stimuli or conditions such as CO₂ inhalation, induced hypercapnia, volitional apnea, induced hypoxia etc. This review provides a comprehensive overview of the majority of available studies on central respiratory control in humans.

Lawrence K, Tachtsidis I, Tong Y, Torricelli A, Urner T, Wabnitz H, Wolf M, Wolf U, Xu S, Yang C, Yodh AG, Yücel MA, Zhou W. Optical imaging and spectroscopy for the study of the human brain: status report

This report is the second part of a comprehensive two-part series aimed at reviewing an extensive and diverse toolkit of novel methods to explore brain health and function. While the first report focused on neurophotonic tools mostly applicable to animal studies, here, we highlight optical spectroscopy and imaging methods relevant to noninvasive human brain studies. We outline current state-of-the-art technologies and software advances, explore the most recent impact of these technologies on neuroscience and clinical applications, identify the areas where innovation is needed, and provide an outlook for the future directions.

NOise Reduction with DIstribution Corrected (NORDIC) PCA improves signal-to-noise in rodent resting-state and optogenetic functional MRI

NOise Reduction with DIstribution Corrected (NORDIC) principal component analysis (PCA) has been shown to selectively suppress thermal noise and improve temporal signal-to-noise ratio (tSNR) in human functional magnetic resonance imaging (fMRI). However, the feasibility to improve rodent fMRI using NORDIC PCA has not been explored. In this study, we developed a rodent fMRI preprocessing pipeline by incorporating NORDIC and evaluated its performance in a range of rodent fMRI applications from resting-state fMRI to task-evoked fMRI using optogenetics. In resting-state fMRI, we demonstrated a significant increase in tSNR by more than 3 times after NORDIC correction with reduced variance and improved task-free relative cerebrovascular reactivity (rCVR) across cortical depth. In optogenetic fMRI, apart from tSNR increase, more activated voxels and a significant decrease in the variance of activated brain signals were observed after NORDIC correction without apparent change in brain morphology. Taken together, our results signified the values of NORDIC correction for better detection of brain activities in rodent fMRI. Clinical Relevance: NORDIC PCA increases temporal signalto- noise ratio in rodent resting-state and task-evoked functional MRI, which can play an important role in improving the image quality for translational medicine and preclinical research, and for guiding future clinical neuroimaging.

Cerebrovascular reactivity deficits in cognitively unimpaired older adults: vasodilatory versus vasoconstrictive responses

Cerebrovascular reactivity (CVR) deficits may index vulnerability to vascular brain injury and cognitive impairment, but findings on age-related changes in CVR have been mixed, and no studies to date have directly compared age-related changes in CVR to hypercapnia versus hypocapnia. The present study compared CVR in 31 cognitively unimpaired older adults (ages 55-87) and 30 healthy younger adults (ages 18-28). Breath control tasks induced CVR to hypocapnia (0.1 Hz paced breathing) and hypercapnia (15s breath holds) during pseudo-continuous arterial spin labeling MRI. Relative to younger adults, cognitively unimpaired older adults displayed lower levels of global CVR under both hypocapnia and hypercapnia. In region-of-interest analyses, older adults exhibited attenuated CVR to hypocapnia in select frontal and temporal regions, and lower CVR to hypercapnia in all cortical, limbic, and subcortical regions examined, relative to younger adults. Results indicate age-related deficits in CVR are detectible even in cognitively unimpaired older adults and are disproportionately related to vasodilatory (hypercapnia) responses relative to vasoconstrictive (hypocapnia) responses. Findings may offer means for early detection of cerebrovascular dysfunction.

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.

The Utility of Cerebrovascular Reactivity MRI in Brain Rehabilitation: A Mechanistic Perspective

Cerebrovascular control and its integration with other physiological systems play a key role in the effective maintenance of homeostasis in brain functioning. Maintenance, restoration, and promotion of such a balance are one of the paramount goals of brain rehabilitation and intervention programs. Cerebrovascular reactivity (CVR), an index of cerebrovascular reserve, plays an important role in chemo-regulation of cerebral blood flow. Improved vascular reactivity and cerebral blood flow are important factors in brain rehabilitation to facilitate desired cognitive and functional outcomes. It is widely accepted that CVR is impaired in aging, hypertension, and cerebrovascular diseases and possibly in neurodegenerative syndromes. However, a multitude of physiological factors influence CVR, and thus a comprehensive understanding of underlying mechanisms are needed. We are currently underinformed on which rehabilitation method will improve CVR, and how this information can inform on a patient's prognosis and diagnosis. Implementation of targeted rehabilitation regimes would be the first step to elucidate whether such regimes can modulate CVR and in the process may assist in improving our understanding for the underlying vascular pathophysiology. As such, the high spatial resolution along with whole brain coverage offered by MRI has opened the door to exciting recent developments in CVR MRI. Yet, several challenges currently preclude its potential as an effective diagnostic and prognostic tool in treatment planning and guidance. Understanding these knowledge gaps will ultimately facilitate a deeper understanding for cerebrovascular physiology and its role in brain function and rehabilitation. Based on the lessons learned from our group's past and ongoing neurorehabilitation studies, we present a systematic review of physiological mechanisms that lead to impaired CVR in aging and disease, and how CVR imaging and its further development in the context of brain rehabilitation can add value to the clinical settings.

Cerebrovascular Reactivity Mapping Using Resting-State BOLD Functional MRI in Healthy Adults and Patients with Moyamoya Disease

Background Cerebrovascular reserve, the potential capacity of brain tissue to receive more blood flow when needed, is a desirable marker in evaluating ischemic risk. However, current measurement methods require acetazolamide injection or hypercapnia challenge, prompting a clinical need for resting-state (RS) blood oxygen level-dependent (BOLD) functional MRI data to measure cerebrovascular reactivity (CVR). Purpose To optimize and evaluate an RS CVR MRI technique and demonstrate its relationship to neurosurgical treatment. Materials and Methods In this HIPAA-compliant study, RS BOLD functional MRI data collected in 170 healthy controls between December 2008 and September 2010 were retrospectively evaluated to identify the optimal frequency range of temporal filtering on the basis of spatial correlation with the reference standard CVR map obtained with CO₂ inhalation. Next, the optimized RS method was applied in a new, prospective cohort of 50 participants with Moyamoya disease who underwent imaging between June 2014 and August 2019. Finally, CVR values were compared between brain hemispheres with and brain hemispheres without revascularization surgery by using Mann-Whitney U test. Results A total of 170 healthy controls (mean age ± standard deviation, 51 years ± 20; 105 women) and 100 brain hemispheres of 50 participants with Moyamoya disease (mean age, 41 years ± 12; 43 women) were evaluated. RS CVR maps based on a temporal filtering frequency of [0, 0.1164 Hz] yielded the highest spatial correlation (r = 0.74) with the CO₂ inhalation CVR results. In patients with Moyamoya disease, 77 middle cerebral arteries (MCAs) had stenosis. RS CVR in the MCA territory was lower in the group that did not undergo surgery (n = 30) than in the group that underwent surgery (n = 47) (mean, 0.407 relative units [ru] ± 0.208 vs 0.532 ru ± 0.182, respectively; P = .006), which is corroborated with the CO₂ inhalation CVR data (mean, 0.242 ru ± 0.273 vs 0.437 ru ± 0.200; P = .003). Conclusion Cerebrovascular reactivity mapping performed by using resting-state blood oxygen level-dependent functional MRI provided a task-free method to measure cerebrovascular reserve and depicted treatment effect of revascularization surgery in patients with Moyamoya disease comparable to that with the reference standard of CO₂ inhalation MRI. © RSNA, 2021 Online supplemental material is available for this article. See also the editorial by Wolf and Ware in this issue.

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