MR thermometry for monitoring tumor ablation

Patient-specific temperature distribution prediction in laser interstitial thermal therapy: single-irradiation data-driven method

Laser interstitial thermal therapy (LITT) is popular for treating brain tumours and epilepsy. The strict control of tissue thermal damage extent is crucial for LITT. Temperature prediction is useful for predicting thermal damage extent. Accurately predictingin vivobrain tissue temperature is challenging due to the temperature dependence and the individual variations in tissue properties. Considering these factors is essential for improving the temperature prediction accuracy.Objective. To present a method for predicting patient-specific tissue temperature distribution within a target lesion area in the brain during LITT.Approach. A magnetic resonance temperature imaging (MRTI) data-driven estimation model was constructed and combined with a modified Pennes bioheat transfer equation (PBHE) to predict patient-specific temperature distribution. In the PBHE for temperature prediction, the individual specificity and temperature dependence of thermal tissue properties and blood perfusion, as well as the individual specificity of optical tissue properties were considered. Only MRTI data during one laser irradiation were required in the method. This enables the prediction of patient-specific temperature distribution and the resulting thermal damage region for subsequent ablations.Main results. Patient-specific temperature prediction was evaluated based on clinical data acquired during LITT in the brain, using intraoperative MRTI data as the reference standard. Our method significantly improved the prediction performance of temperature distribution and thermal damage region. The average root mean square error was decreased by 69.54%, the average intraclass correlation coefficient was increased by 37.5%, the average Dice similarity coefficient was increased by 43.14% for thermal damage region prediction.Significance. The proposed method can predict temperature distribution and thermal damage region at an individual patient level during LITT, providing a promising approach to assist in patient-specific treatment planning for LITT in the brain.

Biomechanical testing of ex vivo porcine tendons following high intensity focused ultrasound thermal ablation

INTRODUCTION

Magnetic resonance-guided focused ultrasound (MRgFUS) has been demonstrated to be able to thermally ablate tendons with the aim to non-invasively disrupt tendon contractures in the clinical setting. However, the biomechanical changes of tendons permitting this disrupting is poorly understood. We aim to obtain a dose-dependent biomechanical response of tendons following magnetic resonance-guided focused ultrasound (MRgFUS) thermal ablation.

METHODS

Ex vivo porcine tendons (n = 72) were embedded in an agar phantom and randomly assigned to 12 groups based on MRgFUS treatment. The treatment time was 10, 20, or 30s, and the applied acoustic power was 25, 50, 75, or 100W. Following each MRgFUS treatment, tendons underwent biomechanical tensile testing on an Instron machine, which calculated stress-strain curves during tendon elongation. Rupture rate, maximum treatment temperature, Young's modulus and ultimate strength were analyzed for each treatment energy.

RESULTS

The study revealed a dose-dependent response, with tendons rupturing in over 50% of cases when energy delivery exceeded 1000J and 100% disruption at energy levels beyond 2000J. The achieved temperatures during MRgFUS were directly proportional to energy delivery. The highest recorded temperature was 56.8°C ± 9.34 (3000J), while the lowest recorded temperate was 18.6°C ± 0.6 (control). The Young's modulus was highest in the control group (47.3 MPa ± 6.5) and lowest in the 3000J group (13.2 MPa ± 5.9). There was no statistically significant difference in ultimate strength between treatment groups.

CONCLUSION

This study establishes crucial thresholds for reliable and repeatable disruption of tendons, laying the groundwork for future in vivo optimization. The findings prompt further exploration of MRgFUS as a non-invasive modality for tendon disruption, offering hope for improved outcomes in patients with musculotendinous contractures.

A review on radiofrequency, laser, and microwave ablations and their thermal monitoring through fiber Bragg gratings

Thermal ablation of tumors aims to apply extreme temperatures inside the target tissue to achieve substantial tumor destruction in a minimally invasive manner. Several techniques are comprised, classified according to the type of energy source. However, the lack of treatment selectivity still needs to be addressed, potentially causing two risks: i) incomplete tumor destruction and recurrence, or conversely, ii) damage of the surrounding healthy tissue. Therefore, the research herein reviewed seeks to develop sensing systems based on fiber Bragg gratings (FBGs) for thermal monitoring inside the lesion during radiofrequency, laser, and microwave ablation. This review shows that, mainly thanks to multiplexing and minimal invasiveness, FBGs provide an optimal sensing solution. Their temperature measurements are the feedback to control the ablation process and allow to investigate different treatments, compare their outcomes, and quantify the impact of factors such as proximity to thermal probe and blood vessels, perfusion, and tissue type.

Noninvasive magnetic resonance-guided focused ultrasound for tendon disruption: an in vivo Animal study

PURPOSE

Surgical resection of the tendon is an effective treatment for severe contracture. Magnetic Resonance-guided Focused Ultrasound (MRgFUS) is a non-invasive ultrasonic therapy which produces a focal increase in temperature, subsequent tissue ablation and disruption. We evaluated MRgFUS as a clinically translatable treatment modality to non-invasively disrupt in vivo porcine tendons.

MATERIAL AND METHODS

In vivo Achilles tendons (n = 28) from 15-20kg Yorkshire pigs (n = 16) were randomly assigned to 4 treatment groups of 600, 900, 1200 and 1500 J. Pretreatment range of motion (ROM) of the ankle joint was measured with the animal under general anesthesia. Following MRgFUS treatment, success of tendon rupture, ROM increase, temperature, thermal dosage, skin burn, and histology analyses were performed.

RESULTS

Rupture success was found to be 29%, 86%, 100% and 100% for treatment energies of 600, 900, 1200 and 1500 J respectfully. ROM difference at 90° flexion showed a statistically significant change in ROM between 900 J and 1200 J from 16° to 27°. There was no statistical significance between other groups, but there was an increase in ROM as more energy was delivered in the treatment. For each of the respective treatment groups, the maximal temperatures were 58.4 °C, 63.3 °C, 67.6 °C, and 69.9 °C. The average areas of thermal dose measured were 24.3mm2, 53.2mm2, 77.8mm2 and 91.6mm2. The average areas of skin necrosis were 5.4mm2, 21.8mm2, 37.2mm2, and 91.4mm2. Histologic analysis confirmed tissue ablation and structural collagen fiber disruption.

CONCLUSIONS

This study demonstrated that MRgFUS is able to disrupt porcine tendons in vivo without skin incisions.

In Vivo Thermal Ablation of Deep Intrahepatic Targets Using a Super-Convergent MRgHIFU Applicator and a Pseudo-Tumor Model

BACKGROUND

HIFU ablation of liver malignancies is particularly challenging due to respiratory motion, high tissue perfusion and the presence of the rib cage. Based on our previous development of a super-convergent phased-array transducer, we aimed to further investigate, in vivo, its applicability to deep intrahepatic targets.

METHODS

In a series of six pigs, a pseudo-tumor model was used as target, visible both on intra-operatory MRI and post-mortem gross pathology. The transcostal MRgHIFU ablation was prescribed coplanar with the pseudo-tumor, either axial or sagittal, but deliberately shifted 7 to 18 mm to the side. No specific means of protection of the ribs were implemented. Post-treatment MRI follow-up was performed at D7, followed by animal necropsy and gross pathology of the liver.

RESULTS

The pseudo-tumor was clearly identified on T1w MR imaging and subsequently allowed the MRgHIFU planning. The peak temperature at the focal point ranged from 58-87 °C. Gross pathology confirmed the presence of the pseudo-tumor and the well-delineated MRgHIFU ablation at the expected locations.

CONCLUSIONS

The specific design of the transducer enabled a reliable workflow. It demonstrated a good safety profile for in vivo transcostal MRgHIFU ablation of deep-liver targets, graded as challenging for standard surgery.

A feasibility study of wireless inductively coupled surface coil for MR-guided high-intensity focused ultrasound ablation of rodents on clinical MRI systems

Recently, to conduct preclinical imaging research on clinical MRI systems has become an attractive alternative to researchers due to its wide availability, cost, and translational application to clinical human studies when compared to dedicated small animal, high-field preclinical MRI. However, insufficient signal-to-noise ratio (SNR) significantly degrades the applicability of those applications which require high SNR, e.g. magnetic resonance guided high-intensity focused ultrasound (MRgHIFU) treatment. This study introduces a wireless inductively coupled surface (WICS) coil design used on a clinical 3 T MRI system for MRgHIFU ablation. To evaluate the SNR improvement and temperature accuracy of WICS coil, the ex vivo experiments were performed on the pork tenderloins (n = 7) and the hind legs of deceased Sprague-Dawley rats (n = 5). To demonstrate the feasibility, the in vivo experiments were performed on the hind leg of Sprague-Dawley rat (n = 1). For all experiments, temperature measurements were performed before and during HIFU ablation. Temperature curves with and without WICS coil were compared to evaluate the temperature precision in ex vivo experiments. The use of WICS coil improves the temperature accuracy from 0.85 to 0.14 °C, demonstrating the feasibility of performing small animal MRgHIFU experiments using clinical 3 T MRI system with WICS coil.

