|
|
||||||||
Experimental Studies |
1 From the Department of Radiology, Division of MRI, Brigham and Womens Hospital, Harvard Medical School, 221 Longwood Ave, Longwood Medical Research Center, 007c, Boston, MA 02115 (N.J.M., R.L.K., F.A.J., K.H.H.), and the Department of Physics and Astronomy, Tufts University, Medford, Mass (N.J.M.). Received September 8, 1999; revision requested October 29; revision received November 29; accepted January 12, 2000. Supported by National Cancer Institute grants CA 46627 and P01 CA 67165 and contract 282-97-0080 from the U.S. Department of Health and Human Services. Address correspondence to N.J.M. (e-mail: njm@bwh.harvard.edu).
| ABSTRACT |
|---|
|
|
|---|
MATERIALS AND METHODS: Sonications were delivered in rabbit thigh muscles at varying powers. Temperature-sensitive MR images obtained during the sonications were used to estimate the temperature and thermal dose. The temperature, thermal dose, and applied power were then correlated to the occurrence of tissue damage observed on postsonication images. An eight-element phased-array transducer was used to produce spatially flat temperature profiles that allowed for averaging to reduce the effects of noise and the voxel size.
RESULTS: The occurrence of tissue damage correlated well with the MR imagingderived temperature and thermal dose measurements but not with the applied power. Tissue damage occurred at all locations with temperatures greater than 50.4°C and thermal doses greater than 31.2 equivalent minutes at 43.0°C. No tissue damage occurred when these values were less than 47.2°C and 4.3 equivalent minutes.
CONCLUSION: MR imaging thermometry and dosimetry provide an index to predict the threshold for tissue damage in vivo. This index offers improved online control over minimally invasive thermal treatments and should allow for more accurate target volume coagulation.
Index terms: Magnetic resonance (MR), thermometry, 44.121412 Ultrasound (US), experimental studies, 44.1298 Ultrasound (US), focused, 44.1298 Ultrasound (US), therapeutic, 44.1298
| INTRODUCTION |
|---|
|
|
|---|
Although MR imaging can depict tissue damage created by thermal therapies, in some cases it may be impossible to do so. For example, it may be hard to distinguish thermally damaged tissue in tumors from the tumor itself (11) and the edema that surrounds it. It also may not be desirable to cause tissue coagulation at all, as in hyperthermia treatments, for example. Thus, having another method beyond posttreatment imaging to determine the extent of tissue damage is necessary for optimal thermal treatments.
Three methods commonly are used to predict whether tissue damage induced by thermal therapies will occur. The first method uses the power output of the treatment device and the exposure length to predict tissue damage (1418). The second method postulates a critical temperature, Tc, above which tissue damage occurs (19). The third method uses the entire temperature history of the tissue to estimate the thermal dose (20).
Temperature and thermal dose measurements have the advantage of being independent of inhomogeneities in the acoustic, optic, or electric properties of tissue and rates of perfusion. Thermal dose measurements have the advantage of being independent of the length and shape of the temporal heating profile. Thermal dose typically has been used in hyperthermia treatments, in which the heating profile is long and stays near 43.0°C, but also has been suggested for use during short-duration, high-temperature, focused ultrasound surgery (21). Others have investigated the use of MR imagingderived dosimetry for predicting the size of lesions produced during focused ultrasound surgery (2224). The purpose of this study was to investigate whether MR imagingderived thermometry and dosimetry of thermal exposures can be used to quantify the threshold of tissue damage in vivo.
| MATERIALS AND METHODS |
|---|
|
|
|---|
The sector-vortex array typically is phased in several modes. In mode zero, each element of the transducer is in phase and is equivalent to a single-element transducer. This mode was used at a power high enough to localize the focus but low enough not to induce tissue damage (10). In mode four, alternating elements of the array are 180° out of phase. In this mode, the ultrasound field is canceled out at the natural focus, and eight peaks are created that are arranged radially about the direction of the ultrasound beam in the shape of a ring. As the ultrasound is delivered, conduction fills in the center of this ring and produces a spatially flat temperature profile. This mode was used in the experiments to create a large focal volume so that the induced temperature rise could be averaged over a section of tissue. A hydrophone scan of the focal plane in mode four is shown in Figure 1. The sonications were 30 seconds in duration. The characteristics of the ultrasound field were measured with a previously described method (1) and are summarized in Table 1.
|
|
The rabbits body temperatures were monitored constantly with copper-constantan thermocouples that were constructed in-house. One thermocouple was inserted rectally. In some experiments, another thermocouple was implanted in the thigh muscle a few centimeters away from the location of the sonications. The temperature in the thigh was used as the baseline temperature in the experiments. When the thermocouple in the thigh was not used, the rectal thermocouple temperature was used as the baseline.
Experimental Setup
The transducer was mounted in a computer-controlled, mechanically driven, MR imagingcompatible positioning system (GE Medical Systems, Milwaukee, Wis). An earlier version of the positioning system has been described (1,2). The rabbits thigh was placed in a hole in a plastic tray positioned above the transducer. The hole was the opening of a bag of deionized, degassed water. The bag rested on a polyvinyl chloride membrane above the transducer that also was submerged in degassed water. To provide proper coupling, a layer of degassed water was poured between the membrane and the plastic bag. A receive-only surface coil with a diameter of 12.7 cm was attached beneath the plastic tray below the thigh to improve the signal-to-noise ratio.
To maintain its body temperature, the rabbit was placed on a plastic mat that circulated heated water. A plastic coil that circulated warm water encircled the water bag to ensure further that the thigh was kept warm during the experiments and to avoid temperature gradients in the tissue. The temperature in the water was monitored with a thermocouple and was kept as close to the body temperature of the rabbit as possible. The whole apparatus was put into a clinical 1.5-T MR imager (Signa; GE Medical Systems). Figure 2 shows a diagram of the experimental setup.
|
times the time in which the phase developed (the echo time of the image). A fast spoiled gradient-echo sequence was used to acquire the phase maps (31). The following parameters were used: 50.6/25 (repetition time msec/echo time msec); flip angle, 30°; bandwidth, 3.1 kHz; field of view, 16 cm; section thickness, 3 mm; matrix size, 256 x 128; imaging time per image, 6.7 seconds. An echo time of 25 msec was used because the optimal signal-to-noise ratio in fast spoiled gradient-echo phase-difference imaging occurs when the echo time is equal to T2*, and the T2* of muscle is around 25 msec (31). The small temperature gradients and the relatively long sonication time allowed for the long imaging time required when the optimal echo time was used. The temperature dependence of the proton-resonant frequency shift in rabbit skeletal muscle is 0.00909 ppm/°C and is linear above the tissue coagulation threshold temperature (24). With an echo time of 25 msec and a magnetic field strength of 1.5 T, the temperature dependence of the sequence was 11.0°C/radian. Twelve sets of magnitude, real, and imaginary images were acquired, which yielded a total imaging time of 80 seconds. A reference phase map was acquired 7 seconds before the application of the ultrasound. The phase-difference images were constructed from this reference image by using a complex phase-subtraction method (31).
To avoid phase wraparound, the images were subtracted pairwise, and then the phase was summed along the time dimension. For the nth phase subtraction 
n:
i was the phase at the ith time point. This method prevented phase wrapping as long as the temperature did not cause a phase change greater than
radians in the time of one image acquisition (
T
34.5°C in 6.7 seconds). The images were acquired in a plane perpendicular to the direction of the ultrasound beam to reduce the effects of averaging of the temperature over the 3-mm section thickness. The position of the focus was found by localizing the maximum voxel value in the region of the focus. The temperature at each voxel was averaged with its eight nearest neighbors to reduce the effects of noise.
Thermal Dose Calculation
As proposed by Chung et al (22), the Sapareto and Dewey equation (20) was used to estimate the thermal dose delivered to the tissue from the time-temperature profiles gathered by means of MR imaging thermometry:
t, and R is a constant that equals 0.50 above 43.0°C and 0.25 below 43.0°C. The temperature profiles of each voxel were interpolated linearly with a step size of
t = 0.1 second.
Postsonication Imaging
After the sonications, T2- and T1-weighted fast spin-echo images were acquired for about 2 hours to determine whether there was a tissue response to the applied sonications. The parameters used in these images are shown in Table 2. The T1-weighted images were acquired before and after injection of a bolus of gadopentetate dimeglumine (Magnevist; Berlex Laboratories, Wayne NJ) at a concentration of 0.3 mmol/kg in each rabbits ear vein. After the injection, images were obtained until the enhancement reached its maximum. Contrast material was administered shortly after the last sonication in each thigh, and again roughly 2 hours after the last sonication.
|
| RESULTS |
|---|
|
|
|---|
|
|
The noise in the image subtractions was calculated by taking four 3 x 3-voxel regions away from the ultrasound focus for each sonication. The regions were chosen to be away from blood vessels, muscle boundaries, and fat. The mean and SD of these voxel regions corresponded to a temperature change of 0.8°C ± 0.3 (n = 3,916). The mean and SD of all of the individual voxels in these regions corresponded to -0.1°C ± 1.0 (n = 35,244). There was no noticeable increase in noise away from the focus when the ultrasound was turned on. At peak temperature rise, the mean SD of a 3 x 3-voxel region centered at the focus was 1.0°C ± 0.8 (n = 91). This corresponded to 8.4% ± 4.8 of the mean temperature rise, which was 12.4°C ± 2.9.
The thermocouple implanted in the thigh was not used during 31 of the 91 sonications analyzed in the experiment. In these cases, the rectal temperature was used for the baseline temperature. The mean difference between the rabbit body temperature measured with the thermocouple in the thigh and with that in the rectum was 0.4°C ± 1.0; the maximum difference was 2.5°C.
Postsonication Images
Figure 5 shows examples of the T2-weighted and contrast-enhanced T1-weighted images of lesions created by the focused ultrasound heating. These lesions were caused by 13.9-, 15.4-, and 21.6-W sonications. The lesions were characterized by a dark region surrounded by a bright ring on T2-weighted images and by a region of no contrast enhancement surrounded by a ring of strong enhancement on the T1-weighted images. The size of the lesions did not vary much during the 2 hours after the delivery of the sonications. However, the signal intensity increased on the T2-weighted images for some of the lesions during this time.
|
|
| DISCUSSION |
|---|
|
|
|---|
The indeterminate range of dose measurements (~27 minutes) and temperature measurements (~3°C) was expected because of the accuracy of the MR imaging temperature estimates and the baseline temperatures in the muscle. Another reason for this indeterminate region may have been because of inhomogeneities in the ultrasound field. Even though the sector-vortex array produced flat temperature profiles after some time, inhomogeneities in the temperature before this point could have caused nonuniform dose profiles. To produce truly flat dose profiles with the sector-vortex array, it would be necessary to switch temporally between different modes (26). However, it is expected that the histologic variations between animals and tissue locations were a substantial source of the variation observed in the study. The results from this study indicate that when all of these uncertainties are considered together, MR imagingderived thermometry and dosimetry provide an index that can be used to estimate the threshold for tissue damage.
MR imagingderived thermometry and dosimetry offer a distinct advantage over previous in vivo methods used to determine the threshold for tissue damage. A whole plane, or volume, can be imaged by means of MR imaging, whereas previous studies used invasive temperature probes with one or a few measurement points that were insensitive to spatial temperature inhomogeneities.
Results of previous studies (3337) on the effects of hyperthermia have shown that the threshold for 100% necrosis in different tissues ranges 50240 minutes. Meshorer et al (33) studied the tissue effects in porcine muscle that occurred after 30 minutes of heating. The acute effects observed 1824 hours after therapy were that minimal necrosis occurred at 40°C43°C, moderate damage occurred at 44°C46°C, and severe damage occurred at temperatures greater than 46°C. These relationships between time and temperature convert to 0.530.0 minutes at 43°C for minimal necrosis, 60240 minutes for moderate damage, and more than 240 minutes for severe damage. Those results agree well with the results found in this study.
According to the hyperthermia literature (21,24,38), this dose threshold of 240 minutes is adequate to coagulate all tissue and thus has been used as an indicator of tissue damage caused by focused ultrasound. The results from this study indicate that 240 minutes is a conservative estimate for predicting tissue damage in muscle tissue. However, more studies, especially in tumor tissue, are required to truly establish the treatment threshold for tumor surgery.
In conclusion, the results of this study show that MR imagingderived thermometry and dosimetry can be used to assume that an adequate thermal exposure is delivered to coagulate the target tissue. With advances in MR imaging thermometry, the threshold presumably could be found with more precision. Knowledge of the threshold exposure necessary to cause tissue damage along with online temperature monitoring could allow for optimal thermal treatments. In addition, use of MR imaging thermometry could allow thermal exposures that are close to the thresholds to be explored for therapeutic purposes. An example of such exposure is increasing the blood vessel permeability for targeted drug or gene therapy (39).
Practical applications: Several ultrasound (4044), laser (3,45,46), and radio-frequency (4,47) thermal surgery devices currently are being tested clinically. The ability to predict the induced tissue damage during treatments with these devices is important for several reasons. The foremost reason is to provide online information to the physician so that the desired target volume is exposed adequately. Second, the thermal imaging and dose calculation can increase safety by providing warning before the surrounding normal tissue is overexposed. Third, the ability to monitor the thermal dose removes the need to overexpose tissue. This reduces the treatment time and the danger of causing thermal buildup in the overlying tissue (4850).
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. Adachi, T. Kaminou, T. Ogawa, T. Kawai, Y. Takaki, K. Sugiura, Y. Ohuchi, and M. Hashimoto Heat Distribution in the Spinal Canal during Radiofrequency Ablation for Vertebral Lesions: Study in Swine Radiology, May 1, 2008; 247(2): 374 - 380. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. D. O'Brien Jr, C. X. Deng, G. R. Harris, B. A. Herman, C. R. Merritt, N. Sanghvi, and J. F. Zachary The Risk of Exposure to Diagnostic Ultrasound in Postnatal Subjects: Thermal Effects J. Ultrasound Med., April 1, 2008; 27(4): 517 - 535. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Kopelman and M. Papa Magnetic Resonance-Guided Focused Ultrasound Surgery for the Noninvasive Curative Ablation of Tumors and Palliative Treatments: A Review Ann. Surg. Oncol., May 1, 2007; 14(5): 1540 - 1550. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. Teichgraber, C. Aube, D. Schmidt, E. Jehle, C. Konig, C. D. Claussen, and P. L. Pereira Percutaneous MR-guided radiofrequency ablation of recurrent sacrococcygeal chordomas. Am. J. Roentgenol., August 1, 2006; 187(2): 571 - 574. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. McDannold, C. M. Tempany, F. M. Fennessy, M. J. So, F. J. Rybicki, E. A. Stewart, F. A. Jolesz, and K. Hynynen Uterine Leiomyomas: MR Imaging-based Thermometry and Thermal Dosimetry during Focused Ultrasound Thermal Ablation. Radiology, July 1, 2006; 240(1): 263 - 272. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. McDannold, S. L. Fossheim, H. Rasmussen, H. Martin, N. Vykhodtseva, and K. Hynynen Heat-activated Liposomal MR Contrast Agent: Initial in Vivo Results in Rabbit Liver and Kidney Radiology, March 1, 2004; 230(3): 743 - 752. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. M. C. Tempany, E. A. Stewart, N. McDannold, B. J. Quade, F. A. Jolesz, and K. Hynynen MR Imaging-guided Focused Ultrasound Surgery of Uterine Leiomyomas: A Feasibility Study Radiology, March 1, 2003; 226(3): 897 - 905. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. E. Huber, J. W. Jenne, R. Rastert, I. Simiantonakis, H.-P. Sinn, H.-J. Strittmatter, D. von Fournier, M. F. Wannenmacher, and J. Debus A New Noninvasive Approach in Breast Cancer Therapy Using Magnetic Resonance Imaging-guided Focused Ultrasound Surgery Cancer Res., December 1, 2001; 61(23): 8441 - 8447. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Hynynen, O. Pomeroy, D. N. Smith, P. E. Huber, N. J. McDannold, J. Kettenbach, J. Baum, S. Singer, and F. A. Jolesz MR Imaging-guided Focused Ultrasound Surgery of Fibroadenomas in the Breast: A Feasibility Study Radiology, April 1, 2001; 219(1): 176 - 185. [Abstract] [Full Text] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| RADIOLOGY | RADIOGRAPHICS | RSNA JOURNALS ONLINE |