Intravoxel incoherent motion (IVIM)-derived perfusion fraction mapping for the visual evaluation of MR-guided high intensity focused ultrasound (MR-HIFU) ablation of uterine fibroids

BACKGROUND

A method for periprocedural contrast agent-free visualization of uterine fibroid perfusion could potentially shorten magnetic resonance-guided high intensity focused ultrasound (MR-HIFU) treatment times and improve outcomes. Our goal was to test feasibility of perfusion fraction mapping by intravoxel incoherent motion (IVIM) modeling using diffusion-weighted MRI as method for visual evaluation of MR-HIFU treatment progression.

METHODS

Conventional and T2-corrected IVIM-derived perfusion fraction maps were retrospectively calculated by applying two fitting methods to diffusion-weighted MRI data (b = 0, 50, 100, 200, 400, 600 and 800 s/mm2 at 1.5 T) from forty-four premenopausal women who underwent MR-HIFU ablation treatment of uterine fibroids. Contrast in perfusion fraction maps between areas with low perfusion fraction and surrounding tissue in the target uterine fibroid immediately following MR-HIFU treatment was evaluated. Additionally, the Dice similarity coefficient (DSC) was calculated between delineated areas with low IVIM-derived perfusion fraction and hypoperfusion based on CE-T1w.

RESULTS

Average perfusion fraction ranged between 0.068 and 0.083 in areas with low perfusion fraction based on visual assessment, and between 0.256 and 0.335 in surrounding tissues (all p < 0.001). DSCs ranged from 0.714 to 0.734 between areas with low perfusion fraction and the CE-T1w derived non-perfused areas, with excellent intraobserver reliability of the delineated areas (ICC 0.97).

CONCLUSION

The MR-HIFU treatment effect in uterine fibroids can be visualized using IVIM perfusion fraction mapping, in moderate concordance with contrast enhanced MRI. IVIM perfusion fraction mapping has therefore the potential to serve as a contrast agent-free imaging method to visualize the MR-HIFU treatment progression in uterine fibroids.

Ultrasound and Microbubbles Mediated Bleomycin Delivery in Feline Oral Squamous Cell Carcinoma—An In Vivo Veterinary Study

To investigate the feasibility and tolerability of ultrasound and microbubbles (USMB)-enhanced chemotherapy delivery for head and neck cancer, we performed a veterinary trial in feline companion animals with oral squamous cell carcinomas. Six cats were treated with a combination of bleomycin and USMB therapy three times, using the Pulse Wave Doppler mode on a clinical ultrasound system and EMA/FDA approved microbubbles. They were evaluated for adverse events, quality of life, tumour response and survival. Furthermore, tumour perfusion was monitored before and after USMB therapy using contrast-enhanced ultrasound (CEUS). USMB treatments were feasible and well tolerated. Among 5 cats treated with optimized US settings, 3 had stable disease at first, but showed disease progression 5 or 11 weeks after first treatment. One cat had progressive disease one week after the first treatment session, maintaining a stable disease thereafter. Eventually, all cats except one showed progressive disease, but each survived longer than the median overall survival time of 44 days reported in literature. CEUS performed immediately before and after USMB therapy suggested an increase in tumour perfusion based on an increase in median area under the curve (AUC) in 6 out of 12 evaluated treatment sessions. In this small hypothesis-generating study, USMB plus chemotherapy was feasible and well-tolerated in a feline companion animal model and showed potential for enhancing tumour perfusion in order to increase drug delivery. This could be a forward step toward clinical translation of USMB therapy to human patients with a clinical need for locally enhanced treatment.

Synthetic CT for the planning of MR-HIFU treatment of bone metastases in pelvic and femoral bones: a feasibility study

Objectives

Visualization of the bone distribution is an important prerequisite for MRI-guided high-intensity focused ultrasound (MRI-HIFU) treatment planning of bone metastases. In this context, we evaluated MRI-based synthetic CT (sCT) imaging for the visualization of cortical bone.

Methods

MR and CT images of nine patients with pelvic and femoral metastases were retrospectively analyzed in this study. The metastatic lesions were osteolytic, osteoblastic or mixed. sCT were generated from pre-treatment or treatment MR images using a UNet-like neural network. sCT was qualitatively and quantitatively compared to CT in the bone (pelvis or femur) containing the metastasis and in a region of interest placed on the metastasis itself, through mean absolute difference (MAD), mean difference (MD), Dice similarity coefficient (DSC), and root mean square surface distance (RMSD).

Results

The dataset consisted of 3 osteolytic, 4 osteoblastic and 2 mixed metastases. For most patients, the general morphology of the bone was well represented in the sCT images and osteolytic, osteoblastic and mixed lesions could be discriminated. Despite an average timespan between MR and CT acquisitions of 61 days, in bone, the average (± standard deviation) MAD was 116 ± 26 HU, MD − 14 ± 66 HU, DSC 0.85 ± 0.05, and RMSD 2.05 ± 0.48 mm and, in the lesion, MAD was 132 ± 62 HU, MD − 31 ± 106 HU, DSC 0.75 ± 0.2, and RMSD 2.73 ± 2.28 mm.

Conclusions

Synthetic CT images adequately depicted the cancellous and cortical bone distribution in the different lesion types, which shows its potential for MRI-HIFU treatment planning.

Key Points

• Synthetic computed tomography was able to depict bone distribution in metastatic lesions.

• Synthetic computed tomography images intrinsically aligned with treatment MR images may have the potential to facilitate MR-HIFU treatment planning of bone metastases, by combining visualization of soft tissues and cancellous and cortical bone.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00330-022-08568-y.

Ultrasound-Mediated Drug Delivery With a Clinical Ultrasound System: In Vitro Evaluation

Chemotherapy efficacy is often reduced by insufficient drug uptake in tumor cells. The combination of ultrasound and microbubbles (USMB) has been shown to improve drug delivery and to enhance the efficacy of several drugs in vitro and in vivo, through effects collectively known as sonopermeation. However, clinical translation of USMB therapy is hampered by the large variety of (non-clinical) US set-ups and US parameters that are used in these studies, which are not easily translated to clinical practice. In order to facilitate clinical translation, the aim of this study was to prove that USMB therapy using a clinical ultrasound system (Philips iU22) in combination with clinically approved microbubbles (SonoVue) leads to efficient in vitro sonopermeation. To this end, we measured the efficacy of USMB therapy for different US probes (S5-1, C5-1 and C9-4) and US parameters in FaDu cells. The US probe with the lowest central frequency (i.e. 1.6 MHz for S5-1) showed the highest USMB-induced intracellular uptake of the fluorescent dye SYTOX™ Green (SG). These SG uptake levels were comparable to or even higher than those obtained with a custom-built US system with optimized US parameters. Moreover, USMB therapy with both the clinical and the custom-built US system increased the cytotoxicity of the hydrophilic drug bleomycin. Our results demonstrate that a clinical US system can be used to perform USMB therapy as efficiently as a single-element transducer set-up with optimized US parameters. Therefore, future trials could be based on these clinical US systems, including validated US parameters, in order to accelerate successful translation of USMB therapy.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Rapid 2D variable flip angle method for accurate and precise T1 measurements over a wide range of T1 values

PURPOSE

