Examining the regional and cerebral depth-dependent BOLD cerebrovascular reactivity response at 7 T

CO2 as an engine for neurofluid flow: Exploring the coupling between vascular reactivity, brain clearance, and changes in tissue properties

The brain relies on an effective clearance mechanism to remove metabolic waste products for the maintenance of homeostasis. Recent studies have focused on elucidating the forces that drive the motion of cerebrospinal fluid (CSF), responsible for removal of these waste products. We demonstrate that vascular responses evoked using controlled manipulations of partial pressure of carbon dioxide (PaCO2) levels, serve as an endogenous driver of CSF clearance from the brain. To demonstrate this, we retrospectively surveyed our database, which consists of brain metastases patients from whom blood oxygen level-dependent (BOLD) images were acquired during targeted hypercapnic and hyperoxic respiratory challenges. We observed a correlation between CSF inflow signal around the fourth ventricle and CO2-induced changes in cerebral blood volume. By contrast, no inflow signal was observed in response to the nonvasoactive hyperoxic stimulus, validating our measurements. Moreover, our results establish a link between the rate of the hemodynamic response (to elevated PaCO2) and peritumoral edema load, which we suspect may affect CSF flow, consequently having implications for brain clearance. Our expanded perspective on the factors involved in neurofluid flow underscores the importance of considering both cerebrovascular responses, as well as the brain mechanical properties, when evaluating CSF dynamics in the context of disease processes.

MRI mapping of hemodynamics in the human spinal cord

Impaired spinal cord vascular function contributes to numerous neurological pathologies, making it important to be able to noninvasively characterize these changes. Here, we propose a functional magnetic resonance imaging (fMRI)-based method to map spinal cord vascular reactivity (SCVR). We used a hypercapnic breath-holding task, monitored with end-tidal CO2 (PETCO2), to evoke a systemic vasodilatory response during concurrent blood oxygenation level-dependent (BOLD) fMRI. SCVR amplitude and hemodynamic delay were mapped at the group level in 27 healthy participants as proof-of-concept of the approach, and then in two highly-sampled participants to probe feasibility/stability of individual SCVR mapping. Across the group and the highly-sampled individuals, a strong ventral SCVR amplitude was initially observed without accounting for local regional variation in the timing of the vasodilatory response. Shifted breathing traces (PETCO2) were used to account for temporal differences in the vasodilatory response across the spinal cord, producing maps of SCVR delay. These delay maps reveal an earlier ventral and later dorsal response and demonstrate distinct gray matter regions concordant with territories of arterial supply. The SCVR fMRI methods described here enable robust mapping of spatiotemporal hemodynamic properties of the human spinal cord. This noninvasive approach has exciting potential to provide early insight into pathology-driven vascular changes in the cord, which may precede and predict future irreversible tissue damage and guide the treatment of several neurological pathologies involving the spine.

Clinical targeting of the cerebral oxygen cascade to improve brain oxygenation in patients with hypoxic-ischaemic brain injury after cardiac arrest

The cerebral oxygen cascade includes three key stages: (a) convective oxygen delivery representing the bulk flow of oxygen to the cerebral vascular bed; (b) diffusion of oxygen from the blood into brain tissue; and (c) cellular utilisation of oxygen for aerobic metabolism. All three stages may become dysfunctional after resuscitation from cardiac arrest and contribute to hypoxic-ischaemic brain injury (HIBI). Improving convective cerebral oxygen delivery by optimising cerebral blood flow has been widely investigated as a strategy to mitigate HIBI. However, clinical trials aimed at optimising convective oxygen delivery have yielded neutral results. Advances in the understanding of HIBI pathophysiology suggest that impairments in the stages of the oxygen cascade pertaining to oxygen diffusion and cellular utilisation of oxygen should also be considered in identifying therapeutic strategies for the clinical management of HIBI patients. Culprit mechanisms for these impairments may include a widening of the diffusion barrier due to peri-vascular oedema and mitochondrial dysfunction. An integrated approach encompassing both intra-parenchymal and non-invasive neuromonitoring techniques may aid in detecting pathophysiologic changes in the oxygen cascade and enable patient-specific management aimed at reducing the severity of HIBI.

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.

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

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

Assessing Cerebrovascular Resistance in Patients With Sickle Cell Disease

In patients with sickle cell disease (SCD) the delivery of oxygen to the brain is compromised by anemia, abnormal rheology, and steno-occlusive vascular disease. Meeting demands for oxygen delivery requires compensatory features of brain perfusion. The cerebral vasculature’s regulatory function and reserves can be assessed by observing the flow response to a vasoactive stimulus. In a traditional approach we measured voxel-wise change in Blood Oxygen-Level Dependent (BOLD) MRI signal as a surrogate of cerebral blood flow (CBF) in response to a linear progressive ramping of end-tidal partial pressure of carbon dioxide (PETCO₂). Cerebrovascular reactivity (CVR) was defined as ΔBOLD/ΔPETCO₂. We used a computer model to fit a virtual sigmoid resistance curve to the progressive CBF response to the stimulus, enabling the calculation of resistance parameters: amplitude, midpoint, range response, resistance sensitivity and vasodilatory reserve. The quality of the resistance sigmoid fit was expressed as the r² of the fit. We tested 35 patients with SCD, as well as 24 healthy subjects to provide an indication of the normal ranges of the resistance parameters. We found that gray matter CVR and resistance amplitude, range, reserve, and sensitivity are reduced in patients with SCD compared to healthy controls, while resistance midpoint was increased. This study is the first to document resistance measures in adult patients with SCD. It is also the first to score these vascular resistance measures in comparison to the normal range. We anticipate these data will complement the current understanding of the cerebral vascular pathophysiology of SCD, identify paths for therapeutic interventions, and provide biomarkers for monitoring the progress of the disease.

Multi-vendor and multisite evaluation of cerebrovascular reactivity mapping using hypercapnia challenge

Cerebrovascular reactivity (CVR), which measures the ability of cerebral blood vessels to dilate or constrict in response to vasoactive stimuli such as CO2 inhalation, is an important index of the brain's vascular health. Quantification of CVR using BOLD MRI with hypercapnia challenge has shown great promises in research and clinical studies. However, in order for it to be used as a potential imaging biomarker in large-scale and multi-site studies, the reliability of CO2-CVR quantification across different MRI acquisition platforms and researchers/raters must be examined. The goal of this report from the MarkVCID small vessel disease biomarkers consortium is to evaluate the reliability of CO2-CVR quantification in three studies. First, the inter-rater reliability of CO2-CVR data processing was evaluated by having raters from 5 MarkVCID sites process the same 30 CVR datasets using a cloud-based CVR data processing pipeline. Second, the inter-scanner reproducibility of CO2-CVR quantification was assessed in 10 young subjects across two scanners of different vendors. Third, test-retest repeatability was evaluated in 20 elderly subjects from 4 sites with a scan interval of less than 2 weeks. In all studies, the CO2 CVR measurements were performed using the fixed inspiration method, where the subjects wore a nose clip and a mouthpiece and breathed room air and 5% CO2 air contained in a Douglas bag alternatively through their mouth. The results showed that the inter-rater CoV of CVR processing was 0.08 ± 0.08% for whole-brain CVR values and ranged from 0.16% to 0.88% in major brain regions, with ICC of absolute agreement above 0.9959 for all brain regions. Inter-scanner CoV was found to be 6.90 ± 5.08% for whole-brain CVR values, and ranged from 4.69% to 12.71% in major brain regions, which are comparable to intra-session CoVs obtained from the same scanners on the same day. ICC of consistency between the two scanners was 0.8498 for whole-brain CVR and ranged from 0.8052 to 0.9185 across major brain regions. In the test-retest evaluation, test-retest CoV across different days was found to be 18.29 ± 17.12% for whole-brain CVR values, and ranged from 16.58% to 19.52% in major brain regions, with ICC of absolute agreement ranged from 0.6480 to 0.7785. These results demonstrated good inter-rater, inter-scanner, and test-retest reliability in healthy volunteers, and suggested that CO2-CVR has suitable instrumental properties for use as an imaging biomarker of cerebrovascular function in multi-site and longitudinal observational studies and clinical trials.

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

