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 (CO2), and cerebrovascular reactivity (CVR), the change in blood flow in response to alterations in CO2, 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 CO2 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 CO2 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 CO2 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
Enlarge this image.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 CO2
Cerebrovascular resistance also responds to CO2, and, in health, the balanced changes in regional vascular resistances in response to a progressive global increase in CO2 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 CO2 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. CO2.
A Conceptual Model for CO2-Induced Redistribution of CBF
To fit the different patterns of BOLD response to CO2 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
Enlarge this image.FIGURE 2. The resistance model; determining the pattern of resistance change from the pattern of blood flow change during a ramp CO2 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 CO2 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 CO2 stimulus of each voxel from their BOLD response patterns (Duffin et al., 2017). These resistance responses to CO2 are sigmoidal, no matter what the BOLD-CO2 response pattern, because of limitations to vasodilation and vasoconstriction at high and low end-tidal partial pressures of CO2 (PETCO2) respectively. The simplified model provides a way of fitting the various BOLD patterns of response to CO2 that has a physiological basis; an interaction among the network of vascular resistances whose responses to CO2 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 CO2 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 PETCO2 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 CO2 stimulus sequence is illustrated in Figure 3 and consists of clamping PETCO2 at the subject’s resting baseline for 2 min, a step increase of 10 mmHg for 2 min followed by 2 min at baseline. PETCO2 is then slowly reduced by hyperventilation over 30 s, followed by a steady rise in PETCO2 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 CO2 targeting approach PETCO2 has been shown to be equivalent to the arterial partial pressure of CO2 (PaCO2) (Ito et al., 2008; Willie et al., 2012).
FIGURE 3
Enlarge this image.FIGURE 3. The test protocol. The CO2 stimulus (blue squares) and the whole brain average BOLD response scaled to 100 nL/s at resting PETCO2 (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 PETCO2 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
PETCO2 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 PETCO2 abruptly decreased at the end of the ramp. The time shift is done to time align the BOLD and PETCO2 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 PETCO2 data series over the range of PETCO2 represented by the ramp portion of the CO2 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 PETCO2.
3. Rref changes sigmoidally with PETCO2, with fixed parameters a, b, c and d according to the following Eq. 1:
Rref=a+b/(1+exp(−PETCO2−c)/d))(1)
Where:
Rref is resistance (pressure/flow, mmHg/nL/s), a function of PETCO2
(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 PETCO2 where slope (sensitivity) is maximum (mmHg);
(d) is the PETCO2 range over which resistance is linear (mmHg).
The model uses the pattern of the BOLD response to PETCO2 as the pattern of the Fvox response to PETCO2. To scale between the measured BOLD and the model flow, the BOLD signal at the resting PETCO2 (BOLDrest) was arbitrarily set to 100 nL/s, and its pattern of response to CO2 is calculated from the change in BOLD with PETCO2 from that at resting PETCO2 as described in Eq. 2.
Fvox=100*BOLD/BOLDrest(2)
Since the model is based on the pattern of BOLD response to CO2, 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 CO2 remains unchanged. Thus, the values of Rart and MAP only affect the scaling of the patterns of the resistance changes with CO2, 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 PETCO2, is determined by calculating Rref from Eq. 1 at each PETCO2 of the ramp stimulus, and then calculating Rvox at each PETCO2 with Eq. 3. The resulting pattern of the Rvox response to the ramp PETCO2 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.
Rvox=(MAP−Fvox*Rart)/(Fvox*(1+Rart/Rref))(3)
Rvox=Start+Amplitude/(1+exp(−PetCO2−Midpoint)/Range))(4)
Figure 4 illustrates the use of the model to convert four patterns of BOLD responses to a ramp PETCO2 stimulus observed in a control subject (resting PETCO2 = 35.5 mmHg) into resistance response patterns.
FIGURE 4
Enlarge this image.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 PETCO2. 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 PETCO2 brings the model resistances at resting PETCO2 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 PETCO2 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 PETCO2 stimulus, where BOLD at resting PETCO2 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.,
Fref=Fvox(5)
With this constraint the two resistances will be identical because they are connected in parallel, i.e.,
Rref=Rvox(6)
(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 PETCO2 using Eq. 7 and since the flow in both branches was scaled to 100 nL/s at resting PETCO2, Rart becomes 30/200 = 0.15 mmHg/nL/s.
Rart=30/(Fref+Fvox)(7)
(3) Using this value of Rart, the Rref response to the ramp PETCO2 stimulus was determined from:
Rref=(MAP−(Fvox+Fref)*Rart)/Fref(8)
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 PETCO2 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 PETCO2 stimulus and were set to a = 0.75 and b = -0.12, respectively. The Rref sigmoid midpoint constant c reflects the PETCO2 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 PETCO2 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 PETCO2 stimulus to which any other voxel resistance sigmoidal response can be compared using the model.
FIGURE 5
Enlarge this image.FIGURE 5. Histograms of the resistance sigmoid parameters observed in the survey of 38 healthy control subjects.
FIGURE 6
Enlarge this image.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 PETCO2 at which reactivity sensitivity is maximum and the range is the span of PETCO2 over which reactivity to PETCO2 is linear, these two measures therefore indicate the common characteristics of the cerebrovascular response to CO2 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
Enlarge this image.FIGURE 7. Histograms of the resistance sigmoid parameters for all voxels in the 10 patients with cerebrovascular disease.
TABLE 1
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 PETCO2 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 PETCO2 than the single line of the CVR, and maps of these r2 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
Enlarge this image.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 PETCO2) 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 PETCO2, scaled from the % changes in BOLD vs. PETCO2 (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 PETCO2.
FIGURE 9
Enlarge this image.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 PETCO2) and 2 (high amplitude, high midpoint). The rest of the caption is the same as for Figure 8.
FIGURE 10
Enlarge this image.FIGURE 10. Example maps for a single axial slice of r2 fit assessments for the subjects shown in Figure 8 (top row) and 9 (bottom row), with histograms of the r2 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 CO2 are considered, based on the relative changes in BOLD responses to CO2. 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 CO2, 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 CO2; the latter supported by the existence of physical limits of vasoconstriction and vasodilation and by the high r2 fitting values obtained over diverse regions. Given that this standardized fitting procedure does account for the variety of observed BOLD-CO2 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 CO2. 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-CO2 response of a voxel is to enable a separation of the innate resistance response to a vasoactive stimulus such as CO2 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 PETCO2, but at an elevated PETCO2. 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-CO2 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-CO2 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 CO2 challenge, it is assumed that cerebral metabolism, neural activation as well as MAP and Rart are unaffected by CO2. 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 CO2 since all cerebral vessel smooth muscles are relaxed by CO2 (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 CO2 so that the model is a reasonable description of the pattern of vascular bed resistance changes with CO2.
Finally, we emphasize that the success of this model analysis depends on controlling the PETCO2 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 RespirActTM controls the target PETCO2 (i.e., prevents it from going below target) independent of ventilation. In this way the PETCO2 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 CO2 (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 CO2 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-CO2 response patterns observed (Fisher et al., 2017).
The model uses a fixed sigmoidal pattern of resistance change with CO2 for the reference branch because the actual resistance change with CO2 is unknown. The resistance sigmoid parameters for all examined voxels are therefore calculated in the model from the measured BOLD-CO2 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 CO2. 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 CO2.
The model analysis depends on the assumptions listed in Section “Methods.” First, that brain metabolism and neural activation is unchanged during the ramp PETCO2 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 PETCO2 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 CO2.
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 PETCO2 as baseline with the ramp from hypocapnia to hypercapnia passing through this PCO2, because the BOLD-PETCO2 relationship depends on the baseline PETCO2 (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 CO2 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 CO2 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 CO2 may not make a significant difference to the calculated resistances because of the low CO2 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-CO2 response pattern of a vascular bed; the innate resistance response to CO2, and the changes in perfusion pressure resulting from the redistribution of flow during a global CO2 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 RespirActTM, 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://[see PDF]/articles/10.3389/fnins.2018.00409/full#supplementary-material
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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].
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
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
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
Ito, S., Mardimae, A., Han, J., Duffin, J., Wells, G., Fedorko, L., et al. (2008). Non-invasive prospective targeting of arterial PCO2 in subjects at rest. J. Physiol. 586, 3675–3682. doi: 10.1113/jphysiol.2008.154716
Liebeskind, D. S. (2003). Collateral circulation. Stroke 34, 2279–2284. doi: 10.1161/01.STR.0000086465.41263.06
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
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
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
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.
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
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
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
Prisman, E., Slessarev, M., Han, J., Poublanc, J., Mardimae, A., Crawley, A., et al. (2008). Comparison of the effects of independently-controlled end-tidal PCO2 and PO2 on blood oxygen level-dependent (BOLD) MRI. J. Magn. Reson. Imaging 27, 185–191. doi: 10.1002/jmri.21102
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
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
Slessarev, M., Han, J., Mardimae, A., Prisman, E., Preiss, D., Volgyesi, G., et al. (2007). Prospective targeting and control of end-tidal CO2 and O2 concentrations. J. Physiol. 581, 1207–1219. doi: 10.1113/jphysiol.2007.129395
Sobczyk, O., Battisti-Charbonney, A., Fierstra, J., Mandell, D. M., Poublanc, J., Crawley, A. P., et al. (2014). A conceptual model for CO2-induced redistribution of cerebral blood flow with experimental confirmation using BOLD MRI. Neuroimage 92, 56–68. doi: 10.1016/j.neuroimage.2014.01.051
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
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.
Symon, L. (1969). The concept of intracerebral steal. Int. Anesthesiol. Clin. 7, 597–615. doi: 10.1097/00004311-196907030-00009
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
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
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 PCO2. Magn. Reson. Med. 45, 1011–1013. doi: 10.1002/mrm.1134
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
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
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
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
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
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
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
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
James Duffin1,2*, Olivia Sobczyk3, Larissa McKetton4, Adrian Crawley4, Julien Poublanc4, Lashmi Venkatraghavan2, Kevin Sam4, W. Alan Mutch5, David Mikulis1,3,4 and Joseph A. Fisher1,2,3
* 1Department of Physiology, University of Toronto, Toronto, ON, Canada
* 2Department of Anaesthesia and Pain Management, University Health Network, University of Toronto, Toronto, ON, Canada
* 3Institute of Medical Science, University of Toronto, Toronto, ON, Canada
* 4Joint Department of Medical Imaging and the Functional Neuroimaging Lab, University Health Network, Toronto, ON, Canada
* 5Department of Anesthesia and Perioperative Medicine, University of Manitoba, Winnipeg, MB, Canada
© 2018. This work is licensed under http://creativecommons.org/licenses/by/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.