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

Closed-loop trans-skull ultrasound hyperthermia leads to improved drug delivery from thermosensitive drugs and promotes changes in vascular transport dynamics in brain tumors

Effective drug delivery in brain tumors remains a major challenge in oncology. Although local hyperthermia and stimuli-responsive delivery systems, such as thermosensitive liposomes, represent promising strategies to locally enhance drug delivery in solid tumors and improve outcomes, their application in intracranial malignancies remains unexplored. We hypothesized that the combined abilities of closed-loop trans-skull Magnetic Resonance Imaging guided Focused Ultrasound (MRgFUS) hyperthermia with those of thermosensitive drugs can alleviate challenges in drug delivery and improve survival in gliomas.

Methods: To conduct our investigations, we first designed a closed loop MR-guided Focused Ultrasound (MRgFUS) system for localized trans-skull hyperthermia (ΔT < 0.5 °C) in rodents and established safety thresholds in healthy mice. To assess the abilities of the developed system and proposed therapeutic strategy for FUS-triggered chemotherapy release we employed thermosensitive liposomal Dox (TSL-Dox) and tested it in two different glioma tumor models (F98 in rats and GL261 in mice). To quantify Dox delivery and changes in the transvascular transport dynamics in the tumor microenvironment we combined fluorescent microscopy, dynamic contrast enhanced MRI (DCE-MRI), and physiologically based pharmacokinetic (PBPK) modeling. Lastly, to assess the therapeutic efficacy of the system and of the proposed therapeutic strategy we performed a survival study in the GL261 glioma bearing mice.

Results: The developed closed-loop trans-skull MRgFUS-hyperthermia system that operated at 1.7 MHz, a frequency that maximized the brain (FUS-focus) to skull temperature ratio in mice, was able to attain and maintain the desired focal temperature within a narrow range. Histological evidence (H&E and Nissl) suggests that focal temperature at 41.5 ± 0.5 °C for 10 min is below the threshold for tissue damage. Quantitative analysis of doxorubicin delivery from TSLs with MRgFUS-hyperthermia demonstrated 3.5-fold improvement in cellular uptake in GL261 glioma mouse tumors (p < 0.001) and 5-fold increase in delivery in F98 glioma rat tumors (p < 0.05), as compared to controls (TSL-Dox-only). Moreover, PBPK modeling of drug transport that was calibrated using the experimental data indicated that thermal stress could lead to significant improvement in the transvascular transport (2.3-fold increase in the vessel diffusion coefficient; P < 0.001), in addition to promoting targeted Dox release. Prospective experimental investigations with DCE-MRI during FUS-hyperthermia, supported these findings and provided evidence that moderate thermal stress (≈41 °C for up to 10 min) can promote acute changes in the vascular transport dynamics in the brain tumor microenvironment (K^trans value for control vs. FUS was 0.0097 and 0.0148 min^-1, respectively; p = 0.026). Crucially, survival analysis demonstrated significant improvement in the survival in the TSL-Dox-FUS group as compared to TSL-Dox-only group (p < 0.05), providing supporting evidence on the therapeutic potential of the proposed strategy.

Conclusions: Our investigations demonstrated that spatially controlled thermal stress can be attained and sustained in the mouse brain, using a trans-skull closed-loop MRgFUS system, and used to promote the effective delivery of chemotherapy in gliomas from thermosensitive drugs. This system also allowed us to conduct mechanistic investigations that resulted in the refinement of our understanding on the role of thermal stress in augmenting mass and drug transport in brain tumors. Overall, our study established a new paradigm for effective drug delivery in brain tumors based on closed-loop ultrasound-mediated thermal stress and thermosensitive drugs.

Novel fractionated ultrashort thermal exposures with MRI-guided focused ultrasound for treating tumors with thermosensitive drugs

Thermosensitive liposomes represent an important paradigm in oncology, where hyperthermia-mediated release coupled with thermal bioeffects enhance the effectiveness of chemotherapy. Their widespread clinical adoption hinges upon performing controlled targeted hyperthermia, and a leading candidate to achieve this is temperature-based magnetic resonance imaging (MRI)-guided focused ultrasound (MRgFUS). However, the current approach to hyperthermia involves exposures lasting tens of minutes to hours, which is not possible to achieve in many circumstances because of blood vessel cooling and respiratory motion. Here, we investigate a novel approach to overcome these limitations: to use fractionated ultrashort (~30 s) thermal exposures (~41° to 45°C) to release doxorubicin from a thermosensitive liposome. This is first demonstrated in a dorsal chamber tumor model using two-photon microscopy. Thermal exposures were then conducted with a rabbit tumor model using a custom MRgFUS system incorporating temperature feedback control. Drug release was confirmed, and longitudinal experiments demonstrated profoundly enhanced tumor growth inhibition and survival.

Focused Ultrasound for Noninvasive, Focal Pharmacologic Neurointervention

A long-standing goal of translational neuroscience is the ability to noninvasively deliver therapeutic agents to specific brain regions with high spatiotemporal resolution. Focused ultrasound (FUS) is an emerging technology that can noninvasively deliver energy up the order of 1 kW/cm² with millimeter and millisecond resolution to any point in the human brain with Food and Drug Administration-approved hardware. Although FUS is clinically utilized primarily for focal ablation in conditions such as essential tremor, recent breakthroughs have enabled the use of FUS for drug delivery at lower intensities (i.e., tens of watts per square centimeter) without ablation of the tissue. In this review, we present strategies for image-guided FUS-mediated pharmacologic neurointerventions. First, we discuss blood–brain barrier opening to deliver therapeutic agents of a variety of sizes to the central nervous system. We then describe the use of ultrasound-sensitive nanoparticles to noninvasively deliver small molecules to millimeter-sized structures including superficial cortical regions and deep gray matter regions within the brain without the need for blood–brain barrier opening. We also consider the safety and potential complications of these techniques, with attention to temporal acuity. Finally, we close with a discussion of different methods for mapping the ultrasound field within the brain and describe future avenues of research in ultrasound-targeted drug therapies.

Improved PRF-based MR thermometry using k-space energy spectrum analysis

PURPOSE

Proton resonance frequency (PRF) thermometry encodes information in the phase of MRI signals. A multiplicative factor converts phase changes into temperature changes, and this factor includes the TE. However, phase variations caused by B₀ and/or B₁ inhomogeneities can effectively change TE in ways that vary from pixel to pixel. This work presents how spatial phase variations affect temperature maps and how to correct for corresponding errors.

METHODS

A method called "k-space energy spectrum analysis" was used to map regions in the object domain to regions in the k-space domain. Focused ultrasound heating experiments were performed in tissue-mimicking gel phantoms under two scenarios: with and without proper shimming. The second scenario, with deliberately de-adjusted shimming, was meant to emulate B₀ inhomogeneities in a controlled manner. The TE errors were mapped and compensated for using k-space energy spectrum analysis, and corrected results were compared with reference results. Furthermore, a volunteer was recruited to help evaluate the magnitude of the errors being corrected.

RESULTS

The in vivo abdominal results showed that the TE and heating errors being corrected can readily exceed 10%. In phantom results, a linear regression between reference and corrected temperatures results provided a slope of 0.971 and R² of 0.9964. Analysis based on the Bland-Altman method provided a bias of -0.0977°C and 95% limits of agreement that were 0.75°C apart.

CONCLUSION

Spatially varying TE errors, such as caused by B₀ and/or B₁ inhomogeneities, can be detected and corrected using the k-space energy spectrum analysis method, for increased accuracy in proton resonance frequency thermometry.

Frequency-domain CBE imaging for ultrasound localization of the HIFU focal spot: a feasibility study

High-intensity focused ultrasound (HIFU) is a well-accepted tool for noninvasive thermal therapy. To control the quality of HIFU treatment, the focal spot generated in tissues must be localized. Ultrasound imaging can monitor heated regions; in particular, the change in backscattered energy (CBE) allows parametric imaging to visualize thermal information in the tissue. Conventional CBE imaging constructed in the spatial domain may be easily affected by noises when the HIFU focal spot is visualized. This study proposes frequency-domain CBE imaging to improve noise tolerance and image contrast in HIFU focal spot monitoring. Phantom experiments were performed in a temperature-controlled environment. HIFU of 2.12 MHz was applied to the phantoms, during which a clinical scanner equipped with a 3-MHz convex array transducer was used to collect raw image data consisting of backscattered signals for B-mode, spatial-, and frequency-domain CBE imaging. Concurrently, temperature changes were measured at the focal spot using a thermocouple for comparison with CBE values by calculating the correlation coefficient r. To further analyze CBE image contrast levels, a contrast factor was introduced, and an independent t-test was performed to calculate the probability value p. Experimental results showed that frequency-domain CBE imaging performed well in thermal distribution visualization, enabling quantitative detection of temperature changes. The CBE value calculated in the frequency domain also correlated strongly with that obtained using the conventional spatial-domain approach (r = 0.97). In particular, compared with the image obtained through the conventional method, the contrast of the CBE image obtained using the method based on frequency-domain analysis increased by 2.5-fold (4 dB; p < 0.05). Frequency-domain computations may constitute a new strategy when ultrasound CBE imaging is used to localize the focal spot in HIFU treatment planning.

Recent technological advancements in thermometry

Thermometry is the key factor for achieving successful thermal therapy. Although invasive thermometry with a probe has been used for more than four decades, this method can only detect the local temperature within the probing volume. Noninvasive temperature imaging using a tomographic technique is ideal for monitoring hot-spot formation in the human body. Among various techniques, such as X-ray computed tomography, microwave tomography, echo sonography, and magnetic resonance (MR) imaging, the proton resonance frequency shift method of MR thermometry is the only method currently available for clinical practice because its temperature sensitivity is consistent in most aqueous tissues and can be easily observed using common clinical scanners. New techniques are being proposed to improve the robustness of this method against tissue motion. MR techniques for fat thermometry were also developed based on relaxation times. One of the latest non-MR techniques to attract attention is photoacoustic imaging.

Elementwise approach for simulating transcranial MRI-guided focused ultrasound thermal ablation

This work explored an elementwise approach to model transcranial MRI-guided focused ultrasound (TcMRgFUS) thermal ablation, a noninvasive approach to neurosurgery. Each element of the phased array transducer was simulated individually and could be simultaneously loaded into computer memory, allowing for rapid (~2.5 s) calculation of the pressure field for different phase offsets used for beam steering and aberration correction. We simulated the pressure distribution for 431 sonications in 32 patients, applied the phase and magnitude values used during treatment, and estimated the resulting temperature rise. We systematically varied the relationship between CT (computerized tomography)-derived skull density and the acoustic attenuation and sound speed to obtain the best agreement between the predictions and MR temperature imaging (MRTI). The optimization was validated with simulations of 396 sonications from 40 additional treatments. After optimization, the predicted and measured heating agreed well (R ²: 0.74 patients 1-32; 0.71 patients 33-72). The dimensions and obliquity of the heating in the simulated temperature maps were correlated with the MRTI (R ²: 0.62, 0.74, respectively), but the measured heating was more spatially diffuse. The energy needed to achieve ablation varied by an order of magnitude (3.3-36.1 kJ). While this elementwise approach required more computation time up front (the combined simulation matrices were approximately 4.6 times higher than a single large simulation), it could be performed in parallel on a computing cluster. It allows for rapid calculation of the three-dimensional heating at the focus for different phase and magnitude values on the array. We also show how this approach can be used to optimize the relationship between CT-derived skull density and acoustic properties. While the relationships found here need further validation in a larger patient population, these results demonstrate the promise of this approach to model TcMRgFUS.

Numerical Simulations of the Nonlinear Interaction of a Bubble Cloud and a High Intensity Focused Ultrasound Field

We studied the effects of a small bubble cloud located at the pre-focal area of a high-intensity focused ultrasound field. Our objective is to show that bubbles can modify the bioeffects of an ultrasound treatment in muscle tissue. We model a three-dimensional ultrasound field in an idealized configuration of real operating conditions. Simulations are performed using a combined method based on the Khokhlov-Zabolotskaya-Kuznetsov equation, describing the ultrasound propagation, and a Rayleigh-Plesset equation, modeling the bubble oscillations. The nonlinear interaction of the ultrasound field and the bubble oscillations is considered. Results with and without bubbles for different void fractions of the cloud and different acoustic powers are compared. The cloud induces scattering, nonlinear distortion, and shielding of ultrasound, which increase the mechanical index in the pre-focal zone, shift the location, reduce the size, and modify the shape of the volume of tissue of high mechanical index values, and lower the pressure at the intended focus considerably. Although some hypothesis and parameters used in the models do not fit the real HIFU situations, the simulation results suggest that the effects caused by a bubble cloud located in the pre-focal area should be considered and monitored to ensure the safety of high-intensity focused ultrasound treatments.

MR-guided Focused Ultrasound Facilitates Sonodynamic Therapy with 5-Aminolevulinic Acid in a Rat Glioma Model

Glioblastoma multiforme (GBM) continues to have a dismal prognosis and significant efforts are being made to develop more effective treatment methods. Sonodynamic therapy (SDT) is an emerging modality for cancer treatment which combines ultrasound with sonosensitizers to produce a localized cytotoxic effect. It has long been known that ultrasound exposure can cause both thermal and non-thermal bioeffects and it remains an open question to what degree does temperature impact the efficacy of SDT. In order to optimize the ultrasound parameters of SDT, transcranial MRI-guided focused ultrasound (MRgFUS) and real-time MRI thermometry were used to monitor the therapy in a rat brain tumor model. Experiments were performed using a C6 intracranial glioma tumor model in 37 male Sprague Dawley rats. Treatments were performed about 7 days following tumor implantation when the tumor reached 1-3 mm in diameter as determined by MRI. 5-aminolevulinic acid (5-ALA) was injected at a dose of 60 mg/kg six hours before sonication. MRgFUS at 1.06 MHz was delivered continuously at an in situ spatial-peak temporal-average intensity of 5.5 W/cm2 for 20 min. MR thermometry was acquired to monitor the temperature change in the brain during sonication. The tumor growth response for animals receiving 5-ALA alone, FUS alone, 5-ALA + FUS and a sham control group were evaluated with MRI every week following treatment. During 20 min of MRgFUS at 5.5 W/cm2, the temperature within the targeted brain tumor was elevated from 32.3 ± 0.5 °C and 37.2 ± 0.7 °C to 33.2 ± 0.9 °C and 38.4 ± 1.1 °C, respectively. Both the tumor growth inhibition and survival were significantly improved in the 5-ALA + FUS group with 32 °C or 37 °C as the starting core body (rectal) temperature. 5-ALA alone and FUS alone did not improve survival. These promising results indicate that relatively low power continuous wave transcranial MRgFUS in conjunction with 5-ALA can produce an inhibitory effect on rat brain tumor growth in the absence of thermal dose. Further investigation of the ultrasound parameters is needed to improve the therapeutic efficacy of MRgFUS and 5-ALA.

Magnetic resonance thermometry and its biological applications - Physical principles and practical considerations

Most parameters that influence the magnetic resonance imaging (MRI) signal experience a temperature dependence. The fact that MRI can be used for non-invasive measurements of temperature and temperature change deep inside the human body has been known for over 30 years. Today, MR temperature imaging is widely used to monitor and evaluate thermal therapies such as radio frequency, microwave, laser, and focused ultrasound therapy. In this paper we cover the physical principles underlying the biological applications of MR temperature imaging and discuss practical considerations and remaining challenges. For biological tissue, the MR signal of interest comes mostly from hydrogen protons of water molecules but also from protons in, e.g., adipose tissue and various metabolites. Most of the discussed methods, such as those using the proton resonance frequency (PRF) shift, T₁, T₂, and diffusion only measure temperature change, but measurements of absolute temperatures are also possible using spectroscopic imaging methods (taking advantage of various metabolite signals as internal references) or various types of contrast agents. Currently, the PRF method is the most used clinically due to good sensitivity, excellent linearity with temperature, and because it is largely independent of tissue type. Because the PRF method does not work in adipose tissues, T₁- and T₂-based methods have recently gained interest for monitoring temperature change in areas with high fat content such as the breast and abdomen. Absolute temperature measurement methods using spectroscopic imaging and contrast agents often offer too low spatial and temporal resolution for accurate monitoring of ablative thermal procedures, but have shown great promise in monitoring the slower and usually less spatially localized temperature change observed during hyperthermia procedures. Much of the current research effort for ablative procedures is aimed at providing faster measurements, larger field-of-view coverage, simultaneous monitoring in aqueous and adipose tissues, and more motion-insensitive acquisitions for better precision measurements in organs such as the heart, liver, and kidneys. For hyperthermia applications, larger coverage, motion insensitivity, and simultaneous aqueous and adipose monitoring are also important, but great effort is also aimed at solving the problem of long-term field drift which gets interpreted as temperature change when using the PRF method.

Microbubble-assisted MRI-guided focused ultrasound for hyperthermia at reduced power levels

PURPOSE

Ultrasound contrast agent microbubbles were combined with magnetic resonance imaging (MRI)-guided focused ultrasound (MRgFUS) as a means to achieve mild hyperthermia at reduced power levels.

METHODS

MRgFUS hyperthermia (42°C for 20 min) was evaluated in rabbit thigh muscle or Vx2 tumors using infusions of microbubbles (Definity, 20 µL/kg) or saline (sham) administered over 5 min. The impact of treatments on drug uptake was assessed with liposomal doxorubicin (Caelyx, 2.5 mg/kg). Applied power levels before and after the injection of microbubbles or saline were compared, and drug uptake was evaluated with fluorometry of tissues harvested 24 hr post-treatment.

RESULTS

MRgFUS hyperthermia in muscle and tumors resulted in accurate temperature control (mean =42.0°C, root mean square error (RMSE) = 0.3°C). The power dropped significantly following the injection of microbubbles in muscle and tumors compared to exposures without microbubbles (-21.9% ± 12.5% vs -5.9% ± 7.8%, p = .009 in muscle; -33.8% ± 9.9% vs -3.0% ± 7.2%, p < .001 in tumors). Cavitation monitoring indicated emission of subharmonic, ultraharmonic, and elevated levels of fourth to sixth harmonic frequencies following microbubble injection. The drug delivery was elevated significantly in muscle with the use of microbubble-assisted relative to conventional heating (0.5 ± 0.5 ng/mg vs 0.20 ± 0.04 ng/mg, p = .05), whereas in tumors similar levels were found (11 ± 3 ng/mg vs 16 ± 4 ng/mg, p = .13).

CONCLUSIONS

The finding that microbubbles reduce the applied power requirements for hyperthermia has considerable clinical implications. The elevated levels of drug found in muscle but not tumor tissue suggest a complex interplay between the heating effects of microbubbles with those of enhanced permeabilization and possible vascular damage.

Combined passive acoustic mapping and magnetic resonance thermometry for monitoring phase-shift nanoemulsion enhanced focused ultrasound therapy

Focused ultrasound (FUS) has the potential to enable precise, image-guided noninvasive surgery for the treatment of cancer in which tumors are identified and destroyed in a single integrated procedure. However, success of the method in highly vascular organs has been limited due to heat losses to perfusion, requiring development of techniques to locally enhance energy absorption and heating. In addition, FUS procedures are conventionally monitored using MRI, which provides excellent anatomical images and can map temperature, but is not capable of capturing the full gamut of available data such as the acoustic emissions generated during this inherently acoustically-driven procedure. Here, we employed phase-shift nanoemulsions (PSNE) embedded in tissue phantoms to promote cavitation and hence temperature rise induced by FUS. In addition, we incorporated passive acoustic mapping (PAM) alongside simultaneous MR thermometry in order to visualize both acoustic emissions and temperature rise, within the bore of a full scale clinical MRI scanner. Focal cavitation of PSNE could be resolved using PAM and resulted in accelerated heating and increased the maximum elevated temperature measured via MR thermometry compared to experiments without nanoemulsions. Over time, the simultaneously acquired acoustic and temperature maps show translation of the focus of activity towards the FUS transducer, and the magnitude of the increase in cavitation and focal shift both increased with nanoemulsion concentration. PAM results were well correlated with MRI thermometry and demonstrated greater sensitivity, with the ability to detect cavitation before enhanced heating was observed. The results suggest that PSNE could be beneficial for enhancement of thermal focused ultrasound therapies and that PAM could be a critical tool for monitoring this process.

MR techniques for guiding high-intensity focused ultrasound (HIFU) treatments

To make full use of the ability of magnetic resonance (MR) to guide high-intensity focused ultrasound (HIFU) treatment, effort has been made to improve techniques for thermometry, motion tracking, and sound beam visualization. For monitoring rapid temperature elevation with proton resonance frequency (PRF) shift, data acquisition and processing can be accelerated with parallel imaging and/or sparse sampling in conjunction with appropriate signal processing methods. Thermometry should be robust against tissue motion, motion-induced magnetic field variation, and susceptibility change. Thus, multibaseline, referenceless, or hybrid techniques have become important. In cases with adipose or bony tissues, for which PRF shift cannot be used, thermometry with relaxation times or signal intensity may be utilized. Motion tracking is crucial not only for thermometry but also for targeting the focus of an ultrasound in moving organs such as the liver, kidney, or heart. Various techniques for motion tracking, such as those based on an anatomical image atlas with optical-flow displacement detection, a navigator echo to seize the diaphragm position, and/or rapid imaging to track vessel positions, have been proposed. Techniques for avoiding the ribcage and near-field heating have also been examined. MR acoustic radiation force imaging (MR-ARFI) is an alternative to thermometry that can identify the location and shape of the focal spot and sound beam path. This technique could be useful for treating heterogeneous tissue regions or performing transcranial therapy. All of these developments, which will be discussed further in this review, expand the applicability of HIFU treatments to a variety of clinical targets while maintaining safety and precision.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 4 J. Magn. Reson. Imaging 2018;47:316-331.

Optothermal profile of an ablation catheter with integrated microcoil for MR-thermometry during Nd:YAG laser interstitial thermal therapies of the liver—an in-vitro experimental and theoretical study

PURPOSE

Flexible microcoils integrated with ablation catheters can improve the temperature accuracy during local MR-thermometry in Nd:YAG laser interstitial thermal therapies. Here, the authors are concerned with obtaining a preliminary confirmation of the clinical utility of the modified catheter. They investigate whether the thin-film substrate and copper tracks of the printed coil inductor affect the symmetry of the thermal profile, and hence of the lesion produced.

METHODS

Transmission spectroscopy in the near infrared was performed to test for the attenuation at 1064 nm through the 25 μm thick Kapton substrate of the microcoil. The radial transmission profile of an infrared high-power, light emitting diode with >80% normalized power at 1064 nm was measured through a cross section of the modified applicator to assess the impact of the copper inductor on the optical profile. The measurements were performed in air, as well as with the applicator surrounded by two types of scattering media; crystals of NaCl and a layer of liver-mimicking gel phantom. A numerical model based on Huygens-Fresnel principle and finite element simulations, using a commercially available package (COMSOL Multiphysics), were employed to compare with the optical measurements. The impact of the modified optical profile on the thermal symmetry was assessed by examining the high resolution microcoil derived thermal maps from a Nd:YAG laser ablation performed on a liver-mimicking gel phantom.

RESULTS

Less than 30% attenuation through the Kapton film was verified. Shadowing behind the copper tracks was observed in air and the measured radial irradiation correlated well with the diffraction pattern calculated numerically using the Huygens-Fresnel principle. Both optical experiments and simulations, demonstrate that shadowing is mitigated by the scattering properties of a turbid medium. The microcoil derived thermal maps at the end of a Nd:YAG laser ablation performed on a gel phantom in a 3 T scanner confirm that the modified irradiation pattern does not disrupt the thermal symmetry, even though, unlike tissue, the gel is minimally scattering.

CONCLUSIONS

The results from this initial assessment indicate that microcoils can be safely integrated with ablation catheters and ensure that the complete necrosis of the liver tumor can still be achieved.

MRI-controlled ultrasound thermal therapy

Advances in medical imaging have enabled the development of new minimally and completely noninvasive therapies that produce a desired biological effect in a target, such as a tumor, with minimal damage to the surrounding tissue. One means of noninvasively achieving bioeffects in tissue is the use of ultrasound to generate heat. Specialized ultrasound transducers can be used to generate focal regions of heating non invasively, without inserting anything into the body or affecting the tissue outside the target region. Ultrasound thermal therapy can be used with magnetic resonance (MR) imaging (MRI) guidance and MRI temperature feedback to automatically control temperature distributions during heating, producing accurate thermal lesions, or maintaining optimal conditions to enhance drug delivery.

High intensity focused ultrasound (HIFU) focal spot localization using harmonic motion imaging (HMI)

Several ultrasound-based imaging modalities have been proposed for image guidance and monitoring of high-intensity focused ultrasound (HIFU) treatment. However, accurate localization and characterization of the effective region of treatment (focal spot) remain important obstacles in the clinical implementation of HIFU ablation. Harmonic motion imaging for focused ultrasound (HMIFU) is a HIFU monitoring technique that utilizes radiation-force-induced localized oscillatory displacement. HMIFU has been shown to correctly identify the formation and extent of HIFU thermal ablation lesions. However a significant problem remains in identifying the location of the HIFU focus, which is necessary for treatment planning. In this study, the induced displacement was employed to localize the HIFU focal spot inside the tissue prior to treatment. Feasibility was shown with two separate systems. The 1D HMIFU system consisted of a HIFU transducer emitting an amplitude-modulated HIFU beam for mechanical excitation and a confocal single-element, pulse-echo transducer for simultaneous RF acquisition. The 2D HIFU system consists of a HIFU phased array, and a co-axial imaging phased array for simultaneous imaging. Initial feasibility was first performed on tissue-mimicking gelatin phantoms and the focal zone was defined as the region corresponding to the -3dB full width at half maximum of the HMI displacement. Using the same parameters, in vitro experiments were performed in canine liver specimens to compare the defined focal zone with the lesion. In vitro measurements showed good agreement between the HMI predicted focal zone and the induced HIFU lesion location. HMIFU was experimentally shown to be capable of predicting and tracking the focal region in both phantoms and in vitro tissues. The accuracy of focal spot localization was evaluated by comparing with the lesion location in post-ablative tissues, with a R(2) = 0.821 at p < 0.002 in the 2D HMI system. We demonstrated the feasibility of using this HMI-based technique to localize the HIFU focal spot without inducing thermal changes during the planning phase. The focal spot localization method has also been applied on ex vivo human breast tissue ablation and can be fully integrated into any HMI system for planning purposes.

MRI compatible head phantom for ultrasound surgery

OBJECTIVE

Develop a magnetic resonance imaging (MRI) compatible head phantom with acoustic attenuation closely matched to the human attenuation, and suitable for testing focused ultrasound surgery protocols.

MATERIALS AND METHODS

Images from an adult brain CT scan were used to segment the skull bone from adjacent cerebral tissue. The segmented model was manufactured in a 3-D printer using (Acrylonitrile Butadiene Styrene) ABS plastic. The cerebral tissue was mimicked by an agar-evaporated milk-silica gel (2% w/v-25% v/v-1.2% w/v) which was molded inside a skull model.

RESULTS

The measured attenuation of the ABS skull was 16 dB/cm MHz. The estimated attenuation coefficient of the gel replicating brain tissue was 0.6 dB/cm MHz. The estimated agar-silica gel's T1 and T2 relaxation times in a 1.5 Tesla magnetic field were 852 ms and 66 ms respectively. The effectiveness of the skull to reduce ultrasonic heating was demonstrated using MRI thermometry.

CONCLUSION

Due to growing interest in using MRI guided focused ultrasound (MRgFUS) for treating brain cancer and its application in sonothrombolysis, the proposed head phantom can be utilized as a very useful tool for evaluating ultrasonic protocols, thus minimizing the need for animal models and cadavers.

Enhanced lesion-to-bubble ratio on ultrasonic Nakagami imaging for monitoring of high-intensity focused ultrasound

OBJECTIVES

This work explored the feasibility of using ultrasonic Nakagami imaging to enhance the contrast between thermal lesions and bubbles induced by high-intensity focused ultrasound (US) in a transparent tissue-mimicking phantom at different acoustic power levels.

METHODS

The term "lesion-to-bubble ratio" was proposed and defined as the ratio of the scattered power from the thermal lesion to the scattered power from the bubbles calculated in the various monitoring of images for high-intensity focused US. Two-dimensional radiofrequency data backscattered from the exposed region were captured by a modified diagnostic US scanner to estimate the Nakagami statistical parameter, m, and reconstruct the ultrasonic B-mode images and Nakagami parameter images. The dynamic changes in the lesion-to-bubble ratio over the US exposure procedure were calculated simultaneously and compared among video photos, B-mode images, and Nakagami images for monitoring of high-intensity focused US.

RESULTS

After a small thermal lesion was induced by high-intensity focused US in the phantom, the lesion-to-bubble ratio values corresponding to the video photo, B-mode image, and Nakagami image were 5.3, 1, and 9.8 dB, respectively. When a large thermal lesion appeared in the phantom, the ratio values increased to 7.2, 3, and 14 dB. During US exposure, the ratio values calculated for the video photo, B-mode image, and Nakagami image began to increase gradually and rose to peak values of 8.3, 2.9, and 14.8 dB at the end of the US exposure.

CONCLUSIONS

This preliminary study on a tissue-mimicking phantom suggests that Nakagami imaging may have a potential use in enhancing the lesion-to-bubble ratio for monitoring high-intensity focused US. Further studies in vivo and in vitro will be needed to evaluate the potential applications for high-intensity focused US.

Simultaneous imaging of in vivo conductivity and susceptibility

PURPOSE

Approaches for quantitative mapping of electric conductivity and magnetic susceptibility using MRI have been developed independently. The purpose of this study is to present a method to simultaneously acquire information on conductivity and susceptibility and to produce images based on these properties.

METHODS

A 3D multiecho gradient-echo sequence was used. Phase evolution during the multiecho was used to produce quantitative susceptibility maps, while the phase value at zero echo time was retrieved, and used to generates quantitative conductivity maps. Electromagnetic simulations were performed to evaluate the phase distribution due to conductivity variations. Phantom and in vivo data were also acquired to assess the quality of images produced.

RESULTS

Simulations demonstrated that phase differences across objects increases with size and conductivity. For an accurate conductivity estimate, the maximum echo time was approximately equal to the true T2* value in order to achieve signal-to-noise ratio maximization. The most accurate susceptibility was obtained when separating phase contribution from conductivity. Phantom and in vivo results showed good quality images representing the electromagnetic properties.

CONCLUSION

A simultaneous quantitative electromagnetic property imaging approach is demonstrated here. The approach not only improves the efficiency of mapping electromagnetic properties, but can also improve the accuracy of susceptibility mapping by separating image phases introduced by conductivity and susceptibility.

The evaluation of steerable ultrasonic catheters for minimally invasive MRI-guided cardiac ablation

PURPOSE

The purpose of this study was to develop steerable MR-compatible ultrasound catheters suitable for minimally invasive MRI-guided cardiac ablation therapies.

METHODS

MRI-compatible ultrasound steerable catheters were developed and tested for their overall tissue heating performance and safety. Ultrasound transducers were mounted on a monodirectional deflectable catheter tip that was made to be MRI-compatible. Catheter safety was assessed on the potential to form hot spots at the distal end of the catheter throughout fast spin echo and thermometry scans. Heating experiments were performed on phantoms and ex vivo porcine cardiac samples.

RESULTS

During catheter safety experiments, a maximum temperature increase of 11.35 ± 0.83°C was evident after a 12-min, 40-s fast spin echo scan with a whole body specific absorption rate (SAR) of 1.9 W/kg and 1.07 ± 0.22°C during thermometry scans (flip angle = 90°; scan time = 12 min, 41 s; whole body SAR = 0.34 W/kg). Temperature elevations induced by the sonication were shown to be on the order of 38.1 ± 5.2°C for phantom experiments and 49.3 ± 9.7°C for ex vivo cardiac samples.

CONCLUSION

Steerable ultrasound catheters have the potential to be safely placed in an MR system with little concern of catheter self-heating and driven to heat surrounding structures to cause ablations. In addition, these catheters have the added benefit of a deflectable tip that allows the treatment of multiple targets from within the bore of the MR scanner.

Multiparametric fat-water separation method for fast chemical-shift imaging guidance of thermal therapies

PURPOSE

A k-means-based classification algorithm is investigated to assess suitability for rapidly separating and classifying fat/water spectral peaks from a fast chemical shift imaging technique for magnetic resonance temperature imaging. Algorithm testing is performed in simulated mathematical phantoms and agar gel phantoms containing mixed fat/water regions.

METHODS

Proton resonance frequencies (PRFs), apparent spin-spin relaxation (T2*) times, and T1-weighted (T1-W) amplitude values were calculated for each voxel using a single-peak autoregressive moving average (ARMA) signal model. These parameters were then used as criteria for k-means sorting, with the results used to determine PRF ranges of each chemical species cluster for further classification. To detect the presence of secondary chemical species, spectral parameters were recalculated when needed using a two-peak ARMA signal model during the subsequent classification steps. Mathematical phantom simulations involved the modulation of signal-to-noise ratios (SNR), maximum PRF shift (MPS) values, analysis window sizes, and frequency expansion factor sizes in order to characterize the algorithm performance across a variety of conditions. In agar, images were collected on a 1.5T clinical MR scanner using acquisition parameters close to simulation, and algorithm performance was assessed by comparing classification results to manually segmented maps of the fat/water regions.

RESULTS

Performance was characterized quantitatively using the Dice Similarity Coefficient (DSC), sensitivity, and specificity. The simulated mathematical phantom experiments demonstrated good fat/water separation depending on conditions, specifically high SNR, moderate MPS value, small analysis window size, and low but nonzero frequency expansion factor size. Physical phantom results demonstrated good identification for both water (0.997 ± 0.001, 0.999 ± 0.001, and 0.986 ± 0.001 for DSC, sensitivity, and specificity, respectively) and fat (0.763 ± 0.006, 0.980 ± 0.004, and 0.941 ± 0.002 for DSC, sensitivity, and specificity, respectively). Temperature uncertainties, based on PRF uncertainties from a 5 × 5-voxel ROI, were 0.342 and 0.351°C for pure and mixed fat/water regions, respectively. Algorithm speed was tested using 25 × 25-voxel and whole image ROIs containing both fat and water, resulting in average processing times per acquisition of 2.00 ± 0.07 s and 146 ± 1 s, respectively, using uncompiled MATLAB scripts running on a shared CPU server with eight Intel Xeon(TM) E5640 quad-core processors (2.66 GHz, 12 MB cache) and 12 GB RAM.

CONCLUSIONS

Results from both the mathematical and physical phantom suggest the k-means-based classification algorithm could be useful for rapid, dynamic imaging in an ROI for thermal interventions. Successful separation of fat/water information would aid in reducing errors from the nontemperature sensitive fat PRF, as well as potentially facilitate using fat as an internal reference for PRF shift thermometry when appropriate. Additionally, the T1-W or R2* signals may be used for monitoring temperature in surrounding adipose tissue.

Outcome of uterine artery embolization versus MR-guided high-intensity focused ultrasound treatment for uterine fibroids: long-term results

OBJECTIVES

To compare the long-term outcome after uterine artery embolization (UAE) versus magnetic resonance-guided high-intensity focused ultrasound (MR-g HIFU) for symptomatic uterine fibroids.

METHODS

Seventy-seven women (median age, 39.3 years; range, 29.2-52.2 years) with symptomatic uterine fibroids, equally eligible for UAE and MR-g HIFU based on our exclusion criteria underwent treatment (UAE, N = 41; MR-g HIFU, N = 36) from 2002 to 2009 at our institution. Symptom severity (SS) and total health-related quality of life (Total HRQoL) scores were assessed by the uterine fibroid symptom and quality of life (UFS-QoL) questionnaire before treatment and at long-term follow-up after UAE (median 61.9 months) and after MR-g HIFU (median: 60.7 months). Re-intervention rates were assessed for each therapy and compared.

RESULTS

Re-intervention was significantly lower after UAE (12.2%) than after MR-g HIFU (66.7%) at long-term follow-up (p<0.001). After UAE changes in SS (50 pre-treatment vs. 6.3 post-treatment) and Total HRQoL (57.8 pre-treatment vs. 100 post-treatment) were significantly better than changes in SS (42.2 pre-treatment vs. 26.6 post-treatment) and Total HRQoL score (66.4 pre-treatment vs. 87.9 post-treatment) after MR-g HIFU (p = 0.019 and 0.049 respectively).

CONCLUSIONS

Improvement of SS and Total HRQoL scores was significantly better after UAE resulting in a significant lower re-intervention rate compared to MR-g HIFU.

A multi-element interstitial ultrasound applicator for the thermal therapy of brain tumors

Interstitial thermal therapy is a minimally invasive treatment modality that has been used clinically for ablating both primary and secondary brain tumors. Here a multi-element interstitial ultrasound applicator is described that allows for increased spatial control during thermal ablation of tumors as compared to existing clinical devices. The device consists of an array of 56 ultrasound elements operating at 6 MHz, oriented on the seven faces of a 3.2 mm flexible catheter. The device was first characterized using the acoustic holography method to examine the functioning of the array. Then experiments were performed to measure heating in tissue-mimicking gel phantoms and ex vivo tissue samples using magnetic resonance imaging-based thermometry. Experimental measurements were compared with results obtained using numerical simulations. Last, simulations were performed to study the feasibility of using the device for thermal ablation in the brain. Experimental results show that the device can be used to induce a temperature rise of greater than 20 °C in ex vivo tissue samples and numerical simulations further demonstrate that tumors with diameters of greater than 30-mm could potentially be treated.

Toward real-time temperature monitoring in fat and aqueous tissue during magnetic resonance-guided high-intensity focused ultrasound using a three-dimensional proton resonance frequency T1 method

PURPOSE

To present a three-dimensional (3D) segmented echoplanar imaging (EPI) pulse sequence implementation that provides simultaneously the proton resonance frequency shift temperature of aqueous tissue and the longitudinal relaxation time (T1 ) of fat during thermal ablation.

METHODS

The hybrid sequence was implemented by combining a 3D segmented flyback EPI sequence, the extended two-point Dixon fat and water separation, and the double flip angle T1 mapping techniques. High-intensity focused ultrasound (HIFU) heating experiments were performed at three different acoustic powers on excised human breast fat embedded in ex vivo porcine muscle. Furthermore, T1 calibrations with temperature in four different excised breast fat samples were performed, yielding an estimate of the average and variation of dT1 /dT across subjects.

RESULTS

The water only images were used to mask the complex original data before computing the proton resonance frequency shift. T1 values were calculated from the fat-only images. The relative temperature coefficients were found in five fat tissue samples from different patients and ranged from 1.2% to 2.6%/°C.

CONCLUSION

The results demonstrate the capability of real-time simultaneous temperature mapping in aqueous tissue and T1 mapping in fat during HIFU ablation, providing a potential tool for treatment monitoring in organs with large fat content, such as the breast.

Feasibility of using Nakagami distribution in evaluating the formation of ultrasound-induced thermal lesions

The acoustic posterior shadowing effects of bubbles influence the accuracy for defining the location and range of ablated thermal lesions during focused ultrasound surgery when using ultrasonic monitoring imaging. This paper explored the feasibility of using Nakagami distribution to evaluate the ablated region induced by focused ultrasound exposures at different acoustic power levels in transparent tissue-mimicking phantoms. The mean value of the Nakagami parameter m was about 0.5 in the cavitation region and increased to around 1 in the ablated region. Nakagami images were not subject to significant shadowing effects of bubbles. Ultrasound-induced thermal lesions observed in the photos and Nakagami images were overshadowed by bubbles in the B-mode images. The lesion size predicted in the Nakagami images was smaller than that predicted in the photos due to the sub resolvable effect of Nakagami imaging at the interface. This preliminary study on tissue-mimicking phantom suggested that the Nakagami parameter m may have the potential use in evaluating the formation of ultrasound-induced thermal lesion when the shadowing effect of bubbles is strong while the thermal lesion was small. Further studies in vivo and in vitro will be needed to evaluate the potential application.

Fat-referenced MR thermometry in the breast and prostate using IDEAL

PURPOSE

To demonstrate a three-echo fat-referenced MR thermometry technique that estimates and corrects for time-varying phase disturbances in heterogeneous tissues.

MATERIALS AND METHODS

Fat protons do not exhibit a temperature-dependent frequency shift. Fat-referenced thermometry methods exploit this insensitivity and use the signal from fat to measure and correct for magnetic field disturbances. In this study, we present a fat-referenced method that uses interpolation of the fat signal to correct for phase disturbances in fat free regions. Phantom and ex vivo tissue cool-down experiments were performed to evaluate the accuracy of this method in the absence of motion. Non-heated in vivo imaging of the breast and prostate was performed to demonstrate measurement robustness in the presence of systemic and motion-induced field disturbances. Measurement accuracy of the method was compared to conventional proton resonance frequency shift MR thermometry.

RESULTS

In the ex vivo porcine tissue experiment, maximum measurement error of the fat-referenced method was reduced 42% from 3.3 to 1.9°C when compared to conventional MR thermometry. In the breasts, measurement errors were reduced by up to 70% from 6.4 to 1.9°C.

CONCLUSION

Ex vivo and in vivo results show that the proposed method reduces measurement errors in the heterogeneous tissue experiments when compared to conventional MR thermometry.

Irreversible change in the T1 temperature dependence with thermal dose using the proton resonance frequency-T1 technique

Denaturation of macromolecules within the tissues is believed to be the major factor contributing to the damage of tissues upon hyperthermia. As a result, the value of the spin-lattice relaxation time T1 of the tissue water, which is related to the translational and rotational rates of water, represents an intrinsic probe for investigating structural changes in tissues at high temperature. Therefore, the goal of this work is to investigate whether the simultaneous measurement of temperature and T1 using a hybrid proton resonance frequency (PRF)-T1 measurement technique can be used to detect irreversible changes in T1 that might be indicative of tissue damage. A new hybrid PRF-T1 sequence was implemented based on the variable flip angle driven-equilibrium single-pulse observation (DESPOT)1 method from a standard three dimensional segmented echo-planar imaging sequence by alternating two flip angles from measurement to measurement. The structural changes of the heated tissue volumes were analyzed based on the derived T1 values and the corresponding PRF temperatures. Using the hybrid PRF-T1 technique, we demonstrate that the change of spin lattice relaxation time T1 is reversible with temperature for low thermal dose (thermal dose ≤ 240 cumulative equivalent minutes [CEM] 43°C) and irreversible with temperature after significant accumulation of thermal dose in ex vivo chicken breast tissue. These results suggest that the hybrid PRF-T1 method may be a potentially powerful tool to investigate the extent and mechanism of heat damage of biological tissues.

Magnetic resonance thermometry at 7T for real-time monitoring and correction of ultrasound induced mild hyperthermia

While Magnetic Resonance Thermometry (MRT) has been extensively utilized for non-invasive temperature measurement, there is limited data on the use of high field (≥7T) scanners for this purpose. MR-guided Focused Ultrasound (MRgFUS) is a promising non-invasive method for localized hyperthermia and drug delivery. MRT based on the temperature sensitivity of the proton resonance frequency (PRF) has been implemented in both a tissue phantom and in vivo in a mouse Met-1 tumor model, using partial parallel imaging (PPI) to speed acquisition. An MRgFUS system capable of delivering a controlled 3D acoustic dose during real time MRT with proportional, integral, and derivative (PID) feedback control was developed and validated. Real-time MRT was validated in a tofu phantom with fluoroptic temperature measurements, and acoustic heating simulations were in good agreement with MR temperature maps. In an in vivo Met-1 mouse tumor, the real-time PID feedback control is capable of maintaining the desired temperature with high accuracy. We found that real time MR control of hyperthermia is feasible at high field, and k-space based PPI techniques may be implemented for increasing temporal resolution while maintaining temperature accuracy on the order of 1°C.

RF field visualization of RF ablation at the Larmor frequency

Radio-frequency ablation (RFA) is an effective minimally invasive treatment for tumors. One primary source of difficulty is monitoring and controlling the ablation region. Currently, RFA is performed at 460 kHz, for which magnetic resonance imaging (MRI) could play a role given its capability for temperature monitoring and tumor visualization. If instead the ablation were to be performed at the MRI Larmor frequency, then the MR capability for B(1) field mapping could be used to directly visualize the radio-frequency (RF) fields created by the ablation currents. Visualizing the RF fields may enable better control of the ablation currents, enabling better control of lesion shape and size and improving repeatability. We demonstrate the feasibility of performing RFAs at 64 MHz and show preliminary results from imaging the RF fields from the ablation. The post-ablation RF fields show an increase in current density in the ablated region, consistent with an increase in conductivity of the ablated tissue.

MR temperature imaging of nanoshell mediated laser ablation

Minimally invasive thermal therapy using high-power diode lasers is an active area of clinical research. Gold nanoshells (AuNS) can be tuned to absorb light in the range used for laser ablation and may facilitate more conformal tumor heating and sparing of normal tissue via enhanced tumor specific heating. This concept was investigated in a xenograft model of prostate cancer (PC-3) using MR temperature imaging (MRTI) in a 1.5T scanner to characterize the spatiotemporal temperature distribution resulting from nanoparticle mediated heating. Tumors with and without intravenously injected AuNS were exposed to an external laser tuned to 808 nm for 180 sec at 4 W/cm(2) under real-time monitoring with proton resonance frequency shift based MRTI. Microscopy indicated that these nanoparticles (140-150 nm) accumulated passively in the tumor and remained close to the tumor microvasculature. MRTI measured a statistically significant (p < 0.001) increase in maximum temperature in the tumor cortex (mean = 21 ± 7°C) in +AuNS tumors versus control tumors. Analysis of the temperature maps helped demonstrate that the overall distribution of temperature within +AuNS tumors was demonstrably higher versus control, and resulted in damage visible on histopathology. This research demonstrates that passive uptake of intravenously injected AuNS in PC-3 xenografts converts the tumor vasculature into a potent heating source for nanoparticle mediated ablation at power levels which do not generate significant damage in normal tissue. When used in conjunction with MRTI, this has implications for development and validation of more conformal delivery of therapy for interstitial laser ablations.

Measurement of temperature changes in cooling dead rats using magnetic resonance thermometry

Magnetic resonance imaging thermometry has been introduced as a technique for measurement of temperature changes in cooling dead rats. Rat pelvic magnetic resonance images were acquired sequentially more than 2h after euthanasia by halothane overdose. A series of temperature difference maps in cooling dead rats was obtained with calculating imaging phase changes induced by the water proton frequency shift caused by temperature changes. Different cooling processes were monitored by the temperature difference maps in the rats. Magnetic resonance imaging thermometry applied in the study of laboratory animals could theoretically reproduce a variety of causes of death with different environmental conditions. Outcomes from experimental animal studies could be translated into a temperature-based time of death estimation in forensics.

