(Stroke. 2003;34:77.)
© 2003 American Heart Association, Inc.
Original Contributions |
From the Departments of Neurology (J.E., M.R., S.M., H.P., T.P.), Electrical Engineering (W.W., H.E.), and Radiology (G.S.), Ruhr University Bochum, Bochum, Germany.
Correspondence to Dr Jens Eyding, Neurologische Klinik der Ruhr-Universität, St Josef Hospital, Gudrunstrasse 56, D-44791 Bochum, Germany. E-mail jeyding{at}web.de
| Abstract |
|---|
|
|
|---|
Methods In a system with a constant concentration of contrast agent, analyzing the decrease in image intensity that occurs with microbubble-destructive imaging modes yields parameters that are considered to correlate with tissue perfusion. This method was first evaluated with a perfusion model that showed that the main resulting parameter "perfusion coefficient" (PC) is a monotonic nonlinear function of flow velocity. Seventeen human volunteers were then scanned according to this method with the use of 2 different contrast agents. Results were correlated with those from perfusion-weighted MRI examinations.
Results The PC did not show significant differences in gray matter areas (ranging from 1.466x10-2 s-1 to 1.641x10-2 s-1) of the brain despite different insonation depths (eg, ipsilateral and contralateral thalamus). In contrast, white matter exhibited significantly lower perfusion values in both imaging modes (PC: 0.604x10-2 s-1 to 0.745x10-2 s-1; P<0.05).
Conclusions CODIM is a promising new tool of imaging parenchymal (brain) perfusion in healthy persons. The method provides semiquantitative and depth-independent perfusion parameters and in this way overcomes the limitations of the perfusion methods using a bolus kinetic. Further investigations must be done to evaluate the potential of the method in patients with perfusion deficits.
Key Words: cerebrovascular circulation magnetic resonance imaging perfusion ultrasonography
| Introduction |
|---|
|
|
|---|
The purpose of this study was to investigate the potential of a new semiquantitative ultrasonic method of imaging parenchymal perfusion in healthy volunteers, referred to as contrast burst depletion imaging (CODIM), with the use of Levovist and Optison as contrast agents. This method makes use of the decrease in contrast agent concentration during destructive imaging to combine the advantages of high-contrast, highsignal-to-noise ratio imaging techniques with the semiquantitative information and easy-to-access information that the replenishment method may also provide. To analyze and validate the proposed method, 3 steps were undertaken: in vitro evaluation with the use of a perfusion phantom, in vivo evaluation with examination of healthy volunteers, and comparison of the method with perfusion-weighted MRI (PWI), ie, the present "gold standard" of brain perfusion measurement in clinical settings.8 The focus of this study was chosen with respect to the diagnosis of perfusion deficits in patients with acute ischemic stroke.
| Subjects and Methods |
|---|
|
|
|---|
Technical Background of CODIM
CODIM is an algorithm that, in combination with microbubble-destructive imaging techniques, estimates perfusion, with some constraints, independent of imaging depth, ie, independent of tissue attenuation. The major constraint is that the contrast imaging method must provide a measure of microbubble concentration. To account for the poor imaging conditions in the brain, a very sensitive contrast agent imaging technique, ie, CBI, was chosen. Nevertheless, the penetration depth is less than that for the morphological B-mode image, and the near field up to 20 to 30 mm is not accessible. A prerequisite for this method is a stable concentration of microbubbles within the circulation during data acquisition achieved by either infusion or bolus application, with a time delay between injection and insonation. In case of contrast agent infusion, the infusion rate must be adjusted carefully to first establish a sufficient microbubble concentration in the blood pool and then to compensate for the slow decay in concentration. According to our experience, a fairly constant microbubble concentration can be assumed after a bolus injection of contrast agent some time after the maximum concentration has been reached because the life span of microbubbles in the blood pool is far longer than the time of data acquisition. For the aforementioned reasons, bolus injection is even more reliable than infusion unless the microbubble concentration is being monitored and the infusion rate is adjusted accordingly. Therefore, we chose the bolus injection for all examinations. After 40 seconds without insonation, imaging using a destructive image mode (eg, CBI) with a high pulse repetition frequency is started. The image intensity in perfused tissue will decrease and approach an equilibrium, which is defined by the destruction by ultrasound and reperfusion of microbubbles (Figure 1). At least 10 images should be acquired, covering a decrease in intensity of 6 to 30 dB in the perfused tissue within the first 5 images, depending on the local perfusion rate. Under these conditions, the number of measurements is greater than the number of parameters that are estimated, and the changes in the measured intensities are also greater than those that might arise from noise and other artifacts. The transmit power and the frame rate must be chosen accordingly. The evaluation of the quantitative data corresponding to the sample volume is based on fitting the measured time-intensity curve to a theoretical model. The concentration of microbubbles [c(1)] in a sample volume at the time of the first insonation is given by the concentration of microbubbles in the blood pool weighted by the blood-to-tissue ratio in the sample volume (Vsample), as follows
|
|
|
For any given sample volume it is assumed that the blood supply does not originate primarily from regions within the same imaging plane. Destruction of microbubbles occurs during or immediately after the insonation and will be observable at high pulse repetition rates. However, compared with the time between 2 frames, the destruction can be considered instantaneous. Between 2 frames, the blood with the destroyed microbubbles is exchanged so that the concentration of microbubbles will approach the concentration in the blood pool c0 again following an exponential curve
|
|
where n=number of insonations,
t=1/frame rate, d=destruction coefficient, and p=perfusion coefficient. Concentration can only be observed at discrete positions in time denoted by the number of insonations n. For concentration as a function of insonations, a closed form solution exists
|
|
Given a sufficient number of measurements with reasonable signal-to-noise ratio, a least squares fit of the model to the measurement yields all parameters of the model. The residual error of the least squares fit is also determined to quantify the reliability of the parameters. The most critical precondition of this type of evaluation is to find a function that converts image intensity to a parameter that is proportional to the concentration of microbubbles. Image intensity is considered to be an arbitrary unit that may be (nonlinearly) rescaled to be proportional to microbubble concentration. The factor of proportionality may still vary within the imaging plane. Maps of c0 will therefore still reflect depth-dependent variations that are typical for ultrasonic imaging. e-d specifies the normalized amount of microbubbles after an insonation. p characterizes an exchange rate, so that its unit is s-1. In the following, the perfusion coefficient p will be termed PC and the destruction coefficient d will be termed DC.
In Vitro Experiments
CODIM was initially investigated in vitro with the use of a perfusion phantom. The phantom consisted of a cylindrical sponge with open pores in the submillimeter range. The sponge had a diameter of 24 mm and a length of approximately 80 mm and was embedded in agar at a mean depth of 50 mm. A rotary pump controlled the flow of contrast agent in the direction of the cylinder axis. According to the size of the phantom, a 7.5-MHz linear array was used for scanning at a frame rate of 2 Hz. The experiments covered a flow velocity range of 0 to 10 mm/s. Transmit power, frame rate, and flow velocity were varied to model different flow and insonation conditions. The experiments confirmed a correlation between flow velocity and the perfusion coefficient as defined in Equation 2. Figure 2 shows an example for Levovist, 100% transmit power, 2 frames per second.
|
In Vivo Examinations
The present study included 17 volunteers with sufficient temporal bone window (9 women and 8 men; mean age, 36.06 years; range, 21 to 51 years). Five volunteers were examined with Levovist alone, 10 with both Levovist and Optison, and 2 with Optison alone (resulting in 15 Levovist and 10 Optison examinations). Exclusion criteria were history or physical signs of cerebrovascular diseases, diseases that might influence brain perfusion, galactosemia, heart insufficiency, known allergy, and pregnancy. Usual extracranial and transcranial ultrasonic examinations had to be normal according to previously published criteria.9 MRI and ultrasonic examinations were performed within 7 days. Informed consent was obtained before subjects entered the study. A positive vote of the local ethics committee was obtained. The procedures followed in this study were in accordance with institutional guidelines. Statistical analysis was performed with the use of ANOVA for repetitive measurements when similarity of values was expected and with the Wilcoxon test for dependent samples when difference of values was expected. The level of significance was defined as 0.05 (2-sided).
Protocol of In Vivo Examinations
A Siemens Sonoline Elegra ultrasound system equipped with a 2.5-MHz phased array transducer (Siemens Medical Systems, Inc, Ultrasound Group) was used for all in vivo examinations, which were performed unilaterally with the transtemporal approach. The field of view was set to an imaging depth of 100 mm and a sector angle of 90 degrees. Ultrasound examinations were performed in CBI mode (a destructive imaging technique that visualizes the ultrasound-induced effects on microbubbles; see above and Reference 4) after intravenous bolus injection of either a galactose- and palmitic acidbased microbubble suspension (Levovist, Schering AG [400 mg/mL], 8 mL, injection speed 2 mL/s) or a perfluropropane-filled agent (Optison, Mallinckrodt GmbH [0.5 mL], injection speed 0.5 mL/s) through an 18-gauge venous access into an antecubital vein. Injection started once the diencephalic plane was correctly adjusted with the third ventricle used as a landmark of orientation. The adjacent thalamus was visualized as a semicircular hypoechogenic structure. Other parenchymal structures were depicted according to previously suggested criteria.1 The ultrasound system was frozen to stop any insonation of the field of view. Data acquisition was started after 40 seconds at a frame rate of 1 Hz and acoustic power set to 63% (mechanical index 1.8). Morphological images were acquired in parallel in the tissue harmonic imaging mode for orientation. Intervals between different examinations in 1 volunteer were 10 minutes. With this protocol, 15 examinations in which Levovist was used and 10 examinations in which Optison was used were undertaken.
Evaluation of Quantitative Parameters With CODIM
Ultrasound data were transferred to a personal computer, where offline evaluation of raw B-mode and CBI data (beam-formed B-mode or power-mode data before any other processing such as scan conversion, dynamic range compression, or thresholding) was performed with the use of dedicated software developed by our group. Manually placed rectangular ROIs were as follows: ipsilateral anterior thalamus, ipsilateral posterior thalamus, ipsilateral head of caudate nucleus, ipsilateral lentiform nucleus, ipsilateral white matter, contralateral anterior thalamus, contralateral posterior thalamus, and contralateral head of caudate nucleus. Defining a ROI for "white matter" is somewhat controversial because white matter and gray matter exhibit identical echogenicity in B-mode (tissue harmonic imaging mode) images. For this reason it is necessary that experienced transcranial B-mode examiners place this specific ROI using the lentiform nucleus and the insonation depth as landmarks of orientation. Size of ROIs varied as previously described,4 ie, for head of the caudate nucleus, anterior thalamus, and posterior thalamus between 20 and 40 mm2 and for lentiform nucleus and white matter between 70 and 100 mm2. For each ROI, perfusion coefficients (PC) and destruction coefficients (DC) were determined as described above.
MRI Examinations
Protocol of MRI Examinations
MR measurements were performed on a 1.5-T clinical whole-body MR scanner (Magnetom Symphony Quantum, Siemens AG) equipped with a gradient overdrive, with the use of the standard head coil. In addition to the ordinary examination protocol (T1, T2), an axial echo-planar PWI sequence was performed. The parameters of the perfusion study with gradient echo were as follows: 7 to 9 slices, slice thickness 5 to 6 mm, interslice gap 20% of slice thickness, echo time 45 ms, field of view 250 mm, and matrix 128x128. Fifty T2*-weighted measurements per slice at intervals of 1.5 seconds were taken. The contrast agent (0.2 mmol/kg Gd-DTPA) was injected at a rate of 5 mL/s.
Evaluation and Postprocessing Image Analysis
PWI data were transferred to a separate personal computer and computed pixel by pixel to create mean transit time, regional cerebral blood volume (rCBV), and regional cerebral blood flow (rCBF) maps with the use of MEDx 3.3 software (Sensor Systems Inc). ROIs were placed into anterior and posterior thalamus, lentiform nucleus, head of the caudate nucleus, white matter, and precentral gyrus of both sides to obtain semiquantitative data (rCBV and rCBF in arbitrary units; regional mean transit time in seconds). ROI sizes were defined in a manner similar to that used for the ultrasound examination, which was guaranteed because the examination was undertaken by the same investigator. The size of precentral gyrus ROIs varied between 20 and 40 mm2.
| Results |
|---|
|
|
|---|
In Vivo Examinations
All in vivo examinations were well tolerated; no side effects were observed for either contrast agent except for local side effects in 4 volunteers using Optison. For practical reasons, 16 of 17 volunteers were examined from the left side; 1 volunteer had a sufficient bone window only on the right side. For this reason data were not referred to as "right side" or "left side" but as "ipsilateral" (ie, in 16/17 cases "left") and "contralateral," which means at the same time "close to" and "distant from" the transducer. All individuals had an adequate bone window that allowed visualization of the third ventricle and the thalamus. Of the 17 volunteers, 5 were examined with Levovist alone, 10 with both Levovist and Optison, and 2 with Optison alone (resulting in 15 Levovist and 10 Optison examinations).
With the use of Levovist, in 2 examinations no depletion curve was obtained in contralateral anterior thalamus and contralateral posterior thalamus, and in 1 examination no contrast enhancement could be achieved in contralateral head of caudate nucleus, ipsilateral head of caudate nucleus, ipsilateral anterior thalamus, ipsilateral posterior thalamus, and ipsilateral white matter. With the use of Optison, the acquisition of usable data failed in 1 examination in contralateral head of caudate nucleus, contralateral anterior thalamus, and contralateral posterior thalamus. In all other ROIs, sufficient data, as shown in Figure 3, could be obtained (a total of 178 of 200 ROIs in 25 examinations).
|
CODIM
Quantitative results of CBI and PWI are summarized in Tables 1 to 4, with relevant data combined for comparison of PC and rCBF values in Tables 1 and 2. Clearly, similar PC values in the following regions of gray matter independent of depth of examination were found: ipsilateral anterior thalamus, ipsilateral posterior thalamus, ipsilateral head of caudate nucleus, ipsilateral lentiform nucleus, contralateral anterior thalamus, contralateral posterior thalamus, and contralateral head of caudate nucleus. The range of mean values of PC in these ROIs was 1.466 to 1.633x10-2 s-1 with Levovist and 1.505 to 1.641x10-2 s-1 with Optison. When we compared PCs in these ROIs statistically (with similarity expected), no significant differences could be found in either group. On the other hand, the PC of ipsilateral white matter was significantly lower (with difference expected) than the PC of all other ROIs with both Levovist and Optison (mean, 0.6042 and 0.7446x10-2 s-1, respectively).
|
|
|
|
For better comparison of these data with PWI data, we calculated (in accordance with established MR criteria10) relative indices by dividing values of identical ROIs from both sides ("mirror regions"). Since contralateral lentiform nucleus values were not obtained by the described ultrasonic method because of shallow insonation depth, we divided ipsilateral lentiform nucleus values by contralateral anterior thalamus values to obtain a relative value for ipsilateral lentiform nucleus. This was only possible because PWI data showed no significantly differing rCBF values between contralateral anterior thalamus and contralateral lentiform nucleus (see below). Data shown in Table 4 indicate that there were no significant differences between different mirror regions.
To evaluate whether the PC depended on the contrast agents used, values of those subjects who were examined with both Levovist and Optison were analyzed. No significant differences between any of the examined ROIs could be found (with the use of both ANOVA and the Wilcoxon test).
DC values, representing acoustic energy in the field of view, showed a great range with a large SD (Table 3). It was not possible to outline similarities or differences between the different ROIs in either direction statistically.
MRI Examinations (PWI)
MRI examinations were also well tolerated without side effects; in all examinations sufficient data for postprocessing analysis were obtained, with a total of 204 of 204 ROIs in 17 examinations. Figure 4 shows typical parameter images for rCBF and rCBV of the adequate plane as well as a time-intensity curve with the ROI placed into the "ipsilateral" posterior thalamus. In accordance with the individual ultrasound data (see above), "right" and "left" were changed to "ipsilateral" and "contralateral."
|
Tables 1 and 2 summarize rCBF values of the PWI examination. Mean values of rCBF in regions of gray matter (ipsilateral anterior thalamus, ipsilateral posterior thalamus, ipsilateral head of caudate nucleus, ipsilateral lentiform nucleus, contralateral anterior thalamus, contralateral posterior thalamus, contralateral head of caudate nucleus, contralateral lentiform nucleus) ranged from 6.797 to 7.869x10-3 and from 8.657 to 9.772x10-3 in the 2 groups (ie, volunteers who had been examined with either Levovist or Optison, respectively); significant differences between these ROIs (with similarity expected) could not be found. Comparison of the calculated mirror region values did not show significant differences (Table 4). Mean values of rCBF for ipsilateral and contralateral white matter in both groups (3.024 versus 3.987x10-3 and 3.161 versus 4.159x10-3, respectively) were significantly lower than for the aforementioned ROIs (with difference expected), and mean values of ipsilateral and contralateral precentral gyrus were significantly higher in both groups (10.36 versus 13.18 and 10.33 versus 12.88x10-3, respectively).
Comparison of CODIM and PWI
Finally, we compared CODIM PC values with PWI rCBF data. Between the calculated relative indices (mirror regions) for the regions, we did not find significant differences in the 2 groups (with both similarity and difference expected; Table 4).
| Discussion |
|---|
|
|
|---|
In the new CODIM method, the "depletion" technique is based on the assumption that insonation decreases the concentration of a specific contrast agent depending on the potential of refilling the insonated area with the agent. With the use of specific algorithms that consider the ultrasound-induced destruction and reperfusion, the semiquantitative parameter PC can be mathematically calculated from the resulting time-intensity curves of the CBI examination. Given that the theoretical model applies, this parameter is, within a certain range, independent of sound intensity within the ROI and therefore independent of depth, which would be a major advantage compared with all other established methods using the bolus method. Furthermore, examination can be shortened to 60 seconds, ie, 40 seconds to establish a constant microbubble concentration in the blood pool and the imaging plane and another 20 seconds for the actual acquisition of 20 frames at 1 Hz.
The aim of the present study was to evaluate the diagnostic potential of this method. As a first step we examined a gel perfusion phantom. Using the flow rate of the model as a referential factor, we demonstrated that different flow rates correlated well with resulting PC values in a nonlinear relation. The fact that PC was non-zero for the no-flow condition is most likely due to 2 factors: the signal-to-noise ratio is poor for very low concentrations of contrast agent, and the agent can still circulate within the sponge. This circulation will occur because, for thermodynamic reasons, the systems will approach equilibrium with respect to temperature and microbubble concentration.
Next we compared the method with PWI in healthy volunteers. Resulting PC values did not differ significantly in any area of similar perfusion (ie, gray matter). The small size of both groups, however, limits statistical strength. PC in white matter was significantly lower by a factor of approximately 0.4 in both the Levovist and Optison groups. Relative values, ie, ratios of values in mirror regions, were approximately 1, indicating again that PC was independent of depth of examination. At the same time, PC was stable when different contrast agents were used, although the stability of the bubbles with respect to ultrasound exposure was quite different. As a probable marker for these differences, the second parameter, DC, showed high variations intraindividually and interindividually.
Finally, comparison with the established clinical "gold standard" PWI revealed that relative values of mirror regions of both methods showed no significant differences, which strongly suggests that both methods are equally independent of depth of examination. Still, it remains striking that PC values had a distinctly smaller range of parameter values interindividually.
Thus, CODIM seems to be a potent new method for semiquantitative assessment of parenchymal cerebral perfusion in healthy volunteers. As expected from former research, acoustic energy was inhomogeneous in the insonation field, as indicated by the DC. Nevertheless, the proposed data acquisition and model-based mathematical analysis could eliminate the depth dependency of former ultrasonic bolus methods (for which wash-in/wash-out method would be a better description), which is a major disadvantage of these methods, and facilitate the diagnostic procedure. The resulting data, although of a small group, may help to create "normal values" of brain perfusion. Furthermore, PC values of white matter were significantly lower than those of gray matter, suggesting the possible ability of CODIM to detect different perfusion rates. Therefore, the method should be tested not only in more healthy volunteers but also in patients with hypoperfused areas (eg, acute stroke), where the relative values could possibly detect regions of ischemia in agreement with PWI criteria.10 Furthermore, it would be of significant interest to compare the replenishment approach with CODIM in healthy volunteers as well as in acute stroke patients to obtain knowledge about the robustness of both methods under clinical conditions.
Some basic limitations of ultrasonic methods also apply for CODIM. Patients must cooperate, since the transducer must be held firmly without loss of the bone window (however, the correct positioning can be monitored during the examination with help of B-mode images on the screen). Second, to date no criteria for independently determining size and location of the chosen ROIs have been described, and therefore experienced examiners are needed to define comparable ROIs. Third, patients must have an adequate temporal bone window.
In conclusion, the present study introduces a new ultrasonic tool of semiquantitative perfusion imaging that was shown to be fairly independent of depth parameters and of physical properties of contrast agents. Moreover, CODIM reduces the time of examination considerably compared with other established methods.
| Acknowledgments |
|---|
Received May 10, 2002; revision received July 30, 2002; accepted August 8, 2002.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
O. Catalano, V. Migaleddu, E. Quaia, and G. Caruso Terminology for Contrast-Enhanced Sonography: A Practical Glossary J. Ultrasound Med., June 1, 2007; 26(6): 717 - 730. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. M. Porter, D. A. B. Smith, and C. K. Holland Acoustic Techniques for Assessing the Optison Destruction Threshold. J. Ultrasound Med., December 1, 2006; 25(12): 1519 - 1529. [Abstract] [Full Text] [PDF] |
||||
![]() |
O. Lucidarme, Y. Kono, J. Corbeil, S.-H. Choi, J.-L. Golmard, J. Varner, and R. F. Mattrey Angiogenesis: Noninvasive Quantitative Assessment with Contrast-enhanced Functional US in Murine Model. Radiology, June 1, 2006; 239(3): 730 - 739. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Eyding, W. Wilkening, C. Krogias, T. Holscher, H. Przuntek, S. Meves, and T. Postert Validation of the Depletion Kinetic in Semiquantitative Ultrasonographic Cerebral Perfusion Imaging Using 2 Different Techniques of Data Acquisition J. Ultrasound Med., August 1, 2004; 23(8): 1035 - 1040. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Kern, F. Perren, K. Schoeneberger, A. Gass, M. Hennerici, and S. Meairs Ultrasound Microbubble Destruction Imaging in Acute Middle Cerebral Artery Stroke Stroke, July 1, 2004; 35(7): 1665 - 1670. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Stroke Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2003 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |