Cerebrovascular Resistance: The Basis of Cerebrovascular Reactivity

The cerebral vascular network regulates blood flow distribution by adjusting vessel diameters, and consequently resistance to flow, in response to metabolic demands (neurovascular coupling) and changes in perfusion pressure (autoregulation). Deliberate changes in carbon dioxide (CO₂) partial pressure may be used to challenge this regulation and assess its performance since CO₂ also acts to change vessel diameter. Cerebrovascular reactivity (CVR), the ratio of cerebral blood flow (CBF) response to CO₂ stimulus is currently used as a performance metric. However, the ability of CVR to reflect the responsiveness of a particular vascular region is confounded by that region’s inclusion in the cerebral vascular network, where all regions respond to the global CO₂ stimulus. Consequently, local CBF responses reflect not only changes in the local vascular resistance but also the effect of changes in local perfusion pressure resulting from redistribution of flow within the network. As a result, the CBF responses to CO₂ take on various non-linear patterns that are not well-described by straight lines. We propose a method using a simple model to convert these CBF response patterns to the pattern of resistance responses that underlie them. The model, which has been used previously to explain the steal phenomenon, consists of two vascular branches in parallel fed by a major artery with a fixed resistance unchanging with CO₂. One branch has a reference resistance with a sigmoidal response to CO₂, representative of a voxel with a robust response. The other branch has a CBF equal to the measured CBF response to CO₂ of any voxel under examination. Using the model to calculate resistance response patterns of the examined branch showed sigmoidal patterns of resistance response, regardless of the measured CBF response patterns. The sigmoid parameters of the resistance response pattern of examined voxels may be mapped to their anatomical location. We show an example for a healthy subject and for a patient with steno-occlusive disease to illustrate. We suggest that these maps provide physiological insight into the regulation of CBF distribution.

Introduction

Cerebral blood flow (CBF) is distributed via an extensive vascular network, where flow to a particular region, is allocated according to a number of regulatory mechanisms that operate via vascular smooth muscle to change vessel diameter and hence the resistance to flow. Cerebrovascular resistance also responds to carbon dioxide (CO₂), and cerebrovascular reactivity (CVR), the change in blood flow in response to alterations in CO₂, is used as a metric to test cerebrovascular regulatory ability (Lythgoe et al., 1999; Vesely et al., 2001). However, the ability of CVR to reflect the responsiveness of a particular vascular region is confounded by that region’s inclusion in the cerebral vascular network, where all regions respond to the global CO₂ stimulus. Consequently, local CBF responses reflect not only changes in the local vascular resistance but also changes in local perfusion pressure resulting from redistribution of flow within the network. As a result, the CBF responses to CO₂ take on various non-linear patterns that are not well-described by straight lines (Fisher et al., 2017). Thus, local changes in flow, do not necessarily reflect local changes in resistance, so that the latter cannot be inferred from the former.

Our aim is to distinguish the CVR status of a region in terms of a physiological parameter, cerebrovascular resistance. We propose using a simple model to convert measured CBF responses to CO₂ to the resistance responses that underlie them. We suggest that anatomical mapping of the resistance response pattern sigmoid parameters provides physiological insight into the regulation of CBF distribution. We note that this presentation is limited to a description of the development of the model, including its limitations, a parameter sensitivity analysis, and a brief description of its application. Although we believe that this model analysis will prove useful, a demonstration of its value must await further experience.

Anatomy and Physiology Background

The blood supply to the brain features several unique anatomical and physiological characteristics that mitigate against a loss of supply (Willie et al., 2014). Anatomically, the circle of Willis acts to redistribute blood flow upon abrupt failure of a systemic supply artery (Zarrinkoob et al., 2015), and long-term compensation involves the development of collateral pial and intracerebral vessels (Moody et al., 1990; Liebeskind, 2003). A simplified schematic of the blood flow pathway to a small vascular bed, such as that within a voxel measured with magnetic resonance imaging (MRI), can be considered as pictured in Figure 1. Blood flows from the aorta via the major extra-cerebral source arteries (e.g., carotid arteries) to a distribution network that is highly redundant due to multiple parallel pathways (e.g., interconnections within the circle of Willis and pial vasculature), then to the voxel vascular bed via smooth muscle-lined vessels.

FIGURE 1

FIGURE 1. A conceptual schematic of cerebral blood flow (CBF) pathways and resistances. Resistances are represented with electrical resistance symbols with variability indicated by an arrow. Mean arterial pressure (MAP) forces flow through the high resistance of the source arteries to anastomose around the circle of Willis (CoW) where they are distributed to vascular territories, via the anterior, middle and posterior cerebral arteries (ACA, MCA, and PCA, respectively) in each hemisphere. The interconnected pial networks supply the brain parenchyma via penetrating arterioles, with flow finally reaching the tissue volumes seen in magnetic resonance imaging (MRI) scans through regulating resistances to venous pressure (Pv).

Physiologically, CBF is largely controlled via vascular smooth muscle. Multiple layers of vascular smooth muscle cells cover the large pial arteries on the surface of the cortex, which branch into penetrating arterioles sheathed in a single layer of vascular smooth muscle cells (Nishimura et al., 2007), and enter the cortical parenchyma where capillary control may occur (Hall et al., 2014; Attwell et al., 2016). Two mechanisms control vascular smooth muscle cells to regulate regional blood flow by manipulating vascular diameter and hence vascular resistance. (i) Autoregulation maintains resting blood flow despite variations in brain perfusion pressure (Tan and Taylor, 2014; Tzeng and Ainslie, 2014), and (ii) neurovascular coupling increases local blood flow in response to increased metabolic demand (Attwell et al., 2011, 2016; Phillips et al., 2016). All of these factors act in concert to change vascular resistance by altering vascular diameter according to circumstances.

The CBF Response to CO₂

Cerebrovascular resistance also responds to CO₂, and, in health, the balanced changes in regional vascular resistances in response to a progressive global increase in CO₂ results in a symmetrical, stereotypical sigmoidal pattern of progressive increase in blood flow in most regions of the brain (Battisti-Charbonney et al., 2011; Bhogal et al., 2014; Sobczyk et al., 2014), as long as arterial blood pressure remains constant (Regan et al., 2014). However, in patients with localized cerebral vascular disease resulting in single or multifocal regions of vascular stenosis, due to for example atherosclerosis or vasculitis, the changes in blood flow with a progressive hypercapnic stimulus become asymmetrical. Healthier vessels draw a disproportionate flow of blood at the expense of the more compromised vessels, and as a result the flow response patterns deviate from a sigmoidal response pattern. Indeed, in a survey using a ramp of a controlled, increasing CO₂ stimulus to explore the full range of cerebrovascular responses, we found that the blood-oxygen-level dependency (BOLD) responses in voxels measured with MRI take on any of four patterns. These may be described as increasing, decreasing, inverted U-shaped and U-shaped (Fisher et al., 2017), with only the increasing pattern (the normal response) as a sigmoidal pattern. With such non-linear responses a more sophisticated method of analysis is needed than a linear fit of BOLD vs. CO₂.

