- 1Department of Medical Biophysics, University of Toronto, Toronto, ON, Canada
- 2Hurvitz Brain Sciences Program, Sunnybrook Research Institute, Toronto, ON, Canada
- 3Dr. Sandra Black Centre for Brain Resilience and Recovery, Toronto, ON, Canada
- 4Institute of Medical Sciences, University of Toronto, Toronto, ON, Canada
- 5Division of Neuroradiology, Joint Department of Medical Imaging, University Health Network, Toronto, ON, Canada
- 6Department of Medical Imaging, University of Toronto, Toronto, ON, Canada
- 7Brain Research Imaging Centre, Centre for Clinical Brain Sciences, UK Dementia Research Institute Centre, The University of Edinburgh, Edinburgh, United Kingdom
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 CO2 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 (R2 = 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.
1. Introduction
The cerebrovascular network is tightly regulated to maintain appropriate cerebral blood flow in the presence of systemic changes in blood pressure (1). Changes in blood pressure, neuronal activity, or arterial blood gases can elicit strong cerebrovascular responses and changes in vascular tone. Increasing the fraction of inhaled CO2 is a common cerebrovascular reactivity (CVR) challenge, as it increases the partial pressure of CO2 in arterial blood. There are two distinct transition phases that occur during a CVR challenge. The transition from a normocapnic to hypercapnic state is facilitated by vasodilation. During hypercapnia there is a decrease in intravascular pH, causing a cascade of chemical reactions in the endothelial and smooth muscle components of the vessel wall that increase the vessel lumen diameter, resulting in an influx of oxygenated blood. Upon cessation of the hypercapnic stimulus (i.e., returning to room air), a reversal effect occurs where blood vessels transition back to a physiological baseline. This process is induced by vasoconstriction, a physiological component of the hypercapnic response, whereby a passive decrease in arterial CO2 occurs once the hypercapnic stimulus ends. A resultant CO2 washout is promoted through a negative feedback loop where a decrease in local CO2 results in a smaller vessel diameter (2). These dynamic structural and physiological responses to CO2 influx and efflux can be influenced by cerebrovascular properties such as tissue type and vascular resistance (3, 4).
CVR can be assessed through dynamic functional neuroimaging to indirectly measure the vascular response to a vasoactive stimulus, providing a means to visualize and quantify cerebrovascular function in-vivo (2). Metrics of CVR have been demonstrated to have utility in the characterization of conditions as diverse as stroke (5, 6), dementia (7–9), small vessel disease (10, 11), obstructive sleep apnea (12–14), and Moyamoya disease (15, 16). Dynamic responses to a CVR challenge can be imaged using blood oxygenation level dependent (BOLD) Magnetic Resonance Imaging (17–19) as well as other nuclear medicine approaches (11, 20). The contrast for BOLD-MRI images arises from changes in the local concentration of oxygenated to deoxygenated blood. What drives this shift in concentration is dependent on multiple characteristics of brain physiology, such as cerebral blood flow and cerebral blood volume. BOLD-MRI indirectly measures the longitudinal changes that occur when these parameters are manipulated by a hypercapnia challenge. This non-invasive imaging technique offers a high contrast-to-noise ratio and an adequate temporal resolution (3, 19, 21–23).
CVR is typically calculated using a general linear model (GLM) that relies on model parameters to yield a single CVR estimate. For example, the end-tidal CO2 (PETCO2) trace is typically convolved with a hemodynamic response function (HRF) in a GLM-CVR analysis. Other modeling approaches account for systematic differences between the PETCO2 and BOLD signal; these examples include a hemodynamic lag term, refined HRF, transfer functions, and non-linear analyses (24–26). However, linear methods often assume that the response pattern of the BOLD and PETCO2 signal are highly consistent. Diseases that affect the cerebrovasculature may challenge this assumption. Furthermore, the conventional CVR measure summarizes the relative change in BOLD signal for multiple challenges in succession despite potential information that could be gained from assessing each challenge individually. Given that cerebrovascular variability may reflect pathophysiology and CBF differences, we look to examine the more granular nuances of the BOLD signal response to CO2 challenges in succession (27, 28).
The primary objective of this study is to develop a non-parametric CVR metric that describes the change in BOLD signal induced by the dilatory response to hypercapnia, while circumventing the need for model-fitting, and allowing the quantification of individual CVR responses. The Theil-Sen Estimator (Sen’s Slope) is a non-parametric regression model and can be used as a CVR metric to characterize the relative change in magnitude of the BOLD signal induced by individual dilatory responses to hypercapnia. Our second objective is to characterize the BOLD signal decrease that occurs during the transition period following a hypercapnia challenge. We hypothesize that the rate of BOLD signal decline during a transition period from hypercapnia to normocapnia will be significantly associated with the non-parametric CVR measure of the following hypercapnia challenge. These metrics are extracted in two clinical populations: individuals with sleep apnea and cognitive concerns. Regional differences in response to hypercapnia are of interest in these clinical populations, hence we conduct a region-of-interest (ROI) analysis. The hippocampus is part of the physiological brain alterations seen in cognitive impairment and/or sleep apnea (29–32). Second, the insula is highly organized in its regulation of the autonomic nervous system, which includes cardiovascular responses to stimuli (33). Third, the thalamus is highly implicated in both sleep disorders and cognitive impairment (34, 35). Fourth, the centrum semiovale is relevant to CVR as it has lower cerebral blood flow compared to grey matter, it has limited collateral arterial blood supply, and is proximal to the draining medullary veins that are implicated in vascular brain dysfunctions (36–38).
2. Methods
2.1. Participants
Participants were recruited at Sunnybrook Research Institute as part of an on-going international study (10) and were either: (1) adults with Alzheimer’s Disease, Mild Cognitive Impairment, or subjective memory concerns seen at a cognitive neurology clinic, or (2) adults with sleep apnea seen at a sleep neurology clinic. A third retrospective cohort comprised individuals with memory complaints (39). All participants gave written informed consent. This research was approved by Sunnybrook Health Sciences Centre Research Ethics Board.
Group 1: Participants underwent clinical testing to ascertain Alzheimer’s disease diagnosis as per the National Institute on Ageing-Alzheimer’s Association (NIA-AA) (40), Mild Cognitive Impairment as per the Albert diagnostic guidelines (41), subcortical Vascular Cognitive Impairment, or subjective complaints. Group 2: Participants were newly diagnosed with sleep apnea based on an apnea hypopnea index ≥15 and at an oxygen desaturation index of 10 de-saturations per hour on diagnostic polysomnography, accompanied by subjective sleepiness, and have never been treated for sleep apnea. Exclusion criteria included: history of major stroke or other central nervous system diseases, use of alpha-blockers, persistent non-sinus arrhythmia, severe pulmonary or cardiac diseases including chronic obstructive pulmonary disease and congestive heart failure, waking saturation of <90%, and history of panic disorder. Group 3: Participants from a retrospective data collection at Sunnybrook Health Sciences Centre (42) were accessed and comprised of older adults with subjective memory complaints; age and sex were matched to Group 1.
For all three groups, participants were excluded based on: contraindications to MRI safety, history of significant head trauma, brain tumours, or hydrocephalus, significant drug or alcohol abuse, and current unstable cardiac disease.
2.2. MRI data collection
MRI data for Group 1 and Group 2 were collected using a 3 Tesla (3 T) Siemens Prisma system scanner (Siemens Medical Solutions, Erlangen, Germany) with a 32-channel head coil reception and body coil transmission. T1-weighted images were collected for registration with functional images and to define structural regions of interests to be segmented. Functional (i.e., BOLD-fMRI) images were collected using a simultaneous multi-slice echo planar imaging (SMS-EPI) pulse sequence (43). T2-weighted FLAIR sequence images were collected for white matter hyperintensity (WMH) quantification.
MRI data for Group 3 were collected using a 3 T Philips Achieva system scanner (Philips Healthcare, Best, Netherlands). T1-weighted and T2-weighted FLAIR images were collected on this scanner. BOLD-fMRI data were collected using a single-shot EPI. MRI parameters are summarized in Table 1.
2.3. Hypercapnia challenge
The CVR set-up used for Group 1 and Group 2, Experimental Design A, was established as part of a collaborative two-site study (10). An enriched-CO2 gas blend contained in a gas cylinder was fed via plastic tubing through the MRI waveguide to the tight-fitting face mask worn by the participant in the MRI. Gas administration was interleaved with room air “baseline” periods to produce a boxcar breathing paradigm. The challenge was administered as follows: 2 min of room air, 3 min of 6% CO2 air, 2 min of room air, 3 min of 6% CO2 air, and 2 min of room air. CO2 and O2 gas traces were monitored using two gas analyzers, O2100C and CO2100C, that are part of the BIOPAC® respiratory gas analysis modules (BIOPAC Systems, California, United States). The hypercapnic challenge for Group 3, Experimental Design B, was administered using a RespirAct Gen3 breathing circuit (Thornhill Research, Toronto, Canada) to reach a targeted 10 mmHg PETCO2 change. Unlike the first breathing paradigm, hypercapnic gas was administered in unequal timing blocks: 45 s rest, 45 s on, 90 s rest, followed by a second dilation challenge of 120 s and a final rest of 180 s. This external data source was retroactively analyzed, and differences in the hypercapnic challenge are not intentional.
2.4. BOLD data processing
BOLD data were pre-processed using parameters appropriate for the image size. Data were motion-corrected using FSL MCFLIRT (FSL, version 5.0, http://fsl.fmrib.ox.ac.uk), spatially smoothed with a FWHM Gaussian Kernel (Group 1 and 2: 3 mm, Group 3: 5 mm), a high-pass filter was applied [Group 1 and 2: 3.18 × 10−4 Hz (500 s), Group 3: 1.60 × 10−3 Hz (300 s)] to remove drift in the signal, and slice-timing correction was performed. An initial set of BOLD volumes were discarded to account for signal equilibrium (Group 1 and 2: 25 volumes, Group 3: 5 volumes).
2.5. Region-of-interest generation
Four pre-defined bilateral ROIs were identified for hypothesis testing in this study, which included the insula, thalamus, hippocampus, and centrum semiovale. The insula was segmented using a MNI152 standard space template and MNI structural atlas (44, 45). The hippocampus was segmented in native T1-space using the convolutional neural network-based tool HippMapp3r (hippmapp3r.readthedocs.io) (46). The thalamus was segmented using FIRST in FSL using the native T1 image. Finally, the bilateral centrum semiovale was defined as contiguous voxels of white matter in an MNI152 standard space axial plane that is 3 mm superior to the top of the lateral ventricles. The standard space masks were transformed into native T1-space. Bilateral ROIs in T1-space are conceptualized (Figure 1A). For each of the ROIs, the masks were down-sampled from T1-space to BOLD-space to generate BOLD-CVR time series data. All registrations were performed using FSL FLIRT (FSL, version 5.0, http://fsl.fmrib.ox.ac.uk).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1199805/fneur-14-1199805-HTML/image_m/fneur-14-1199805-g001.jpg)
Figure 1. From left to right: sagittal (posterior to anterior), coronal (right to left), and axial (right to left) views of T1-space ROIs (A). From left to right: sagittal (posterior to anterior), coronal (right to left), and axial (right to left) views of COPE-CVR map (B).
2.6. Isolation of transition phases during CVR challenge
Physiologically distinct periods were isolated according to the experimental breathing paradigm of the group (refer to Section 2.3) and are summarized in Table 2. The first and second non-parametric CVR metrics (CVR-1 and CVR-2) arise from the transition periods that describe the BOLD signal response to the first and second administration of CO2-rich air, respectively. The hypercapnia to normocapnia transition period describes the BOLD response as it returns to room air condition, specifically the first half of the rest period. These analytical time windows were applied consistently for all CVR data.
2.7. Transition from normocapnia to hypercapnia
We used Sen’s Slope as a non-parametric metric to quantify the percent change in BOLD signal when transitioning from normo-to-hypercapnia conditions; this approach circumvents the need to fit the BOLD signal to the PETCO2 signal and is robust to outlier measurements that may occur due to imaging noise. CVR-1 describes the response to the first CO2 challenge, while CVR-2 describes the response to the second challenge. BOLD signal differences were calculated between every pair of values in the time series excerpt, resulting in a total of n(n−1)/2 pairs of data. These values were then divided by the interval range corresponding to the calculated signal difference. The Sen’s Slope was based on the median rate of BOLD change from among each of the calculated values (Equation 1), followed by multiplying by the number of data points within the dilation period, dividing the mean baseline BOLD signal prior to the hypercapnia challenge, and yielding a unit-less BOLD percent signal change (Equation 2). BOLD percent signal change is then converted to a CVR metric by dividing by the change in PETCO2 (Equation 3).
Where xi and xj are BOLD signal values at times i and j (j > i).
To compare this approach against a conventional method, we calculated the coefficient of parameter estimate (COPE) values based on the FSL feat tool.1 This CVR estimator uses the voxel with the highest agreement (i.e., COPE value) between the signal and experimental paradigm to produce a CVR estimate (Figure 1B). The COPE-CVR estimate is representative of the entire hypercapnia challenge, rather than individual responses.
2.8. Transition from hypercapnia to normocapnia
The BOLD signal will tend to decrease after each hypercapnia challenge. This signal recovery is related to vasoconstriction and the transition back to normocapnia. We approximate this response as the linear rate of change in a regression model (lm function in R, http://www.R-project.org) where the BOLD signal is modeled as a function of time starting from the transition period (Equation 4). The resultant parameter estimate (𝜏) was further normalized by the change in PETCO2 to account for the influence of CO2 on the BOLD signal (Equation 5).
2.9. Simulation
We evaluated non-parametric and parametric CVR estimates through simulated hypercapnia events. The BOLD signal response to CO2 inhalation can vary depending on several factors; to simulate this we considered two (fast and slow) responses. A fast response was approximated by a sigmoidal curve for the BOLD signal change, while a slow response was simulated as a slower ramp up. Increments of simulated noise were added to each response to produce a total of 4,000 normocapnia-to-hypercapnia simulations. To characterize temporal noise characteristics, a temporal signal-to-noise (tSNR) ratio was calculated by dividing the mean baseline signal by its standard deviation (47). The tSNR was used to evaluate the influence that incremental increases of white noise added to the simulated time series would have on the CVR estimates. Both non-parametric (Sen’s) and parametric (linear regression) slopes for the simulated BOLD signal change were determined. The regressor in the parametric approach was a convolution of the gamma-HRF and the experimental paradigm. An Akaike Information Criterion (AIC) estimate was computed for both models to compare which model had the most likelihood of producing simulated scenarios.
2.10. Statistical analysis
All statistics were performed in R (version 1.2.1335). We tested our hypothesis on the empirical data by conducting two regression models for each ROI. In the first regression model (CVR Model-1), we assess the effect of the transition rate on the non-parametric CVR estimate of the second hypercapnia challenge (CVR-2), while adjusting for age, sex, clinical population, and CVR experimental design differences. For the ‘clinical population’ confounder, Group 1 and Group 3 were pooled because they are a similar clinical population (e.g., cognitive concerns), while Group 2 participants were grouped separately. For the ‘experimental design’ confounder, Group 1 and Group 2 were pooled together due to similarities in experimental design. In the second regression model (CVR Model-2) we account for similarities between hypercapnia challenges in succession by orthogonalizing the transition rate relative to CVR-2 and considering the first CVR estimate as an additional explanatory variable. A p-value was considered significant if it reached p ≤ 0.0125, to account for 4 ROIs. Standardized parameter estimates and adjusted-R2 values are reported for each predictor. Pearson correlation estimates between transition rate, CVR-1, CVR-2, and COPE-CVR are reported for the thalamus to demonstrate the utility of characterizing individual temporal features of a CVR challenge.
3. Results
3.1. Simulation
Simulation results revealed a marginal benefit of using the parametric method for fast CVR responses (Figure 2A), as reflected by lower AIC values and less variability compared to the non-parametric values. For the simulated slow response (Figure 2B), however, the non-parametric method had a marginally lower AIC than the parametric model. As expected, the fast responses had lower AIC values compared to the slow response (Figure 2C). There was a general trend of lower AIC as the tSNR increased for all simulations.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1199805/fneur-14-1199805-HTML/image_m/fneur-14-1199805-g002.jpg)
Figure 2. Example of a simulated fast (A) and simulated slow (B) response. AIC values for non-parametric and parametric models (C) are shown for a range of tSNR levels. Semi-transparent points represent individual simulation AIC results, while solid points represent the average AIC value for simulations across tSNR values. AIC was lower for the parametric model in the fast response simulation, and lower for the non-parametric model in the slow response simulation.
3.2. Empirical CVR
Cohort characteristics for the 54 participants in this study are provided in Table 3. Age, sex, years of education, Montreal Cognitive Assessment (MoCA) scores, and white matter hyperintensity (WMH) volumes are reported.
For CVR Model-1, the rate of BOLD signal decrease during the transition to normocapnia was significantly associated with the relative change in BOLD signal during a following hypercapnia event for each ROI, after accounting for age, sex, experimental design, and group as covariates (p < 0.008 for all ROIs) (Figure 3A). The strength of association was highest for the hippocampus (R2 = 0.51, p < 0.001), while the insula had the weakest association (R2 = 0.30, p > 0.0125) as reflected by the adjusted R-squared in Table 4. We observed that no other model confounders were significantly associated with the 2nd CVR metric in any ROI (p > 0.0125). In CVR Model-2, the two CVR metrics describing the first and second challenge were significantly associated with each other across all regions (p < 0.001) as expected due to the homogeneity between challenges. However, the adjusted transition rate was also significant (p < 0.008) for the thalamus, hippocampus, and centrum semiovale, but not for the insula (p = 0.38). Additionally, the CVR experimental design covariate was associated with dilation for the insula (p = 0.009), but no other covariates were significantly associated for the ROIs (p > 0.0125).
Pearson correlation estimates between transition rate, CVR-1, CVR-2 and COPE-CVR (Figure 3B) demonstrate COPE-CVR is strongly correlated to the transition rate and non-parametric CVR estimates in Groups 1 and 3, but not in Group 2. Because this was not the primary outcome of this study, further analysis was not conducted.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1199805/fneur-14-1199805-HTML/image_m/fneur-14-1199805-g003.jpg)
Figure 3. Scatterplots (A) demonstrating linear relationships between transition rate and CVR-2. The colored lines correspond to group assignment. Correlation plots (B) between transition rate, CVR-1, CVR-2, and COPE-CVR in the thalamus. For the correlation plots, the top-right values represent the R2 value between two variables, the middle-diagonal row is a histogram of values for a given variable, and the left-bottom plots represent scatterplots between two variables.
4. Discussion
In this study we investigated BOLD-CVR responses in clinical samples and found the transition rate from hyper-to-normocapnia was significantly associated with the non-parametric CVR estimate of the subsequent challenge. The association between these two responses remained robust for all but one ROI, after refining the regression model to account for the first CVR response. These results demonstrated that the Sen’s Slope non-parametric approach was conducive for analyzing the intra-session effect of a CVR experiment. Simulated CVR responses demonstrated that the use of a non-parametric CVR dilation metric is comparable to the more conventional parametric method.
Simulated CVR responses supported the use of the Sen’s Slope as a viable alternative to a parametric CVR estimate, which is appealing when estimating individual CVR responses. The parametric CVR metric tended to show better performance (i.e., lower AIC) compared to the non-parametric estimate for fast CVR response scenarios. Meanwhile, the non-parametric estimation of BOLD response had a lower AIC value for the slower simulated responses. These simulation results are consistent with various CVR studies which reported delayed responses in cerebrovascular and white matter diseases populations (3, 9, 48–50). Parametric models consider additional covariates, such as hemodynamic response functions, to account for these delayed responses (51). In contrast, the non-parametric method does not rely on assumptions concerning the underlying BOLD signal, allowing the method to be used for different brain regions, populations, and/or experimental design. Assumptions of linearity between the BOLD-CVR response and the hypercapnia stimuli may be problematic in some circumstances (27). This assumption is circumvented by taking the delta of the CO2 signal. Shifting the PETCO2 time course, for example, can address discrepancies in CO2 arrival but may not account for faster and slower regional or between-patient responses. It has been demonstrated that the underlying sigmoidal relationship between vascular resistance and PaCO2 can often result in various non-linear responses to CO2, especially in diseased populations (28). By using a non-parametric approach, we can estimate BOLD signal change without having to assume how underlying physiology may affect the signal. Collectively, these simulation results again suggest that a non-parametric, model-free approach can be utilized to characterize a single vasodilatory event.
By separating temporal features of the BOLD-CVR experiment, we were able to assess the relationship between transition periods that occur during a hypercapnic challenge. The rise in BOLD signal during the administration of CO2-rich air can be used as an indirect measure of the physiological changes that occur due to CO2 inhalation, such as an increase in cerebral blood flow. The transition rate to normocapnia, which describes the rate at which the BOLD signal returns to baseline, provides an indirect measure of the physiological response to CO2 washout. We found evidence to suggest that the rate of BOLD signal decrease during a participant’s recovery from hypercapnia air has an influence on the BOLD signal’s rise induced by the following hypercapnia challenge. For CVR Model-1, the transition to normocapnia significantly influenced BOLD signal increase during the next challenge in all ROIs after accounting for age, sex, group, and experimental paradigm. This relationship revealed that a rapid recovery was linearly associated with a greater dilatory response during a subsequent administration of CO2 stimulus. CVR Model-2 further supported our primary finding by showing that the orthogonalized transition rate remained a significant feature to explain variance in the subsequent CVR estimate. Several explanations are worthy of discussion.
One reason for this observation may be that it is physiologically optimal for these transitions to and from hypercapnia to agree to avoid unbalanced and/or exaggerated hemodynamic responses. Furthermore, recovery of the basal CBF level may contribute to a carry-over influence on the subsequent hypercapnia BOLD response. We hypothesize that the physiological underpinnings of signal recovery influence factors related to vessel dilation, such as intracellular hydrogen concentration, baseline vascular tension and resistance (27, 52). The regulation of intracellular hydrogen has been hypothesized by Duffin et al. to be a key driver of the dilatory response to hypercapnia (53). This factor could potentially be influenced by the transition period to normocapnia prior to a subsequent hypercapnic challenge. For example, a short transition rate would clear PaCO2 quickly, minimizing the disturbance of PaCO2 on hydrogen concentration. This would allow for the vessel to return to its physiological baseline state and dilation to occur uninhibited by the effects of the previous challenge. In the case of a slower transition rate, the vessel may not return to baseline before the subsequent hypercapnia challenge occurs again. This would decrease vascular tension and resistance, increasing the baseline vascular diameter prior to dilation, and result in smaller fractional BOLD signal changes in response to CO2. While this phenomenon has not been observed in the cerebral vasculature, it has been studied in other vessels, such as the brachial artery. The relationship was observed in a conventional assessment of endothelial function by cuff-occlusion, which constricts blood vessels, before applying flow-mediated dilation. Aizawa et al. (54) and Harbin et al. (55) found evidence of an inverse relationship between vessel constriction and subsequent dilation during this assessment. These findings support the notion of common mechanisms across different vascular beds.
It is interesting to note that regional differences exist for these temporal associations as well. In both dilation models, the strongest association between the two periods occurred in the hippocampus. Hippocampal arteries are small in diameter and, depending on the vascularization pattern of the hippocampus, may have limited collateral supply for arterial blood compared to more vascularized grey matter regions. Thus, the lingering effects of vessel constriction and decreased capacity for vessels to dilate may be magnified. This implies that CVR features from within low-flow brain regions may be more dependent on the outflow of CO2. This may also explain why a region such as the centrum semiovale, a white matter region with limited arterial supply, also exhibits a strong relationship between temporal features of a CVR challenge. Conversely, the insula showed the lowest association between the two periods in CVR Model-1 and showed no association between the two features in CVR Model-2. The control of vascular resistance in the insula may be a factor due to a higher density of arterioles that can quickly react to local changes in arterial gas concentrations. The rapid ability to change vessel diameter suggests that the residual effects of hypercapnia (i.e., increased CBF) are less likely to impact a subsequent hypercapnia challenge. Also, the insula is in the middle cerebral artery territory, which is in physical proximity to the chemo-sensitive receptors of the carotid bulb. Hence, the insula may contain a higher density of arterial blood gas receptors that also contribute to quick reactivity to CO2. The regional differences observed in this study warrant further investigations into CVR regional heterogeneity.
We observed the influence of experimental design on dilation in the insula for CVR Model-2. This finding indicates that certain brain regions may be more sensitive to the delivery of CO2. Tancredi and Hoge (21) demonstrated that differences in the administration of CO2 can result in significantly different changes in cerebral blood flow. Our results reveal that this effect may be dependent on tissue type or vascularization. This observation suggests that careful considerations are needed if BOLD-CVR were to be pooled across sites, groups, and/or studies, if the CVR paradigms differ.
When briefly comparing the proposed non-parametric estimates and transition rate to a traditional CVR measure in the thalamus, we observed differences in the strength of correlation amongst clinical populations. For groups with cognitive concerns, the correlation between COPE-CVR and the newly derived measures (CVR-1, CVR-2, Transition Rate) is strong. Conversely, the sleep apnea group did not exhibit a similar pattern. These results suggest that traditional CVR measurements can benefit from alternative methods of measuring the BOLD signal response to gas manipulations.
There are limitations in this study that are discussed herein. While BOLD-CVR may be a suitable choice for this analysis, it does not provide a direct physiological measurement of vessel dilation or constriction. Therefore, the hemodynamic features explored in this study are only a proxy measure of the underlying physiological changes that occur during hypercapnia. Future work could investigate these features using a more direct measurement of vascular changes. We also recognize that CO2 variability was not captured as our primary focus was situated on the changes in BOLD signal. Thus, future work might include accounting for this variability and incorporating dynamic CO2 changes or further mitigation of CO2 variability through CO2 level stabilization during scans. More advanced BOLD imaging with multi-echo and high frame rate (simultaneous multi-slice) acquisition, cerebral blood flow ASL, or hybrid BOLD/ASL imaging would provide greater confidence that the CO2 challenges in succession are intricately linked. In the case of tissue segmentation for the different ROIs, the white matter region was manually segmented while other ROIs were segmented using automated software, which could confound CVR results. Furthermore, white matter hyperintensities were not removed from the centrum semiovale region nor was the ROI corrected for potential partial voluming effects that can occur during the transformation from T1-space to BOLD-space, which could affect the CVR response for this ROI (23). Future work would include optimizing the current masks to reduce potential influences of noise while also investigating additional regions of interest. Another limitation is the choice of non-parametric CVR estimation. As noted previously, the Sen’s slope is not without its limitations. Namely, fewer assumptions in this non-parametric method compared to parametric means that a noise term is not explicitly modeled. In this instance, future work may also include a comparison between linear models such as GLM as a typical method and our non-parametric model using Sen’s Slope for BOLD-CVR signal changes. It would also be intriguing to model the transition to normocapnia using a more sophisticated model to capture the non-linear recovery response.
A strength of this study was the opportunity to characterize CVR across clinical populations and experimental paradigms. This method attempts to address the ambiguity of BOLD-CVR data interpretation in diseased populations by providing an alternative method to characterize responses to hypercapnia. We did not observe group differences between the sleep apnea and cognitive impairment groups, indicating that these two clinical populations may either have similar physiological changes occurring in the brain or that the relationship between normo- and hypercapnia CVR physiological components (such as vasodilation and vasoconstriction) are preserved across different clinical populations. The impact of disease burden could be explored in future studies by assessing inter-group differences and including an age and sex-matched healthy control group.
This study attempted to develop an alternative method to quantify the dynamic BOLD-MRI responses to hypercapnia in clinical populations. We introduced a non-parametric measure of CVR and a transition to normocapnia rate. More work is needed to refine these measures, but our results indicate that the BOLD signal induced by a physiological response to hypercapnia can be influenced by a preceding hypercapnia challenge in various regions across the brain. Future work could compare how this response differs between healthy and non-healthy populations.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by Sunnybrook Health Sciences Centre Research Ethics Board. The patients/participants provided their written informed consent to participate in this study.
Author contributions
K-JM-R contributed to study design, data analysis, interpretation, and writing. XJ, JP, RS, and AC contributed to data collection, interpretation, and writing. JR and DM contributed to study design, analysis, and edits. MT, JW, AL, SB, and BM contributed to study design, data collection, harmonization of study design and data collection, interpretation, and edits. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by funding from Fondation Leducq (16CVD05).
Acknowledgments
The authors would like to thank the MR imaging technologists, Ruby Endre and Garry Detzler, for their contributions in data acquisition. This research was supported by Foundation Leducq.
Conflict of interest
DM contributed to the development of the RespirAct Gen3 breathing circuit (Thornhill Research, Toronto, Canada).
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Footnotes
References
1. Johnson, PC. Autoregulation of blood flow. Circ Res. (1986) 59:483–95. doi: 10.1161/01.res.59.5.483
2. MacKay, CM, Skow, RJ, Tymko, MM, Boulet, LM, Davenport, MH, Steinback, CD, et al. Central respiratory chemosensitivity and cerebrovascular CO2 reactivity: a rebreathing demonstration illustrating integrative human physiology. Adv Physiol Educ. (2016) 40:79–92. doi: 10.1152/advan.00048.2015
3. Thomas, BP, Liu, P, Park, DC, van Osch, MJP, and Lu, H. Cerebrovascular reactivity in the brain white matter: magnitude, temporal characteristics, and age effects. J Cereb Blood Flow Metab. (2014) 34:242–7. doi: 10.1038/jcbfm.2013.194
4. Bhogal, AA, Philippens, MEP, Siero, JCW, Fisher, JA, Petersen, ET, Luijten, PR, et al. Examining the regional and cerebral depth-dependent BOLD cerebrovascular reactivity response at 7T. NeuroImage. (2015) 114:239–48. doi: 10.1016/j.neuroimage.2015.04.014
5. Papassin, J, Heck, O, Condamine, E, Pietras, J, Detante, O, and Krainik, A. Impaired cerebrovascular reactivity is associated with recurrent stroke in patients with severe intracranial arterial stenosis: a C02 BOLD fMRI study. J Neuroradiol. (2021) 48:339–45. doi: 10.1016/j.neurad.2020.04.005
6. Raut, RV, Nair, VA, Sattin, JA, and Prabhakaran, V. Hypercapnic evaluation of vascular reactivity in healthy aging and acute stroke via functional MRI. Neuroimage Clin. (2016) 12:173–9. doi: 10.1016/j.nicl.2016.06.016
7. Vicenzini, E, Ricciardi, MC, Altieri, M, Puccinelli, F, Bonaffini, N, di Piero, V, et al. Cerebrovascular reactivity in degenerative and vascular dementia: a transcranial Doppler study. Eur Neurol. (2007) 58:84–9. doi: 10.1159/000103642
8. Cantin, S, Villien, M, Moreaud, O, Tropres, I, Keignart, S, Chipon, E, et al. Impaired cerebral vasoreactivity to CO2 in Alzheimer’s disease using BOLD fMRI. NeuroImage. (2011) 58:579–87. doi: 10.1016/j.neuroimage.2011.06.070
9. Holmes, KR, Tang-Wai, D, Sam, K, McKetton, L, Poublanc, J, Crawley, AP, et al. Slowed temporal and parietal cerebrovascular response in patients with Alzheimer’s disease. Can J Neurol Sci. (2020) 47:366–73. doi: 10.1017/cjn.2020.30
10. Thrippleton, MJ, Shi, Y, Blair, G, Hamilton, I, Waiter, G, Schwarzbauer, C, et al. Cerebrovascular reactivity measurement in cerebral small vessel disease: rationale and reproducibility of a protocol for MRI acquisition and image processing. Int J Stroke. (2018) 13:195–206. doi: 10.1177/1747493017730740
11. Chen, D-W, Zheng, J, Shi, J, Yin, YW, Song, C, Yang, F, et al. Assessment of the cerebral hemodynamic benefits of carotid artery stenting for patients with preoperative hemodynamic impairment using cerebral single photon emission computed tomography (SPECT) and carbon dioxide inhalation. Med Sci Monit. (2018) 24:5398–404. doi: 10.12659/MSM.909401
12. Li, N, Liu, Y, Zhao, Y, Wu, X, Tong, J, and Hua, Y. Cerebrovascular reactivity in young and old patients with obstructive sleep apnea. Sleep Med. (2020) 73:125–9. doi: 10.1016/j.sleep.2020.04.029
13. Ponsaing, LB, Lindberg, U, Rostrup, E, Iversen, HK, Larsson, HBW, and Jennum, P. Impaired cerebrovascular reactivity in obstructive sleep apnea: a case-control study. Sleep Med. (2018) 43:7–13. doi: 10.1016/j.sleep.2017.10.010
14. Wu, P-H, Rodríguez-Soto, AE, Rodgers, ZB, Englund, EK, Wiemken, A, Langham, MC, et al. MRI evaluation of cerebrovascular reactivity in obstructive sleep apnea. J Cereb Blood Flow Metab. (2020) 40:1328–37. doi: 10.1177/0271678X19862182
15. Liu, P, Liu, G, Pinho, MC, Lin, Z, Thomas, BP, Rundle, M, et al. Cerebrovascular reactivity mapping using resting-state BOLD functional MRI in healthy adults and patients with moyamoya disease. Radiology. (2021) 299:419–25. doi: 10.1148/radiol.2021203568
16. Mazerolle, EL, Ma, Y, Sinclair, D, and Pike, GB. Impact of abnormal cerebrovascular reactivity on BOLD fMRI: a preliminary investigation of moyamoya disease. Clin Physiol Funct Imaging. (2018) 38:87–92. doi: 10.1111/cpf.12387
17. Goode, SD, Krishan, S, Alexakis, C, Mahajan, R, and Auer, DP. Precision of cerebrovascular reactivity assessment with use of different quantification methods for hypercapnia functional MR imaging. AJNR Am J Neuroradiol. (2009) 30:972–7. doi: 10.3174/ajnr.A1496
18. Ziyeh, S, Rick, J, Reinhard, M, Hetzel, A, Mader, I, and Speck, O. Blood oxygen level-dependent MRI of cerebral CO2 reactivity in severe carotid stenosis and occlusion. Stroke. (2005) 36:751–6. doi: 10.1161/01.STR.0000157593.03470.3d
19. Duffin, J, Sobczyk, O, Crawley, AP, Poublanc, J, Mikulis, DJ, and Fisher, JA. The dynamics of cerebrovascular reactivity shown with transfer function analysis. NeuroImage. (2015) 114:207–16. doi: 10.1016/j.neuroimage.2015.04.029
20. Fierstra, J, van Niftrik, C, Warnock, G, Wegener, S, Piccirelli, M, Pangalu, A, et al. Staging hemodynamic failure with blood oxygen-level-dependent functional magnetic resonance imaging cerebrovascular reactivity: a comparison versus gold standard (15O-)H2O-positron emission tomography. Stroke. (2018) 49:621–9. doi: 10.1161/STROKEAHA.117.020010
21. Tancredi, FB, and Hoge, RD. Comparison of cerebral vascular reactivity measures obtained using breath-holding and CO2 inhalation. J Cereb Blood Flow Metab. (2013) 33:1066–74. doi: 10.1038/jcbfm.2013.48
22. Blockley, NP, Driver, ID, Francis, ST, Fisher, JA, and Gowland, PA. An improved method for acquiring cerebrovascular reactivity maps. Magn Reson Med. (2011) 65:1278–86. doi: 10.1002/mrm.22719
23. Sam, K, Peltenburg, B, Conklin, J, Sobczyk, O, Poublanc, J, Crawley, AP, et al. Cerebrovascular reactivity and white matter integrity. Neurology. (2016) 87:2333–9. doi: 10.1212/WNL.0000000000003373
24. Poublanc, J, Crawley, AP, Sobczyk, O, Montandon, G, Sam, K, Mandell, DM, et al. Measuring cerebrovascular reactivity: the dynamic response to a step hypercapnic stimulus. J Cereb Blood Flow Metab. (2015) 35:1746–56. doi: 10.1038/jcbfm.2015.114
25. Prokopiou, PC, Pattinson, KTS, Wise, RG, and Mitsis, GD. Modeling of dynamic cerebrovascular reactivity to spontaneous and externally induced CO2 fluctuations in the human brain using BOLD-fMRI. NeuroImage. (2019) 186:533–48. doi: 10.1016/j.neuroimage.2018.10.084
26. Champagne, AA, Bhogal, AA, Coverdale, NS, Mark, CI, and Cook, DJ. A novel perspective to calibrate temporal delays in cerebrovascular reactivity using hypercapnic and hyperoxic respiratory challenges. NeuroImage. (2019) 187:154–65. doi: 10.1016/j.neuroimage.2017.11.044
27. Duffin, J, Sobczyk, O, Crawley, A, Poublanc, J, Venkatraghavan, L, Sam, K, et al. The role of vascular resistance in BOLD responses to progressive hypercapnia. Hum Brain Mapp. (2017) 38:5590–602. doi: 10.1002/hbm.23751
28. Fisher, JA, Sobczyk, O, Crawley, A, Poublanc, J, Dufort, P, Venkatraghavan, L, et al. Assessing cerebrovascular reactivity by the pattern of response to progressive hypercapnia. Hum Brain Mapp. (2017) 38:3415–27. doi: 10.1002/hbm.23598
29. Reas, ET, Hagler, DJ, White, NS, Kuperman, JM, Bartsch, H, Wierenga, CE, et al. Microstructural brain changes track cognitive decline in mild cognitive impairment. Neuroimage Clin. (2018) 20:883–91. doi: 10.1016/j.nicl.2018.09.027
30. Bettio, LEB, Rajendran, L, and Gil-Mohapel, J. The effects of aging in the hippocampus and cognitive decline. Neurosci Biobehav Rev. (2017) 79:66–86. doi: 10.1016/j.neubiorev.2017.04.030
31. BARTLETT, D, RAE, C, THOMPSON, C, BYTH, K, JOFFE, D, ENRIGHT, T, et al. Hippocampal area metabolites relate to severity and cognitive function in obstructive sleep apnea. Sleep Med. (2004) 5:593–6. doi: 10.1016/j.sleep.2004.08.004
32. Morrell, MJ, McRobbie, DW, Quest, RA, Cummin, AR, Ghiassi, R, and Corfield, DR. Changes in brain morphology associated with obstructive sleep apnea. Sleep Med. (2003) 4:451–4. doi: 10.1016/S1389-945703.00159-X
33. Oppenheimer, SM, Gelb, A, Girvin, JP, and Hachinski, VC. Cardiovascular effects of human insular cortex stimulation. Neurology. (1992) 42:1727–32. doi: 10.1212/WNL.42.9.1727
34. Prilipko, O, Huynh, N, Thomason, ME, Kushida, CA, and Guilleminault, C. An fMRI study of cerebrovascular reactivity and perfusion in obstructive sleep apnea patients before and after CPAP treatment. Sleep Med. (2014) 15:892–8. doi: 10.1016/j.sleep.2014.04.004
35. Valdés Hernández, MC, Maconick, LC, Muñoz Maniega, S, Wang, X, Wiseman, S, Armitage, PA, et al. A comparison of location of acute symptomatic vs. “silent” small vessel lesions. Int J Stroke. (2015) 10:1044–50. doi: 10.1111/ijs.12558
36. Zhang, K, Zhou, Y, Zhang, W, Li, Q, Sun, J, and Lou, M. MRI-visible perivascular spaces in basal ganglia but not centrum semiovale or hippocampus were related to deep medullary veins changes. J Cereb Blood Flow Metab. (2022) 42:136–44. doi: 10.1177/0271678X211038138
37. Keith, J, Gao, F-Q, Noor, R, Kiss, A, Balasubramaniam, G, Au, K, et al. Collagenosis of the deep medullary veins: an underrecognized pathologic correlate of white matter hyperintensities and periventricular infarction? J Neuropathol Exp Neurol. (2017) 76:299–312. doi: 10.1093/jnen/nlx0009
38. Moody, DM, Brown, WR, Challa, VR, and Anderson, RL. Periventricular venous collagenosis: association with leukoaraiosis. Radiology. (1995) 194:469–76. doi: 10.1148/radiology.194.2.7824728
39. Atwi, S, Shao, H, Crane, DE, da Costa, L, Aviv, RI, Mikulis, DJ, et al. BOLD-based cerebrovascular reactivity vascular transfer function isolates amplitude and timing responses to better characterize cerebral small vessel disease. NMR Biomed. (2019) 32:e4064. doi: 10.1002/nbm.4064
40. Hyman, BT, Phelps, CH, Beach, TG, Bigio, EH, Cairns, NJ, Carrillo, MC, et al. National Institute on Aging–Alzheimer’s association guidelines for the neuropathologic assessment of Alzheimer’s disease. Alzheimers Dement. (2012) 8:1–13. doi: 10.1016/j.jalz.2011.10.007
41. Albert, MS, DeKosky, ST, Dickson, D, Dubois, B, Feldman, HH, Fox, NC, et al. The diagnosis of mild cognitive impairment due to Alzheimer’s disease: recommendations from the National Institute on Aging-Alzheimer’s association workgroups on diagnostic guidelines for Alzheimer’s disease. Alzheimers Dement. (2011) 7:270–9. doi: 10.1016/j.jalz.2011.03.008
42. Sam, K, Crawley, AP, Conklin, J, Poublanc, J, Sobczyk, O, Mandell, DM, et al. Development of white matter hyperintensity is preceded by reduced cerebrovascular reactivity. Ann Neurol. (2016) 80:277–85. doi: 10.1002/ana.24712
43. Feinberg, DA, and Setsompop, K. Ultra-fast MRI of the human brain with simultaneous multi-slice imaging. J Magn Reson. (2013) 229:90–100. doi: 10.1016/j.jmr.2013.02.002
44. Collins, DL, Holmes, CJ, Peters, TM, and Evans, AC. Automatic 3-D model-based neuroanatomical segmentation. Hum Brain Mapp. (1995) 3:190–208. doi: 10.1002/hbm.460030304
45. Mazziotta, J, Toga, A, Evans, A, Fox, P, Lancaster, J, Zilles, K, et al. A probabilistic atlas and reference system for the human brain: international consortium for brain mapping (ICBM). Philos Trans R Soc Lond Ser B Biol Sci. (2001) 356:1293–322. doi: 10.1098/rstb.2001.0915
46. Goubran, M, Ntiri, E, Akhavein, H, Holmes, M, Nestor, S, Ramirez, J, et al. Hippocampal segmentation for atrophied brains using three-dimensional convolutional neural networks. Hum Brain Mapp. (2020) 41:291–308. doi: 10.1002/hbm.24811
47. Posse, S, Ackley, E, Mutihac, R, Rick, J, Shane, M, Murray-Krezan, C, et al. Enhancement of temporal resolution and BOLD sensitivity in real-time fMRI using multi-slab echo-volumar imaging. NeuroImage. (2012) 61:115–30. doi: 10.1016/j.neuroimage.2012.02.059
48. Leung, J, Duffin, J, Fisher, JA, and Kassner, A. MRI-based cerebrovascular reactivity using transfer function analysis reveals temporal group differences between patients with sickle cell disease and healthy controls. Neuroimage Clin. (2016) 12:624–30. doi: 10.1016/j.nicl.2016.09.009
49. McKetton, L, Sobczyk, O, Duffin, J, Poublanc, J, Sam, K, Crawley, AP, et al. The aging brain and cerebrovascular reactivity. NeuroImage. (2018) 181:132–41. doi: 10.1016/j.neuroimage.2018.07.007
50. Champagne, AA, and Bhogal, AA. Insights into cerebral tissue-specific response to respiratory challenges at 7T: evidence for combined blood flow and CO2-mediated effects. Front Physiol. (2021) 12:601369. doi: 10.3389/fphys.2021.601369
51. Yao, J, Yang, HC, Wang, JH, Liang, Z, Talavage, TM, Tamer, GG, et al. A novel method of quantifying hemodynamic delays to improve hemodynamic response, and CVR estimates in CO2 challenge fMRI. J Cereb Blood Flow Metab. (2021) 41:1886–98. doi: 10.1177/0271678X20978582
52. Halani, S, Kwinta, JB, Golestani, AM, Khatamian, YB, and Chen, JJ. Comparing cerebrovascular reactivity measured using BOLD and cerebral blood flow MRI: the effect of basal vascular tension on vasodilatory and vasoconstrictive reactivity. NeuroImage. (2015) 110:110–23. doi: 10.1016/j.neuroimage.2015.01.050
53. Duffin, J, Mikulis, DJ, and Fisher, JA. Control of cerebral blood flow by blood gases. Front Physiol. (2021) 12:640075. doi: 10.3389/fphys.2021.640075
54. Aizawa, K, Elyas, S, Adingupu, DD, Casanova, F, Gooding, KM, Strain, WD, et al. Reactivity to low-flow as a potential determinant for brachial artery flow-mediated vasodilatation. Physiol Rep. (2016) 4:e12808. doi: 10.14814/phy2.12808
Keywords: cerebrovascular reactivity, inter-regional heterogeneity, non-parametric, sleep apnea, small vessel disease, blood oxygenation level dependent functional magnetic resonance imaging
Citation: Marchena-Romero K-J, Ji X, Sommer R, Centen A, Ramirez J, Poulin JM, Mikulis D, Thrippleton M, Wardlaw J, Lim A, Black SE and MacIntosh BJ (2023) Examining temporal features of BOLD-based cerebrovascular reactivity in clinical populations. Front. Neurol. 14:1199805. doi: 10.3389/fneur.2023.1199805
Edited by:
Hongyu An, Washington University in St. Louis, United StatesReviewed by:
Chunwei Ying, Washington University in St. Louis, United StatesPeiying Liu, University of Maryland, United States
Copyright © 2023 Marchena-Romero, Ji, Sommer, Centen, Ramirez, Poulin, Mikulis, Thrippleton, Wardlaw, Lim, Black and MacIntosh. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bradley J. MacIntosh, brad.macintosh@utoronto.ca