Regional deep hyperthermia: quantitative evaluation of predicted and direct measured temperature distributions in patients with high-risk extremity soft-tissue sarcoma

Robust, planning-based targeted locoregional tumour heating in small animals

Objective.To improve hyperthermia in clinical practice, pre-clinical hyperthermia research is essential to investigate hyperthermia effects and assess novel treatment strategies. Translating pre-clinical hyperthermia findings into clinically viable protocols requires laboratory animal treatment techniques similar to clinical hyperthermia techniques. The ALBA micro8 electromagnetic heating system (Med-logix SRL, Rome, Italy) has recently been developed to provide the targeted locoregional tumour heating currently lacking for pre-clinical research. This study evaluates the heat focusing properties of this device and its ability to induce robust locoregional tumour heating under realistic physiological conditions using simulations.Approach.Simulations were performed using the Plan2Heat treatment planning package (Amsterdam UMC, the Netherlands). First, the specific absorption rate (SAR) focus was characterised using a homogeneous phantom. Hereafter, a digital mouse model was used for the characterisation of heating robustness in a mouse. Device settings were optimised for treatment of a pancreas tumour and tested for varying circumstances. The impact of uncertainties in tissue property and perfusion values was evaluated using polynomial chaos expansion. Treatment quality and robustness were evaluated based on SAR and temperature distributions.Main results.The SAR distributions within the phantom are well-focused and can be adjusted to target any specific location. The focus size (full-width half-maximum) is a spheroid with diameters 9 mm (radially) and 20 mm (axially). The mouse model simulations show strong robustness against respiratory motion and intestine and stomach filling (∆T90≤0.14°C).Mouse positioning errors in the cranial-caudal direction lead to∆T90≤0.23°C. Uncertainties in tissue property and perfusion values were found to impact the treatment plan up to 0.56 °C (SD), with a variation onT90of 0.32 °C (1 SD).Significance.Our work shows that the pre-clinical phased-array system can provide adequate and robust locoregional heating of deep-seated target regions in mice. Using our software, robust treatment plans can be generated for pre-clinical hyperthermia research.

Quantification of tissue property and perfusion uncertainties in hyperthermia treatment planning: Multianalysis using polynomial chaos expansion

INTRODUCTION

Hyperthermia treatment planning (HTP) tools can guide treatment delivery, particularly with locoregional radiative phased array systems. Uncertainties in tissue and perfusion property values presently lead to quantitative inaccuracy of HTP, leading to sub-optimal treatment. Assessment of these uncertainties would allow for better judgement of the reliability of treatment plans and improve their value for treatment guidance. However, systematically investigating the impact of all uncertainties on treatment plans is a complex, high-dimensional problem and too computationally expensive for traditional Monte Carlo approaches. This study aims to systematically quantify the treatment-plan impact of tissue property uncertainties by investigating their individual contribution to, and combined impact on predicted temperature distributions.

METHODS

A novel Polynomial Chaos Expansion (PCE)-based HTP uncertainty quantification was developed and applied for locoregional hyperthermia of modelled tumours in the pancreatic head, prostate, rectum, and cervix. Patient models were based on the Duke and Ella digital human models. Using Plan2Heat, treatment plans were created to optimise tumour temperature (represented by T90) for treatment using the Alba4D system. For all 25-34 modelled tissues, the impact of tissue property uncertainties was analysed individually i.e., electrical and thermal conductivity, permittivity, density, specific heat capacity and perfusion. Next, combined analyses were performed on the top 30 uncertainties with the largest impact.

RESULTS

Uncertainties in thermal conductivity and heat capacity were found to have negligible impact on the predicted temperature ( < 1 × 10-10 °C), density and permittivity uncertainties had a small impact (< 0.3 °C). Uncertainties in electrical conductivity and perfusion can lead to large variations in predicted temperature. However, variations in muscle properties result in the largest impact at locations that could limit treatment quality, with a standard deviation up to almost 6 °C (pancreas) and 3.5 °C (prostate) for perfusion and electrical conductivity, respectively. The combined influence of all significant uncertainties leads to large variations with a standard deviation up to 9.0, 3.6, 3.7 and 4.1 °C for the pancreatic, prostate, rectal and cervical cases, respectively.

CONCLUSION

Uncertainties in tissue and perfusion property values can have a large impact on predicted temperatures from hyperthermia treatment planning. PCE-based analysis helps to identify all major uncertainties, their impact and judge the reliability of treatment plans.

Field and Temperature Shaping for Microwave Hyperthermia: Recent Treatment Planning Tools to Enhance SAR-Based Procedures