Low-Power Magnetic Resonance-Guided Focused Ultrasound Tumor Ablation upon Controlled Accumulation of Magnetic Nanoparticles by Cascade-Activated DNA Cross-Linkers

Magnetic resonance-guided focused ultrasound (MRgFUS) is a promising non-invasive surgical technique with spatial specificity and minimal off-target effects. Despite the expanding clinical applications, the major obstacles associated with MRgFUS still lie in low magnetic resonance imaging (MRI) sensitivity and safety issues. High ultrasound power is required to resist the energy attenuation during the delivery to the tumor site and may cause damage to the surrounding healthy tissues. Herein, a surface modification strategy is developed to simultaneously strengthen MRI and ultrasound ablation of MRgFUS by prolonging Fe₃O₄ nanoparticles' blood circulation and tumor-environment-triggered accumulation and retention at the tumor site. Specifically, reactive oxygen species-labile methoxy polyethylene glycol and pH-responsive DNA cross-linkers are modified on the surface of Fe₃O₄ nanoparticles, which can transform nanoparticles into aggregations through the cascade responsive reactions at the tumor site. Notably, DNA is selected as the pH-responsive cross-linker because of its superior biocompatibility as well as the fast and sensitive response to the weak acidity of 6.5-6.8, corresponding to the extracellular pH of tumor tissues. Due to the significantly enhanced delivery and retention amount of Fe₃O₄ nanoparticles at the tumor site, the MRI sensitivity was enhanced by 1.7-fold. In addition, the ultrasound power was lowered by 35% to reach a sufficient thermal ablation effect. Overall, this investigation demonstrates a feasible resolution to promote the MRgFUS treatment by enhancing the therapeutic efficacy and reducing the side effects, which will be helpful to guide the clinical practice in the future.

A Novel Concept of a Phased-Array HIFU Transducer Optimized for MR-Guided Hepatic Ablation: Embodiment and First In-Vivo Studies

Purpose

High-intensity focused ultrasound (HIFU) is challenging in the liver due to the respiratory motion and risks of near-/far-field burns, particularly on the ribs. We implemented a novel design of a HIFU phased-array transducer, dedicated to transcostal hepatic thermo-ablation. Due to its large acoustic window and strong focusing, the transducer should perform safely for this application.

Material and Methods

The new HIFU transducer is composed of 256 elements distributed on 5 concentric segments of a specific radius (either 100, 111, or 125 mm). It has been optimally shaped to fit the abdominal wall. The shape and size of the acoustic elements were optimized for the largest emitting surface and the lowest symmetry. Calibration tests have been conducted on tissue-mimicking gels under 3-T magnetic resonance (MR) guidance. In-vivo MR-guided HIFU treatment was conducted in two pigs, aiming to create thermal ablation deep in the liver without significant side effects. Imaging follow-up was performed at D0 and D7. Sacrifice and post-mortem macroscopic examination occurred at D7, with the ablated tissue being fixed for pathology.

Results

The device showed −3-dB focusing capacities in a volume of 27 × 46 × 50 mm³ as compared with the numerical simulation volume of 18 × 48 × 60 mm³. The shape of the focal area was in millimeter-range agreement with the numerical simulations. No interference was detected between the HIFU sonication and the MR acquisition. In vivo, the temperature elevation in perivascular liver parenchyma reached 28°C above physiological temperature, within one breath-hold. The lesion was visible on Gd contrast-enhanced MRI sequences and post-mortem examination. The non-perfused volume was found in pig #1 and pig #2 of 8/11, 6/8, and 7/7 mm along the LR, AP, and HF directions, respectively. No rib burns or other near-field side effects were visually observed on post-mortem gross examination. High-resolution contrast-enhanced 3D MRI indicated a minor lesion on the sternum.

Conclusion

The performance of this new HIFU transducer has been demonstrated in vitro and in vivo. The transducer meets the requirement to perform thermal lesions in deep tissues, without the need for rib-sparing means.

Biopsy Artifact in Laser Interstitial Thermal Therapy: A Technical Note

Purpose

The safety and effectiveness of laser interstitial thermal therapy (LITT) relies critically on the ability to continuously monitor the ablation based on real-time temperature mapping using magnetic resonance thermometry (MRT). This technique uses gradient recalled echo (GRE) sequences that are especially sensitive to susceptibility effects from air and blood. LITT for brain tumors is often preceded by a biopsy and is anecdotally associated with artifact during ablation. Thus, we reviewed our experience and describe the qualitative signal dropout that can interfere with ablation.

Methods

We retrospectively reviewed all LITT cases performed in our intraoperative MRI suite for tumors between 2017 and 2020. We identified a total of 17 LITT cases. Cases were reviewed for age, sex, pathology, presence of artifact, operative technique, and presence of blood/air on post-operative scans.

Results

We identified six cases that were preceded by biopsy, all six had artifact present during ablation, and all six were noted to have air/blood on their post-operative MRI or CT scans. In two of those cases, the artifactual signal dropout qualitatively interfered with thermal damage thresholds at the borders of the tumor. There was no artifact in the 11 non-biopsy cases and no obvious blood or air was noted on the post-ablation scans.

Conclusion

Additional consideration should be given to pre-LITT biopsies. The presence of air/blood caused an artifactual signal dropout effect in cases with biopsy that was severe enough to interfere with ablation in a significant number of those cases. Additional studies are needed to identify modifying strategies.

Interleaved water and fat MR thermometry for monitoring high intensity focused ultrasound ablation of bone lesions

PURPOSE

To demonstrate that interleaved MR thermometry can monitor temperature in water and fat with adequate temporal resolution. This is relevant for high intensity focused uUltrasounds (HIFU) treatment of bone lesions, which are often found near aqueous tissues, as muscle, or embedded in adipose tissues, as subcutaneous fat and bone marrow.

METHODS

Proton resonance frequency shift (PRFS)-based thermometry scans and T1 -based 2D variable flip angle (2D-VFA) thermometry scans were acquired alternatingly over time. Temperature in water was monitored using PRFS thermometry, and in fat by 2D-VFA thermometry with slice profile effect correction. The feasibility of interleaved water/fat temperature monitoring was studied ex vivo in porcine bone during MR-HIFU sonication. Precision and stability of measurements in vivo were evaluated in a healthy volunteer under non-heating conditions.

RESULTS

The method allowed observing temperature change over time in muscle and fat, including bone marrow, during MR-HIFU sonication, with a temporal resolution of 6.1 s. In vivo, the apparent temperature change was stable on the time scale of the experiment: In 7 min the systematic drift was <0.042°C/min in muscle (PRFS after drift correction) and <0.096°C/min in bone marrow (2D-VFA). The SD of the temperature change averaged over time was 0.98°C (PRFS) and 2.7°C (2D-VFA).

CONCLUSIONS

Interleaved MR thermometry allows temperature measurements in water and fat with a temporal resolution high enough for monitoring HIFU ablation. Specifically, combined fat and water thermometry provides uninterrupted information on temperature changes in tissue close to the bone cortex.

Photoacoustic computed tomography of mechanical HIFU-induced vascular injury

Mechanical high-intensity focused ultrasound (HIFU) has been used for cancer treatment and drug delivery. Existing monitoring methods for mechanical HIFU therapies such as MRI and ultrasound imaging often suffer from high cost, poor spatial-temporal resolution, and/or low sensitivity to tissue's hemodynamic changes. Evaluating vascular injury during mechanical HIFU treatment, therefore, remains challenging. Photoacoustic computed tomography (PACT) is a promising tool to meet this need. Intrinsically sensitive to optical absorption, PACT provides high-resolution imaging of blood vessels using hemoglobin as the endogenous contrast. In this study, we have developed an integrated HIFU-PACT system for detecting vascular rupture in mechanical HIFU treatment. We have demonstrated singular value decomposition for enhancing hemorrhage detection. We have validated the HIFU-PACT performance on phantoms and in vivo animal tumor models. We expect that PACT-HIFU will find practical applications in oncology research using small animal models.

Rapid Impedance Spectroscopy for Monitoring Tissue Impedance, Temperature, and Treatment Outcome During Electroporation-Based Therapies

OBJECTIVE

