Magnetic resonance shear wave elastography using transient acoustic radiation force excitations and sinusoidal displacement encoding.

A magnetic resonance (MR) shear wave elastography technique that uses transient acoustic radiation force impulses from a focused ultrasound (FUS) transducer and a sinusoidal-shaped MR displacement encoding strategy is presented. Using this encoding strategy, an analytic expression for calculating the shear wave speed in a heterogeneous medium was derived. Green's function-based simulations were used to evaluate the feasibility of calculating shear wave speed maps using the analytic expression. Accuracy of simulation technique was confirmed experimentally in a homogeneous gelatin phantom. The elastography measurement was compared to harmonic MR elastography in a homogeneous phantom experiment and the measured shear wave speed values differed by less than 14%. This new transient elastography approach was able to map the position and shape of inclusions sized from 8.5 to 14 mm in an inclusion phantom experiment. These preliminary results demonstrate the feasibility of using a straightforward analytic expression to generate shear wave speed maps from MR images where sinusoidal-shaped motion encoding gradients are used to encode the displacement-time history of a transiently propagating wave-packet. This new measurement technique may be particularly well suited for performing elastography before, during, and after MR-guided FUS therapies since the same device used for therapy is also used as an excitation source for elastography.

Development and characterization of a tissue mimicking psyllium husk gelatin phantom for ultrasound and magnetic resonance imaging.

Purpose: To develop and characterize a tissue-mimicking phantom that enables the direct comparison of magnetic resonance (MR) and ultrasound (US) imaging techniques useful for monitoring high-intensity focused ultrasound (HIFU) treatments. With no additions, gelatin phantoms produce little if any scattering required for US imaging. This study characterizes the MR and US image characteristics as a function of psyllium husk concentration, which was added to increase US scattering.Methods: Gelatin phantoms were constructed with varying concentrations of psyllium husk. The effects of psyllium husk concentration on US B-mode and MR imaging were evaluated at nine different concentrations. T1, T2, and T2* MR maps were acquired. Acoustic properties (attenuation and speed of sound) were measured at frequencies of 0.6, 1.0, 1.8, and 3.0 MHz using a through-transmission technique. Phantom elastic properties were evaluated for both time and temperature dependence.Results: Ultrasound image echogenicity increased with increasing psyllium husk concentration while quality of gradient-recalled echo MR images decreased with increasing concentration. For all phantoms, the measured speed of sound ranged between 1567-1569 m/s and the attenuation ranged between 0.42-0.44 dB/(cm·MHz). Measured T1 ranged from 974-1051 ms. The T2 and T2* values ranged from 97-108 ms and 48-88 ms, respectively, with both showing a decreasing trend with increased psyllium husk concentration. Phantom stiffness, measured using US shear-wave speed measurements, increased with age and decreased with increasing temperature.Conclusions: The presented dual-use tissue-mimicking phantom is easy to manufacture and can be used to compare and evaluate US-guided and MR-guided HIFU imaging protocols.

Efficient shear wave elastography using transient acoustic radiation force excitations and MR displacement encoding.

PURPOSE: To present a novel MR shear wave elastography (MR-SWE) method that efficiently measures the speed of propagating wave packets generated using acoustic radiation force (ARF) impulses. METHODS: ARF impulses from a focused ultrasound (FUS) transducer were applied sequentially to a preselected set of positions and motion encoded MRI was used to acquire volumetric images of the propagating shear wavefront emanating from each point. The wavefront position at multiple propagation times was encoded in the MR phase image using a train of motion encoding gradient lobes. Generating a transient propagating wavefront at multiple spatial positions and sampling each at multiple time-points allowed for shear wave speed maps to be efficiently created. MR-SWE was evaluated in tissue mimicking phantoms and ex vivo bovine liver tissue before and after ablation. RESULTS: MR-SWE maps, covering an in-plane area of ~5 × 5 cm, were acquired in 12 s for a single slice and 144 s for a volumetric scan. MR-SWE detected inclusions of differing stiffness in a phantom experiment. In bovine liver, mean shear wave speed significantly increased from 1.65 ± 0.18 m/s in normal to 2.52 ± 0.18 m/s in ablated region (n = 581 pixels; P-value < 0.001). CONCLUSION: MR-SWE is an elastography technique that enables precise targeting and excitation of the desired tissue of interest. MR-SWE may be particularly well suited for treatment planning and endpoint assessment of MR-guided FUS procedures because the same device used for therapy can be used as an excitation source for tissue stiffness quantification.

Validation of hybrid angular spectrum acoustic and thermal modelling in phantoms.

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

Phase aberration simulation study of MRgFUS breast treatments.

