Thermal diffusivity and perfusion constants from in vivo MR-guided focussed ultrasound treatments: a feasibility study

Magnetic resonance imaging of blood perfusion rate based on Helmholtz decomposition of heat flux

Objective.Thermal property (TP) maps of human tissues are useful for tumor treatment and diagnosis. In particular, the blood perfusion rate is significantly different for tumors and healthy tissues. Noninvasive techniques that reconstruct TPs from the temperature measured via magnetic resonance imaging (MRI) by solving an inverse bioheat transfer problem have been developed. A few conventional methods can reconstruct spatially varying TP distributions, but they have several limitations. First, most methods require the numerical Laplacian computation of the temperature, and hence they are sensitive to noise. In addition, some methods require the division of a region of interest (ROI) into sub-regions with homogeneous TPs using prior anatomical information, and they assume an unmeasurable initial temperature distribution. We propose a novel robust reconstruction method without the division of an ROI or the assumption of an initial temperature distribution.Approach.The proposed method estimates blood perfusion rate maps from relative temperature changes. This method avoids the computation of the Laplacian by using integral representations of the Helmholtz decomposition of the heat flux.Main Result.We compare the reconstruction results of the conventional and proposed methods using numerical simulations. The results indicate the robustness of the proposed method.Significance.This study suggests the feasibility of thermal property mapping with MRI using the robust proposed method.

Comparison of computer simulations and clinical treatment results of magnetic resonance-guided focused ultrasound surgery (MRgFUS) of uterine fibroids

PURPOSE

Magnetic resonance-guided focused ultrasound surgery (MRgFUS) can be used to noninvasively treat symptomatic uterine fibroids by heating with focused ultrasound sonications while monitoring the temperature with magnetic resonance (MR) thermometry. While prior studies have compared focused ultrasound simulations to clinical results, studies involving uterine fibroids remain scarce. In our study, we perform such a comparison to assess the suitability of simulations for treatment planning.

METHODS

Sonications (N = 67) were simulated retrospectively using acoustic and thermal models based on the Rayleigh integral and Pennes bioheat equation followed by MR-thermometry simulation in seven patients who underwent MRgFUS treatment for uterine fibroids. The spatial accuracy of simulated focus location was assessed by evaluating displacements of the centers of mass of the thermal dose distributions between simulated and treatment MR thermometry slices. Temperature-time curves and sizes of 240 equivalent minutes at 43°C (240EM43 ) volumes between treatment and simulation were compared.

RESULTS

The simulated focus location showed errors of 2.7 ± 4.1, -0.7 ± 2.0, and 1.3 ± 1.2 mm (mean ± SD) in the anterior-posterior, foot-head, and right-left directions for a fibroid absorption coefficient of 4.9 Np m-1  MHz-1 and perfusion parameter of 1.89 kg m-3 s-1 . Linear regression of 240EM43 volumes of 67 sonications of patient treatments and simulations utilizing these parameters yielded a slope of 1.04 and a correlation coefficient of 0.54. The temperature rise ratio of simulation to treatment near the end of sonication was 0.47 ± 0.22, 1.28 ± 0.60, and 1.49 ± 0.71 for 66 sonications simulated utilizing fibroid absorption coefficient of 1.2, 4.9, and 8.6 Np m-1  MHz-1 , respectively, and the aforementioned perfusion value. The impact of perfusion on peak temperature rise is minimal between 1.89 and 10 kg m-3 s-1 , but became more substantial when utilizing a value of 100 kg m-3 s-1 .

CONCLUSIONS

The results of this study suggest that perfusion, while in some cases having a substantial impact on thermal dose volumes, has less impact than ultrasound absorption for predicting peak temperature elevation at least when using perfusion parameter values up to 10 kg m-3 s-1 for this particular array geometry, frequencies, and tissue target which is good for clinicians to be aware of. The results suggest that simulations show promise in treatment planning, particularly in terms of spatial accuracy. However, in order to use simulations to predict temperature rise due to a sonication, knowledge of the patient-specific tissue parameters, in particular the absorption coefficient is important. Currently, spatially varying patient-specific tissue parameter values are not available during treatment, so simulations can only be used for planning purposes to estimate sonication performance on average.