A Conceptual Model for CO₂-Induced Redistribution of CBF

To fit the different patterns of BOLD response to CO₂ we used an approach based on the conceptual model shown in Figure 2. This model takes into account the significant resistance of the major blood vessels (Rart) that conduct blood flow to the brain (Faraci and Heistad, 1990; Warnert et al., 2016), and can explain the paradoxical reduction in flow observed during a vasodilating stimulus (Sobczyk et al., 2014) first described as cerebrovascular steal (Brawley, 1968; Symon, 1969; Conklin et al., 2010) or a reverse Robin Hood phenomenon (Alexandrov et al., 2007).

FIGURE 2

FIGURE 2. The resistance model; determining the pattern of resistance change from the pattern of blood flow change during a ramp CO₂ stimulus. MAP, mean arterial blood pressure; Pv, venous blood pressure; Rart, major arteries resistance; Pbranch, perfusion pressure at branch; Rref, standard reference branch resistance; Fref, standard reference branch blood flow; Rvox, examined voxel branch resistance; Fvox, examined voxel branch blood flow; Ftotal, total blood flow through Rart.

A global CO₂ stimulus causes vasodilation in both branches, and consequently a decrease in branch pressure because of the upstream resistance of the major cerebral arteries. If an unhealthy branch cannot vasodilate, and the other healthy branch can reduce its resistance to increase its flow, a reduction in branch pressure results and consequently the flow in the unhealthy branch decreases and steal results. It should be noted that the venous pressure is assumed to be zero in this model. This simplified model therefore not only explains the steal phenomenon but also shows that any mismatch of resistance responses between the two branches will alter the distribution of flow between them. Thus, CBF is distributed to brain regions such that flow is preferentially routed to regions of lowest resistance.

Using this resistance model of the interactions between the two branches it is possible to calculate the patterns of the resistance responses to a ramp CO₂ stimulus of each voxel from their BOLD response patterns (Duffin et al., 2017). These resistance responses to CO₂ are sigmoidal, no matter what the BOLD-CO₂ response pattern, because of limitations to vasodilation and vasoconstriction at high and low end-tidal partial pressures of CO₂ (PETCO₂) respectively. The simplified model provides a way of fitting the various BOLD patterns of response to CO₂ that has a physiological basis; an interaction among the network of vascular resistances whose responses to CO₂ have a sigmoidal pattern.

A Standardized Resistance Response Pattern

This application of the model can be further developed. As previously described (Duffin et al., 2017) its use was limited to voxel BOLD response observations in individual subjects, choosing pairs of BOLD patterns of response to CO₂ for the model vascular beds, one from a voxel with a robust response serving as a standard comparator for any other voxel BOLD response pattern. This use of a robust BOLD response pattern as a standard comparator for other voxel BOLD response patterns in an individual subject led to the idea of using a resistance sigmoid response pattern as a standard comparator instead. A standard resistance sigmoid pattern can thus be used in one vascular bed of the model as the comparator for BOLD response patterns observed in any voxel.

With such a standard comparator, the BOLD response patterns in any voxel can be examined to determine the model resistance pattern of response for that voxel. That pattern is sigmoidal, and since it has been derived for each voxel with the model using the same standard comparator resistance response pattern, all voxel resistance pattern sigmoid parameters can be compared between voxels in that subject and individualized resistance parameter maps generated. An extension of the standard comparator concept is to use the same standard reference resistance response pattern for interrogating the resistance of any voxel in any subject. Anatomical maps of the voxel resistance sigmoidal response parameters allow comparisons based on anatomical location within and between subjects, analogous to the CVR analysis described by Sobczyk et al. (2015).

Methods

Experimental Protocol

Subjects

The example data used to illustrate the resistance maps was drawn from our database of healthy control subjects and patients with cerebrovascular disease that had undergone a standardized testing protocol as a result of participation in a Research Ethics Board (REB) approved study at our institution. The data used here was from 38 (25F) control subjects aged 18–76, mean (SD) = 41.6 (16.7), and 10 patients with clinical symptoms and/or known cerebrovascular disease, (4F) aged 23–72, mean (SD) = 45.4 (16.5). Each control subject was in good health, a non-smoker, not taking medication, and had no history of cardiovascular, respiratory or cerebrovascular disease, had no structural lesions on their anatomical scan, and did not have hypertension or diabetes. The patients were chosen to reflect a wide variety of cerebrovascular dysfunction. These studies conformed to the standards set by the latest revision of the Declaration of Helsinki and were approved by the REB of the University Health Network and Health Canada. All subjects in the database were competent and gave written informed consent.

Data Acquisition

During testing, subjects breathe via a face mask, connected to a sequential gas delivery breathing circuit (Somogyi et al., 2005). Targeting PETCO₂ was via sequential gas delivery (Slessarev et al., 2007; Fisher et al., 2016) implemented by a computerized gas blender (RespirAct, Thornhill Research, Inc., Toronto, ON, Canada). The CO₂ stimulus sequence is illustrated in Figure 3 and consists of clamping PETCO₂ at the subject’s resting baseline for 2 min, a step increase of 10 mmHg for 2 min followed by 2 min at baseline. PETCO₂ is then slowly reduced by hyperventilation over 30 s, followed by a steady rise in PETCO₂ at a rate of 0.1 mmHg/s to 15 mmHg or more above baseline over 4.0 min, and a return to baseline for 2 min. With this CO₂ targeting approach PETCO₂ has been shown to be equivalent to the arterial partial pressure of CO₂ (PaCO₂) (Ito et al., 2008; Willie et al., 2012).

FIGURE 3

FIGURE 3. The test protocol. The CO₂ stimulus (blue squares) and the whole brain average BOLD response scaled to 100 nL/s at resting PETCO₂ (magenta squares) in a control subject. The step change is used for measuring the speed of response in a separate analysis, and the ramp portion shown between the dashed lines is used in assessing the resistance response.

The BOLD acquisitions with echo planar imaging gradient echo (TR/TE = 2400/30 ms, 3.5 mm isotropic voxels, field of view 24 cm × 24 cm, 39 slices, slice thickness 3.5 mm, matrix size 64 × 64, number of frames = 405, flip angle = 70°) were obtained using a 3.0-Tesla Signa HDx scanner with an 8-channel phased-array receiver coil (GE Healthcare, Milwaukee, WI, United States). The acquired MRI and PETCO₂ data were analyzed using AFNI software (Cox, 1996). BOLD images were then volume registered and slice-time corrected and co-registered to an axial 3-D T1-weighted Inversion-Recovery prepared Fast Spoiled Gradient-Echo (IR-FSPGR) volume (TI/TR/TE = 450/8/3 ms, matrix size 256 × 256, field of view 22 cm × 22 cm, slice thickness = 1 mm, and flip angle = 15°) that was acquired during the same scan session (Saad et al., 2009).