The aim of the article is to provide a summary of the work carried out in the framework of a research project funded by the Italian Ministry of Research. The main goal of the activity was to introduce multiple tools for reliable, affordable, and high-performance microwave hyperthermia for cancer therapy. The proposed methodologies and approaches target microwave diagnostics, accurate in vivo electromagnetic parameters estimation, and improvement in treatment planning using a single device. This article provides an overview of the proposed and tested techniques and shows their complementarity and interconnection. To highlight the approach, we also present a novel combination of specific absorption rate optimization via convex programming with a temperature-based refinement method implemented to mitigate the effect of thermal boundary conditions on the final temperature map. To this purpose, numerical tests were carried out for both simple and anatomically detailed 3D scenarios for the head and neck region. These preliminary results show the potential of the combined technique and improvements in the temperature coverage of the tumor target with respect to the case wherein no refinement is adopted.

Antenna Arrangement in UWB Helmet Brain Applicators for Deep Microwave Hyperthermia

Deep microwave hyperthermia applicators are typically designed as narrow-band conformal antenna arrays with equally spaced elements, arranged in one or more rings. This solution, while adequate for most body regions, might be sub-optimal for brain treatments. The introduction of ultra-wide-band semi-spherical applicators, with elements arranged around the head and not necessarily aligned, has the potential to enhance the selective thermal dose delivery in this challenging anatomical region. However, the additional degrees of freedom in this design make the problem non-trivial. We address this by treating the antenna arrangement as a global SAR-based optimization process aiming at maximizing target coverage and hot-spot suppression in a given patient. To enable the quick evaluation of a certain arrangement, we propose a novel E-field interpolation technique which calculates the field generated by an antenna at any location around the scalp from a limited number of initial simulations. We evaluate the approximation error against full array simulations. We demonstrate the design technique in the optimization of a helmet applicator for the treatment of a medulloblastoma in a paediatric patient. The optimized applicator achieves 0.3 °C higher T90 than a conventional ring applicator with the same number of elements.

Using MRI to measure position and anatomy changes and assess their impact on the accuracy of hyperthermia treatment planning for cervical cancer

PURPOSE

We studied the differences between planning and treatment position, their impact on the accuracy of hyperthermia treatment planning (HTP) predictions, and the relevance of including true treatment anatomy and position in HTP based on magnetic resonance (MR) images.

MATERIALS AND METHODS

All volunteers were scanned with an MR-compatible hyperthermia device, including a filled waterbolus, to replicate the treatment setup. In the planning setup, the volunteers were scanned without the device to reproduce the imaging in the current HTP. First, we used rigid registration to investigate the patient position displacements between the planning and treatment setup. Second, we performed HTP for the planning anatomy at both positions and the treatment mimicking anatomy to study the effects of positioning and anatomy on the quality of the simulated hyperthermia treatment. Treatment quality was evaluated using SAR-based parameters.

RESULTS

We found an average displacement of 2 cm between planning and treatment positions. These displacements caused average absolute differences of ∼12% for TC25 and 10.4%-15.9% in THQ. Furthermore, we found that including the accurate treatment position and anatomy in treatment planning led to an improvement of 2% in TC25 and 4.6%-10.6% in THQ.

CONCLUSIONS

This study showed that precise patient position and anatomy are relevant since these affect the accuracy of HTP predictions. The major part of improved accuracy is related to implementing the correct position of the patient in the applicator. Hence, our study shows a clear incentive to accurately match the patient position in HTP with the actual treatment.

Application of photogrammetry reconstruction for hyperthermia quality control measurements

PURPOSE

Hyperthermia is a cancer treatment in which the target region is heated to temperatures of 40-44 °C usually applying external electromagnetic field sources. The behavior of the hyperthermia applicators (antennas) in clinical practice should be periodically checked with phantom experiments to verify the applicator's performance over time. The purpose of this study was to investigate the application of photogrammetry reconstructions of 3D applicator position in these quality control procedure measurements.

METHODS

Photogrammetry reconstruction was applied at superficial hyperthermia scenario using the Lucite cone applicator (LCA) and phased-array heating in the head and neck region using the HYPERcollar3D. Wire-frame models of the entire measurement setups were created from multiple-view images and used for recreation of the setup inside 3D electromagnetic field simulation software. We evaluated applicator relation (R_a) between measured and simulated absolute specific absorption rate (SAR) for manually created and photogrammetry reconstructed simulation setups.

RESULTS