Visualization of thermal washout due to spatiotemporally heterogenous perfusion in the application of a model-based control algorithm for MR-HIFU mediated hyperthermia

PURPOSE

This article will report results from the in-vivo application of a previously published model-predictive control algorithm for MR-HIFU hyperthermia. The purpose of the investigation was to test the controller's in-vivo performance and behavior in the presence of heterogeneous perfusion.

MATERIALS AND METHODS

Hyperthermia at 42°C was induced and maintained for up to 30 min in a circular section of a thermometry slice in the biceps femoris of German landrace pigs (n=5) using a commercial MR-HIFU system and a recently developed MPC algorithm. The heating power allocation was correlated with heat sink maps and contrast-enhanced MRI images. The temporal change in perfusion was estimated based on the power required to maintain hyperthermia.

RESULTS

The controller performed well throughout the treatments with an absolute average tracking error of 0.27 ± 0.15 °C and an average difference of 1.25 ± 0.22 °C between T10 and T90. The MPC algorithm allocates additional heating power to sub-volumes with elevated heat sink effects, which are colocalized with blood vessels visible on contrast-enhanced MRI. The perfusion appeared to have increased by at least a factor of ∼1.86 on average.

CONCLUSIONS

The MPC controller generates temperature distributions with a narrow spectrum of voxel temperatures inside the target ROI despite the presence of spatiotemporally heterogeneous perfusion due to the rapid thermometry feedback available with MR-HIFU and the flexible allocation of heating power. The visualization of spatiotemporally heterogeneous perfusion presents new research opportunities for the investigation of stimulated perfusion in hypoxic tumor regions.

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 tissue preparation to characterize uterine fibroid tissue properties for thermal therapies

PURPOSE

Treating uterine fibroids with less invasive therapies such as magnetic resonance-guided focused ultrasound (MRgFUS) is an attractive alternative to surgery. Treatment planning can improve MRgFUS procedures and reduce treatment times, but the tissue properties that currently inform treatment planning tools are not adequate. This study aims to develop an ex vivo uterine fibroid model that can emulate the in vivo environment allowing for characterization of the uterus and fibroid MR, acoustic, and thermal tissue properties while maintaining viability for the necessary postsurgical histopathological assessments.

METHODS

Women undergoing a hysterectomy due to fibroid-related symptoms were invited to undergo a preoperative pelvic MRI and to permit postoperative testing of their uterine specimen. Patients that declined or could not be scheduled for a pre-operative MRI were still able to allow post-operative testing of their excised tissue. Following surgical removal of the uterus, nonmorcellated tissues were reperfused with a Krebs-Henseleit buffer solution. An MR-compatible perfusion system was designed to maintain tissue viability inside the MR suite during scanning. MR imaging protocols utilized preoperatively were repeated on whole sample, reperfused ex vivo uterus specimens. Thermal properties including thermal diffusivity and thermal conductivity of the uterus and fibroids were determined using an invasive needle sensor device in 50% of the specimens. Acoustic property measurements (density, speed of sound and attenuation) were obtained for approximately 20% of the tissue samples using both through-transmission and radiation force balance techniques. Differences between fibroid and uterus and in vivo and ex vivo measurements were evaluated with a two-tailed Student t test.

RESULTS

Fourteen patients participated in the study and measurements were obtained from 22 unique fibroids. Of the 16 fibroids available for preoperative MRI testing, 69% demonstrated classic hypo-intensity relative to the myometrium, with the remainder presenting with iso- (25%) or hyper-intensity (6%). While thermal diffusivity was not significantly different between fibroid and myometrium tissues (0.217 ± 0.047 and 0.204 ± 0.039 mm² /s, respectively), the acoustic attenuation in fibroid tissue was significantly higher than myometrium (0.092 ± 0.021 and 0.052 ± 0.023 Np/cm/MHz, respectively). When comparing in vivo with ex vivo MRI T1 and T2 measurements in fibroids and myometrium tissue, the only difference was found in the fibroid T2 property (P < 0.05). Finally, the developed perfusion protocol successfully maintained tissue viability in ex vivo tissues as evaluated through histological analysis.