To perform dynamic T₁ mapping using a 2D variable flip angle (VFA) method, a correction for the slice profile effect is needed. In this work we investigated the impact of flip angle selection and excitation RF pulse profile on the performance of slice profile correction when applied to T₁ mapping over a range of T₁ values.

METHODS

A correction of the slice profile effect is proposed, based on Bloch simulation of steady-state signals. With this correction, Monte Carlo simulations were performed to assess the accuracy and precision of 2D VFA T₁ mapping in the presence of noise, for RF pulses with time-bandwidth products of 2, 3 and 10 and with flip angle pairs in the range [1°-90°]. To evaluate its performance over a wide range of T₁ , maximum errors were calculated for six T₁ values between 50 ms and 1250 ms. The method was demonstrated using in vitro and in vivo experiments.

RESULTS

Without corrections, 2D VFA severely underestimates T₁ . Slice profile errors were effectively reduced with the correction based on simulations, both in vitro and in vivo. The precision and accuracy of the method depend on the nominal T₁ values, the FA pair, and the RF pulse shape. FA pairs leading to <5% errors in T₁ can be identified for the common RF shapes, for T₁ values between 50 ms and 1250 ms.

CONCLUSIONS

2D VFA T₁ mapping with Bloch-simulation-based correction can deliver T₁ estimates that are accurate and precise to within 5% over a wide T₁ range.

Combining radiotherapy and focused ultrasound for pain palliation of cancer induced bone pain; a stage I/IIa study according to the IDEAL framework

•

Combined treatment of EBRT and MR-HIFU is feasible and well tolerated by patients.

•

Clinical outcomes of combined treatment of EBRT and MR-HIFU are promising.

•

Superiority of combined treatment over standard EBRT needs to be evaluated in a comparative study.

Background

Cancer induced bone pain (CIBP) strongly interferes with patient’s quality of life. Currently, the standard of care includes external beam radiotherapy (EBRT), resulting in pain relief in approximately 60% of patients. Magnetic Resonance guided High Intensity Focused Ultrasound (MR-HIFU) is a promising treatment modality for CIBP.

Methods

A single arm, R-IDEAL stage I/IIa study was conducted. Patients presenting at the department of radiation oncology with symptomatic bone metastases in the appendicular skeleton, as well as in the sacrum and sternum were eligible for inclusion. All participants underwent EBRT, followed by MR-HIFU within 4 days. Safety and feasibility were assessed, and pain scores were monitored for 4 weeks after completing the combined treatment.

Results

Six patients were enrolled. Median age was 67 years, median lesion diameter was 56,5 mm. In all patients it was logistically possible to plan and perform the MR-HIFU treatment within 4 days after EBRT. All patients tolerated the combined procedure well. Pain response was reported by 5 out of 6 patients at 7 days after completion of the combined treatment, and stabilized on 60% at 4 weeks follow up. No treatment related serious adverse events occurred.

Conclusion

This is the first study to combine EBRT with MR-HIFU. Our results show that combined EBRT and MR-HIFU in first-line treatment of CIBP is safe and feasible, and is well tolerated by patients. Superiority over standard EBRT, in terms of (time to) pain relief and quality of life need to be evaluated in comparative (randomized) study.

Deep correction of breathing-related artifacts in real-time MR-thermometry

Real-time MR-imaging has been clinically adapted for monitoring thermal therapies since it can provide on-the-fly temperature maps simultaneously with anatomical information. However, proton resonance frequency based thermometry of moving targets remains challenging since temperature artifacts are induced by the respiratory as well as physiological motion. If left uncorrected, these artifacts lead to severe errors in temperature estimates and impair therapy guidance. In this study, we evaluated deep learning for on-line correction of motion related errors in abdominal MR-thermometry. For this, a convolutional neural network (CNN) was designed to learn the apparent temperature perturbation from images acquired during a preparative learning stage prior to hyperthermia. The input of the designed CNN is the most recent magnitude image and no surrogate of motion is needed. During the subsequent hyperthermia procedure, the recent magnitude image is used as an input for the CNN-model in order to generate an on-line correction for the current temperature map. The method's artifact suppression performance was evaluated on 12 free breathing volunteers and was found robust and artifact-free in all examined cases. Furthermore, thermometric precision and accuracy was assessed for in vivo ablation using high intensity focused ultrasound. All calculations involved at the different stages of the proposed workflow were designed to be compatible with the clinical time constraints of a therapeutic procedure.

Phase I feasibility study of Magnetic Resonance guided High Intensity Focused Ultrasound-induced hyperthermia, Lyso-Thermosensitive Liposomal Doxorubicin and cyclophosphamide in de novo stage IV breast cancer patients: study protocol of the i-GO study

INTRODUCTION

In breast cancer, local tumour control is thought to be optimised by administering higher local levels of cytotoxic chemotherapy, in particular doxorubicin. However, systemic administration of higher dosages of doxorubicin is hampered by its toxic side effects. In this study, we aim to increase doxorubicin deposition in the primary breast tumour without changing systemic doxorubicin concentration and thus without interfering with systemic efficacy and toxicity. This is to be achieved by combining Lyso-Thermosensitive Liposomal Doxorubicin (LTLD, ThermoDox, Celsion Corporation, Lawrenceville, NJ, USA) with mild local hyperthermia, induced by Magnetic Resonance guided High Intensity Focused Ultrasound (MR-HIFU). When heated above 39.5°C, LTLD releases a high concentration of doxorubicin intravascularly within seconds. In the absence of hyperthermia, LTLD leads to a similar biodistribution and antitumour efficacy compared with conventional doxorubicin.

METHODS AND ANALYSIS

This is a single-arm phase I study in 12 chemotherapy-naïve patients with de novo stage IV HER2-negative breast cancer. Previous endocrine treatment is allowed. Study treatment consists of up to six cycles of LTLD at 21-day intervals, administered during MR-HIFU-induced hyperthermia to the primary tumour. We will aim for 60 min of hyperthermia at 40°C-42°C using a dedicated MR-HIFU breast system (Profound Medical, Mississauga, Canada). Afterwards, intravenous cyclophosphamide will be administered. Primary endpoints are safety, tolerability and feasibility. The secondary endpoint is efficacy, assessed by radiological response.This approach could lead to optimal loco-regional control with less extensive or even no surgery, in de novo stage IV patients and in stage II/III patients allocated to receive neoadjuvant chemotherapy.

ETHICS AND DISSEMINATION

This study has obtained ethical approval by the Medical Research Ethics Committee Utrecht (Protocol NL67422.041.18, METC number 18-702). Informed consent will be obtained from all patients before study participation. Results will be published in an academic peer-reviewed journal.

TRIAL REGISTRATION NUMBERS

NCT03749850, EudraCT 2015-005582-23.

A Doxorubicin-Glucuronide Prodrug Released from Nanogels Activated by High-Intensity Focused Ultrasound Liberated β-Glucuronidase