Electroporation-based therapies (EBTs) employ high voltage pulsed electric fields (PEFs) to permeabilize tumor tissue; this results in changes in electrical properties detectable using electrical impedance spectroscopy (EIS). Currently, commercial potentiostats for EIS are limited by impedance spectrum acquisition time (  ∼ 10 s); this timeframe is much larger than pulse periods used with EBTs (  ∼ 1 s). In this study, we utilize rapid EIS techniques to develop a methodology for characterizing electroporation (EP) and thermal effects associated with high-frequency irreversible EP (H-FIRE) in real-time by monitoring inter-burst impedance changes.

METHODS

A charge-balanced, bipolar rectangular chirp signal is proposed for rapid EIS. Validation of rapid EIS measurements against a commercial potentiostat was conducted in potato tissue using flat-plate electrodes and thereafter for the measurement of impedance changes throughout IRE treatment. Flat-plate electrodes were then utilized to uniformly heat potato tissue; throughout high-voltage H-FIRE treatment, low-voltage inter-burst impedance measurements were used to continually monitor impedance change and to identify a frequency at which thermal effects are delineated from EP effects.

RESULTS

Inter-burst impedance measurements (1.8 kHz - 4.93 MHz) were accomplished at 216 discrete frequencies. Impedance measurements at frequencies above  ∼ 1 MHz served to delineate thermal and EP effects in measured impedance.

CONCLUSION

We demonstrate rapid-capture ( 1 s) EIS which enables monitoring of inter-burst impedance in real-time. For the first time, we show impedance analysis at high frequencies can delineate thermal effects from EP effects in measured impedance.

SIGNIFICANCE

The proposed waveform demonstrates the potential to perform inter-burst EIS using PEFs compatible with existing pulse generator topologies.

Techniques for Temperature Monitoring of Myocardial Tissue Undergoing Radiofrequency Ablation Treatments: An Overview

Cardiac radiofrequency ablation (RFA) has received substantial attention for the treatment of multiple arrhythmias. In this scenario, there is an ever-growing demand for monitoring the temperature trend inside the tissue as it may allow an accurate control of the treatment effects, with a consequent improvement of the clinical outcomes. There are many methods for monitoring temperature in tissues undergoing RFA, which can be divided into invasive and non-invasive. This paper aims to provide an overview of the currently available techniques for temperature detection in this clinical scenario. Firstly, we describe the heat generation during RFA, then we report the principle of work of the most popular thermometric techniques and their features. Finally, we introduce their main applications in the field of cardiac RFA to explore the applicability in clinical settings of each method.

Feature-based automated segmentation of ablation zones by fuzzy c-mean clustering during low-dose computed tomography

PURPOSE

Intra-procedural monitoring and post-procedural follow-up is necessary for a successful ablation treatment. An imaging technique which can assess the ablation geometry accurately is beneficial to monitor and evaluate treatment. In this study, we developed an automated ablation segmentation technique for serial low-dose, noisy ablation computed tomography (CT) or contrast-enhanced CT (CECT).

METHODS

Low-dose, noisy temporal CT and CECT volumes were acquired during microwave ablation on normal porcine liver (four with non-contrast CT and eight with CECT). Highly constrained backprojection (HYPR) processing was used to recover ablation zone information compromised by low-dose noise. First-order statistic features and normalized fractional Brownian features (NBF) were used to segment ablation zones by fuzzy c-mean clustering. After clustering, the segmented ablation zone was refined by cyclic morphological processing. Automatic and manual segmentations were compared to gross pathology with Dice's coefficient (morphological similarity), while cross-sectional dimensions were compared by percent difference.

RESULTS

Automatic and manual segmentations of the ablation zone were very similar to gross pathology (Dice Coefficients: Auto.-Path. = 0.84 ± 0.02; Manu.-Path. = 0.76 ± 0.03, P = 0.11). The differences in ablation area, major diameter and minor diameter were 17.9 ± 3.2%, 11.1 ± 3.2% and 16.2 ± 3.4%, respectively, when comparing automatic segmentation to gross pathology, which were lower than the differences of 32.9 ± 16.8%, 13.0 ± 9.8% and 21.8 ± 5.8% when comparing manual segmentation to gross pathology. Manual segmentations tended to overestimate gross pathology when ablation area was less than 15 cm2 , but the automated segmentation tended to underestimate gross pathology when ablation zone is larger than 20 cm2 .

CONCLUSION

Fuzzy c-means clustering may be used to aid automatic segmentation of ablation zones without prior information or user input, making serial CT/CECT has more potential to assess treatments intra-procedurally.

Fiber Bragg Grating Sensors for Performance Evaluation of Fast Magnetic Resonance Thermometry on Synthetic Phantom

The increasing recognition of minimally invasive thermal treatment of tumors motivate the development of accurate thermometry approaches for guaranteeing the therapeutic efficacy and safety. Magnetic Resonance Thermometry Imaging (MRTI) is nowadays considered the gold-standard in thermometry for tumor thermal therapy, and assessment of its performances is required for clinical applications. This study evaluates the accuracy of fast MRTI on a synthetic phantom, using dense ultra-short Fiber Bragg Grating (FBG) array, as a reference. Fast MRTI is achieved with a multi-slice gradient-echo echo-planar imaging (GRE-EPI) sequence, allowing monitoring the temperature increase induced with a 980 nm laser source. The temperature distributions measured with 1 mm-spatial resolution with both FBGs and MRTI were compared. The root mean squared error (RMSE) value obtained by comparing temperature profiles showed a maximum error of 1.2 °C. The Bland-Altman analysis revealed a mean of difference of 0.1 °C and limits of agreement 1.5/−1.3 °C. FBG sensors allowed to extensively assess the performances of the GRE-EPI sequence, in addition to the information on the MRTI precision estimated by considering the signal-to-noise ratio of the images (0.4 °C). Overall, the results obtained for the GRE-EPI fully satisfy the accuracy (~2 °C) required for proper temperature monitoring during thermal therapies.

Systematic review of pre-clinical and clinical devices for magnetic resonance-guided radiofrequency hyperthermia

Clinical trials have demonstrated the therapeutic benefits of adding radiofrequency (RF) hyperthermia (HT) as an adjuvant to radio- and chemotherapy. However, maximum utilization of these benefits is hampered by the current inability to maintain the temperature within the desired range. RF HT treatment quality is usually monitored by invasive temperature sensors, which provide limited data sampling and are prone to infection risks. Magnetic resonance (MR) temperature imaging has been developed to overcome these hurdles by allowing noninvasive 3D temperature monitoring in the target and normal tissues. To exploit this feature, several approaches for inserting the RF heating devices into the MR scanner have been proposed over the years. In this review, we summarize the status quo in MR-guided RF HT devices and analyze trends in these hybrid hardware configurations. In addition, we discuss the various approaches, extract best practices and identify gaps regarding the experimental validation procedures for MR - RF HT, aimed at converging to a common standard in this process.

Temperature Dependence of Tissue Thermal Parameters Should Be Considered in the Thermal Lesion Prediction in High-Intensity Focused Ultrasound Surgery

This study considers the temperature-dependent thermal parameters (specific heat capacity, thermal diffusivity and thermal conductivity) used when predicting the temperature rise of tissue exposed to high-intensity focused ultrasound (HIFU). Numerical analysis was performed using the Khokhlov-Zabolotskaya-Kuznetsov equation coupled with a bioheat transfer function. The thermal parameters were set as the functions of temperature using experimental data. The results revealed that, for liver tissue exposed to HIFU with a focal intensity of 3000 W/cm² for 10 s, the predicted focal temperature rise was 23% lower and the thermal lesion area 41% smaller than those predicted without considering the temperature dependence. The prediction was validated by experimental observations on thermal lesions visualized in a tissue-mimicking phantom. The present results suggest that temperature-dependent thermal parameters should be considered in the prediction of HIFU-induced temperature rise to avoid lowering ultrasonic output settings for HIFU surgery. The aim of the present study was to investigate how significantly the temperature dependence of the thermal parameters affects the thermal dose imposed on the tissue by a typical clinical HIFU device. A numerical simulation was performed using a thermo-acoustic algorithm coupling the non-linear Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation (Meaney et al. 1998; Filonenko and Khokhlova 2001) and a bio-heat transfer (BHT) equation (Pennes 1948). Thermal parameters of liver tissue were modeled in the present study as functions of temperature and were incorporated into the BHT equation to compensate for the variations in thermal parameters with temperature. Experimental validation was achieved by comparing the predictions with the thermal lesions formed in the tissue-mimicking phantoms.

A thermometry software tool for monitoring laser-induced interstitial thermotherapy

