Study of laser ablation in the in vivo rabbit brain with MR thermometry

Experimental assessment of microwave ablation computational modeling with MR thermometry

PURPOSE

Computational models are widely used during the design and characterization of microwave ablation (MWA) devices, and have been proposed for pretreatment planning. Our objective was to assess three-dimensional (3D) transient temperature and ablation profiles predicted by MWA computational models with temperature profiles measured experimentally using magnetic resonance (MR) thermometry in ex vivo bovine liver.

MATERIALS AND METHODS

We performed MWA in ex vivo tissue under MR guidance using a custom, 2.45 GHz water-cooled applicator. MR thermometry data were acquired for 2 min prior to heating, during 5-10 min microwave exposures, and for 3 min following heating. Fiber-optic temperature sensors were used to validate the accuracy of MR temperature measurements. A total of 13 ablation experiments were conducted using 30-50 W applied power at the applicator input. MWA computational models were implemented using the finite element method, and incorporated temperature-dependent changes in tissue physical properties. Model-predicted ablation zone extents were compared against MRI-derived Arrhenius thermal damage maps using the Dice similarity coefficient (DSC).

RESULTS

Prior to heating, the observed standard deviation of MR temperature data was in the range of 0.3-0.7°C. Mean absolute error between MR temperature measurements and fiber-optic temperature probes during heating was in the range of 0.5-2.8°C. The mean DSC between model-predicted ablation zones and MRI-derived Arrhenius thermal damage maps for 13 experimental set-ups was 0.95. When comparing simulated and experimentally (i.e. using MRI) measured temperatures, the mean absolute error (MAE %) relative to maximum temperature change was in the range 5%-8.5%.

CONCLUSION

We developed a system for characterizing 3D transient temperature and ablation profiles with MR thermometry during MWA in ex vivo liver tissue, and applied the system for experimental validation of MWA computational models.

Optical Fiber Temperature Sensors and Their Biomedical Applications

The use of sensors in the real world is on the rise, providing information on medical diagnostics for healthcare and improving quality of life. Optical fiber sensors, as a result of their unique properties (small dimensions, capability of multiplexing, chemical inertness, and immunity to electromagnetic fields) have found wide applications, ranging from structural health monitoring to biomedical and point-of-care instrumentation. Furthermore, these sensors usually have good linearity, rapid response for real-time monitoring, and high sensitivity to external perturbations. Optical fiber sensors, thus, present several features that make them extremely attractive for a wide variety of applications, especially biomedical applications. This paper reviews achievements in the area of temperature optical fiber sensors, different configurations of the sensors reported over the last five years, and application of this technology in biomedical applications.

Improved MR thermometry for laser interstitial thermotherapy

OBJECTIVES

To develop, test and evaluate improved 2D and 3D protocols for proton resonance frequency shift magnetic resonance temperature imaging (MRTI) of laser interstitial thermal therapy (LITT). The objective was to develop improved MRTI protocols in terms of temperature measurement precision and volume coverage compared to the 2D MRTI protocol currently used with a commercially available LITT system.

METHODS

Four different 2D protocols and four different 3D protocols were investigated. The 2D protocols used multi-echo readouts to prolong the total MR sampling time and hence the MRTI precision, without prolonging the total acquisition time. The 3D protocols provided volumetric thermometry by acquiring a slab of 12 contiguous slices in the same acquisition time as the 2D protocols. The study only considered readily available pulse sequences (Cartesian 2D and 3D gradient recalled echo and echo planar imaging [EPI]) and methods (partial Fourier and parallel imaging) to ensure wide availability and rapid clinical implementation across vendors and field strengths. In vivo volunteer studies were performed to investigate and compare MRTI precision and image quality. Phantom experiments with LITT heating were performed to investigate and compare MRTI precision and accuracy. Different coil setups were used in the in vivo studies to assess precision differences between using local (such as flex and head coils) and non-local (i.e., body coil) receive coils. Studies were performed at both 1.5 T and 3 T.

RESULTS

The improved 2D protocols provide up to a factor of two improvement in the MRTI precision in the same acquisition time, compared to the currently used clinical protocol. The 3D echo planar imaging protocols provide comparable precision as the currently used 2D clinical protocol, but over a substantially larger field of view, without increasing the acquisition time. As expected, local receive coils perform substantially better than the body coil, and 3 T provides better MRTI accuracy and precision than 1.5 T. 3D data can be zero-filled interpolated in all three dimensions (as opposed to just two dimensions for 2D data), reducing partial volume effects and measuring higher maximum temperature rises.

CONCLUSIONS

With the presented protocols substantially improved MRTI precision (for 2D imaging) or greatly improved field of view coverage (for 3D imaging) can be achieved in the same acquisition time as the currently used protocol. Only widely available pulse sequences and acquisition methods were investigated, which should ensure quick translation to the clinic. Lasers Surg. Med. 51:286-300, 2019. © 2019 Wiley Periodicals, Inc.

Stereotactic Laser Interstitial Thermal Therapy for Recurrent High-Grade Gliomas

BACKGROUND

The value of maximal safe cytoreductive surgery in recurrent high-grade gliomas (HGGs) is gaining wider acceptance. However, patients may harbor recurrent tumors that may be difficult to access with open surgery. Laser interstitial thermal therapy (LITT) is emerging as a technique for treating a variety of brain pathologies, including primary and metastatic tumors, radiation necrosis, and epilepsy.

OBJECTIVE

To review the role of LITT in the treatment of recurrent HGGs, for which current treatments have limited efficacy, and to discuss the possible role of LITT in the disruption of the blood-brain barrier to increase delivery of chemotherapy locoregionally.

METHODS

A MEDLINE search was performed to identify 17 articles potentially appropriate for review. Of these 17, 6 reported currently commercially available systems and as well as magnetic resonance thermometry to monitor the ablation and, thus, were thought to be most appropriate for this review. These studies were then reviewed for complications associated with LITT. Ablation volume, tumor coverage, and treatment times were also reviewed.

