Simultaneous magnetic resonance phase and magnitude temperature maps in muscle

Simultaneous measurement of temperature and velocity maps by inversion recovery tagging method

A new method, called the inversion recovery (IR) tagging method, for simultaneous measurement of temperature and velocity maps of flowing fluid has been developed. The present method employs a set of tagging pulses which acts as an inversion pulse of the conventional IR method, based on the temperature dependence of the spin-lattice relaxation of water proton in a fluid, and has the advantage of being able to compensate the reduction of the NMR signal intensity due to flow motion and to reduce the total time to measure these maps. First, the accuracy of the temperature measurement of stagnant doped water in a differentially heated cell using the conventional IR method, as the basic sequence of the IR tagging method, has been evaluated. The accuracy was within 10% of the temperature difference DeltaT = 17.2 degrees C and the measurable temperature resolution was within +/-0.5 degrees C. Then temperature and velocity maps of the flowing doped-water through a cooled pipe were measured simultaneously by the IR tagging method, and the accuracy of temperature measurement was evaluated. The accuracy obtained using the present method was within 15% of the temperature difference DeltaT = 15 degrees C.

Temperature mapping using the water proton chemical shift: self-referenced method with echo-planar spectroscopic imaging

An echo-planar spectroscopic imaging method of temperature mapping is proposed. This method is sufficiently faster than the so-called 3D magnetic resonance spectroscopic imaging (3D-MRSI) method and does not require image subtractions, unlike the conventional phase mapping method when an internal reference signal is detectable. The water proton chemical shift measured by using the tissue lipid as an internal reference clearly visualized the temperature change in a porcine liver sample in vitro. It was also demonstrated that the internally referenced echo-planar spectroscopic imaging method could markedly reduce a temperature error caused by a simple, translational motion between scans compared with the phase-mapping method.

MRI phase mapping of temperature distributions induced in food by microwave heating

A range of temperature-sensitive MRI parameters of water (T2, T1, diffusion coefficient, and chemical shift) were evaluated to map in three dimensions the non-uniform temperature distributions induced by microwave heating in both model and real food systems. Phase mapping was found to be the most robust method, and evaluations of possible experimental errors were based on semi-quantitative studies of homogeneous and heterogeneous systems. The MRI protocol provides complementary phase and magnitude data, which are related to the sample temperature and structural heterogeneity, respectively. Used together, they relate the temperature changes to the differential thermal properties of the various components within a heterogeneous sample. The potential applications of this technique to microwave and other forms of heating is discussed.

Reliability of water proton chemical shift temperature calibration for focused ultrasound ablation therapy

Our purpose in this work was to assess the reliability of the calibration coefficient for magnetic resonance water proton chemical shift temperature mapping. Over a six month period, the calibration coefficient was measured 15 times in several different phantoms. A highly linear relationship between water proton chemical shift and temperature change was found. The average temperature calibration coefficient determined from all studies was 0.009+/-0.001 ppm/degrees C. Four of the 15 studies were conducted on the same day using the same phantom. The average temperature calibration coefficient of these four studies was 0.0096+/-0.0001 ppm/degrees C.

Thermal dosimetry of a focused ultrasound beam in vivo by magnetic resonance imaging

Magnetic resonance imaging (MRI) thermometry has been utilized for in vivo evaluation of thermal exposure induced by a focused ultrasound beam. A simulation study of the focused ultrasound beam was conducted to select imaging parameters for reducing the error due to the spatial and temporal averaging of MRI. Temperature imaging based on the proton resonance frequency shift was utilized to obtain the temperature distribution during sonication in the skeletal muscle of eight rabbits. MRI-derived temperature information was then used to calculate the thermal dose distribution induced by the sonication and to estimate the coagulated tissue volume. The tissue changes were also evaluated directly by taking the T2-weighted and the contrast agent enhanced T1-weighted MR images. Errors in the temperature and thermal dose measurements were found to be minimal using the following parameters: slice thickness = 3 mm, voxel dimension = 0.6 mm, and scan time per image = 3.4 s. The estimated dimensions of the coagulated tissue volume were in good agreement with the tissue damages seen on the contrast agent enhanced T1-weighted images. The tissue damage seen on the histology was closely matched to the ones seen on the T2-weighted images. This study showed that MRI thermometry has significant potential for both monitoring the thermal exposure and evaluating the tissue damage. This would allow real-time control of the sonication parameters to optimize clinical treatments.

