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

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.

Hemodynamic imaging parameters in brain metastases patients - Agreement between multi-delay ASL and hypercapnic BOLD

Arterial spin labeling (ASL) MRI is a routine clinical imaging technique that provides quantitative cerebral blood flow (CBF) information. A related technique is blood oxygenation level-dependent (BOLD) MRI during hypercapnia, which can assess cerebrovascular reactivity (CVR). ASL is weighted towards arteries, whereas BOLD is weighted towards veins. Their associated parameters in heterogeneous tissue types or under different hemodynamic conditions remains unclear. Baseline multi-delay ASL MRI and BOLD MRI during hypercapnia were performed in fourteen patients with brain metastases. In the ROI analysis, the CBF and CVR values were positively correlated in regions showing sufficient reserve capacity (i.e. non-steal regions, rrm = 0.792). Additionally, longer hemodynamic lag times were related to lower baseline CBF (rrm = -0.822) and longer arterial arrival time (AAT; rrm = 0.712). In contrast, in regions exhibiting vascular steal an inverse relationship was found with higher baseline CBF related to more negative CVR (rrm = -0.273). These associations were confirmed in voxelwise analyses. The relationship between CBF, AAT and CVR measures seems to be dependent on the vascular status of the underlying tissue. Healthy tissue relationships do not hold in tissues experiencing impaired or exhausted autoregulation. CVR metrics can possibly identify at-risk areas before perfusion deficiencies become visible on ASL MRI, specifically within vascular steal regions.

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.

Generating dynamic carbon-dioxide traces from respiration-belt recordings: Feasibility using neural networks and application in functional magnetic resonance imaging

INTRODUCTION

In the context of functional magnetic resonance imaging (fMRI), carbon dioxide (CO2) is a well-known vasodilator that has been widely used to monitor and interrogate vascular physiology. Moreover, spontaneous fluctuations in end-tidal carbon dioxide (PETCO2) reflects changes in arterial CO2 and has been demonstrated as the largest physiological noise source for denoising the low-frequency range of the resting-state fMRI (rs-fMRI) signal. However, the majority of rs-fMRI studies do not involve CO2 recordings, and most often only heart rate and respiration are recorded. While the intrinsic link between these latter metrics and CO2 led to suggested possible analytical models, they have not been widely applied.

METHODS

In this proof-of-concept study, we propose a deep-learning (DL) approach to reconstruct CO2 and PETCO2 data from respiration waveforms in the resting state.

RESULTS

We demonstrate that the one-to-one mapping between respiration and CO2 recordings can be well predicted using fully convolutional networks (FCNs), achieving a Pearson correlation coefficient (r) of 0.946 ± 0.056 with the ground truth CO2. Moreover, dynamic PETCO2 can be successfully derived from the predicted CO2, achieving r of 0.512 ± 0.269 with the ground truth. Importantly, the FCN-based methods outperform previously proposed analytical methods. In addition, we provide guidelines for quality assurance of respiration recordings for the purposes of CO2 prediction.

DISCUSSION

Our results demonstrate that dynamic CO2 can be obtained from respiration-volume using neural networks, complementing the still few reports in DL of physiological fMRI signals, and paving the way for further research in DL based bio-signal processing.

Hemodynamic and metabolic changes during hypercapnia with normoxia and hyperoxia using pCASL and TRUST MRI in healthy adults

Blood oxygenation level-dependent (BOLD) or arterial spin labeling (ASL) MRI with hypercapnic stimuli allow for measuring cerebrovascular reactivity (CVR). Hypercapnic stimuli are also employed in calibrated BOLD functional MRI for quantifying neuronally-evoked changes in cerebral oxygen metabolism (CMRO2). It is often assumed that hypercapnic stimuli (with or without hyperoxia) are iso-metabolic; increasing arterial CO2 or O2 does not affect CMRO2. We evaluated the null hypothesis that two common hypercapnic stimuli, 'CO2 in air' and carbogen, are iso-metabolic. TRUST and ASL MRI were used to measure the cerebral venous oxygenation and cerebral blood flow (CBF), from which the oxygen extraction fraction (OEF) and CMRO2 were calculated for room-air, 'CO2 in air' and carbogen. As expected, CBF significantly increased (9.9% ± 9.3% and 12.1% ± 8.8% for 'CO2 in air' and carbogen, respectively). CMRO2 decreased for 'CO2 in air' (-13.4% ± 13.0%, p < 0.01) compared to room-air, while the CMRO2 during carbogen did not significantly change. Our findings indicate that 'CO2 in air' is not iso-metabolic, while carbogen appears to elicit a mixed effect; the CMRO2 reduction during hypercapnia is mitigated when including hyperoxia. These findings can be important for interpreting measurements using hypercapnic or hypercapnic-hyperoxic (carbogen) stimuli.

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.