RESULTS

Sixty-four lesions in 63 patients with recurrent HGGs were treated with LITT. Frontal (n = 34), temporal (n = 14), and parietal (n = 16) were the most common locations. Permanent neurological deficits were seen in 7 patients (12%), vascular injuries occurred in 2 patients (3%), and wound infection was observed in 1 patient (2%). Ablation coverage of the lesions ranged from 78% to 100%.

CONCLUSION

Although experience using LITT for recurrent HGGs is growing, current evidence is insufficient to offer a recommendation about its role in the treatment paradigm for recurrent HGGs.

ABBREVIATIONS

BBB, blood-brain barrierFDA, US Food and Drug AdministrationGBM, glioblastoma multiformeHGG, high-grade gliomaLITT, laser interstitial thermal therapy.

A methodology for thermal dose model parameter development using perioperative MRI

Post-treatment imaging is the principal method for evaluating thermal lesions following image-guided thermal ablation procedures. While real-time temperature feedback using magnetic resonance temperature imaging (MRTI) is a complementary tool that can be used to optimise lesion size throughout the procedure, a thermal dose model is needed to convert temperature-time histories to estimates of thermal damage. However, existing models rely on empirical parameters derived from laboratory experiments that are not direct indicators of post-treatment radiologic appearance. In this work, we investigate a technique that uses perioperative MR data to find novel thermal dose model parameters that are tailored to the appearance of the thermal lesion on post-treatment contrast-enhanced imaging. Perioperative MR data were analysed for five patients receiving magnetic resonance-guided laser-induced thermal therapy (MRgLITT) for brain metastases. The characteristic enhancing ring was manually segmented on post-treatment T1-weighted imaging and registered into the MRTI geometry. Post-treatment appearance was modelled using a coupled Arrhenius-logistic model and non-linear optimisation techniques were used to find the maximum-likelihood kinetic parameters and dose thresholds that characterise the inner and outer boundary of the enhancing ring. The parameter values and thresholds were consistent with previous investigations, while the average difference between the predicted and segmented boundaries was on the order of one pixel (1 mm). The areas predicted using the optimised model parameters were also within 1 mm of those predicted by clinically utilised dose models. This technique makes clinically acquired data available for investigating new thermal dose model parameters driven by clinically relevant endpoints.

Magnetic Resonance Thermometry and Laser Interstitial Thermal Therapy for Brain Tumors

Recent technological advancements in intraoperative imaging are shaping the way for a new era in brain tumor surgery. Magnetic resonance thermometry has provided intraoperative real-time imaging feedback for safe and effective application of laser interstitial thermal therapy (LITT) in neuro-oncology. Thermal ablation has also established itself as a surgical option in epilepsy surgery and is currently used in spine oncology with promising results. This article reviews the principles and rationale as well as the clinical application of LITT for brain tumors. It also discusses the technical nuances of the current commercially available systems.

Image-guided thermal ablation with MR-based thermometry

Thermal ablation techniques such as radiofrequency, microwave, high intensity focused ultrasound (HIFU) and laser have been used as minimally invasive strategies for the treatment of variety of cancers. MR thermometry methods are readily available for monitoring thermal distribution and deposition in real time, leading to decrease of incidents of normal tissue damage around targeted lesion. HIFU and laser-induced thermal therapy (LITT) are the two widely accepted tumor ablation techniques because of their compatibility with MR systems. MRI provides multiple temperature dependent parameters for thermal imaging, such as signal intensity, T1, T2, diffusion coefficient, magnetization transfer, proton resonance frequency shift (PRFS, including phase imaging and spectroscopy) as well as frequency shift of temperature sensitive contrast agents. Absolute temperature mapping techniques, including both spectroscopic imaging using metabolites as a reference and phase imaging using fat as a reference, are immune to susceptibility effects and are not dependent on phase differences. These techniques are intrinsically more reliable than relative temperature measurement by phase mapping methods. If the limitation of low temporal and spatial resolution could be overcome, these methods may be preferred for MR-guided thermal ablation systems. As of today, the most popular MR thermal imaging method applied in tumor thermal ablation surgery is, however, still PRFS based phase mapping technique, which only provides relative temperature change and is prone to motion artifacts.

ANALYTICAL SOLUTIONS IN THE MODELING OF THE LOCAL RF ABLATION

Coupled mathematical models for the radiofrequency (RF) ablation performed in biomedical sciences have been developed based on the bioheat transfer theory. The heat exchange problem is important to be analytically studied in order to control the size of the necrosis zone caused by RF ablation. This lesion size in the tissue may be predicted by the knowledge of the internal tissue temperature. We propose an analytical solution for the Pennes heat transfer equation in bi- and tri-region domains, applicable to the RF ablation of cancerigeneous tissue — a clinical relevant problem. The model consists of two partial differential equations describing the spatio-temporal interactions between the electric and thermic effects. The aim is to find simple algebraic expressions of analytical solutions that may allow to generate quantitative results which in turn may be interpreted (including uncertainties). The dependence of the temperature as function of the electrothermal parameters in both diseased and surrounding healthy tissues is pointed out. Two cases, namely the tumor–tissue and tumor–tissue–skin systems, are graphically computed, and important findings include the fact that the presence of tissue with smaller value parameters protects somehow healthy cells. Moreover, the graphical representations are conducted to highlight the link of the profile of current density distribution in the physiological problem with the (neither oval nor circular) shape of the temperature isoclinic lines.

Tissue Ablation Dynamics During Magnetic Resonance-Guided, Laser-Induced Thermal Therapy

BACKGROUND

Magnetic resonance-guided, laser-induced thermal therapy is a real-time magnetic resonance thermometry-guided, minimally invasive procedure used in the treatment of intracranial tumors, epilepsy, and pain. Little is known about its dynamics and the effects of various pathologies on overall ablation.

OBJECTIVE