Brain stress testing using blood oxygenation level-dependent (BOLD) MRI to evaluate changes in cerebrovascular reactivity (CVR) is of growing interest for evaluating white matter integrity. However, even under healthy conditions, the white matter BOLD-CVR response differs notably from that observed in the gray matter. In addition to actual arterial vascular control, the venous draining topology may influence the WM-CVR response leading to signal delays and dispersions. These types of alterations in hemodynamic parameters are sometimes linked with pathology, but may also arise from differences in normal venous architecture. In this work, high-resolution T2*weighted anatomical images combined with BOLD imaging during a hypercapnic breathing protocol were acquired using a 7 tesla MRI system. Hemodynamic parameters including base CVR, hemodynamic lag, lag-corrected CVR, response onset and signal dispersion, and finally ΔCVR (corrected CVR minus base CVR) were calculated in 8 subjects. Parameter maps were spatially normalized and correlated against an MNI-registered white matter medullary vein atlas. Moderate correlations (Pearson's rho) were observed between medullary vessel frequency (MVF) and ΔCVR (0.52; 0.58 for total WM), MVF and hemodynamic lag (0.42; 0.54 for total WM), MVF and signal dispersion (0.44; 0.53 for total WM), and finally MVF and signal onset (0.43; 0.52 for total WM). Results indicate that, when assessed in the context of the WM venous architecture, changes in the response shape may only be partially reflective of the actual vascular reactivity response occurring further upstream by control vessels. This finding may have implications when attributing diseases mechanisms and/or progression to presumed impaired WM BOLD-CVR.

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

Elevated carbon dioxide (CO2) in breathing air is widely used as a vasoactive stimulus to assess cerebrovascular functions under hypercapnia (i.e., "stress test" for the brain). Blood-oxygen-level-dependent (BOLD) is a contrast mechanism used in functional magnetic resonance imaging (fMRI). BOLD is used to study CO2-induced cerebrovascular reactivity (CVR), which is defined as the voxel-wise percentage BOLD signal change per mmHg change in the arterial partial pressure of CO2 (PaCO2). Besides the CVR, two additional important parameters reflecting the cerebrovascular functions are the arrival time of arterial CO2 at each voxel, and the waveform of the local BOLD signal. In this study, we developed a novel analytical method to accurately calculate the arrival time of elevated CO2 at each voxel using the systemic low frequency oscillations (sLFO: 0.01-0.1 Hz) extracted from the CO2 challenge data. In addition, 26 candidate hemodynamic response functions (HRF) were used to quantitatively describe the temporal brain reactions to a CO2 stimulus. We demonstrated that our approach improved the traditional method by allowing us to accurately map three perfusion-related parameters: the relative arrival time of blood, the hemodynamic response function, and CVR during a CO2 challenge.

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.

Compromised resting cerebral metabolism after sport-related concussion: A calibrated MRI study

Altered resting cerebral blood flow (CBF₀) in the acute phase post-concussion may contribute to neurobehavioral deficiencies, often reported weeks after the injury. However, in addition to changes in CBF₀, little is known about other physiological mechanisms that may be disturbed within the cerebrovasculature. The aim of this study was to assess whether changes in baseline perfusion following sport-related concussion (SRC) were co-localized with changes in cerebral metabolic demand. Forty-two subjects (15 SRC patients 8.0 ± 4.6 days post-injury and 27 age-matched healthy control athletes) were studied cross-sectionally. CBF₀, cerebrovascular reactivity (CVR), resting oxygen extraction (OEF₀) and cerebral metabolic rate of oxygen consumption (CMRO_2|0) were measured using a combination of hypercapnic and hyperoxic breathing protocols, and the biophysical model developed in calibrated MRI. Blood oxygenation level dependent and perfusion data were acquired simultaneously using a dual-echo arterial spin labelling sequence. SRC patients showed significant decreases in CBF₀ spread across the grey-matter (P < 0.05, corrected), and these differences were also confounded by the effects of baseline end-tidal CO₂ (P < 0.0001). Lower perfusion was co-localized with reductions in regional CMRO_2|0 (P = 0.006) post-SRC, despite finding no group-differences in OEF₀ (P = 0.800). Higher CVR within voxels showing differences in CBF was also observed in the SRC group (P = 0.001), compared to controls. Reductions in metabolic demand despite no significant changes in OEF₀ suggests that hypoperfusion post-SRC may reflect compromised metabolic function after the injury. These results provide novel insight about the possible pathophysiological mechanisms underlying concussion that may affect the clinical recovery of athletes after sport-related head injuries.

Cerebrovascular Reactivity: Purpose, Optimizing Methods, and Limitations to Interpretation - A Personal 20-Year Odyssey of (Re)searching

The brain is a neurovascular organ. A stimulus-response approach is effective in interrogating the physiology of its vasculature. Ideally, the stimulus is standardized across patients, and in a single patient over time. We developed a standard stimulus and attempted to measure, classify, and interpret the many forms of responses. Over the past 20 years, our work has delivered nuanced insights into normal cerebral vascular physiology, as well as adaptive physiological responses in the presence of disease. The trajectory of our understanding did not follow a logical linear progression; rather, it emerged as a coalescence of new, old, and previously dismissed, ideas that had accumulated over time. In this essay, we review what we believe were our most valuable - and sometimes controversial insights during our two decades-long journey.

Effects of head-down tilt bed rest plus elevated CO2 on cognitive performance

Microgravity and elevated CO2 levels are two important environmental spaceflight stressors that can adversely affect astronaut cognitive performance and jeopardize mission success. This study investigated the effects of 6° head-down tilt bed rest (HDBR) with (n = 11 participants, 30-day HDBR) and without (n = 8 participants, 60-day HDBR) elevated ambient (3.73 mmHg) CO2 concentrations on cognitive performance. Participants of both groups performed all 10 tests of NASA's Cognition battery and a brief alertness and mood survey repeatedly before, during, and after the HDBR period. Test scores were adjusted for practice and stimulus set effects. Concentrating on the first 30 days of HDBR, a modest but statistically significant slowing across a range of cognitive domains was found in both groups (controls: -0.37 SD; 95% CI -0.48, -0.27; adjusted P < 0.0001; CO2: -0.25 SD; 95% CI -0.34, -0.16; adjusted P < 0.001), most prominently for sensorimotor speed. These changes were observed early during HDBR and did not further deteriorate or improve with increasing time in HDBR. The study found similar cognitive effects of HDBR irrespective of CO2 levels, suggesting that elevated CO2 neither ameliorated nor worsened the HDBR effects. In both groups, cognitive performance after 15 days of recovery was statistically indistinguishable from pre-HDBR performance. However, subjects undergoing 60 days of HDBR rated themselves as feeling more sleepy, tired, physically exhausted, stressed, and unhealthy during recovery compared to their 30-day counterparts.NEW AND NOTEWORTHY This study investigated the effects of prolonged head-down tilt bed rest with and without elevated (3.73 mmHg) levels of ambient CO2 on cognitive performance across a range of cognitive domains and is one of the few studies investigating combined effects of environmental stressors prevalent in spaceflight. The study showed moderate declines in cognitive speed induced by head-down tilt bed rest and suggests that exposure to elevated levels of ambient CO2 did not modify this effect.

Effects of Spaceflight Stressors on Brain Volume, Microstructure, and Intracranial Fluid Distribution

Astronauts are exposed to elevated CO2 levels onboard the International Space Station. Here, we investigated structural brain changes in 11 participants following 30-days of head-down tilt bed rest (HDBR) combined with 0.5% ambient CO2 (HDBR + CO2) as a spaceflight analog. We contrasted brain changes observed in the HDBR + CO2 group with those of a previous HDBR sample not exposed to elevated CO2. Both groups exhibited a global upward shift of the brain and concomitant intracranial free water (FW) redistribution. Greater gray matter changes were seen in the HDBR + CO2 group in some regions. The HDBR + CO2 group showed significantly greater FW decrements in the posterior cerebellum and the cerebrum than the HDBR group. In comparison to the HDBR group, the HDBR + CO2 group exhibited greater diffusivity increases. In half of the participants, the HDBR + CO2 intervention resulted in signs of Spaceflight Associated Neuro-ocular Syndrome (SANS), a constellation of ocular structural and functional changes seen in astronauts. We therefore conducted an exploratory comparison compared between subjects that did and did not develop SANS and found asymmetric lateral ventricle enlargement in the SANS group. These results enhance our understanding of the underlying mechanisms of spaceflight-induced brain changes, which is critical for promoting astronaut health and performance.

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.

Insights Into Cerebral Tissue-Specific Response to Respiratory Challenges at 7T: Evidence for Combined Blood Flow and CO2-Mediated Effects

Cerebrovascular reactivity (CVR) mapping is finding increasing clinical applications as a non-invasive probe for vascular health. Further analysis extracting temporal delay information from the CVR response provide additional insight that reflect arterial transit time, blood redistribution, and vascular response speed. Untangling these factors can help better understand the (patho)physiology and improve diagnosis/prognosis associated with vascular impairments. Here, we use hypercapnic (HC) and hyperoxic (HO) challenges to gather insight about factors driving temporal delays between gray-matter (GM) and white-matter (WM). Blood Oxygen Level Dependent (BOLD) datasets were acquired at 7T in nine healthy subjects throughout BLOCK- and RAMP-HC paradigms. In a subset of seven participants, a combined HC+HO block, referred as the “BOOST” protocol, was also acquired. Tissue-based differences in Rapid Interpolation at Progressive Time Delays (RIPTiDe) were compared across stimulus to explore dynamic (BLOCK-HC) versus progressive (RAMP-HC) changes in CO₂, as well as the effect of bolus arrival time on CVR delays (BLOCK-HC versus BOOST). While GM delays were similar between the BLOCK- (21.80 ± 10.17 s) and RAMP-HC (24.29 ± 14.64 s), longer WM lag times were observed during the RAMP-HC (42.66 ± 17.79 s), compared to the BLOCK-HC (34.15 ± 10.72 s), suggesting that the progressive stimulus may predispose WM vasculature to longer delays due to the smaller arterial content of CO₂ delivered to WM tissues, which in turn, decreases intravascular CO₂ gradients modulating CO₂ diffusion into WM tissues. This was supported by a maintained ∼10 s offset in GM (11.66 ± 9.54 s) versus WM (21.40 ± 11.17 s) BOOST-delays with respect to the BLOCK-HC, suggesting that the vasoactive effect of CO₂ remains constant and that shortening of BOOST delays was be driven by blood arrival reflected through the non-vasodilatory HO contrast. These findings support that differences in temporal and magnitude aspects of CVR between vascular networks reflect a component of CO₂ sensitivity, in addition to redistribution and steal blood flow effects. Moreover, these results emphasize that the addition of a BOOST paradigm may provide clinical insights into whether vascular diseases causing changes in CVR do so by way of severe blood flow redistribution effects, alterations in vascular properties associated with CO₂ diffusion, or changes in blood arrival time.

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.

Persistent alterations in cerebrovascular reactivity in response to hypercapnia following pediatric mild traumatic brain injury

Much attention has been paid to the effects of mild traumatic brain injury (mTBI) on cerebrovascular reactivity in adult populations, yet it remains understudied in pediatric injury. In this study, 30 adolescents (12-18 years old) with pediatric mTBI (pmTBI) and 35 age- and sex-matched healthy controls (HC) underwent clinical and neuroimaging assessments during sub-acute (6.9 ± 2.2 days) and early chronic (120.4 ± 11.7 days) phases of injury. Relative to controls, pmTBI reported greater initial post-concussion symptoms, headache, pain, and anxiety, resolving by four months post-injury. Patients reported increased sleep issues and exhibited deficits in processing speed and attention across both visits. In grey-white matter interface areas throughout the brain, pmTBI displayed increased maximal fit/amplitude of a time-shifted end-tidal CO₂ regressor to blood oxygen-level dependent response relative to HC, as well as increased latency to maximal fit. The alterations persisted through the early chronic phase of injury, with maximal fit being associated with complaints of ongoing sleep disturbances during post hoc analyses but not cognitive measures of processing speed or attention. Collectively, these findings suggest that deficits in the speed and degree of cerebrovascular reactivity may persist longer than current conceptualizations about clinical recovery within 30 days.

Functional engagement of white matter in resting-state brain networks

The topological characteristics of functional networks, derived from measurements of resting-state connectivity in gray matter (GM), are associated with individual cognitive abilities or specific dysfunctions. However, blood oxygen level-dependent (BOLD) signals in white matter (WM) are usually ignored or even regressed out as nuisance factors in the data analyses that underlie network models. Recent studies have demonstrated reliable detection of WM BOLD signals and imply these reflect associated neural activities. Here we evaluate quantitatively the contributions of individual WM voxels to the identification of functional networks, which we term their engagement (or conceptually, their importance). We quantify the engagement by measuring the reductions of connectivity, produced by ignoring the signal fluctuations within each WM voxel, with respect to both the entire network (global) or a single GM node (local). We observed highly reproducible spatial distributions of global engagement maps, as well as a trend toward increased relevance of deep WM voxels at delayed times. Local engagement maps exhibit homogeneous spatial distributions with respect to internal nodes that constitute a well-recognized sub-functional network, but inhomogeneous distributions with respect to other nodes. WM voxels show distinct distributions of engagement depending on their anatomical locations. These findings demonstrate the important role of WM in network modeling, thus supporting the need for changes of conventional views that WM signal variations represent only physiological noise.

Predicting PET Cerebrovascular Reserve with Deep Learning by Using Baseline MRI: A Pilot Investigation of a Drug-Free Brain Stress Test

Background Cerebrovascular reserve (CVR) may be measured by using an acetazolamide test to clinically evaluate patients with cerebrovascular disease. However, acetazolamide use may be contraindicated and/or undesirable in certain clinical settings. Purpose To predict CVR images generated from acetazolamide vasodilation with a deep learning network by using only images before acetazolamide administration. Materials and Methods Simultaneous oxygen 15 (¹⁵O)-labeled water PET/MRI before and after acetazolamide injection were retrospectively analyzed for patients with Moyamoya disease and healthy control participants from April 2017 to May 2019. Inputs to deep learning models were perfusion-based images (arterial spin labeling [ASL]), structural scans (T2 fluid-attenuated inversion-recovery, T1), and brain location. Two models, that is, ¹⁵O-labeled water PET cerebral blood flow (CBF) and MRI (PET-plus-MRI model) before acetazolamide administration and only MRI (MRI-only model) before acetazolamide administration, were trained and tested with sixfold cross-validation. The models learned to predict a voxelwise relative CBF change (rΔCBF) map by using rΔCBF measured with PET due to acetazolamide as ground truth. Quantitative analysis included image quality metrics (peak signal-to-noise ratio, root mean square error, and structural similarity index), as well as comparison between the various methods by using correlation and Bland-Altman analyses. Identification of vascular territories with impaired rΔCBF was evaluated by using receiver operating characteristic metrics. Results Thirty-six participants were included: 24 patients with Moyamoya disease (mean age ± standard deviation, 41 years ± 12; 17 women) and 12 age-matched healthy control participants (mean age, 39 years ± 16; nine women). The rΔCBF maps predicted by both deep learning models demonstrated better image quality metrics than did ASL (all P < .001 in patients) and higher correlation coefficient with PET than with ASL (PET-plus-MRI model, 0.704; MRI-only model, 0.690 vs ASL, 0.432; both P < .001 in patients). Both models also achieved high diagnostic performance in identifying territories with impaired rΔCBF (area under receiver operating characteristic curve, 0.95 for PET-plus-MRI model [95% confidence interval: 0.90, 0.99] and 0.95 for MRI-only model [95% confidence interval: 0.91, 0.98]). Conclusion By using only images before acetazolamide administration, PET-plus-MRI and MRI-only deep learning models predicted cerebrovascular reserve images without the need for vasodilator injection. © RSNA, 2020 Online supplemental material is available for this article.

Reproducibility of a ramping protocol to measure cerebral vascular reactivity using functional magnetic resonance imaging

Though individual differences in arterial carbon dioxide and oxygen levels inherently exist, the degree of their influence on cerebral vascular reactivity (CVR) is less clear. We examined the reproducibility of BOLD signal changes to an iso-oxic ramping P_et CO₂ protocol. CVR changes were induced by altering P_et CO₂ while holding P_et O₂ constant using a computer-controlled sequential gas delivery (SGD) device. Two MRI scans, each including a linear change in P_et CO₂ , were performed using a 3-Tesla (3T) scanner. This ramp sequence consisted of 1 min at 30 mmHg followed by 4 min period during where P_et CO₂ was linearly increased from 30 to 50 mmHg, 1 min at 51 mmHg, and concluded with 4 min at baseline. The protocol was repeated at a separate visit with 3 days between visits (minimum). Intraclass correlation coefficients (ICC) and coefficients of variation (CV) were used to verify reproducibility. Eleven subjects (6 females; mean age 26.5 ± 5.7 years) completed the full testing protocol. Good reproducibility was observed for the within-visit ramp sequence (Visit 1: ICC = 0.82, CV = 6.5%; Visit 2: ICC = 0.74, CV = 6.4%). Similarly, ramp sequence were reproducible between visits (Scan 1: ICC = 0.74, CV = 6.5%; Scan 2: ICC = 0.66, CV = 6.1%). Establishing reproducible methodologies for measuring BOLD signal changes in response to P_et CO₂ alterations using a ramp protocol will allow researchers to study CVR functionality. Finally, adding a ramping protocol to CVR studies could provide information about changes in CVR over a broad range of P_et CO₂ .

Neural Correlates of Vestibular Processing During a Spaceflight Analog With Elevated Carbon Dioxide (CO2): A Pilot Study

Astronauts return to Earth from spaceflight missions with impaired mobility and balance; recovery can last weeks postflight. This is due in large part to the altered vestibular signaling and sensory reweighting that occurs in microgravity. The neural mechanisms of spaceflight-induced vestibular changes are not well understood. Head-down-tilt bed rest (HDBR) is a common spaceflight analog environment that allows for study of body unloading, fluid shifts, and other consequences of spaceflight. Subjects in this context still show vestibular changes despite being in Earth's gravitational environment, potentially due to sensory reweighting. Previously, we found evidence of sensory reweighting and reduced neural efficiency for vestibular processing in subjects who underwent a 70-day HDBR intervention. Here we extend this work by evaluating the impact of HDBR paired with elevated carbon dioxide (CO₂) to mimic International Space Station conditions on vestibular neural processing. Eleven participants (6 males, 34 ± 8 years) completed 30 days of HDBR combined with 0.5% atmospheric CO₂ (HDBR + CO₂). Participants underwent six functional magnetic resonance imaging (fMRI) sessions pre-, during, and post- HDBR + CO₂ while we measured brain activity in response to pneumatic skull taps (a validated method of vestibular stimulation). We also measured mobility and balance performance several times before and after the intervention. We found support for adaptive neural changes within the vestibular system during bed rest that subsequently recovered in several cortical and cerebellar regions. Further, there were multiple brain regions where greater pre- to post- deactivation was associated with reduced pre- to post- balance declines. That is, increased deactivation of certain brain regions associated with better balance post-HDBR + CO₂. We also found that, compared to HDBR alone (n = 13 males; 29 ± 3 years) HDBR + CO₂ is associated with greater increases in activation of multiple frontal, parietal, and temporal regions during vestibular stimulation. This suggests interactive or additive effects of bed rest and elevated CO₂. Finally, we found stronger correlations between pre- to post- HDBR + CO₂ brain changes and dependence on the visual system during balance for subjects who developed signs of Spaceflight-Associated Neuro-ocular Syndrome (SANS). Together, these findings have clear implications for understanding the neural mechanisms of bed rest and spaceflight-related changes in vestibular processing, as well as adaptation to altered sensory inputs.

Longitudinal Neuroimaging in Pediatric Traumatic Brain Injury: Current State and Consideration of Factors That Influence Recovery

Traumatic brain injury (TBI) is a leading cause of death and disability for children and adolescents in the U.S. and other developed and developing countries. Injury to the immature brain varies greatly from that of the mature, adult brain due to numerous developmental, pre-injury, and injury-related factors that work together to influence the trajectory of recovery during the course of typical brain development. Substantial damage to brain structure often underlies subsequent functional limitations that persist for years following pediatric TBI. Advances in neuroimaging have established an important role in the acute management of pediatric TBI, and magnetic resonance imaging (MRI) techniques have a particular relevance for the sequential assessment of long-term consequences from injuries sustained to the developing brain. The present paper will discuss the various factors that influence recovery and review the findings from the present neuroimaging literature to assess altered development and long-term outcome following pediatric TBI. Four MR-based neuroimaging modalities have been used to examine recovery from pediatric TBI longitudinally: (1) T1-weighted structural MRI is sensitive to morphological changes in gray matter volume and cortical thickness, (2) diffusion-weighted MRI is sensitive to changes in the microstructural integrity of white matter, (3) MR spectroscopy provides a sensitive assessment of metabolic and neurochemical alterations in the brain, and (4) functional MRI provides insight into the functional changes that occur as a result of structural damage and typical developmental processes. As reviewed in this paper, 13 cohorts have contributed to only 20 studies published to date using neuroimaging to examine longitudinal changes after TBI in pediatric patients. The results of these studies demonstrate considerable heterogeneity in post-injury outcome; however, the existing literature consistently shows that alterations in brain structure, function, and metabolism can persist for an extended period of time post-injury. With larger sample sizes and multi-site cooperation, future studies will be able to further examine potential moderators of outcome, such as the developmental, pre-injury, and injury-related factors discussed in the present review.

Validating faster DENSE measurements of cardiac-induced brain tissue expansion as a potential tool for investigating cerebral microvascular pulsations

Displacement Encoding with Stimulated Echoes (DENSE) has recently shown potential for measuring cardiac-induced cerebral volumetric strain in the human brain. As such, it may provide a powerful tool for investigating the cerebral small vessels. However, further development and validation are necessary. This study aims, first, to validate a retrospectively-gated implementation of the DENSE method for assessing brain tissue pulsations as a physiological marker, and second, to use the acquired measurements to explore intracranial volume dynamics. We acquired repeated measurements of cerebral volumetric strain in 8 healthy subjects, and internally validated these measurements by comparing them to spinal CSF stroke volumes obtained in the same scan session. Peak volumetric strain was found to be highly repeatable between scan sessions. First/second measured peak volumetric strains were: (6.4 ​± ​1.7)x10^-4/(6.7 ​± ​1.6)x10^-4 for whole brain, (9.5 ​± ​2.5)x10^-4/(9.6 ​± ​2.4)x10^-4 for grey matter, and (4.4 ​± ​1.7)x10^-4/(4.1 ​± ​0.8)x10^-4 for white matter. Grey matter showed significantly higher peak strain (p ​< ​0.001) and earlier time-to-peak strain (p ​< ​0.02) than white matter. An approximately linear relationship was found between CSF and brain tissue volume pulsations over the cardiac cycle (mean slope and R² of 0.88 ​± ​0.23 and 0.89 ​± ​0.07, respectively). The close similarity between CSF and brain tissue volume pulsations implies limited contributions from large intracranial vessel pulsations, providing further evidence for venous compression as an additional mechanism for maintaining stable intracranial pressures over the cardiac cycle. Cerebral pulsatility showed consistent inter-subject peak values in healthy subjects, and was strongly correlated to CSF stroke volumes. These results strengthen the potential of brain tissue volumetric strain as a means for investigating the intracranial dynamics of the ageing brain in normal or diseased states.

Head Down Tilt Bed Rest Plus Elevated CO2 as a Spaceflight Analog: Effects on Cognitive and Sensorimotor Performance

Long duration head down tilt bed rest (HDBR) has been widely used as a spaceflight analog environment to understand the effects of microgravity on human physiology and performance. Reports have indicated that crewmembers onboard the International Space Station (ISS) experience symptoms of elevated CO₂ such as headaches at lower levels of CO₂ than levels at which symptoms begin to appear on Earth. This suggests there may be combinatorial effects of elevated CO₂ and the other physiological effects of microgravity including headward fluid shifts and body unloading. The purpose of the current study was to investigate these effects by evaluating the impact of 30 days of 6° HDBR and 0.5% CO₂ (HDBR + CO₂) on mission relevant cognitive and sensorimotor performance. We found a facilitation of processing speed and a decrement in functional mobility for subjects undergoing HDBR + CO₂ relative to our previous study of HDBR in ambient air. In addition, nearly half of the participants in this study developed signs of Spaceflight Associated Neuro-ocular Syndrome (SANS), a constellation of ocular structural and functional changes seen in approximately one third of long duration astronauts. This allowed us the unique opportunity to compare the two subgroups. We found that participants who exhibited signs of SANS became more visually dependent and shifted their speed-accuracy tradeoff, such that they were slower but more accurate than those that did not incur ocular changes. These small subgroup findings suggest that SANS may have an impact on mission relevant performance inflight via sensory reweighting.

Characterization of the hemodynamic response function in white matter tracts for event-related fMRI

Accurate estimates of the BOLD hemodynamic response function (HRF) are crucial for the interpretation and analysis of event-related functional MRI data. To date, however, there have been no comprehensive measurements of the HRF in white matter (WM) despite increasing evidence that BOLD signals in WM change after a stimulus. We performed an event-related cognitive task (Stroop color-word interference) to measure the HRF in selected human WM pathways. The task was chosen in order to produce robust, distributed centers of activity throughout the cortex. To measure the HRF in WM, fiber tracts were reconstructed between each pair of activated cortical areas. We observed clear task-specific HRFs with reduced magnitudes, delayed onsets and prolonged initial dips in WM tracts compared with activated grey matter, thus calling for significant changes to current standard models for accurately characterizing the HRFs in WM and for modifications of standard methods of analysis of functional imaging data.

Neuroimaging of vascular reserve in patients with cerebrovascular diseases

Cerebrovascular reactivity, defined broadly as the ability of brain parenchyma to adjust cerebral blood flow in response to altered metabolic demand or a vasoactive stimulus, is being measured with increasing frequency and may have a use for portending new or recurrent stroke risk in patients with cerebrovascular disease. The purpose of this review is to outline (i) the physiological basis of variations in cerebrovascular reactivity, (ii) available approaches for measuring cerebrovascular reactivity in research and clinical settings, and (iii) clinically-relevant cerebrovascular reactivity findings in the context of patients with cerebrovascular disease, including atherosclerotic arterial steno-occlusion, non-atherosclerotic arterial steno-occlusion, anemia, and aging. Literature references summarizing safety considerations for these procedures and future directions for standardizing protocols and post-processing procedures across centers are presented in the specific context of major unmet needs in the setting of cerebrovascular disease.

Modeling of dynamic cerebrovascular reactivity to spontaneous and externally induced CO2 fluctuations in the human brain using BOLD-fMRI

In this work, we investigate the regional characteristics of the dynamic interactions between arterial CO₂ and BOLD (dynamic cerebrovascular reactivity - dCVR) during normal breathing and hypercapnic, externally induced step CO₂ challenges. To obtain dCVR curves at each voxel, we use a custom set of basis functions based on the Laguerre and gamma basis sets. This allows us to obtain robust dCVR estimates both in larger regions of interest (ROIs), as well as in individual voxels. We also implement classification schemes to identify brain regions with similar dCVR characteristics. Our results reveal considerable variability of dCVR across different brain regions, as well as during different experimental conditions (normal breathing and hypercapnic challenges), suggesting a differential response of cerebral vasculature to spontaneous CO₂ fluctuations and larger, externally induced CO₂ changes that are possibly associated with the underlying differences in mean arterial CO₂ levels. The clustering results suggest that anatomically distinct brain regions are characterized by different dCVR curves that in some cases do not exhibit the standard, positive valued curves that have been previously reported. They also reveal a consistent set of dCVR cluster shapes for resting and forcing conditions, which exhibit different distribution patterns across brain voxels.

Cerebral arterial time constant calculated from the middle and posterior cerebral arteries in healthy subjects

The cerebral arterial blood volume changes (∆C_aBV) during a single cardiac cycle can be estimated using transcranial Doppler ultrasonography (TCD) by assuming pulsatile blood inflow, constant, and pulsatile flow forward from large cerebral arteries to resistive arterioles [continuous flow forward (CFF) and pulsatile flow forward (PFF)]. In this way, two alternative methods of cerebral arterial compliance (C_a) estimation are possible. Recently, we proposed a TCD-derived index, named the time constant of the cerebral arterial bed (τ), which is a product of C_a and cerebrovascular resistance and is independent of the diameter of the insonated vessel. In this study, we aim to examine whether the τ estimated by either the CFF or the PFF model differs when calculated from the middle cerebral artery (MCA) and the posterior cerebral artery (PCA). The arterial blood pressure and TCD cerebral blood flow velocity (CBFV_a) in the MCA and in the PCA were non-invasively measured in 32 young, healthy volunteers (median age: 24, minimum age: 18, maximum age: 31). The τ was calculated using both the PFF and CFF models from the MCA and the PCA and compared using a non-parametric Wilcoxon signed-rank test. Results are presented as medians (25th-75th percentiles). The cerebrovascular time constant estimated in both arteries using the PFF model was shorter than when using the CFF model (ms): [64.83 (41.22-104.93) vs. 178.60 (160.40-216.70), p < 0.001 in the MCA, and 44.04 (17.15-81.17) vs. 183.50 (153.65-204.10), p < 0.001 in the PCA, respectively]. The τ obtained using the PFF model was significantly longer from the MCA than from the PCA, p = 0.004. No difference was found in the τ when calculated using the CFF model. Longer τ from the MCA might be related to the higher C_a of the MCA than that of the PCA. Our results demonstrate MCA-PCA differences in the τ, but only when the PFF model was applied.

Effect sizes of BOLD CVR, resting-state signal fluctuations and time delay measures for the assessment of hemodynamic impairment in carotid occlusion patients

Background and purpose

The BOLD signal amplitude as a response to a hypercapnia stimulus is commonly used to assess cerebrovascular reserve. Despite recent advances, the implementation remains cumbersome and alternative ways to assess hemodynamic impairment are desirable. Resting-state BOLD signal fluctuations (rsBOLD) have been proposed however data on its sensitivity and dependence on baseline venous cerebral blood volume (vCBV) is limited. The primary aim of this study was to compare the effect sizes of resting-state and hypercapnia induced BOLD signal changes in the detection of hemodynamic impairment. The second aim of the study was to assess the dependence of BOLD signal variability on vCBV.

Materials and methods

Fifteen patients with internal carotid artery occlusive disease and 15 matched healthy controls were included in this study. The BOLD signal was derived from a dual-echo gradient-echo echo-planar sequence during hypercapnia (HC) and hyperoxia (HO) gas modulations. BOLD (fractional) amplitude of low frequency fluctuations ((f)ALFF) was compared to HC-BOLD, BOLD response delays derived from time delay analysis and ΔBOLD in response to progressively increasing HC. Effect sizes (i.e. the standard mean difference between patients and controls) were calculated. HO-BOLD was used to estimate vCBV, and its contribution to the variability in rsBOLD signal was evaluated.

Results

The effect sizes of ALFF and fALFF (0.61 and 0.72) were lower than the effect sizes related to hypercapnia-based hemodynamic assessment analysis; 1.62, 1.56 and 0.90 for HC-BOLD, BOLD response delays and ΔBOLD in response to progressively increasing HC. A moderate relation was found between (f)ALFF and HC-BOLD in controls (R² of 0.61 and 0.42), but this relation decreased in patients (R² of 0.33 and 0.15). (f)ALFF did not differ between patients and controls whereas HC-BOLD did (p < 0.005). The ΔBOLD response to progressively increasing HC was significantly different in between patients and controls for ΔEtCO₂ values ≥ 2 mmHg (at +2  mmHg F(1, 18) = 5.85, p = 0.026). Up to 31% and 53% of the variance in the ALFF and HC-BOLD spatial distribution could be explained by HO-BOLD.

Conclusion

ALFF and fALFF demonstrated a moderate effect size to detect hemodynamic impairment whereas the effect size was large for methods employing a hypercapnia-based vascular stress stimulus. Based on our analysis of BOLD signal change as a response to a progressively increasing hypercapnia stimulus we can argue that a hypercapnia stimulus of at least 2 mmHg above baseline EtCO₂ is necessary to evaluate hemodynamic impairment. We also demonstrated that a substantial amount of information imbedded in the rsBOLD and HC-BOLD was explained by HO-BOLD. HO-BOLD can serve as a proxy for vCBV and this thus indicates that one should be careful when adopting these techniques in disease cases with compromised CBV.

Contribution of systemic vascular effects to fMRI activity in white matter

To investigate a potential contribution of systemic physiology to recently reported BOLD fMRI signals in white matter, we compared photo-plethysmography (PPG) and whole-brain fMRI signals recorded simultaneously during long resting-state scans from an overnight sleep study. We found that intermittent drops in the amplitude of the PPG signal exhibited strong and widespread correlations with the fMRI signal, both in white matter (WM) and in gray matter (GM). The WM signal pattern resembled that seen in previous resting-state fMRI studies and closely tracked the location of medullary veins. Its temporal cross-correlation with the PPG amplitude was bipolar, with an early negative value. In GM, the correlation was consistently positive. Consistent with previous studies comparing physiological signals with fMRI, these findings point to a systemic vascular contribution to WM fMRI signals. The PPG drops are interpreted as systemic vasoconstrictive events, possibly related to intermittent increases in sympathetic tone related to fluctuations in arousal state. The counter-intuitive polarity of the WM signal is explained by long blood transit times in the medullary vasculature of WM, which cause blood oxygenation loss and a substantial timing mismatch between blood volume and blood oxygenation effects. A similar mechanism may explain previous findings of negative WM signals around large draining veins during both task- and resting-state fMRI.

Cerebrovascular reactivity (CVR) MRI with CO2 challenge: A technical review

Cerebrovascular reactivity (CVR) is an indicator of cerebrovascular reserve and provides important information about vascular health in a range of brain conditions and diseases. Unlike steady-state vascular parameters, such as cerebral blood flow (CBF) and cerebral blood volume (CBV), CVR measures the ability of cerebral vessels to dilate or constrict in response to challenges or maneuvers. Therefore, CVR mapping requires a physiological challenge while monitoring the corresponding hemodynamic changes in the brain. The present review primarily focuses on methods that use CO2 inhalation as a physiological challenge while monitoring changes in hemodynamic MRI signals. CO2 inhalation has been increasingly used in CVR mapping in recent literature due to its potency in causing vasodilation, rapid onset and cessation of the effect, as well as advances in MRI-compatible gas delivery apparatus. In this review, we first discuss the physiological basis of CVR mapping using CO2 inhalation. We then review the methodological aspects of CVR mapping, including gas delivery apparatus, the timing paradigm of the breathing challenge, the MRI imaging sequence, and data analysis. In addition, we review alternative approaches for CVR mapping that do not require CO2 inhalation.

Clinical vascular imaging in the brain at 7T

Stroke and related cerebrovascular diseases are a major cause of mortality and disability. Even at standard-field-strengths (1.5T), MRI is by far the most sensitive imaging technique to detect acute brain infarctions and to characterize incidental cerebrovascular lesions, such as white matter hyperintensities, lacunes and microbleeds. Arterial time-of-flight (TOF) MR angiography (MRA) can depict luminal narrowing or occlusion of the major brain feeding arteries, and this without the need for contrast administration. Compared to 1.5T MRA, the use of high-field strength (3T) and even more so ultra-high-field strengths (7T), enables the visualization of the lumen of much smaller intracranial vessels, while adding a contrast agent to TOF MRA at 7T may enable the visualization of even more distal arteries in addition to veins and venules. Moreover, with 3T and 7T, the arterial vessel walls beyond the circle of Willis become visible with high-resolution vessel wall imaging. In addition, with 7T MRI, the brain parenchyma can now be visualized on a submillimeter scale. As a result, high-resolution imaging studies of the brain and its blood supply at 7T have generated new concepts of different cerebrovascular diseases. In the current article, we will discuss emerging clinical applications and future directions of vascular imaging in the brain at 7T MRI.

Effects of −12° head-down tilt with and without elevated levels of CO2 on cognitive performance: the SPACECOT study

Microgravity and elevated levels of CO2 are two common environmental stressors in spaceflight that may affect cognitive performance of astronauts. In this randomized, double-blind, crossover trial (SPACECOT), 6 healthy males (mean ± SD age: 41 ± 5 yr) were exposed to 0.04% (ambient air) and 0.5% CO2 concentrations during 26.5-h periods of −12° head-down tilt (HDT) bed rest with a 1-wk washout period between exposures. Subjects performed the 10 tests of the Cognition Test Battery before and on average 0.1, 5.2, and 21.0 h after the initiation of HDT bed rest. HDT in ambient air induced a change in response strategy, with increased response speed (+0.19 SD; P = 0.0254) at the expense of accuracy (−0.19 SD; P = 0.2867), resulting in comparable cognitive efficiency. The observed effects were small and statistically significant for cognitive speed only. However, even small declines in accuracy can potentially cause errors during mission-critical tasks in spaceflight. Unexpectedly, exposure to 0.5% CO2 reversed the response strategy changes observed under HDT in ambient air. This was possibly related to hypercapnia-induced cerebrovascular reactivity that favors cortical regions in general and the frontal cortex in particular, or to the CNS arousing properties of mildly to moderately increased CO2 levels. There were no statistically significant time-in-CO2 effects for any cognitive outcome. The small sample size and the small effect sizes are major limitations of this study and its findings. The results should not be generalized beyond the group of investigated subjects until they are confirmed by adequately powered follow-up studies.

NEW & NOTEWORTHY Simulating microgravity with exposure to 21 h of −12° head-down tilt bed rest caused a change in response strategy on a range of cognitive tests, with a statistically significant increase in response speed at the expense of accuracy. Cognitive efficiency was not affected. The observed speed-accuracy tradeoff was small but may nevertheless be important for mission-critical tasks in spaceflight. Importantly, the change in response strategy was reversed by increasing CO2 concentrations to 0.5%.

Vascular reactivity in small cerebral perforating arteries with 7 T phase contrast MRI - A proof of concept study

Existing cerebrovascular reactivity (CVR) techniques assess flow reactivity in either the largest cerebral vessels or at the level of the parenchyma. We examined the ability of 2D phase contrast MRI at 7 T to measure CVR in small cerebral perforating arteries.

Blood flow velocity in perforators was measured in 10 healthy volunteers (mean age 26 years) using a 7 T MR scanner, using phase contrast acquisitions in the semioval center (CSO), the basal ganglia (BG) and the middle cerebral artery (MCA). Changes in flow velocity in response to a hypercapnic breathing challenge were assessed, and expressed as the percentual increase of flow velocity as a function of the increase in end tidal partial pressure of CO₂.

The hypercapnic challenge increased (fit ± standard error) flow velocity by 0.7 ± 0.3%/mmHg in the CSO (P < 0.01). Moreover, the number of detected perforators (mean [range]) increased from 63 [27–88] to 108 [61–178] (P < 0.001). In the BG, the hypercapnic challenge increased flow velocity by 1.6 ± 0.5%/mmHg (P < 0.001), and the number of detected perforators increased from 48 [24–66] to 63 [32–91] (P < 0.01). The flow in the MCA increased by 5.2 ± 1.4%/mmHg (P < 0.01).

Small vessel specific reactivity can now be measured in perforators of the CSO and BG, using 2D phase contrast at 7 T.

•

We show that 2D phase contrast at 7T MRI is capable of measuring reactivity in small cerebral perforating arteries.

•

Reactivity to hypercapnia was measured in perforating arteries of the semi-oval center and the basal ganglia.

•

Both blood flow velocity and the number of detected perforating arteries increased during hypercapnia.

•

The proposed method bridges the gap between current reactivity measurements in parenchyma and large arteries.

Patient-Specific Alterations in CO2 Cerebrovascular Responsiveness in Acute and Sub-Acute Sports-Related Concussion

Background

Preliminary studies suggest that sports-related concussion (SRC) is associated with alterations in cerebral blood flow (CBF) regulation. Here, we use advanced magnetic resonance imaging (MRI) techniques to measure CBF and cerebrovascular responsiveness (CVR) in individual SRC patients and healthy control subjects.

Methods

15 SRC patients (mean age = 16.3, range 14–20 years) and 27 healthy control subjects (mean age = 17.6, range 13–21 years) underwent anatomical MRI, pseudo-continuous arterial spin labeling (pCASL) MRI and model-based prospective end-tidal targeting (MPET) of CO₂ during blood oxygenation level-dependent (BOLD) MRI. Group differences in global mean resting CBF were examined. Voxel-by-voxel group and individual differences in regional CVR were examined using statistical parametric mapping (SPM). Leave-one-out receiver operating characteristic curve analysis was used to evaluate the utility of brain MRI CO₂ stress testing biomarkers to correctly discriminate between SRC patients and healthy control subjects.

Results

All studies were tolerated with no complications. Traumatic structural findings were identified in one SRC patient. No significant group differences in global mean resting CBF were observed. There were no significant differences in the CO₂ stimulus and O₂ targeting during BOLD MRI. Significant group and patient-specific differences in CVR were observed with SRC patients demonstrating a predominant pattern of increased CVR. Leave-one-out ROC analysis for voxels demonstrating a significant increase in CVR was found to reliably discriminate between SRC patients and healthy control subjects (AUC of 0.879, p = 0.0001). The optimal cutoff for increased CVR declarative for SRC was 1,899 voxels resulting in a sensitivity of 0.867 and a specificity of 0.778 for this specific ROC analysis. There was no correlation between abnormal voxel counts and Postconcussion Symptom Scale scores among SRC patients.

Conclusion

Acute and subacute SRCs are associated with alterations in CVR that can be reliably detected by brain MRI CO₂ stress testing in individual patients.

A novel perspective to calibrate temporal delays in cerebrovascular reactivity using hypercapnic and hyperoxic respiratory challenges

Redistribution of blood flow across different brain regions, arising from the vasoactive nature of hypercapnia, can introduce errors when examining cerebrovascular reactivity (CVR) response delays. In this study, we propose a novel analysis method to characterize hemodynamic delays in the blood oxygen level dependent (BOLD) response to hypercapnia, and hyperoxia, as a way to provide insight into transient differences in vascular reactivity between cortical regions, and across tissue depths. A pseudo-continuous arterial spin labeling sequence was used to acquire BOLD and cerebral blood flow simultaneously in 19 healthy adults (12 F; 20 ± 2 years) during boxcar CO₂ and O₂ gas inhalation paradigms. Despite showing distinct differences in hypercapnia-induced response delay times (P < 0.05; Bonferroni corrected), grey matter regions showed homogenous hemodynamic latencies (P > 0.05) once calibrated for bolus arrival time derived using non-vasoactive hyperoxic gas challenges. Longer hypercapnic temporal delays were observed as the depth of the white matter tissue increased, although no significant differences in response lag were found during hyperoxia across tissue depth, or between grey and white matter. Furthermore, calibration of hypercapnic delays using hyperoxia revealed that deeper white matter layers may be more prone to dynamic redistribution of blood flow, which introduces response lag times ranging between 1 and 3 s in healthy subjects. These findings suggest that the combination of hypercapnic and hyperoxic gas-inhalation MRI can be used to distinguish between differences in CVR that arise as a result of delayed stimulus arrival time (due to the local architecture of the cerebrovasculature), or preferential blood flow distribution. Calibrated response delays to hypercapnia provide important insights into cerebrovascular physiology, and may be used to correct response delays associated with vascular impairment.

Iterative analysis of cerebrovascular reactivity dynamic response by temporal decomposition

OBJECTIVE

To improve quantitative cerebrovascular reactivity (CVR) measurements and CO ₂ arrival times, we present an iterative analysis capable of decomposing different temporal components of the dynamic carbon dioxide- Blood Oxygen-Level Dependent (CO ₂-BOLD) relationship.

EXPERIMENTAL DESIGN

Decomposition of the dynamic parameters included a redefinition of the voxel-wise CO ₂ arrival time, and a separation from the vascular response to a stepwise increase in CO ₂ (Delay to signal Plateau - DTP) and a decrease in CO ₂ (Delay to signal Baseline -DTB). Twenty-five (normal) datasets, obtained from BOLD MRI combined with a standardized pseudo-square wave CO ₂ change, were co-registered to generate reference atlases for the aforementioned dynamic processes to score the voxel-by-voxel deviation probability from normal range. This analysis is further illustrated in two subjects with unilateral carotid artery occlusion using these reference atlases.

PRINCIPAL OBSERVATIONS

We have found that our redefined CO ₂ arrival time resulted in the best data fit. Additionally, excluding both dynamic BOLD phases (DTP and DTB) resulted in a static CVR, that is maximal response, defined as CVR calculated only over a normocapnic and hypercapnic calibrated plateau.

CONCLUSION

Decomposition and novel iterative modeling of different temporal components of the dynamic CO ₂-BOLD relationship improves quantitative CVR measurements.

The role of vascular resistance in BOLD responses to progressive hypercapnia

The ability of the cerebral vasculature to regulate vascular diameter, hence resistance and cerebral blood flow (CBF), in response to metabolic demands (neurovascular coupling), and perfusion pressure changes (autoregulation) may be assessed by measuring the CBF response to carbon dioxide (CO₂ ). In healthy individuals, the CBF response to a ramp CO₂ stimulus from hypocapnia to hypercapnia is assumed sigmoidal or linear. However, other response patterns commonly occur, especially in individuals with cerebrovascular disease, and these remain unexplained. CBF responses to CO₂ in a vascular region are determined by the combined effects of the innate vascular responses to CO₂ and the local perfusion pressure; the latter ensuing from pressure-flow interactions within the cerebral vascular network. We modeled this situation as two vascular beds perfused in parallel from a fixed resistance source. Our premise is that all vascular beds have a sigmoidal reduction of resistance in response to a progressive rise in CO₂ . Surrogate CBF data to test the model was provided by magnetic resonance imaging of blood oxygen level-dependent (BOLD) signals. The model successfully generated all the various BOLD-CO₂ response patterns, providing a physiological explanation of CBF distribution as relative differences in the network of vascular bed resistance responses to CO₂ . Hum Brain Mapp 38:5590-5602, 2017. © 2017 Wiley Periodicals, Inc.

Assessing cerebrovascular reactivity by the pattern of response to progressive hypercapnia

Cerebral blood flow responds to a carbon dioxide challenge, and is often assessed as cerebrovascular reactivity, assuming a linear response over a limited stimulus range or a sigmoidal response over a wider range. However, these assumed response patterns may not necessarily apply to regions with pathophysiology. Deviations from sigmoidal responses are hypothesised to result from upstream flow limitations causing competition for blood flow between downstream regions, particularly with vasodilatory stimulation; flow is preferentially distributed to regions with more reactive vessels. Under these conditions, linear or sigmoidal fitting may not fairly describe the relationship between stimulus and flow. To assess the range of response patterns and their prevalence a survey of healthy control subjects and patients with cerebrovascular disease was conducted. We used a ramp carbon dioxide challenge from hypo- to hypercapnia as the stimulus, and magnetic resonance imaging to measure the flow responses. We categorized BOLD response patterns into four types based on the signs of their linear slopes in the hypo- and hypercapnic ranges, color coded and mapped them onto their respective anatomical scans. We suggest that these type maps complement maps of linear cerebrovascular reactivity by providing a better indication of the actual response patterns. Hum Brain Mapp 38:3415-3427, 2017. © 2017 Wiley Periodicals, Inc.

7-T MRI in Cerebrovascular Diseases: Challenges to Overcome and Initial Results

Magnetic resonance imaging (MRI) plays a key role in the investigation of cerebrovascular diseases. Compared with computed tomography (CT) and digital subtraction angiography (DSA), its advantages in diagnosing cerebrovascular pathology include its superior tissue contrast, its ability to visualize blood vessels without the use of a contrast agent, and its use of magnetic fields and radiofrequency pulses instead of ionizing radiation. In recent years, ultrahigh field MRI at 7 tesla (7 T) has shown promise in the diagnosis of many cerebrovascular diseases. The increased signal-to-noise ratio (SNR; 2.3x and 4.7x increase compared with 3 and 1.5 T, respectively) and contrast-to-noise ratio (CNR) at this higher field strength can be exploited to obtain a higher spatial resolution and higher lesion conspicuousness, enabling assessment of smaller brain structures and lesions. Cerebrovascular diseases can be assessed at different tissue levels; for instance, changes of the arteries feeding the brain can be visualized to determine the cause of ischemic stroke, regional changes in brain perfusion can be mapped to predict outcome after revascularization, and tissue damage, including old and recent ischemic infarcts, can be evaluated as a marker of ischemic burden. For the purpose of this review, we will discriminate 3 levels of assessment of cerebrovascular diseases using MRI: Pipes, Perfusion, and Parenchyma (3 Ps). The term Pipes refers to the brain-feeding arteries from the heart and aortic arch, upwards to the carotid arteries, vertebral arteries, circle of Willis, and smaller intracranial arterial branches. Perfusion is the amount of blood arriving at the brain tissue level, and includes the vascular reserve and perfusion territories. Parenchyma refers to the acute and chronic burden of brain tissue damage, which includes larger infarcts, smaller microinfarcts, and small vessel disease manifestations such as white matter lesions, lacunar infarcts, and microbleeds. In this review, we will describe the key developments in the last decade of 7-T MRI of cerebrovascular diseases, subdivided for these 3 levels of assessment.

Cerebral microvascular blood flow and CO2 reactivity in pulmonary arterial hypertension

Hypocapnia and endothelial dysfunction might impair microvascular cerebral blood flow (CBF_micr) and cerebrovascular reactivity to CO₂ (CVR_CO2). Pulmonary arterial hypertension (PAH) is characteristically associated with chronic alveolar hyperventilation and microvascular endothelial dysfunction. We therefore determined CBF_micr (pre-frontal blood flow index (BFI) by the indocyanine green-near infrared spectroscopy methodology) during hypocapnia and hypercapnia in 25 PAH patients and 10 gender- and age-matched controls. Cerebral BFI was lower in patients than controls at similar transcutaneous PCO₂ (PtcCO₂) levels in both testing conditions. In fact, while BFI increased from hypocapnia to hypercapnia in all controls, it failed to increase in 17/25 (68%) patients. Thus, BFI increased to a lesser extent from hypo to hypercapnia ("Δ") in patients, i.e., they showed lower Δ BFI/Δ PtcCO₂ ratios than controls. In conclusion, CBF_micr and CVR_CO2 are lessened in clinically stable, mildly-impaired patients with PAH. These abnormalities might be associated with relevant clinical outcomes (hyperventilation and dyspnea, cognition, cerebrovascular disease) being potentially amenable to pharmacological treatment.

Longitudinal Brain Magnetic Resonance Imaging CO2 Stress Testing in Individual Adolescent Sports-Related Concussion Patients: A Pilot Study

Background

Advanced neuroimaging studies in concussion have been limited to detecting group differences between concussion patients and healthy controls. In this small pilot study, we used brain magnetic resonance imaging (MRI) CO₂ stress testing to longitudinally assess cerebrovascular responsiveness (CVR) in individual sports-related concussion (SRC) patients.

Methods

Six SRC patients (three males and three females; mean age = 15.7, range = 15–17 years) underwent longitudinal brain MRI CO₂ stress testing using blood oxygen level-dependent (BOLD) MRI and model-based prospective end-tidal CO₂ targeting under isoxic conditions. First-level and second-level comparisons were undertaken using statistical parametric mapping (SPM) to score the scans and compare them to an atlas of 24 healthy control subjects.

Results

All tests were well tolerated and without any serious adverse events. Anatomical MRI was normal in all study participants. The CO₂ stimulus was consistent between the SRC patients and control subjects and within SRC patients across the longitudinal study. Individual SRC patients demonstrated both quantitative and qualitative patient-specific alterations in CVR (p < 0.005) that correlated strongly with clinical findings, and that persisted beyond clinical recovery.

Conclusion

Standardized brain MRI CO₂ stress testing is capable of providing a longitudinal assessment of CVR in individual SRC patients. Consequently, larger prospective studies are needed to examine the utility of brain MRI CO₂ stress testing as a clinical tool to help guide the evaluation, classification, and longitudinal management of SRC patients.