The poor pharmacokinetics and selectivity of low-molecular-weight anticancer drugs contribute to the relatively low effectiveness of chemotherapy treatments. To improve the pharmacokinetics and selectivity of these treatments, the combination of a doxorubicin-glucuronide prodrug (DOX-propGA3) nanogel formulation and the liberation of endogenous β-glucuronidase from cells exposed to high-intensity focused ultrasound (HIFU) were investigated in vitro. First, a DOX-propGA3-polymer was synthesized. Subsequently, DOX-propGA3-nanogels were formed from this polymer dissolved in water using inverse mini-emulsion photopolymerization. In the presence of bovine β-glucuronidase, the DOX-propGA3 in the nanogels was quantitatively converted into the chemotherapeutic drug doxorubicin. Exposure of cells to HIFU efficiently induced liberation of endogenous β-glucuronidase, which in turn converted the prodrug released from the DOX-propGA3-nanogels into doxorubicin. β-glucuronidase liberated from cells exposed to HIFU increased the cytotoxicity of DOX-propGA3-nanogels to a similar extend as bovine β-glucuronidase, whereas in the absence of either bovine β-glucuronidase or β-glucuronidase liberated from cells exposed to HIFU, the DOX-propGA3-nanogels hardly showed cytotoxicity. Overall, DOX-propGA3-nanogels systems might help to further improve the outcome of HIFU-related anticancer therapy.

Use of multiparametric MRI to characterize uterine fibroid tissue types

BACKGROUND

Although the biological characteristics of uterine fibroids (UF) have implications for therapy choice and effectiveness, there is limited MRI data about these characteristics. Currently, the Funaki classification and Scaled Signal Intensity (SSI) are used to predict treatment outcome but both screening-tools appear to be suboptimal. Therefore, multiparametric and quantitative MRI was studied to evaluate various biological characteristics of UF.

METHODS

87 patients with UF underwent an MRI-examination. Differences between UF tissues and myometrium were investigated using T2-mapping, Apparent Diffusion Coefficient (ADC) maps with different b-value combinations, contrast-enhanced T1-weighted and T2-weighted imaging. Additionally, the Funaki classification and SSI were calculated.

RESULTS

Significant differences between myometrium and UF tissue in T2-mapping (p = 0.001), long-TE ADC low b-values (p = 0.002), ADC all b-values (p < 0.001) and high b-values (p < 0.001) were found. Significant differences between Funaki type 3 versus type 1 and 2 were observed in SSI (p < 0.001) and T2-values (p < 0.001). Significant correlations were found between SSI and T2-mapping (p < 0.001; ρ_s = 0.82), ADC all b-values (p = 0.004; ρ_s = 0.31), ADC high b-values (p < 0.001; ρ_s = 0.44) and long-TE ADC low b-values (p = 0.004; ρ_s = 0.31).

CONCLUSIONS

Quantitative MR-data allowed us to distinguish UF tissue from myometrium and to discriminate different UF tissue types and may, therefore, be a useful tool to predict treatment outcome/determine optimal treatment modality.

Tumor Drug Distribution after Local Drug Delivery by Hyperthermia, In Vivo

Tumor drug distribution and concentration are important factors for effective tumor treatment. A promising method to enhance the distribution and the concentration of the drug in the tumor is to encapsulate the drug in a temperature sensitive liposome. The aim of this study was to investigate the tumor drug distribution after treatment with various injected doses of different liposomal formulations of doxorubicin, ThermoDox (temperature sensitive liposomes) and DOXIL (non-temperature sensitive liposomes), and free doxorubicin at macroscopic and microscopic levels. Only ThermoDox treatment was combined with hyperthermia. Experiments were performed in mice bearing a human fibrosarcoma. At low and intermediate doses, the largest growth delay was obtained with ThermoDox, and at the largest dose, the largest growth delay was obtained with DOXIL. On histology, tumor areas with increased doxorubicin concentration correlated with decreased cell proliferation, and substantial variations in doxorubicin heterogeneity were observed. ThermoDox treatment resulted in higher tissue drug levels than DOXIL and free doxorubicin for the same dose. A relation with the distance to the vasculature was shown, but vessel perfusion was not always sufficient to determine doxorubicin delivery. Our results indicate that tumor drug distribution is an important factor for effective tumor treatment and that its dependence on delivery formulation merits further systemic investigation.

Field drift correction of proton resonance frequency shift temperature mapping with multichannel fast alternating nonselective free induction decay readouts

PURPOSE

To demonstrate that proton resonance frequency shift MR thermometry (PRFS-MRT) acquisition with nonselective free induction decay (FID), combined with coil sensitivity profiles, allows spatially resolved B₀ drift-corrected thermometry.

METHODS

Phantom experiments were performed at 1.5T and 3T. Acquisition of PRFS-MRT and FID were performed during MR-guided high-intensity focused ultrasound heating. The phase of the FIDs was used to estimate the change in angular frequency δω_drift per coil element. Two correction methods were investigated: (1) using the average δω_drift over all coil elements (0th-order) and (2) using coil sensitivity profiles for spatially resolved correction. Optical probes were used for independent temperature verification. In-vivo feasibility of the methods was evaluated in the leg of 1 healthy volunteer at 1.5T.

RESULTS

In 30 minutes, B₀ drift led to an apparent temperature change of up to -18°C and -98°C at 1.5T and 3T, respectively. In the sonicated area, both corrections had a median error of 0.19°C at 1.5T and -0.54°C at 3T. At 1.5T, the measured median error with respect to the optical probe was -1.28°C with the 0th-order correction and improved to 0.43°C with the spatially resolved correction. In vivo, without correction the spatiotemporal median of the apparent temperature was at -4.3°C and interquartile range (IQR) of 9.31°C. The 0th-order correction had a median of 0.75°C and IQR of 0.96°C. The spatially resolved method had the lowest median at 0.33°C and IQR of 0.80°C.

CONCLUSION

FID phase information from individual receive coil elements allows spatially resolved B₀ drift correction in PRFS-based MRT.

Investigation of the influence of B0 drift on the performance of the PLANET method and an algorithm for drift correction

Purpose

The PLANET method was designed to simultaneously reconstruct maps of T₁ and T₂, the off‐resonance, the RF phase, and the banding free signal magnitude. The method requires a stationary B₀ field over the course of a phase‐cycled balanced SSFP acquisition. In this work we investigated the influence of B₀ drift on the performance of the PLANET method for single‐component and two‐component signal models, and we propose a strategy for drift correction.

Methods

The complex phase‐cycled balanced SSFP signal was modeled with and without frequency drift. The behavior of the signal influenced by drift was mathematically interpreted as a sum of drift‐dependent displacement of the data points along an ellipse and drift‐dependent rotation around the origin. The influence of drift on parameter estimates was investigated experimentally on a phantom and on the brain of healthy volunteers and was verified by numerical simulations. A drift correction algorithm was proposed and tested on a phantom and in vivo.

Results

Drift can be assumed to be linear over the typical duration of a PLANET acquisition. In a phantom (a single‐component signal model), drift induced errors of 4% and 8% in the estimated T₁ and T₂ values. In the brain, where multiple components are present, drift only had a minor effect. For both single‐component and two‐component signal models, drift‐induced errors were successfully corrected by applying the proposed drift correction algorithm.

Conclusion

We have demonstrated theoretically and experimentally the sensitivity of the PLANET method to B₀ drift and have proposed a drift correction method.

Hyperthermia-triggered release of hypoxic cell radiosensitizers from temperature-sensitive liposomes improves radiotherapy efficacy in vitro