PURPOSE: This simulation study evaluates the effects of phase aberration in breast MR-guided focused ultrasound (MRgFUS) ablation treatments performed with a phased-array transducer positioned laterally to the breast. A quantification of these effects in terms of thermal dose delivery and the potential benefits of phase correction is demonstrated in four heterogeneous breast numerical models. METHODS: To evaluate the effects of varying breast tissue properties on the quality of the focus, four female volunteers with confirmed benign fibroadenomas were imaged using 3T MRI. These images were segmented into numerical models with six tissue types, with each tissue type assigned standard acoustic properties from the literature. Simulations for a single-plane 16-point raster-scan treatment trajectory centered in a fibroadenoma in each modeled breast were performed for a breast-specific MRgFUS system. At each of the 16 points, pressure patterns both with and without applying a phase correction technique were determined with the hybrid-angular spectrum method. Corrected phase patterns were obtained using a simulation-based phase aberration correction technique to adjust each element's transmit phase to obtain maximized constructive interference at the desired focus. Thermal simulations were performed for both the corrected and uncorrected pressure patterns using a finite-difference implementation of the Pennes bioheat equation. The effect of phase correction was evaluated through comparison of thermal dose accumulation both within and outside a defined treatment volume. Treatment results using corrected and uncorrected phase aberration simulations were compared by evaluating the power required to achieve a 20 °C temperature rise at the first treatment location. The extent of the volumes that received a minimum thermal dose of 240 CEM at 43 °C inside the intended treatment volume as well as the volume in the remaining breast tissues was also evaluated in the form of a dose volume ratio (DVR), a DVR percent change between corrected and uncorrected phases, and an additional metric that measured phase spread. RESULTS: With phase aberration correction applied, there was an improvement in the focus for all breast anatomies as quantified by a reduction in power required (13%-102%) to reach 20 °C when compared to uncorrected simulations. Also, the DVR percent change increased by 5%-77% in seven out of eight cases, indicating an improvement to the treatment as measured by a reduction in thermal dose deposited to the nontreatment tissues. Breast compositions with a higher degree of heterogeneity along the ultrasound beam path showed greater reductions in thermal dose delivered outside of the treatment volume with correction applied than beam trajectories that propagated through more homogeneous breast compositions. An increasing linear trend was observed between the DVR percent change and the phase-spread metric (R(2) = 0.68). CONCLUSIONS: These results indicate that performing phase aberration correction for breast MRgFUS treatments is beneficial for the small-aperture transducer (14.4 × 9.8 cm) evaluated in this work. While all breast anatomies could benefit from phase aberration correction, greater benefits are observed in more heterogeneous anatomies.

A simulation technique for 3D MR-guided acoustic radiation force imaging.

PURPOSE: In magnetic resonance-guided focused ultrasound (MRgFUS) therapies, the in situ characterization of the focal spot location and quality is critical. MR acoustic radiation force imaging (MR-ARFI) is a technique that measures the tissue displacement caused by the radiation force exerted by the ultrasound beam. This work presents a new technique to model the displacements caused by the radiation force of an ultrasound beam in a homogeneous tissue model. METHODS: When a steady-state point-source force acts internally in an infinite homogeneous medium, the displacement of the material in all directions is given by the Somigliana elastostatic tensor. The radiation force field, which is caused by absorption and reflection of the incident ultrasound intensity pattern, will be spatially distributed, and the tensor formulation takes the form of a convolution of a 3D Green's function with the force field. The dynamic accumulation of MR phase during the ultrasound pulse can be theoretically accounted for through a time-of-arrival weighting of the Green's function. This theoretical model was evaluated experimentally in gelatin phantoms of varied stiffness (125-, 175-, and 250-bloom). The acoustic and mechanical properties of the phantoms used as parameters of the model were measured using independent techniques. Displacements at focal depths of 30- and 45-mm in the phantoms were measured by a 3D spin echo MR-ARFI segmented-EPI sequence. RESULTS: The simulated displacements agreed with the MR-ARFI measured displacements for all bloom values and focal depths with a normalized RMS difference of 0.055 (range 0.028-0.12). The displacement magnitude decreased and the displacement pattern broadened with increased bloom value for both focal depths, as predicted by the theory. CONCLUSIONS: A new technique that models the displacements caused by the radiation force of an ultrasound beam in a homogeneous tissue model theory has been rigorously validated through comparison with experimentally obtained 3D displacement data in homogeneous gelatin phantoms using a 3D MR-ARFI sequence. The agreement of the experimentally measured and simulated results demonstrates the potential to use MR-ARFI displacement data in MRgFUS therapies.