The purpose of this study was to develop a thermometry software tool for temperature monitoring during laser-induced interstitial thermotherapy (LITT). C++ programming language and several libraries including DICOM Toolkit, Grassroots DICOM library, Insight Segmentation and Registration Toolkit, Visualization Toolkit and Quasar Toolkit were used. The software's graphical user interface creates windows displaying the temperature map and the coagulation extent in the tissue, determined by the magnetic resonance imaging (MRI) thermometry with the echo planar imaging sequence and a numerical simulation based on the radiation and heat transfer in biological tissues, respectively. The software was evaluated applying the MRI-guided LITT to ex vivo pig liver and simultaneously measuring the temperature through a fiber-optic thermometer as reference. Using the software, the temperature distribution determined by the MRI method was compared with the coagulation extent simulation. An agreement was shown between the MRI temperature map and the simulated coagulation extent. Furthermore, the MRI-based and simulated temperatures agreed with the measured one - a correlation coefficient of 0.9993 and 0.9996 was obtained, respectively. The precision of the MRI temperature amounted to 2.4°C. In conclusion, the software tool developed in the present study can be applied for monitoring and controlling the LITT procedure in ex vivo tissues.

Toward Image Data-Driven Predictive Modeling for Guiding Thermal Ablative Therapy

OBJECTIVE

Accurate prospective modeling of microwave ablation (MWA) procedures can provide powerful planning and navigational information to physicians. However, patient-specific tissue properties are generally unavailable and can vary based on factors such as relative perfusion and state of disease. Therefore, a need exists for modeling frameworks that account for variations in tissue properties.

METHODS

In this study, we establish an inverse modeling approach to reconstruct a set of tissue properties that best fit the model-predicted and observed ablation zone extents in a series of phantoms of varying fat content. We then create a model of these tissue properties as a function of fat content and perform a comprehensive leave-one-out evaluation of the predictive property model. Furthermore, we validate the inverse-model predictions in a separate series of phantoms that include co-recorded temperature data.

RESULTS

This model-based approach yielded thermal profiles in close agreement with experimental measurements in the series of validation phantoms (average root-mean-square error of 4.8 °C). The model-predicted ablation zones showed compelling overlap with observed ablations in both the series of validation phantoms (93.4 ± 2.2%) and the leave-one-out cross validation study (86.6 ± 5.3%). These results demonstrate an average improvement of 17.3% in predicted ablation zone overlap when comparing the presented property-model to properties derived from phantom component volume fractions.

CONCLUSION

These results demonstrate accurate model-predicted ablation estimates based on image-driven determination of tissue properties.

SIGNIFICANCE

The work demonstrates, as a proof-of-concept, that physical modeling parameters can be linked with quantitative medical imaging to improve the utility of predictive procedural modeling for MWA.

Percutaneous MR-guided prostate cancer cryoablation technical updates and literature review

Prostate cancer (PCa) is the most common malignant tumor in males. The benefits in terms of overall reduction in specific mortality due to the widespread use of Prostate Specific Antigen (PSA) screening and the advancements in the curative treatments (radical prostatectomy or radiotherapy) appear to have reached a plateau. There remains, however, the questions of overdiagnosis and overtreatment of such patients. Currently, the main challenge in the treatment of patients with clinically organ-confined PCa is to offer an oncologically efficient treatment with as little morbidity as possible. Amongst the arising novel curative techniques for PCa, cryoablation (CA) is the most established one, which is also included in the NICE and AUA guidelines. CA is commonly performed under ultrasound guidance with the inherent limitations associated with this technique. The recent advancements in MRI have significantly improved the accuracy of detecting and characterizing a clinically significant PCa. This, alongside the development of wide bore interventional MR scanners, has opened the pathway for in bore PCa treatment. Under MRI guidance, PCa CA can be used either as a standard whole gland treatment or as a tumor targeted one. With MR-fluoroscopy, needle guidance capability, multiplanar and real-time visualization of the iceball, MRI eliminates the inherent limitations of ultrasound guidance and can potentially lead to a lower rate of local complications. The aim of this review article is to provide an overview about PCa CA with a more specific insight on MR guided PCa CA; the limitations, challenges and applications of this novel technique will be discussed.

Can we ablate liver lesions close to large portal and hepatic veins with MR-guided HIFU? An experimental study in a porcine model

OBJECTIVES

Invasive treatment of tumors adjacent to large hepatic vessels is a continuous clinical challenge. The primary aim of this study was to examine the feasibility of ablating liver tissue adjacent to large hepatic and portal veins with magnetic resonance imaging-guided high-intensity focused ultrasound (MRgHIFU). The secondary aim was to compare sonication data for ablations performed adjacent to hepatic veins (HV) versus portal veins (PV).

MATERIALS AND METHODS

MRgHIFU ablations were performed in six male land swine under general anesthesia. Ablation cells of either 4 or 8 mm diameter were planned in clusters (two/animal) adjacent either to HV (n = 6) or to PV (n = 6), with diameter ≥ 5 mm. Ablations were made using 200 W and 1.2 MHz. Post-procedure evaluation was made on contrast-enhanced MRI (T1w CE-MRI), histopathology, and ablation data from the HIFU system.

RESULTS

A total of 153 ablations in 81 cells and 12 clusters were performed. There were visible lesions with non-perfused volumes in all animals on T1w CE-MRI images. Histopathology showed hemorrhage and necrosis in all 12 clusters, with a median shortest distance to vessel wall of 0.4 mm (range 0-2.7 mm). Edema and endothelial swelling were observed without vessel wall rupture. In 8-mm ablations (n = 125), heat sink was detected more often for HV (43%) than for PV (19%; p = 0.04).

CONCLUSIONS

Ablations yielding coagulative necrosis of liver tissue can be performed adjacent to large hepatic vessels while keeping the vessel walls intact. This indicates that perivascular tumor ablation in the liver is feasible using MRgHIFU.

KEY POINTS

• High-intensity focused ultrasound ablation is a non-invasive treatment modality that can be used for treatment of liver tumors. • This study shows that ablations of liver tissue can be performed adjacent to large hepatic vessels in an experimental setting. • Liver tumors close to large vessels can potentially be treated using this modality.

Weighting the Passive Acoustic Mapping Technique With the Phase Coherence Factor for Passive Ultrasound Imaging of Ultrasound-Induced Cavitation

Ultrasound (US) cavitation is currently being explored for low-invasive therapy techniques applied to a wide panel of pathologies. Because of the random behavior of cavitation, a real-time spatial monitoring system may be required. For this purpose, the US passive imaging techniques have been recently investigated. In particular, the passive acoustic mapping (PAM) beamforming method enables the reconstruction of cavitation activity maps by beamforming acoustic signals passively recorded by an array transducer. In this paper, an optimized version of PAM, PAM weighted with a phase coherence factor (PAM-PCF), is considered. A general validation process is developed including simulations on a point source and experiments on a wire. Furthermore, using a focused regulated US-induced cavitation generator, reproducible cavitation experiments are conducted in water and in agar gel. The spatial behavior of a bubble cavitation cloud is determined using the PAM-PCF beamforming method to localize the focal cavitation point in two perpendicular imaging planes.

Focused ultrasound for the treatment of bone metastases: effectiveness and feasibility

Background

To evaluate the effectiveness and feasibility of high-intensity focused ultrasound (HIFU) for the treatment of bone metastases.

Methods

A single-center prospective study was made involving 17 consecutive patients with symptomatic bone metastases. Patients were treated by Focused Ultrasound (FUs) performed with magnetic resonance (MR) guidance. Surgical treatment or radiotherapy treatment was not indicated for patients who underwent FUs. Lesions were located in the appendicular and axial skeleton and consisted of secondary symptomatic lesions. The clinical course of pain was evaluated using the Visual Analog Scale (VAS) before treatment, at 1 week, and at 1 month after treatment and the Oral Morphine Equivalent Daily Dose (OMEDD) was also recorded. We used Wilcoxon signed rank test to assess change in patient pain (R CRAN software V 3.1.1).

Results

We observed a significant decrease in the pain felt by patients between pre- procedure and 1 week post-procedure (p = 2.9.10-4), and pre-procedure and 1 month post-procedure (p = 3.10-4). The proportion of responders according to the International Bone Metastases Consensus Working Party was: Partial Response 50% (8/16) and Complete Response 37.5% (6/16).

Conclusions

HIFU under MR-guidance seems to be an effective and safe procedure in the treatment of symptomatic bone lesions for patients suffering from metastatic disease. A significant decrease of patient pain was observed.

Trial registration

NCT01091883. Registered 24 March 2010. Level of evidence: Level 3.

Quantitative MR thermometry based on phase-drift correction PRF shift method at 0.35 T

Background

Noninvasive magnetic resonance thermometry (MRT) at low-field using proton resonance frequency shift (PRFS) is a promising technique for monitoring ablation temperature, since low-field MR scanners with open-configuration are more suitable for interventional procedures than closed systems. In this study, phase-drift correction PRFS with first-order polynomial fitting method was proposed to investigate the feasibility and accuracy of quantitative MR thermography during hyperthermia procedures in a 0.35 T open MR scanner.

Methods

Unheated phantom and ex vivo porcine liver experiments were performed to evaluate the optimal polynomial order for phase-drift correction PRFS. The temperature estimation approach was tested in brain temperature experiments of three healthy volunteers at room temperature, and in ex vivo porcine liver microwave ablation experiments. The output power of the microwave generator was set at 40 W for 330 s. In the unheated experiments, the temperature root mean square error (RMSE) in the inner region of interest was calculated to assess the best-fitting order for polynomial fit. For ablation experiments, relative temperature difference profile measured by the phase-drift correction PRFS was compared with the temperature changes recorded by fiber optic temperature probe around the microwave ablation antenna within the target thermal region.

Results

The phase-drift correction PRFS using first-order polynomial fitting could achieve the smallest temperature RMSE in unheated phantom, ex vivo porcine liver and in vivo human brain experiments. In the ex vivo porcine liver microwave ablation procedure, the temperature error between MRT and fiber optic probe of all but six temperature points were less than 2 °C. Overall, the RMSE of all temperature points was 1.49 °C.

Conclusions

Both in vivo and ex vivo experiments showed that MR thermometry based on the phase-drift correction PRFS with first-order polynomial fitting could be applied to monitor temperature changes during microwave ablation in a low-field open-configuration whole-body MR scanner.

Passive acoustic mapping of cavitation using eigenspace-based robust Capon beamformer in ultrasound therapy

Pulse-echo imaging technique can only play a role when high intensity focused ultrasound (HIFU) is turned off due to the interference between the primary HIFU signal and the transmission pulse. Passive acoustic mapping (PAM) has been proposed as a tool for true real-time monitoring of HIFU therapy. However, the most-used PAM algorithm based on time exposure acoustic (TEA) limits the quality of cavitation image. Recently, robust Capon beamformer (RCB) has been used in PAM to provide improved resolution and reduced artifacts over TEA-based PAM, but the presented results have not been satisfactory. In the present study, we applied an eigenspace-based RCB (EISRCB) method to further improve the PAM image quality. The optimal weighting vector of the proposed method was found by projecting the RCB weighting vector onto the desired vector subspace constructed from the eigenstructure of the covariance matrix. The performance of the proposed PAM was validated by both simulations and in vitro histotripsy experiments. The results suggested that the proposed PAM significantly outperformed the conventionally used TEA and RCB-based PAM. The comparison results between pulse-echo images of the residual bubbles and cavitation images showed the potential of our proposed PAM in accurate localization of cavitation activity during HIFU therapy.

Acoustic field characterization of a clinical magnetic resonance-guided high-intensity focused ultrasound system inside the magnet bore

PURPOSE

With the expanding clinical application of magnetic resonance-guided high-intensity focused ultrasound (MR-HIFU), acoustic field characterization of MR-HIFU systems is needed for facilitating regulatory approval and ensuring consistent and safe power output of HIFU transducers. However, the established acoustic field measurement techniques typically use equipment that cannot be used in a magnetic resonance imaging (MRI) suite, thus posing a challenge to the development and execution of HIFU acoustic field characterization techniques. In this study, we developed and characterized a technique for HIFU acoustic field calibration within the MRI magnet bore, and validated the technique with standard hydrophone measurements outside of the MRI suite.

METHODS

A clinical Philips MR-HIFU system (Sonalleve V2, Philips, Vantaa, Finland) was used to assess the proposed technique. A fiber-optic hydrophone with a long fiber was inserted through a 24-gauge angiocatheter and fixed inside a water tank that was placed on the HIFU patient table above the acoustic window. The long fiber allowed the hydrophone control unit to be placed outside of the magnet room. The location of the fiber tip was traced on MR images, and the HIFU focal point was positioned at the fiber tip using the MR-HIFU therapy planning software. To perform acoustic field mapping inside the magnet, the HIFU focus was positioned relative to the fiber tip using an MRI-compatible 5-axis robotic transducer positioning system embedded in the HIFU patient table. To perform validation measurements of the acoustic fields, the HIFU table was moved out of the MRI suite, and a standard laboratory hydrophone measurement setup was used to perform acoustic field measurements outside the magnetic field.

RESULTS

The pressure field scans along and across the acoustic beam path obtained inside the MRI bore were in good agreement with those obtained outside of the MRI suite. At the HIFU focus with varying nominal acoustic powers of 10-500 W, the peak positive pressure and peak negative pressure measured inside the magnet bore were 3.87-68.67 MPa and 3.56-12.06 MPa, respectively, while outside the MRI suite the corresponding pressures were 3.27-67.32 MPa and 3.06-12.39 MPa, respectively. There was no statistically significant difference (P > 0.05) between measurements inside the magnet bore and outside the MRI suite for the p⁺ and p^- at any acoustic power level. The spatial-peak pulse-average intensities (I_SPPA ) for these powers were 312-17816 W/cm² and 220-15698 W/cm² for measurements inside and outside the magnet room, respectively. In addition, when the scanning step size of the HIFU focus was increased from 100 μm to 500 μm, the execution time for scanning a 4 × 4 mm² area decreased from 210 min to 10 min, the peak positive pressure decreased by 14%, the peak negative pressure decreased by 5%, and the lateral full width at half maximum dimension of pressure profiles increased from 1.15 mm to 1.55 mm.

CONCLUSIONS

The proposed hydrophone measurement technique offers a convenient and reliable method for characterizing the acoustic fields of clinical MR-HIFU systems inside the magnet bore. The technique was validated for use by measurements outside the MRI suite using a standard hydrophone calibration technique. This technique can be a useful tool in MR-HIFU quality assurance and acoustic field assessment.

Monitoring Microwave Ablation of Ex Vivo Bovine Liver Using Ultrasonic Attenuation Imaging

Thermal ablation of soft tissue changes the tissue microstructure and, consequently, induces changes in its acoustic properties. Although B-mode ultrasound provides high-resolution and high-frame-rate images of ablative therapeutic procedures, it is not particularly effective at delineating boundaries of ablated regions because of poor contrast in echogenicity between ablated and surrounding normal tissue. Quantitative ultrasound techniques can provide quantitative estimates of acoustic properties, such as backscatter and attenuation coefficients, and differentiate ablated and unablated regions more effectively, with the potential for monitoring minimally invasive thermal therapies. In this study, a previously introduced attenuation estimation method was used to create quantitative attenuation coefficient maps for 11 microwave ablation procedures performed on refrigerated ex vivo bovine liver. The attenuation images correlate well with the pathologic images of the ablated region. The mean attenuation coefficient for regions of interest drawn inside and outside the ablated zones were 0.9 (±0.2) and 0.45 (±0.15) dB/cm/MHz, respectively. These estimates agree with reported values in the literature and establish the usefulness of non-invasive attenuation imaging for monitoring therapeutic procedures in the liver.

Theoretical investigation of transgastric and intraductal approaches for ultrasound-based thermal therapy of the pancreas

Background

The goal of this study was to theoretically investigate the feasibility of intraductal and transgastric approaches to ultrasound-based thermal therapy of pancreatic tumors, and to evaluate possible treatment strategies.

Methods

This study considered ultrasound applicators with 1.2 mm outer diameter tubular transducers, which are inserted into the tissue to be treated by an endoscopic approach, either via insertion through the gastric wall (transgastric) or within the pancreatic duct lumen (intraductal). 8 patient-specific, 3D, transient, biothermal and acoustic finite element models were generated to model hyperthermia (n = 2) and ablation (n = 6), using sectored (210°–270°, n = 4) and 360° (n = 4) transducers for treatment of 3.3–17.0 cm³ tumors in the head (n = 5), body (n = 2), and tail (n = 1) of the pancreas. A parametric study was performed to determine appropriate treatment parameters as a function of tissue attenuation, blood perfusion rates, and distance to sensitive anatomy.

Results

Parametric studies indicated that pancreatic tumors up to 2.5 or 2.7 cm diameter can be ablated within 10 min with the transgastric and intraductal approaches, respectively. Patient-specific simulations demonstrated that 67.1–83.3% of the volumes of four sample 3.3–11.4 cm³ tumors could be ablated within 3–10 min using transgastric or intraductal approaches. 55.3–60.0% of the volume of a large 17.0 cm³ tumor could be ablated using multiple applicator positions within 20–30 min with either transgastric or intraductal approaches. 89.9–94.7% of the volume of two 4.4–11.4 cm³ tumors could be treated with intraductal hyperthermia. Sectored applicators are effective in directing acoustic output away from and preserving sensitive structures. When acoustic energy is directed towards sensitive structures, applicators should be placed at least 13.9–14.8 mm from major vessels like the aorta, 9.4–12.0 mm from other vessels, depending on the vessel size and flow rate, and 14 mm from the duodenum.

Conclusions

