Abstract
BACKGROUND AND PURPOSE: Similar to digital subtraction angiography, dynamic spin labeling angiography (DSLA) provides time-resolved measurements of the influx of blood into the cerebral vascular tree. We determined whether DSLA may help in assessing the degree of stenosis and whether it provides information about intracerebral collateralization and allows us to monitor the hemodynamic effects of vascular interventions.
METHODS: We developed a segmented DSLA sequence that allowed the formation of images representing inflow delays in 41-ms increments. Thirty patients with unilateral carotid artery stenosis and 10 control subjects underwent DSLA. Arrival times of the labeled arterial blood bolus were measured in the carotid siphon (CS) and the middle cerebral artery (MCA) on both sides, and the corresponding side-to-side arrival time differences (ATDs) were calculated. ATDs before and after carotid endarterectomy or percutaneous angioplasty were studied in 10 patients.
RESULTS: The degree of stenosis was significantly correlated with ATD in the cerebral vessels. Receiver operating characteristic analysis yielded a cutoff CS ATD of 110 ms to separate stenoses <70% from those ≥70%, with a sensitivity of 90% and a specificity of 67%. In one third of patients, ATD was higher in the MCA than in the CS; this finding suggested an absence of collateralization. Most patients had reduced ATD in the MCA. The degree of ATD reduction was regarded as a quantitative measure of collateralization. Successful intervention resulted in normalized ATDs.
CONCLUSION: DSLA is a promising method that allowed us to noninvasively quantify the hemodynamic effect of extracranial carotid stenosis and the resulting intracranial collateralization.
Current CT- and MR imaging-based angiographic techniques provide high spatial resolution. Combined with Doppler ultrasonography, they allow us to precisely determine the degree of extracranial internal carotid artery (ICA) stenoses (1). The results are closely correlated with those of digital subtraction angiography DSA. However, both MR and CT angiography are limited by a relatively low temporal resolution on the order of several seconds. With this limitation, it is nearly impossible to directly assess the collateralization of an ischemic area from primarily unaffected vascular territories. However, indirect information can be obtained. Supplementary to morphologic measurement of the site, length and lumen of a stenosis, contrast-enhanced dynamic MR or CT imaging provides physiologic and hemodynamic information. Regional tissue perfusion deficits associated with occlusive cerebrovascular disease can be demonstrated by using these techniques. Impaired autoregulation was found in patients with carotid artery stenosis by analyzing cerebral blood flow maps before and after acetazolamide stimulation (2). Another study showed a significant relation between hemispheric mean transit time and the degree of collateralization, as seen on DSA (3). In the circle of Willis, phase-contrast MR angiography (MRA) can be used to depict the direction of blood flow and its changes over time without the use of contrast agents (4). Nevertheless, the lack of adequate information on collateralization is an important reason why DSA is still used in patients with extracranial stenoses of the carotid or vertebral arteries despite the invasiveness of this technique and the occurrence of postprocedural diffusion disturbances in up to 10% (5).
In the present study, we addressed whether MR spin labeling is useful for evaluating hemodynamic changes associated with stenosing cerebrovascular disease. With a number of MR techniques, spin labeling is used to visualize blood flow noninvasively. In spin labeling, the magnetization of inflowing blood is prepared by using inversion or saturation pulses, which make it distinguishable from static tissue. The altered magnetization decays within a few seconds, but this time suffices to quantify blood flow in tissue (6). Most studies have focused on flow on the capillary level to assess tissue perfusion (e.g., with functional MR imaging) or to determine local organ perfusion, such as that of the kidney or the heart. Various techniques are used to specifically suppress the blood signal intensity from the arteries (e.g., by applying flow-sensitive gradients or by saturating the distal end of the labeled blood bolus in combination with delayed image acquisition) (7). Although spin labeling was show to be suitable for visualizing vessels on projection angiograms by using pulsed (8, 9) or continuous inversion (10), it is rarely used for angiography because time-of-flight or contrast-enhanced MRAs are of higher quality and, more important, 3D at identical acquisition times.
We present a DSLA technique characterized by a particularly high temporal resolution. The basic approach is to vary the inversion delay in pulsed spin-labeling angiography, as described in the early 1990s (11, 12). To our knowledge, no group has investigated the use of this technique for analyzing changes in cerebral blood flow dynamics in patients with extracranial vascular stenoses. We optimized DSLA for clinical use by introducing segmented data acquisition to shorten the acquisition time at the expense of temporal resolution. The selection of a combination of acquired lines, bandwidth, and segmentation allowed us to adjust the temporal and spatial resolution, as well as the acquisition time, over wide ranges.
The purpose of the present study was to determine whether this technique is suitable for determining the degree of stenosis in patients with extracranial carotid artery stenoses by precise analysis of arrival time differences (ATDs). In addition and probably more important, the investigation addressed the issue of whether DSLA provides information about intracerebral collateralization and whether it allows us to monitor the hemodynamic effects of interventional vascular procedures.
Methods
DSLA Technique
DSLA is based on the same principle as DSA. In DSA, contrast medium is injected directly into an artery through a catheter, and its distribution is visualized on a series of radiographs acquired in rapid succession. For analysis, the nonenhanced fluoroscopic image is digitally subtracted from the subsequent images to eliminate the background signal intensity of static tissue. In spin labeling angiography, the contrast medium is replaced by magnetic labeling of the blood, which is done at a defined point in time (pulsed spin labeling) (13, 14). As with DSA, control images without labeling are obtained and subtracted, leaving only the blood signal intensity. A quasi-dynamic study can be done by varying the delay (i.e., TI) between labeling and image acquisition (12).
We used the flow alternating inversion-recovery (FAIR) labeling scheme (15, 16). With FAIR, label and control images were acquired after global and section-selective inversion, respectively. This scheme had the advantage of labeling all blood in the transmitter coil and precluding subtraction errors due to magnetization transfer. We implemented hyperbolic secant inversion pulses that provided a good section profile and made the inversion almost insensitive to B1 inhomogeneities. Figure 1 illustrates the technique of DSLA. Instead of acquiring data in a single temporal window with a varying postlabeling delay, multiple time steps were sampled in cine-like fashion by using gradient-echo data acquisition with a small excitation flip angle. This principle of simultaneously sampling multiple inversion times was first introduced by Look and Locker for fast T1 determination (17). It is clear from Figure 1 that the readout length and the number of segments N determine the temporal resolution. After the acquisition of the N identical k-space lines for all phases after global inversion, time was given to allow for decay and complete washout of labeled blood from the section. Subsequently, the same readout train acquired N k-space lines of the control images following section selective inversion. With M lines to measure for a complete image, acquisition of a data set was finished after 2·M/N inversion pulses.
DSLA sequence scheme. After the section selective/nonselective labeling pulse, multiphase data acquisition started in a cine-like fashion. Any segmentable gradient-echo readout may be used. Figure depicts FLASH readout with threefold segmentation and flow compensation in the section and read directions. ADC indicates analog-to-digital converter; RF, radio frequency; S, section direction gradient.
However, DSLA and DSA have several differences. First and most important, magnetic labeling is subject to relaxation and can only be traced for a short period of 1 to 2 seconds. Second, the excitation pulses used for image acquisition accelerate relaxation and thereby destroy the labeling effect. Imaging with a flip angle of 90° completely excites longitudinal magnetization in the imaging plane and fully exploits the labeling effect, but does not allow for repetitive acquisition in multiple phases. By using a smaller excitation angle, data can be acquired in multiple phases after a single blood inversion. However, with a decreased measuring signal intensity compared with a single phase measurement. Nevertheless, the overall signal-to-noise ratio of such an image series is higher than that of several images acquired individually (18). Third, to achieve high temporal resolution with DSLA, only some k-space lines of each image in the series are measured after one magnetization preparation. The dynamic study thus comprises data acquired in a large number of repeated labeling experiments. This procedure requires blood flow dynamics to be the same after each magnetization, which is why ECG triggering should be used.
Patients and Subjects
The study was approved by the institutional review board, and informed consent was obtained from all participants. As a control group, 10 healthy subjects (three men, seven women; age range, 53–80 years; mean age, 67 years) were examined. The absence of stenotic lesions of extracranial and intracranial brain arteries was confirmed with extracranial and transcranial ultrasonography. Thirty patients (16 men, 14 women; age range 50–80 years; mean age, 67 years)) with unilateral extracranial stenoses of the ICA were examined. The degree of stenosis was determined (C.W., C.Z.) on DSA in 18 patients by using the North American Symptomatic Carotid Endarterectomy Trial (NASCET) method (19). In the other 12 patients, two independent examiners (J.M.V., S.J.S.) determined the degree with extracranial duplex sonography, according to standard European Carotid Surgery Trial (ECST) criteria (20). These results were converted into NASCET values (21). The average interval between assessment of the degree of stenosis and MR imaging was 4 days, with a maximum of 22 days.
ICA stenosis was symptomatic in 22 patients and asymptomatic in eight. Among the patients with symptomatic stenosis, 17 had had a cerebral infarction, and five had presented with a TIA within 4 weeks preceding the examination. In two symptomatic patients, cerebral infarction was older than 28 days. Ten of the patients underwent repeat imaging after either carotid endarterectomy (n = 6) or angioplasty with stent placement (n = 4). The maximum interval between carotid endarterectomy or stent angioplasty and imaging was 11 days; the average was 8 days. Patients with additional extracranial or intracranial vascular stenoses on the contralateral side were not included in the study.
Study Protocol
All examinations were performed by using a 1.5-T whole-body MR imager (Magnetom Vision, Siemens Erlangen, Germany) with the transmit/receive head coil. Three-lead ECG (Bruker, Rheinstetten, Germany) with active amplification was used for triggering. Proton density- and T2-weighted axial images of the entire brain were acquired by using a double-echo turbo spin-echo sequence (TR/TE1/TE2 = 3800/22/90, turbo factor of 5, flip angle of 90°), followed by time-of-flight MRA of the intracranial vessels by using a gradient-echo sequence (TR/TE = 22/5.3, flip angle of 20°).
Intracerebral flow was studied by using a matrix of 192 × 256 (6/8 230-mm field of view), resulting in 128 inversion pulses. With an interval of 2500 ms between two inversion pulses, the minimum acquisition time was 5 minutes 20 seconds. A triply segmented flow-compensated gradient-echo sequence (FLASH) was used with a TE of 6.1 ms, a bandwidth of 195 Hz per pixel, and a temporal resolution of 41 ms or 25 images per second. This was adequate for visualizing the rapid blood flow in the cerebral arteries. Use of a flip angle of 10° allowed us to track the labeled blood on a series of 40 images, covering an interval of 0 to 1640 ms.
Two ECG-triggered DSLA acquisitions with different section thicknesses (40 and 30 mm) were performed, one at the level of the circle of Willis and one at the level of the carotid siphon (CS). Both acquisitions were performed in the transverse orientation depicting the vessels in head-to-foot projection. By using FAIR, blood entering the imaging section from above (e.g., in the sagittal and transverse sinuses) was likewise visualized, but this was not a problem in a transverse projection. It was important for labeling to be performed distal to the stenosis, as the label was not imaged while passing through it. Consequently, possible passage delays did not contribute to the arrival time, and instead slower flow or late filling through the circle of Willis was solely responsible.
Analysis
Movie sequences showing cerebral blood arrival were generated from each examination for a qualitative analysis of cerebral blood flow and visual assessment of collateralization, which was rated in the following way: High collateral flow was assigned to patients with bright connecting vessels and a distinct directionality of blood flow to the affected hemisphere. Low collateral flow meant less prominent collaterals. Assessment was rated inconclusive in case of low vessel signal intensity or if the temporal resolution was insufficient to determine the flow direction.
For quantitative analysis, arrival of magnetized blood was measured in control subjects and in patients on both sides at two sites: proximal to the origin of the ophthalmic artery in the DSLA obtained at the level of the CS and at the bifurcation of the MCA in the DSLA obtained at the level of the circle of Willis. Figure 2 depicts such signal intensity time courses measured in the MCA of a patient. They were normalized to their respective maximum for comparison. For precise determination of arterial arrival, the time at which the signal intensity reached 50% of the maximal amplitude was measured on both sides. The temporal resolution was virtually increased by linear interpolation between points, thus signal intensities were used to a certain extent to improve ATD accuracy. Arrival times were determined on a 10-ms time scale instead of the 41 ms originating from the acquisition. In the group of patients, absolute arrival times and ipsilateral and contralateral ATDs were determined. In the control group, right-side measurements were subtracted from left-side measurements, yielding the ATD.
Determination of ATDs. Left, Original signal intensity-time courses measured in the left and right MCAs. Right, Normalized curves (divided by the respective maximum value). By using linear interpolation between the measured points, the intersection at which the signals reached 50% was determined. Difference was referred to as the ATD (arrows).
Statistical analysis was performed by using SPSS (SPSS Inc., Chicago, IL). The degree of stenosis and the ATDs were tested for statistically significant correlations by using a nonparametric test according to Kendall. Groups were defined (controls, stenoses <70% and ≥70%) and compared for statistically significant ATD differences by using a Mann-Whitney nonparametric test. Significance was assumed at P < .05. A receiver operating characteristic (ROC) analysis was performed to assess sensitivity and specificity of group discrimination on the basis of calculated ATDs. The optimal cutoff value for maximizing the product of sensitivity and specificity and that for 100% specificity were determined.
Results
In contrast to the analysis of ATDs, evaluation of absolute arrival times yielded no meaningful results. Tables 1 and 2 summarize the data on control subjects and patients, respectively. Listed are the degree of stenosis determined on DSA or duplex ultrasonography and ATDs of arterial blood in the CS and the MCA. In two patients with 95% stenosis, no blood signal intensity could be detected in the CS.
Quantitative results in the control group
Quantitative results in patients
Control group
The average ATD between the right and left side in the control group was 16 ms at the level of the CS and 19 ms at the level of the MCA. The greatest differences between the right and left side were 30 and 40 ms, respectively.
Patients
Qualitative analysis of DSLA.— DSLA, in particular, allowed for good qualitative assessment of collateralization via the circle of Willis, i.e., whether collateral flow occurs mainly through the anterior communicating artery from the contralateral hemisphere or rather from the ipsilateral posterior territory via the posterior communicating artery. Figure 3 shows images from a patient with pronounced collateral flow via the circle of Willis. In one patient, retrograde flow in the ophthalmic artery and leptomeningeal collaterals were clearly depicted in the DSLA image series.
Patient 30, with an 85% stenosis of the left ICA. Dynamic angiograms of the circle of Willis in foot-to-head projections at 60, 100, 140, 220, 300, and 580 ms after labeling in a-f. In a, the right ICA (upper arrow) and basilar artery (lower arrow) fill first. In b, Collateral flow into the left MCA via the left posterior communicating artery and the anterior communicating artery is shown. In d, Left ICA fills. In f, because of the finite length of the labeled bolus, all vessels but the left ICA and MCA contain unlabeled blood at this late phase.
Quantitative Analysis of ATDs.— The degree of stenosis was clearly related to the ATDs calculated for the CS and for the MCA (Fig 4). The correlation is statistically significant for both sites of measurement (P < .01). For further evaluation, the subjects were subdivided into three groups, according to their degree of stenosis: 0% (controls), <70% and ≥70%. ATDs in the CS and MCA were significantly different between the control group and patients with stenoses ≥70%, as well as between patients with ≥70% and those with <70% stenosis. The only factor that significantly differed between control subjects and patients with stenoses <70% was ATD in the CS. By using ROC analysis, cutoff values were determined to predict the classification of each patient into either of the groups with highest specificity and sensitivity. Because the ATD values at the CS level were of highest significance regarding the degree of stenosis, only these were used in the ROC analysis. Figure 5 depicts the ROC curves. The optimal cutoff value separating controls and stenoses <70% was 35 ms (sensitivity, 80%; specificity, 100%). The cutoff separating stenoses <70% from those ≥70% was 110 ms (sensitivity, 67%; specificity, 90%). Transit time delays of more than 150 ms in the CS are indicative for high-grade stenoses (sensitivity, 56%; specificity, 100%).
ATDs calculated for the CS (left) and MCA (right) increase significantly with the degree of stenosis (P < .01). Values at 60% and 75% stenosis are those of individual patients.
ROC curves for ATDs in the CS. Left, Distinction of control subjects and patients with 30–60% stenosis. Right, Distinction of control subjects and patients with <70% stenoses from those with 70–95% stenosis. ATD of >150 ms proved a high-grade stenosis (specificity, 100%; sensitivity, 56%).
Figure 6 shows a scatterplot of the individual ATDs proximal and distal to the circle of Willis of all patients, grouping them according to the visually rated poststenotic collateral flow. The degree of collateralization visible in the circle of Willis was correlated with the ATD CS/ATD MCA ratio, with a higher degree of collateralization associated with reduced distal ATDs. No interrelation could be found between the patients’ clinical picture (asymptomatic, TIA, or infarction) and the degree of collateralization.
Scatterplot of the ATDs calculated for the MCA versus those for the CS. Patients were subdivided by the degree of collateral flow visible on DSLA. Degree of collateralization was determined by subjectively assessing flow in the anterior and posterior communicating arteries. Patients with pronounced collateral flow are in the group in whom MCA ATD was less than CS ATD.
DSLA Before and After Treatment of Carotid Stenosis
Ten patients with high-grade extracranial carotid stenoses underwent DSLA before and after treatment (carotid endarterectomy in six, stent angioplasty in fur). Figure 7 presents the image series and signal intensity-time courses in the CS a patient before and after endarterectomy. One patient (patient 32) had occlusion after the intervention. All other patients had dramatically reduced delays after treatment, but three had a CS ATD of more than 60 ms. Postinterventional ATDs were still significantly elevated compared with those of controls.
In a, DSLA difference images of the circle of Willis in patient 40, with a 95% stenosis before and after endarterectomy. In b, Signal intensity-time courses in the CS. After intervention, arrival time on the affected side was the same as on that the unaffected side (about 125 ms). Amplitudes were not calibrated; hence, higher signal intensity was not associated with higher blood flow or blood volume in a vessel.
Discussion
The phenomenon of delayed arterial arrival associated with vascular stenoses was first described in 1960s, initially in studies by using conventional angiography (22) or scintigraphy (23, 24), whereas most later studies used contrast-enhanced dynamic CT (25), echo-enhanced transcranial Doppler sonography (26), or perfusion MR imaging (3, 27), to determine tissue perfusion in the presence of stenosis. All these studies did, however, require the administration of a contrast medium. This method might alter physiologic blood flow conditions and cannot be repeated arbitrarily. DSLA avoids these problems because it does not need contrast agents. Conventional CT angiography or MRA provide little information about the dynamics of poststenotic blood flow. The temporal resolution achieved with two-dimensional projection techniques is on the order of 0.5–1 second and 2–6 seconds per volume dataset when 3D acquisition with temporal and spatial interpolation is used (28, 29). Both techniques incorporate a compromise between temporal and spatial resolution, but the image frame rate still does not suffice for precise depiction of the hemodynamic situation in the circle of Willis or in the deeper CS. DSLA requires markedly longer acquisition times, but it provides 10–20-fold higher temporal resolution. As opposed to DSA, the method allows for arbitrary projection directions. DSLA is not limited to being performed as a projection technique but also allows for acquiring 3D datasets. However, the increased acquisition time is no longer acceptable in clinical examinations; this limitation suggests that temporal or spatial resolution must be reduced here as well.
This study showed that the high temporal resolution of the spin labeling method makes the technique especially suitable for monitoring hemodynamic effects of stenosing cerebrovascular disease. To our knowledge, the DSLA technique presented herein is the first fully noninvasive MR imaging procedure that allows quantification of ATDs in patients with extracranial carotid stenosis. Moreover, DSLA enabled us to show the compensation of such ATDs by adequate distal collateral flow and thus provided important supplementary information on hemodynamic changes secondary to vascular stenoses.
Our results demonstrated that the degree of stenosis was best correlated with ATD in the CS. Although the proportion of stenoses <70% was small, the results suggested that an ATD exceeding 110 ms was indicative of a stenosis of at least 70%, with a sensitivity of 67% and a specificity of 90%, or even higher specificity if control subjects and subjects with low-grade stenoses are merged in one group. All stenoses with ATD of at least 150 ms were of high grade. This additional possibility for grading vascular stenoses by using MR imaging is clinically relevant. On the basis of MRA alone, even contrast-enhanced MRA, the assessment of the degree of high-grade carotid stenoses is sometimes difficult and may lead to a misclassification of patients (30–32). Combining MRA with DSLA might significantly improve the diagnostic accuracy in such patients and furthermore reduce the necessity for conventional invasive angiography. Moreover, the technique can be used to assess the effect of treatment in patients after surgery or stent placement for ICA stenosis and to identify recurrent stenosis.
In two patients with 95% stenosis, no blood signal intensity was detected in the siphon. In both, the amount of blood was so small or it arrived so late that its signal intensity could not be distinguished from background noise; this finding suggested near-occlusion of the ICA. Large ATD variations were observed within the groups (e.g., in patient 28 who had a 30% stenosis and ATDs of more than 100 ms in the MCA and CS). This patient had no anatomic abnormalities (e.g., vessel loop) at time-of-flight MRA. However, he underwent duplex ultrasonography more than 2 weeks before DSLA. Further narrowing of the vascular lumen in the meantime to a high-degree stenosis is improbable, unless a nonsymptomatic thromboembolic event happened.
In the patients who underwent treatment, we still observed a significant transit time delay after the intervention. This result was consistent with others’ observations of a considerable variability in flow velocity and its changes in the cerebral arteries of treated patients (33–35).
In addition to information on the degree of stenosis, DLSA provides information on hemodynamic changes in the brain. In 33% of patients, ATDs in the MCA are equal to or greater than those in the CS, suggesting an absence of collateralization. However, most patients had a reduced ATD distal to the circle of Willis, most probably due to collateral flow. Whether these patients have a decreased risk of cerebral ischemia remains to be determined in further studies. Of interest, patients with stenoses ≥95% had relatively early MCA filling and thus small ATDs comparable to those of patients with intermediate-grade ICA stenosis. This observation was attributed to a well-developed collateralization and may explain why patients with near-occlusion hardly benefit from intervention (36). Although thromboembolic events secondary to atherosclerotic wall processes play a role in most patients with stroke, a hemodynamic infarct pattern can be detected in up to 50% of patients with high-grade ICA stenosis (37). For this subset of patients, it seems to be important to have methods that allow us to quantify the hemodynamic effect of ICA stenosis. In acute ischemia, this could be of great diagnostic value. DSLA can potentially help in visualizing the complex postischemic vascular reaction.
A significant limitation of our method was the inability to assess the degree of bilateral carotid stenoses, which occurred in more than 30% of patients. Determination of absolute blood arrival times in these patients is equivalent to measuring flow velocity in the vessels. However, the analysis of these absolute values gave no statistically significant results in our study. This shows that a differential method is more sensitive, as it is less dependent on section positioning. Nevertheless DSLA depicts collateralization in patients with bilateral disease as well. Having said that, in practice, patients with unilateral carotid artery stenoses can be identified with a high level of certainty by using conventional MRA, which has a high negative predictive value (1, 38, 39). However, a restriction is that DSLA cannot be performed after the administration of contrast agent.
The greatest disadvantage of all spin labeling techniques is their poor signal-to-noise ratio, which markedly improves at high magnetic field strength (e.g., 3T). The gain for spin-labeling techniques is twofold: Magnetization in thermal equilibrium, and thus the signal intensity is higher, and decay of the magnetic label is slower. T1 relaxation time of blood is about 1650 ms at 3T compared with about 1200 ms at 1.5 T. The wider availability of MR imagers with high field strength provides new options for the use of spin labeling techniques. The FLASH readout we used to perform DSLA does not yield optimal results in terms of signal-to-noise ratio. Double that ratio could be achieved by using segmented spiral readout at identical acquisition time and temporal resolution (40).
Conclusion
DSLA is an auspicious method for noninvasively quantifying the hemodynamic effect of an extracranial carotid artery stenosis and the resulting intracranial collateralization. It cannot replace morphologic MRA, but it may serve as a supplementary examination to enhance specificity in grading carotid artery stenoses, to monitor patients after intervention, and to provide additional functional information.
Footnotes
Supported by the German Research Foundation (Deutsche Forschungsgesellschaft, DFG), grant GRK 238/2.
References
- Received May 28, 2004.
- Accepted after revision October 6, 2004.
- Copyright © American Society of Neuroradiology