MR temperature mapping of focused ultrasound surgery

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.

Real-time interactive MRI on a conventional scanner

A real-time interactive MRI system capable of localizing coronary arteries and imaging arrhythmic hearts in real-time is described. Non-2DFT acquisition strategies such as spiral-interleaf, spiral-ring, and circular echo-planar imaging provide short scan times on a conventional scanner. Real-time gridding reconstruction at 8-20 images/s is achieved by distributing the reconstruction on general-purpose UNIX workstations. An X-windows application provides interactive control. A six-interleaf spiral sequence is used for cardiac imaging and can acquire six images/s. A sliding window reconstruction achieves display rates of 16-20 images/s. This allows cardiac images to be acquired in real-time, with minimal motion and flow artifacts, and without breath holding or cardiac gating. Abdominal images are acquired at over 2.5 images/s with spiral-ring or circular echo-planar sequences. Reconstruction rates are 8-10 images/s. Rapid localization in the abdomen is demonstrated with the spiral-ring acquisition, whereas peristaltic motion in the small bowel is well visualized using the circular echo-planar sequence.

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.

Determination of laser-induced temperature distributions using echo-shifted TurboFLASH

An echo-shifted TurboFLASH sequence implemented on a clinical whole body MR scanner was used to determine thermal changes in tissue. With this snapshot-like data acquisition, temperature-related phase shifts were measured with a temporal resolution of 1.3 s. For different types of tissue (postmortem porcine brain, liver, and muscle) the temperature coefficients of the proton chemical shift were recorded during uniform heating of the specimen in a water bath. The specific temperature-dependent frequency shifts appeared similar to the proton chemical shift of free water (-0.01 ppm/degrees C). With this method, laser-induced ablation in postmortem porcine brain was monitored by temperature mapping. Comparison of the induced temperature profiles measured with NiCrNi-thermocouples with the MR calculated profiles demonstrated excellent temperature sensitivity and accuracy for this method of MR thermometry, with a maximum deviation of the determined temperatures of only 1.8 degrees C. This investigation was designed as a feasibility study for this rapid version of the phase mapping method, and no in vivo studies were performed.

Influence of ablated tissue on the formation of high-intensity focused ultrasound lesions

In order to ablate tumours using high-intensity focused ultrasound (HIFU) it is necessary to irradiate the tumour with a confluent array of single ultrasound exposures. We have identified a phenomenon that we term lesion-to-lesion interaction, which occurs when the spatial separation of individual exposures is such that an existing lesions appears to affect the formation of a subsequent lesion. This article investigates the implications of this phenomenon for strategies to ablate large tissue volumes in the treatment of hepatic metastases. Experiments on pig and rat livers have been carried out using a focused ultrasound system with a frequency of 1.7 MHz, an in situ spatially averaged focal intensity (ISAL) of 133-658 W cm-2 (ISP of 239-1185 W cm-2) and an exposure duration of 5-15 s. The results show that there is interaction between lesions that spatial exposure separations that depend on the intensities and exposure durations used. As a result, either subsequent lesions form closer to the ultrasound source (if the focal peak of the ultrasound beam is placed deep inside the liver tissue) or their length is reduced (if the focal peak is near the liver surface). An explanation is suggested for this effect and a strategy for its avoidance during in vivo HIFU treatment is discussed.

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.

2D locally focused MRI: applications to dynamic and spectroscopic imaging

Conventional magnetic resonance images have uniform spatial resolution across the entire field of view. A method of creating MR images with user-specified spatial resolution along one dimension of the field of view was described recently by the authors. This paper presents the 2D generalization of this technique, which allows the user to specify arbitrary spatial resolution in arbitrary 2D regions. These images are reconstructed from signals that sparsely sample the k-space representation of the image. Therefore, locally focused images can be acquired in less time than that required by Fourier imaging with uniformaly high resolution. In this paper the authors show how to increase the temporal resolution of dynamic imaging (e.g., interventional imaging) by using high resolution in areas of expected change and lower resolution elsewhere. Alternatively, by matching the local spatial resolution to the expected edge content of the image, it is possible to avoid the localized truncation artifacts that mark Fourier images reconstructed from the same number of signals. For example, the authors show how proton spectroscopic images of the head may be improved by using high resolution in the neighborhood of scalp lipids that might otherwise cause truncation artifacts.