Treatment envelope evaluation in transcranial magnetic resonance-guided focused ultrasound utilizing 3D MR thermometry.

BACKGROUND: Current clinical targets for transcranial magnetic resonance-guided focused ultrasound (tcMRgFUS) are all located close to the geometric center of the skull convexity, which minimizes challenges related to focusing the ultrasound through the skull bone. Non-central targets will have to be reached to treat a wider variety of neurological disorders and solid tumors. Treatment envelope studies utilizing two-dimensional (2D) magnetic resonance (MR) thermometry have previously been performed to determine the regions in which therapeutic levels of FUS can currently be delivered. Since 2D MR thermometry was used, very limited information about unintended heating in near-field tissue/bone interfaces could be deduced. METHODS: In this paper, we present a proof-of-concept treatment envelope study with three-dimensional (3D) MR thermometry monitoring of FUS heatings performed in a phantom and a lamb model. While the moderate-sized transducer used was not designed for transcranial geometries, the 3D temperature maps enable monitoring of the entire sonication field of view, including both the focal spot and near-field tissue/bone interfaces, for full characterization of all heating that may occur. 3D MR thermometry is achieved by a combination of k-space subsampling and a previously described temporally constrained reconstruction method. RESULTS: We present two different types of treatment envelopes. The first is based only on the focal spot heating-the type that can be derived from 2D MR thermometry. The second type is based on the relative near-field heating and is calculated as the ratio between the focal spot heating and the near-field heating. This utilizes the full 3D MR thermometry data achieved in this study. CONCLUSIONS: It is shown that 3D MR thermometry can be used to improve the safety assessment in treatment envelope evaluations. Using a non-optimal transducer, it is shown that some regions where therapeutic levels of FUS can be delivered, as suggested by the first type of envelope, are not necessarily safely treated due to the amount of unintended near-field heating occurring. The results presented in this study highlight the need for 3D MR thermometry in tcMRgFUS.

Adaptive model-predictive controller for magnetic resonance guided focused ultrasound therapy.

PURPOSE: Minimising treatment time and protecting healthy tissues are conflicting goals that play major roles in making magnetic resonance image-guided focused ultrasound (MRgFUS) therapies clinically practical. We have developed and tested in vivo an adaptive model-predictive controller (AMPC) that reduces treatment time, ensures safety and efficacy, and provides flexibility in treatment set-up. MATERIALS AND METHODS: The controller realises time savings by modelling the heated treatment cell's future temperatures and thermal dose accumulation in order to anticipate the optimal time to switch to the next cell. Selected tissues are safeguarded by a configurable temperature constraint. Simulations quantified the time savings realised by each controller feature as well as the trade-offs between competing safety and treatment time parameters. In vivo experiments in rabbit thighs established the controller's effectiveness and reliability. RESULTS: In all in vivo experiments the target thermal dose of at least 240 CEM43 was delivered everywhere in the treatment volume. The controller's temperature safety limit reliably activated and constrained all protected tissues to <9 CEM43. Simulations demonstrated the path independence of the controller, and that a path which successively proceeds to the hottest untreated neighbouring cell leads to significant time savings, e.g. when compared to a concentric spiral path. Use of the AMPC produced a compounding time-saving effect; reducing the treatment cells' heating times concurrently reduced heating of normal tissues, which eliminated cooling periods. CONCLUSIONS: Adaptive model-predictive control can automatically deliver safe, effective MRgFUS treatments while significantly reducing treatment times.

The accuracy and precision of two non-invasive, magnetic resonance-guided focused ultrasound-based thermal diffusivity estimation methods.

PURPOSE: The use of correct tissue thermal diffusivity values is necessary for making accurate thermal modelling predictions during magnetic resonance-guided focused ultrasound (MRgFUS) treatment planning. This study evaluates the accuracy and precision of two non-invasive thermal diffusivity estimation methods, a Gaussian temperature method and a Gaussian specific absorption rate (SAR) method. MATERIALS AND METHODS: Both methods utilise MRgFUS temperature data obtained during cooling following a short (<25 s) heating pulse. The Gaussian SAR method can also use temperatures obtained during heating. Experiments were performed at low heating levels (ΔT∼10 °C) in ex vivo pork muscle and in vivo rabbit back muscle. The non-invasive MRgFUS thermal diffusivity estimates were compared with measurements from two standard invasive methods. RESULTS: Both non-invasive methods accurately estimated thermal diffusivity when using MR temperature cooling data (overall ex vivo error <6%, in vivo <12%). Including heating data in the Gaussian SAR method further reduced errors (ex vivo error <2%, in vivo <3%). The significantly lower standard deviation values (p < 0.03) of the Gaussian SAR method indicated that it had better precision than the Gaussian temperature method. CONCLUSIONS: With repeated sonications, either MR-based method could provide accurate thermal diffusivity values for MRgFUS therapies. Fitting to more data simultaneously likely made the Gaussian SAR method less susceptible to noise, and using heating data helped it converge more consistently to the FUS fitting parameters and thermal diffusivity. These effects led to the improved precision of the Gaussian SAR method.