This study demonstrated the feasibility of generating shaped or conformal ablative or hyperthermic temperature distributions within pancreatic tumors using transgastric or intraductal ultrasound.

Endoluminal ultrasound applicators for MR-guided thermal ablation of pancreatic tumors: Preliminary design and evaluation in a porcine pancreas model

PURPOSE

Endoluminal ultrasound may serve as a minimally invasive option for delivering thermal ablation to pancreatic tumors adjacent to the stomach or duodenum. The objective of this study was to explore the basic feasibility of this treatment strategy through the design, characterization, and evaluation of proof-of-concept endoluminal ultrasound applicators capable of placement in the gastrointestinal (GI) lumen for volumetric pancreas ablation under MR guidance.

METHODS

Two variants of the endoluminal applicator, each containing a distinct array of two independently powered transducers (10 × 10 mm 3.2 MHz planar; or 8 × 10 × 20 mm radius of curvature 3.3 MHz curvilinear geometries) at the distal end of a meter long flexible catheter assembly, were designed and fabricated. Transducers and circulatory water flow for acoustic coupling and luminal cooling were contained by a low-profile polyester balloon covering the transducer assembly fixture. Each applicator incorporated miniature spiral MR coils and mechanical features (guiding tips and hinges) to facilitate tracking and insertion through the GI tract under MRI guidance. Acoustic characterization of each device was performed using radiation force balance and hydrophone measurements. Device delivery into the upper GI tract, adjacent to the pancreas, and heating characteristics for treatment of pancreatic tissue were evaluated in MR-guided ex vivo and in vivo porcine experiments. MR guidance was utilized for anatomical target identification, tracking/positioning of the applicator, and MR temperature imaging (MRTI) for PRF-based multislice thermometry, implemented in the real-time RTHawk software environment.

RESULTS

Force balance and hydrophone measurements indicated efficiencies of 48.8% and 47.8% and -3 dB intensity beam-widths of 3.2 and 1.2 mm for the planar and curvilinear transducers, respectively. Ex vivo studies on whole-porcine carcasses revealed capabilities of producing ablative temperature rise (ΔT > 15 °C) contours in pancreatic tissue 4-40 mm long and 4-28 mm wide for the planar transducer applicator (1-13 min sonication duration, ∼4 W/cm(2) applied acoustic intensity). Curvilinear transducers produced more selective heating, with a narrower ΔT > 15 °C contour length and width of up to 1-24 mm and 2-7 mm, respectively (1-7 min sonication duration, ∼4 W/cm(2) applied acoustic intensity). Active tracking of the miniature spiral coils was achieved using a Hadamard encoding tracking sequence, enabling real-time determination of each coil's coordinates and automated prescription of imaging planes for thermometry. In vivo MRTI-guided heating trials in three pigs demonstrated capability of ∼20 °C temperature elevation in pancreatic tissue at 2 cm depths from the applicator, with 5-7 W/cm(2) applied intensity and 6-16 min sonication duration. Dimensions of thermal lesions in the pancreas ranged from 12 to 28 mm, 3 to 10 mm, and 5 to 10 mm in length, width, and depth, respectively, as verified through histological analysis of tissue sections. Multiple-baseline reconstruction and respiratory-gated acquisition were demonstrated to be effective strategies in suppressing motion artifacts for clear evolution of temperature profiles during MRTI in the in vivo studies.

CONCLUSIONS

This study demonstrates the technical feasibility of generating volumetric ablation in pancreatic tissue using endoluminal ultrasound applicators positioned in the stomach lumen. MR guidance facilitates target identification, device tracking/positioning, and treatment monitoring through real-time multislice PRF-based thermometry.

Evaluation of a tissue-mimicking thermochromic phantom for radiofrequency ablation

PURPOSE

This work describes the characterization and evaluation of a tissue-mimicking thermochromic phantom (TMTCP) for direct visualization and quantitative determination of temperatures during radiofrequency ablation (RFA).

METHODS

TMTCP material was prepared using polyacrylamide gel and thermochromic ink that permanently changes color from white to magenta when heated. Color vs temperature calibration was generated in matlab by extracting RGB color values from digital photographs of phantom standards heated in a water bath at 25-75 °C. RGB and temperature values were plotted prior to curve fitting in mathematica using logistic functions of form f(t) = a + b/(1 + e((c(t-d)))), where a, b, c, and d are coefficients and t denotes temperature. To quantify temperatures based on TMTCP color, phantom samples were heated to temperatures blinded to the investigators, and two methods were evaluated: (1) visual comparison of sample color to the calibration series and (2) in silico analysis using the inverse of the logistic functions to convert sample photograph RGB values to absolute temperatures. For evaluation of TMTCP performance with RFA, temperatures in phantom samples and in a bovine liver were measured radially from an RF electrode during heating using fiber-optic temperature probes. Heating and cooling rates as well as the area under the temperature vs time curves were compared. Finally, temperature isotherms were generated computationally based on color change in bisected phantoms following RFA and compared to temperature probe measurements.

RESULTS

TMTCP heating resulted in incremental, permanent color changes between 40 and 64 °C. Visual and computational temperature estimation methods were accurate to within 1.4 and 1.9 °C between 48 and 67 °C, respectively. Temperature estimates were most accurate between 52 and 62 °C, resulting in differences from actual temperatures of 0.6 and 1.6 °C for visual and computational methods, respectively. Temperature measurements during RFA using fiber-optic probes matched closely with maximum temperatures predicted by color changes in the TMTCP. Heating rate and cooling rate, as well as the area under the temperature vs time curve were similar for TMTCP and ex vivo liver.

CONCLUSIONS

The TMTCP formulated for use with RFA can be used to provide quantitative temperature information in mild hyperthermic (40-45 °C), subablative (45-50 °C), and ablative (>50 °C) temperature ranges. Accurate visual or computational estimates of absolute temperatures and ablation zone geometry can be made with high spatial resolution based on TMTCP color. As such, the TMTCP can be used to assess RFA heating characteristics in a controlled, predictable environment.

Feasibility of real-time MR thermal dose mapping for predicting radiofrequency ablation outcome in the myocardium in vivo

BACKGROUND

Clinical treatment of cardiac arrhythmia by radiofrequency ablation (RFA) currently lacks quantitative and precise visualization of lesion formation in the myocardium during the procedure. This study aims at evaluating thermal dose (TD) imaging obtained from real-time magnetic resonance (MR) thermometry on the heart as a relevant indicator of the thermal lesion extent.

METHODS

MR temperature mapping based on the Proton Resonance Frequency Shift (PRFS) method was performed at 1.5 T on the heart, with 4 to 5 slices acquired per heartbeat. Respiratory motion was compensated using navigator-based slice tracking. Residual in-plane motion and related magnetic susceptibility artifacts were corrected online. The standard deviation of temperature was measured on healthy volunteers (N = 5) in both ventricles. On animals, the MR-compatible catheter was positioned and visualized in the left ventricle (LV) using a bSSFP pulse sequence with active catheter tracking. Twelve MR-guided RFA were performed on three sheep in vivo at various locations in left ventricle (LV). The dimensions of the thermal lesions measured on thermal dose images, on 3D T1-weighted (T1-w) images acquired immediately after the ablation and at gross pathology were correlated.

RESULTS

MR thermometry uncertainty was 1.5 °C on average over more than 96% of the pixels covering the left and right ventricles, on each volunteer. On animals, catheter repositioning in the LV with active slice tracking was successfully performed and each ablation could be monitored in real-time by MR thermometry and thermal dosimetry. Thermal lesion dimensions on TD maps were found to be highly correlated with those observed on post-ablation T1-w images (R = 0.87) that also correlated (R = 0.89) with measurements at gross pathology.

CONCLUSIONS

Quantitative TD mapping from real-time rapid CMR thermometry during catheter-based RFA is feasible. It provides a direct assessment of the lesion extent in the myocardium with precision in the range of one millimeter. Real-time MR thermometry and thermal dosimetry may improve safety and efficacy of the RFA procedure by offering a reliable indicator of therapy outcome during the procedure.

Combination of principal component analysis and optical-flow motion compensation for improved cardiac MR thermometry

The use of magnetic resonance (MR) thermometry for the monitoring of thermal ablation is rapidly expanding. However, this technique remains challenging for the monitoring of the treatment of cardiac arrhythmia by radiofrequency ablation due to the heart displacement with respiration and contraction. Recent studies have addressed this problem by compensating in-plane motion in real-time with optical-flow based tracking technique. However, these algorithms are sensitive to local variation of signal intensity on magnitude images associated with tissue heating. In this study, an optical-flow algorithm was combined with a principal component analysis method to reduce the impact of such effects. The proposed method was integrated to a fully automatic cardiac MR thermometry pipeline, compatible with a future clinical workflow. It was evaluated on nine healthy volunteers under free breathing conditions, on a phantom and in vivo on the left ventricle of a sheep. The results showed that local intensity changes in magnitude images had lower impact on motion estimation with the proposed method. Using this strategy, the temperature mapping accuracy was significantly improved.