Ultrafast 1D MR thermometry using phase or frequency mapping

OBJECT

To develop an ultrafast MRI-based temperature monitoring method for application during rapid ultrasound exposures in moving organs.

MATERIALS AND METHODS

A slice selective 90° - 180° pair of RF pulses was used to solicit an echo from a column, which was then sampled with a train of gradient echoes. In a gel phantom, phase changes of each echo were compared to standard gradient-echo thermometry, and temperature monitoring was tested during focused ultrasound sonications. Signal-to-noise ratio (SNR) performance was evaluated in vivo in a rabbit brain, and feasibility was tested in a human heart.

RESULTS

The correlation between each echo in the acquisition and MRI-based temperature measurements was good (R = 0.98 ± 0.03). A temperature sampling rate of 19 Hz was achieved at 3T in the gel phantom. It was possible to acquire the water frequency in the beating heart muscle with 5-Hz sampling rate during a breath hold.

CONCLUSION

Ultrafast thermometry via phase or frequency monitoring along single columns was demonstrated. With a temporal resolution around 50 ms, it may be possible to monitor focal heating produced by short ultrasound pulses.

Combining two-dimensional spatially selective RF excitation, parallel imaging, and UNFOLD for accelerated MR thermometry imaging

MR thermometry can be a very challenging application, as good resolution may be needed along spatial, temporal, and temperature axes. Given that the heated foci produced during thermal therapies are typically much smaller than the anatomy being imaged, much of the imaged field-of-view is not actually being heated and may not require temperature monitoring. In this work, many-fold improvements were obtained in terms of temporal resolution and/or 3D spatial coverage by sacrificing some of the in-plane spatial coverage. To do so, three fast-imaging approaches were jointly implemented with a spoiled gradient echo sequence: (1) two-dimensional spatially selective RF excitation, (2) unaliasing by Fourier encoding the overlaps using the temporal dimension (UNFOLD), and (3) parallel imaging. The sequence was tested during experiments with focused ultrasound heating in ex vivo tissue and a tissue-mimicking phantom. Temperature maps were estimated from phase-difference images based on the water proton resonance frequency shift. Results were compared to those obtained from a spoiled gradient echo sequence sequence, using a t-test. Temporal resolution was increased by 24-fold, with temperature uncertainty less than 1°C, while maintaining accurate temperature measurements (mean difference between measurements, as observed in gel = 0.1°C ± 0.6; R = 0.98; P > 0.05).

Multipathway sequences for MR thermometry

MR-based thermometry is a valuable adjunct to thermal ablation therapies as it helps to determine when lethal doses are reached at the target and whether surrounding tissues are safe from damage. When the targeted lesion is mobile, MR data can further be used for motion-tracking purposes. The present work introduces pulse sequence modifications that enable significant improvements in terms of both temperature-to-noise-ratio properties and target-tracking abilities. Instead of sampling a single magnetization pathway as in typical MR thermometry sequences, the pulse-sequence design introduced here involves sampling at least one additional pathway. Image reconstruction changes associated with the proposed sampling scheme are also described. The method was implemented on two commonly used MR thermometry sequences: the gradient-echo and the interleaved echo-planar imaging sequences. Data from the extra pathway enabled temperature-to-noise-ratio improvements by up to 35%, without increasing scan time. Potentially of greater significance is that the sampled pathways featured very different contrast for blood vessels, facilitating their detection and use as internal landmarks for tracking purposes. Through improved temperature-to-noise-ratio and lesion-tracking abilities, the proposed pulse-sequence design may facilitate the use of MR-monitored thermal ablations as an effective treatment option even in mobile organs such as the liver and kidneys.

High intensity focused ultrasound in clinical tumor ablation

Recent advances in high intensity focused ultrasound (HIFU), which was developed in the 1940s as a viable thermal tissue ablation approach, have increased its popularity. In clinics, HIFU has been applied to treat a variety of solid malignant tumors in a well-defined volume, including the pancreas, liver, prostate, breast, uterine fibroids, and soft-tissue sarcomas. In comparison to conventional tumor/cancer treatment modalities, such as open surgery, radio- and chemo-therapy, HIFU has the advantages of non-invasion, non-ionization, and fewer complications after treatment. Over 100 000 cases have been treated throughout the world with great success. The fundamental principles of HIFU ablation are coagulative thermal necrosis due to the absorption of ultrasound energy during transmission in tissue and the induced cavitation damage. This paper reviews the clinical outcomes of HIFU ablation for applicable cancers, and then summarizes the recommendations for a satisfactory HIFU treatment according to clinical experience. In addition, the current challenges in HIFU for engineers and physicians are also included. More recent horizons have broadened the application of HIFU in tumor treatment, such as HIFU-mediated drug delivery, vessel occlusion, and soft tissue erosion (“histotripsy”). In summary, HIFU is likely to play a significant role in the future oncology practice.

MRI-Guided Electrophysiology Intervention

Catheter ablation is a first-line treatment for many cardiac arrhythmias and is generally performed under X-ray fluoroscopy guidance. However, current techniques for ablating complex arrhythmias such as atrial fibrillation and ventricular tachycardia are associated with sub-optimal success rates and prolonged radiation exposure. Pre-procedure 3-D magnetic resonance imaging (MRI) has improved understanding of the anatomic basis of complex arrhythmias and is being used for planning and guidance of ablation procedures. A particular strength of MRI compared to other imaging modalities is the ability to visualize ablation lesions. Post-procedure MRI is now being applied to assess ablation lesion location and permanence with the goal of identifying factors leading to procedure success and failure. In the future, intra-procedure real-time MRI, together with the ability to image complex 3-D arrhythmogenic anatomy and target additional ablation to regions of incomplete lesion formation, may allow for more successful treatment of even complex arrhythmias without exposure to ionizing radiation. Development of clinical grade MRI-compatible electrophysiology devices is required to transition intra-procedure MRI from preclinical studies to more routine use in patients.

MR imaging-guided interventions in the genitourinary tract: an evolving concept

MR imaging-guided interventions are well established in routine patient care in many parts of the world. There are many approaches, depending on magnet design and clinical need, based on MR imaging providing excellent inherent tissue contrast without ionizing radiation risk for patients. MR imaging-guided minimally invasive therapeutic procedures have advantages over conventional surgical procedures. In the genitourinary tract, MR imaging guidance has a role in tumor detection, localization, and staging and can provide accurate image guidance for minimally invasive procedures. The advent of molecular and metabolic imaging and use of higher strength magnets likely will improve diagnostic accuracy and allow targeted therapy to maximize disease control and minimize side effects.

An internal reference model-based PRF temperature mapping method with Cramer-Rao lower bound noise performance analysis

The conventional phase difference method for MR thermometry suffers from disturbances caused by the presence of lipid protons, motion-induced error, and field drift. A signal model is presented with multi-echo gradient echo (GRE) sequence using a fat signal as an internal reference to overcome these problems. The internal reference signal model is fit to the water and fat signals by the extended Prony algorithm and the Levenberg-Marquardt algorithm to estimate the chemical shifts between water and fat which contain temperature information. A noise analysis of the signal model was conducted using the Cramer-Rao lower bound to evaluate the noise performance of various algorithms, the effects of imaging parameters, and the influence of the water:fat signal ratio in a sample on the temperature estimate. Comparison of the calculated temperature map and thermocouple temperature measurements shows that the maximum temperature estimation error is 0.614 degrees C, with a standard deviation of 0.06 degrees C, confirming the feasibility of this model-based temperature mapping method. The influence of sample water:fat signal ratio on the accuracy of the temperature estimate is evaluated in a water-fat mixed phantom experiment with an optimal ratio of approximately 0.66:1.

