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

Computational predictions of magnetic resonance acoustic radiation force imaging for breast cancer focused ultrasound therapy

PURPOSE

In magnetic resonance-guided focused ultrasound (MRgFUS) breast therapies, the focal location must be characterized to guide successful treatment. Focal characterization is difficult because heterogeneous breast tissues introduce phase aberrations that blur and shift the focus and traditional guidance methods do not work in adipose tissues. The purpose of this work is to evaluate numerical simulations of MRgFUS that predict the focal location. Those simulations are compared to clinical magnetic resonance acoustic radiation force imaging (MR-ARFI) data collected during in vivo treatment of breast tumors.

METHODS

The focal location was evaluated before MRgFUS treatment with MR-ARFI in five patients. The hybrid angular spectrum method (HAS) was applied to simulate pressure fields which were converted to forces, then convolved with a 3D Green's function (with time-of-arrival weighting) to produce a simulation of the MR-ARFI tissue displacement.

RESULTS

The focal locations found by the simulations and the MR-ARFI measurements were on average separated by 3.7 mm (SD: 0.9 mm). Characterization of the focal zone spatial distributions had a normalized root mean squared difference of 8.1% (SD: 2.5%). The displacement magnitudes of the simulations underestimated the MR-ARFI measurements by 82% (SD: 5.6%).

CONCLUSIONS

The agreement between MR-ARFI measurements and simulations demonstrates that HAS can predict the in vivo focal location in heterogeneous tissues, though accurate patient-specific properties are needed to improve predictions of tissue displacement magnitude. Tools developed in this study could be used to streamline MRgFUS treatment planning and optimization, for biomechanical property estimation, and in developing phase aberration correction techniques.

Effectiveness of fat suppression methods and influence on proton-resonance frequency shift (PRFS) MR thermometry

PURPOSE

To evaluate the effectiveness of fat suppression techniques experimentally and illustrate their influence on the accuracy of PRFS MR-thermometry.

METHODS

The residual magnitudes of the main fat peaks are measured using a water-fat decomposition algorithm in an oil phantom and in vivo in swine bone marrow, either with spectral fat saturation (FS), water excitation (WE) or fast water excitation (FWE), as implemented on 1.5 T whole-body clinical MRIs. Thermometry experiments in tissue-mimicking oil-water phantoms (10 and 30 % fat) allow determining temperature errors in PRFS MR-thermometry with no fat suppression, FS and WE, compared against reference fiber optic thermometry.

RESULTS

WE attenuates the signal of the main methylene fat peak more than FS (2 % and 22 % amplitude attenuation in the oil phantom, respectively), while the olefinic and glycerol peaks surrounding the water peak remain unaltered with both FS and WE. Within the 37 °C to 60 °C temperature range explored, FS and WE strongly attenuate temperature errors compared to PRFS without fat suppression. The residual fat signal after FS and WE leads to errors in PRFS thermometry, that increase with the fat content and oscillate with TE and temperature. In our tests limited to a single MR provider, fat suppression with WE appears to suppress fat signal more effectively.

CONCLUSIONS

We propose a protocol to quantify the remaining fraction of each spectral fat peak after fat suppression. In PRFS thermometry, despite spectral fat suppression, the remnant fat signal leads to temperature underestimation or overestimation depending on TE, fat fraction and temperature range. Fat suppression techniques should be evaluated specifically for quantitative MRI methods such as PRFS thermometry.

In vivo simultaneous proton resonance frequency shift thermometry and single reference variable flip angle T1 measurements

PURPOSE

The single reference variable flip angle sequence with a multi-echo stack of stars acquisition (SR-VFA-SoS) simultaneously measures temperature change using proton resonance frequency (PRF) shift and T1-based thermometry methods. This work evaluates SR-VFA-SoS thermometry in MR-guided focused ultrasound in an in vivo rabbit model.

METHODS

Simultaneous PRF shift thermometry and T1-based thermometry were obtained in a New Zealand white rabbit model (n = 7) during MR-guided focused ultrasound surgery using the SR-VFA-SoS sequence at 3 T. Distinct locations in muscle (n = 16), fat (n = 12), or the interface of both tissues (n = 23) were heated. The T1-temperature coefficient of fat was determined using least-squares fitting of inversion recovery-based T1 maps of untreated fat harvested from the animal and was applied to the in vivo measured heat-induced T1 changes to create temperature maps.

RESULTS

Using k-space weighted image contrast reconstruction, temporal resolution of 1.71 s was achieved for simultaneous thermometry at 1.5 × 1.5 × 2 mm voxel resolution. PRF shift thermometry was not sensitive to heating in fat. T1 changes were observed in fat at the ultrasound focus. The mean T1-temperature coefficient for fat was determined to be 1.9%/°C ± 0.2%/°C. Precision was 0.76°C ± 0.18°C for PRF shift thermometry in muscle and 1.93°C ± 0.60°C for T1-based thermometry in fat. Sonications in muscle showed an increase in T1 of 2.4%/°C ± 0.9%/°C.

CONCLUSION

The SR-VFA-SoS sequence was shown to simultaneously measure temperature change using PRF shift and T1-based methods in an in vivo model, providing thermometry for both aqueous and fat tissues.

An All-in-One Nanoheater and Optical Thermometer Fabricated from Fractal Nanoparticle Assemblies

We designed and optimized a dual-functional photothermal agent that performs as a nanoheater and real-time optical thermometer by leveraging gold nanoparticle (AuNP) self-assembly and anti-Stokes thermometry. We engineered colloidally stable fractal AuNP clusters with well-defined nanogaps to absorb strongly in the near-infrared and enhance anti-Stokes vibrational modes via surface-enhanced Raman scattering (SERS) for electromagnetic (EM) hotspot-localized thermometry during plasmonic heating. Photothermal characterization and simulations of a range of AuNP building block sizes demonstrated that 40 nm AuNPs are optimum for combined plasmonic heating and SERS due to the high probability of in resonance chains within assemblies. We explored the relationship between the far-field of our AuNP clusters and the near-field enhancement of anti-Stokes modes in the context of SERS thermometry, setting out design considerations for applying SERS thermometry. Finally, using a single near-infrared (NIR) laser source, we demonstrated plasmonic heating of a colloidal system with simultaneous accurate temperature measurement from EM hotspots via the thermal information encoded in the anti-Stokes mode of surface-bound Raman reporter molecules. Ultimately, our approach could enable real-time noninvasive temperature feedback from plasmonic nanoparticles within tumor tissue environments to guide safe and effective temperature increases during cancer photothermal therapy.

Complex-valued neural networks to speed-up MR thermometry during hyperthermia using Fourier PD and PDUNet

Hyperthermia (HT) in combination with radio- and/or chemotherapy has become an accepted cancer treatment for distinct solid tumour entities. In HT, tumour tissue is exogenously heated to temperatures between 39 and 43 °C for 60 min. Temperature monitoring can be performed non-invasively using dynamic magnetic resonance imaging (MRI). However, the slow nature of MRI leads to motion artefacts in the images due to the movements of patients during image acquisition. By discarding parts of the data, the speed of the acquisition can be increased - known as undersampling. However, due to the invalidation of the Nyquist criterion, the acquired images might be blurry and can also produce aliasing artefacts. The aim of this work was, therefore, to reconstruct highly undersampled MR thermometry acquisitions with better resolution and with fewer artefacts compared to conventional methods. The use of deep learning in the medical field has emerged in recent times, and various studies have shown that deep learning has the potential to solve inverse problems such as MR image reconstruction. However, most of the published work only focuses on the magnitude images, while the phase images are ignored, which are fundamental requirements for MR thermometry. This work, for the first time, presents deep learning-based solutions for reconstructing undersampled MR thermometry data. Two different deep learning models have been employed here, the Fourier Primal-Dual network and the Fourier Primal-Dual UNet, to reconstruct highly undersampled complex images of MR thermometry. MR images of 44 patients with different sarcoma types who received HT treatment in combination with radiotherapy and/or chemotherapy were used in this study. The method reduced the temperature difference between the undersampled MRIs and the fully sampled MRIs from 1.3 to 0.6 °C in full volume and 0.49 °C to 0.06 °C in the tumour region for a theoretical acceleration factor of 10.

EPI proton resonant frequency temperature mapping at 0.5T in the brain: Comparison to single-echo gradient recalled echo

PURPOSE

Evaluate the use of both single-echo gradient recalled echo (SE-GRE) and EPI approaches to creating temperature maps on a mid-field head-only scanner, both in vivo and on a tissue mimicking gel.

METHODS

Three 2D protocols were investigated (an SE-GRE, single-shot EPI, and an averaged single-shot EPI). The protocols used either a gradient recalled acquisition or an echo planar acquisition, with EPI parameters optimized for the longerT 2 * $$ {\mathrm{T}}_2^{\ast } $$ at lower field-strengths. Phantom experiments were conducted to evaluate temperature tracking while cooling, comparing protocol to measurements from an optical fiber thermometer. Studies were performed on a 0.5T head only MR scanner. Temperature stability maps were produced in vivo for the various protocols to evaluate precision.

RESULTS