Fast lipid-suppressed MR temperature mapping with echo-shifted gradient-echo imaging and spectral-spatial excitation

The water proton resonance frequency (PRF) is temperature dependent and can thus be used for magnetic resonance (MR) thermometry. Since lipid proton resonance frequencies do not depend on temperature, fat suppression is essential for PRF-based temperature mapping. The efficacy of echo-shifted (TE > TR) gradient-echo imaging with spectral-spatial excitation is demonstrated, resulting in accurate and rapid, lipid-suppressed, MR thermometry. The method was validated on phantoms, fatty duck liver, and rat thigh, demonstrating improvements in both the speed and precision of temperature mapping. Heating of a rat thigh with focused ultrasound was monitored in vivo with an accuracy of 0.37 degree C and a time resolution of 438 msec.

Ultrasound technology for hyperthermia

Hyperthermia (HT) is used in the clinical management of cancer and benign disease. Numerous biological and clinical investigations have demonstrated that HT in the 41-45 degrees C range can significantly enhance clinical responses to radiation therapy, and has potential for enhancing other therapies, such as chemotherapy, immunotherapy and gene therapy. Furthermore, high-temperature hyperthermia (greater than 50 degrees C) alone is being used for selective tissue destruction as an alternative to conventional invasive surgery. The degree of thermal enhancement of these therapies is strongly dependent on the ability to localize and maintain therapeutic temperature elevations. Due to the often heterogeneous and dynamic properties of tissues, most notably blood perfusion and the presence of thermally significant blood vessels, therapeutic temperature elevations are difficult to spatially and temporally control during these forms of HT therapy. However, ultrasound technology has significant advantages that allow for a higher degree of spatial and dynamic control of the heating compared to other commonly utilized heating modalities. These advantages include a favorable range of energy penetration characteristics in soft tissue and the ability to shape the energy deposition patterns. Thus, heating systems have been developed for interstitial, intracavitary, or external approaches that utilize properties such as multiple transducer arrays, phased arrays, focused beams, mechanical and/or electrical scanning, dynamic frequency control and transducers of various shapes and sizes. This article provides a general review of a selection of ultrasound hyperthermia systems that are either in clinical use or currently under development, that utilize these advantages as a means to better localize and control HT for the aforementioned therapies.

The feasibility of elastographic visualization of HIFU-induced thermal lesions in soft tissues. Image-guided high-intensity focused ultrasound

The potential for visualizing high-intensity focused ultrasound (HIFU)-induced thermal lesions in biological soft tissues in vitro using elastography was investigated. Thermal lesions were created in rabbit paraspinal skeletal muscle in vivo. The rabbits were sacrificed 60 h following the treatment and lesioned tissues were excised. The tissues were cast in a block of clear gel and elastographic images of the lesions were acquired. Gross pathology of the tissue samples confirmed the characteristics of the lesions.

Comparison of four magnetic resonance methods for mapping small temperature changes

Non-invasive detection of small temperature changes (< 1 degree C) is pivotal to the further advance of regional hyperthermia as a treatment modality for deep-seated tumours. Magnetic resonance (MR) thermography methods are considered to be a promising approach. Four methods exploiting temperature-dependent parameters were evaluated in phantom experiments. The investigated temperature indicators were spin-lattice relaxation time T1, diffusion coefficient D, shift of water proton resonance frequency (water PRF) and resonance frequency shift of the methoxy group of the praseodymium complex (Pr probe). The respective pulse sequences employed to detect temperature-dependent signal changes were the multiple readout single inversion recovery (T One by Multiple Read Out Pulses; TOMROP), the pulsed gradient spin echo (PGSE), the fast low-angle shot (FLASH) with phase difference reconstruction, and the classical chemical shift imaging (CSI). Applying these sequences, experiments were performed in two separate and consecutive steps. In the first step, calibration curves were recorded for all four methods. In the second step, applying these calibration data, maps of temperature changes were generated and verified. With the equal total acquisition time of approximately 4 min for all four methods, the uncertainties of temperature changes derived from the calibration curves were less than 1 degree C (Pr probe 0.11 degrees C, water PRF 0.22 degrees C, D 0.48 degrees C and T1 0.93 degrees C). The corresponding maps of temperature changes exhibited slightly higher errors but still in the range or less than 1 degree C (0.97 degrees C, 0.41 degrees C, 0.70 degrees C, 1.06 degrees C respectively). The calibration results indicate the Pr probe method to be most sensitive and accurate. However, this advantage could only be partially transferred to the thermographic maps because of the coarse 16 x 16 matrix of the classical CSI sequence. Therefore, at present the water PRF method appears to be most suitable for MR monitoring of small temperature changes during hyperthermia treatment.