Dynamic changes of integrated backscatter, attenuation coefficient and bubble activities during high-intensity focused ultrasound (HIFU) treatment

This paper simultaneously investigated the transient characteristics of integrated backscatter (IBS), attenuation coefficient and bubble activities as time traces before, during and after HIFU treatment, with different HIFU parameters (acoustic power and duty cycle) in both transparent tissue-mimicking phantoms and freshly excised bovine livers. These dynamic changes of acoustic parameters and bubble activities were correlated with the visualization of lesion development selected from photos, conventional B-mode ultrasound images and differential IBS images over the whole procedure of HIFU treatment. Two-dimensional radiofrequency (RF) data were acquired by a modified diagnostic ultrasound scanner to estimate the changes of mean IBS and attenuation coefficient averaged in the lesion region, and to construct the differential IBS images and B-mode ultrasound images simultaneously. Bubble activities over the whole procedure of HIFU treatment were investigated by the passive cavitation detection (PCD) method and the changes in subharmonic and broadband noise were correlated with the transient characteristics of IBS and attenuation coefficient. When HIFU was switched on, IBS and attenuation coefficient increased with the appearance of bubble clouds in the B-mode and differential IBS image. At the same time, the level of subharmonic and broadband noise rose abruptly. Then, there was an initial decrease in the attenuation coefficient, followed by an increase when at lower HIFU power. As the lesion appeared, IBS and attenuation coefficient both increased rapidly to a value twice that of normal. Then the changes in IBS and attenuation coefficient showed more complex patterns, but still showed a slower trend of increases with lesion development. Violent bubble activities were visible in the gel and were evident as strongly echogenic regions in the differential IBS images and B-mode images simultaneously. This was detected by a dramatic high level of subharmonic and broadband noise at the same time. These bubble activities caused fluctuations in IBS and attenuation coefficient during HIFU treatment. After HIFU, IBS and attenuation coefficient decreased gradually accompanied by the fadeout of bright hyperechoic spot in the B-mode and differential IBS image, but were still higher than normal when they were stable. The increases of IBS and attenuation coefficient were greater when using higher acoustic power or a higher duty cycle of the therapeutic emission. These experiments indicated that the bubble activities had the dominant effects on the transient characteristics of IBS and attenuation. This should be taken into consideration when using the dynamic acoustic-property changes for the potentially real-time monitoring imaging of HIFU treatment.

Improved volumetric MR-HIFU ablation by robust binary feedback control

Volumetric high-intensity focused ultrasound (HIFU) guided by multiplane magnetic resonance (MR) thermometry has been shown to be a safe and efficient method to thermally ablate large tissue volumes. However, the induced temperature rise and thermal lesions show significant variability, depending on exposure parameters, such as power and timing, as well as unknown tissue parameters. In this study, a simple and robust feedback-control method that relies on rapid MR thermometry to control the HIFU exposure during heating is introduced. The binary feedback algorithm adjusts the durations of the concentric ablation circles within the target volume to reach an optimal temperature. The efficacy of the binary feedback control was evaluated by performing 90 ablations in vivo and comparing the results with simulations. Feedback control of the sonications improved the reproducibility of the induced lesion size. The standard deviation of the diameter was reduced by factors of 1.9, 7.2, 5.0, and 3.4 for 4-, 8-, 12-, and 16-mm lesions, respectively. Energy efficiency was also improved, as the binary feedback method required less energy to create the desired lesion. These results show that binary feedback improves the quality of volumetric ablation by consistently producing thermal lesions of expected size while reducing the required energy as well.

Volumetric HIFU ablation under 3D guidance of rapid MRI thermometry

A volumetric sonication method is proposed that produces volume ablations by steering the focal point along a predetermined trajectory consisting of multiple concentric outward-moving circles. This method was tested in vivo on pig thigh muscle (32 ablations in nine animals). Trajectory diameters were 4, 12, and 16 mm with sonication duration depending on the trajectory size and ranging from 20 to 73 s. Despite the larger trajectories requiring more energy to reach necrosis within the desired volume, the ablated volume per unit applied energy increased with trajectory size, indicating improved treatment efficiency for larger trajectories. The higher amounts of energy required for the larger trajectories also increased the risk of off-focus heating, especially along the beam axis in the near field. To avoid related adverse effects, rapid volumetric multiplane MR thermometry was introduced for simultaneous monitoring of the temperature and thermal dose evolution along the beam axis and in the near field, as well as in the target region with a total coverage of six slices acquired every 3 s. An excellent correlation was observed between the thermal dose and both the nonperfused (R=0.929 for the diameter and R=0.964 for the length) and oedematous (R=0.913 for the diameter and R=0.939 for the length) volumes as seen in contrast-enhanced T1-weighted difference images and T2-weighted postsonication images, respectively. Histology confirmed the presence of a homogeneous necrosis inside the heated volumes. These results show that volumetric high-intensity focused ultrasound (HIFU) sonication allows for efficiently creating large thermal lesions while reducing treatment duration and also that the rapid multiplane MR thermometry improves the safety of the therapeutic procedure by monitoring temperature evolution both inside as well as outside the targeted volume.

Computational aspects in high intensity ultrasonic surgery planning

Therapeutic ultrasound treatment planning is discussed and computational aspects regarding it are reviewed. Nonlinear ultrasound simulations were solved with a combined frequency domain Rayleigh and KZK model. Ultrasonic simulations were combined with thermal simulations and were used to compute heating of muscle tissue in vivo for four different focused ultrasound transducers. The simulations were compared with measurements and good agreement was found for large F-number transducers. However, at F# 1.9 the simulated rate of temperature rise was approximately a factor of 2 higher than the measured ones. The power levels used with the F# 1 transducer were too low to show any nonlinearity. The simulations were used to investigate the importance of nonlinarities generated in the coupling water, and also the importance of including skin in the simulations. Ignoring either of these in the model would lead to larger errors. Most notably, the nonlinearities generated in the water can enhance the focal temperature by more than 100%. The simulations also demonstrated that pulsed high power sonications may provide an opportunity to significantly (up to a factor of 3) reduce the treatment time. In conclusion, nonlinear propagation can play an important role in shaping the energy distribution during a focused ultrasound treatment and it should not be ignored in planning. However, the current simulation methods are accurate only with relatively large F-numbers and better models need to be developed for sharply focused transducers.

Temperature mapping of thermal ablation using MRI

MRI is a unique tool for minimally invasive thermal ablation in that it can provide both targeting, monitoring and control during the procedure. Monitoring is achieved by using MRI temperature mapping. In this review the relevant physics is explained as a background to the state-of-the-art methods for computing temperature maps as well as the more cutting edge methods. The review covers both methods to monitor heating and cooling of tissue and explains temperature mapping using Proton Resonance Frequency shift, T1 mapping, diffusion mapping, R2* mapping and thermal models.

MR thermometry-based feedback control of efficacy and safety in minimum-time thermal therapies: Phantom andin-vivoevaluations

The experimental validation of a model-based, thermal therapy control system which automatically and simultaneously achieves the specified efficacy and safety objectives of the treatment is reported. MR-thermometry measurements are used in real-time to control the power of a stationary, focused ultrasound transducer in order to achieve the desired treatment outcome in minimum time without violating the imposed safety constraints. Treatment efficacy is quantified in terms of the thermal dose delivered to the target. Normal tissue safety is ensured by automatically maintaining normal tissue temperature below the imposed limit in the user-specified locations. To reflect hardware limitations, constraints on the maximum applied power are also imposed. At the pretreatment stage, MR imaging and thermometry are used to localize the treatment target and identify thermal and actuation models. The results of phantom and canine experiments demonstrate that spatially-distributed, real-time MR temperature measurements enhance one's ability to robustly achieve the desired treatment outcome in minimum time without violating safety constraints. Post-treatment evaluation of the outcome using T2-weighted images of canine muscle showed good spatial correlation between the sonicated area and thermally damaged tissue.