The use of an EPI protocol for thermometry improved temperature precision in a temperature control phantom and provided an 18% improvement in temperature measurement precision in vivo. Temperature tracking using a fast (<2 s) update rate EPI thermometry sequence provided a similar precision to the slower SE-GRE protocol.

CONCLUSION

While SE-GRE PRF thermometry shows good performance, EPI methods offer improved tracking precision or update rate, making them a better option for thermometry in the brain at mid-field.

Flexible, high spatial and temporal resolution absolute thermometry of MRI phantoms using ethylene glycol

PURPOSE

Quantitative MRI (qMRI) parameters such as relaxation rates and diffusion parameters are typically temperature dependent. Therefore, using phantoms to evaluate the accuracy of qMRI pulse sequences requires accurate knowledge of the absolute temperature throughout the phantom. This work aims to evaluate the use of ethylene glycol (EG) together with a multi-echo gradient recalled echo (ME-GRE) pulse sequence for MR spectroscopic-based measurement of absolute temperature in phantoms.

METHODS

We develop and test a simple MR spectroscopic imaging approach to rapidly, automatically, accurately, and precisely measure absolute temperature at multiple locations throughout a qMRI phantom. A series of EG cubes are located throughout the phantom and imaged with a ME-GRE pulse sequence over a wide range of temperatures (˜7°C-37°C) and on multiple scanners. The resulting images were automatically processed to isolate the EG. The measured spectral peak spacing was calibrated to temperature using fiber optic probes. The accuracy and precision of the measurements were evaluated between scanners over a range of temperatures.

RESULTS

The spectral peak spacing of EG can be used to predict temperature with an error of 0.6°C ± 0.3°C with an increase of 0.1°C when evaluated on different scanners. An automatic reconstruction approach without manual input is described, highlighting the feasibility of online implementation on a clinical MRI scanner.

CONCLUSION

The described MR spectroscopic imaging approach is easy to implement and provides robust, automatic, and accurate measurement of absolute temperature throughout a phantom.

Overview of Therapeutic Ultrasound Applications and Safety Considerations: 2024 Update

A 2012 review of therapeutic ultrasound was published to educate researchers and physicians on potential applications and concerns for unintended bioeffects (doi: 10.7863/jum.2012.31.4.623). This review serves as an update to the parent article, highlighting advances in therapeutic ultrasound over the past 12 years. In addition to general mechanisms for bioeffects produced by therapeutic ultrasound, current applications, and the pre-clinical and clinical stages are outlined. An overview is provided for image guidance methods to monitor and assess treatment progress. Finally, other topics relevant for the translation of therapeutic ultrasound are discussed, including computational modeling, tissue-mimicking phantoms, and quality assurance protocols.

Focused Ultrasound: Noninvasive Image-Guided Therapy

ABSTRACT

Invasive open surgery used to be compulsory to access tumor mass to perform excision or resection. Development of minimally invasive laparoscopic procedures followed, as well as catheter-based approaches, such as stenting, endovascular surgery, chemoembolization, brachytherapy, which minimize side effects and reduce the risks to patients. Completely noninvasive procedures bring further benefits in terms of reducing risk, procedure time, recovery time, potential of infection, or other side effects. Focusing ultrasound waves from the outside of the body specifically at the disease site has proven to be a safe noninvasive approach to localized ablative hyperthermia, mechanical ablation, and targeted drug delivery. Focused ultrasound as a medical intervention was proposed decades ago, but it only became feasible to plan, guide, monitor, and control the treatment procedures with advanced radiological imaging capabilities. The purpose of this review is to describe the imaging capabilities and approaches to perform these tasks, with the emphasis on magnetic resonance imaging and ultrasound. Some procedures already are in clinical practice, with more at the clinical trial stage. Imaging is fully integrated in the workflow and includes the following: (1) planning, with definition of the target regions and adjacent organs at risk; (2) real-time treatment monitoring via thermometry imaging, cavitation feedback, and motion control, to assure targeting and safety to adjacent normal tissues; and (3) evaluation of treatment efficacy, via assessment of ablation and physiological parameters, such as blood supply. This review also focuses on sonosensitive microparticles and nanoparticles, such as microbubbles injected in the bloodstream. They enable ultrasound energy deposition down to the microvascular level, induce vascular inflammation and shutdown, accelerate clot dissolution, and perform targeted drug delivery interventions, including focal gene delivery. Especially exciting is the ability to perform noninvasive drug delivery via opening of the blood-brain barrier at the desired areas within the brain. Overall, focused ultrasound under image guidance is rapidly developing, to become a choice noninvasive interventional radiology tool to treat disease and cure patients.

A high magnetic resonance imaging (MRI) contrast agar/silica-based phantom for evaluating focused ultrasound (FUS) protocols

PURPOSE

Thermal ablation therapies require tissue mimicking phantoms for evaluating novel systems. Herein, an agar phantom exhibiting high magnetic resonance imaging (MRI) contrast to noise ratio (CNR) was developed for testing focused ultrasound (FUS) protocols.

METHODS

Four agar based phantoms (6 % w/v) were fabricated with varied silica concentrations (0, 2, 4, or 6 % w/v) and subjected to FUS inside a 3 T MRI. T2-Weighted Fast Spin Echo (T2-W FSE) images were acquired after sonications to assess the effect of varied silica on CNR of inflicted lesions. The highest CNR phantom was sonicated and its proton resonance frequency (PRF) coefficient, thermal dose denaturation threshold and ability to sustain good lesion CNR 0-44 min post exposures were assessed.

RESULTS

T2-W median lesion CNR between 1.5-453.5 was observed, exponentially increasing with increased silica concentration. High CNR was achieved with 4 % w/v silica, with the PRF coefficient of the phantom calculated at -0.00954 ppm/°C. The thermal dose denaturation threshold was revealed at 2 × 106 CEM43°C by comparing thermal dose maps with T2-W FSE lesion hyperenhancement. Progressive lesion CNR loss was observed, with CNR lost 28 min after sonications.

CONCLUSIONS

The proposed phantom possesses excellent T2-W contrast of inflicted lesions while exhibiting a tissue like PRF coefficient and can thus constitute an inexpensive reusable tool for validating FUS systems and protocols.

Evaluation of a deformable image registration algorithm for image-guided thermal ablation of liver tumors on clinically acquired MR-temperature maps

BACKGROUND

Quantitative real-time MRI-based temperature mapping techniques are hampered by abdominal motion. Intrascan motion can be reduced by rapid acquisition sequences such as 2D echo planar imaging (EPI), and inter-scan organ displacement can be compensated by image processing such as optical flow (OF) algorithms. However, motion field estimation can be seriously affected by local variation of signal intensity on magnitude images inherent to tissue heating, potentially leading to erroneous temperature estimates.

PURPOSE

This study aims to characterize, in the context of clinical MRI-guided microwave ablation (MWA), a novel deformable image registration (DIR) algorithm that enhances the generation of thermal maps aligned to a reference position, a critical step for calculating cumulative thermal dose and, consequently, for the real-time evaluation of interventional procedure progress.

METHODS

A retrospective image analysis was performed on 11 patients that underwent MWA of a liver tumor (primary or metastasis). Ablation duration was set to 9 ± 2 min with a 14-gauge large antenna. A stack of 13-20 contiguous slices was acquired dynamically (350 repetitions) at 1.5T using a single-shot EPI sequence. Evaluation was first performed on motion-free datasets (5 gated acquisitions using a cushion positioned in the patient abdomen) then with ones with motion (8 fixed-frequency acquisitions at 0.5 Hz). Temperature, thermal dose and lesion size were computed using three workflows: (i) standard phase subtraction (gold standard), (ii) conventional OF motion compensation, (iii) PCA-based OF motion compensation. The impact of flow field, temperature and lesion volume estimation were compared using averaged endpoint error (AEE), NRMSE and bland Altman plot, respectively.

RESULTS

Intensity signal decreases (close to 50%) were observed in the vicinity of the probe during MW energy delivery. Both motion correction algorithms reduce the NRMSE of magnitude images throughout the acquisition (p < 0.005) and achieve similar results between them. Gated acquisition results. Conventional OF produced erroneous vector fields compared to the PCA-based OF, leading to higher maximal EE (3 mm vs. 1 mm) and temperature errors up to 15°C-20°C. PCA-based OF algorithm significantly reduces the NRMSE of temperature (p < 0.005). The conventional OF method underestimated the final size of lesions with a bias of 0.93 cm3 while the PCA-based OF reported a bias of 0.5 cm3. Fixed frequency acquisition results. The temperature estimation without motion correction led to strong fluctuations or loss of temperature measurement while the proposed PCA-based OF recovered both a stable and precise measurement with null bias.

CONCLUSION

The deformable image registration algorithm is less sensitive to local variations of the signal. Volumetric temperature imaging without gating (20 slices/2 s) could be performed with the same accuracy, and offer trade-offs in acquisition time or volume coverage. Such a strategy is expected to increase procedure safety by monitoring large volumes more rapidly for MR-guided thermotherapy on mobile organs.

Transcranial MRI-guided Histotripsy Targeting Using MR-thermometry and MR-ARFI

OBJECTIVE