Data Analysis

PETCO₂ data was re-sampled at the TR of the BOLD scan, and time-shifted to the point of coincidence with the mean brain BOLD signal where PETCO₂ abruptly decreased at the end of the ramp. The time shift is done to time align the BOLD and PETCO₂ data that were recorded on different computers. The amount of shift therefore has no physiological significance as the data acquisitions were not synchronized.

Cerebrovascular reactivity was calculated from the slope of a linear least-squares fit of the BOLD signal data series to the PETCO₂ data series over the range of PETCO₂ represented by the ramp portion of the CO₂ sequence on a voxel-by-voxel basis. The CVR was color-coded to a spectrum of colors corresponding to the direction (positive or negative) and the magnitude of the slope.

Resistance Model Description

The resistance model pictured in Figure 2 is labeled with the variables used in the analysis. Note that the model constants and variables of pressure, flow and resistance although given their respective units, they do not correspond with actual measurements in the brain. Pressure units were chosen as mmHg, and nL/s were arbitrarily chosen as the unit of model flow, making resistance units mmHg/nL/s and the appropriate axes in the graphs are labeled accordingly. The following assumptions apply:

1.   Mean arterial pressure (MAP) and Rart are constants (MAP is arbitrarily set to 100 mmHg and Pv to zero).

2.   Pbranch, Rref, Fref, Rvox, Fvox, and Ftotal are variables that change with PETCO₂.

3.   Rref changes sigmoidally with PETCO₂, with fixed parameters a, b, c and d according to the following Eq. 1:

$\text{Rref}=\text{a}+\text{b}/(1+\exp (-{\text{PETCO}}_2-\text{c})/\text{d}))\mathrm{\left(}\mathrm{1}\mathrm{\right)}$

Where:

Rref is resistance (pressure/flow, mmHg/nL/s), a function of PETCO₂

(a)   is the maximum resistance in hypocapnia (mmHg/nL/s);

(b)   is the sigmoid amplitude = minimum–maximum resistance (mmHg/nL/s);

(c)   is the sigmoid midpoint PETCO₂ where slope (sensitivity) is maximum (mmHg);

(d)   is the PETCO₂ range over which resistance is linear (mmHg).

The model uses the pattern of the BOLD response to PETCO₂ as the pattern of the Fvox response to PETCO₂. To scale between the measured BOLD and the model flow, the BOLD signal at the resting PETCO₂ (BOLDrest) was arbitrarily set to 100 nL/s, and its pattern of response to CO₂ is calculated from the change in BOLD with PETCO₂ from that at resting PETCO₂ as described in Eq. 2.

$\text{Fvox}=100*\text{BOLD}/\text{BOLDrest}\mathrm{\left(}\mathrm{2}\mathrm{\right)}$

Since the model is based on the pattern of BOLD response to CO₂, which is a relative change, this convention scales the model parameters to values that are convenient. This choice and that of MAP = 100 mmHg, determine both Rart and the Rref sigmoid parameters, as well as the calculated Rvox resistance sigmoid parameters. However, although these choices change the absolute values in the model, the pattern of Rvox resistance change in response to CO₂ remains unchanged. Thus, the values of Rart and MAP only affect the scaling of the patterns of the resistance changes with CO₂, and alternative choices could have been used (see Supplementary File).

Calculating Voxel Resistance Response Patterns

Once the BOLD response is converted to Fvox, the pattern of the Rvox response to PETCO₂, is determined by calculating Rref from Eq. 1 at each PETCO₂ of the ramp stimulus, and then calculating Rvox at each PETCO₂ with Eq. 3. The resulting pattern of the Rvox response to the ramp PETCO₂ stimulus is then fitted with the sigmoid Eq. 4 using a constrained Levenburg–Marquardt least absolute residual algorithm (LabVIEW, National Instruments, Austin, TX, United States). The Rvox sigmoid parameters derived from the model can then be color coded and mapped onto anatomical maps.

$\text{Rvox}=(\text{MAP}-\text{Fvox}*\text{Rart})/(\text{Fvox}*(1+\text{Rart}/\text{Rref}))\mathrm{\left(}\mathrm{3}\mathrm{\right)}$

$\begin{array}{c}\text{Rvox}=\text{Start}+\text{Amplitude}/(1+\exp (-{\text{PetCO}}_2\\ -\text{Midpoint})/\text{Range)})\end{array}\mathrm{\left(}\mathrm{4}\mathrm{\right)}$

Figure 4 illustrates the use of the model to convert four patterns of BOLD responses to a ramp PETCO₂ stimulus observed in a control subject (resting PETCO₂ = 35.5 mmHg) into resistance response patterns.

FIGURE 4

FIGURE 4. Examples of the four types of response patterns for model flow (left) and the calculated resistance response patterns (right) with the reference resistance (black line). The four voxel branch flow patterns (color coded solid lines) are increasing (A, red), inverted U-shaped (B, light blue), declining (C, dark blue) and U-shaped (D, orange), with the accompanying reference branch flow patterns shown as dashed lines. The dotted vertical line indicates the subject’s resting PETCO₂. Note that despite the varying shapes of the BOLD curves, the resistance curves are all sigmoidal in shape as expected except type C, which is unresponsive. The choice of BOLD scaling to model flow of 100 nL/s at the resting PETCO₂ brings the model resistances at resting PETCO₂ to approximately the same value for all patterns.

Determining the Model Constants

The constants of the model were established from a survey of the patterns of the mean BOLD response to a ramp PETCO₂ stimulus in brain regions we judged representative of a functionally healthy response (CVR between 0.5 and 0.6%/mmHg) in each of 38 healthy control subjects. To find the constants we proceeded as follows:

(1) In each subject, both branch blood flows of the model, Fref and Fvox, were set equal to a scaled version of the mean healthy BOLD pattern of response to the ramp PETCO₂ stimulus, where BOLD at resting PETCO₂ was set equal to a model flow of 100 nL/s as previously described in Eq. (2). Thus, each branch of the model had the same flow, i.e.,

$\text{Fref}=\text{Fvox}\mathrm{\left(}\mathrm{5}\mathrm{\right)}$

With this constraint the two resistances will be identical because they are connected in parallel, i.e.,

$\text{Rref}=\text{Rvox}\mathrm{\left(}\mathrm{6}\mathrm{\right)}$

(2) To determine Rart, it was assumed that a 30 mmHg pressure drop occurred from MAP to Pbranch, based on the findings detailed in Faraci and Heistad (1990). Rart was calculated from the Fref and Fvox scaled values at the resting PETCO₂ using Eq. 7 and since the flow in both branches was scaled to 100 nL/s at resting PETCO₂, Rart becomes 30/200 = 0.15 mmHg/nL/s.