Normal BOLD Response to a Step CO2 Stimulus After Correction for Partial Volume Averaging

Cerebrovascular reactivity (CVR) is defined as the change in cerebral blood flow induced by a change in a vasoactive stimulus. CVR using BOLD MRI in combination with changes in end-tidal CO₂ is a very useful method for assessing vascular performance. In recent years, this technique has benefited from an advanced gas delivery method where end-tidal CO₂ can be targeted, measured very precisely, and validated against arterial blood gas sampling (Ito et al., 2008). This has enabled more precise comparison of an individual patient against a normative atlas of healthy subjects. However, expected control ranges for CVR metrics have not been reported in the literature. In this work, we calculate and report the range of control values for the magnitude (mCVR), the steady state amplitude (ssCVR), and the speed (TAU) of the BOLD response to a standard step stimulus, as well as the time delay (TD) as observed in a cohort of 45 healthy controls. These CVR metrics maps were corrected for partial volume averaging for brain tissue types using a linear regression method to enable more accurate quantitation of CVR metrics. In brief, this method uses adjacent voxel CVR metrics in combination with their tissue composition to write the corresponding set of linear equations for estimating CVR metrics of gray matter (GM), white matter (WM), and cerebrospinal fluid (CSF). After partial volume correction, mCVR and ssCVR increase as expected in gray matter, respectively, by 25 and 19%, and decrease as expected in white matter by 33 and 13%. In contrast, TAU and TD decrease in gray matter by 33 and 13%. TAU increase in white matter by 24%, but TD surprisingly decreased by 9%. This correction enables more accurate voxel-wise tissue composition providing greater precision when reporting gray and white matter CVR values.

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.

Changes in volumetric and metabolic parameters relate to differences in exposure to sub-concussive head impacts

Structural and calibrated magnetic resonance imaging data were acquired on 44 collegiate football players prior to the season (PRE), following the first four weeks in-season (PTC) and one month after the last game (POST). Exposure data collected from g-Force accelerometers mounted to the helmet of each player were used to split participants into HIGH (N = 22) and LOW (N = 22) exposure groups, based on the frequency of impacts sustained by each athlete. Significant decreases in grey-matter volume specific to the HIGH group were documented at POST (P = 0.009), compared to baseline. Changes in resting cerebral blood flow (CBF0), corrected for partial volume effects, were observed within the HIGH group, throughout the season (P < 0.0001), suggesting that alterations in perfusion may follow exposure to sub-concussive collisions. Co-localized significant increases in cerebral metabolic rate of oxygen consumption (CMRO2|0) mid-season were also documented in the HIGH group, with respect to both PRE- and POST values. No physiological changes were observed in the LOW group. Therefore, cerebral metabolic demand may be elevated in players with greater exposure to head impacts. These results provide novel insight into the effects of sub-concussive collisions on brain structure and cerebrovascular physiology and emphasize the importance of multi-modal imaging for a complete characterization of cerebral health.

Multi-modal normalization of resting-state using local physiology reduces changes in functional connectivity patterns observed in mTBI patients

Blood oxygenation level dependent (BOLD) resting-state functional magnetic resonance imaging (rs-fMRI) may serve as a sensitive marker to identify possible changes in the architecture of large-scale networks following mild traumatic brain injury (mTBI). Differences in functional connectivity (FC) measurements derived from BOLD rs-fMRI may however be confounded by changes in local cerebrovascular physiology and neurovascular coupling mechanisms, without changes in the underlying neuronally driven connectivity of networks. In this study, multi-modal neuroimaging data including BOLD rs-fMRI, baseline cerebral blood flow (CBF0) and cerebrovascular reactivity (CVR; acquired using a hypercapnic gas breathing challenge) were collected in 23 subjects with reported mTBI (14.6±14.9 months post-injury) and 27 age-matched healthy controls. Despite no group differences in CVR within the networks of interest (P > 0.05, corrected), significantly higher CBF0 was documented in the mTBI subjects (P < 0.05, corrected), relative to the controls. A normalization method designed to account for differences in CBF0 post-mTBI was introduced to evaluate the effects of such an approach on reported group differences in network connectivity. Inclusion of regional perfusion measurements in the computation of correlation coefficients within and across large-scale networks narrowed the differences in FC between the groups, suggesting that this approach may elucidate unique changes in connectivity post-mTBI while accounting for shared variance with CBF0. Altogether, our results provide a strong paradigm supporting the need to account for changes in physiological modulators of BOLD in order to expand our understanding of the effects of brain injury on large-scale FC of cortical networks.

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

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