Transcranial magnetic resonance imaging (MRI)-guided histotripsy has been demonstrated to treat various locations in in vivo swine brain through a human skull. To ensure that the histotripsy treatment is delivered to the intended target location, accurate pre-treatment targeting is necessary. In this work, we investigate the feasibility of MR-thermometry and MR-acoustic radiation force imaging (MR-ARFI) to perform pre-treatment targeting of histotripsy in ex vivo bovine brain through a human skull.

METHODS

A 700 kHz, 128-element MR-compatible histotripsy array was used to generate histotripsy and tone-burst sonications. The array's electronic drivers were modified to also generate low-amplitude tone-burst sonications to perform MR-thermometry and MR-ARFI-based targeting. Twelve ex vivo bovine brains were treated with histotripsy at 35 MPa, 75 MPa and through a skull at 36 MPa. Before treating the tissue, both MR-ARFI and MR-thermometry were used to estimate the lesion location. Finally, the location of the histotripsy lesion was compared with the focus estimated by MR-thermometry and MR-ARFI.

RESULTS

MR-thermometry and MR-ARFI were able to successfully perform pre-treatment targeting of histotripsy using the modified histotripsy array driver. Histotripsy focus was estimated with mean absolute errors along the transverse/longitudinal axis of 2.06/2.95 mm and 2.13/2.51 mm for MR-ARFI and MR-thermometry, respectively. The presence of the human skull reduced the pressure at the focal region, but it did not compromise the targeting accuracy of either of the two methods with a mean absolute error of 1.10/2.91 mm and 1.29/2.91 mm for MR-ARFI and MR-thermometry, respectively.

CONCLUSION

This study demonstrated that transcranial histotripsy pre-treatment targeting is feasible with MR-thermometry and MR-ARFI.

MR-Guided Transrectal Focal Laser Ablation for Localized Low- and Intermediate-Risk Prostate Cancer: Initial Outcomes Using an Integrated Laser Ablation System

PURPOSE

To investigate the feasibility and safety of magnetic resonance (MR) imaging-guided focal laser ablation (FLA) in localized, International Society of Urological Pathology (ISUP) Grade 1-3, prostate cancer (PCa) using an integrated system.

MATERIALS AND METHODS

Ten consecutive males (mean age, 66 years [SD ± 7]) with low-to-intermediate-risk PCa were prospectively included (April 2022-May 2023) and treated with MR imaging-guided FLA using an integrated system for laser energy control and MR thermometry monitoring. Primary end points were technical success, procedure-related adverse events (AEs) following Society of Interventional Radiology (SIR) classification, and 12-month local tumor progression-free survival, defined as no evident residual/recurrent disease on follow-up imaging or histopathology at the treatment site. Secondary end points included MR imaging-derived percentual volumetric tumor coverage, prostate-specific antigen (PSA) levels, and sexual and urinary function response measured by the Sexual Health Inventory for Men (SHIM) and International Prostate Symptom Score (IPSS) index questionnaires, respectively.

RESULTS

Technical success was achieved in all 10 (100%) patients (ISUP Grade 1, n = 1; Grade 2, n = 8; and Grade 3, n = 1). Three AEs were observed: urinary tract infection (n = 2; SIR Grade 2) and acute urinary retention (n = 1; SIR Grade 3). Cumulative 12-month local tumor progression-free survival was 80% (8/10 patients). Median tumor coverage was 100% (IQR, 95%-100%). Compared with baseline, the mean PSA level decreased, but did not reach statistical significance (6.6 vs 4.4 ng/mL; P = .06), and mean urinary (8.6 vs 7.3; P = 0.60) and sexual function (11.3 vs 10.5; P = 1.00) scores were nonsignificantly altered at 12-month follow-up.

CONCLUSIONS

MR imaging-guided FLA in patients with low-to-intermediate-risk PCa using an integrated system was feasible and safe and indicated promising short-term oncologic and functional outcomes.

Magnetic Resonance Acoustic Radiation Force Imaging (MR-ARFI)

This review covers the theoretical background, pulse sequence considerations, practical implementations, and multitudes of applications of magnetic resonance acoustic radiation force imaging (MR-ARFI) described to date. MR-ARFI is an approach to encode tissue displacement caused by the acoustic radiation force of a focused ultrasound field into the phase of a MR image. The displacement encoding is done with motion encoding gradients (MEG) which have traditionally been added to spin echo-type and gradient recalled echo-type pulse sequences. Many different types of MEG (monopolar, bipolar, tripolar etc.) have been described and pros and cons are discussed. We further review studies investigating the safety of MR-ARFI, as well as approaches to simulate the MR-ARFI displacement. Lastly, MR-ARFI applications such as for focal spot localization, tissue stiffness interrogation following thermal ablation, trans-skull aberration correction, and simultaneous MR-ARFI and MR thermometry are discussed. EVIDENCE LEVEL: N/A TECHNICAL EFFICACY: Stage 1.

Image-Based Monitoring of Thermal Ablation

Thermal therapy is a commonly used local treatment technique in clinical practice. Monitoring the treatment process is essential for ensuring its success. In this review, we analyze recent image-based methods for thermal therapy monitoring, focusing particularly on their feasibility for synchronous or immediate postoperative monitoring. This includes thermography and other techniques that track the physical changes in tissue during thermal ablation. Potential directions and challenges for further clinical applications are also summarized.

Millisecond-level transient heating and temperature monitoring technique for ultrasound-induced thermal strain imaging

Background: Ultrasound-induced thermal strain imaging (US-TSI) is a promising ultrasound imaging modality that has been demonstrated in preclinical studies to identify a lipid-rich necrotic core of an atherosclerotic plaque. However, human physiological motion, e.g., cardiac pulsation, poses challenges in implementing US-TSI applications, where achieving a millisecond-level temperature rise by delivering acoustic energy from a compact US-TSI probe is a key requirement. This study aims to develop a transient ultrasound heating and thermocouple monitoring technique at the millisecond level for US-TSI applications. Methods: We designed, prototyped, and characterized a novel US-TSI probe that includes a high-power, 3.5 MHz heating transducer with symmetrical dual 1D concave array. Additionally, millisecond-level temperature monitoring was demonstrated with fast-response thermocouples in laser- and ultrasound- induced thermal tests. Subsequently, we demonstrated the prototyped US-TSI probe can produce a desired temperature rise in a millisecond-short time window in vitro phantom and in vivo animal tests. Results: The prototyped US-TSI probe delivered zero-to-peak acoustic pressure up to 6.2 MPa with a 90 VPP input voltage. Both laser- and ultrasound- induced thermal tests verified that the selected thermocouples can monitor temperature change within 50 ms. The fast-response thermocouple confirmed the transient heating ability of the US-TSI probe, achieving a 3.9 °C temperature rise after a 25 ms heating duration (50% duty cycle) in the gel phantom and a 2.0 °C temperature rise after a 50 ms heating duration (50% duty cycle) in a pig model. Conclusions: We successfully demonstrated a millisecond-level transient heating and temperature monitoring technique utilizing the novel US-TSI probe and fast-response thermocouples. The reported transient ultrasound heating and thermocouple monitoring technique is promising for future in vivo human subject studies in US-TSI or other ultrasound-related thermal investigations.

Real-time multislice MR-thermometry of the prostate: Assessment of feasibility, accuracy and sources of biases in patients

PURPOSE

The primary purpose of this study was to evaluate the accuracy of an MR-thermometry sequence for monitoring prostate temperature. The secondary purposes were to analyze clinical and technical factors that may affect accuracy and testing the method in a realistic setting, with MR-guided Laser ablation on an ex vivo muscle sample.

MATERIALS AND METHODS

An ex vivo muscle sample was subjected to Laser ablation while using a two-dimensional multislice segmented echo planar imaging sequence for MR thermometry. The MR thermometry measurements were compared with invasive sensor temperature readings to assess accuracy. Subsequently, 56 men with a median age of 70 years (age range: 53-84 years) who underwent prostate MRI examinations at 1.5- (n = 27) or 3 T (n = 24) were prospectively included. For each patient, the proportion of 'noisy voxels' (i.e., those with a temporal standard deviation of temperature [SD(T)] > 2 °C) in the prostate was calculated. The impact of clinical and technical factors on the proportion of noisy voxels was also examined.

RESULTS

MR-thermometry showed excellent correlation with invasive sensors during MR-guided Laser ablation on the ex vivo muscle sample. The median proportion of noisy voxels per patient in the entire cohort was 1 % (Q1, 0.2; Q3, 4.9; range: 0-90.4). No significant differences in median proportion of noisy voxels were observed between examinations performed at 1.5 T and those at 3 T (P = 0.89 before and after adjustment). No clinical or technical factors significantly influenced the proportion of noisy voxels.

CONCLUSION

Two-dimensional real time multislice MR-thermometry is feasible and accurate for monitoring prostate temperature in patients.

Nanoparticles based image-guided thermal therapy and temperature feedback