$\text{Rart}=30/(\text{Fref}+\text{Fvox})\mathrm{\left(}\mathrm{7}\mathrm{\right)}$

(3) Using this value of Rart, the Rref response to the ramp PETCO₂ stimulus was determined from:

$\text{Rref}=(\text{MAP}-(\text{Fvox}+\text{Fref})*\text{Rart})/\text{Fref}\mathrm{\left(}\mathrm{8}\mathrm{\right)}$

where: Fvox and Fref were given by Eqs (2) and (5).

Finally, the sigmoid parameters a, b, c, and d associated with Rref were estimated using a constrained Levenburg–Marquardt least absolute residual algorithm (LabVIEW, National Instruments, Austin, TX, United States).

Determining the Resistance Sigmoid Fitting Constraints

To investigate CVR status in regions that may not be functionally healthy, a survey of 10 patients with various types of cerebrovascular disease was undertaken. We assumed that the model parameters for Rref remained unchanged, but Fref and Fvox were no longer the same. With Fvox given by Eq. (2), the local voxel BOLD response, and Rref given by Eq. (1), we calculated Rvox from Eq. (3). The four sigmoid parameters of Eq. (4) could then be estimated using the least absolute residual fitting algorithm for Rvox mentioned above. This analysis was used to determine the voxel resistance sigmoids for all voxels in each patient. Histograms were made for each parameter and each subject to discern the range of each parameter that would typically be encountered in a patient, and these determined the sigmoid fitting constraints.

Results

The Model Constants

Histograms of the values for the resistance sigmoid parameters a, b, c, and d found in the survey of the patterns of the mean BOLD response to a ramp PETCO₂ stimulus in brain regions we judged representative of a functionally healthy response are presented in Figure 5. The Rref sigmoid parameters, were set to the median values found in the survey of 38 control subjects. Figure 6 shows the sigmoids for all subjects in the survey and the Rref sigmoid based on the median values of the population sigmoid parameters. Rref sigmoid parameters a and b determine the reference resistance and its overall extent of change during the ramp PETCO₂ stimulus and were set to a = 0.75 and b = -0.12, respectively. The Rref sigmoid midpoint constant c reflects the PETCO₂ where sensitivity is a maximum, and was set to 40 mmHg, the median value found in the survey. The Rref range constant d reflects the range of PETCO₂ over which the Rref sigmoid can be considered linear, and thus affects the shape of the sigmoid response; it was set to 4.5 mmHg, the median value we observed in our survey. These choices establish the pattern of the reference standard resistance sigmoid response to a ramp PETCO₂ stimulus to which any other voxel resistance sigmoidal response can be compared using the model.

FIGURE 5

FIGURE 5. Histograms of the resistance sigmoid parameters observed in the survey of 38 healthy control subjects.

FIGURE 6

FIGURE 6. The resistance sigmoids of the healthy regions for all 38 healthy control subjects surveyed. The dotted line shows the reference sigmoid chosen for the model from the median values of the cohort sigmoid parameters.

The Resistance Sigmoid Fitting Constraints

The resistance sigmoid parameter histograms as shown in Figure 7 show the variation in parameter values over all voxels from the survey of 10 patients with various types of cerebrovascular disease. They indicate that the most common midpoint was approximately 40 mmHg, and the most common range was approximately 5 mmHg. Since the midpoint is the PETCO₂ at which reactivity sensitivity is maximum and the range is the span of PETCO₂ over which reactivity to PETCO₂ is linear, these two measures therefore indicate the common characteristics of the cerebrovascular response to CO₂ in this group of 10 patients. As well as setting the display scale limits, the skewedness of these parameter histograms was also reflected in the map color scales for start, amplitude and midpoint deviation, whose colors were distributed in a square law fashion in order to highlight small changes for the most common values. The parameter constraints and display scale limits listed in Table 1 were chosen based on these histograms.

FIGURE 7

FIGURE 7. Histograms of the resistance sigmoid parameters for all voxels in the 10 patients with cerebrovascular disease.

TABLE 1

TABLE 1. Constraints and scales for the resistance sigmoid parameters.

Example Resistance Parameter Maps

Figures 8, 9 compare CVR maps and resistance parameter maps for a healthy control subject and a patient with steno-occlusive disease. These examples show that maps of the voxel resistance parameters describing the resistance sigmoid response to a ramp of PETCO₂ provide further information about the physiology of blood flow distribution in the brain. Figure 10 shows the that the model analysis provides a better fit for the BOLD pattern of response to PETCO₂ than the single line of the CVR, and maps of these r² assessments, especially for the resistance sigmoid fit, provide not only assurance that the assumption of a sigmoid was correct, but also identify regions where noisy data and little or no response occurs.

FIGURE 8

FIGURE 8. Example maps for a single axial slice and their color scales from a healthy control subject with analysis graphs and resistance sigmoids from two example voxels. The maps show the locations of example voxels 1 (high amplitude, midpoint at resting PETCO₂) and 2 (low amplitude, high midpoint). The analysis graph shows the model fitting process, using the reference resistance (black solid line) and its calculated flow (black dashed line): model voxel branch flow pattern of response to PETCO₂, scaled from the % changes in BOLD vs. PETCO₂ (magenta points), is converted to resistance (blue points), then fitted with a sigmoid (blue line). The flow patterns of response for the voxel branch and reference resistance branch are then calculated from the resistance response patterns (voxel, blue line and reference black dashed line) The resistance sigmoids graphs show the relation of the fitted voxel sigmoid resistances (blue line) to the reference resistance sigmoid (black line), with their respective midpoints indicated by the vertical lines. The dashed vertical line shows the subject’s resting PETCO₂.

FIGURE 9

FIGURE 9. Example maps for a single axial slice and their color scales, with analysis graphs and resistance sigmoids from two example voxels for a patient with a near occlusion of the right carotid and bilateral foci of stenosis in the vertebral arteries. The maps show the locations of example voxels 1 (low amplitude, midpoint at resting PETCO₂) and 2 (high amplitude, high midpoint). The rest of the caption is the same as for Figure 8.

FIGURE 10

FIGURE 10. Example maps for a single axial slice of r² fit assessments for the subjects shown in Figure 8 (top row) and 9 (bottom row), with histograms of the r² fits for all voxels.

Discussion

General

First, we must emphasize that this model is not an anatomical one. The resistances are not actual, nor are the model flows; only the patterns of model flow and resistance changes with CO₂ are considered, based on the relative changes in BOLD responses to CO₂. Thus, it could be viewed as an exercise in pattern fitting, a way of fitting the variety of patterns observed in the BOLD response to CO₂, especially in patients, that has a common physiological basis. The model is simply a way of providing a standardized method based on the physiological assumptions of network interactions affecting a voxel’s perfusion pressure and a sigmoidal change in resistance with CO₂; the latter supported by the existence of physical limits of vasoconstriction and vasodilation and by the high r² fitting values obtained over diverse regions. Given that this standardized fitting procedure does account for the variety of observed BOLD-CO₂ response patterns, we suggest that the resistance sigmoid parameters of the examined voxels describe underlying physiological effects.