Investigating the acoustic release of doxorubicin from targeted micelles.

The main problem associated with the administration of anti-cancer medication is that the drug is delivered throughout the body causing undesirable side effects. Therefore, it is important to synthesize drug carriers capable of minimizing the adverse side effects of chemotherapy by preferentially targeting tumor cells both actively (e.g. a folate receptor) and using external stimulus (e.g. ultrasound). In this paper, we report the synthesis of Pluronic P105 micelles with a folate targeting moiety (with a yield of 48%) containing doxorubicin (Dox). We applied low frequency ultrasound as an external stimulus and measured the amount of release of Dox from these folated micelles. The results showed that the percent drug release increases as the power intensity of ultrasound increases. The maximum amount of release (14%) was measured at 5.4 W/cm(2). A power density threshold at approximately 0.55 W/cm(2) exists below which no statistically significant release was observed. This lower threshold suggests that cavitation plays an important role in triggering drug release from targeted micelles.

Thermal sensitivity of endothelial cells on synthetic vascular graft material.

PURPOSE: The goal is to identify thermal exposures capable of reducing or eliminating cell survival on expanded polytetrafluoroethylene (ePTFE), in an effort to develop a mild hyperthermia treatment of neointimal hyperplasia in ePTFE vascular grafts. MATERIALS AND METHODS: Viable and dead bovine aortic endothelial cells were quantified following different thermal exposure conditions: cells on collagen-coated ePTFE sheets or tissue culture polystyrene dishes were heated at 42° and 45°C to determine their thermal sensitivity on different surfaces, and cells cultured on collagen-coated ePTFE sheets were heated at 43-50°C for various durations, followed by incubation at 37°C for 0 and 20 h, respectively. Significant cell death was set to be 50%. Two types of cell death, apoptosis and necrosis, were distinguished by cell morphology and membrane integrity assessments. RESULTS: The attachment and survival of cells on ePTFE sheets were more sensitive to inhibition by mild heating than those on tissue culture dishes. Exposure to 45°C for 90 min and 50°C for 30 min caused significant necrotic cell death on ePTFE (65% and 75%, respectively). A 37°C/20-h incubation following 30-min exposures at 47° and 50°C increased total cell death (necrosis + apoptosis) from 20% to 50% and 75% to 100%, respectively. CONCLUSION: Cells grown on ePTFE were more susceptible to mild hyperthermia-induced death, compared to those on tissue culture dishes. Significant cell death on ePTFE mainly via apoptosis can be achieved by optimising temperature and duration of exposure.

Extension of the angular spectrum method to calculate pressure from a spherically curved acoustic source.

The angular spectrum method is an accurate and computationally efficient method for modeling acoustic wave propagation. The use of the typical 2D fast Fourier transform algorithm makes this a fast technique but it requires that the source pressure (or velocity) be specified on a plane. Here the angular spectrum method is extended to calculate pressure from a spherical transducer-as used extensively in applications such as magnetic resonance-guided focused ultrasound surgery-to a plane. The approach, called the Ring-Bessel technique, decomposes the curved source into circular rings of increasing radii, each ring a different distance from the intermediate plane, and calculates the angular spectrum of each ring using a Fourier series. Each angular spectrum is then propagated to the intermediate plane where all the propagated angular spectra are summed to obtain the pressure on the plane; subsequent plane-to-plane propagation can be achieved using the traditional angular spectrum method. Since the Ring-Bessel calculations are carried out in the frequency domain, it reduces calculation times by a factor of approximately 24 compared to the Rayleigh-Sommerfeld method and about 82 compared to the Field II technique, while maintaining accuracies of better than 96% as judged by those methods for cases of both solid and phased-array transducers.

Microbubble Generation in Phase-Shift Nanoemulsions used as Anticancer Drug Carriers.