Optimization of spoiled gradient-echo phase imaging for in vivo localization of a focused ultrasound beam

The parameters of a spoiled gradient-echo (SPGR) pulse sequence have been optimized for in vivo localization of a focused ultrasound beam. Temperature elevation was measured by using the proton resonance frequency shift technique, and the phase difference signal-to-noise ratio (SNR delta phi) was estimated in skeletal muscle and kidney cortex in 10 rabbits. Optimized parameters included the echo time equivalent to T2* of the tissue, the longest repetition time possible with a 20-s sonication, and the flip angle equivalent to the Ernst angle. Optimal SPGR phase imaging can detect a sonication beam with a peak phase difference of 0.55 radian, which corresponds to a temperature elevation of 7.3 degrees C. The sonication beam can be localized within one voxel (0.6 x 0.6 x 5 mm3) at power levels that are below the threshold for thermal damage of the tissue.

Sonochemically produced fluorocarbon microspheres: a new class of magnetic resonance imaging agent

With the intent of increasing the signal-to-noise ratio (SNR) of fluorine magnetic resonance imaging and enabling new applications, we have developed a novel class of agents based on protein encapsulation of fluorocarbons. Microspheres formed by high-intensity ultrasound have a gaussian size distribution with an average diameter of 2.5 microns. As with conventional emulsions, these microspheres target the reticuloendothelial system. However, our sonochemically produced microspheres, because of a high encapsulation efficiency, show increases in the SNR of up to 300% compared to commercially available emulsions. We also demonstrate an increase in the circulation lifetime of the microspheres with the bloodstream by more than 30-fold with a chemical modification of the outer surface of the microsphere. Finally, by encapsulating mixtures of fluorocarbons that undergo solid/liquid phase transitions, we can map temperature in the reticuloendothelial system, with signal changes of approximately 20-fold over a 5 degrees C range.

Simultaneous magnetic resonance phase and magnitude temperature maps in muscle

Noninvasive magnetic resonance temperature maps that are used to monitor thermal ablation of tissue are described. In magnetic resonance images, thermally induced proton nuclear magnetic resonance frequency shifts, and changes in the longitudinal relaxation time produce both phase and magnitude changes in the MR signal. Temperature maps with improved sensitivity are derived from the complex-difference nuclear magnetic resonance signal. Bovine muscle specimens were heated with focused ultrasound to model thermal surgery and create a known thermal distribution to test the method. Resulting MR images acquired in 2 s produce temperature maps with 1 min resolution and 2 degrees C temperature sensitivity. The temperature sensitivity was increased by extending the acquisition to 5 s, by decreasing the receiver bandwidth, and increasing the echo time.

Potential adverse effects of high-intensity focused ultrasound exposure on blood vessels in vivo

The aim of the study was to evaluate the potential adverse effects of high intensity ultrasound exposure on blood vessels during noninvasive focused ultrasound surgery. A hydraulic MR-compatible positioning device was used to manipulate a focused ultrasound transducer (frequency 1.49 MHz, f-number = 0.8) in an MRI scanner. The system was used to sonicate a branch of the femoral artery and vein of 19 rabbits (26 thighs) in vivo at intensity levels above the threshold for transient cavitation; i.e., between 4400 and 8800 W cm-2 with multiple 1 s pulses stepped across the vessels (step size = 0.7 mm). The vessels were located and followed by MR angiography. In 13 rabbits, x-ray angiograms were also performed after the animals were euthanized. The results demonstrated that the 1 s high-intensity exposures caused the arteries to constrict at all exposure levels tested. At the intensity of 5800 W cm-2 and above, the MRI angiogram immediately after the sonications showed no flow. The x-ray angiograms (1-2 h later) showed that the blood vessels were open, but constricted to about 50% or less of their diameter. Both the MR and x-ray angiograms showed that the vessel diameters relaxed toward their initial diameter during the first week after sonication. In five cases, hemorrhage or vessel rupture was caused by the sonication. This study demonstrates that short, high-intensity focused ultrasound exposure can cause vessel spasm and hemorrhage when transient cavitation is present. This condition should be avoided during noninvasive focused ultrasound surgery.

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.

A precise and fast temperature mapping using water proton chemical shift