Hypoxia is a characteristic feature of solid tumors and an important cause of resistance to radiotherapy. Hypoxic cell radiosensitizers have been shown to increase radiotherapy efficacy, but dose-limiting side effects prevent their widespread use in the clinic. We propose the encapsulation of hypoxic cell radiosensitizers in temperature-sensitive liposomes (TSL) to target the radiosensitizers specifically to tumors and to avoid unwanted accumulation in healthy tissues. The main objective of the present study is to develop and characterize TSL loaded with the radiosensitizer pimonidazole (PMZ) and to evaluate the in vitro efficacy of free PMZ and PMZ encapsulated in TSL in combination with hyperthermia and radiotherapy. PMZ was actively loaded into TSL at different drug/lipid ratios, and the physicochemical characteristics and the stability of the resulting TSL-PMZ were evaluated. PMZ release was determined at 37 °C and 42 °C in HEPES buffer saline and fetal bovine serum. The concentration-dependent radiosensitizing effect of PMZ was investigated by exposing FaDu cells to different PMZ concentrations under hypoxic conditions followed by exposure to ionizing irradiation. The efficacy of TSL-PMZ in combination with hyperthermia and radiotherapy was determined in vitro, assessing cell survival and DNA damage by means of the clonogenic assay and histone H2AX phosphorylation, respectively. All TSL-PMZ formulations showed high encapsulation efficiencies and were stable for 30 d upon storage at 4 °C and 20 °C. Fast PMZ release was observed at 42 °C, regardless of the drug/lipid ratio. Increasing the PMZ concentration significantly enhanced the effect of ionizing irradiation. Pre-heated TSL-PMZ in combination with radiotherapy caused a 14.3-fold increase in cell death as compared to radiotherapy treatment alone. In conclusion, our results indicate that TSL-PMZ in combination with hyperthermia can assist in improving the efficacy of radiotherapy under hypoxic conditions.

Assessment of 3D motion modeling performance for dose accumulation mapping on the MR-linac by simultaneous multislice MRI

Hybrid MR-linac systems enable intrafraction motion monitoring during radiation therapy. Since time-resolved 3D MRI is still challenging, various motion models have been developed that rely on time-resolved 2D imaging. Continuous validation of these models is important for accurate dose accumulation mapping. In this study we used 2D simultaneous multislice (SMS) imaging to improve the PCA-based motion modeling method developed previously (Stemkens et al 2016 Phys. Med. Biol. 61 5335-55). From the additional simultaneously acquired slices, several independent motion models could be generated, which allowed for an assessment of the sensitivity of the motion model to the location of the time-resolved 2D slices. Additionally, the best model could be chosen at every time-point, increasing the method's robustness. Imaging experiments were performed in six healthy volunteers using three simultaneous slices, which generated three independent models per volunteer. For each model the motion traces of the liver tip and both kidneys were estimated. We found that the location of the 2D slices influenced the model's error in five volunteers significantly with a p -value  <0.05, and that selecting the best model at every time-point can improve the method. This allows for more accurate and robust motion characterization in MR-guided radiotherapy.

Respiratory- and cardiac-triggered three-dimensional sheath inked rapid acquisition with refocused echoes imaging (SHINKEI) of the abdomen for magnetic resonance neurography of the celiac plexus

The visualisation of the celiac plexus using respiratory- and cardiac-triggered three-dimensional (3D) sheath inked rapid acquisition with refocused echoes imaging (SHINKEI) was evaluated. After ethical approval and written informed consent, eight volunteers (age 27 ± 5 years, mean ± standard deviation) were scanned at 1.5 and 3 T. Displacement of the celiac ganglia due to aortic pulsatility was studied on axial single-slice breath-hold balanced turbo field-echo cine sequences in five volunteers and found to be 3.0 ± 0.5 mm (left) and 3.1 ± 0.4 mm (right). Respiratory- and cardiac-triggered 3D SHINKEI images were compared to respiratory- and cardiac-triggered fat-suppressed 3D T2-weighted turbo spin-echo and respiratory-triggered 3D SHINKEI in all volunteers. Visibility of the celiac ganglia was rated by three radiologists as visible or non-visible. On 3D SHINKEI with double-triggering at 1.5 T, the left and right ganglia were seen by all observers in 7/8 and 8/8 volunteers, respectively. At 3 T, this was the case for 6/8 and 7/8 volunteers, respectively. The nerve-to-muscle signal ratio increased from 1.9 ± 0.5 on fat-suppressed 3D T2-weighted turbo spin-echo to 4.7 ± 0.8 with 3D SHINKEI. Anatomical validation was performed in a human cadaver. An expert in anatomy confirmed that the hyperintense structure visible on ex vivo 3D SHINKEI scans was the celiac plexus. In conclusion, double-triggering allowed visualisation of the celiac plexus using 3D SHINKEI at both 1.5 T and 3 T.

Enabling free-breathing background suppressed renal pCASL using fat imaging and retrospective motion correction

Purpose

For free‐breathing renal perfusion imaging using arterial spin labeling (ASL), retrospective image realignment has been found essential to reduce subtraction artifacts and, independently, background suppression has been demonstrated to reduce physiologic noise. However, negative results on ASL precision and accuracy have been reported for the combination of both. In this study, the effect of background suppression ‐level in combination with image registration on free‐breathing renal ASL signal quality, with registration either on ASL‐images themselves or guided by additionally acquired fat‐images, was investigated. The results from free‐breathing acquisitions were compared with the reference paced‐breathing motion compensation strategy.

Methods

Pseudocontinuous ASL (pCASL) data with additional fat‐images were acquired from 10 subjects at 1.5T with varying background suppression levels during free‐breathing and paced‐breathing. Images were registered using the ASL‐images themselves (ASLReg) or using their corresponding fat‐images (FatReg). Temporal signal‐to‐noise ratio (tSNR) served to evaluate precision and perfusion weighted signal (PWS) to assess accuracy.

Results

In combination with image registration, background suppression significantly improved tSNR by 50% (P < .05). For heavy suppression, ASLReg and FatReg showed similar performance in terms of tSNR and PWS. Background suppression with two inversion pulses induced a small, nonsignificant (P > .05) PWS reduction, but increased PWS accuracy. When applying heavy background suppression, free‐breathing acquisitions resulted in similar ASL‐quality to paced‐breathing acquisitions.

Conclusion

Background suppression was found beneficial for free‐breathing renal pCASL precision without compromising accuracy, despite motion challenges. In combination with ASLReg or FatReg, background suppression enabled clinically viable free‐breathing renal pCASL.

On the accuracy and precision of PLANET for multiparametric MRI using phase-cycled bSSFP imaging

Purpose

In this work we demonstrate how sequence parameter settings influence the accuracy and precision in T₁, T₂, and off‐resonance maps obtained with the PLANET method for a single‐component signal model. In addition, the performance of the method for the particular case of a two‐component relaxation model for white matter tissue was assessed.

Methods

Numerical simulations were performed to investigate the influence of sequence parameter settings on the accuracy and precision in the estimated parameters for a single‐component model, as well as for a two‐component white matter model. Phantom and in vivo experiments were performed for validation. In addition, the effects of Gibbs ringing were investigated.

Results