To determine the relationship between thermal energy delivery and the time to maximal estimated thermal damage and whether differences exist between various intracranial pathologies.

METHODS

We used real-time ablation data from 28 patients across 5 unique intracranial pathologies. All ablations were performed using the Visualase Thermal Therapy System (Medtronic, Inc, Minneapolis, Minnesota), which uses a 980-nm diffusing tip diode laser. The thermal damage area was plotted against time for each ablation. We then estimated the duration of time required to reach 50% (t50) and 97% (t97) of maximal damage. Comparisons were then made between different intracranial pathologies.

RESULTS

The duration required to reach maximal thermal damage estimate (TDE) among all ablations was 159 ± 62 seconds, and the t50 and t97 were 43 ± 21 and 136 ± 57 seconds, respectively, where t97 was reached at an average of 23 seconds before the maximal TDE. The t97 was shorter in the recurrent metastasis/radiation necrosis and epilepsy groups compared with the previously untreated glioblastoma multiforme group.

CONCLUSION

The optimal duration can be estimated by the t97, which can be achieved in less than 3 minutes and differs across ablation targets. TDE expansion decelerates with prolonged ablation. Future studies are needed to examine the radiographic and clinical outcomes as well as the effects of ablation power, irrigation speed, and the effect of previous therapies on ablation dynamics.

Renaissance of laser interstitial thermal ablation

Laser interstitial thermal therapy (LITT) is a minimally invasive technique for treating intracranial tumors, originally introduced in 1983. Its use in neurosurgical procedures was historically limited by early technical difficulties related to the monitoring and control of the extent of thermal damage. The development of magnetic resonance thermography and its application to LITT have allowed for real-time thermal imaging and feedback control during laser energy delivery, allowing for precise and accurate provision of tissue hyperthermia. Improvements in laser probe design, surgical stereotactic targeting hardware, and computer monitoring software have accelerated acceptance and clinical utilization of LITT as a neurosurgical treatment alternative. Current commercially available LITT systems have been used for the treatment of neurosurgical soft-tissue lesions, including difficult to access brain tumors, malignant gliomas, and radiosurgery-resistant metastases, as well as for the ablation of such lesions as epileptogenic foci and radiation necrosis. In this review, the authors aim to critically analyze the literature to describe the advent of LITT as a neurosurgical, laser excision tool, including its development, use, indications, and efficacy as it relates to neurosurgical applications.

Noninvasive, in vivo determination of uterine fibroid thermal conductivity in MRI-guided high intensity focused ultrasound therapy

PURPOSE

To estimate the local thermal conductivity of uterine fibroid in vivo at a high temperature range (60-80°C) typically encountered in magnetic resonance imaging-guided high-intensity focused ultrasound (MRgHIFU) surgery. The thermal conductivity of uterine fibroids in vivo is unknown and knowledge about tissue thermal conductivity may aid in effective delivery of thermal energy for ablation.

MATERIALS AND METHODS

All subjects (nine women) provided written informed consent to participate in this Institutional Review Board-approved study. A total of 10 fibroids were treated using MRgHIFU surgery with real-time temperature monitoring during both heating and cooling periods. The local thermal conductivity was determined by analyzing the spatiotemporal spread of temperature during the cooling period.

RESULTS

The thermal conductivity of MRgHIFU-treated uterine fibroids was 0.47 ± 0.07 W·m(-1) ·K(-1) (range: 0.25∼0.67 W·m(-1) ·K(-1) ) which is slightly lower than the reported value for skeletal muscle at temperatures of <40°C (0.52 to 0.62 W·m(-1) ·K(-1) ).

CONCLUSION

It is possible to estimate the thermal conductivity of uterine fibroids in vivo from the spatiotemporal spread of temperature around the HIFU focus during the cooling period.

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.

Correlation between the temperature dependence of intrinsic MR parameters and thermal dose measured by a rapid chemical shift imaging technique

In order to investigate simultaneous MR temperature imaging and direct validation of tissue damage during thermal therapy, temperature-dependent signal changes in proton resonance frequency (PRF) shifts, R(2)* values, and T1-weighted amplitudes are measured from one technique in ex vivo tissue. Using a multigradient echo acquisition and the Stieglitz-McBride algorithm, the temperature sensitivity coefficients of these parameters are measured in each tissue at high spatiotemporal resolutions (1.6 x 1.6 x 4 mm 3,≤ 5sec) at the range of 25-61 °C. Non-linear changes in MR parameters are examined and correlated with an Arrhenius rate dose model of thermal damage. Using logistic regression, the probability of changes in these parameters is calculated as a function of thermal dose to determine if changes correspond to thermal damage. Temperature sensitivity of R(2)* and, in some cases, T1-weighted amplitudes are statistically different before and after thermal damage occurred. Significant changes in the slopes of R(2)* as a function of temperature are observed. Logistic regression analysis shows that these changes could be accurately predicted using the Arrhenius rate dose model (Ω = 1.01 ± 0.03), thereby showing that the changes in R(2)* could be direct markers of protein denaturation. Overall, by using a chemical shift imaging technique with simultaneous temperature estimation, R(2)* mapping and T1-W imaging, it is shown that changes in the sensitivity of R(2)* and, to a lesser degree, T1-W amplitudes are measured in ex vivo tissue when thermal damage is expected to occur. These changes could possibly be used for direct validation of thermal damage in contrast to model-based predictions.

Threshold-based prediction of the coagulation zone in sequential temperature mapping in MR-guided radiofrequency ablation of liver tumours

OBJECTIVE

To evaluate different cut-off temperature levels for a threshold-based prediction of the coagulation zone in magnetic resonance (MR)-guided radiofrequency (RF) ablation of liver tumours.

METHODS

Temperature-sensitive measurements were acquired during RF ablation of 24 patients with primary (6) and secondary liver lesions (18) using a wide-bore 1.5 T MR sytem and compared with the post-interventional coagulation zone. Temperature measurements using the proton resonance frequency shift method were performed directly subsequent to energy application. The temperature maps were registered on the contrast-enhanced follow-up MR images acquired 4 weeks after treatment. Areas with temperatures above 50°, 55° and 60°C were segmented and compared with the coagulation zones. Sensitivity and positive predictive value were calculated.

