A novel method of quantifying hemodynamic delays to improve hemodynamic response, and CVR estimates in CO2 challenge fMRI

Effects of clamping end-tidal CO2 on neurofluidic low-frequency oscillations

In recent years, low-frequency oscillations (LFOs) (0.01-0.1 Hz) have been a subject of interest in resting-state functional magnetic resonance imaging research. They are believed to have many possible driving mechanisms, from both regional and global sources. Internal fluctuations in the partial pressure of CO2 (PCO2) has long been thought of as one of these major driving forces, but its exact contributions compared with other mechanisms have yet to be fully understood. This study examined the effects of end-tidal PCO2 (PetCO2) oscillations on LF cerebral hemodynamics and cerebrospinal fluid (CSF) dynamics under "clamped PetCO2" and "free-breathing" conditions. Under clamped PetCO2, a participant's PetCO2 levels were fixed to their baseline average, whereas PetCO2 was not controlled in free breathing. Under clamped PetCO2, the fractional amplitude of hemodynamic LFOs in the occipital and sensorimotor cortex and temporal lobes were found to be significantly reduced. Additionally, the fractional amplitude of CSF LFOs, measured at the fourth ventricle, was found to be reduced by almost one-half. However, the spatiotemporal distributions of blood and CSF delay times, as measured by cross-correlation in the LF domain, were not significantly altered between conditions. This study demonstrates that, while PCO2 oscillations significantly mediate LFOs, especially those observed in the CSF, other mechanisms are able to maintain LFOs, with high correlation, even in their absence.

Assessment of MR blood-oxygen-level-dependent (BOLD) cerebrovascular reactivity under general anesthesia in children with moyamoya

BACKGROUND

Moyamoya is a progressive, non-atherosclerotic cerebral arteriopathy that may present in childhood and currently has no cure. Early diagnosis is critical to prevent a lifelong risk of neurological morbidity. Blood-oxygen-level-dependent (BOLD) MRI cerebrovascular reactivity (CVR) imaging provides a non-invasive, in vivo measure of autoregulatory capacity and cerebrovascular reserve. However, non-compliant or younger children require general anesthesia to achieve BOLD-CVR imaging.

OBJECTIVE

To determine the same-day repeatability of BOLD-CVR imaging under general anesthesia in children with moyamoya.

MATERIALS AND METHODS

Twenty-eight examination pairs were included (mean patient age = 7.3 ± 4.0 years). Positive and negatively reacting voxels were averaged over signals and counted over brain tissue and vascular territory. The intraclass correlation coefficient (ICC), Wilcoxon signed-rank test, and Bland-Altman plots were used to assess the variability between the scans.

RESULTS

There was excellent-to-good (≥ 0.59) within-day repeatability in 18 out of 28 paired studies (64.3%). Wilcoxon signed-rank tests demonstrated no significant difference in the grey and white matter CVR estimates, between repeat scans (all p-values > 0.05). Bland-Altman plots of differences in mean magnitude of positive and negative and fractional positive and negative CVR estimates illustrated a reasonable degree of agreement between repeat scans and no systematic bias.

CONCLUSION

BOLD-CVR imaging provides repeatable assessment of cerebrovascular reserve in children with moyamoya imaged under general anesthesia.

Resting-State fMRI: Emerging Concepts for Future Clinical Application

Resting-state functional magnetic resonance imaging (rsfMRI) has been developed as a method of investigating spontaneous neural activity. Based on its low-frequency signal synchronization, rsfMRI has made it possible to identify multiple macroscopic structures termed resting-state networks (RSNs) on a single scan of less than 10 minutes. It is easy to implement even in clinical practice, in which assigning tasks to patients can be challenging. These advantages have accelerated the adoption and growth of rsfMRI. Recently, studies on the global rsfMRI signal have attracted increasing attention. Because it primarily arises from physiological events, less attention has hitherto been paid to the global signal than to the local network (i.e., RSN) component. However, the global signal is not a mere nuisance or a subsidiary component. On the contrary, it is quantitatively the dominant component that accounts for most of the variance in the rsfMRI signal throughout the brain and provides rich information on local hemodynamics that can serve as an individual-level diagnostic biomarker. Moreover, spatiotemporal analyses of the global signal have revealed that it is closely and fundamentally associated with the organization of RSNs, thus challenging the basic assumptions made in conventional rsfMRI analyses and views on RSNs. This review introduces new concepts emerging from rsfMRI spatiotemporal analyses focusing on the global signal and discusses how they may contribute to future clinical medicine. EVIDENCE LEVEL: 5 TECHNICAL EFFICACY: Stage 1.