A clinical study of thermal monitoring techniques of ultrasound-guided microwave ablation for hepatocellular carcinoma in high-risk locations

To confirm the safety and effectiveness of the minimally invasive thermal monitor technique on percutaneous ultrasound-guided microwave ablation (MWA) for hepatocellular carcinoma (HCC) in high-risk locations, a total of 189 patients with 226 HCC nodules in high-risk locations were treated with MWA. The real-time temperature of the tissue between the lesion margin and the vital structures was monitored by inserting a 21G thermal monitoring needle. The major indexes of technical success, technique effectiveness, local tumour progression and complications were observed during the follow-up period. Technical success was acquired in all patients. Technique effectiveness was achieved with one session in 119 lesions based on contrast-enhanced ultrasound (CEUS) 3-5 days after treatment. An additional 95 lesions achieved technique effectiveness at the second session. Within the follow-up period of 6-58 months (median 38 months), the 1-, 2-, 3-, and 4-year local tumour progression rate was 11.1%, 18.1%, 19.1%, and 19.9%, respectively. There were no major complications in all the patients except for the common side effects. These results indicate that the thermal monitor technique can be applied to prevent major complications in vulnerable structures and allow percutaneous MWA to achieve satisfactory technique effectiveness in the treatment of HCC in high-risk locations.

Laser thermal therapy monitoring using complex differential variance in optical coherence tomography

Conventional thermal therapy monitoring techniques based on temperature are often invasive, limited by point sampling, and are indirect measures of tissue injury, while techniques such as magnetic resonance and ultrasound thermometry are limited by their spatial resolution.  The visualization of the thermal coagulation zone at high spatial resolution is particularly critical to the precise delivery of thermal energy to epithelial lesions. In this work, an integrated thulium laser thermal therapy monitoring system was developed based on complex differential variance (CDV), which enables the 2D visualization of the dynamics of the thermal coagulation process at high spatial and temporal resolution with an optical frequency domain imaging system. With proper calibration to correct for noise, the CDV-based technique was shown to accurately delineate the thermal coagulation zone, which is marked by the transition from high CDV upon heating to a significantly reduced CDV once the tissue is coagulated, in 3 different tissue types ex vivo: skin, retina, and esophagus. The ability to delineate thermal lesions in multiple tissue types at high resolution opens up the possibility of performing microscopic image-guided procedures in a vast array of epithelial applications ranging from dermatology, ophthalmology, to gastroenterology and beyond.

Analysis of iodinated contrast delivered during thermal ablation: is material trapped in the ablation zone?

Intra-procedural contrast-enhanced CT (CECT) has been proposed to evaluate treatment efficacy of thermal ablation. We hypothesized that contrast material delivered concurrently with thermal ablation may become trapped in the ablation zone, and set out to determine whether such an effect would impact ablation visualization. CECT images were acquired during microwave ablation in normal porcine liver with: (A) normal blood perfusion and no iodinated contrast, (B) normal perfusion and iodinated contrast infusion or (C) no blood perfusion and residual iodinated contrast. Changes in CT attenuation were analyzed from before, during and after ablation to evaluate whether contrast was trapped inside of the ablation zone. Visualization was compared between groups using post-ablation contrast-to-noise ratio (CNR). Attenuation gradients were calculated at the ablation boundary and background to quantitate ablation conspicuity. In Group A, attenuation decreased during ablation due to thermal expansion of tissue water and water vaporization. The ablation zone was difficult to visualize (CNR  =  1.57  ±  0.73, boundary gradient  =  0.7  ±  0.4 HU mm(-1)), leading to ablation diameter underestimation compared to gross pathology. Group B ablations saw attenuation increase, suggesting that iodine was trapped inside the ablation zone. However, because the normally perfused liver increased even more, Group B ablations were more visible than Group A (CNR  =  2.04  ±  0.84, boundary gradient  =  6.3  ±  1.1 HU mm(-1)) and allowed accurate estimation of the ablation zone dimensions compared to gross pathology. Substantial water vaporization led to substantial attenuation changes in Group C, though the ablation zone boundary was not highly visible (boundary gradient  =  3.9  ±  1.1 HU mm(-1)). Our results demonstrate that despite iodinated contrast being trapped in the ablation zone, ablation visibility was highest when contrast is delivered intra-procedurally. Therefore, CECT may be feasible for real-time thermal ablation monitoring.

Sum-of-harmonics method for improved narrowband and broadband signal quantification during passive monitoring of ultrasound therapies

Passive Acoustic Mapping (PAM) enables real-time monitoring of ultrasound therapies by beamforming acoustic emissions emanating from the ultrasound focus. Reconstruction of the narrowband or broadband acoustic emissions component enables mapping of different physical phenomena, with narrowband emissions arising from non-linear propagation and scattering, non-inertial cavitation or tissue boiling, and broadband (generally, of significantly lower amplitude) indicating inertial cavitation. Currently, accurate classification of the received signals based on pre-defined frequency-domain comb filters cannot be guaranteed because varying levels of leakage occur as a function of signal amplitude and the choice of windowing function. This work presents a time-domain parametric model aimed at enabling accurate estimation of the amplitude of time-varying narrowband components in the presence of broadband signals. Conversely, the method makes it possible to recover a weak broadband signal in the presence of a dominant harmonic or other narrowband component. Compared to conventional comb filtering, the proposed sum-of-harmonics method enables PAM of cavitation sources that better reflect their physical location and extent.

Low-Pressure Burst-Mode Focused Ultrasound Wave Reconstruction and Mapping for Blood-Brain Barrier Opening: A Preclinical Examination

Burst-mode focused ultrasound (FUS) exposure has been shown to induce transient blood-brain barrier (BBB) opening for potential CNS drug delivery. FUS-BBB opening requires imaging guidance during the intervention, yet current imaging technology only enables postoperative outcome confirmation. In this study, we propose an approach to visualize short-burst low-pressure focal beam distribution that allows to be applied in FUS-BBB opening intervention on small animals. A backscattered acoustic-wave reconstruction method based on synchronization among focused ultrasound emission, diagnostic ultrasound receiving and passively beamformed processing were developed. We observed that focal beam could be successfully visualized for in vitro FUS exposure with 0.5–2 MHz without involvement of microbubbles. The detectable level of FUS exposure was 0.467 MPa in pressure and 0.05 ms in burst length. The signal intensity (SI) of the reconstructions was linearly correlated with the FUS exposure level both in-vitro (r² = 0.9878) and in-vivo (r² = 0.9943), and SI level of the reconstructed focal beam also correlated with the success and level of BBB-opening. The proposed approach provides a feasible way to perform real-time and closed-loop control of FUS-based brain drug delivery.

Fast temperature estimation from undersampled k-space with fully-sampled center for MR guided microwave ablation

PURPOSE

This study aims to accelerate MR temperature imaging using the proton resonance frequency (PRF) shift method for real time temperature monitoring during thermal ablation.

MATERIALS AND METHODS

The proposed method estimates temperature changes from undersampled k-space with a fully sampled center. This proposed algorithm is based on the hybrid multi-baseline and referenceless treatment image model and can be seen as an extension of the conventional k-space-based hybrid thermometry. The parameters of hybrid model are acquired by utilizing information from low resolution images which are obtained from fully-sampled centers of k-space. Registration is used to correct temperature errors due to the displacement of the subject. Phantom heating simulations, motion simulations, phantom heating and in-vivo experiments were performed to investigate the efficiency of the proposed method. SPIRiT and the conventional k-space estimation reconstruction thermometry were implemented for comparison using the same sampling pattern.

RESULTS

The phantom heating simulations showed that the proposed method results in lower RMSEs than the conventional k-space hybrid thermometry and SPIRiT at various reduction factors tested. The motion simulations indicated the robustness of the proposed method to displacement of the subject. Phantom heating experiment further demonstrated the ability of the method to reconstruct temperature maps with less computation time and higher accuracy (RMSEs lower than 0.4°C) at a net reduction factor of 3.5 in the presence of large noise caused by a microwave needle. In-vivo experiments validated the feasibility of the proposed method to estimate temperature changes from undersampled k-space (net reduction factor 4.3) in presence of respiratory motion and complicated anatomical structure, while reducing computation time as much as 10-fold compared with the conventional k-space method.

CONCLUSION

The proposed method accelerates the PRF-shift MR thermometry and provides more accurate temperature maps in presence of motion with relatively short computation time, which may make real time imaging for MR-guided microwave ablation possible.