RESULTS

No major complications occurred and all tumours were completely treated. No tumour recurrence was observed at the follow-up examination after 4 weeks. Two patients with secondary liver lesions showed local tumour recurrence after 4 and 7 months. The 60°C threshold level achieved the highest positive predictive value (87.7 ± 9.9) and the best prediction of the coagulation zone.

CONCLUSIONS

For a threshold-based prediction of the coagulation zone, the 60°C cut-off level achieved the best prediction of the coagulation zone among the tested levels.

KEY POINTS

• Temperature monitoring can be used to survey MR-guided radiofrequency ablation • The developing ablation zone can be estimated based on post-interventional temperature measurements • A 60°C threshold level can be used to predict the ablation zone • The 50°C and 55°C temperature zones tend to overestimate the ablation zone.

Thresholds for thermal damage to normal tissues: an update

The purpose of this review is to summarise a literature survey on thermal thresholds for tissue damage. This review covers published literature for the consecutive years from 2002-2009. The first review on this subject was published in 2003. It included an extensive discussion of how to use thermal dosimetric principles to normalise all time-temperature data histories to a common format. This review utilises those same principles to address sensitivity of a variety of tissues, but with particular emphasis on brain and testis. The review includes new data on tissues that were not included in the original review. Several important observations have come from this review. First, a large proportion of the papers examined for this review were discarded because time-temperature history at the site of thermal damage assessment was not recorded. It is strongly recommended that future research on this subject include such data. Second, very little data is available examining chronic consequences of thermal exposure. On a related point, the time of assessment of damage after exposure is critically important for assessing whether damage is transient or permanent. Additionally, virtually no data are available for repeated thermal exposures which may occur in certain recreational or occupational activities. For purposes of regulatory guidelines, both acute and lasting effects of thermal damage should be considered.

Clinical evaluation of MR temperature monitoring of laser-induced thermotherapy in human liver using the proton-resonance-frequency method and predictive models of cell death

PURPOSE

To assess the feasibility, precision, and accuracy of real-time temperature mapping (TMap) during laser-induced thermotherapy (LITT) for clinical practice in patients liver with a gradient echo (GRE) sequence using the proton resonance frequency (PRF) method.

MATERIALS AND METHODS

LITT was performed on 34 lesions in 18 patients with simultaneous real-time visualization of relative temperature changes. Correlative contrast-enhanced T1-weighted magnetic resonance (MR) images of the liver were acquired after treatment using the same slice positions and angulations as TMap images acquired during LITT. For each slice, TMap and follow-up images were registered for comparison. Afterwards, segmentation based on temperature (T) >52°C on TMap and based on necrosis seen on follow-up images was performed. These segmented structures were overlaid and divided into zones where the TMap was found to either over- or underestimate necrosis on the postcontrast images. Regions with T>52°C after 20 minutes were defined as necrotic tissue based on data received from two different thermal dose models.

RESULTS

The average intersecting region of TMap and necrotic zone was 87% ± 5%, the overestimated 13% ± 4%, and the underestimated 13% ± 5%.

CONCLUSION

This study demonstrates that MR temperature mapping appears reasonably capable of predicting tissue necrosis on the basis of indicating regions having greater temperatures than 52°C and could be used to monitor and adjust the thermal therapy appropriately during treatment.

Modeling of laser coagulation of tissue with MRI temperature monitoring

Light energy from a laser source that is delivered into body tissue via a fiber-optic probe with minimal invasiveness has been used to ablate solid tumors. This thermal coagulation process can be guided and monitored accurately by continuous magnetic resonance imaging (MRI) since the laser energy delivery system does not interfere with MRI. This report deals with mathematical modeling and analysis of laser coagulation of tissue. This model is intended for "real-time" analysis of magnetic resonance images obtained during the coagulation process to guide clinical treatment. A mathematical model is developed to simulate the thermal response of tissue to a laser light heating source. For fast simulation, an approximate solution of the thermal model is used to predict the dynamics of temperature distribution and tissue damage induced by a laser energy line source. The validity of these simulations is tested by comparison with MRI-based temperature data acquired from in vivo experiments in rabbits. The model-simulated temperature distribution and predicted lesion dynamics correspond closely with MRI-based data. These results demonstrate the potential for using this combination of fast modeling and MRI technologies during laser heating of tissue for online prediction of tumor lesion size during laser heating.

Quantitative comparison of thermal dose models in normal canine brain

PURPOSE

Minimally invasive thermal ablative therapies as alternatives to conventional surgical management of solid tumors and other pathologies is increasing owing to the potential benefits of performing these procedures in an outpatient setting with reduced complications and comorbidity. Magnetic resonance temperature imaging (MRTI) measurement allows existing thermal dose models to use the spatiotemporal temperature history to estimate the thermal damage to tissue. However, the various thermal dose models presented in the literature employ different parameters and thresholds, affecting the reliability of thermal dosimetry. In this study, the authors quantitatively compared three thermal dose models (Arrhenius rate process, CEM43, and threshold temperature) using the dice similarity coefficient (DSC).

METHODS

The DSC was used to compare the spatial overlap between the region of thermal damage as predicted by the models for in vivo normal canine brain during thermal therapy to the region of thermal damage as revealed by contrast-enhanced T1-weighted images acquired immediately after therapy (< 20 min). The outer edge of the hyperintense rim of the ablation region was used as the surrogate marker for the limits of thermal coagulation. The DSC was also used to investigate the impact of varying the thresholds on each models' ability to predict the zone of thermal necrosis.

RESULTS

At previously reported thresholds, the authors found that all three models showed good agreement (defined as DSC > 0.7) with post-treatment imaging. All three models examined across the range of commonly applied thresholds consistently showed highly accurate spatial overlap, low variability, and little dependence on temperature uncertainty. DSC values corresponding to cited thresholds were not significantly different from peak DSC values.