An Investigation of Cerebral Vascular Functional Properties in Middle-to-Old Age Community People With High Vascular Risk Profiles

BACKGROUND

Vascular degeneration is an important cause of brain damage in aging. Assessing the functional properties of the cerebral vascular system may aid early diagnosis and prevention.

PURPOSE

To investigate the relationships between potential vascular functional markers and vascular risks, brain parenchymal damage, and cognition.

STUDY TYPE

Retrospective.

SUBJECTS

Two hundred two general community subjects (42-80 years, males/females: 127/75).

FIELD STRENGTH/SEQUENCE

3 T, spin echo T1W/T2W/FLAIR, resting-state functional MRI with an echo-planar sequence (rsfMRI), pseudo-continuous arterial spin labeling (pCASL) with a three-dimensional gradient-spin echo sequence.

ASSESSMENT

Cerebral blood flow (CBF) in gray matter calculated using pCASL, blood transit times calculated using rsfMRI, and the SD of internal carotid arteries signal (ICAstd ) calculated using rsfMRI; visual assessment for lacunes; quantification of white matter hyperintensity volume; permutation test for quality control; collection of demographic and clinical data, Montreal Cognitive Assessment, Mini-Mental State Examination.

STATISTICAL TESTS

Kolmogorov-Smirnov test; Spearman rank correlation analysis; Multivariable linear regression analysis controlling for covariates; The level of statistical significance was set at P < 0.05.

RESULTS

Age was negatively associated with ICAstd (β = -0.180). Diabetes was associated with longer blood transit time from large arteries to capillary bed (β = 0.185, adjusted for age, sex, and intracranial volume). Larger ICAstd was associated with less presence of lacunes (odds ratio: 0.418, adjusted for age and sex). Higher gray matter CBF (β = 0.154) and larger ICAstd (β = 0.136) were associated with better MoCA scores (adjusted for age, sex, and education).

DATA CONCLUSION

Prolonged blood transit time, decreased ICAstd , and diminished CBF were associated with vascular dysfunction and cognitive impairment. They may serve as vascular functional markers in future studies.

EVIDENCE LEVEL

3 TECHNICAL EFFICACY: Stage 3.

Resting-state networks representation of the global phenomena

Resting-state functional magnetic resonance imaging (rsfMRI) has been widely applied to investigate spontaneous neural activity, often based on its macroscopic organization that is termed resting-state networks (RSNs). Although the neurophysiological mechanisms underlying the RSN organization remain largely unknown, accumulating evidence points to a substantial contribution from the global signals to their structured synchronization. This study further explored the phenomenon by taking advantage of the inter- and intra-subject variations of the time delay and correlation coefficient of the signal timeseries in each region using the global mean signal as the reference signal. Consistent with the hypothesis based on the empirical and theoretical findings, the time lag and correlation, which have consistently been proven to represent local hemodynamic status, were shown to organize networks equivalent to RSNs. The results not only provide further evidence that the local hemodynamic status could be the direct source of the RSNs' spatial patterns but also explain how the regional variations in the hemodynamics, combined with the changes in the global events' power spectrum, lead to the observations. While the findings pose challenges to interpretations of rsfMRI studies, they further support the view that rsfMRI can offer detailed information related to global neurophysiological phenomena as well as local hemodynamics that would have great potential as biomarkers.

Examining temporal features of BOLD-based cerebrovascular reactivity in clinical populations

Background

Conventional cerebrovascular reactivity (CVR) estimation has demonstrated that many brain diseases and/or conditions are associated with altered CVR. Despite the clinical potential of CVR, characterization of temporal features of a CVR challenge remains uncommon. This work is motivated by the need to develop CVR parameters that characterize individual temporal features of a CVR challenge.

Methods

Data were collected from 54 adults and recruited based on these criteria: (1) Alzheimer's disease diagnosis or subcortical Vascular Cognitive Impairment, (2) sleep apnea, and (3) subjective cognitive impairment concerns. We investigated signal changes in blood oxygenation level dependent (BOLD) contrast images with respect to hypercapnic and normocapnic CVR transition periods during a gas manipulation paradigm. We developed a model-free, non-parametric CVR metric after considering a range of responses through simulations to characterize BOLD signal changes that occur when transitioning from normocapnia to hypercapnia. The non-parametric CVR measure was used to examine regional differences across the insula, hippocampus, thalamus, and centrum semiovale. We also examined the BOLD signal transition from hypercapnia back to normocapnia.