Experimental evaluation of two simple thermal models using transient temperature analysis

Thermal models are used to predict temperature distributions of heated tissues during thermal therapies. Recent interest in short duration high temperature therapeutic procedures necessitates the accurate modelling of transient temperature profiles in heated tissues. Blood flow plays an important role in tissue heat transfer and the resultant temperature distribution. This work examines the transient predictions of two simple mathematical models of heat transfer by blood flow (the bioheat transfer equation model and the effective thermal conductivity equation model) and compares their predictions to measured transient temperature data. Large differences between the two models are predicted in the tissue temperature distribution as a function of blood flow for a short heat pulse. In the experiments a hot water needle, approximately 30 degrees C above ambient, delivered a 20 s heating pulse to an excised fixed porcine kidney that was used as a flow model. Temperature profiles of a thermocouple that primarily traversed the kidney cortex were examined. Kidney locations with large vessels were avoided in the temperature profile analysis by examination of the vessel geometry using high resolution computed tomography angiography and the detection of the characteristic large vessel localized cooling or heating patterns in steady-state temperature profiles. It was found that for regions without large vessels, predictions of the Pennes bioheat transfer equation were in much better agreement with the experimental data when compared to predictions of the scalar effective thermal conductivity equation model. For example, at a location r approximately 2 mm away from the source, the measured delay time was 10.6 +/- 0.5 s compared to predictions of 9.4 s and 5.4 s of the BHTE and ETCE models, respectively. However, for the majority of measured locations, localized cooling and heating effects were detected close to large vessels when the kidney was perfused. Finally, it is shown that increasing flow in regions without large vessels minimally perturbs temperature profiles for short exposure times; regions with large vessels still have a significant effect.

MRI thermometry in phantoms by use of the proton resonance frequency shift method: application to interstitial laser thermotherapy

In this work the temperature dependence of the proton resonance frequency was assessed in agarose gel with a high melting temperature (95 degrees C) and in porcine liver in vitro at temperatures relevant to thermotherapy (25-80 degrees C). Furthermore, an optically tissue-like agarose gel phantom was developed and evaluated for use in MRI. The phantom was used to visualize temperature distributions from a diffusing laser fibre by means of the proton resonance frequency shift method. An approximately linear relationship (0.0085 ppm degrees C(-1)) between proton resonance frequency shift and temperature change was found for agarose gel, whereas deviations from a linear relationship were observed for porcine liver. The optically tissue-like agarose gel allowed reliable MRI temperature monitoring, and the MR relaxation times (T1 and T2) and the optical properties were found to be independently alterable. Temperature distributions around a diffusing laser fibre, during irradiation and subsequent cooling, were assessed with high spatial resolution (voxel size = 4.3 mm3) and with random uncertainties ranging from 0.3 degrees C to 1.4 degrees C (1 SD) with a 40 s scan time.

Monitoring and visualization techniques for MR-guided laser ablations in an open MR system

Our purpose was to develop temperature-sensitive MR sequences and image-processing techniques to assess their potential of monitoring interstitial laser therapy (ILT) in brain tumors (n = 3) and liver tumors (n = 7). ILT lasted 2 to 26 minutes, whereas images from T1-weighted fast-spin-echo (FSE) or spoiled gradient-recalled (SPGR) sequences were acquired within 5 to 13 seconds. Pixel subtraction and visualization of T1-weighted images or optical flow computation was done within less than 110 msec. Alternating phase-mapping of real and imaginary components of SPGR sequences was performed within 220 msec. Pixel subtraction of T1-weighted images identified thermal changes in liver and brain tumors but could not evaluate the temperature values as chemical shift-based imaging, which was, however, more affected by susceptibility effects and motion. Optical flow computation displayed the predicted course of thermal changes and revealed that the rate of heat deposition can be anisotropic, which may be related to heterogeneous tumor structure and/or vascularization.

NMR temperature measurements using a paramagnetic lanthanide complex