The choice of model flow scaling to BOLD and the choices of MAP and Rart do change the absolute values of resistance of the model branches but not the resistance responses in terms of the pattern of change with CO₂. The choice of reference resistance sigmoid parameters determines the calculated parameters for the voxels, but, since they are all related to the same reference, maps of the voxel resistance sigmoid parameters are relative maps; only the color scaling changes with different reference resistance parameters. The effect of such model parameter variations are examined in a Supplementary File.

One way of viewing this model approach to the interpretation of the BOLD-CO₂ response of a voxel is to enable a separation of the innate resistance response to a vasoactive stimulus such as CO₂ from the confounding effects of local perfusion pressure changes resulting from the network redistribution of flow during a global vasodilatory stimulus. The resulting resistance sigmoid characterizes the ability of a region (i.e., voxel or vascular bed) to respond to local regulatory demands, such as those made by neurovascular coupling. For this reason, we suggest that the resistance parameter maps provide an insight into vascular physiology and pathophysiology not discernable from conventional CVR maps, and may therefore prove useful to understanding pathophysiological changes in disease.

Indeed, the examples shown for a control subject (Figure 8) and patient (Figure 9) identify characteristics that are not discernible from the CVR maps, such as voxel 2 in the control subject that has a low amplitude as might be expected from the CVR map but a midpoint, where the responsiveness is a maximum, not at resting PETCO₂, but at an elevated PETCO₂. By contrast, voxel 2 in the patient also has an elevated midpoint in a region of high CVR and amplitude. These observations require a physiological explanation and as further experience with these maps is gained valuable insights may accrue.

Model Assumptions

The assumptions made for this model are very similar to those presented in a previous paper (Duffin et al., 2017). The model uses single variable resistances to represent the equivalent total resistances of vascular capillary beds, a conceptual simplification of the actual vascular anatomy. The major difference between the previous and present models is the adoption of a universal standard resistance response sigmoid for the reference branch here, rather than the BOLD-CO₂ response of a healthy voxel for the individual subject under study in the previous model. This change came about after accumulating experience with the previous model analysis. We noted that the reference branch resistance responses calculated from the BOLD-CO₂ responses of healthy voxels in many subjects had similar sigmoid parameters. Consequently, it became apparent that it would be possible to adopt a standard resistance response sigmoid that could be used in the model for all subjects, rather than one that was unique to each subject.

The standard resistance response used for the reference branch is intended to simulate the effect of the vascular network on the perfusion pressure experienced by the examined voxel branch of the model during a global stimulus. Previous studies of CBF have shown that redistribution during a global vasodilatory stimulus depends on the relative CVR; high CVR regions affect regions of low CVR (Sobczyk et al., 2014), resulting in cerebrovascular steal in extreme cases (Conklin et al., 2010). The reference branch resistance response does not produce the actual change in perfusion pressure experienced by the examined voxel branch but is a way of exposing every examined voxel branch to the same challenge, in all subjects.

During the ramp CO₂ challenge, it is assumed that cerebral metabolism, neural activation as well as MAP and Rart are unaffected by CO₂. Subjects are supine during these measurements and are instructed to relax but may invoke neural activation in response to their surroundings to an unknown extent. Whether hypercapnia changes cerebral metabolism is debateable (Yablonskiy, 2011), but it can affect MAP (Battisti-Charbonney et al., 2011), and if so remains a confound (Regan et al., 2014). Until MRI compatible devices capable of measuring MAP continuously, only intermittent MAP measures are available, and in our experience most subjects do not show a significant rise in MAP (>10 mmHg) during these experiments. For those subjects that do, the model results must be interpreted with considerable caution. Rart is fixed for the model but may change with very high CO₂ since all cerebral vessel smooth muscles are relaxed by CO₂ (Willie et al., 2014). However, the small vessels of vascular beds are the most sensitive (Wei et al., 1980) and dominate the resistance changes with CO₂ so that the model is a reasonable description of the pattern of vascular bed resistance changes with CO₂.

Finally, we emphasize that the success of this model analysis depends on controlling the PETCO₂ stimulus so that a slow ramp increase occurs over a wide stimulus range from hypocapnia to hypercapnia, with the hypocapnic range requiring the cooperation of subjects to hyperventilate while the RespirAct^TM controls the target PETCO₂ (i.e., prevents it from going below target) independent of ventilation. In this way the PETCO₂ stimulus should reach the limits of vasodilation and vasoconstriction that define the sigmoidal response. It is assumed that the speed of response does not influence the relations between flow distribution and CO₂ (Blockley et al., 2011; Duffin et al., 2015; Poublanc et al., 2015) because the ramp stimulus is sufficiently slow. If this assumption is incorrect then any very slow response regions would show an increased resistance sigmoid midpoint.

Limitations

The first point of discussion must consider the validity of this model analysis. It should be understood that the model analysis is dependent on the assumption that every voxel examined experiences perfusion pressure changes due to the vasodilatory response of other brain regions. Should no such perfusion pressure changes occur during the CO₂ stress, then the model premise is violated and the predicted sigmoidal resistance is incorrect. We are nevertheless confident that such competition exists (Brawley, 1968; Symon, 1969; Sobczyk et al., 2014) and accounts for the wide variety of BOLD-CO₂ response patterns observed (Fisher et al., 2017).

The model uses a fixed sigmoidal pattern of resistance change with CO₂ for the reference branch because the actual resistance change with CO₂ is unknown. The resistance sigmoid parameters for all examined voxels are therefore calculated in the model from the measured BOLD-CO₂ response using the same standard resistance response in the reference branch. In this way the calculated resistance response for the examined voxel fits the pattern of its measured BOLD response to CO₂. Consequently, the calculated sigmoid parameters for the examined voxel resistance are not measures of the actual resistance response but relative measures. Thus, we are not able to compare voxel absolute resistances, only comparisons between the sigmoid parameters describing the pattern of resistance response to a ramp of CO₂.

The model analysis depends on the assumptions listed in Section “Methods.” First, that brain metabolism and neural activation is unchanged during the ramp PETCO₂ challenge. Although, this assumption has been challenged by results from steady-state experiments (Xu et al., 2011) its validity is unknown for slow changes in PETCO₂ such as the ramp challenge we used. We suggest that the effect of brain metabolism and neural activity on the BOLD signal that occurs over the duration of the ramp stimulus is very small compared to the vascular effect of CO₂.