We found a displacement of 7.9 mm for the LCA and 8.2 mm for the HYPERcollar3D setups when comparing manually created and photogrammetry reconstructed applicator models placements. R_a improved from 1.24 to 1.18 for the LCA and from 1.17 to 1.07 for the HYPERcollar3D when using photogrammetry reconstructed simulation setups.

CONCLUSION

Photogrammetry reconstruction technique holds promise to improve measurement setup reconstruction and agreement between measured and simulated absolute SAR.

Experimental Validation of the MRcollar: An MR Compatible Applicator for Deep Heating in the Head and Neck Region

Clinical effectiveness of hyperthermia treatments, in which tumor tissue is artificially heated to 40-44 °C for 60-90 min, can be hampered by a lack of accurate temperature monitoring. The need for noninvasive temperature monitoring in the head and neck region (H&N) and the potential of MR thermometry prompt us to design an MR compatible hyperthermia applicator: the MRcollar. In this work, we validate the design, numerical model, and MR performance of the MRcollar. The MRcollar antennas have low reflection coefficients (<-15 dB) and the intended low interaction between the individual antenna modules (<-32 dB). A 10 °C increase in 3 min was reached in a muscle-equivalent phantom, such that the specifications from the European Society for Hyperthermic Oncology were easily reached. The MRcollar had a minimal effect on MR image quality and a five-fold improvement in SNR was achieved using the integrated coils of the MRcollar, compared to the body coil. The feasibility of using the MRcollar in an MR environment was shown by a synchronous heating experiment. The match between the predicted SAR and measured SAR using MR thermometry satisfied the gamma criteria [distance-to-agreement = 5 mm, dose-difference = 7%]. All experiments combined show that the MRcollar delivers on the needs for MR-hyperthermia in the H&N and is ready for in vivo investigation.

ESHO benchmarks for computational modeling and optimization in hyperthermia therapy

BACKGROUND

The success of cancer hyperthermia (HT) treatments is strongly dependent on the temperatures achieved in the tumor and healthy tissues as it correlates with treatment efficacy and safety, respectively. Hyperthermia treatment planning (HTP) simulations have become pivotal for treatment optimization due to the possibility for pretreatment planning, optimization and decision making, as well as real-time treatment guidance.

MATERIALS AND METHODS

The same computational methods deployed in HTP are also used for in silico studies. These are of great relevance for the development of new HT devices and treatment approaches. To aid this work, 3 D patient models have been recently developed and made available for the HT community. Unfortunately, there is no consensus regarding tissue properties, simulation settings, and benchmark applicators, which significantly influence the clinical relevance of computational outcomes.

RESULTS AND DISCUSSION