Multi-parametric analysis reveals metabolic and vascular effects driving differences in BOLD-based cerebrovascular reactivity associated with a history of sport concussion

Objective: Identify alterations in cerebrovascular reactivity (CVR) based on the history of sport-related concussion (SRC). Further explore possible mechanisms underlying differences in vascular physiology using hemodynamic parameters modeled using calibrated magnetic resonance imaging (MRI). Method: End-tidal targeting and dual-echo MRI were combined to probe hypercapnic and hyperoxic challenges in athletes with (n = 32) and without (n = 31) a history of SRC. Concurrent blood oxygenation level dependent (BOLD) and arterial spin labeling (ASL) data were used to compute BOLD-CVR, ASL-CVR, and other physiological parameters including resting oxygen extraction fraction (OEF0) and cerebral blood volume (CBV0). Multiple linear and logistic regressions were then used to identify dominant parameters driving group-differences in BOLD-CVR. Results: Robust evidence for elevated BOLD-CVR were found in athletes with SRC history spreading over parts of the cortical hemispheres. Follow-up analyses showed co-localized differences in ASL-CVR (representing modulation of cerebral blood flow) and hemodynamic factors representing static vascular (i.e., CBV0) and metabolic (i.e., OEF0) effects suggesting that group-based differences in BOLD-CVR may be driven by a mixed effect from factors with vascular and metabolic origins. Conclusion: These results emphasize that while BOLD-CVR offers promises as a surrogate non-specific biomarker for cerebrovascular health following SRC, multiple hemodynamic parameters can affect its relative measurements. Abbreviations: [dHb]: concentration of deoxyhemoglobin; AFNI: Analysis of Functional NeuroImages ( https://afni.nimh.nih.gov ); ASL: arterial spin labeling; BIG: position group: defensive and offensive linemen; BIG-SKILL: position group: full backs, linebackers, running backs, tight-ends; BOLD: blood oxygen level dependent; CBF: cerebral blood flow; CMRO2: cerebral metabolic rate of oxygen consumption; CTL: group of control subjects; CVR: cerebrovascular reactivity; fMRI: functional magnetic resonance imaging; FSL: FMRIB software library ( https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/ ); HC: hypercapnia; HO: hyperoxia; HX: group with history of concussion; M: maximal theoretical BOLD signal upon complete removal of venous dHb; pCASL: pseudo-continuous arterial spin labeling; PETCO2: end-tidal carbon dioxide; PETO2: end-tidal oxygen; SCAT: sport-concussion assessment tool; SKILL: position group: defensive backs, kickers, quarterbacks, safeties, wide-receivers; SRC: sport-related concussion.

Resting CMRO2 fluctuations show persistent network hyper-connectivity following exposure to sub-concussive collisions

Exposure to head impacts may alter brain connectivity within cortical hubs such as the default-mode network (DMN). However, studies have yet to consider the confounding effects of altered resting cerebral blood flow (CBF0) and cerebrovascular reactivity (CVR) on changes in connectivity following sub-concussive impacts. Here, 23 Canadian collegiate football players were followed during a season using calibrated resting-state MRI and helmet accelerometers to examine the interplay between the neural and vascular factors that determine functional connectivity (FC). Connectivity-based analyses using blood oxygen level dependent (BOLD) and cerebral metabolic rate of oxygen consumption (CMRO2) mapping were used to study the DMN longitudinally. Network-specific decreases in CBF0 were observed one month following the season, while impaired CVR was documented at both mid-season and one month following the season, compared to pre-season baseline. Alterations in CBF0 and BOLD-based CVR throughout the season suggest that neurophysiological markers may show different susceptibility timelines following head impacts. DMN connectivity was increased throughout the season, independent of changes in cerebrovascular physiology, suggesting that alterations in FC following sub-concussive impacts are robust and independent of changes in brain hemodynamics. No significant correlations between impact kinematics and DMN connectivity changes were documented in this study. Altogether, these findings create a strong paradigm for future studies to examine the underlying neural and vascular mechanisms associated with increases in network connectivity following repeated exposure to sub-concussive collisions, in an effort to improve management of head impacts in contact sports.

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.