One of the common criticisms of studies such as this one is concerns with the fact that BOLD is only a surrogate rather than a direct measure of CBF. These concerns have been discussed at length in previous papers (Duffin et al., 2015, 2017; Fisher et al., 2017) and the interested reader can find them there. In brief, first, although BOLD deviates from a linear relationship to CBF at very high flows (Hoge et al., 1999), the increases in BOLD we observed were <10%, which is within the linear BOLD-CBF relation according to the Davis [dHb] dilution model (Davis et al., 1998). Second, BOLD measures are affected by cerebral blood volume (CBV), which also changes with hypercapnia. However, the fractional change in CBV relative to baseline is approximated by the fractional change in CBF relative to baseline raised to the power of 0.2 (Chen and Pike, 2010; Mark and Pike, 2012), and therefore likely to have little effect on the sigmoidal responses we measured. Third, BOLD may also reflect changes in cerebral metabolic rate, and changes in oxygen tension (Bulte et al., 2007; Prisman et al., 2008). With respect to metabolic rate the assumption that hypercapnia does not change cerebral metabolic rate is debateable (Yezhuvath et al., 2009; Yablonskiy, 2011), and must therefore be regarded as a caveat for these experiments. We maintained normoxic isoxia during these tests to avoid changes in oxygen tension. Finally, we used the subject’s resting PETCO₂ as baseline with the ramp from hypocapnia to hypercapnia passing through this PCO₂, because the BOLD-PETCO₂ relationship depends on the baseline PETCO₂ (Cohen et al., 2002; Sobczyk et al., 2014; Halani et al., 2015).

At present BOLD is a widely available method that provides a sufficient spatial and temporal resolution to be useful, and methods to obtain similar resolutions with direct flow measurement methods are not widely available. Moreover, where a direct comparison of BOLD and flow measures such as arterial spin labeling (ASL) are made, the authors conclude that ASL MRI confirms that, even in patients with steno-occlusive disease, the BOLD MRI signal response to hypercapnia predominantly reflects changes in CBF (Mandell et al., 2008), and this conclusion is supported by recent experiments comparing BOLD and positive emission tomography (Fierstra et al., 2018). Therefore, we suggest that, for this initial exploration of a method to derive cerebrovascular resistance response patterns, measures of BOLD response patterns to a CO₂ stimulus ramp are appropriate. When improvements in CBF measurement become available the method can be applied to them.

That the flow resistance of the supply artery (Rart) is fixed and does not change with CO₂ requires investigation. Current experiments find that hypercapnia does affect the larger supply arteries (Willie et al., 2012), but are confounded by a concomitant rise in MAP (Regan et al., 2014). This assumption has the advantage of simplifying the model, but we suggest that a more complex model incorporating changes of Rart with CO₂ may not make a significant difference to the calculated resistances because of the low CO₂ sensitivity of the major vessels; however, this assumption remains to be tested.

Conclusion

The development of an analysis method to generate the resistance sigmoid parameter maps presented here evolved from a simple model of interaction between two vascular beds, which was used as a way of separating the two effects on the BOLD-CO₂ response pattern of a vascular bed; the innate resistance response to CO₂, and the changes in perfusion pressure resulting from the redistribution of flow during a global CO₂ stimulus. Since the response of a vascular bed to a vasoactive stimulus is not usually to a global stimulus, we suggest that the innate resistance response to a vasoactive stimulus may describe the response to a local stimulus such as the metabolic demand of neurovascular coupling. In that case, the resistance sigmoid parameter maps should provide considerable insight into the underlying physiology and pathophysiology of the control of the CBF and prove to be of clinical value. This conclusion therefore merits further testing as experience with resistance sigmoid parameter maps accumulate in patients with a wide range of neurological and cerebrovascular diseases.

Author Contributions

OS, LM, and KS data collection. OS, LM, JP, and JD data analysis. JD, AC, WM, LV, DM, and JF concept and design. JD, WM, DM, and JF manuscript drafting. JD, OS, LM, AC, JP, LV, KS, WM, DM, and JF manuscript revision and approval.

Conflict of Interest Statement

JF and JD are senior scientists at Thornhill Research, Inc., (TRI), a spin-off company affiliated with the University Health Network that developed the RespirAct^TM, a non-commercial research tool assembled by TRI to enable cerebrovascular reactivity studies. DM is a shareholder in TRI. OS is a part time employee of TRI.

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.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2018.00409/full#supplementary-material

References

Alexandrov, A. V., Sharma, V. K., Lao, A. Y., Tsivgoulis, G., Malkoff, M. D., and Alexandrov, A. W. (2007). Reversed robin hood syndrome in acute ischemic stroke patients. Stroke 38, 3045–3048. doi: 10.1111/j.1552-6569.2006.00088.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Attwell, D., Buchan, A. M., Charpak, S., Lauritzen, M., Macvicar, B. A., and Newman, E. A. (2011). Glial and neuronal control of brain blood flow. Nature 468, 232–243. doi: 10.1038/nature09613

PubMed Abstract | CrossRef Full Text | Google Scholar

Attwell, D., Mishra, A., Hall, C. N., O’Farrell, F. M., and Dalkara, T. (2016). What is a pericyte? J. Cereb. Blood Flow Metab. 36, 451–455. doi: 10.1177/0271678X15610340

PubMed Abstract | CrossRef Full Text | Google Scholar

Battisti-Charbonney, A., Fisher, J., and Duffin, J. (2011). The cerebrovascular response to carbon dioxide in humans. J. Physiol. 589(Pt 12), 3039–3048. doi: 10.1113/jphysiol.2011.206052

CrossRef Full Text | Google Scholar

Bhogal, A., Siero, J. C., Fisher, J. A., Froeling, M., Luijten, P., Philippens, M., et al. (2014). Investigating the non-linearity of the BOLD cerebrovascular reactivity response to targeted hypo/hypercapnia at 7T. Neuroimage 98, 296–305. doi: 10.1016/j.neuroimage.2014.05.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Blockley, N. P., Driver, I. D., Francis, S. T., Fisher, J. A., and Gowland, P. A. (2011). An improved method for acquiring cerebrovascular reactivity maps. Magn. Reson. Med. 65, 1278–1286. doi: 10.1002/mrm.22719

PubMed Abstract | CrossRef Full Text | Google Scholar

Brawley, B. W. (1968). The pathophysiology of intracerebral steal following carbon dioxide inhalation, an experimental study. Scand. J. Clin. Lab. Invest. 22(Suppl. 102), XIII. doi: 10.3109/00365516809169045

PubMed Abstract | CrossRef Full Text | Google Scholar

Bulte, D. P., Chiarelli, P. A., Wise, R. G., and Jezzard, P. (2007). Cerebral perfusion response to hyperoxia. J. Cereb. Blood Flow Metab. 27, 69–75. doi: 10.1038/sj.jcbfm.9600319

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J. J., and Pike, G. B. (2010). MRI measurement of the BOLD-specific flow–volume relationship during hypercapnia and hypocapnia in humans. Neuroimage 53, 383–391. doi: 10.1016/j.neuroimage.2010.07.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Cohen, E. R., Ugurbil, K., and Kim, S. G. (2002). Effect of basal conditions on the magnitude and dynamics of the blood oxygenation level-dependent fMRI response. J. Cereb. Blood Flow Metab. 22, 1042–1053. doi: 10.1097/00004647-200209000-00002