Nanoparticles have emerged as versatile tools in the realm of thermal therapy, offering precise control and feedback mechanisms for targeted treatments. This review explores the intersection of nanotechnology and thermal therapy, focusing on the utilization of nanoparticles for image-guided interventions and temperature monitoring. Starting with an exploration of local temperature dynamics compared to whole-body responses, we delve into the landscape of nanomaterials and their pivotal role in nanomedicine. Various physical stimuli employed in therapy and imaging are scrutinized, laying the foundation for nanothermal therapies and the accompanying challenges. A comprehensive analysis of nanomaterial architecture ensues, delineating the functionalities of magnetic, plasmonic, and luminescent nanomaterials within the context of thermal therapies. Nano-design intricacies, including core-shell structures and monodisperse properties, are dissected for their impact on therapeutic efficacy. Furthermore, considerations in designing in vivo nanomaterials, such as hydrodynamic radii and core sizes at sub-tissue levels, are elucidated. The review then delves into specific modalities of thermally induced therapy, including magnetically induced hyperthermia and luminescent-based thermal treatments. Magnetic hyperthermia treatment is explored alongside its imaging and relaxometric properties, emphasizing the implications of imaging formulations on biotransformation and biodistribution. This review also provides an overview of the magnetic hyperthermia treatment using magnetic nanoparticles to induce localized heat in tissues. Similarly, optical and thermal imaging techniques utilizing luminescent nanomaterials are discussed, highlighting their potential for light-induced thermal therapy and cellular-level temperature monitoring. Finally, the application landscape of diagnosis and photothermal therapy (PTT) is surveyed, encompassing diverse areas such as cancer treatment, drug delivery, antibacterial therapy, and immunotherapy. The utility of nanothermometers in elucidating thermal relaxation dynamics as a diagnostic tool is underscored, alongside discussions on PTT hyperthermia protocols. Moreover, the advancements in nanoparticle magnetic imaging and implications of imaging formulations especially in creating positive MRI contrast agents are also presented. This comprehensive review offers insights into the evolving landscape of nanoparticle-based image-guided thermal therapies, promising advancements in precision medicine and targeted interventions, underscoring the importance of continued research in optimization for the full potential of magnetic hyperthermia to improve its efficacy and clinical translation.

Patient-Specific Prediction of Immediate Phase Lung Microwave Ablation Zone Size

PURPOSE

To investigate the effect of patient and tumor-specific characteristics on the size of immediate phase lung microwave ablation (MWA) zone and establish a prediction model.

MATERIALS AND METHODS

This institutional review board (IRB)-approved, Health Insurance Portability and Accountability Act (HIPAA)-compliant cohort included 164 lesions from 99 patients who underwent computed tomography (CT)-guided lung MWA, and the 2-dimensional elliptical ground-glass opacity ablation zone was measured. Duration, maximum temperature, tumor depth, presence of emphysema, history of ipsilateral lung ablation, surgery, and radiotherapy were recorded. K-fold cross validation with k = 5 and Least Absolute Shrinkage and Selection Operator were used to build prediction models for the major and minor axes and area of the ablation zone.

RESULTS

The median of immediate phase ablation duration was 2 minutes (interquartile range, 1.5-4.25 minutes) with 65 W of power for all ablations. The mean major and minor axes and area of ablation zone were 3.1 cm (SD ± 0.6), 2.0 cm (SD ± 0.5), and 5.1 cm2 (SD ± 2.1), respectively. The major and minor axes and area of immediate phase ablation zone dimensions were significantly associated with duration (P < .001, P < .001, and P < .001, respectively), maximum temperature (P < .001, P < .001, and P < .001, respectively), tumor depth (P = .387, P < .001, and P < .001, respectively), history of ipsilateral lung ablation (P = .008, P = .286, and P = .076, respectively), and lung radiotherapy (P = .001, P = .042, and P = .015, respectively). The prediction model showed R2 values for major and minor axes and area of the ablation zone to be 0.50, 0.45, and 0.53, respectively.

CONCLUSIONS

Duration of ablation, maximum temperature, tumor depth, history of ipsilateral lung ablation, surgery, and radiotherapy were significantly associated with the ablation zone dimensions and size and can be used to build the prediction model to approximate the immediate phase lung MWA zone.

MR Imaging-Guided Focused Ultrasound-Clinical Applications in Managing Malignant Gliomas

Malignant gliomas (MGs) are the most common primary brain tumors in adults. Despite recent advances in understanding the biology and potential therapeutic vulnerabilities of MGs, treatment options remain limited as the delivery of drugs is often impeded by the blood-brain barrier (BBB), and safe, complete surgical resection may not always be possible, especially for deep-seated tumors. In this review, the authors highlight emerging applications for MR imaging-guided focused ultrasound (MRgFUS) as a noninvasive treatment modality for MGs. Specifically, the authors discuss MRgFUS's potential role in direct tumor cell killing, opening the BBB, and modulating antitumor immunity.

Real-time non-invasive control of ultrasound hyperthermia using high-frequency ultrasonic backscattered energy inex vivotissue andin vivoanimal studies

Objective.A reliable, calibrated, non-invasive thermometry method is essential for thermal therapies to monitor and control the treatment. Ultrasound (US) is an effective thermometry modality due to its relatively high sensitivity to temperature changes, and fast data acquisition and processing capabilities.Approach.In this work, the change in backscattered energy (CBE) was used to control the tissue temperature non-invasively using a real-time proportional-integral-derivative (PID) controller. A clinical high-frequency US scanner was used to acquire radio-frequency echo data fromex vivoporcine tissue samples andin vivomice hind leg tissue while the tissue was treated with mild hyperthermia by a focused US applicator. The PID controller maintained the focal temperature at approximately 40 °C for about 4 min.Main results.The results show that the US thermometry based on CBE estimated by a high-frequency US scanner can produce 2D temperature maps of a localized heating region and to estimate the focal temperature during mild hyperthermia treatments. The CBE estimated temperature varied by an average of ±0.85 °C and ±0.97 °C, compared to a calibrated thermocouple, inex vivoandin vivostudies, respectively. The mean absolute deviations of CBE thermometry during the controlled hyperthermia treatment were ±0.45 °C and ±0.54 °C inex vivoandin vivo,respectively.Significance.It is concluded that non-invasive US thermometry via backscattered energies at high frequencies can be used for real-time monitoring and control of hyperthermia treatments with acceptable accuracy. This provides a foundation for an US mediated drug delivery system.

Magnetic Resonance Thermometry of Focused Ultrasound Using a Preclinical Focused Ultrasound Robotic System at 3T

AIM

Focused ultrasound (FUS) therapies are often performed within magnetic resonance imaging (MRI) systems providing thermometry-based temperature monitoring. Herein, MRI thermometry was assessed for FUS sonications executed using a preclinical system on agar-based phantoms at 1.5T and 3T MRI scanners, using the proton resonance frequency shift technique.

MATERIALS AND METHODS

Sonications were executed at 1.5T and 3T to assess the system and observe variations in magnetic resonance (MR) thermometry temperature measurements. MR thermometry was assessed at 3T, for identical sonications on three agar-based phantoms doped with varied silica and evaporated milk concentrations, and for sonications executed at varied acoustic power of 1.5-45 W. Moreover, echo time (TE) values of 5-20 ms were used to assess the effect on the signal-to-noise ratio (SNR) and temperature change sensitivity.

RESULTS

Clearer thermal maps with a 2.5-fold higher temporal resolution were produced for sonications at 3T compared to 1.5T, despite employment of similar thermometry sequences. At 3T, temperature changes between 41°C and 50°C were recorded for the three phantoms produced with varied silica and evaporated milk, with the addition of 2% w/v silica resulting in a 20% increase in temperature change. The lowest acoustic power that produced reliable beam detection within a voxel was 1.5 W. A TE of 10 ms resulted in the highest temperature sensitivity with adequate SNR.

CONCLUSIONS

MR thermometry performed at 3T achieved short temporal resolution with temperature dependencies exhibited with the sonication and imaging parameters. Present data could be used in preclinical MRI-guided FUS feasibility studies to enhance MR thermometry.

Therapeutic ultrasound transducer technology and monitoring techniques: a review with clinical examples

The exponential growth of therapeutic ultrasound applications demonstrates the power of the technology to leverage the combinations of transducer technology and treatment monitoring techniques to effectively control the preferred bioeffect to elicit the desired clinical effect.Objective: This review provides an overview of the most commonly used bioeffects in therapeutic ultrasound and describes existing transducer technologies and monitoring techniques to ensure treatment safety and efficacy.Methods and materials: Literature reviews were conducted to identify key choices that essential in terms of transducer design, treatment parameters and procedure monitoring for therapeutic ultrasound applications. Effective combinations of these options are illustrated through descriptions of several clinical indications, including uterine fibroids, prostate disease, liver cancer, and brain cancer, that have been successful in leveraging therapeutic ultrasound to provide effective patient treatments.Results: Despite technological constraints, there are multiple ways to achieve a desired bioeffect with therapeutic ultrasound in a target tissue. Visualizations of the interplay of monitoring modality, bioeffect, and applied acoustic parameters are presented that demonstrate the interconnectedness of the field of therapeutic ultrasound. While the clinical indications explored in this review are at different points in the clinical evaluation path, based on the ever expanding research being conducted in preclinical realms, it is clear that additional clinical applications of therapeutic ultrasound that utilize a myriad of bioeffects will continue to grow and improve in the coming years.Conclusions: Therapeutic ultrasound will continue to improve in the next decades as the combination of transducer technology and treatment monitoring techniques will continue to evolve and be translated in clinical settings, leading to more personalized and efficient therapeutic ultrasound mediated therapies.