By making a proper choice for sequence parameter settings, accurate and precise parameter estimation can be achieved for a single‐component signal model over a wide range of relaxation times at realistic SNR levels. Due to the presence of a second myelin‐related signal component in white matter, an underestimation of approximately 30% in T₁ and T₂ was observed, predicted by simulations and confirmed by measurements. Gibbs ringing artifacts correction improved the precision and accuracy of the parameter estimates.

Conclusion

For a single‐component signal model there is a broad “sweet spot” of sequence parameter combinations for which a high accuracy and precision in the parameter estimates is achieved over a wide range of relaxation times. For a multicomponent signal model, the single‐component PLANET reconstruction results in systematic errors in the parameter estimates as expected.

Anatomically plausible models and quality assurance criteria for online mono- and multi-modal medical image registration

Medical imaging is currently employed in the diagnosis, planning, delivery and response monitoring of cancer treatments. Due to physiological motion and/or treatment response, the shape and location of the pathology and organs-at-risk may change over time. Establishing their location within the acquired images is therefore paramount for an accurate treatment delivery and monitoring. A feasible solution for tracking anatomical changes during an image-guided cancer treatment is provided by image registration algorithms. Such methods are, however, often built upon elements originating from the computer vision/graphics domain. Since the original design of such elements did not take into consideration the material properties of particular biological tissues, the anatomical plausibility of the estimated deformations may not be guaranteed. In the current work we adapt two existing variational registration algorithms, namely Horn-Schunck and EVolution, to online soft tissue tracking. This is achieved by enforcing an incompressibility constraint on the estimated deformations during the registration process. The existing and the modified registration methods were comparatively tested against several quality assurance criteria on abdominal in vivo MR and CT data. These criteria included: the Dice similarity coefficient (DSC), the Jaccard index, the target registration error (TRE) and three additional criteria evaluating the anatomical plausibility of the estimated deformations. Results demonstrated that both the original and the modified registration methods have similar registration capabilities in high-contrast areas, with DSC and Jaccard index values predominantly in the 0.8-0.9 range and an average TRE of 1.6-2.0 mm. In contrast-devoid regions of the liver and kidneys, however, the three additional quality assurance criteria have indicated a considerable improvement of the anatomical plausibility of the deformations estimated by the incompressibility-constrained methods. Moreover, the proposed registration models maintain the potential of the original methods for online image-based guidance of cancer treatments.

A framework for continuous target tracking during MR-guided high intensity focused ultrasound thermal ablations in the abdomen

Background

During lengthy magnetic resonance-guided high intensity focused ultrasound (MRg-HIFU) thermal ablations in abdominal organs, the therapeutic work-flow is frequently hampered by various types of physiological motion occurring at different time-scales. If left un-addressed this can lead to an incomplete therapy and/or to tissue damage of organs-at-risk. While previous studies focus on correction schemes for displacements occurring at a particular time-scale within the work-flow of an MRg-HIFU therapy, in the current work we propose a motion correction strategy encompassing the entire work-flow.

Methods

The proposed motion compensation framework consists of several linked components, each being adapted to motion occurring at a particular time-scale. While respiration was addressed through a fast correction scheme, long term organ drifts were compensated using a strategy operating on time-scales of several minutes. The framework relies on a periodic examination of the treated area via MR scans which are then registered to a reference scan acquired at the beginning of the therapy. The resulting displacements were used for both on-the-fly re-optimization of the interventional plan and to ensure the spatial fidelity between the different steps of the therapeutic work-flow. The approach was validated in three complementary studies: an experiment conducted on a phantom undergoing a known motion pattern, a study performed on the abdomen of 10 healthy volunteers and during 3 in-vivo MRg-HIFU ablations on porcine liver.

Results

Results have shown that, during lengthy MRg-HIFU thermal therapies, the human liver and kidney can manifest displacements that exceed acceptable therapeutic margins. Also, it was demonstrated that the proposed framework is capable of providing motion estimates with sub-voxel precision and accuracy. Finally, the 3 successful animal studies demonstrate the compatibility of the proposed approach with the work-flow of an MRg-HIFU intervention under clinical conditions.

Conclusions

In the current study we proposed an image-based motion compensation framework dedicated to MRg-HIFU thermal ablations in the abdomen, providing the possibility to re-optimize the therapy plan on-the-fly with the patient on the interventional table. Moreover, we have demonstrated that even under clinical conditions, the proposed approach is fully capable of continuously ensuring the spatial fidelity between the different phases of the therapeutic work-flow.

Dynamic Fluorescence Microscopy of Cellular Uptake of Intercalating Model Drugs by Ultrasound-Activated Microbubbles

PURPOSE

The combination of ultrasound and microbubbles can facilitate cellular uptake of (model) drugs via transient permeabilization of the cell membrane. By using fluorescent molecules, this process can be studied conveniently with confocal fluorescence microscopy. This study aimed to investigate the relation between cellular uptake and fluorescence intensity increase of intercalating model drugs.

PROCEDURES

SYTOX Green, an intercalating fluorescent dye that displays >500-fold fluorescence enhancement upon binding to nucleic acids, was used as a model drug for ultrasound-induced cellular uptake. SYTOX Green uptake was monitored in high spatiotemporal resolution to qualitatively assess the relation between uptake and fluorescence intensity in individual cells. In addition, the kinetics of fluorescence enhancement were studied as a function of experimental parameters, in particular, laser duty cycle (DC), SYTOX Green concentration and cell line.

RESULTS

Ultrasound-induced intracellular SYTOX Green uptake resulted in local fluorescence enhancement, spreading throughout the cell and ultimately accumulating in the nucleus during the 9-min acquisition. The temporal evolution of SYTOX Green fluorescence was substantially influenced by laser duty cycle: continuous laser (100 % DC) induced a 6.4-fold higher photobleaching compared to pulsed laser (3.3 % DC), thus overestimating the fluorescence kinetics. A positive correlation of fluorescence kinetics and SYTOX Green concentration was found, increasing from 0.6 × 10-3 to 2.2 × 10-3 s-1 for 1 and 20 μM, respectively. Finally, C6 cells displayed a 2.4-fold higher fluorescence rate constant than FaDu cells.

CONCLUSIONS

These data show that the temporal behavior of intracellular SYTOX Green fluorescence enhancement depends substantially on nuclear accumulation and not just on cellular uptake. In addition, it is strongly influenced by the experimental conditions, such as the laser duty cycle, SYTOX Green concentration, and cell line.

Pharmacological and physical vessel modulation strategies to improve EPR-mediated drug targeting to tumors

The performance of nanomedicine formulations depends on the Enhanced Permeability and Retention (EPR) effect. Prototypic nanomedicine-based drug delivery systems, such as liposomes, polymers and micelles, aim to exploit the EPR effect to accumulate at pathological sites, to thereby improve the balance between drug efficacy and toxicity. Thus far, however, tumor-targeted nanomedicines have not yet managed to achieve convincing therapeutic results, at least not in large cohorts of patients. This is likely mostly due to high inter- and intra-patient heterogeneity in EPR. Besides developing (imaging) biomarkers to monitor and predict EPR, another strategy to address this heterogeneity is the establishment of vessel modulation strategies to homogenize and improve EPR. Over the years, several pharmacological and physical co-treatments have been evaluated to improve EPR-mediated tumor targeting. These include pharmacological strategies, such as vessel permeabilization, normalization, disruption and promotion, as well as physical EPR enhancement via hyperthermia, radiotherapy, sonoporation and phototherapy. In the present manuscript, we summarize exemplary studies showing that pharmacological and physical vessel modulation strategies can be used to improve tumor-targeted drug delivery, and we discuss how these advanced combination regimens can be optimally employed to enhance the (pre-) clinical performance of tumor-targeted nanomedicines.