PubMed Abstract | CrossRef Full Text | Google Scholar

Conklin, J., Fierstra, J., Crawley, A. P., Han, J. S., Poublanc, J., Mandell, D. M., et al. (2010). Impaired cerebrovascular reactivity with steal phenomenon is associated with increased diffusion in white matter of patients with Moyamoya disease. Stroke 41, 1610–1616. doi: 10.1161/STROKEAHA.110.579540

PubMed Abstract | CrossRef Full Text | Google Scholar

Cox, R. W. (1996). AFNI: software for analysis and visualization of functional magnetic resonance neuroimages. Comput. Biomed. Res. 29, 162–173. doi: 10.1006/cbmr.1996.0014

CrossRef Full Text | Google Scholar

Davis, T. L., Kwong, K. K., Weisskoff, R. M., and Rosen, B. R. (1998). Calibrated functional MRI: mapping the dynamics of oxidative metabolism. Proc. Natl. Acad. Sci. U.S.A. 95, 1834–1839. doi: 10.1073/pnas.95.4.1834

PubMed Abstract | CrossRef Full Text | Google Scholar

Duffin, J., Sobczyk, O., Crawley, A., Poublanc, J., Venkatraghavan, L., Sam, K., et al. (2017). The role of vascular resistance in BOLD responses to progressive hypercapnia. Hum. Brain Mapp. 38, 5590–5602. doi: 10.1002/hbm.23751

PubMed Abstract | CrossRef Full Text | Google Scholar

Duffin, J., Sobczyk, O., Crawley, A. P., Poublanc, J., Mikulis, D. J., and Fisher, J. A. (2015). The dynamics of cerebrovascular reactivity shown with transfer function analysis. Neuroimage 114, 207–216. doi: 10.1016/j.neuroimage.2015.04.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Faraci, F., and Heistad, D. (1990). Regulation of large cerebral arteries and cerebral microvascular pressure. Circ. Res. 66, 8–17. doi: 10.1161/01.RES.66.1.8

CrossRef Full Text | Google Scholar

Fierstra, J., van Niftrik, C., Warnock, G., Wegener, S., Piccirelli, M., Pangalu, A., et al. (2018). Staging hemodynamic failure with blood oxygen-level–dependent functional magnetic resonance imaging cerebrovascular reactivity. Stroke 49, 621–629. doi: 10.1161/strokeaha.117.020010

PubMed Abstract | CrossRef Full Text | Google Scholar