CONCLUSIONS

Thus, the authors conclude that the all three thermal dose models can be used as a reliable surrogate for postcontrast tissue damage verification imaging in rapid ablation procedures and can also be used to enhance the capability of MRTI to control thermal therapy in real time.

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.

Quantitative MRI-based temperature mapping based on the proton resonant frequency shift: Review of validation studies

MRI-based temperature imaging that exploits the temperature-sensitive water proton resonant frequency shift is currently the only available method for reliable quantification of temperature changes in vivo. Extensive pre-clinical work has been performed to validate this method for guiding thermal therapies. That work has shown the method to be useful for all stages of the thermal therapy, from resolving heating below the threshold for damage to ensuring that the thermal exposure is sufficient within the target volume and protecting surrounding critical structures and to accurately predicting the extent of the ablated volume. In this paper, these validation studies will be reviewed. In addition, clinical studies that have shown this method feasible in human treatments will be overviewed.

Autoregressive moving average modeling for spectral parameter estimation from a multigradient echo chemical shift acquisition

The authors investigated the performance of the iterative Steiglitz-McBride (SM) algorithm on an autoregressive moving average (ARMA) model of signals from a fast, sparsely sampled, multiecho, chemical shift imaging (CSI) acquisition using simulation, phantom, ex vivo, and in vivo experiments with a focus on its potential usage in magnetic resonance (MR)-guided interventions. The ARMA signal model facilitated a rapid calculation of the chemical shift, apparent spin-spin relaxation time (T2*), and complex amplitudes of a multipeak system from a limited number of echoes (< or equal 16). Numerical simulations of one- and two-peak systems were used to assess the accuracy and uncertainty in the calculated spectral parameters as a function of acquisition and tissue parameters. The measured uncertainties from simulation were compared to the theoretical Cramer-Rao lower bound (CRLB) for the acquisition. Measurements made in phantoms were used to validate the T2* estimates and to validate uncertainty estimates made from the CRLB. We demonstrated application to real-time MR-guided interventions ex vivo by using the technique to monitor a percutaneous ethanol injection into a bovine liver and in vivo to monitor a laser-induced thermal therapy treatment in a canine brain. Simulation results showed that the chemical shift and amplitude uncertainties reached their respective CRLB at a signal-to-noise ratio (SNR) > or =5 for echo train lengths (ETLs) > or =4 using a fixed echo spacing of 3.3 ms. T2* estimates from the signal model possessed higher uncertainties but reached the CRLB at larger SNRs and/or ETLs. Highly accurate estimates for the chemical shift (<0.01 ppm) and amplitude (<1.0%) were obtained with > or =4 echoes and for T2*(<1.0%) with > or =7 echoes. We conclude that, over a reasonable range of SNR, the SM algorithm is a robust estimator of spectral parameters from fast CSI acquisitions that acquire < or =16 echoes for one- and two-peak systems. Preliminary ex vivo and in vivo experiments corroborated the results from simulation experiments and further indicate the potential of this technique for MR-guided interventional procedures with high spatiotemporal resolution approximately 1.6 x 1.6 x 4 mm3 in < or =5 s.

Dynamics of MRI-Guided thermal ablation of VX2 tumor in paraspinal muscle of rabbits

This study combines fast magnetic resonance imaging (MRI) and model simulation of tissue thermal ablation for monitoring and predicting the dynamics of lesion size for tumor destruction. In vivo experiments were conducted using radiofrequency (RF) thermal ablation in paraspinal muscle of rabbit with a VX2 tumor. Before ablation, turbo-spin echo (TSE) images visualized the 3-D tumor (necrotic core and tumor periphery) and surrounding normal tissue. MR gradient-recalled echo (GRE) phase and magnitude images were acquired repeatedly in 3.3 s at 30-s intervals during and after thermal ablation to follow tissue temperature distribution dynamics and lesion development in tumor and surrounding normal tissue. Final lesion sizes estimated from GRE magnitude, post-ablation TSE, and stained histologic images were compared. Model simulations of temperature distribution and lesion development dynamics closely corresponded to the experimental data from MR images in tumor and normal tissue. The combined use of MR image monitoring and model simulation has the potential for improving pretreatment planning and real-time prediction of lesion-size dynamics for guidance of thermal ablation of tumors.

The risk of exposure to diagnostic ultrasound in postnatal subjects: thermal effects

This review evaluates the thermal mechanism for ultrasound-induced biological effects in postnatal subjects. The focus is the evaluation of damage versus temperature increase. A view of ultrasound-induced temperature increase is presented, based on thermodynamic Arrhenius analyses. The hyperthermia and other literature revealed data that allowed for an estimate of a temperature increase threshold of tissue damage for very short exposure times. This evaluation yielded an exposure time extension of the 1997 American Institute of Ultrasound in Medicine Conclusions Regarding Heat statement (American Institute of Ultrasound in Medicine, Laurel, MD) to 0.1 second for nonfetal tissue, where, at this exposure time, the temperature increase threshold of tissue damage was estimated to be about 18 degrees C. The output display standard was also evaluated for soft tissue and bone cases, and it was concluded that the current thermal indices could be improved to reduce the deviations and scatter of computed maximum temperature rises.

MRI-guided Thermal Ablation Therapy: Model and Parameter Estimates to Predict Cell Death from MR Thermometry Images