Early health technology assessment of magnetic resonance-guided high intensity focused ultrasound ablation for the treatment of early-stage breast cancer

BACKGROUND

Magnetic resonance-guided high intensity focused ultrasound (MR-HIFU) ablation is in development for minimally invasive treatment of breast cancer. Cost-effectiveness has not been assessed yet. An early health technology assessment was performed to estimate costs of MR-HIFU ablation, compared to breast conserving treatment (BCT).

METHODS

An MR-HIFU treatment model using the dedicated MR-HIFU breast system (Sonalleve, Philips Healthcare) was developed. Input parameters (treatment steps and duration) were based on the analysis of questionnaire data from an expert panel. MR-HIFU experts assessed face validity of the model. Data collected by questionnaires were compared to published data of an MR-HIFU breast feasibility study. Treatment costs for tumours of 1 to 3 cm were calculated.

RESULTS

The model structure was considered of acceptable face validity by consulted experts, and questionnaire data and published data were comparable. Costs of MR-HIFU ablation were higher than BCT costs. MR-HIFU best-case scenario costs exceeded BCT costs with approximately €1000. Cooling times and breathing correction contributed most to treatment costs.

CONCLUSIONS

MR-HIFU ablation is currently not a cost-effective alternative for BCT. MR-HIFU experience is limited, increasing uncertainty of estimations. The potential for cost-effectiveness increases if future research reduces treatment durations and might substantiate equal or improved results.

Thermal ablation of a confluent lesion in the porcine kidney with a clinically available MR-HIFU system

The incidence of small renal masses (SRMs) sized  <4 cm has increased over the decades (as co-findings/or due to introduction of cross sectional imaging). Currently, partial nephrectomy (PN) or watchful waiting is advised in these patients. Ultimately, 80-90% of these SRMs require surgical treatment and PN is associated with a 15% complication rate. In this aging population, with possible comorbidities and poor health condition, both PN and watchful waiting are non-ideal treatment options. This resulted in an increased need for early, non-invasive treatment strategies such as MR-guided high intensity focused ultrasound (MR-HIFU). (i) To investigate the feasibility of creating a confluent lesion in the kidney using respiratory-gated MR-HIFU under clinical conditions in a pre-clinical study and (ii) to evaluate the reproducibility of the MR-HIFU ablation strategy. Healthy pigs (n  =  10) under general anesthesia were positioned on a clinical MR-HIFU system with integrated cooling. A honeycomb pattern of seven overlapping ablation cells (4  ×  4  ×  10 mm³, 450 W, <30 s) was ablated successively in the cortex of the porcine kidney. Both MR thermometry and acoustic energy delivery were respiratory gated using a pencil beam navigator on the contralateral kidney. The non-perfused volume (NPV) was visualized after the last sonication by contrast-enhanced (CE) T ₁-weighted MR (T ₁ w) imaging. Cell viability staining was performed to visualize the extent of necrosis.

RESULTS

a median NPV of 0.62 ml was observed on CE-T ₁ w images (IQR 0.58-1.57 ml, range 0.33-2.75 ml). Cell viability staining showed a median damaged volume of 0.59 ml (IQR 0.24-1.35 ml, range 0-4.1 ml). Overlooking of the false rib, shivering of the pig, and too large depth combined with a large heat-sink effect resulted in insufficient heating in 4 cases. The NPV and necrosed volume were confluent in all cases in which an ablated volume could be observed. Our results demonstrated the feasibility of creating a confluent volume of ablated kidney cortical tissue in vivo with MR-HIFU on a clinically available system using respiratory gating and near-field cooling and showed its reproducibility.

Fluid filling of the digestive tract for improved proton resonance frequency shift-based MR thermometry in the pancreas

PURPOSE

To demonstrate that fluid filling of the digestive tract improves the performance of respiratory motion-compensated proton resonance frequency shift (PRFS)-based magnetic resonance (MR) thermometry in the pancreas.

MATERIALS AND METHODS

In seven volunteers (without heating), we evaluated PRFS thermometry in the pancreas with and without filling of the surrounding digestive tract. All data acquisition was performed at 1.5T, then all datasets were analyzed and compared with three different PRFS respiratory motion-compensated thermometry methods: gating, multibaseline, and referenceless. The temperature precision of the different methods was evaluated by assessing temperature standard deviation over time, while a simulation experiment was used to study the accuracy of the methods.

RESULTS

Without fluid intake, errors in temperature precision in the pancreas up to 10°C were observed for all evaluated methods. After liquid intake, temperature precision improved to median values between 1.8 and 2.9°C. The simulations showed that gating had the lowest accuracy, with errors up to 7°C. Multibaseline and referenceless thermometry performed better, with a median error in the pancreas between -3 and +3°C after fluid intake, for all volunteers.

CONCLUSION

Preparation of the digestive tract near the pancreas by filling it with fluid improved MR thermometry precision and accuracy for all common respiratory motion-compensated methods evaluated. These improvements are attributed to reducing field inhomogeneity in the pancreas.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2018;47:692-701.

PLANET: An ellipse fitting approach for simultaneous T1 and T2 mapping using phase-cycled balanced steady-state free precession

Purpose

To demonstrate the feasibility of a novel, ellipse fitting approach, named PLANET, for simultaneous estimation of relaxation times T₁ and T₂ from a single 3D phase‐cycled balanced steady‐state free precession (bSSFP) sequence.

Methods

A method is presented in which the elliptical signal model is used to describe the phase‐cycled bSSFP steady‐state signal. The fitting of the model to the acquired data is reformulated into a linear convex problem, which is solved directly by a linear least squares method, specific to ellipses. Subsequently, the relaxation times T₁ and T₂, the banding free magnitude, and the off‐resonance are calculated from the fitting results.

Results

Maps of T₁ and T₂, as well as an off‐resonance and a banding free magnitude can be simultaneously, quickly, and robustly estimated from a single 3D phase‐cycled bSSFP sequence. The feasibility of the method was demonstrated in a phantom and in the brain of healthy volunteers on a clinical MR scanner. The results were in good agreement for the phantom, but a systematic underestimation of T₁ was observed in the brain.

Conclusion

The presented method allows for accurate mapping of relaxation times and off‐resonance, and for the reconstruction of banding free magnitude images at realistic signal‐to‐noise ratios. Magn Reson Med 79:711–722, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution‐NonCommercial‐NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non‐commercial and no modifications or adaptations are made.

Short and long time MR signal behavior of randomly distributed water and fat-numerical simulations