Herein, we propose a comprehensive set of applicator benchmarks, efficacy and safety optimization algorithms, simulation settings and clinical parameters, to establish benchmarks for method comparison and code verification, to provide guidance, and in view of the 2021 ESHO Grand Challenge (Details on the ESHO grand challenge on HTP will be provided at https://www.esho.info/).

CONCLUSION

We aim to establish guidelines to promote standardization within the hyperthermia community such that novel approaches can quickly prove their benefit as quickly as possible in clinically relevant simulation scenarios. This paper is primarily focused on radiofrequency and microwave hyperthermia but, since 3 D simulation studies on heating with ultrasound are now a reality, guidance as well as a benchmark for ultrasound-based hyperthermia are also included.

Suitability of eigenvalue beam-forming for discrete multi-frequency hyperthermia treatment planning

PURPOSE

Thermal dose delivery in microwave hyperthermia for cancer treatment is expected to benefit from the introduction of ultra-wideband (UWB)-phased array applicators. A full exploitation of the combination of different frequencies to improve the deposition pattern is, however, a nontrivial problem. It is unclear whether the cost functions used for hyperthermia treatment planning (HTP) optimization in the single-frequency setting can be meaningfully extended to the UWB case.

METHOD

We discuss the ability of the eigenvalue (EV) and a novel implementation of iterative-EV (i-EV) beam-forming methods to fully exploit the available frequency spectrum when a discrete set of simultaneous operating frequencies is available for treatment. We show that the quadratic power deposition ratio solved by the methods can be maximized by only one frequency in the set, therefore rendering EV inadequate for UWB treatment planning. We further investigate whether this represents a limitation in two realistic test cases, comparing the thermal distributions resulting from EV and i-EV to those obtained by optimizing for other nonlinear cost functions that allow for multi-frequency.

RESULTS

The classical EV-based single-frequency HTP yields systematically lower target SAR deposition and temperature values than nonlinear HTP. In a larynx target, the proposed single-frequency i-EV scheme is able to compensate for this and reach temperatures comparable to those given by global nonlinear optimization. In a meninges target, the multi-frequency setting outperforms the single-frequency one, achieving better target coverage and 0 . 5 ∘ C higher T 90 in the tumor than single-frequency-based HTP.

CONCLUSIONS

Classical EV performs poorly in terms of resulting target temperatures. The proposed single-frequency i-EV scheme can be a viable option depending on the patient and tumor to be treated, as long as the proper operating frequency can be selected across a UWB range. Multi-frequency HTP can bring a considerable benefit in regions typically difficult to treat such as the brain.

Dielectric Properties of Human Active Liver, Kidney and Spleen Compared to Those of Respective Inactive Tissues, Porcine Tissues and the Data Provided by a Database in the Frequency Range of 10 Hz to 100 MHz

OBJECTIVE

The purpose of this work is to study whether the active state and species of biological tissues can influence changes in their dielectric properties.

METHODS

In this paper, the dielectric properties of liver, kidney and spleen tissues from human active, human inactive and animal tissues are measured in the frequency range of 10 Hz to 100 MHz. The four- and two-electrode methods are used to measure dielectric properties at different frequencies. Statistical analysis and the pattern recognition method are used to compare the dielectric properties of human active tissues, human inactive tissues, animal tissues and data provided by the IFAC database.

RESULTS

The results show that the dielectric properties of human active tissues are significantly different from those of human inactive tissues and animal tissues, resulting in a great difference between the dielectric properties provided by the IFAC database and those of human active tissues. The dielectric properties of human active tissues can be identified by the pattern recognition method based on principal component analysis, which further proves that the dielectric properties of human active tissues cannot be replaced.

CONCLUSION

The dielectric properties of biological tissues are closely related to the activity and species of tissues. The dielectric properties of human active tissues cannot be replaced by those of human cadaver tissues or animal tissues.

SIGNIFICANCE

The significance of this study is suggesting that the IFAC database should be updated with the dielectric properties of human active tissues to provide accurate data for bioelectromagnetics research.

Contactless Thermometry by MRI and MRS: Advanced Methods for Thermotherapy and Biomaterials

Control of temperature variation is of primordial importance in particular areas of biomedicine. In this context, medical treatments such as hyperthermia and cryotherapy, and also the development and use of hydrogel-based biomaterials, are of particular concern. To enable accurate temperature measurement without perturbing or even destroying the biological tissue or material to be monitored, contactless thermometry methods are preferred. Among these, the most suitable are based on magnetic resonance imaging and spectroscopy (MRI, MRS). Here, we address the latest developments in this field as well as their current and anticipated practical applications. We highlight recent progress aimed at rendering MR thermometry faster and more reproducible, versatile, and sophisticated and provide our perspective on how these new techniques broaden the range of applications in medical treatments and biomaterial development by enabling insight into finer details of thermal behavior. Thus, these methods facilitate optimization of clinical and industrial heating and cooling protocols.

Quantitative, Multi-institutional Evaluation of MR Thermometry Accuracy for Deep-Pelvic MR-Hyperthermia Systems Operating in Multi-vendor MR-systems Using a New Anthropomorphic Phantom

Clinical outcome of hyperthermia depends on the achieved target temperature, therefore target conformal heating is essential. Currently, invasive temperature probe measurements are the gold standard for temperature monitoring, however, they only provide limited sparse data. In contrast, magnetic resonance thermometry (MRT) provides unique capabilities to non-invasively measure the 3D-temperature. This study investigates MRT accuracy for MR-hyperthermia hybrid systems located at five European institutions while heating a centric or eccentric target in anthropomorphic phantoms with pelvic and spine structures. Scatter plots, root mean square error (RMSE) and Bland-Altman analysis were used to quantify accuracy of MRT compared to high resistance thermistor probe measurements. For all institutions, a linear relation between MRT and thermistor probes measurements was found with R² (mean ± standard deviation) of 0.97 ± 0.03 and 0.97 ± 0.02, respectively for centric and eccentric heating targets. The RMSE was found to be 0.52 ± 0.31 °C and 0.30 ± 0.20 °C, respectively. The Bland-Altman evaluation showed a mean difference of 0.46 ± 0.20 °C and 0.13 ± 0.08 °C, respectively. This first multi-institutional evaluation of MR-hyperthermia hybrid systems indicates comparable device performance and good agreement between MRT and thermistor probes measurements. This forms the basis to standardize treatments in multi-institution studies of MR-guided hyperthermia and to elucidate thermal dose-effect relations.