2D ultrasound thermometry during thermal ablation with high-intensity focused ultrasound

The clinical use of high intensity focused ultrasound (HIFU) therapy for noninvasive tissue ablation has recently gained momentum. Guidance is provided by either magnetic resonance imaging (MRI) or conventional B-mode ultrasound imaging, each with its own advantages and disadvantages. The main limitation of ultrasound imaging is its inability to provide temperature measurements over the ranges corresponding to the target temperatures during ablative thermal therapies (between 55 °C and 70 °C). Here, variations in ultrasound backscattered energy (ΔBSE) were used to monitor temperature increases in liver tissue up to an absolute value of 90 °C during and after HIFU treatment. In vitro experimental measurements were performed in 47 bovine liver samples using a toroidal HIFU transducer operating at 2.5 MHz to increase the temperature of tissues. An ultrasound imaging probe working at 7.5 MHz was placed in the center of the HIFU transducer to monitor the backscattered signals. The free-field acoustic power was set to 9 W, 12 W or 16 W in the different experiments. HIFU sonications were performed for 240 s using a duty cycle of 83 % to allow ultrasound imaging and raw radiofrequency data acquisition during exposures. Measurements showed a linear relationship between ΔBSE (in dB) and temperature (r = 0.94, p < 0.001) over a temperature range from 37 °C to 90 °C, with a high reliability of temperature measurements below 75 °C. Monitoring can be performed at the frame rate of ultrasound imaging scanners with an accuracy within an acceptable threshold of 5 °C, given the temperatures targeted during thermal ablations. If the maximum temperature reached is below 70 °C, ΔBSE is also a reliable approach for estimating the temperature during cooling. Histological analysis shown the impact of the treatment on the spatial arrangement of cells that can explain the observed variation of ΔBSE. These results demonstrate the ability of ΔBSE measurements to estimate temperature in ultrasound images within an effective therapeutic range. This method can be implemented clinically and potentially applied to other thermal-based therapies.

Technical advances in motion-robust MR thermometry

Proton resonance frequency shift (PRFS) MR thermometry is the most common method used in clinical thermal treatments because of its fast acquisition and high sensitivity to temperature. However, motion is the biggest obstacle in PRFS MR thermometry for monitoring thermal treatment in moving organs. This challenge arises because of the introduction of phase errors into the PRFS calculation through multiple methods, such as image misregistration, susceptibility changes in the magnetic field, and intraframe motion during MRI acquisition. Various approaches for motion correction have been developed for real-time, motion-robust, and volumetric MR thermometry. However, current technologies have inherent trade-offs among volume coverage, processing time, and temperature accuracy. These tradeoffs should be considered and chosen according to the thermal treatment application. In hyperthermia treatment, precise temperature measurements are of increased importance rather than the requirement for exceedingly high temporal resolution. In contrast, ablation procedures require robust temporal resolution to accurately capture a rapid temperature rise. This paper presents a comprehensive review of current cutting-edge MRI techniques for motion-robust MR thermometry, and recommends which techniques are better suited for each thermal treatment. We expect that this study will help discern the selection of motion-robust MR thermometry strategies and inspire the development of motion-robust volumetric MR thermometry for practical use in clinics.

Fusion and Validation Method for Laser Interstitial Thermal Therapy Simulation Model and MRI

In this paper, we present a simulation program based on Monte Carlo simulation and bio-heat transfer models for laser interstitial thermal therapy (LITT). Additionally, we present a data fusion strategy that synchronizes MRI data with simulation results. The simulation model's validity was checked and proven by using this method to combine simulation calculations with MRI-measured temperature data from cases of brain lesions. For LITT treatment planning, the alignment fusion method and higher temporal-spatial resolution simulation model developed in this paper provide a clear three-dimensional visualization of the ablation temperature field displayed in situ on MRI images. This enables the preoperative planning of LITT to be conducted efficiently.

Estimation of the Proton Resonance Frequency Coefficient in Agar-based Phantoms

AIM

Agar-based phantoms are popular in high intensity focused ultrasound (HIFU) studies, with magnetic resonance imaging (MRI) preferred for guidance since it provides temperature monitoring by proton resonance frequency (PRF) shift magnetic resonance (MR) thermometry. MR thermometry monitoring depends on several factors, thus, herein, the PRF coefficient of agar phantoms was estimated.

MATERIALS AND METHODS

Seven phantoms were developed with varied agar (2, 4, or 6% w/v) or constant agar (6% w/v) and varied silica concentrations (2, 4, 6, or 8% w/v) to assess the effect of the concentration on the PRF coefficient. Each phantom was sonicated using varied acoustical power for a 30 s duration in both a laboratory setting and inside a 3T MRI scanner. PRF coefficients were estimated through linear trends between phase shift acquired using gradient sequences and thermocouple-based temperatures changes.

RESULTS

Linear regression (R 2 = 0.9707-0.9991) demonstrated a proportional dependency of phase shift with temperature change, resulting in PRF coefficients between -0.00336 ± 0.00029 and -0.00934 ± 0.00050 ppm/°C for the various phantom recipes. Weak negative linear correlations of the PRF coefficient were observed with increased agar. With silica concentrations, the negative linear correlation was strong. For all phantoms, calibrated PRF coefficients resulted in 1.01-3.01-fold higher temperature changes compared to the values calculated using a literature PRF coefficient.

CONCLUSIONS

Phantoms developed with a 6% w/v agar concentration and doped with 0%-8% w/v silica best resemble tissue PRF coefficients and should be preferred in HIFU studies. The estimated PRF coefficients can result in enhanced MR thermometry monitoring and evaluation of HIFU protocols.

Nanomaterial-Enabled Photothermal Heating and Its Use for Cancer Therapy via Localized Hyperthermia

Photothermal therapy (PTT), which employs nanoscale transducers delivered into a tumor to locally generate heat upon irradiation with near-infrared light, shows great potential in killing cancer cells through hyperthermia. The efficacy of such a treatment is determined by a number of factors, including the amount, distribution, and dissipation of the generated heat, as well as the type of cancer cell involved. The amount of heat generated is largely controlled by the number of transducers accumulated inside the tumor, the absorption coefficient and photothermal conversion efficiency of the transducer, and the irradiance of the light. The efficacy of treatment depends on the distribution of the transducers in the tumor and the penetration depth of the light. The vascularity and tissue thermal conduction both affect the dissipation of heat and thereby the distribution of temperature. The successful implementation of PTT in the clinic setting critically depends on techniques for real-time monitoring and management of temperature.

MR Imaging-Guided Prostate Cancer Therapies

Prostate cancer is the most common malignancy diagnosed in men. MR imaging-guided therapies for prostate cancer have become an increasingly common treatment alternative to traditional whole-gland therapies, such as radical prostatectomy or radiation therapy. This is especially true in men with localized, low- to intermediate-risk prostate cancer. Although long-term oncologic data remain limited, the authors describe several MR imaging-guided therapeutic options for the treatment of prostate cancer, including cryoablation, laser ablation, transrectal high-intensity focused ultrasound, and transurethral ultrasound ablation.

MR Thermometry during Transcranial MR Imaging-Guided Focused Ultrasound Procedures: A Review

Interest in transcranial MR imaging-guided focused ultrasound procedures has recently grown. These incisionless procedures enable precise focal ablation of brain tissue using real-time monitoring by MR thermometry. This article will provide an updated review on clinically applicable technical underpinnings and considerations of proton resonance frequency MR thermometry, the most common clinically used MR thermometry sequence.

Spatial patterning of energy metabolism during tissue morphogenesis

Biophysical signaling organizes forces to drive tissue morphogenesis, a process co-opted during disease progression. The systematic buildup of forces at the tissue scale is energetically demanding. Just as mechanical forces, gene expression, and concentrations of morphogens vary spatially across a developing tissue, there might similarly be spatial variations in energy consumption. Recent studies have started to uncover the connections between spatial patterns of mechanical forces and spatial patterns of energy metabolism. Here, we define and review the concept of energy metabolism during tissue morphogenesis. We highlight experiments showing spatial variations in energy metabolism across several model systems, categorized by morphogenetic motif, including convergent extension, branching, and migration. Finally, we discuss approaches to further enable quantitative measurements of energy production and consumption during morphogenesis.

Aqueous Dispersion of Manganese-Zinc Ferrite Nanoparticles Protected by PEG as a T2 MRI Temperature Contrast Agent

Mixed manganese-zinc ferrite nanoparticles coated with PEG were studied for their potential usefulness in MRI thermometry as temperature-sensitive contrast agents. Particles in the form of an 8.5 nm core coated with a 3.5 nm layer of PEG were fabricated using a newly developed, one-step method. The composition of Mn0.48Zn0.46Fe2.06O4 was found to have a strong thermal dependence of magnetization in the temperature range between 5 and 50 °C. Nanoparticles suspended in an agar gel mimicking animal tissue and showing non-significant impact on cell viability in the biological test were studied with NMR and MRI over the same temperature range. For the concentration of 0.017 mg/mL of Fe, the spin-spin relaxation time T2 increased from 3.1 to 8.3 ms, while longitudinal relaxation time T1 shows a moderate decrease from 149.0 to 125.1 ms. A temperature map of the phantom exposed to the radial temperature gradient obtained by heating it with an 808 nm laser was calculated from T2 weighted spin-echo differential MR images. Analysis of temperature maps yields thermal/spatial resolution of 3.2 °C at the distance of 2.9 mm. The experimental relaxation rate R2 data of water protons were compared with those obtained from calculations using a theoretical model incorporating the motion averaging regime.