The MR time-signal behavior of water has been reported to be different on short and long time scales for systems of randomly distributed perturbers in water in the static dephasing regime. Up to now, the signal of the perturbers in such systems has not been taken into consideration. Water-fat emulsions are macroscopically homogeneous systems and can be considered as microscopically randomly distributed perturbing fat spheres embedded in water. In such water-fat systems, the signal of the perturber, fat, cannot be ignored. Since water and fat are within the same system, the fat signal behavior may show similarities with water, with differences in short and long time scales. This could complicate fat-referenced MR thermometry (MRT) methods such as multi-gradient echo-based (MGE) MRT. Simulations were performed using a numerical phantom comprising spherical fat objects embedded in a spherical water medium. To characterize the fat signal, the theoretical signal description of water was fitted to the simulated fat signal. The simulated signals were sampled as an MGE signal and MGE MRT was used to calculate temperatures. The sampling was done with and without delay, to investigate the effect on the temperature error of the time ranges in which the signal was sampled. It was confirmed that the fat signal behavior was similar to that of water and consisted of two regimes. The separation between the short and long time scales was approximately at 55 ms for fat, as compared with 8.9 ms for water. Without delayed signal sampling, the MGE MRT temperature error was about 2.5°C. With delayed sampling such that both the water and the fat signals were either in the short or in the long time scale the error was reduced to 0.2°C.

Spatiotemporal control of gene expression in bone-marrow derived cells of the tumor microenvironment induced by MRI guided focused ultrasound

The tumor microenvironment is an interesting target for anticancer therapies but modifying this compartment is challenging. Here, we demonstrate the feasibility of a gene therapy strategy that combined targeting to bone marrow-derived tumor microenvironment using genetically modified bone-marrow derived cells and control of transgene expression by local hyperthermia through a thermo-inducible promoter. Chimera were obtained by engraftment of bone marrow from transgenic mice expressing reporter genes under transcriptional control of heat shock promoter and inoculated sub-cutaneously with tumors cells. Heat shocks were applied at the tumor site using a water bath or magnetic resonance guided high intensity focused ultrasound device. Reporter gene expression was followed by bioluminescence and fluorescence imaging and immunohistochemistry. Bone marrow-derived cells expressing reporter genes were identified to be mainly tumor-associated macrophages. We thus provide the proof of concept for a gene therapy strategy that allows for spatiotemporal control of transgenes expression by macrophages targeted to the tumor microenvironment.

Procedural sedation and analgesia for respiratory-gated MR-HIFU in the liver: a feasibility study

BACKGROUND

Previous studies demonstrated both pre-clinically and clinically the feasibility of magnetic resonance-guided high-intensity focused ultrasound (MR-HIFU) ablations in the liver. To overcome the associated problem of respiratory motion of the ablation area, general anesthesia (GA) and mechanical ventilation was used in conjunction with either respiratory-gated energy delivery or energy delivery during induced apnea. However, clinical procedures requiring GA are generally associated with increased mortality, morbidity, and complication rate compared to procedural sedation and analgesia (PSA). Furthermore, PSA is associated with faster recovery and an increased eligibility for non- and mini-invasive interventions.

METHODS

In this study, we investigate both in an animal model and on a small patient group the kinetics of the diaphragm during free-breathing, when a tailored remifentanil/propofol-based PSA protocol inducing partial respiratory depression is used. Subsequently, we demonstrate in an animal study the compatibility of the resulting respiratory pattern of the PSA protocol with a gated HIFU ablation in the liver by direct comparison with gated ablations conducted under GA. Wilcoxon signed-rank tests were performed for statistical analysis of non-perfused and necrosed tissue volumes. Duty cycles (ratio or percentage of the breathing cycle with the diaphragm in its resting position, such that acoustic energy delivery with MR-HIFU was allowed) were statistically compared for both GA and PSA using student's t tests.

RESULTS

In both animal and human experiments, the breathing frequency was decreased below 9/min, while maintaining stable vital functions. Furthermore an end-exhalation resting phase was induced by this PSA protocol during which the diaphragm is virtually immobile. Median non-perfused volumes, non-viable volumes based on NADH staining, and duty cycles were larger under PSA than under GA or equal.

CONCLUSIONS

We conclude that MR-HIFU ablations of the liver under PSA are feasible and potentially increase the non-invasive nature of this type of intervention.

DCE-MRI and IVIM-MRI of rabbit Vx2 tumors treated with MR-HIFU-induced mild hyperthermia

Background

The purpose of this study is to investigate whether changes could be detected in dynamic contrast-enhanced (DCE) and intra-voxel incoherent motion (IVIM) MR parameters upon MR-guided high-intensity focused ultrasound (MR-HIFU)-induced hyperthermia in a rabbit Vx2 tumor model.

Methods

Five Vx2 tumor-bearing New Zealand white rabbits were treated with hyperthermia using a clinical MR-HIFU system. Data were acquired before and after hyperthermia. For the DCE analysis, the extended Tofts model was used. For the IVIM analysis, a Bayesian approach was used. Maps were reconstructed of the DCE parameters (K^trans, k_ep, and v_p) and IVIM parameters (D_t, f_p, and D_p). Individual parameter histograms and two-dimensional cross-correlation histograms were constructed to analyze changes in the parameters after hyperthermia. Changes in median values were tested for statistical significance with the Mann-Whitney U test.

Results

The MR temperature measurements confirmed that mild hyperthermia (40 to 42 °C) was successfully achieved in all rabbits. One rabbit died during treatment and was excluded from the analysis. In the remaining four rabbits, an increase in D_t was observed. In three rabbits, an increase in K^trans was observed, while in the other rabbits, all three DCE parameter values decreased. Mixed changes were seen for v_p and f_p.

Conclusions

Changes in DCE and IVIM parameters were detected after hyperthermia and were variable between the rabbits. DCE- and IVIM-MRI may be promising tools to assess tumor responses to hyperthermia. Further research in a larger number of subjects is necessary in order to assess their value for treatment response monitoring.

Sonochemotherapy: from bench to bedside

The combination of microbubbles and ultrasound has emerged as a promising method for local drug delivery. Microbubbles can be locally activated by a targeted ultrasound beam, which can result in several bio-effects. For drug delivery, microbubble-assisted ultrasound is used to increase vascular- and plasma membrane permeability for facilitating drug extravasation and the cellular uptake of drugs in the treated region, respectively. In the case of drug-loaded microbubbles, these two mechanisms can be combined with local release of the drug following destruction of the microbubble. The use of microbubble-assisted ultrasound to deliver chemotherapeutic agents is also referred to as sonochemotherapy. In this review, the basic principles of sonochemotherapy are discussed, including aspects such as the type of (drug-loaded) microbubbles used, the routes of administration used in vivo, ultrasound devices and parameters, treatment schedules and safety issues. Finally, the clinical translation of sonochemotherapy is discussed, including the first clinical study using sonochemotherapy.

MRI methods for the evaluation of high intensity focused ultrasound tumor treatment: Current status and future needs

Thermal ablation with high intensity focused ultrasound (HIFU) is an emerging noninvasive technique for the treatment of solid tumors. HIFU treatment of malignant tumors requires accurate treatment planning, monitoring and evaluation, which can be facilitated by performing the procedure in an MR-guided HIFU system. The MR-based evaluation of HIFU treatment is most often restricted to contrast-enhanced T1 -weighted imaging, while it has been shown that the non-perfused volume may not reflect the extent of nonviable tumor tissue after HIFU treatment. There are multiple studies in which more advanced MRI methods were assessed for their suitability for the evaluation of HIFU treatment. While several of these methods seem promising regarding their sensitivity to HIFU-induced tissue changes, there is still ample room for improvement of MRI protocols for HIFU treatment evaluation. In this review article, we describe the major acute and delayed effects of HIFU treatment. For each effect, the MRI methods that have been-or could be-used to detect the associated tissue changes are described. In addition, the potential value of multiparametric MRI for the evaluation of HIFU treatment is discussed. The review ends with a discussion on future directions for the MRI-based evaluation of HIFU treatment.