Cluster analysis of DCE-MRI data identifies regional tracer-kinetic changes after tumor treatment with high intensity focused ultrasound

Evaluation of high intensity focused ultrasound (HIFU) treatment with MRI is generally based on assessment of the non-perfused volume from contrast-enhanced T1-weighted images. However, the vascular status of tissue surrounding the non-perfused volume has not been extensively investigated with MRI. In this study, cluster analysis of the transfer constant K(trans) and extravascular extracellular volume fraction ve , derived from dynamic contrast-enhanced MRI (DCE-MRI) data, was performed in tumor tissue surrounding the non-perfused volume to identify tumor subregions with distinct contrast agent uptake kinetics. DCE-MRI was performed in CT26.WT colon carcinoma-bearing BALB/c mice before (n = 12), directly after (n = 12) and 3 days after (n = 6) partial tumor treatment with HIFU. In addition, a non-treated control group (n = 6) was included. The non-perfused volume was identified based on the level of contrast enhancement. Quantitative comparison between non-perfused tumor fractions and non-viable tumor fractions derived from NADH-diaphorase histology showed a stronger agreement between these fractions 3 days after treatment (R(2) to line of identity = 0.91) compared with directly after treatment (R(2) = 0.74). Next, k-means clustering with four clusters was applied to K(trans) and ve parameter values of all significantly enhanced pixels. The fraction of pixels within two clusters, characterized by a low K(trans) and either a low or high ve , significantly increased after HIFU. Changes in composition of these clusters were considered to be HIFU induced. Qualitative H&E histology showed that HIFU-induced alterations in these clusters may be associated with hemorrhage and structural tissue disruption. Combined microvasculature and hypoxia staining suggested that these tissue changes may affect blood vessel functionality and thereby tumor oxygenation. In conclusion, it was demonstrated that, in addition to assessment of the non-perfused tumor volume, the presented methodology gives further insight into HIFU-induced effects on tumor vascular status. This method may aid in assessment of the consequences of vascular alterations for the fate of the tissue.

Investigation of factors affecting hypothermic pelvic tissue cooling using bio-heat simulation based on MRI-segmented anatomic models

This study applied a simulation method to map the temperature distribution based on magnetic resonance imaging (MRI) of individual patients, and investigated the influence of different pelvic tissue types as well as the choice of thermal property parameters on the efficiency of endorectal cooling balloon (ECB). MR images of four subjects with different prostate sizes and pelvic tissue compositions, including fatty tissue and venous plexus, were analyzed. The MR images acquired using endorectal coil provided a realistic geometry of deformed prostate that resembled the anatomy in the presence of ECB. A single slice with the largest two-dimensional (2D) cross-sectional area of the prostate gland was selected for analysis. The rectal wall, prostate gland, peri-rectal fatty tissue, peri-prostatic fatty tissue, peri-prostatic venous plexus, and urinary bladder were manually segmented. Pennes' bioheat thermal model was used to simulate the temperature distribution dynamics, by using an in-house finite element mesh based solver written in MATLAB. The results showed that prostate size and periprostatic venous plexus were two major factors affecting ECB cooling efficiency. For cases with negligible amount of venous plexus and small prostate, the average temperature in the prostate and neurovascular bundles could be cooled down to 25 °C within 30 min. For cases with abundant venous plexus and large prostate, the temperature could not reach 25 °C at the end of 3 h cooling. Large prostate made the cooling difficult to propagate through. The impact of fatty tissue on cooling effect was small. The filling of bladder with warm urine during the ECB cooling procedure did not affect the temperature in the prostate or NVB. In addition to the 2D simulation, in one case a 3D pelvic model was constructed for volumetric simulation. It was found that the 2D slice with the largest cross-sectional area of prostate had the most abundant venous plexus, and was the most difficult slice to cool, thus it may provide a conservative prediction of the cooling effect. This feasibility study demonstrated that the simulation tool could potentially be used for adjusting the setting of ECB for individual patients during hypothermic radical prostatectomy. Further studies using MR thermometry are required to validate the in silico results obtained using simulation.

Correcting heat-induced chemical shift distortions in proton resonance frequency-shift thermometry

PURPOSE

To reconstruct proton resonance frequency-shift temperature maps free of chemical shift distortions.

THEORY AND METHODS

Tissue heating created by thermal therapies such as focused ultrasound surgery results in a change in proton resonance frequency that causes geometric distortions in the image and calculated temperature maps, in the same manner as other chemical shift and off-resonance distortions if left uncorrected. We propose an online-compatible algorithm to correct these distortions in 2DFT and echo-planar imaging acquisitions, which is based on a k-space signal model that accounts for proton resonance frequency change-induced phase shifts both up to and during the readout. The method was evaluated with simulations, gel phantoms, and in vivo temperature maps from brain, soft tissue tumor, and uterine fibroid focused ultrasound surgery treatments.

RESULTS

Without chemical shift correction, peak temperature and thermal dose measurements were spatially offset by approximately 1 mm in vivo. Spatial shifts increased as readout bandwidth decreased, as shown by up to 4-fold greater temperature hot spot asymmetry in uncorrected temperature maps. In most cases, the computation times to correct maps at peak heat were less than 10 ms, without parallelization.

CONCLUSION

Heat-induced proton resonance frequency changes create chemical shift distortions in temperature maps resulting from MR-guided focused ultrasound surgery ablations, but the distortions can be corrected using an online-compatible algorithm. Magn Reson Med 76:172-182, 2016. © 2015 Wiley Periodicals, Inc.

Magnetic Resonance-Monitored Coaxial Electrochemical Ablation--Preliminary Evaluation of Technical Feasibility

PURPOSE

To evaluate the technical feasibility of a coaxial electrode configuration to rapidly create a mechanically defined electrochemical ablation zone monitored by magnetic resonance (MR) imaging in real time.

MATERIALS AND METHODS

A direct current generator supplied the nitinol cathode cage and central platinum anode for coaxial electrochemical ablation. Safety and efficacy were evaluated by measuring local pH, temperature, and current scatter in saline solutions. Ablation zone diameters of 3-6 cm (n = 72) were created on ex vivo bovine liver and verified by gross pathology. Feasibility of MR monitoring was evaluated using 8 swine livers to create ablations of 3 cm (n = 12), 4 cm (n = 4), and 5 cm (n = 4) verified by histology.

RESULTS

Local pH was 3.2 at the anode and 13.8 at the cathode. Current scatter was negligible. Ablation progress increased relative to local ion concentration, and MR signal changes corresponded to histologic findings. In the ex vivo model, the times to achieve complete ablation were 15 minutes, 20 minutes, 35 minutes, and 40 minutes for diameters of 3 cm, 4 cm, 5 cm, and 6 cm, respectively. Ablation times for the in situ model were 15 minutes, 35 minutes, and 50 minutes for 3 cm, 4 cm, and 5 cm, respectively.

CONCLUSIONS

The coaxial configuration mechanically defined the electrochemical ablation zone with times similar to comparably sized thermal ablations. MR compatibility allowed for real-time monitoring of ablation progress.

Spatiotemporal filtering of MR-temperature artifacts arising from bowel motion during transurethral MR-HIFU

PURPOSE

Transurethral MR-HIFU is a minimally invasive image-guided treatment for localized prostate cancer that enables precise targeting of tissue within the gland. The treatment is performed within a clinical MRI to obtain real-time MR thermometry used as an active feedback to control the spatial heating pattern in the prostate and to monitor for potential damage to surrounding tissues. This requires that the MR thermometry measurements are an accurate representation of the true tissue temperature. The proton resonance frequency shift thermometry method used is sensitive to tissue motion and changes in the local magnetic susceptibility that can be caused by the motion of air bubbles in the rectum, which can impact the performance of transurethral MR-HIFU in these regions of the gland.

METHODS

A method is proposed for filtering of temperature artifacts based on the temporal variance of the temperature, using empirical and dynamic positional knowledge of the ultrasonic heating beam, and an estimation of the measurement noise. A two-step correction strategy is introduced which eliminates artifact-detected temperature variations while keeping the noise level low through spatial averaging.

RESULTS

The filter has been evaluated by postprocessing data from five human transurethral ultrasound treatments. The two-step correction process led to reduced final temperature standard deviation in the prostate and rectum areas where the artifact was located, without negatively affecting areas distal to the artifact. The performance of the filter was also found to be consistent across all six of the data sets evaluated. The evaluation of the detection criterion parameter M determined that a value of M = 3 achieves a conservative filter with minimal loss of spatial resolution during the process.

CONCLUSIONS

The filter was able to remove most artifacts due to the presence of moving air bubbles in the rectum during transurethral MR-HIFU. A quantitative estimation of the filter capabilities shows a systematic improvement in the standard deviation of the corrected temperature maps in the rectum zone as well as in the entire acquired slice.