Organic Fluorine Compounds and Their Uses as Molecular MR-Based Temperature Sensors

The interest in fluorinated substances has increased significantly in recent decades due to their diverse properties and possible uses. An important analytical method in this context is NMR spectroscopy, which provides information on the structure as well as on intermolecular interactions or generally on changes in the environment of the nucleus under consideration. A physical quantity that is of great importance in most studies is temperature. However, this is not always easy, e. g. in shielded systems or within an organism. However, the application potential in chemical reactors or in medical diagnosis and therapy is very high and for this reason 13 fluorinated organic compound were chosen for a first 19 F NMR signal temperature sensitivity examination for determination of local temperatures in solution. Polyfluorinated molecules with separate 19 F MR signals are particularly suitable for temperature determination. Those can be serve as internal error-correcting thermometers without the need of a reference substance. Under these conditions, a 19 F MR signal shift of up to 0.03 ppm/K was detectable. Fluorine position and chemical environment were very important for the temperature sensitivity.

Straightforward Magnetic Resonance Temperature Measurements Combined with High Frame Rate and Magnetic Susceptibility Correction

Proton resonance frequency shift (PRFS) is an MRI-based simple temperature mapping method that exhibits higher spatial and temporal resolution than temperature mapping methods based on T1 relaxation time and diffusion. PRFS temperature measurements are validated against fiber-optic thermal sensors (FOSs). However, the use of FOSs may introduce temperature errors, leading to both underestimation and overestimation of PRFS measurements, primarily due to material susceptibility changes caused by the thermal sensors. In this study, we demonstrated susceptibility-corrected PRFS (scPRFS) with a high frame rate and accuracy for suitably distributed temperatures. A single-echo-based background removal technique was employed for phase variation correction, primarily owing to magnetic susceptibility, which enabled fast temperature mapping. The scPRFS was used to validate the temperature fidelity by comparing the temperatures of fiber-optic sensors and conventional PRFS through phantom-mimicked human and ex vivo experiments. This study demonstrates that scPRFS measurements in agar-gel are in good agreement with the thermal sensor readings, with a root mean square error (RMSE) of 0.33-0.36 °C in the phantom model and 0.12-0.16 °C in the ex vivo experiment. These results highlight the potential of scPRFS for precise thermal monitoring and ablation in both low- and high-temperature non-invasive therapies.

Modeling and in vivo experimental validation of 1,064 nm laser interstitial thermal therapy on brain tissue

Introduction

Laser interstitial thermal therapy (LITT) at 1064 nm is widely used to treat epilepsy and brain tumors; however, no numerical model exists that can predict the ablation region with careful in vivo validation.

Methods

In this study, we proposed a model with a system of finite element methods simulating heat transfer inside the brain tissue, radiative transfer from the applicator into the brain tissue, and a model for tissue damage.

Results

To speed up the computation for practical applications, we also validated P1-approximation as an efficient and fast method for calculating radiative transfer by comparing it with Monte Carlo simulation. Finally, we validated the proposed numerical model in vivo on six healthy canines and eight human patients with epilepsy and found strong agreement between the predicted temperature profile and ablation area and the magnetic resonance imaging-measured results.

Discussion

Our results demonstrate the feasibility and reliability of the model in predicting the ablation area of 1,064 nm LITT, which is important for presurgical planning when using LITT.

Magnetic Resonance Acoustic Radiation Force Imaging (MR-ARFI) for the monitoring of High Intensity Focused Ultrasound (HIFU) ablation in anisotropic tissue

OBJECTIVE

We introduce a non-invasive MR-Acoustic Radiation Force Imaging (ARFI)-based elastography method that provides both the local shear modulus and temperature maps for the monitoring of High Intensity Focused Ultrasound (HIFU) therapy.

MATERIALS AND METHODS

To take tissue anisotropy into account, the local shear modulus μ is determined in selected radial directions around the focal spot by fitting the phase profiles to a linear viscoelastic model, including tissue-specific mechanical relaxation time τ. MR-ARFI was evaluated on a calibrated phantom, then applied to the monitoring of HIFU in a gel phantom, ex vivo and in vivo porcine muscle tissue, in parallel with MR-thermometry.

RESULTS

As expected, the shear modulus polar maps reflected the isotropy of phantoms and the anisotropy of muscle. In the HIFU monitoring experiments, both the shear modulus polar map and the thermometry map were updated with every pair of MR-ARFI phase images acquired with opposite MR-ARFI-encoding. The shear modulus was found to decrease (phantom and ex vivo) or increase (in vivo) during heating, before remaining steady during the cooling phase. The mechanical relaxation time, estimated pre- and post-HIFU, was found to vary in muscle tissue.

DISCUSSION

MR-ARFI allowed for monitoring of viscoelasticity changes around the HIFU focal spot even in anisotropic muscle tissue.

Advanced software for MRgFUS treatment planning

BACKGROUND AND OBJECTIVES

Herein, a user-friendly software platform for 3-dimensional Focused Ultrasound treatment planning based on Magnetic Resonance Imaging (MRI) images is presented.

METHODS

The software directly retrieves and loads MRI images. Various design tools can be used on the MRI images to define the treatment area and the sonication parameters. Based on the treatment plan, the software controls the robotic motion and motion pattern of Magnetic Resonance guided Focused Ultrasound (MRgFUS) robotic systems to execute the treatment procedure. Real-time treatment monitoring is achieved through MRI images and thermometry. The software's functionality and performance were evaluated in both laboratory and MRI environments. Different treatment plans were designed on MRI images and sonications were executed on agar-based phantoms and polymer films.

RESULTS

Magnetic Resonance (MR) thermometry maps were acquired in the agar-based phantoms. An exceptional agreement was observed between the software-planned treatment area and the lesions produced on the polymer films.

CONCLUSIONS

The developed software was successfully integrated with the MRI and robotic system controls for performing accurate treatment planning and real-time monitoring during sonications. The software provides an extremely user-friendly interface, while in the future it could be enhanced by providing dynamic modulation of the ultrasonic parameters during the treatment process.

19 F VT NMR: Novel Tm3+ and Ce3+ Complexes Provide New Insight into Temperature Measurement Using Molecular Sensors

In the past few decades, magnetic resonance spectroscopy (MRS) and MR imaging (MRI) have developed into a powerful non-invasive tool for medical diagnostic and therapy. Especially 19 F MR shows promising potential because of the properties of the fluorine atom and the negligible background signals in the MR spectra. The detection of temperature in a living organism is quite difficult, and usually external thermometers or fibers are used. Temperature determination via MRS needs temperature-sensitive contrast agents. This article reports first results of solvent and structural influences on the temperature sensitivity of 19 F NMR signals of chosen molecules. By using this chemical shift sensitivity, a local temperature can be determined with a high precision. Based on this preliminary study, we synthesized five metal complexes and compared the results of all variable temperature measurements. It is shown that the highest 19 F MR signal temperature dependence is detectable for a fluorine nucleus in a Tm3+ -complex.

Magnetic resonance fingerprinting based thermometry (MRFT): application toex vivoimaging near DBS leads

The purpose of this study is to demonstrate the first work ofT1-based magnetic resonance thermometry using magnetic resonance fingerprinting (dubbed MRFT). We compared temperature estimation of MRFT with proton resonance frequency shift (PRFS) thermometry onex vivobovine muscle. We demonstrated MRFT's feasibility in predicting temperature onex vivobovine muscles with deep brain stimulation (DBS) lead.B0maps generated from MRFT were compared with gold standardB0maps near the DBS lead. MRFT and PRFS estimated temperatures were compared in the presence of motion. All experiments were performed on a 3 Tesla whole-body GE Premier system with a 21-channel receive head coil (GE Healthcare, Milwaukee, WI). Four fluoroptic probes were used to measure the temperature at the center of a cold muscle (probe 1), the room temperature water bottle (probe 2), and the center and periphery of the heated muscle (probes 3 and 4). We selected regions of interest (ROIs) around the location of the probes and used simple linear regression to generate the temperature sensitivity calibration equations that convertT1maps and Δsmaps to temperature maps. We then repeated the same setup and compared MRFT and PRFS thermometry temperature estimation with gold standard probe measurements. For the MRFT experiment on DBS lead, we taped the probe to the tip of the DBS lead and used a turbo spin echo sequence to induce heating near the lead. We selected ROIs around the tip of the lead to compare MRFT temperature estimation with probe measurements and compared with PRFS temperature estimation. Vendor-suppliedB0mapping sequence was acquired to compare with MRFT-generatedB0maps. We found strong linear relationships (R2> 0.958) betweenT1and temperature and Δsand temperatures in our temperature sensitivity calibration experiment. MRFT and PRFS thermometry both accurately predict temperature (RMSE < 1.55 °C) compared to probe measurements. MRFT estimated temperature near DBS lead has a similar trend as the probe temperature. BothB0maps show inhomogeneities around the lead. MRFT estimated temperature is less sensitive to motion.