Cavitation-enhanced back projection for acoustic rib detection and attenuation mapping

High-intensity focused ultrasound allows for minimally invasive, highly localized cancer therapies that can complement surgical procedures or chemotherapy. For high-intensity focused ultrasound interventions in the upper abdomen, the thoracic cage obstructs and aberrates the ultrasonic beam, causing undesired heating of healthy tissue. When a phased array therapeutic transducer is used, such complications can be minimized by applying an apodization law based on analysis of beam path obstructions. In this work, a rib detection method based on cavitation-enhanced ultrasonic reflections is introduced and validated on a porcine tissue sample containing ribs. Apodization laws obtained for different transducer positions were approximately 90% similar to those obtained using image analysis. Additionally, the proposed method provides information on attenuation between transducer elements and the focus. This principle was confirmed experimentally on a polymer phantom. The proposed methods could, in principle, be implemented in real time for determination of the optimal shot position in intercostal high-intensity focused ultrasound therapy.

Influence of water and fat heterogeneity on fat-referenced MR thermometry

PURPOSE

To investigate the effect of the aqueous and fatty tissue magnetic susceptibility distribution on absolute and relative temperature measurements as obtained directly from the water/fat (w/f) frequency difference.

METHODS

Absolute thermometry was investigated using spherical phantoms filled with pork and margarine, which were scanned in three orthogonal orientations. To evaluate relative fat referencing, multigradient echo scans were acquired before and after heating pork tissue via high-intensity focused ultrasound (HIFU). Simulations were performed to estimate the errors that can be expected in human breast tissue.

RESULTS

The sphere experiment showed susceptibility-related errors of 8.4 °C and 0.2 °C for pork and margarine, respectively. For relative fat referencing measurements, fat showed pronounced phase changes of opposite polarity to aqueous tissue. The apparent mean temperature for a numerical breast model assumed to be 37 °C was 47.2 ± 21.6 °C. Simulations of relative fat referencing for a HIFU sonication (ΔT = 29.7 °C) yielded a maximum temperature error of 6.6 °C compared with 2.5 °C without fat referencing.

CONCLUSION

Variations in the observed frequency difference between water and fat are largely due to variations in the w/f spatial distribution. This effect may lead to considerable errors in absolute MR thermometry. Additionally, fat referencing may exacerbate rather than correct for proton resonance frequency shift-temperature measurement errors.

Understanding ultrasound induced sonoporation: definitions and underlying mechanisms

In the past two decades, research has underlined the potential of ultrasound and microbubbles to enhance drug delivery. However, there is less consensus on the biophysical and biological mechanisms leading to this enhanced delivery. Sonoporation, i.e. the formation of temporary pores in the cell membrane, as well as enhanced endocytosis is reported. Because of the variety of ultrasound settings used and corresponding microbubble behavior, a clear overview is missing. Therefore, in this review, the mechanisms contributing to sonoporation are categorized according to three ultrasound settings: i) low intensity ultrasound leading to stable cavitation of microbubbles, ii) high intensity ultrasound leading to inertial cavitation with microbubble collapse, and iii) ultrasound application in the absence of microbubbles. Using low intensity ultrasound, the endocytotic uptake of several drugs could be stimulated, while short but intense ultrasound pulses can be applied to induce pore formation and the direct cytoplasmic uptake of drugs. Ultrasound intensities may be adapted to create pore sizes correlating with drug size. Small molecules are able to diffuse passively through small pores created by low intensity ultrasound treatment. However, delivery of larger drugs such as nanoparticles and gene complexes, will require higher ultrasound intensities in order to allow direct cytoplasmic entry.

Magnetic resonance-guided high-intensity focused ultrasound (MR-HIFU) ablation of liver tumours

Recent decades have seen a paradigm shift in the treatment of liver tumours from invasive surgical procedures to minimally invasive image-guided ablation techniques. Magnetic resonance-guided high-intensity focused ultrasound (MR-HIFU) is a novel, completely non-invasive ablation technique that has the potential to change the field of liver tumour ablation. The image guidance, using MR imaging and MR temperature mapping, provides excellent planning images and real-time temperature information during the ablation procedure. However, before clinical implementation of MR-HIFU for liver tumour ablation is feasible, several organ-specific challenges have to be addressed. In this review we discuss the MR-HIFU ablation technique, the liver-specific challenges for MR-HIFU tumour ablation, and the proposed solutions for clinical translation.

Feasibility of fast MR-thermometry during cardiac radiofrequency ablation

Online MR temperature monitoring during radiofrequency (RF) ablation of cardiac arrhythmias may improve the efficacy and safety of the treatment. MR thermometry at 1.5 T using the proton resonance frequency (PRF) method was assessed in 10 healthy volunteers under normal breathing conditions, using a multi-slice, ECG-gated, echo planar imaging (EPI) sequence in combination with slice tracking. Temperature images were post-processed to remove residual motion-related artifacts. Using an MR-compatible steerable catheter and electromagnetic noise filter, RF ablation was performed in the ventricles of two sheep in vivo. The standard deviation of the temperature evolution in time (TSD) was computed. Temperature mapping of the left ventricle was achieved at an update rate of approximately 1 Hz with a mean TSD of 3.6 ± 0.9 °C. TSD measurements at the septum showed a higher precision (2.8 ± 0.9 °C) than at the myocardial regions at the heart-lung and heart-liver interfaces (4.1 ± 0.9 °C). Temperature rose maximally by 9 °C and 16 °C during 5 W and 10 W RF applications, respectively, for 60 s each. Tissue temperature can be monitored at an update rate of approximately 1 Hz in five slices. Typical temperature changes observed during clinical RF application can be monitored with an acceptable level of precision.

Robust real-time-constrained estimation of respiratory motion for interventional MRI on mobile organs

Real-time magnetic resonance imaging is a promising tool for image-guided interventions. For applications such as thermotherapy on moving organs, a precise image-based compensation of motion is required in real time to allow quantitative analysis, retrocontrol of the interventional device, or determination of the therapy endpoint. Reduced field-of-view imaging represents a promising way to improve spatial and/or temporal resolution. However, it introduces new challenges for target motion estimation, since structures near the target may appear transiently due to the respiratory motion and the limited spatial coverage. In this paper, a new image-based motion estimation method is proposed combining a global motion estimation with a novel optical flow approach extending the initial Horn and Schunck (H&S) method by an additional regularization term. This term integrates the displacement of physiological landmarks into the variational formulation of the optical flow problem. This allowed for a better control of the optical flow in presence of transient structures. The method was compared to the same registration pipeline employing the H&S approach on a synthetic dataset and in vivo image sequences. Compared to the H&S approach, a significant improvement (p<0.05) of the Dice's similarity criterion computed between the reference and the registered organ positions was achieved.