Solid tumors and other pathologies can be treated using laser thermal ablation under interventional magnetic resonance imaging (iMRI) guidance. A model was developed to predict cell death from magnetic resonance (MR) thermometry measurements based on the temperature-time history, and validated using in vivo rabbit brain data. To align post-ablation T2-weighted spin-echo MR lesion images to gradient-echo MR images, from which temperature is derived, a registration method was used that aligned fiducials placed near the thermal lesion. The outer boundary of the hyperintense rim in the post-ablation MR lesion image was used as the boundary for cell death, as verified from histology. Model parameters were simultaneously estimated using an iterative optimization algorithm applied to every interesting voxel in 328 images from multiple experiments having various temperature histories. For a necrotic region of 766 voxels across all lesions, the model provided a voxel specificity and sensitivity of 98.1 and 78.5%, respectively. Mislabeled voxels were typically within one voxel from the segmented necrotic boundary with median distances of 0.77 and 0.22 mm for false positives (FP) and false negatives (FN), respectively. As compared to the critical temperature cell death model and the generalized Arrhenius model, our model typically predicted fewer FP and FN. This is good evidence that iMRI temperature maps can be used with our model to predict therapeutic regions in real-time during treatment.

MRI-guided laser thermal ablation: model to predict cell death from MR thermometry images for real-time therapy monitoring

Solid tumors and other pathologies can be treated using laser thermal ablation under interventional magnetic resonance imaging (iMRI) guidance. We developed a model to predict cell death from MR thermometry measurements and applied it to in vivo rabbit brain data. To align post-ablation T2-weighted spin-echo MR lesion images to gradient echo MR images, from which temperature is derived, we used a registration method that aligned fiducials placed near the thermal lesion. We used the outer boundary of the hyperintense rim in the post-ablation MR lesion image as the boundary for cell death, as verified from histology. Model parameters were simultaneously estimated using an iterative optimization algorithm applied to every interesting pixel in 328 images from multiple experiments having various temperature histories. For a necrotic region of 766 voxels across all lesions, the model gave a voxel specificity and sensitivity of 98.1% and 78.4%, respectively. Median distance between the segmented necrotic boundary and the mislabeled voxels was within one MR voxel. Furthermore, our model predicted fewer errors as compared to the critical temperature cell death model. This is good evidence that iMRI temperature maps can be used with our model to predict therapeutic regions in real-time.

Magnetic resonance imaging and model prediction for thermal ablation of tissue

PURPOSE

To monitor and predict tissue temperature distributions and lesion boundaries during thermal ablation by combining MRI and thermal modeling methods.

MATERIALS AND METHODS

Radiofrequency (RF) ablation was conducted in the paraspinal muscles of rabbits with MRI monitoring. A gradient-recalled echo (GRE) sequence via a 1.5T MRI system provided tissue temperature distribution from the phase images and lesion progression from changes in magnitude images. Post-ablation GRE estimates of lesion size were compared with post-ablation T2-weighted turbo-spin-echo (TSE) images and hematoxylin and eosin (H&E)-stained histological slices. A three-dimensional (3D) thermal model was used to simulate and predict tissue temperature and lesion size dynamics.

RESULTS

The lesion area estimated from repeated GRE images remained constant during the post-heating period when the temperature of the lesion boundary was less than a critical temperature. The final lesion areas estimated from multi-slice (M/S) GRE, TSE, and histological slices were not statistically different. The model-simulated tissue temperature distribution and lesion area closely corresponded to the GRE-based MR measurements throughout the imaging experiment.

CONCLUSION

For normal tissue in vivo, the dynamics of tissue temperature distribution and lesion size during RF thermal ablation can be 1) monitored with GRE phase and magnitude images, and 2) simulated for prediction with a thermal model.

Ultrasound-biophysics mechanisms

Ultrasonic biophysics is the study of mechanisms responsible for how ultrasound and biological materials interact. Ultrasound-induced bioeffect or risk studies focus on issues related to the effects of ultrasound on biological materials. On the other hand, when biological materials affect the ultrasonic wave, this can be viewed as the basis for diagnostic ultrasound. Thus, an understanding of the interaction of ultrasound with tissue provides the scientific basis for image production and risk assessment. Relative to the bioeffect or risk studies, that is, the biophysical mechanisms by which ultrasound affects biological materials, ultrasound-induced bioeffects are generally separated into thermal and non-thermal mechanisms. Ultrasonic dosimetry is concerned with the quantitative determination of ultrasonic energy interaction with biological materials. Whenever ultrasonic energy is propagated into an attenuating material such as tissue, the amplitude of the wave decreases with distance. This attenuation is due to either absorption or scattering. Absorption is a mechanism that represents that portion of ultrasonic wave that is converted into heat, and scattering can be thought of as that portion of the wave, which changes direction. Because the medium can absorb energy to produce heat, a temperature rise may occur as long as the rate of heat production is greater than the rate of heat removal. Current interest with thermally mediated ultrasound-induced bioeffects has focused on the thermal isoeffect concept. The non-thermal mechanism that has received the most attention is acoustically generated cavitation wherein ultrasonic energy by cavitation bubbles is concentrated. Acoustic cavitation, in a broad sense, refers to ultrasonically induced bubble activity occurring in a biological material that contains pre-existing gaseous inclusions. Cavitation-related mechanisms include radiation force, microstreaming, shock waves, free radicals, microjets and strain. It is more challenging to deduce the causes of mechanical effects in tissues that do not contain gas bodies. These ultrasonic biophysics mechanisms will be discussed in the context of diagnostic ultrasound exposure risk concerns.

Microbubble contrast agent with focused ultrasound to create brain lesions at low power levels: MR imaging and histologic study in rabbits

PURPOSE

To evaluate magnetic resonance (MR) imaging-based thermometry for predicting the onset and spatial extent of lesions produced by focused ultrasound combined with a microbubble contrast agent (Optison; GE Healthcare, Milwaukee, Wis) and to compare the resulting induced temperature increase and threshold for damage with those in studies performed without the agent.

MATERIALS AND METHODS

The experiments were approved by the animal care committee. Fifty-three locations in the brains of 15 rabbits were sonicated with various exposure parameters by using a 1.5-MHz focused ultrasound transducer. MR imaging was used to map the temperature rise and, along with light microscopy, to examine the lesions. Diameters of isotherms created from thermometry were compared with the resulting lesions by using Bland-Altman analysis and linear regression. The minimum acoustic power necessary for lesion creation was determined, and the apparent temperature threshold for damage was calculated with probit analysis. These thresholds were compared with prior work performed without the contrast agent. The heating induced with the microbubbles was compared with that in sonications performed without them by using a t test.