Magnetic Resonance Imaging of Focused Ultrasound Radiation Force Strain Fields for Discrimination of Solid and Liquid Phases

OBJECTIVE

Focused ultrasound (FUS) has become a non-invasive option for some surgical procedures, including tumor ablation and thalamotomy. Extension of magnetic resonance (MR) imaging-guided focused ultrasound for ablation of slowly perfused cerebrovascular lesions requires a novel treatment monitoring method that does not rely on thermometry or high-frequency Doppler methods. The goal of this study was to evaluate the sensitivity and specificity of strain estimates based on MR acoustic radiation force imaging (MR-ARFI) for differentiation of solids and liquids.

METHODS

Strain fields were estimated in gelatin-based tissue-mimicking focused ultrasound phantoms on the basis of apparent displacement fields measured by MR-ARFI. MR-ARFI and diffusion-weighted imaging (DWI) measurements were made before and after FUS-induced heating to evaluate the performance of displacement, strain and apparent diffusion coefficient (ADC) measurements for the discrimination of solid and liquid phases.

RESULTS

As revealed by receiver operating characteristic analyses, axial normal strain and shear strain components performed significantly better than axial displacement measurements alone when predicting whether a gelatin had melted. Additional measurements must be made to estimate certain strain components, so this trade-off must be considered when developing clinical strategies. ADC had the best overall performance, but DWI is vulnerable to signal dropouts and susceptibility artifacts near cerebrovascular lesions, so this metric may have limited clinical applicability.

CONCLUSION

Strain components based on MR-ARFI apparent displacement measurements perform better than apparent displacement measurements alone at discriminating between solids and liquids. These methods are applicable to FUS treatment monitoring and evaluation of mechanical tissue properties in vivo.

Simultaneous proton resonance frequency T1 - MR shear wave elastography for MR-guided focused ultrasound multiparametric treatment monitoring

PURPOSE

To develop an efficient MRI pulse sequence to simultaneously measure multiple parameters that have been shown to correlate with tissue nonviability following thermal therapies.

METHODS

A 3D segmented EPI pulse sequence was used to simultaneously measure proton resonance frequency shift (PRFS) MR thermometry (MRT), T1 relaxation time, and shear wave velocity induced by focused ultrasound (FUS) push pulses. Experiments were performed in tissue mimicking gelatin phantoms and ex vivo bovine liver. Using a carefully designed FUS triggering scheme, a heating duty cycle of approximately 65% was achieved by interleaving FUS ablation pulses with FUS push pulses to induce shear waves in the tissue.

RESULTS

In phantom studies, temperature increases measured with PRFS MRT and increases in T1 correlated with decreased shear wave velocity, consistent with material softening with increasing temperature. During ablation in ex vivo liver, temperature increase measured with PRFS MRT initially correlated with increasing T1 and decreasing shear wave velocity, and after tissue coagulation with decreasing T1 and increasing shear wave velocity. This is consistent with a previously described hysteresis in T1 versus PRFS curves and increased tissue stiffness with tissue coagulation.

CONCLUSION

An efficient approach for simultaneous and dynamic measurements of PRSF, T1 , and shear wave velocity during treatment is presented. This approach holds promise for providing co-registered dynamic measures of multiple parameters, which correlates to tissue nonviability during and following thermal therapies, such as FUS.

Magnetic Resonance Thermometry Targeting for Magnetic Resonance-Guided Histotripsy Treatments

OBJECTIVE

The potential of transcranial magnetic resonance (MR)-guided histotripsy for brain applications has been described in prior in vivo studies in the swine brain through an excised human skull. The safety and accuracy of transcranial MR-guided histotripsy (tcMRgHt) rely on pre-treatment targeting guidance. In the work described here, we investigated the feasibility and accuracy of using ultrasound-induced low-temperature heating and MR thermometry for histotripsy pre-treatment targeting in ex vivo bovine brain.

METHODS

A 15-element, 750-kHz MRI-compatible ultrasound transducer with modified drivers that can deliver both low-temperature heating and histotripsy acoustic pulses was used to treat seven bovine brain samples. The samples were first heated to an approximately 1.6°C temperature increase at the focus, and MR thermometry was used to localize the target. Once the targeting was confirmed, a histotripsy lesion was generated at the focus and visualized on post-histotripsy MR images.

DISCUSSION

The accuracy of MR thermometry targeting was evaluated with the mean/standard deviation of the difference between the locus of peak heating identified by MR thermometry and the center of mass of the post-treatment histotripsy lesion, which was 0.59/0.31 mm and 1.31/0.93 mm in the transverse and longitudinal directions, respectively.

CONCLUSION

This study determined that MR thermometry could provide reliable pre-treatment targeting for transcranial MR-guided histotripsy treatment.

Advanced Radio Frequency Applicators for Thermal Magnetic Resonance Theranostics of Brain Tumors

Thermal Magnetic Resonance (ThermalMR) is a theranostic concept that combines diagnostic magnetic resonance imaging (MRI) with targeted thermal therapy in the hyperthermia (HT) range using a radiofrequency (RF) applicator in an integrated system. ThermalMR adds a therapeutic dimension to a diagnostic MRI device. Focused, targeted RF heating of deep-seated brain tumors, accurate non-invasive temperature monitoring and high-resolution MRI are specific requirements of ThermalMR that can be addressed with novel concepts in RF applicator design. This work examines hybrid RF applicator arrays combining loop and self-grounded bow-tie (SGBT) dipole antennas for ThermalMR of brain tumors, at magnetic field strengths of 7.0 T, 9.4 T and 10.5 T. These high-density RF arrays improve the feasible transmission channel count, and provide additional degrees of freedom for RF shimming not afforded by using dipole antennas only, for superior thermal therapy and MRI diagnostics. These improvements are especially relevant for ThermalMR theranostics of deep-seated brain tumors because of the small surface area of the head. ThermalMR RF applicators with the hybrid loop+SGBT dipole design outperformed applicators using dipole-only and loop-only designs, with superior MRI performance and targeted RF heating. Array variants with a horse-shoe configuration covering an arc (270°) around the head avoiding the eyes performed better than designs with 360° coverage, with a 1.3 °C higher temperature rise inside the tumor while sparing healthy tissue. Our EMF and temperature simulations performed on a virtual patient with a clinically realistic intracranial tumor provide a technical foundation for implementation of advanced RF applicators tailored for ThermalMR theranostics of brain tumors.

Improved MR temperature imaging at 0.5 T using view-sharing accelerated multiecho thermometry for MR-guided laser interstitial thermal therapy

The aim of the current study was to improve temperature-monitoring precision using multiecho proton resonance frequency shift-based thermometry with view-sharing acceleration for MR-guided laser interstitial thermal therapy (MRgLITT) on a 0.5-T low-field MR system. Both precision and speed of the temperature measurement for clinical MRgLITT treatments suffer at low field, due to reduced image signal-to-noise ratio (SNR), decreased temperature-induced phase changes, and limited RF receiver channels. In this work, a bipolar multiecho gradient-recalled echo sequence with a temperature-to-noise ratio optimal weighted echo combination is applied to improve the temperature precision. A view-sharing-based approach is utilized to accelerate signal acquisitions while preserving image SNRs. The method was evaluated using ex vivo (pork and pig brain) LITT heating experiments and in vivo (human brain) nonheating experiments on a high-performance 0.5-T scanner. In terms of results, (1) after echo combination, multiecho thermometry (i.e., ~7.5-40.5 ms, 7 TEs) provides ~1.5-1.9 times higher temperature precision than the no echo combination case (i.e., TE7 = 40.5 ms) within the same readout bandwidth. Additionally, echo registration is necessary for the bipolar multiecho sequence; (2) for a threefold acceleration, the view-sharing approach with variable-density subsampling shows around 1.8 times lower temperature errors than the GRAPPA method. Particularly for view-sharing, variable-density subsampling performs better than Interleave subsampling; and (3) ex vivo heating and in vivo nonheating experiments demonstrated that the temperature accuracy was less than 0.5 ° C $$ {}^{{}^{\circ}}\mathrm{C} $$ and that the temperature precision was less than 0.6 ° C $$ {}^{{}^{\circ}}\mathrm{C} $$ using the proposed 0.5-T thermometry. It was concluded that view-sharing accelerated multiecho thermometry is a practical temperature measurement approach for MRgLITT at 0.5 T.

MR-guided ultrasound-stimulated microbubble therapy enhances radiation-induced tumor response

