DOI: 10.1148/rg.276075097
RadioGraphics 2007;27:1809-1818
© RSNA, 2007
From the RSNA Refresher Courses
MR Imaging in Hyperthermia1
James R. MacFall, PhD and
Brian J. Soher, PhD
1 From the Department of Radiology, Duke University Medical Center, Box 3808, Erwin Road, Durham, NC 27710. Presented as a refresher course at the 2007 RSNA Annual Meeting. Received May 4, 2007; revision requested May 16 and received June 21; accepted June 25. Supported by grant P01-CA42745 from the National Institutes of Health. All authors have no financial relationships to disclose.
Address correspondence to J.R.M. (e-mail: james.macfall{at}duke.edu).
There is growing clinical evidence that the combination of radiation therapy and hyperthermia, when delivered at moderate temperatures (40°–45°C) for sustained times (30–90 minutes), is of benefit with regard to palliative relief of cancer, tumor response, local control, and survival. Adequate measurement of the temperature distribution achieved with the hyperthermia is a key element in successful therapy. Thermal dosimetry, even invasive dosimetry, is a complex topic when applied to the heterogeneous tissue of a tumor and associated organ systems. Imaging in hyperthermia therapy is performed primarily for estimation and control of temperature. Magnetic resonance (MR) imaging has unique parameter dependences that make it possible to monitor hyperthermia therapy by detection of proton resonant frequency changes or diffusion coefficient changes. In addition, MR imaging can be used to assess vascular parameters that not only allow selection of suitable patients for therapy but may also allow demonstration of response to therapy. Finally, as the use of thermally sensitive liposomes for delivery of chemotherapeutic agents is developed, MR imaging may allow determination of local drug dose.
© RSNA, 2007
There is growing clinical evidence that the combination of radiation therapy and hyperthermia, when delivered at moderate temperatures (40°–45°C) for sustained times (30–90 minutes), is of benefit with regard to palliative relief, tumor response, local control, and survival (1–3). Positive phase III clinical trials comparing radiation therapy with or without hyperthermia have been reported for several sites, including recurrent breast cancer on the chest wall (4), melanoma (5), advanced head and neck cancer (6), esophageal cancer (7), cervical cancer (8), and glioblastoma multiforme (9). In some trials where invasive temperatures were measured during treatment, retrospective analysis has shown that descriptors of temperatures achieved were significantly correlated with beneficial outcome (10–14). Other phase III trials showed no significant correlation between time-temperature relationships and outcome, despite the positive benefit seen by adding hyperthermia to radiation therapy (6,9). There were important limitations in the design of these clinical trials: (a) thermometry data were sparse, being based on invasive measurements, and (b) there was no a priori attempt to control the thermal dose delivered.
With respect to the first point above, a number of research groups have demonstrated that magnetic resonance (MR) imaging is effective and accurate for noninvasively assessing temperature changes in tissues associated with absorption of nonionizing radiation (15–21). Concerning the second point, two phase III trials (22–24) have demonstrated that controlled delivery of thermal dose is clearly related to treatment outcome. These trials did not have the limitations listed above. Nevertheless, extensive invasive thermometry and real-time manual manipulation of treatment settings were performed, which was inefficient and lacked the accuracy of data that reflected the true temperature distribution in the tumor (25). Extensive invasive thermometry cannot be done in a routine fashion, particularly in deep-seated tumors, and is a confounding factor in any trial designed to demonstrate clinical benefit from thermal therapy. The positive results of these two clinical trials clearly provide compelling rationale for continued development of noninvasive thermometry in hyperthermia.
Hyperthermia can be an effective adjunctive therapy, but adequate attention must be paid to ensuring that the target temperatures have been achieved for an adequate period of time.
Thermal dosimetry, even invasive dosimetry, is a complex topic when applied to the heterogeneous tissue of a tumor and associated organ systems. The dose is generally calculated by using invasive thermal probes such as the Luxtron (Mountain View, Calif) fiberoptic probe at five to 10 locations. The temperatures from the probes are evaluated, for example, every minute as a function of the current T90 (the temperature exceeded by 90% of the measured temperature points), and thermal doses are summed together for each minute of treatment to obtain the total dose in cumulative equivalent minutes at 43°C for the T90 temperature (CEM43°CT90) (1). As mentioned above, this methodology is challenged by the low density of temperature samples that can practically be achieved. Thus, there is a natural desire to try to achieve sampling densities closer to the pixel density of a digital image of the tumor volume. The development of image-based thermometry methods has the potential to automate and simplify hyperthermia treatment.
Regardless of the potential direct benefit of hyperthermia therapy alone or in conjunction with other therapies, another emerging role for tissue heating and thermal control is the use of temperature-sensitive liposomes to encapsulate chemotherapeutic agents such as doxorubicin to allow high-dose delivery to a tumor (http://clinicaltrials.gov/show/NCT00441376). Preclinical (26) and clinical work with doxorubicin-containing thermally sensitive liposomes also points to a strong need for improved thermal dosimetry. The doxorubicin-containing thermally sensitive liposomes show optimal release kinetics in the range between 40° and 42°C, and deviation above or below this range leads to less drug release and delivery. One concept being tested is temporal modulation of the temperature distribution across the tumor volume in order to achieve more uniform drug delivery. The implementation of volumetric thermal dosimetry and control is essential for clinical implementation of such strategies.
Also, as a mature imaging modality, MR imaging offers opportunities for more than just temperature measurement. It can provide data for treatment planning, dynamic control of treatment delivery (27), and posttreatment assessment of tissue damage (28).
Thus, while there are a variety of other imaging methods such as computed tomography, positron emission tomography, and ultrasonography that are used to stage and monitor cancer therapy, MR imaging has emerged as providing unique advantages for hyperthermia therapy. Hence, this article concentrates on MR imaging and reviews the methods for measurement of temperature with MR imaging during hyperthermia therapy.
It is this key feature of MR imaging, the ability to both monitor temperature throughout a volume and obtain useful morphometric and functional information from tumor and normal tissues, that makes it the modality of choice to image regional hyperthermia therapy delivery.
Proton Resonant Frequency Shift Technique
The fundamental thermal sensitivity of certain MR imaging parameters has been recognized for some time (29–32), but it was not thought to be a feasible method for hyperthermia thermometry because of implied incompatibilities between heating devices and the MR system. There were also concerns raised about hysteresis effects caused by tissue changes during heating (33).
The work by Delannoy et al (34) at the National Institutes of Health in 1990 was a landmark accomplishment because it showed for the first time that a clinically relevant heating device could be operated in the bore of a clinical MR imaging unit. While this work was done by using the apparent diffusion coefficient (ADC), Ishihara et al (35) performed key measurements to show the temperature dependence of the proton resonant frequency shift (PRFS) of water, which was known in the NMR literature but had not been developed for use in imaging. Also, MR temperature imaging has been demonstrated for thermal ablative procedures (28,36), where the temperatures needed to achieve tissue destruction are relatively high compared with hyperthermic temperatures.
Further developments (15–20,32,37–40) have led to the delivery of clinical hyperthermia treatments in human patients while acquiring temperature data by using MR imaging (21,41,42).
When a subject or water phantom is placed in an MR imaging unit of magnetic field strength B0, the water protons acquire a resonant frequency (fL) described by the Larmor expression: fL =
B0, where
is the gyromagnetic ratio. The PRFS method relies on the fact that this resonant frequency of protons will shift by 0.01 ppm/°C.
This is a result of the bond length between the hydrogen and the oxygen in the water getting longer as temperature increases; thus, the chemical shielding of the proton is reduced, giving increased chemical shift. For a 1.5-T imaging system where fL = 63.85 MHz, the proton in a water molecule will thus change its resonant frequency by 0.6385 Hz/°C.
If MR images are made at two different temperatures, the familiar magnitude image will not show a change, but we can also make an image where the intensity of each pixel is proportional to the image phase angle
, which depends on the pixel frequency fp (which is close to fL but shifted by temperature and field inhomogeneity) and the echo time TE, as
= fpTE. With an echo time of 20 msec, the image phase will thus change by 4.6° of image phase for each degree Celsius of temperature change at 1.5 T. It is important to remember that the method can measure only temperature change and the size of the effect depends on the echo time and magnetic field strength. While the image phase has a limited range of –180° to +180° of phase (if the real and imaginary parts of the image are available), this does not limit the method in practice if temperature updates are performed often enough that the phase change between updates does not reach the limit. The total temperature change over time is then given by summing the individual changes.
As an example of the PRFS method, Figure 1 shows a cylindrical gel phantom that is inside of a radiofrequency heating device that has four H-shaped antennae on the outside with a water coupling path between the antennae and the phantom. The phantom can be heated by energizing the antennae with, typically, 50 W of 140-MHz radiofrequency energy. This causes a temperature rise of about 1°C per minute, since the phantom is designed to have conductivity and other material parameters adjusted to approximate those of tissue. The phantom has thin catheter tubes along its length that allow MR-compatible temperature probes (Luxtron) to be inserted to measure the temperature directly. Figure 2 is a graph of the image phase as a function of the probe temperature as the phantom is heated. The inset is a phase image obtained along the axis of the phantom (coronal) that shows the positions of the temperature probes.

View larger version (100K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 1. Photograph shows a gel-filled phantom (large white arrow) inside a cylinder with four radiofrequency antennae on the outside. The small white arrow indicates one of the H-shaped antennae. There is a water-filled space (black arrow) between the antennae and the phantom.
|
|

View larger version (24K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 2. Graph shows the MR image phase as a function of directly measured temperature at six different locations in a gel phantom surrounded by a water bolus, as shown in Figure 1. The inset is a phase image obtained along the axis of the phantom (coronal view); the small white inhomogeneities are the probe locations.
|
|
A major difficulty with the PRFS method is that of image phase "drift" over the time of a series of temperature measurements that may span a hyperthermia therapy delivery session of between 1 and 2 hours. The image phase may drift enough to make the PRFS method unusable. An obvious approach is to provide reference materials in the image that do not change phase with temperature in order to monitor and correct the phase drift. Mineral oils provide a suitable such material, since the resonant frequency of oil changes a negligible amount when heated. A difficulty with such a scheme is that the references must be outside of the region where the temperature is to be estimated. Just measuring the phase drift outside of the object does not necessarily determine the drift inside the object. However, it turns out that making simple assumptions such as that the drift is slowly varying over space can provide very good estimates of the drift in places where it is not measured.
Use of such reference materials also illustrates the potential of measuring absolute temperature, since the fat-water frequency difference should be a function of absolute temperature, which may be especially valuable in regions that have mixtures of fat and water such as breast tissue.
Figure 3 shows an axial image obtained across a cylindrical radiofrequency heating device and phantom similar to that in Figure 1. The phantom contains temperature probes for direct temperature measurement; the water coupling bolus contains tubes filled with a mineral oil reference material. The graph in Figure 4a shows the phase temperature change over a period of 65 minutes for a number of ROIs in the phantom and references. During the first 35 minutes, the phantom remained at constant temperature. Radiofrequency heating was then initiated for the next 20 minutes, followed by 10 minutes with the radiofrequency heating stopped. Although the phantom stayed at a constant temperature during the initial unheated period, the calculated temperature drifted downward uniformly for all regions. The drift downward in the gel phantom during heating was masked by the actual temperature rise occurring.

View larger version (41K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 3a. (a) Axial MR image of the phantom in Figure 1 shows the positions of the catheters that contain the direct temperature measurement probes (small arrows). The large arrows indicate the locations of the mineral oil–filled tubes that were placed outside the phantom in the water bolus. The lines in the water bolus are the folds of the plastic membrane that contains the water, and the bright inhomogeneities on the outside are due to the conductive antennae for heating placed on the outside of the outer cylinder. (b) Image obtained by summing all of the temperature changes measured with the PRFS method, which are indicated by an orange-toned intensity scale, shows the circularly symmetric heating pattern at the end of the heating.
|
|

View larger version (113K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 3b. (a) Axial MR image of the phantom in Figure 1 shows the positions of the catheters that contain the direct temperature measurement probes (small arrows). The large arrows indicate the locations of the mineral oil–filled tubes that were placed outside the phantom in the water bolus. The lines in the water bolus are the folds of the plastic membrane that contains the water, and the bright inhomogeneities on the outside are due to the conductive antennae for heating placed on the outside of the outer cylinder. (b) Image obtained by summing all of the temperature changes measured with the PRFS method, which are indicated by an orange-toned intensity scale, shows the circularly symmetric heating pattern at the end of the heating.
|
|

View larger version (18K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 4a. Graphs of the temperature change at the different regions of interest (ROIs) in Figure 3b calculated with the PRFS method over time. During the 65 minutes shown, the phantom remained unheated and at a constant temperature for the first 35 minutes, followed by 20 minutes of heating and then 10 minutes of cooling. (a) Graph shows data that have not been corrected for phase drift, which is clearly demonstrated by the PRFS measurements in the oil references (the constantly decreasing values). (b) Graph shows the data corrected for drift by using spatial interpolation and the oil references. The data are seen to be comparable with the direct temperature measurements (black lines).
|
|

View larger version (19K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 4b. Graphs of the temperature change at the different regions of interest (ROIs) in Figure 3b calculated with the PRFS method over time. During the 65 minutes shown, the phantom remained unheated and at a constant temperature for the first 35 minutes, followed by 20 minutes of heating and then 10 minutes of cooling. (a) Graph shows data that have not been corrected for phase drift, which is clearly demonstrated by the PRFS measurements in the oil references (the constantly decreasing values). (b) Graph shows the data corrected for drift by using spatial interpolation and the oil references. The data are seen to be comparable with the direct temperature measurements (black lines).
|
|
The graph in Figure 4b, which was obtained by using the ROIs to sample the oil reference phase values and a low-order interpolation spatially, shows the temperature measured with the temperature probes (black lines) and the temperature estimated with the PRFS technique, corrected for drift. After correction, the gel was shown to have a constant temperature during the nonheated period and to have temperature changes consistent with the measured temperature changes during the heat periods. There is an apparent overcorrection in the locations near the edge of the phantom during the period just before the heating. This may indicate the degree of validity of the spatial interpolation method or assumptions. The orange-toned image in Figure 3b is a cumulative temperature change map formed at the end of the heating period and showing the typical circular heating pattern in the phantom.
An example of this method for a subject who underwent hyperthermia therapy in a clinical trial is shown in Figure 5 (the data were obtained under an institutional review board–approved protocol with informed consent). The subjects tumor (soft-tissue sarcoma of the calf) is visible in the upper right quadrant of Figure 5a. The graph shows the PRFS temperature estimates in the ROIs defined by the color blocks in Figure 5c and 5d. ROI 1 is in the tumor, near the invasive fiberoptic probe (not visible), ROI 2 is in normal muscle, and ROIs 3 and 4 are in the oil reference materials. The phases and amplitude of the applied radiofrequency heating were adjusted to deliver the highest heating to the tumor region.

View larger version (104K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 5a. Estimates of temperature change (starting at body temperature) during hyperthermia therapy in a patient with a soft-tissue sarcoma of the lower leg. (a) Axial image of the calf and the tumor shows the anatomy. The leg is surrounded by a water bolus and the heating antennae are on the outside of the cylinder, as shown in Figure 1. (b) Graph shows the temperature changes during the 70 minutes of therapy. The different lines correspond to the ROIs, which are shown as color blocks in c and d. 1 = ROI in the tumor near the invasive fiberoptic temperature probe (black line on graph), 2 = ROI in normal muscle, 3 and 4 = ROIs in the oil references. The drift correction has been applied to the graph, as demonstrated by the near-constant values of the oil references. (c, d) Anatomy images obtained at the beginning (c) and near the end (d) of therapy show the regional temperature changes by means of an overlaid orange-toned temperature scale.
|
|

View larger version (19K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 5b. Estimates of temperature change (starting at body temperature) during hyperthermia therapy in a patient with a soft-tissue sarcoma of the lower leg. (a) Axial image of the calf and the tumor shows the anatomy. The leg is surrounded by a water bolus and the heating antennae are on the outside of the cylinder, as shown in Figure 1. (b) Graph shows the temperature changes during the 70 minutes of therapy. The different lines correspond to the ROIs, which are shown as color blocks in c and d. 1 = ROI in the tumor near the invasive fiberoptic temperature probe (black line on graph), 2 = ROI in normal muscle, 3 and 4 = ROIs in the oil references. The drift correction has been applied to the graph, as demonstrated by the near-constant values of the oil references. (c, d) Anatomy images obtained at the beginning (c) and near the end (d) of therapy show the regional temperature changes by means of an overlaid orange-toned temperature scale.
|
|

View larger version (55K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 5c. Estimates of temperature change (starting at body temperature) during hyperthermia therapy in a patient with a soft-tissue sarcoma of the lower leg. (a) Axial image of the calf and the tumor shows the anatomy. The leg is surrounded by a water bolus and the heating antennae are on the outside of the cylinder, as shown in Figure 1. (b) Graph shows the temperature changes during the 70 minutes of therapy. The different lines correspond to the ROIs, which are shown as color blocks in c and d. 1 = ROI in the tumor near the invasive fiberoptic temperature probe (black line on graph), 2 = ROI in normal muscle, 3 and 4 = ROIs in the oil references. The drift correction has been applied to the graph, as demonstrated by the near-constant values of the oil references. (c, d) Anatomy images obtained at the beginning (c) and near the end (d) of therapy show the regional temperature changes by means of an overlaid orange-toned temperature scale.
|
|

View larger version (60K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 5d. Estimates of temperature change (starting at body temperature) during hyperthermia therapy in a patient with a soft-tissue sarcoma of the lower leg. (a) Axial image of the calf and the tumor shows the anatomy. The leg is surrounded by a water bolus and the heating antennae are on the outside of the cylinder, as shown in Figure 1. (b) Graph shows the temperature changes during the 70 minutes of therapy. The different lines correspond to the ROIs, which are shown as color blocks in c and d. 1 = ROI in the tumor near the invasive fiberoptic temperature probe (black line on graph), 2 = ROI in normal muscle, 3 and 4 = ROIs in the oil references. The drift correction has been applied to the graph, as demonstrated by the near-constant values of the oil references. (c, d) Anatomy images obtained at the beginning (c) and near the end (d) of therapy show the regional temperature changes by means of an overlaid orange-toned temperature scale.
|
|
The graph in Figure 5b shows that the invasive measure and ROI 1 (both in the tumor) agree quite well and that the drift correction method results in a constant value for the oil reference materials. The tumor heated up rapidly and showed modulation of the temperature due to blood flow changes alternated with radiofrequency power changes as the body attempted to cool the tumor region. The orange-toned images in Figure 5c and 5d show the temperature changes at the beginning (virtually zero) and after 60 minutes of heating.
The water self-diffusion coefficient, which is called the ADC in heterogeneous tissue, is temperature sensitive. The ADC of water (D, usually expressed in units of millimeters squared per second) in tissue changes with temperature through the following simple relationship:
where Ea is the water activation energy, k is the Boltzmann constant, and T is the absolute temperature. This relationship can be restated in terms of temperature change (
T) as follows:
Thus, the fractional change in temperature can be calculated from the fractional change in the diffusion coefficient. The value of Ea is about 0.21 eV for water, so this method yields about a 2.7% change in D for each degree of temperature change.
The ADC can be measured with most modern MR imaging systems by using echo-planar imaging pulse sequences. The process involves acquiring two images, each with different sensitivity to diffusion. The diffusion sensitivity is related to the strength and timing of the gradients and is denoted as the b value of the pulse sequence. The best signal-to-noise ratio in measuring an ADC value results from obtaining one image with diffusion weighting b = 0 and a second one with b = 1/D. If the first image intensity is denoted S(0) and the second is S(b), then the ADC is calculated from the natural logarithm of the image ratio:
 |
Echo-planar imaging is a fast technique that is challenged by low resolution, low signal-to-noise ratio, and image distortion. In addition, the need to sensitize the images to diffusion leads to long echo times, which further reduce signal-to-noise ratio for tissues that have short T2 relaxation times. Fortunately, many tumors have longer T2 values, so this method can still work in such pathologic tissues. As an example of the noise and linearity of the method, Figure 6 shows the ADC fractional change versus temperature change in canine tissues being heated with radiofrequency energy.
Although the ADC method thus appears to have adequate linearity with temperature, a more complicated analysis shows that it has lower sensitivity to temperature change than the PRFS method. As an example, for an echo-planar diffusion measurement at 1.5 T on a system with a maximum gradient strength of 23 mT/m with repetition time = 6 seconds, 4-mm section thickness, 90-kHz bandwidth, 24-cm field of view, 128 x 128 matrix, and b = 750 sec/mm2, the minimum echo time value is 109 msec, leading to a temperature noise (calculated from the image and diffusion noise) of 1.2°C. To measure temperature to within 1°C, longer imaging times or larger voxels would have to be used.
For the same object, by using the PRFS technique and a gradient-echo acquisition with repetition time = 34 msec, echo time = 20 msec, 16-kHz bandwidth, and a 256 x 256 matrix, the temperature noise was only 0.26°C. The PRFS method in this example thus has twice the image spatial resolution, a shorter 8-second imaging time compared with the longer 12-second imaging time of the ADC measurement (since two images must be acquired), and, most important, much lower temperature noise (by about a factor of 5!). Hence, the PRFS method has become the preferred method of MR temperature measurement.
To have successful hyperthermia therapy, the tumor must be "heatable." If a tumor has large sections of highly vascular tissue, it may be difficult to achieve local heating in these sections of the tumor without overheating poorly or normally vascularized tissue. Hence it is of value, when deciding whether a particular tumor is a good candidate for hyperthermia, to perform contrast material–enhanced or, even better, dynamic contrast-enhanced MR imaging to evaluate tumor vascularity (17,19,43). In addition to allowing the basic morphometry and vascularity of the tumor to be assessed, the vascular parameters (permeability, extravascular space) that can be obtained from an analysis of the dynamic contrast-enhanced data are potentially markers of tissue response to therapy.
As an example of the value of contrast-enhanced MR imaging, Figure 7 shows contrast-enhanced MR images of two different tumors along with color-coded temperature maps obtained at the end of heat therapy, which show the hottest areas in red and the cooler areas in blue (note the somewhat smaller magnification compared with that of the MR images). The top MR image (Fig 7a) shows a strongly enhancing tumor rim consistent with high vascularity and a necrotic core that does not enhance. The thermal map (Fig 7b) shows that only cooler temperatures were achieved by the end of the heating, indicating that the tissue cooling from the vascular flow was enough to reduce the heating—even in the necrotic region.

View larger version (164K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 7a. MR images and temperature maps of two tumors with higher (a, b) and lower (c, d) vascularity. (a, c) Contrast-enhanced MR images show two sarcomas of the leg. (b, d) Corresponding temperature maps show temperature change from the beginning of heating by means of a color scale, which ranges from blue (0°–2°C) to red (>6°C). Note that the image magnification is different between the MR images and the temperature maps.
|
|

View larger version (83K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 7b. MR images and temperature maps of two tumors with higher (a, b) and lower (c, d) vascularity. (a, c) Contrast-enhanced MR images show two sarcomas of the leg. (b, d) Corresponding temperature maps show temperature change from the beginning of heating by means of a color scale, which ranges from blue (0°–2°C) to red (>6°C). Note that the image magnification is different between the MR images and the temperature maps.
|
|

View larger version (141K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 7c. MR images and temperature maps of two tumors with higher (a, b) and lower (c, d) vascularity. (a, c) Contrast-enhanced MR images show two sarcomas of the leg. (b, d) Corresponding temperature maps show temperature change from the beginning of heating by means of a color scale, which ranges from blue (0°–2°C) to red (>6°C). Note that the image magnification is different between the MR images and the temperature maps.
|
|

View larger version (72K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 7d. MR images and temperature maps of two tumors with higher (a, b) and lower (c, d) vascularity. (a, c) Contrast-enhanced MR images show two sarcomas of the leg. (b, d) Corresponding temperature maps show temperature change from the beginning of heating by means of a color scale, which ranges from blue (0°–2°C) to red (>6°C). Note that the image magnification is different between the MR images and the temperature maps.
|
|
The tumor shown in the lower MR image (Fig 7c) does not have such strong enhancement of the rim and has a larger necrotic region. The thermal map (Fig 7d) shows that higher temperatures were achieved in the vascular regions of this tumor compared with those achieved in the more vascular tumor and also shows a hot spot in the necrotic region.
Careful evaluation of subjects as to whether their tumors are heatable has been an important part of conducting successful clinical trials of hyperthermia, since older, less successful trials may have entered subjects who were not suitable for such therapy.
As mentioned in the Introduction, an exciting new area of thermal therapy is to use heat-sensitive liposomes to encapsulate chemotherapeutic agents to allow targeted therapy through hyperthermia. These liposomes can have their properties adjusted to release their contents at 42°C, and thus local hyperthermia will cause most of the drug to be released in the heated region. By encapsulating MR contrast agents in the same liposomes, it may be possible to make MR images of the regions where the agent is released. By using quantitative MR imaging, the concentration of the MR contrast agent can be calculated, which may allow accurate assessment of the concentration of the drug deposited throughout the tumor (26,44–48).
Imaging in hyperthermia therapy is performed primarily to estimate and control temperature. MR imaging has unique parameter dependences that make this possible by detection of proton resonant frequency changes or diffusion coefficient changes. In addition, MR imaging can be used to assess vascular parameters that not only allow suitable selection of subjects for therapy but may also allow detection of response to therapy. Finally, as the use of thermally sensitive liposomes matures for delivery of chemotherapeutic agents, imaging may allow determination of local drug dose.
The authors gratefully acknowledge valuable discussions with and images from Mark Dewhirst, DVM, PhD, Paul Stauffer, MSEE, CCE, Oana Craciunescu, PhD, and Zeljko Vujaskovic, MD, PhD.
Abbreviations: ADC = apparent diffusion coefficient, PRFS = proton resonant frequency shift, ROI = region of interest
- Dewhirst M, Jones E, Samulski RJ, Vujaskovic Z, Li C, Prosnitz L. Hyperthermia. In: Kufe DW, et al, eds. Cancer medicine 6. Hamilton, Ontario: Decker, 2003; 623–636.
- Falk MH, Issels RD. Hyperthermia in oncology. Int J Hyperthermia 2001;17(1):1–18.[Medline]
- Wust P, Hildebrandt B, Sreenivasa G, et al. Hyperthermia in combined treatment of cancer. Lancet Oncol 2002;3(8):487–497.[CrossRef][Medline]
- Vernon CC, Hand JW, Field SB, et al. Radiotherapy with or without hyperthermia in the treatment of superficial localized breast cancer: results from five randomized controlled trials. International Collaborative Hyperthermia Group. Int J Radiat Oncol Biol Phys 1996;35(4):731–744.[CrossRef][Medline]
- Overgaard J, Gonzalez Gonzalez D, Hulshof MC, et al. Hyperthermia as an adjuvant to radiation therapy of recurrent or metastatic malignant melanoma: a multicentre randomized trial by the European Society for Hyperthermic Oncology. Int J Hyperthermia 1996;12(1):3–20.[Medline]
- Valdagni R, Amichetti M. Report of long-term follow-up in a randomized trial comparing radiation therapy and radiation therapy plus hyperthermia to metastatic lymph nodes in stage IV head and neck patients. Int J Radiat Oncol Biol Phys 1994;28(1):163–169.[Medline]
- Kitamura K, Kuwano H, Watanabe M, et al. Prospective randomized study of hyperthermia combined with chemoradiotherapy for esophageal carcinoma. J Surg Oncol 1995;60(1):55–58.[Medline]
- van der Zee J, González González D, van Rhoon GC, van Dijk JD, van Putten WL, Hart AA. Comparison of radiotherapy alone with radiotherapy plus hyperthermia in locally advanced pelvic tumours: a prospective, randomised, multicentre trial. Dutch Deep Hyperthermia Group. Lancet 2000;355(9210):1119–1125.[CrossRef][Medline]
- Sneed PK, Stauffer PR, McDermott MW, et al. Survival benefit of hyperthermia in a prospective randomized trial of brachytherapy boost +/– hyperthermia for glioblastoma multiforme. Int J Radiat Oncol Biol Phys 1998;40(2):287–295.[CrossRef][Medline]
- Oleson JR, Samulski TV, Leopold KA, et al. Sensitivity of hyperthermia trial outcomes to temperature and time: implications for thermal goals of treatment. Int J Radiat Oncol Biol Phys 1993; 25(2):289–297.[Medline]
- Kapp DS, Cox RS. Thermal treatment parameters are most predictive of outcome in patients with single tumor nodules per treatment field in recurrent adenocarcinoma of the breast. Int J Radiat Oncol Biol Phys 1995;33(4):887–899.[CrossRef][Medline]
- Hand JW, Machin D, Vernon CC, Whaley JB. Analysis of thermal parameters obtained during phase III trials of hyperthermia as an adjunct to radiotherapy in the treatment of breast carcinoma. Int J Hyperthermia 1997;13(4):343–364.[Medline]
- Sherar M, Liu FF, Pintilie M, et al. Relationship between thermal dose and outcome in thermoradiotherapy treatments for superficial recurrences of breast cancer: data from a phase III trial. Int J Radiat Oncol Biol Phys 1997;39(2):371–380.[CrossRef][Medline]
- Seegenschmiedt MH, Martus P, Fietkau R, Iro H, Brady LW, Sauer R. Multivariate analysis of prognostic parameters using interstitial thermoradiotherapy (IHT-IRT): tumor and treatment variables predict outcome. Int J Radiat Oncol Biol Phys 1994;29(5):1049–1063.[Medline]
- Clegg ST, Das SK, Zhang Y, Macfall J, Fullar E, Samulski TV. Verification of a hyperthermia model method using MR thermometry. Int J Hyperthermia 1995;11(3):409–424.[Medline]
- Craciunescu OI, Samulski TV, MacFall JR, Clegg ST. Perturbations in hyperthermia temperature distributions associated with counter-current flow: numerical simulations and empirical verification. IEEE Trans Biomed Eng 2000;47(4):435–443.[CrossRef][Medline]
- Craciunescu OI, Das SK, McCauley RL, MacFall JR, Samulski TV. 3D numerical reconstruction of the hyperthermia induced temperature distribution in human sarcomas using DE-MRI measured tissue perfusion: validation against non-invasive MR temperature measurements. Int J Hyperthermia 2001;17(3):221–239.[CrossRef][Medline]
- MacFall JR, Prescott DM, Charles HC, Samulski TV. 1H MRI phase thermometry in vivo in canine brain, muscle, and tumor tissue. Med Phys 1996; 23(10):1775–1782.[CrossRef][Medline]
- Craciunescu OI, Raaymakers BW, Kotte AN, Das SK, Samulski TV, Lagendijk JJ. Discretizing large traceable vessels and using DE-MRI perfusion maps yields numerical temperature contours that match the MR noninvasive measurements. Med Phys 2001;28(11):2289–2296.[CrossRef][Medline]
- Samulski TV, Clegg ST, Das S, MacFall J, Prescott DM. Application of new technology in clinical hyperthermia. Int J Hyperthermia 1994;10(3): 389–394.[Medline]
- Carter DL, MacFall JR, Clegg ST, et al. Magnetic resonance thermometry during hyperthermia for human high-grade sarcoma. Int J Radiat Oncol Biol Phys 1998;40(4):815–822.[CrossRef][Medline]
- Jones E, Thrall D, Dewhirst MW, Vujaskovic Z. Prospective thermal dosimetry: the key to hyperthermias future. Int J Hyperthermia 2006;22(3): 247–253.[CrossRef][Medline]
- Thrall DE, LaRue SM, Yu D, et al. Thermal dose is related to duration of local control in canine sarcomas treated with thermoradiotherapy. Clin Cancer Res 2005;11(14):5206–5214.[Abstract/Free Full Text]
- Jones EL, Oleson JR, Prosnitz LR, et al. Randomized trial of hyperthermia and radiation for superficial tumors. J Clin Oncol 2005;23(13):3079–3085.[Abstract/Free Full Text]
- Edelstein-Keshet L, Dewhirst MW, Oleson JR, Samulski TV. Characterization of tumour temperature distributions in hyperthermia based on assumed mathematical forms. Int J Hyperthermia 1989;5(6):757–777.[Medline]
- Chen Q, Tong S, Dewhirst MW, Yuan F. Targeting tumor microvessels using doxorubicin encapsulated in a novel thermosensitive liposome. Mol Cancer Ther 2004;3(10):1311–1317.[Abstract/Free Full Text]
- Hutchinson E, Dahleh M, Hynynen K. The feasibility of MRI feedback control for intracavitary phased array hyperthermia treatments. Int J Hyperthermia 1998;14(1):39–56.[Medline]
- Cline HE, Hynynen K, Schneider E, et al. Simultaneous magnetic resonance phase and magnitude temperature maps in muscle. Magn Reson Med 1996;35(3):309–315.[Medline]
- De Poorter J. Noninvasive MRI thermometry with the proton resonance frequency method: study of susceptibility effects. Magn Reson Med 1995; 34(3):359–367.[Medline]
- Ishihara Y, Calderon A, Watanabe H, et al. A precise and fast temperature mapping using water proton chemical shift. Magn Reson Med 1995; 34(6):814–823.[Medline]
- Parker DL. Applications of NMR imaging in hyperthermia: an evaluation of the potential for localized tissue heating and noninvasive temperature monitoring. IEEE Trans Biomed Eng 1984;31(1): 161–167.[CrossRef][Medline]
- Samulski TV, MacFall J, Zhang Y, Grant W, Charles C. Non-invasive thermometry using magnetic resonance diffusion imaging: potential for application in hyperthermic oncology. Int J Hyperthermia 1992;8(6):819–829.[Medline]
- Cetas TC. Will thermometric tomography become practical for hyperthermia treatment monitoring? Cancer Res 1984;44(10 suppl):4805s–4808s.[Abstract/Free Full Text]
- Delannoy J, LeBihan D, Hoult DI, Levin RL. Hyperthermia system combined with a magnetic resonance imaging unit. Med Phys 1990;17(5): 855–860.[CrossRef][Medline]
- Ishihara Y, Watanabe H, Okamoto K, Kanamatsu T, Tsukada Y. Temperature monitoring of internal body heating induced by decoupling pulses in animal (13)C-MRS experiments. Magn Reson Med 2000;43(6):796–803.[CrossRef][Medline]
- Jolesz FA, Bleier AR, Jakab P, Ruenzel PW, Huttl K, Jako GJ. MR imaging of laser-tissue interactions. Radiology 1988;168(1):249–253.[Abstract/Free Full Text]
- Samulski TV, Fessenden P, Lee ER, Kapp DS, Tanabe E, McEuen A. Spiral microstrip hyperthermia applicators: technical design and clinical performance. Int J Radiat Oncol Biol Phys 1990; 18(1):233–242.[Medline]
- Zhang Y, Joines WT, Oleson JR. Prediction of heating patterns of a microwave interstitial antenna array at various insertion depths. Int J Hyperthermia 1991;7(1):197–207.[Medline]
- Prescott DM, Charles HC, Sostman HD, et al. Therapy monitoring in human and canine soft tissue sarcomas using magnetic resonance imaging and spectroscopy. Int J Radiat Oncol Biol Phys 1994;28(2):415–423.[Medline]
- Das SK, Clegg ST, Samulski TV. Computational techniques for fast hyperthermia temperature optimization. Med Phys 1999;26(2):319–328.[CrossRef][Medline]
- Gellermann J, Wlodarczyk W, Ganter H, et al. A practical approach to thermography in a hyperthermia/magnetic resonance hybrid system: validation in a heterogeneous phantom. Int J Radiat Oncol Biol Phys 2005;61(1):267–277.[CrossRef][Medline]
- Wust P, Gellermann J, Seebass M, et al. Part-body hyperthermia with a radiofrequency multiantenna applicator under online control in a 1.5 T MR-tomograph [in German]. Rofo 2004;176(3):363–374.[Medline]
- Craciunescu OI, Das SK, Poulson JM, Samulski TV. Three-dimensional tumor perfusion reconstruction using fractal interpolation functions. IEEE Trans Biomed Eng 2001;48(4):462–473.[CrossRef][Medline]
- Viglianti BL, Ponce AM, Michelich CR, et al. Chemodosimetry of in vivo tumor liposomal drug concentration using MRI. Magn Reson Med 2006;56(5):1011–1018.[CrossRef][Medline]
- Ponce AM, Vujaskovic Z, Yuan F, Needham D, Dewhirst MW. Hyperthermia mediated liposomal drug delivery. Int J Hyperthermia 2006;22(3): 205–213.[CrossRef][Medline]
- Hauck ML, LaRue SM, Petros WP, et al. Phase I trial of doxorubicin-containing low temperature sensitive liposomes in spontaneous canine tumors. Clin Cancer Res 2006;12(13):4004–4010.[Abstract/Free Full Text]
- Viglianti BL, Abraham SA, Michelich CR, et al. In vivo monitoring of tissue pharmacokinetics of liposome/drug using MRI: illustration of targeted delivery. Magn Reson Med 2004;51(6):1153–1162.[CrossRef][Medline]
- Needham D, Dewhirst MW. The development and testing of a new temperature-sensitive drug delivery system for the treatment of solid tumors. Adv Drug Deliv Rev 2001;53(3):285–305.[CrossRef][Medline]