NMR thermometry has previously suffered from poor thermal resolution owing to the relatively weak dependence of chemical shift on temperature in diamagnetic molecules. In contrast, the shifts of nuclear spins near a paramagnetic center exhibit strong temperature dependencies. The chemical shifts of the thulium 1,4,7, 10-tetraazacyclododecane-1,4,7,10-tetrakis(methylene phosphonate) complex (TmDOTP5-) have been studied as a function of temperature, pH, and Ca2+ concentration over ranges which may be encountered in vivo. The results demonstrate that the 1H and 31P shifts in TmDOTP5- are highly sensitive to temperature and may be used for NMR thermometry with excellent accuracy and resolution. A new technique is also described which permits simultaneous measurements of temperature and pH changes from the shifts of multiple TmDOTP5- spectral lines.

Tissue temperature monitoring with multiple gradient-echo imaging sequences

The inherent sensitivity of multiple gradient-echo sequences to the chemical shift is exploited to rapidly map muscle water frequency shifts caused by ultrasonic heating. The use of multiple echoes is shown to offer several advantages over single gradient-echo approaches previously proposed for temperature measurement. An increase in the effective bandwidth significantly reduces aliasing problems observed with single gradient-echo methods in high temperature applications. Of greater significance is the improved immunity to intrascan motion found for multi-echo versus single echo gradient methods, making the former more attractive for clinical applications. Finally, a sensitivity to the presence of multiple spectral components unavailable with single gradient-echo methods is obtained.

Time and temperature dependence of MR parameters during thermal coagulation of ex vivo rabbit muscle

Detailed measurements of the T1-weighted, T2-weighted, and MT-weighted signal were performed for ex vivo muscle samples heated to various temperatures for different times. Consistent, monotonic increases in signal intensity were observed with progressive thermal coagulation, corresponding to an increase in T2 relaxation time and an increase in MT-weighted signal for temperatures above 60 degrees C. The relationship for T1 relaxation was more complex, showing a decrease in T1 relaxation from 40 to 60 degrees C and an increase above 60 degrees C. These techniques provide a more direct measure of tissue thermal coagulation than that provided by MR thermometry and suggest MR imaging strategies for the optimization and monitoring of thermal coagulation therapy protocols that create thermal damage in target tissues.

Temperature monitoring of interstitial thermal tissue coagulation using MR phase images

The temperature-dependent water proton frequency shift was investigated for temperature monitoring of interstitial thermal coagulation. A procedure for on-line temperature calculation was developed, and errors due to temperature-dependent susceptibility were investigated by finite element analysis and reference measurements. The temperature coefficient of magnetic susceptibility and proton chemical shift were determined for brain tissue and other substances. With the proposed procedure, the location of isotherms could be well visualized during laser-induced interstitial coagulation in vitro and in vivo. Systematic errors caused by magnetic susceptibility changes with temperature depend strongly on the characteristics of the heat source and can exceed susceptibility effects caused by physiologic tissue changes. For the laser applicators discussed here, however, a first order compensation for this effect was found to be satisfactory, because it reduces the absolute error to the range of +/- 1 degrees C. The proposed method represents a very promising approach for monitoring of the interstitial thermal coagulation.

Three-dimensional monitoring of small temperature changes for therapeutic hyperthermia using MR

Radiofrequency hyperthermia of deep-seated pelvic tumors requires noninvasive monitoring of temperature distributions in patients. Methods of MR thermography were reported to be a promising tool in solving this problem. However, to be truly useful for monitoring hyperthermia treatments, MR thermography should be able to cover the entire pelvis in acquisition times no longer than for a breath-hold (< or = 15 seconds) and to resolve small temperature differences (< 1 degrees C). Three methods exploiting the temperature dependence of spin-lattice relaxation time (T1), of self-diffusion coefficient (D), and of chemical shift of proton resonance frequency (PRF) were applied in phantom experiments; the pulse sequences were the T1-weighted gradient echo, the pulsed diffusion gradient spin echo made faster through the keyhole technique, and the gradient echo with the phase reconstruction, respectively. The high planar resolution was compromised, and instead, coarse and more isotropic voxels were used. Experiments were performed in two consecutive steps, thus imitating a possible scenario for monitoring hyperthermia. In the first step, calibration curves were recorded, which were then used in the second step to obtain maps of temperature changes. The results show a clear superiority of the PRF method, followed by the D and the T1 methods. The uncertainty of temperature changes predicted both from calibration curves and from maps was less than 1 degrees C only with the PRF and the D-based methods.