The paper describes droplet-to-bubble transition in block copolymer stabilized perfluoropentane nanoemulsions. Three physical factors that trigger droplet-to-bubble transition in liquid emulsions and gels were evaluated, namely heat, ultrasound, and injections through fine-gauge needles. Among those listed, ultrasound irradiation was found the most efficient factor. Possible mechanisms of bubble generation and growth discussed in the paper include liquid-to-gas transition inside the individual bubble; bubble coalescence; and diffusion of dissolved air and/or perfluoropentane from small bubbles into larger bubbles (i.e., Oswald ripening). The last two factors result in irreversibility of the droplet-to-bubble transition. In gel matrices, ultrasound-induced droplet-to-bubble transition was substantially inhibited but was catalyzed by large (hundred micron) pre-existing bubbles irradiated by low frequency (hundred kilohertz) ultrasound. The dependence of the droplet-to-bubble transition on initial bubble size is theoretically treated and the role of increase of surface area in promoting bubble coalescence is discussed. Therapeutic implications of observed effects are discussed.

Drug-loaded nano/microbubbles for combining ultrasonography and targeted chemotherapy.

A new class of multifunctional nanoparticles that combine properties of polymeric drug carriers, ultrasound imaging contrast agents, and enhancers of ultrasound-mediated drug delivery has been developed. At room temperature, the developed systems comprise perfluorocarbon nanodroplets stabilized by the walls made of biodegradable block copolymers. Upon heating to physiological temperatures, the nanodroplets convert into nano/microbubbles. The phase state of the systems and bubble size may be controlled by the copolymer/perfluorocarbon volume ratio. Upon intravenous injections, a long-lasting, strong and selective ultrasound contrast is observed in the tumor volume indicating nanobubble extravasation through the defective tumor microvasculature, suggesting their coalescence into larger, highly echogenic microbubbles in the tumor tissue. Under the action of tumor-directed ultrasound, microbubbles cavitate and collapse resulting in a release of the encapsulated drug and dramatically enhanced intracellular drug uptake by the tumor cells. This effect is tumor-selective; no accumulation of echogenic microbubbles is observed in other organs. Effective chemotherapy of the MDA MB231 breast cancer tumors has been achieved using this technique.

Release of doxorubicin from unstabilized and stabilized micelles under the action of ultrasound.

Polymeric micelles are being investigated as chemotherapy drug delivery carriers using ultrasound as a trigger mechanism. The aim of this paper is to measure the release of Doxorubicin (Dox) from the core of unstabilized Pluronic P105 micelles, Pluronic P105 micelles stabilized with an interpenetrating network of N,N-diethylacrylamide, and micelles of poly(ethylene oxide)-b-poly (N-isopropylacrylamide)-b-poly(oligolactylmethacrylate) with stabilized cores. An ultrasonic exposure chamber with fluorescence detection was used to measure the release of the antineoplastic agent from both stabilized and unstabilized micelles. The release of Dox at 37 degrees C from unstabilized Pluronic appears to be several times higher than release from the more stabilized and crosslinked copolymers at the same temperature. Although there is a difference in the amount of release between the different compounds, the onset of release occurs at about the same ultrasonic power density for all carriers investigated in this study. The threshold of drug release for all the compounds correlates to the emergence of subharmonic peaks detected in the acoustic spectra. We hypothesize that shearing events caused by cavitating bubbles play an important role in the acoustically activated release of chemotherapy agents delivered from various polymeric drug delivery vehicles.

The role of cavitation in acoustically activated drug delivery.

Pluronic P105 micelles are potential candidates as chemotherapy drug delivery vehicles using ultrasonic stimulation as a release trigger. Acoustic power has been previously shown to release two anthracycline agents from these polymeric carriers. In this study, an ultrasonic exposure chamber with fluorescence detection was used to examine the mechanism of doxorubicin release from P105 micelles. Acoustic spectra were collected and analyzed, at the same spatial position as fluorescence data, to probe the role of cavitation in drug release. Our study showed a strong correlation between percent drug release and subharmonic acoustic emissions, and we attribute the drug release to collapse cavitation that perturbs the structure of the micelle and releases drug.

Drug delivery in pluronic micelles: effect of high-frequency ultrasound on drug release from micelles and intracellular uptake.

The effect of high-frequency ultrasound on doxorubicin (DOX) release from Pluronic micelles and intracellular DOX uptake was studied for promyelocytic leukemia HL-60 cells, ovarian carcinoma drug-sensitive and multidrug-resistant (MDR) cells (A2780 and A2780/ADR, respectively), and breast cancer MCF-7 cells. Cavitation events initiated by high-frequency ultrasound were recorded by radical trapping. The onset of transient cavitation and DOX release from micelles were observed at much higher power densities than at low-frequency ultrasound (20-100 kHz). Even a short (15-30 s) exposure to high-frequency ultrasound significantly enhanced the intracellular DOX uptake from PBS, RPMI 1640, and Pluronic micelles. The mechanisms of the observed effects are discussed.