Results

We found a linear association between isolated temporal features of successive CO₂ challenges. Our study concluded that the transition rate from hypercapnia to normocapnia was significantly associated with the second CVR response across all regions of interest (p < 0.001), and this association was highest in the hippocampus (R² = 0.57, p < 0.0125).

Conclusion

This study demonstrates that it is feasible to examine individual responses associated with normocapnic and hypercapnic transition periods of a BOLD-based CVR experiment. Studying these features can provide insight on between-subject differences in CVR.

The many layers of BOLD. The effect of hypercapnic and hyperoxic stimuli on macro- and micro-vascular compartments quantified by CVR, M, and CBV across cortical depth

Ultra-high field functional magnetic resonance imaging (fMRI) offers the spatial resolution to measure neuronal activity at the scale of cortical layers. However, cortical depth dependent vascularization differences, such as a higher prevalence of macro-vascular compartments near the pial surface, have a confounding effect on depth-resolved blood-oxygen-level dependent (BOLD) fMRI signals. In the current study, we use hypercapnic and hyperoxic breathing conditions to quantify the influence of all venous vascular and micro-vascular compartments on laminar BOLD fMRI, as measured with gradient-echo (GE) and spin-echo (SE) scan sequences, respectively. We find that all venous vascular and micro-vascular compartments are capable of comparable theoretical maximum signal intensities, as represented by the M-value parameter. However, the capacity for vessel dilation, as reflected by the cerebrovascular reactivity (CVR), is approximately two and a half times larger for all venous vascular compartments combined compared to the micro-vasculature at superficial layers. Finally, there is roughly a 35% difference in estimates of CBV changes between all venous vascular and micro-vascular compartments, although this relative difference was approximately uniform across cortical depth. Thus, our results suggest that fMRI BOLD signal differences across cortical depth are likely caused by differences in dilation properties between macro- and micro-vascular compartments.

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.

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.

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.

Measuring Cerebrovascular Reactivity: Sixteen Avoidable Pitfalls

An increase in arterial PCO₂ is the most common stressor used to increase cerebral blood flow for assessing cerebral vascular reactivity (CVR). That CO₂ is readily obtained, inexpensive, easy to administer, and safe to inhale belies the difficulties in extracting scientifically and clinically relevant information from the resulting flow responses. Over the past two decades, we have studied more than 2,000 individuals, most with cervical and cerebral vascular pathology using CO₂ as the vasoactive agent and blood oxygen-level-dependent magnetic resonance imaging signal as the flow surrogate. The ability to deliver different forms of precise hypercapnic stimuli enabled systematic exploration of the blood flow-related signal changes. We learned the effect on CVR of particular aspects of the stimulus such as the arterial partial pressure of oxygen, the baseline PCO₂, and the magnitude, rate, and pattern of its change. Similarly, we learned to interpret aspects of the flow response such as its magnitude, and the speed and direction of change. Finally, we were able to test whether the response falls into a normal range. Here, we present a review of our accumulated insight as 16 “lessons learned.” We hope many of these insights are sufficiently general to apply to a range of types of CO₂-based vasoactive stimuli and perfusion metrics used for CVR.

The Role of Cerebrovascular-Reactivity Mapping in Functional MRI: Calibrated fMRI and Resting-State fMRI

Task and resting-state functional MRI (fMRI) is primarily based on the same blood-oxygenation level-dependent (BOLD) phenomenon that MRI-based cerebrovascular reactivity (CVR) mapping has most commonly relied upon. This technique is finding an ever-increasing role in neuroscience and clinical research as well as treatment planning. The estimation of CVR has unique applications in and associations with fMRI. In particular, CVR estimation is part of a family of techniques called calibrated BOLD fMRI, the purpose of which is to allow the mapping of cerebral oxidative metabolism (CMRO2) using a combination of BOLD and cerebral-blood flow (CBF) measurements. Moreover, CVR has recently been shown to be a major source of vascular bias in computing resting-state functional connectivity, in much the same way that it is used to neutralize the vascular contribution in calibrated fMRI. Furthermore, due to the obvious challenges in estimating CVR using gas challenges, a rapidly growing field of study is the estimation of CVR without any form of challenge, including the use of resting-state fMRI for that purpose. This review addresses all of these aspects in which CVR interacts with fMRI and the role of CVR in calibrated fMRI, provides an overview of the physiological biases and assumptions underlying hypercapnia-based CVR and calibrated fMRI, and provides a view into the future of non-invasive CVR measurement.