Temperature monitoring of ultrasonically heated muscle with RARE chemical shift imaging

The ability to monitor tissue temperature in ultrasonically heated rabbit muscle is demonstrated using a chemical shift imaging approach based on the rapid acquisition with relaxation enhancement (RARE) fast imaging method [Hennig et al., Magn. Reson. Med. 3, 823-833 (1986)] applied in a line scan format. A three echo sequence with a 16 Hz spectral resolution with 64 ms echo readouts and 78 ms echo spacings is shown capable of measuring relevantly small water frequency shifts in phantoms. Applied to the in vivo model of ultrasonically heated rabbit muscle, water resonance frequencies at the ultrasonic focal point were found to be linearly related to temperature with a slope of -0.007 +/- 0.001 ppm/degree C (N = 6 studies). Measurements of the frequency shift in unheated tissue located away from the ultrasonically heated tissue varied by approximately 0.011 ppm over the course of the experiments, leading to an estimated temperature accuracy of +/- 1.6 degrees C in vivo.

Phase imaging on a .2-T MR scanner: application to temperature monitoring during ablation procedures

Proton phase shift imaging methods with keyholing were developed to rapidly monitor temperature during MR-guided radiofrequency (RF) interventional procedures on a .2-T open configuration scanner. Temperature calibration was performed on thermally controlled gel phantom and ex vivo bovine liver samples. Keyholing methods were implemented for rapid imaging and tested both in simulation experiments and in the gel phantom. Phase drifts from extraneous sources were monitored and compensated for using reference phantoms. Sequence parameters TE, TR, and flip angle (FA) were optimized for maximum temperature sensitivity and minimum noise. Reduction of phase noise from coupling of the magnetic field to external perturbations using navigator-echo-based correction schemes were also investigated. The extraneous phase drifts from the magnet could be minimized by keeping the electromagnet on continuously. Navigator echo corrected keyholed FLASH sequences (TE = 30 msec, TR = 60 msec, FA = 40 degrees, 64 x 128 matrix) were used to monitor the RF lesioning process in gel phantoms yielding images every 4 seconds with a temperature sensitivity of .015 ppm/degree C. RF ablation in the bovine tissue was monitored using navigator-echo-corrected keyholed fast low angle shot (FLASH) sequences (TE = 30 msec, TR = 100 msec, FA = 40 degrees, 128 x 256 matrix) with a temporal resolution of 13 seconds and a temperature sensitivity of .007 ppm/degree C. The results indicate that monitoring of an RF ablation procedure by mapping temperature with sufficient temporal resolution is possible using phase images of FLASH sequences on a .2-T open scanner.

TmDOTP5-: a substance for NMR temperature measurements in vivo

The chemical shifts of 31P and 1H in thulium 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetrakis(methylene phosphonate) (TmDOTP5-) are approximately two orders of magnitude more sensitive to temperature than are water proton and 19F shifts. In the physiologically relevant pH range, the 31P and 1H chemical shifts of TmDOTP5- are linear functions of temperature between 25 and 47 degrees C. The results indicate that using TmDOTP5- can provide measurements of temperature in vivo that are significantly more accurate than methods based on water and fluorocarbon chemical shifts.

Noninvasive arterial occlusion using MRI-guided focused ultrasound

The purpose of this work was to test the hypothesis that reproducible and sustainable arterial occlusion can be induced by focused ultrasound energy deposition noninvasively within deep tissue. An MRI-compatible focused ultrasound transducer was used to sonicate a branch of the renal artery (diameter about 0.6 mm) in vivo (nine rabbits). An intravenous MRI contrast agent bolus was injected about 30 min and up to 7 days after the sonication. After follow-up, in vitro magnification x-ray angiograms were obtained and the kidneys were fixed in formaldehyde for histologic study. The ultrasound pulses resulted in complete cessation of blood flow, as shown by the gradient echo images. In seven of the nine rabbits, a wedge-shaped unenhanced area was seen at the part of the kidney that was perfused by the vessel after the contrast agent injection. This area extended laterally (outside of the sonicated volume) to the cortical surface of the kidney. The x-ray angiograms showed that the artery was completely occluded. Postmortem histologic evaluation showed an infarcted tissue volume corresponding to the wedge shape seen in the images. This study showed that appropriately focused ultrasound can be used to close arteries noninvasively. This finding has significant clinical potential.