RESULTS

The MR imaging-mapped temperature distributions matched the shape of the lesions. The diameters of isotherms correlated well with diameters measured at contrast material-enhanced MR imaging (mean difference between measurements, 0.0 mm +/- 0.5; R = 0.93). The temperature increase with microbubbles was statistically larger (P < .01) than for sonications performed without microbubbles. In some locations (mostly continuous wave exposures), damage was observed along the ultrasound beam path. The time-averaged acoustic power damage threshold was reduced by 91% for 10-second exposures when compared with earlier studies performed without microbubbles. The probability of producing lesions was 50% at a temperature increase of 5.9 degrees C, 5.5 degrees C lower than was observed earlier without the agent.

CONCLUSION

MR imaging-based temperature measurements appeared to correlate with focused ultrasound-induced lesions in the brain when microbubbles were present, even though the temperature appeared to be below the threshold for thermal damage.

Uterine leiomyomas: MR imaging-based thermometry and thermal dosimetry during focused ultrasound thermal ablation

PURPOSE

To retrospectively evaluate magnetic resonance (MR) imaging-based thermometry and thermal dosimetry during focused ultrasound treatments of uterine leiomyomas (ie, fibroids).

MATERIALS AND METHODS

All patients gave written informed consent for the focused ultrasound treatments and the current HIPAA-compliant retrospective study, both of which were institutional review board approved. Thermometry performed during the treatments of 64 fibroids in 50 women (mean age, 46.6 years +/- 4.5 [standard deviation]) was used to create thermal dose maps. The areas that reached dose values of 240 and 18 equivalent minutes at 43 degrees C were compared with the nonperfused regions measured on contrast material-enhanced MR images by using the Bland-Altman method. Volume changes in treated fibroids after 6 months were compared with volume changes in nontreated fibroids and with MR-based thermal dose estimates.

RESULTS

While the thermal dose estimates were shown to have a clear relationship with resulting nonperfused regions, the nonperfused areas were, on average, larger than the dose estimates (means of 1.9 +/- 0.7 and 1.2 +/- 0.4 times as large for areas that reached 240- and 18-minute threshold dose values, respectively). Good correlation was observed for smaller treatment volumes at the lower dose threshold (mean ratio, 1.0 +/- 0.3), but for larger treatment volumes, the nonperfused region extended to locations within the fibroid that clearly were not heated. Variations in peak temperature increase were as large as a factor of two, both between patients and within individual treatments. On average, the fibroid volume reduction at 6 months increased as the ablated volume estimated by using the thermal dose increased.

CONCLUSION

Study results showed good correlation between thermal dose estimates and resulting nonperfused areas for smaller ablated volumes. For larger treatment volumes, nonperfused areas could extend within the fibroid to unheated areas.

MR imaging-controlled focused ultrasound ablation: a noninvasive image-guided surgery

The history of MR-guided FUS demonstrates the need for merging advanced therapy technology with advanced imaging. Without the ability of MR imaging to localize the tumor margins and without the temperature-sensitive imaging that provides the closed-loop control of energy deposition, this method is inadequate for most clinical applications. Given these limitations,high-intensity focused ultrasound initially appeared to have a narrow application area and was not able to compete with other surgical or ablation methods. Today, MR imaging-guided FUS has become a safe and effective means of performing probe-delivered thermal ablations and minimally invasive surgery. Moreover, it has the potential to replace treatments that use ionizing radiation such as radiosurgery and brachytherapy. Although the cost of integrating"big ticket" MR imaging systems with complex and expensive phased arrays is high, this expenditure will largely be offset by eliminating hospitalization and anesthesia and by reducing complications. In effect, an investment in this emerging technology will ultimately redound to the benefit of the health care delivery system and, most important, to the patient. The FUS system provides a safe, repeatable treatment approach for benign tumors (eg, uterine fibroid and breast fibroadenoma) that do not require an aggressive approach. MR-guided FUS can also be used for debulking cancerous tissue. It has already been tested as a breast cancer treatment; its application for other malignancies in the brain, liver, and prostate is under development. MR-guided FUS offers an attractive alternative to conventional surgery because it incorporates intraoperative MR imaging, which provides far more precise target definition than is possible with the surgeon's direct visualization of the lesion. MR-guided FUS is undeniably the most promising interventional MR imaging method in the field of image-guided therapy today. It is applicable not only in the thermal coagulative treatment of tumors but also in several other medical situations for which invasive surgery or radiation may not be treatment options. The use of FUS for treating vascular malformation or functional disorders of the brain is also exciting. It is uniquely applicable for image-guided therapy using targeted drug delivery methods and gene therapy. Further advances in this technology will no doubt improve energy deposition and reduce treatment times. In the near future, FUS will offer a viable alternative to conventional surgery and radiation therapy; in the longer-term, it may also enable a host of targeted treatment methods aimed at eradicating or arresting heretofore intractable diseases such as certain brain malignancies and forms of epilepsy.

MRI guided and monitored focused ultrasound thermal ablation methods: a review of progress

This paper reviews the current status in using magnetic resonance imaging (MRI) to guide and monitor thermal coagulation of tumours using focused ultrasound. The patient treatment procedure with a second generation phased array system will be described. Several clinical trials have found that patient treatments are feasible and that MRI thermometry allows noninvasive monitoring of clinical treatments. Overall, this emerging modality holds significant potential for non-invasive tumour treatment of both benign and malignant tumours.

MRI-guided interstitial ultrasound thermal therapy of the prostate: A feasibility study in the canine model

