Abstract
BACKGROUND AND PURPOSE: With the advances of magnetic resonance technology, the CBF, oxygen extraction fraction, and cerebral metabolic rate of oxygen can be measured in MRI. Our aim was to measure the CBF, oxygen extraction fraction, and cerebral metabolic rate of oxygen use in patients with different severities of middle cerebral artery stenosis or acute stroke by using the arterial spin-labeling and susceptibility-weighted imaging techniques.
MATERIALS AND METHODS: Fifty-seven patients with MCA stenosis or acute stroke were recruited and classified into 4 groups: mild MCA stenosis (group 1), severe MCA stenosis (group 2), occluded MCA (group 3), and acute stroke (group 4). Arterial spin-labeling and SWI sequences were used to acquire CBF, oxygen extraction fraction, and cerebral metabolic rate of oxygen.
RESULTS: The oxygen extraction fraction in hemispheres with mild MCA stenosis (group 1) was remarkably higher than that in the contralateral hemisphere. In addition, hemispheres with severe MCA stenosis (group 2) had significantly lower CBF and a significantly higher oxygen extraction fraction than the contralateral hemisphere. Hemispheres with occluded MCA (group 3) or acute stroke (group 4) had a significantly lower CBF and cerebral metabolic rate of oxygen and a significantly higher oxygen extraction fraction than the contralateral hemisphere.
CONCLUSIONS: The oxygen extraction fraction gradually increased in groups 1–3. When this offset a decrease in CBF, the cerebral metabolic rate of oxygen remained at a normal level. An occluded MCA led to reduction in both the CBF and cerebral metabolic rate of oxygen. Moreover, the oxygen extraction fraction and cerebral metabolic rate of oxygen significantly increased and decreased, respectively, in the occluded MCA region during acute stroke.
ABBREVIATIONS:
- ASL
- arterial spin-labeling
- Hct
- hematocrit
- CMRO2
- cerebral metabolic rate of oxygen
- HR
- high-resolution
- OEF
- oxygen extraction fraction
- rCBF
- relative cerebral blood flow
- rCMRO2
- relative cerebral metabolic rate of oxygen
- rOEF
- relative oxygen extraction fraction
- Yv
- venous oxygen saturation
A reduction in cerebral blood flow in brain tissue is typically accompanied by a compensatory increase in the oxygen extraction fraction (OEF) to maintain normal neuronal function.1 However, the risk of stroke is greatly increased once the maximum OEF is achieved.2 The presence of increased OEF in stroke is an independent predictor of subsequent stroke in patients.3,4 Consequently, both CBF and OEF are important indicators for stroke; however, neither parameter sufficiently predicts the risk of this condition. It is possible to use both parameters and the arterial oxygen content to derive cerebral metabolic rate of oxygen (CMRO2) use, which is of critical importance in the occurrence of stroke.1
Positron-emission tomography provides the most accurate in vivo OEF and CMRO2 measurements; however, it is expensive and requires administration of radioactive isotopes. Both OEF and CBF can be determined using MR susceptometry and arterial spin-labeling (ASL).5,6 MR susceptometry uses gradient-echo phase maps to estimate oxygen saturation in segments of the jugular and gray matter veins.7,8 This technique requires blood vessels to have a substantially greater length than the diameter.7,8 Susceptibility-weighted imaging is a type of MR susceptometry that can obtain measurements of venous oxygen saturation by using the phase difference between the venous blood and surrounding tissue.9,10 ASL is a noninvasive technique that can measure CBF.11,12 The gray matter CMRO2 can be obtained from local CBF and OEF measurements.
Most previous PET studies have focused on OEF and CMRO2 in patients with an occluded carotid artery; few studies have investigated the stenotic or occluded middle cerebral artery.13,14 To the best of our knowledge, no studies have assessed CBF, OEF, and CMRO2 relative to the degree of MCA stenosis. To address this question, we measured CBF, OEF, and CMRO2 in patients with different degrees of MCA stenosis or stroke via ASL and susceptibility-weighted phase imaging.
Materials and Methods
Subjects
This prospective study was approved by the institutional review board, and written informed consent was obtained from all subjects or their guardians. Ten healthy volunteers (5 men, 5 women) with a mean age of 55.9 years (age range, 40–67 years) were enrolled in the study. All volunteers had normal findings on neurologic examination and brain MR imaging, no stenotic middle cerebral or carotid artery, and no history of neurologic disease.
Eighty-four patients with MCA stenosis or acute stroke were recruited in the study and underwent clinical evaluation, laboratory testing, MR imaging examination, and color Doppler sonography to assess the cervical vasculature. Inclusion criteria were as follows: 1) unilateral stenosis of the M1 segment of the MCA or unilateral acute MCA stroke, 2) no stenosis in the internal carotid and/or common carotid artery, 3) no extracranial-intracranial vessel bypass operation, and 4) no motion artifacts in the susceptibility-weighted or ASL images. Periodic stripe-like artifacts in the phase-code direction of the images were identified as motion artifacts. One and 5 patients were excluded due to motion artifacts and bilateral MCA stenosis, respectively. Six patients were excluded because no sonographic examination of the cervical vasculature had been performed, and 10 patients were excluded because arteriosclerosis was found in the internal carotid and/or common carotid artery by color Doppler ultrasongraphy. Five patients were excluded because 5 veins were not detected in the SWIs in a unilateral hemisphere. Subsequently, 57 patients (29 men, 28 women) with a mean age of 55.9 years (age range, 37–75 years) were included in the study. Nine patients had a history of hypertension, 6 had diabetes mellitus, and 4 had both comorbidities. Fourteen patients had an acute unilateral MCA stroke. The time between symptom onset and MR imaging examination in patients with acute stroke was between 12 and 65 hours, with a median of 32 hours. All subjects and patients underwent measurements for hemoglobin (grams per deciliter) and hematocrit (Hct) levels within 3 days of MR imaging examinations.
Patient Categorization
Classification of the degree of MCA stenosis was based on the diameter of the vessel at the point of maximal narrowing relative to the normal distal lumen diameter, which was obtained from high-resolution (HR) T2-weighted images. For data analysis, MCA stenosis was graded as mild (<50% diameter reduction) or severe (>50% diameter reduction) stenosis or occlusion.
Patients were divided into 4 groups: 3 without and 1 with stroke (Table 1). Patients without stroke were classified according to the degree of MCA narrowing: mild MCA stenosis (group 1), severe MCA stenosis (group 2), and occluded MCA (group 3). All patients had normal findings on brain MR imaging, except on the MR angiography, with an acceptable level of white matter abnormalities (ie, pencil-thin lining of periventricular hyperintensity and foci of high signal intensity in the deep white matter) for subjects older than 60 years of age. The white matter abnormalities were not high signals in the diffusion-weighted images. The lesions in these patients with acute stroke were located in the region of the MCA. The number, age, and clinical characteristics of the various subgroups are listed in Table 1.
MR Imaging Measurements
All MR imaging examinations were performed on a 3T whole-body scanner (Signa HDxt; GE Healthcare, Milwaukee, Wisconsin) with an 8-channel head coil. The MR imaging acquisition protocols included axial T2-weighted imaging, diffusion-weighted imaging, 3D time-of-flight MRA, HR-T2WI, SWI, and pulsed-continuous ASL. The total scan time for each patient was approximately 18 minutes.
T2WI parameters were as follows: TR, 9600 ms; TE, 117 ms; TI, 2400 ms; section thickness, 6 and 1.0 mm gap between sections; FOV, 24 × 24 cm2; matrix size, 288 × 256; and NEX, 2. Axial single-shot spin-echo echo-planar sequences were used to acquire DWI with the following parameters: TR/TE, 4000/70 ms; b factors, 0 and 1000 s/mm2; matrix size, 128 × 128; FOV, 24 × 24 cm2; NEX, 4; section thickness, 6 mm; and section gap. 1 mm.
3D TOF MRA was performed with the following parameters: TR/TE, 20/3.4 ms; flip angle, 20°; FOV, 22 × 22 cm2; matrix size, 384 × 192; and section thickness, 1.6 mm. The HR-T2WI sequence was perpendicular to the stenotic M1 segment of the MCA shown on 3D TOF MRA. HR-T2WI sequence parameters were as follows: TR/TE, 2800/60 ms; FOV, 12 × 12 cm2; matrix size, 512 × 256; section thickness, 2 mm; section gap, 0.5 mm; and NEX, 4. Figure 1 shows representative MRA and HR-T2WI findings.
A 3D single-shot fast spin-echo with an in-plane spiral readout acquisition was performed to generate the quantitative CBF maps for ASL. The ASL pulse parameters were the following: point, 512; arm, 8; postlabeling delay, 2025 ms; NEX, 3; FOV, 22 × 22 cm2; spatial resolution, 3.4 × 3.4 × 2 mm3; section thickness, 2.0 mm; scan locations, 44; and bandwidth, 62.5. The acquisition time was 4 minutes 26 seconds.
For axial 3D SWI, a fully flow-compensated 3D fast low-angle shot sequence was used with the following parameters to obtain the OEF: TR/TE, 55.9/12.1 ms; NEX, 1; matrix, 448 × 384; flip angle, 15°; and FOV, 22 × 22 cm2. The acquisition time was 6 minutes 50 seconds.
After acquisition, the magnitude and phase images were reconstructed. The phase images were processed with high-pass filtering to remove unwanted phase shifts (field variations) from the main magnet.9
Theory
The principle of MR susceptometry is based on the measurement of differences in susceptibility between a candidate cerebral vein and its surroundings as obtained from axial phase images.15 The cerebral veins can be approximated as infinitely long cylinders parallel to a static magnetic field. The differences in susceptibility between the vein and tissue manifest as9: where γ is the gyromagnetic ratio, TE is the echo time, B0 is the main magnetic strength, and θ is the angle between the vein and B0. Meanwhile, Δχvein-tissue is the susceptibility shift between the vein and tissue, as follows9,16: where Δχdo (deoxygenated blood) is 0.27 ppm,15 Hct is the blood hematocrit, and Yv is the venous oxygen saturation. Under normal conditions, arterial oxygen saturation is nearly 100%. Consequently, OEF can be approximated as OEF = 1 − Yv.9
CBF was calculated by using the following equation17: where the partition coefficient, λ, is 0.9; the time of saturation performed before imaging, Tsat, is 2.0 seconds; the longitudinal relaxation time (T1) of gray matter for 3T, T1gm,18 is 1.2 seconds; the postlabel delay, W, is 2.0 seconds; the longitudinal relaxation time (T1) of blood for 3T, T1B,18 is 1.4 seconds; the labeling efficiency, which is a combination of both label pseudocontinuous ASL (0.8) and background suppression (0.75), α, is 0.6 (0.8 × 0.75); the labeling time, τ, is 1.5 seconds; ΔS is the difference between the reference and labeled image signal intensities; and SO, the signal intensity of the proton density–weighted reference image.
Oxygen consumption manifests as16: where Ya is the arterial oxygen saturation approximated as 1, CBF is brain-blood perfusion, and Ca is the oxygen concentration of blood per 100 mL. Ca is determined by using the hemoglobin level and O2 carrying capacity per gram of hemoglobin. For a healthy person, this is 55.6 μmol O2.16 In this study, we obtained the CMRO2 by measuring Yv in the gray matter candidate veins by using SWI and determining the local CBF in the region supplied by the MCA by using ASL.
Image Processing
Susceptibility-weighted and ASL images were transferred to the workstation (AW4.4, GE Healthcare) for postprocessing. SWI was reconstructed as axis magnitude and phase images. High-pass filtering (32 × 32 Hanning) was performed on the phase images to remove variance in the background field.
The phase value measurements of the candidate vein (ϕV) and surrounding tissue (ϕT) were based on methods proposed by Fan et al16 and Fujima et al.19 ϕV and ϕT were measured in ROIs by using a square ROI (20 × 20 pixels) containing the gray matter candidate vein placed on the phase image (Fig 2). When the length of the vessel exceeds 6× its diameter, it can be fitted to the long-cylinder model in MR susceptometry.8 Consequently, we selected veins that had a length of at least 6× their diameter. The orientation of the vein in the gray matter was verified as parallel to B0 by tracing the vessel through several adjacent axial, sagittal, and coronal sections. The sagittal and coronal sections were reconstructed from the phase axis section. When there was a tilt angle of the vein, it was calculated as follows8: where x1 and y1 are the coordinates of the vein measured for phase in the axis section, x2 and y2 are the coordinates of the vein in the fifth axis section above or below the phase axis section, and Δz is the section thickness × 5. Subsequently, the mean and SD phase values were calculated in the ROI.
The definition of ϕV was the darkest pixel or the average of at least 4 dark pixels of similar phase value in the largest vessels (Fig 2).16,19 Estimation of ϕT was by using the average phase value up to the threshold value (mean ± 2 SDs) in the ROI.16,19 For 1 candidate vein, ϕV and ϕT were measured 3 times by an experimenter who was blinded to the vascular results and the values were averaged. The Δϕvein-tissue value was calculated as ϕV − ϕT. We randomly chose 5 candidate veins on each side of the region in the cerebral cortex supplied by the MCA for each subject, to obtain 5 average Δϕvein-tissue values. The veins in the stroke regions were not included in the measurement. Data were excluded if 5 candidate veins were not available. Data from 5 patients were excluded due to this.
The ASL data were postprocessed by using an automated reconstruction script to generate CBF maps. This included motion correction, pair-wise subtraction between control and labeled images, and the generation of mean-difference images and quantitative CBF maps.
Regional perfusion in both MCAs was measured by ROIs based on a published template for healthy subjects and patients (Fig 3).20 The ROIs of the regional CBF in each hemisphere were 2100–4000 mm2. The regional CBF of tissue supplied by the MCA was quantified in both hemispheres from the top to bottom sections (15–20 sections), and the data were averaged. The regional CMRO2 was calculated for the tissue supplied by the MCA by using Yv and the local CBF values.
All patients with acute stroke had lesions located in the region around the MCA. The stroke regions were not included in the ROIs of CBF measurement.
Statistical Analysis
All data were analyzed by using SPSS 19.0 (IBM, Armonk, New York). Paired t tests were used to compare CBF, OEF, and CMRO2 between hemispheres in each subject. The ratios of the mean contralateral/lesion hemispheric CBF, OEF, and CMRO2 values were calculated and defined as the relative CBF (rCBF), OEF (rOEF), and CMRO2 (rCMRO2) in the patient groups. Analysis of variance and Student-Newman-Keuls tests were conducted to compare rCBF, rOEF, and rCMRO2 among the groups. A P value of < .05 was statistically significant.
Results
The largest diameter that an ischemic lesion demonstrated on DWI was 1.8–5.4 cm (median value, 2.6 cm).
The diameter of the vessel at the point of maximal narrowing relative to the normal distal lumen diameter was 0.675 ± 0.105 cm, 0.371 ± 0.078 cm, 0 cm, and 0.298 ± 0.254 cm in groups 1, 2, 3, and 4, respectively. The tilt angle of all veins lay in the 10°–34° range (mean values, 25.8°± 5.1°).
For the control subjects, the mean gray matter Yv, CBF, OEF, and CMRO2 of the hemispheres supplied by the MCA was 0.629 ± 0.032 (0.579–0.674), 55.7 ± 9.4 (31.9–67.4) mL × 100 g−1 × min−1, 37.1% ± 3.2% (32.6%–42.1%), and 145.8 ± 25.1 (100.6–193.4) μmol × 100 g−1 × min−1, respectively. There were no significant differences in CBF, OEF, and CMRO2 between the hemispheres (Table 2).
In group 1, significant differences between hemispheres were observed in OEF, but not in CBF or CMRO2. The OEF in hemispheres with mild MCA stenosis was significantly higher than that in the contralateral hemispheres. Additionally, in group 2, there were significant differences in CBF and OEF, but not in CMRO2 (Table 3). Hemispheres with severe MCA stenosis had a significantly lower CBF and significantly higher OEF than the contralateral hemispheres. Significant differences between hemispheres were observed for CBF, OEF, and CMRO2 in groups 3 and 4 (Table 3). In addition, hemispheres with occluded MCAs had a significantly lower CBF and CMRO2 and significantly higher OEF than the contralateral hemispheres.
There were significant differences in rCBF, rCMRO2 and rOEF in the 4 patient groups (Table 4). rCBF and rCMRO2 values gradually decreased among groups, but rOEF values gradually increased from group 1 to 3, coincident with an increase in the degree of angiostenosis. Table 5 lists the Student-Newman-Keuls results comparing the lesion-to-normal ratios in the 4 patient groups.
Discussion
In the present study, the mean venous oxygen saturation, Yv, (0.629 ± 0.032) of healthy subjects was similar to those in earlier reports (0.64 ± 0.04,6 0.638 ± 0.03,7 0.628 ± 0.05321) that were measured in the sagittal sinus by using quantitative susceptibility mapping.
Table 6 shows that the range of mean OEF measured by PET has been reported as 35%–43%.22⇓⇓⇓⇓–27 The OEF (37.1% ± 3.2%) calculated in this study was within this range. With a catheter-derived arteriovenous approach, OEF at the right jugular bulb has been calculated as 36% ± 5%,23 which is very close to the OEF calculated in this study.
CMRO2 measured by PET has been reported as ranging from 127 ± 17.4 to 156 ± 22.3 μmol × 100 g−1 × min−1 (Table 6).22⇓⇓⇓⇓–27 The measurement location was typically global brain regions in these studies.22⇓⇓–25 The CMRO2 in gray matter was measured in 2 studies, and values of 156 ± 22.3 or 184.9 ± 20.1 μmol × 100 g−1 × min−1 were reported.25,27 The ROI in this study incorporated both gray matter and white matter supplied by the MCA. Gray matter CBF is 2.8 times that of white matter,25 that is, CBF in white matter is significantly lower than in gray matter. Therefore, CBF in the ROI is lower than in the gray matter. The CBF decreased, leading to a decrease in CMRO2 in the ROI. This means that CMRO2 found in this study (145.8 ± 25.1 μmol × 100 g−1 × min−1) is lower than gray matter CMRO2.25,27 Furthermore, quantitative susceptibility mapping–based CMRO2 measurements within left and right cortical gray matter regions associated with MCA territories have been reported as 149.3 ± 27.1 and 154.8 ± 37.7 μmol × 100 g−1 × min−1.28 This finding indicates that the measurement method and CMRO2 values in this study are consistent with quantitative susceptibility mapping–based CMRO2 measurements,28 suggesting the robustness of these methods for the measurement of OEF and CMRO2.
The original phase-based regional oxygen metabolism method used a multiecho sequence.16 Multiecho phase values fit a straight line that is a phase against TE, and a new phase for a TE is extrapolated from the straight line. The difference between the new phase and the measured phase is due to noise. Multiecho sequences remove the effect of noise on phase measurements. The SD in Yv estimates due to noise is <5% of oxygen saturation, with a signal-to-noise ratio is as low as 10,16 whereas an ordinary phase image has an average SNR in the range of 40–70. This means that thermal noise within this SNR range leads to proportionately small error values. In addition, multiecho sequences require triple the calculations of single-echo sequences, which makes it difficult to apply the former in clinical practice.
We reconstructed the magnitude and phase images by using a single-echo sequence. The mean Yv (0.629 ± 0.032) in this study is similar to the value (0.638 ± 0.03) measured in the sagittal sinus by using a dual-echo sequence in healthy subjects.7 Furthermore, the OEF and CMRO2 values of the healthy subjects in this study are in agreement with the values determined by using PET or MR imaging.22⇓⇓⇓⇓⇓–28 These results demonstrate that the present method of OEF measurement is accurate; therefore, OEF measured by a single-echo sequence is simple and useful.
An earlier study used standard approximate values for Hct and hemoglobin16,28,29; however, Hct and hemoglobin were calculated directly from blood measurements in this study.
At present, there is controversy over the value of Δχdo.16 The 2 reported values (0.27 and 0.18 ppm) lead to a difference of approximately 13% in SvO2 estimates.7 We adopted 0.27 ppm in this study on the basis of recently published research.15
The susceptibility reproduced by the long-cylinder approximation has excellent accuracy and precision when the length of a vein exceeds 6× its diameter8; therefore, a vein with a length of at least 6× its diameter can be used in the long-cylinder model in MR susceptometry. On the basis of our observations, this standard is easy to meet in the veins of the brain.
All vessels were not parallel to the B0 field, and the tilt angle was not equal to zero. Veins thread tortuously through the brain sulcus. Therefore, significant lengths that are parallel to the B0 field are unlikely.
A 32 × 32 Hanning filter was used to enhance the phase image in this study. This resulted in <1% Yv bias.16 Therefore, the effects on the quantification of the vein-tissue phase caused by the 32 × 32 Hanning filter can be safely ignored.
To the best of our knowledge, no previous studies have investigated changes in CBF, OEF, and CMRO2 based on the severity of MCA stenosis. With the help of HR-T2WI, we measured and classified the degree of MCA stenosis, which has not been published previously. Moreover, ASL and SWI are noninvasive examinations compared with PET, making these techniques attractive for a clinical setting.
As shown in Table 3, CBF significantly decreased in groups 2 and 3. In addition, OEF increased from groups 2 through 4, coincident with the severity of the MCA stenosis. This increase is a compensatory mechanism to counteract the decrease in CBF and stabilize cerebral oxygen metabolism, indicating a self-regulatory mechanism in the brain.
Tables 4 and 5 show that a small increase in rOEF was present in patients with severe MCA stenosis (group 2) compared with patients with mild MCA stenosis (group 1). Most interesting, rOEF was significantly higher in patients with occlusive MCA (group 3) compared with patients with mild MCA stenosis. This result demonstrates that collateral circulation can compensate for the shortage in CBF in patients with mild MCA stenosis, but not in those with severe MCA stenosis. Moreover, in mild or severe MCA stenosis, the CMRO2 of the hemisphere with MCA stenosis was similar to or slightly lower than the contralateral CMRO2. This difference indicates that an increase in the OEF can balance a decrease in CBF in the presence of mild or severe MCA stenosis, ultimately maintaining brain oxygen availability. Conversely, when an occlusion was present in the MCA, the CMRO2 of the affected hemisphere was significantly lower than that of the contralateral hemisphere. These findings indicate that the presence of an occlusion prevents an increase in OEF to offset the decrease in CBF, leading to a decrease in CMRO2. This decrease in CMRO2 can disrupt normal cellular function or metabolism and ultimately lead to cell death. Furthermore, a decrease in CMRO2 indicates that the OEF has reached a maximal or nearly maximal level. Consequently, the risk of stroke is higher in patients with severe stenosis or an occluded MCA than in the presence of mild stenosis.
The CBF in stroke regions was significantly lower than that in the normal brain. There can be significant measurement errors if the stroke area is incorporated in the ROI; therefore, we avoided the stroke region when measuring CBF in patients with acute stroke. We observed that the OEF and CMRO2 significantly increased and decreased, respectively, in the territory of the occluded MCA during acute stroke. This observation is consistent with the findings of previous studies by using PET.30⇓⇓–33 Tables 3⇑–5 show that the CBF, OEF, and CMRO2 values of patients with acute stroke were similar to those of patients in group 3, but not group 2. Conversely, the degree of the angiostenosis in group 4 was similar to that of group 2, but not group 3. Table 5 shows that there was a significant difference in rCBF between groups 2 and 4. The OEF may have reached its maximum value because of the markedly decreased CBF in group 4, which consequently led to the acute stroke. Most patients in group 3 had neurologic symptoms, such as headache. Research has revealed that there is a high incidence of stroke in patients with an increased OEF due to an occluded carotid artery.34 As a result, patients with an increased OEF in group 3 may have been close to having a stroke.
This study has several limitations. First, although we simplified the method as much as possible, it remains complex and difficult to master for clinician doctors. Second, the controversy regarding the value of Δχdo remains. Third, this method can only obtain measurements for the entire MCA territory, not the regional OEF in the infarction area or penumbra.
Conclusions
The study obtained 2 primary results: 1) The OEF gradually increased from group 1 through 3. When the increased OEF compensated for a decreased CBF, CMRO2 remained at a normal level; however, an occluded MCA led to a reduction in both CBF and CMRO2. 2) The OEF and CMRO2 were significantly increased and decreased, respectively, in the territory around the occluded MCA in acute stroke.
REFERENCES
- Received April 20, 2015.
- Accepted after revision August 9, 2015.
- © 2016 by American Journal of Neuroradiology