A new temperature measurement procedure using phase mapping was developed that makes use of the temperature dependence of the water proton chemical shift. Highly accurate and fast measurements were obtained during phantom and in vivo experiments. In the pure water phantom experiments, an accuracy of more than +/- 0.5 degrees C was obtained within a few seconds/slice using a field echo pulse sequence (TR/TE = 115/13 ms, matrix = 128 x 128, number of slices = 5). The temperature dependence of the water proton chemical shift was found to be almost the same for different materials with a chemical composition similar to living tissues (water, glucide, protein). Using this method, the temperature change inside a cat's brain was obtained with an accuracy of more than +/- 1 degree C and an in-plane resolution of 0.6 x 0.6 mm. The temperature measurement error was affected by several factors in the living system (B0 shifts caused by position shifts of the sample, blood flow, etc.), the position shift effect being the most serious.

The estimation of local brain temperature by in vivo 1H magnetic resonance spectroscopy

Brain temperature may be important for investigating pathology and cerebroprotective effects of pharmaceuticals and hypothermia. Two methods for estimating temperature using 1H magnetic resonance spectroscopy are described: a partially water-suppressed binomial sequence and non-water-suppressed point-resolved spectroscopy. Relative to N-acetylaspartate (Naa), water chemical shift (delta H2O-Naa) in piglet brain depended linearly on temperature from 30 degrees to 40 degrees C: temperature was 286.9-94.0 delta H2O-Naa degrees C. Thalamic temperature in six normal infants was 38.1 degrees +/- 0.4 degree C indicating that local brain temperature could be estimated with adequate sensitivity for studying pathologic and therapeutic changes.

MR monitoring of focused ultrasonic surgery of renal cortex: experimental and simulation studies

The aim of the study was to test the hypothesis that magnetic resonance (MR) imaging-guided and -monitored noninvasive ultrasonic surgery can be performed in highly perfused tissues from outside the body. A simulation study was performed to evaluate the optimal sonication parameters. An MR-compatible positioning device was then used to manipulate a focused ultrasound transducer in an MR imager, which was used to sonicate kidneys of five rabbits at various power levels and different durations. Temperature elevation during sonication was monitored with a T1-weighted spoiled gradient-echo sequence. The simulation study demonstrated that a sharply focused transducer and relatively short sonication times (30 seconds or less) are necessary to prevent damage to the overlying skin and muscle tissue, which have a much lower blood perfusion rate than kidney. The experiments showed that the imaging sequence was sensitive enough to show temperature elevation during sonication, thereby indicating the location of the beam focus. Histologic evaluations showed that kidney necrosis could be consistently induced without damage to overlying skin and muscle. The study demonstrated that highly perfused tissues such as the renal cortex can be coagulated from outside the body with focused ultrasound and that MR imaging can be used to guide and monitor this surgery.

Histologic effects of high intensity pulsed ultrasound exposure with subharmonic emission in rabbit brain in vivo

In this study, the threshold for subharmonic emission during in vivo sonication of rabbit brain was investigated. In addition, the histologic effects of pulsed sonication above this threshold were studied. Two spherically curved focused ultrasound transducers with a diameter of 80 mm and a radius of curvature of 70 mm were used in the sonications. The operating frequencies of the transducers were 0.936 and 1.72 MHz. The sonication duration was varied between 0.001 and 1 s and the repetition frequency between 0.1 and 5 Hz. The threshold for subharmonic emission at the frequency of 0.936 MHz was found to be approximately 2000 W cm-2 and 3600 W cm-2 for pulse durations of 1 s and 0.001 s, respectively. The threshold was approximately 1.5-fold as high at a frequency of 1.72 MHz. However, there was considerable variation from experiment to experiment. The multiple pulse experiments at a frequency of 1.72 MHz and an intensity of 7000 W cm-2 showed that the histologic effects ranged from no observable damage of the tissue, to blood-brain barrier breakage, to local haemorrhagia, to local destruction of the tissue, to gross hemorrhage resulting in the death of the animal. The severity of the tissue damage increased as the pulse duration, number of pulses and their repetition frequency increased. The results indicate that the end point of the tissue damage may be controlled by selecting the sonication parameters. Such control over tissue effects can have several different applications when brain disorders are treated.