High intensity focused ultrasound (HIFU) systems have been approved for therapeutic ultrasound delivery to cause tissue ablation or induced hyperthermia. Microbubble agents have also been used in combination with sonication exposures. These require temperature feedback and monitoring to prevent unstable cavitation and prevent excess tissue heating. Previous work has utilized lower power and pressure to oscillate microbubbles and transfer energy to endothelial cells in the absence of thermally induced damage that can radiosensitize tumors. This work investigated whether reduced acoustic power and pressure on a commercial available MR-integrated HIFU system could result in enhanced radiation-induced tumor response after exposure to ultrasound-stimulated microbubbles (USMB) therapy. A commercially available MR-integrated HIFU system was used with a hyperthermia system calibration provided by the manufacturer. The ultrasound transducer was calibrated to reach a peak negative pressure of - 750 kPa. Thirty male New Zealand white rabbits bearing human derived PC3 tumors were grouped to receive no treatment, 14 min of USMB, 8 Gy of radiation in a separate irradiation cabinet, or combined treatments. In vivo temperature changes were collected using MR thermometry at the tumor center and far-field muscle region. Tissues specimens were collected 24 h post radiation therapy. Tumor cell death was measured and compared to untreated controls through hematoxylin and eosin staining and immunohistochemical analysis. The desired peak negative pressure of - 750 kPa used for previous USMB occurred at approximately an input power of 5 W. Temperature changes were limited to under 4 °C in ten of twelve rabbits monitored. The median temperature in the far-field muscle region of the leg was 2.50 °C for groups receiving USMB alone or in combination with radiation. Finally, statistically significant tumor cell death was demonstrated using immunohistochemical analysis in the combined therapy group compared to untreated controls. A commercial MR-guided therapy HIFU system was able to effectively treat PC3 tumors in a rabbit model using USMB therapy in combination with radiation exposures. Future work could find the use of reduced power and pressure levels in a commercial MR-guided therapy system to mechanically stimulate microbubbles and damage endothelial cells without requiring high thermal doses to elicit an antitumor response.

Validation of a drift-corrected 3D MR temperature imaging sequence for breast MR-guided focused ultrasound treatments

Real-time temperature monitoring is critical to the success of thermally ablative therapies. This work validates a 3D thermometry sequence with k-space field drift correction designed for use in magnetic resonance-guided focused ultrasound treatments for breast cancer. Fiberoptic probes were embedded in tissue-mimicking phantoms, and temperature change measurements from the probes were compared with the magnetic resonance temperature imaging measurements following heating with focused ultrasound. Precision and accuracy of measurements were also evaluated in free-breathing healthy volunteers (N = 3) under a non-heating condition. MR temperature measurements agreed closely with those of fiberoptic probes, with a 95% confidence interval of measurement difference from -2.0 °C to 1.4 °C. Field drift-corrected measurements in vivo had a precision of 1.1 ± 0.7 °C and were accurate within 1.3 ± 0.9 °C across the three volunteers. The field drift correction method improved precision and accuracy by an average of 46 and 42%, respectively, when compared to the uncorrected data. This temperature imaging sequence can provide accurate measurements of temperature change in aqueous tissues in the breast and support the use of this sequence in clinical investigations of focused ultrasound treatments for breast cancer.

Monitoring MR-guided high intensity focused ultrasound therapy using transient supersonic shear wave MR-elastography

Objective.The aim of the paper is to propose an all-in-one method based on magnetic resonance-supersonic shear wave imaging (MR-SSI) and proton resonance frequency shift (PRFS) to monitor high intensity focused ultrasound (HIFU) thermal ablations.Approach.Mechanical properties have been shown to be related to tissue damage induced by thermal ablations. Monitoring elasticity in addition to temperature changes may help in ensuring the efficacy and the accuracy of HIFU therapies. For this purpose, an MR-SSI method has been developed where the ultrasonic transducer is used for both mechanical wave generation and thermal ablation. Transient quasi-planar shear waves are generated using the acoustic radiation force, and their propagation is monitored in motion-sensitized phase MR images. Using a single-shot gradient-echo echo-planar-imaging sequence, MR images can be acquired at a sufficiently high temporal resolution to provide an update of PRFS thermometry and MR-SSI elastography maps in real time.Main results.The proposed method was first validated on a calibrated elasticity phantom, in which both the possibility to detect inclusions with different stiffness and repeatability were demonstrated. The standard deviation between the 8 performed measurements was 2% on the background of the phantom and 11%, at most, on the inclusions. A second experiment consisted in performing a HIFU heating in a gelatin phantom. The temperature increase was estimated to be 9 °C and the shear modulus was found to decrease from 2.9 to 1.8 kPa, reflecting the gel softening around the HIFU focus, whereas it remained steady in non-heated areas.Significance.The proposed MR-SSI technique allows monitoring HIFU ablations using thermometry and elastography simultaneously, without the need for an additional external mechanical exciter such as those used in MR elastography.

Multi-echo MR thermometry in the upper leg at 7 T using near-harmonic 2D reconstruction for initialization

PURPOSE

The aim of this work is the development of a thermometry method to measure temperature increases in vivo, with a precision and accuracy sufficient for validation against thermal simulations. Such an MR thermometry model would be a valuable tool to get an indication on one of the major safety concerns in MR imaging: the tissue heating occurring due to radiofrequency (RF) exposure. To prevent excessive temperature rise, RF power deposition, expressed as specific absorption rate, cannot exceed predefined thresholds. Using these thresholds, MRI has demonstrated an extensive history of safe usage. Nevertheless, MR thermometry would be a valuable tool to address some of the unmet needs in the area of RF safety assessment, such as validation of specific absorption rate and thermal simulations, investigation of local peak temperatures during scanning, or temperature-based safety guidelines.

METHODS

The harmonic initialized model-based multi-echo approach is proposed. The method combines a previously published model-based multi-echo water/fat separated approach with an also previously published near-harmonic 2D reconstruction method. The method is tested on the human thigh with a multi-transmit array at 7 T, in three volunteers, and for several RF shims.

RESULTS

Precision and accuracy are improved considerably compared to a previous fat-referenced method (precision: 0.09 vs. 0.19°C). Comparison of measured temperature rise distributions to subject-specific simulated counterparts show good relative agreement for multiple RF shim settings.

CONCLUSION

The high precision shows promising potential for validation purposes and other RF safety applications.

Fast, interleaved, Look-Locker-based T1 mapping with a variable averaging approach: Towards temperature mapping at low magnetic field

Proton resonance frequency shift (PRFS) is currently the gold standard method for magnetic resonance thermometry. However, the linearity between the temperature-dependent phase accumulation and the static magnetic field B₀ confines its use to rather high-field scanners. Applications such as thermal therapies could naturally benefit from lower field MRI settings through leveraging increased accessibility, a lower physical and economical footprint, and further consideration of the technical challenges associated with the integration of heating systems into conventional clinical scanners. T 1 -based thermometry has been proposed as an alternative to the gold standard; however, because of longer acquisition times, it has found clinical use solely with adipose tissue where PRFS fails. At low field, the enhanced T 1 dispersion, combined with reduced relaxation times, make T 1 mapping an appealing candidate. Here, an interleaved Look-Locker-based T 1 mapping sequence was proposed for temperature quantification at 0.1 T. A variable averaging scheme was introduced, to maximize the signal-to-noise ratio throughout T 1 recovery. In calibrated samples, an average T 1 accuracy of 85% ± 4% was achieved in 10 min, compared with the 77% ± 7% obtained using a standard averaging scheme. Temperature maps between 29.0 and 41.7°C were eventually reconstructed, with a precision of 3.0 ± 1.1°C and an accuracy of 1.5 ± 1.0°C. Accounting for longer thermal treatments and less strict temperature constraints, applications such as MR-guided mild hyperthermia treatments at low field could be envisioned.

Motion-robust, multi-slice, real-time MR thermometry for MR-guided thermal therapy in abdominal organs

PURPOSE

To develop an effective and practical reconstruction pipeline to achieve motion-robust, multi-slice, real-time MR thermometry for monitoring thermal therapy in abdominal organs.

METHODS

The application includes a fast spiral magnetic resonance imaging (MRI) pulse sequence and a real-time reconstruction pipeline based on multi-baseline proton resonance frequency shift (PRFS) method with visualization of temperature imaging. The pipeline supports multi-slice acquisition with minimal reconstruction lag. Simulations with a virtual motion phantom were performed to investigate the influence of the number of baselines and respiratory rate on the accuracy of temperature measurement. Phantom experiments with ultrasound heating were performed using a custom-made motion phantom to evaluate the performance of the pipeline. Lastly, experiments in healthy volunteers (N = 2) without heating were performed to evaluate the accuracy and stability of MR thermometry in abdominal organs (liver and kidney).

RESULTS

The multi-baseline approach with greater than 25 baselines resulted in minimal temperature errors in the simulation. Phantom experiments demonstrated a 713 ms update time for 3-slice acquisitions. Temperature maps with 30 baselines showed clear temperature distributions caused by ultrasound heating in the respiratory phantom. Finally, the pipeline was evaluated with physiologic motions in healthy volunteers without heating, which demonstrated the accuracy (root mean square error [RMSE]) of 1.23  ±   0.18 °C (liver) and 1.21  ±   0.17 °C (kidney) and precision of 1.13  ±   0.11 °C (liver) and 1.16  ±   0.15 °C (kidney) using 32 baselines.

CONCLUSIONS

The proposed real-time acquisition and reconstruction pipeline allows motion-robust, multi-slice, real-time temperature monitoring within the abdomen during free breathing.