The feasibility of MRI-guided interstitial ultrasound thermal therapy of the prostate was evaluated in an in vivo canine prostate model. MRI compatible, multielement interstitial ultrasound applicators were developed using 1.5 mm diameter cylindrical piezoceramic transducers (7 to 8 MHz) sectored to provide 180 degrees of angular directional heating. Two in vivo experiments were performed in canine prostate. The first using two interstitial ultrasound applicators, the second using three ultrasound applicators in conjunction with rectal and urethral cooling. In both experiments, the applicators were inserted transperineally into the prostate with the energy directed ventrally, away from the rectum. Electrical power levels of 5-17 W per element (approximately 1.6-5.4 W acoustic output power) were applied for heating periods of 18 and 48 min. Phase-sensitive gradient-echo MR imaging was used to monitor the thermal treatment in real-time on a 0.5 T interventional MRI system. Contrast-enhanced T1-weighted images and vital-stained serial tissue sections were obtained to assess thermal damage and correlate to real-time thermal contour plots and calculated thermal doses. Results from these studies indicated a large volume of ablated (nonstained) tissue within the prostate, extending 1.2 to 2.0 cm from the applicators to the periphery of the gland, with the dorsal margin of coagulation well-defined by the applicator placement and directionality. The shape of the lesions correlated well to the hypointense regions visible in the contrast-enhanced T1-weighted images, and were also in good agreement with the contours of the 52 degrees C threshold temperature and t43 > 240 min. This study demonstrates the feasibility of using directional interstitial ultrasound in conjunction with MRI thermal imaging to monitor and possibly control thermal coagulation within a targeted tissue volume while potentially protecting surrounding tissue, such as rectum, from thermal damage.

Noninvasive thermal ablation of hepatocellular carcinoma by using magnetic resonance imaging-guided focused ultrasound

A number of minimally invasive methods have been tested for the thermal ablation of liver tumors as an alternative to surgical resection. The use of focused ultrasound transducers to ablate deep tumors offers the first completely noninvasive alternative to these techniques. By increasing the flexibility of this technology with modern phased-array transducer design and by combining it with magnetic resonance imaging for targeting and online guidance, a powerful tool results with the potential to offer treatment to a larger population of patients, to reduce trauma to the patient, and to reduce the cost of treatment. In this article, we review previous work with focused ultrasound in the liver and recent experimental results with magnetic resonance imaging guidance.

Interventional and intraoperative MR: review and update of techniques and clinical experience

The concept of interventional magnetic resonance imaging (MRI) is based on the integration of diagnostic and therapeutic procedures, favored by the combination of the excellent morphological and functional imaging characteristics of MRI. The spectrum of MRI-assisted interventions ranges from biopsies and intraoperative guidance to thermal ablation modalities and vascular interventions. The most relevant recently published experimental and clinical results are discussed. In the future, interventional MRI is expected to play an important role in interventional radiology, minimal invasive therapy and guidance of surgical procedures. However, the associated high costs require a careful evaluation of its potentials in order to ensure cost-effective medical care.

Thermal modeling of lesion growth with radiofrequency ablation devices

BACKGROUND

Temperature is a frequently used parameter to describe the predicted size of lesions computed by computational models. In many cases, however, temperature correlates poorly with lesion size. Although many studies have been conducted to characterize the relationship between time-temperature exposure of tissue heating to cell damage, to date these relationships have not been employed in a finite element model.

METHODS

We present an axisymmetric two-dimensional finite element model that calculates cell damage in tissues and compare lesion sizes using common tissue damage and iso-temperature contour definitions. The model accounts for both temperature-dependent changes in the electrical conductivity of tissue as well as tissue damage-dependent changes in local tissue perfusion. The data is validated using excised porcine liver tissues.

RESULTS

The data demonstrate the size of thermal lesions is grossly overestimated when calculated using traditional temperature isocontours of 42 degrees C and 47 degrees C. The computational model results predicted lesion dimensions that were within 5% of the experimental measurements.

CONCLUSION

When modeling radiofrequency ablation problems, temperature isotherms may not be representative of actual tissue damage patterns.

MRI investigation of the threshold for thermally induced blood-brain barrier disruption and brain tissue damage in the rabbit brain

The ability of MRI-derived thermometry to predict thermally induced tissue changes in the brain was tested, and the thermal thresholds for blood-brain barrier (BBB) disruption and brain tissue damage were estimated. In addition, the ability of standard MRI to detect threshold-level effects was confirmed. These safety thresholds are being investigated to provide guidelines for clinical thermal ablation studies in the brain. MRI-monitored focused ultrasound heating was delivered to 63 locations in 26 rabbits. Tissue changes were detected in T(2)-weighted imaging and T(1)-weighted imaging (with and without contrast) and with light microscopy. The probability for tissue damage as a function of the accumulated thermal dose, the peak temperature achieved, the applied acoustic energy, and the peak acoustic power was estimated with probit regression. The discriminative abilities of these parameters were compared using the areas under the receiver operator characteristic (ROC) curves. In MRI, BBB disruption was observed in contrast-enhanced T(1)-weighted imaging shortly after the ultrasound exposures, sometimes accompanied by changes in T(2)-weighted imaging. Two days later, changes in T(2)-weighted imaging were observed, sometimes accompanied by changes in T(1)-weighted imaging. In histology, tissue damage was seen at every location where MRI changes were observed, ranging from small (diameter <1.0 mm) areas of tissue necrosis to severe vascular damage and associated hemorrhagic infarct. In one location, small (diameter: 0.8 mm) damage was not detected in MRI. The thermal dose and peak temperature thresholds were between 12.3-40.1 equivalent min at 43 degrees C and 48.0-50.8 degrees C, respectively, and values of 17.5 equivalent min at 43 degrees C and 48.4 degrees C were estimated to result in tissue damage with 50% probability. Thermal dose and peak temperature were significantly better predictors than the applied acoustic energy and peak acoustic power (P < 0.01). BBB disruption was always accompanied by tissue damage. The temperature information was better than the applied acoustic power or energy for predicting the damage than the ultrasound parameters. MRI was sensitive in detecting threshold-level damage.