CONCLUSIONS

This study developed an MR-compatible extracorporeal perfusion technique that effectively maintains tissue viability, allowing for the direct measurement of patient-specific MR, thermal, and acoustic property values for both fibroid and myometrium tissues. These measured tissue property values will enable further development and validation of treatment planning models that can be utilized during MRgFUS uterine fibroid treatments.

Validation of hybrid angular spectrum acoustic and thermal modelling in phantoms

In focused ultrasound (FUS) thermal ablation of diseased tissue, acoustic beam and thermal simulations enable treatment planning and optimization. In this study, a treatment-planning methodology that uses the hybrid angular spectrum (HAS) method and the Pennes' bioheat equation (PBHE) is experimentally validated in homogeneous tissue-mimicking phantoms. Simulated three-dimensional temperature profiles are compared to volumetric MR thermometry imaging (MRTI) of FUS sonications in the phantoms, whose acoustic and thermal properties are independently measured. Additionally, Monte Carlo (MC) uncertainty analysis is performed to quantify the effect of tissue property uncertainties on simulation results. The mean error between simulated and experimental spatiotemporal peak temperature rise was +0.33°C (+6.9%). Despite this error, the experimental temperature rise fell within the expected uncertainty of the simulation, as determined by the MC analysis. The average errors of the simulated transverse and longitudinal full width half maximum (FWHM) of the profiles were -1.9% and 7.5%, respectively. A linear regression and local sensitivity analysis revealed that simulated temperature amplitude is more sensitive to uncertainties in simulation inputs than in the profile width and shape. Acoustic power, acoustic attenuation and thermal conductivity had the greatest impact on peak temperature rise uncertainty; thermal conductivity and volumetric heat capacity had the greatest impact on FWHM uncertainty. This study validates that using the HAS and PBHE method can adequately predict temperature profiles from single sonications in homogeneous media. Further, it informs the need to accurately measure or predict patient-specific properties for improved treatment planning of ablative FUS surgeries.

Experimental assessment of phase aberration correction for breast MRgFUS therapy

PURPOSE

This study validates that phase aberrations in breast magnetic resonance-guided focussed ultrasound (MRgFUS) therapies can be corrected in a clinically relevant time frame to generate more intense, smaller and more spatially accurate foci.

MATERIALS AND METHODS

Hybrid angular spectrum (HAS) ultrasound calculations in an magnetic resonance imaging (MRI)-based tissue model, were used to compute phase aberration corrections for improved experimental MRgFUS heating in four heterogeneous breast-mimicking phantoms (n = 18 total locations). Magnetic resonance(MR) temperature imaging was used to evaluate the maximum temperature rise, focus volume and focus accuracy for uncorrected and phase aberration-corrected sonications. Thermal simulations assessed the effectiveness of the phase aberration correction implementation.

RESULTS

In 13 of 18 locations, the maximum temperature rise increased by an average of 30%, focus volume was reduced by 40% and focus accuracy improved from 4.6 to 3.6 mm. Mixed results were observed in five of the 18 locations, with focus accuracy improving from 6.1 to 2.5 mm and the maximum temperature rise decreasing by 8% and focus volume increasing by 10%. Overall, the study demonstrated significant improvements (p < 0.005) in maximum temperature rise, focus volume and focus accuracy. Simulations predicted greater improvements than observed experimentally, suggesting potential for improvement in implementing the technique. The complete phase aberration correction procedure, including model generation, segmentation and phase aberration computations, required less than 45 min per sonication location.

CONCLUSION

The significant improvements demonstrated in this study i.e., focus intensity, size and accuracy from phase aberration correction have the potential to improve the efficacy, time-efficiency and safety of breast MRgFUS therapies.