Fisher, J. A., Iscoe, S., and Duffin, J. (2016). Sequential gas delivery provides precise control of alveolar gas exchange. Respir. Physiol. Neurobiol. 225, 60–69. doi: 10.1016/j.resp.2016.01.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Fisher, J. A., Sobczyk, O., Crawley, A., Poublanc, J., Dufort, P., Venkatraghavan, L., et al. (2017). Assessing cerebrovascular reactivity by the pattern of response to progressive hypercapnia. Hum. Brain Mapp. doi: 10.1002/hbm.23598 [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Halani, S., Kwinta, J. B., Golestani, A. M., Khatamian, Y. B., and Chen, J. J. (2015). Comparing cerebrovascular reactivity measured using BOLD and cerebral blood flow MRI: the effect of basal vascular tension on vasodilatory and vasoconstrictive reactivity. Neuroimage 110, 110–123. doi: 10.1016/j.neuroimage.2015.01.050

PubMed Abstract | CrossRef Full Text | Google Scholar

Hall, C. N., Reynell, C., Gesslein, B., Hamilton, N. B., Mishra, A., Sutherland, B. A., et al. (2014). Capillary pericytes regulate cerebral blood flow in health and disease. Nature 508, 55–60. doi: 10.1038/nature13165

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoge, R. D., Atkinson, J., Gill, B., Crelier, G. R., Marrett, S., and Pike, G. B. (1999). Investigation of BOLD signal dependence on cerebral blood flow and oxygen consumption: the deoxyhemoglobin dilution model. Magn. Reson. Med. 42, 849–863. doi: 10.1002/(SICI)1522-2594(199911)42:5<849::AID-MRM4>3.0.CO;2-Z

PubMed Abstract | CrossRef Full Text | Google Scholar

Ito, S., Mardimae, A., Han, J., Duffin, J., Wells, G., Fedorko, L., et al. (2008). Non-invasive prospective targeting of arterial PCO₂ in subjects at rest. J. Physiol. 586, 3675–3682. doi: 10.1113/jphysiol.2008.154716

CrossRef Full Text | Google Scholar

Liebeskind, D. S. (2003). Collateral circulation. Stroke 34, 2279–2284. doi: 10.1161/01.STR.0000086465.41263.06

PubMed Abstract | CrossRef Full Text | Google Scholar

Lythgoe, D. J., Williams, S. C., Cullinane, M., and Markus, H. S. (1999). Mapping of cerebrovascular reactivity using BOLD magnetic resonance imaging. Magn. Reson. Imaging 17, 495–502. doi: 10.1016/S0730-725X(98)00211-2

CrossRef Full Text | Google Scholar

Mandell, D. M., Han, J. S., Poublanc, J., Crawley, A. P., Stainsby, J. A., Fisher, J. A., et al. (2008). Mapping cerebrovascular reactivity using blood oxygen level-dependent MRI in Patients with arterial steno-occlusive disease: comparison with arterial spin labeling MRI. Stroke 39, 2021–2028. doi: 10.1161/STROKEAHA.107.506709

PubMed Abstract | CrossRef Full Text | Google Scholar

Mark, C. I., and Pike, G. B. (2012). Indication of BOLD-specific venous flow-volume changes from precisely controlled hyperoxic vs. hypercapnic calibration. J. Cereb. Blood Flow Metab. 32, 709–719. doi: 10.1038/jcbfm.2011.174

PubMed Abstract | CrossRef Full Text | Google Scholar

Moody, D. M., Bell, M. A., and Challa, V. R. (1990). Features of the cerebral vascular pattern that predict vulnerability to perfusion or oxygenation deficiency: an anatomic study. Am. J. Neuroradiol. 11, 431–439.

Google Scholar

Nishimura, N., Schaffer, C. B., Friedman, B., Lyden, P. D., and Kleinfeld, D. (2007). Penetrating arterioles are a bottleneck in the perfusion of neocortex. Proc. Natl. Acad. Sci. U.S.A. 104, 365–370. doi: 10.1073/pnas.0609551104

PubMed Abstract | CrossRef Full Text | Google Scholar

Phillips, A. A., Chan, F. H., Zheng, M. M., Krassioukov, A. V., and Ainslie, P. N. (2016). Neurovascular coupling in humans: physiology, methodological advances and clinical implications. J. Cereb. Blood Flow Metab. 36, 647–664. doi: 10.1177/0271678X15617954

PubMed Abstract | CrossRef Full Text | Google Scholar

Poublanc, J., Crawley, A. P., Sobczyk, O., Montandon, G., Sam, K., Mandell, D. M., et al. (2015). Measuring cerebrovascular reactivity: the dynamic response to a step hypercapnic stimulus. J. Cereb. Blood Flow Metab. 35, 1–11. doi: 10.1038/jcbfm.2015.114

PubMed Abstract | CrossRef Full Text | Google Scholar

Prisman, E., Slessarev, M., Han, J., Poublanc, J., Mardimae, A., Crawley, A., et al. (2008). Comparison of the effects of independently-controlled end-tidal PCO₂ and PO₂ on blood oxygen level-dependent (BOLD) MRI. J. Magn. Reson. Imaging 27, 185–191. doi: 10.1002/jmri.21102

PubMed Abstract | CrossRef Full Text | Google Scholar

Regan, R. E., Fisher, J. A., and Duffin, J. (2014). Factors affecting the determination of cerebrovascular reactivity. Brain Behav. 4, 775–788. doi: 10.1002/brb3.275

PubMed Abstract | CrossRef Full Text | Google Scholar

Saad, Z. S., Glen, D. R., Chen, G., Beauchamp, M. S., Desai, R., and Cox, R. W. (2009). A new method for improving functional-to-structural MRI alignment using local Pearson correlation. Neuroimage 44, 839–848. doi: 10.1016/j.neuroimage.2008.09.037

PubMed Abstract | CrossRef Full Text | Google Scholar

Slessarev, M., Han, J., Mardimae, A., Prisman, E., Preiss, D., Volgyesi, G., et al. (2007). Prospective targeting and control of end-tidal CO₂ and O₂ concentrations. J. Physiol. 581, 1207–1219. doi: 10.1113/jphysiol.2007.129395

CrossRef Full Text | Google Scholar

Sobczyk, O., Battisti-Charbonney, A., Fierstra, J., Mandell, D. M., Poublanc, J., Crawley, A. P., et al. (2014). A conceptual model for CO₂-induced redistribution of cerebral blood flow with experimental confirmation using BOLD MRI. Neuroimage 92, 56–68. doi: 10.1016/j.neuroimage.2014.01.051

PubMed Abstract | CrossRef Full Text | Google Scholar

Sobczyk, O., Battisti-Charbonney, A., Poublanc, J., Crawley, A. P., Sam, K., Fierstra, J., et al. (2015). Assessing cerebrovascular reactivity abnormality by comparison to a reference atlas. J. Cereb. Blood Flow Metab. 35, 213–220. doi: 10.1038/jcbfm.2014.184

PubMed Abstract | CrossRef Full Text | Google Scholar

Somogyi, R. B., Vesely, A. E., Preiss, D., Prisman, E., Volgyesi, G., Azami, T., et al. (2005). Precise control of end-tidal carbon dioxide levels using sequential rebreathing circuits. Anaesth. Intensive Care 33, 726–732.

PubMed Abstract | Google Scholar

Symon, L. (1969). The concept of intracerebral steal. Int. Anesthesiol. Clin. 7, 597–615. doi: 10.1097/00004311-196907030-00009

CrossRef Full Text | Google Scholar

Tan, C. O., and Taylor, J. A. (2014). Integrative physiological and computational approaches to understand autonomic control of cerebral autoregulation. Exp. Physiol. 99, 3–15. doi: 10.1113/expphysiol.2013.072355

PubMed Abstract | CrossRef Full Text | Google Scholar

Tzeng, Y. C., and Ainslie, P. N. (2014). Blood pressure regulation IX: cerebral autoregulation under blood pressure challenges. Eur. J. Appl. Physiol. 114, 545–559. doi: 10.1007/s00421-013-2667-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Vesely, A., Sasano, H., Volgyesi, G., Somogyi, R., Tesler, J., Fedorko, L., et al. (2001). MRI mapping of cerebrovascular reactivity using square wave changes in end-tidal PCO₂. Magn. Reson. Med. 45, 1011–1013. doi: 10.1002/mrm.1134

PubMed Abstract | CrossRef Full Text | Google Scholar

Warnert, E. A., Hart, E. C., Hall, J. E., Murphy, K., and Wise, R. G. (2016). The major cerebral arteries proximal to the Circle of Willis contribute to cerebrovascular resistance in humans. J. Cereb. Blood Flow Metab. 36, 1384–1395. doi: 10.1177/0271678X15617952

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, E. P., Kontos, H. A., and Patterson, J. L. Jr. (1980). Dependence of pial arteriolar response to hypercapnia on vessel size. Am. J. Physiol. 238, 697–703. doi: 10.1152/ajpheart.1980.238.5.H697

PubMed Abstract | CrossRef Full Text | Google Scholar

Willie, C. K., Macleod, D. B., Shaw, A. D., Smith, K. J., Tzeng, Y. C., Eves, N. D., et al. (2012). Regional brain blood flow in man during acute changes in arterial blood gases. J. Physiol. 590, 3261–3275. doi: 10.1113/jphysiol.2012.228551

PubMed Abstract | CrossRef Full Text | Google Scholar

Willie, C. K., Tzeng, Y. C., Fisher, J. A., and Ainslie, P. N. (2014). Integrative regulation of human brain blood flow. J. Physiol. 592, 841–859. doi: 10.1113/jphysiol.2013.268953

CrossRef Full Text | Google Scholar

Xu, F., Uh, J., Brier, M. R., Hart, J. Jr., Yezhuvath, U. S., Gu, H., et al. (2011). The influence of carbon dioxide on brain activity and metabolism in conscious humans. J. Cereb. Blood Flow Metab. 31, 58–67. doi: 10.1038/jcbfm.2010.153

PubMed Abstract | CrossRef Full Text | Google Scholar

Yablonskiy, D. A. (2011). Cerebral metabolic rate in hypercapnia: controversy continues. J. Cereb. Blood Flow Metab. 31, 1502–1503. doi: 10.1038/jcbfm.2011.32

PubMed Abstract | CrossRef Full Text | Google Scholar

Yezhuvath, U. S., Lewis-Amezcua, K., Varghese, R., Xiao, G., and Lu, H. (2009). On the assessment of cerebrovascular reactivity using hypercapnia BOLD MRI. NMR Biomed. 22, 779–786. doi: 10.1002/nbm.1392

PubMed Abstract | CrossRef Full Text | Google Scholar

Zarrinkoob, L., Ambarki, K., Wahlin, A., Birgander, R., Eklund, A., and Malm, J. (2015). Blood flow distribution in cerebral arteries. J. Cereb. Blood Flow Metab. 35, 648–654. doi: 10.1038/jcbfm.2014.241

PubMed Abstract | CrossRef Full Text | Google Scholar