Development of an MR-Guided Focused Ultrasound (MRgFUS) Lesioning Approach for the Fornix in the Rat Brain.

OBJECTIVE: High-intensity magnetic resonance-guided focused ultrasound (MRgFUS) is a non-invasive therapy to lesion brain tissue, used clinically in patients and pre-clinically in several animal models. Challenges with focused ablation in rodent brains can include skull and near-field heating and accurately targeting small and deep brain structures. We overcame these challenges by creating a novel method consisting of a craniectomy skull preparation, a high-frequency transducer (3 MHz) with a small ultrasound focal spot, a transducer positioning system with an added manual adjustment of ∼0.1 mm targeting accuracy, and MR acoustic radiation force imaging for confirmation of focal spot placement. METHODS: The study consisted of two main parts. First, two skull preparation approaches were compared. A skull thinning approach (n = 7 lesions) was compared to a craniectomy approach (n = 22 lesions), which confirmed a craniectomy was necessary to decrease skull and near-field heating. Second, the two transducer positioning systems were compared with the fornix chosen as a subcortical ablation target. We evaluated the accuracy of targeting using histologic methods from a high-frequency transducer with a small ultrasound focal spot and MR acoustic radiation force imaging. RESULTS: Comparing a motorized adjustment system (∼1 mm precision, n = 17 lesions) to the motorized system with an added micromanipulator (∼0.1 mm precision, n = 14 lesions), we saw an increase in the accuracy of targeting the fornix by 133%. CONCLUSIONS: The described work allows for repeatable and accurate targeting of small and deep structures in the rodent brain, such as the fornix, enabling the investigation of neurological disorders in chronic disease models.

REhabilitation Approaches in CHildren with cerebellar mutism syndrome (REACH): An international cross-disciplinary survey study.

PURPOSE: Pediatric cerebellar mutism syndrome (pCMS) can occur following resection of a posterior fossa tumor and, although some symptoms are transient, many result in long-lasting neurological deficits. A multi-disciplinary rehabilitation approach is often used in cases of pCMS; however, there have been no clinical trials to determine gold standards in rehabilitation practice in this population, which remains a research priority. The purpose of this study was to identify and compare intervention practices used in pCMS throughout the disciplines of occupational and physical therapy, speech-language pathology, and neuropsychology across geographic regions. METHODS: A 55-question e-survey was created by an international multidisciplinary research group made up of members of the Posterior Fossa Society and sent to rehabilitation professionals in pediatric neuro-oncology centers in the US, Canada, and Europe. RESULTS: Although some differences in the type of intervention used in pCMS were identified within each discipline, many of the targeted interventions including dose, frequency, and intensity were similar within disciplines across geographic regions. In addition, there were common themes identified across disciplines regarding challenges in the rehabilitation of this population. CONCLUSION: These results provide a foundation of current practices on which to build future intervention-based clinical trials.

Evaluation of acoustic-thermal simulations of in vivo magnetic resonance guided focused ultrasound ablative therapy.

PURPOSE: To evaluate numerical simulations of focused ultrasound (FUS) with a rabbit model, comparing simulated heating characteristics with magnetic resonance temperature imaging (MRTI) data collected during in vivo treatment. METHODS: A rabbit model was treated with FUS sonications in the biceps femoris with 3D MRTI collected. Acoustic and thermal properties of the rabbit muscle were determined experimentally. Numerical models of the rabbits were created, and tissue-type-specific properties were assigned. FUS simulations were performed using both the hybrid angular spectrum (HAS) method and k-Wave. Simulated power deposition patterns were converted to temperature maps using a Pennes' bioheat equation-based thermal solver. Agreement of pressure between the simulation techniques and temperature between the simulation and experimental heating was evaluated. Contributions of scattering and absorption attenuation were considered. RESULTS: Simulated peak pressures derived using the HAS method exceeded the simulated peak pressures from k-Wave by 1.6 ± 2.7%. The location and FWHM of the peak pressure calculated from HAS and k-Wave showed good agreement. When muscle acoustic absorption value in the simulations was adjusted to approximately 54% of the measured attenuation, the average root-mean-squared error between simulated and experimental spatial-average temperature profiles was 0.046 ± 0.019 °C/W. Mean distance between simulated and experimental COTMs was 3.25 ± 1.37 mm. Transverse FWHMs of simulated sonications were smaller than in in vivo sonications. Longitudinal FWHMs were similar. CONCLUSIONS: Presented results demonstrate agreement between HAS and k-Wave simulations and that FUS simulations can accurately predict focal position and heating for in vivo applications in soft tissue.

A Non-Contrast Multi-Parametric MRI Biomarker for Assessment of MR-Guided Focused Ultrasound Thermal Therapies.

OBJECTIVE: We present the development of a non-contrast multi-parametric magnetic resonance (MPMR) imaging biomarker to assess treatment outcomes for magnetic resonance-guided focused ultrasound (MRgFUS) ablations of localized tumors. Images obtained immediately following MRgFUS ablation were inputs for voxel-wise supervised learning classifiers, trained using registered histology as a label for thermal necrosis. METHODS: VX2 tumors in New Zealand white rabbits quadriceps were thermally ablated using an MRgFUS system under 3 T MRI guidance. Animals were re-imaged three days post-ablation and euthanized. Histological necrosis labels were created by 3D registration between MR images and digitized H&E segmentations of thermal necrosis to enable voxel-wise classification of necrosis. Supervised MPMR classifier inputs included maximum temperature rise, cumulative thermal dose (CTD), post-FUS differences in T2-weighted images, and apparent diffusion coefficient, or ADC, maps. A logistic regression, support vector machine, and random forest classifier were trained in red a leave-one-out strategy in test data from four subjects. RESULTS: In the validation dataset, the MPMR classifiers achieved higher recall and Dice than a clinically adopted 240 cumulative equivalent minutes at 43 °C (CEM 43) threshold (0.43) in all subjects. The average Dice scores of overlap with the registered histological label for the logistic regression (0.63) and support vector machine (0.63) MPMR classifiers were within 6% of the acute contrast-enhanced non-perfused volume (0.67). CONCLUSIONS: Voxel-wise registration of MPMR data to histological outcomes facilitated supervised learning of an accurate non-contrast MR biomarker for MRgFUS ablations in a rabbit VX2 tumor model.

Development of an MR-guided focused ultrasound (MRgFUS) lesioning approach for small and deep structures in the rat brain.

Objective: High-intensity magnetic resonance-guided focused ultrasound (MRgFUS) is a noninvasive therapy to lesion brain tissue, used clinically in patients and preclinically in several animal models. Challenges with focused ablation in rodent brains can include skull and near-field heating and accurately targeting small and deep brain structures. We overcame these challenges by creating a novel method consisting of a craniectomy skull preparation, a high-frequency transducer (3 MHz) with a small ultrasound focal spot, a transducer positioning system with an added manual adjustment of ∼0.1 mm targeting accuracy, and MR acoustic radiation force imaging for confirmation of focal spot placement. Methods: The study consisted of two main parts. First, two skull preparation approaches were compared. A skull thinning approach (n=7 lesions) was compared to a craniectomy approach (n=22 lesions), which confirmed a craniectomy was necessary to decrease skull and near-field heating. Second, the two transducer positioning systems were compared with the fornix chosen as a subcortical ablation target. We evaluated the accuracy of targeting using a high-frequency transducer with a small ultrasound focal spot and MR acoustic radiation force imaging. Results: Comparing a motorized adjustment system (∼1 mm precision, n=17 lesions) to the motorized system with an added micromanipulator (∼0.1 mm precision, n=14 lesions), we saw an increase in the accuracy of targeting the fornix by 133%. The described work allows for repeatable and accurate targeting of small and deep structures in the rodent brain, such as the fornix, enabling the investigation of neurological disorders in chronic disease models.

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

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

Experimental validation of acoustic and thermal modeling in heterogeneous phantoms using the hybrid angular spectrum method.

PURPOSE: The aim was to quantitatively validate the hybrid angular spectrum (HAS) algorithm, a rapid wave propagation technique for heterogeneous media, with both pressure and temperature measurements. METHODS: Heterogeneous tissue-mimicking phantoms were used to evaluate the accuracy of the HAS acoustic modeling algorithm in predicting pressure and thermal patterns. Acoustic properties of the phantom components were measured by a through-transmission technique while thermal properties were measured with a commercial probe. Numerical models of each heterogeneous phantom were segmented from 3D MR images. Cylindrical phantoms 30-mm thick were placed in the pre-focal field of a focused ultrasound beam and 2D pressure measurements obtained with a scanning hydrophone. Peak pressure, full width at half maximum, and normalized root mean squared difference (RMSDn) between the measured and simulated patterns were compared. MR-guided sonications were performed on 150-mm phantoms to obtain MR temperature measurements. Using HAS-predicted power density patterns, temperature simulations were performed. Experimental and simulated temperature patterns were directly compared using peak and mean temperature plots, RMSDn metrics, and accuracy of heating localization. RESULTS: The average difference between simulated and hydrophone-measured peak pressures was 9.0% with an RMSDn of 11.4%. Comparison of the experimental MRI-derived and simulated temperature patterns showed RMSDn values of 10.2% and 11.1% and distance differences between the centers of thermal mass of 2.0 and 2.2 mm. CONCLUSIONS: These results show that the computationally rapid hybrid angular spectrum method can predict pressure and temperature patterns in heterogeneous models, including uncertainties in property values and other parameters, to within approximately 10%.

Histology to 3D in vivo MR registration for volumetric evaluation of MRgFUS treatment assessment biomarkers.

Advances in imaging and early cancer detection have increased interest in magnetic resonance (MR) guided focused ultrasound (MRgFUS) technologies for cancer treatment. MRgFUS ablation treatments could reduce surgical risks, preserve organ tissue and function, and improve patient quality of life. However, surgical resection and histological analysis remain the gold standard to assess cancer treatment response. For non-invasive ablation therapies such as MRgFUS, the treatment response must be determined through MR imaging biomarkers. However, current MR biomarkers are inconclusive and have not been rigorously evaluated against histology via accurate registration. Existing registration methods rely on anatomical features to directly register in vivo MR and histology. For MRgFUS applications in anatomies such as liver, kidney, or breast, anatomical features that are not caused by the treatment are often insufficient to drive direct registration. We present a novel MR to histology registration workflow that utilizes intermediate imaging and does not rely on anatomical MR features being visible in histology. The presented workflow yields an overall registration accuracy of 1.00 ± 0.13 mm. The developed registration pipeline is used to evaluate a common MRgFUS treatment assessment biomarker against histology. Evaluating MR biomarkers against histology using this registration pipeline will facilitate validating novel MRgFUS biomarkers to improve treatment assessment without surgical intervention. While the presented registration technique has been evaluated in a MRgFUS ablation treatment model, this technique could be potentially applied in any tissue to evaluate a variety of therapeutic options.

AAPM Task Group 241: A medical physicist's guide to MRI-guided focused ultrasound body systems.

Magnetic resonance-guided focused ultrasound (MRgFUS) is a completely non-invasive technology that has been approved by FDA to treat several diseases. This report, prepared by the American Association of Physicist in Medicine (AAPM) Task Group 241, provides background on MRgFUS technology with a focus on clinical body MRgFUS systems. The report addresses the issues of interest to the medical physics community, specific to the body MRgFUS system configuration, and provides recommendations on how to successfully implement and maintain a clinical MRgFUS program. The following sections describe the key features of typical MRgFUS systems and clinical workflow and provide key points and best practices for the medical physicist. Commonly used terms, metrics and physics are defined and sources of uncertainty that affect MRgFUS procedures are described. Finally, safety and quality assurance procedures are explained, the recommended role of the medical physicist in MRgFUS procedures is described, and regulatory requirements for planning clinical trials are detailed. Although this report is limited in scope to clinical body MRgFUS systems that are approved or currently undergoing clinical trials in the United States, much of the material presented is also applicable to systems designed for other applications.

A T1-based correction method for proton resonance frequency shift thermometry in breast tissue.

PURPOSE: Develop and evaluate the effectiveness of a T1-based correction method for errors in proton resonant frequency shift thermometry due to non-local field effects caused by heating in fatty breast tissues. METHODS: Computational models of human breast tissue were created by segmenting MRI data from a healthy human volunteer. MR-guided focused ultrasound (MRgFUS) heating and MR thermometry measurements were simulated in several locations in the heterogeneous segmented breast models. A T1-based correction method for PRF thermometry errors was applied and the maximum positive and negative errors and the root mean squared error (RMSE) in a region around each heating location was evaluated with and without correction. The method uses T1 measurements to estimate the temperature change in fatty tissues and correct for their influence. Experimental data from a heating study in cadaver breast tissue were analyzed, and the expected PRFS error computed. RESULTS: The simulated MR thermometry had maximum single voxel errors ranging between 10% and 18% when no correction was applied. Applying the correction led to a considerable improvement, lowering the maximum error range to 2%-5%. The 5th to 95th percentile interval of the temperature error distribution was also lowered with correction, from approximately 3.5 to 1°C. This correction worked even when T1 times were uniformly raised or lowered by 5%-10%. The experimental data showed predicted errors of 15%. CONCLUSIONS: This simulation study demonstrates that the T1-based correction method reduces MR thermometry errors due to non-local effects from heating in fatty tissues, potentially improving the accuracy of thermometry measurements during MRgFUS treatments. The presented correction method is reliant on having a patient-specific 3D model of the breast, and may be limited by the accuracy of the fat temperatures which in turn may be limited by noise or bias present in the T1 measurements.

Design and evaluation of an open-source, conformable skin-cooling system for body magnetic resonance guided focused ultrasound treatments.

PURPOSE: Magnetic resonance guided focused ultrasound (MRgFUS) treatment of tumors uses inter-sonication delays to allow heat to dissipate from the skin and other near-field tissues. Despite inter-sonication delays, treatment of tumors close to the skin risks skin burns. This work has designed and evaluated an open-source, conformable, skin-cooling system for body MRgFUS treatments to reduce skin burns and enable ablation closer to the skin. METHODS: A MR-compatible skin cooling system is described that features a conformable skin-cooling pad assembly with feedback control allowing continuous flow and pressure maintenance during the procedure. System performance was evaluated with hydrophone, phantom and in vivo porcine studies. Sonications were performed 10 and 5 mm from the skin surface under both control and forced convective skin-cooling conditions. 3D MR temperature imaging was acquired in real time and the accumulated thermal dose volume was measured. Gross analysis of the skin post-sonication was further performed. Device conformability was demonstrated at several body locations. RESULTS: Hydrophone studies demonstrated no beam aberration, but a 5-12% reduction of the peak pressure due to the presence of the skin-cooling pad assembly in the acoustic near field. Phantom evaluation demonstrated there is no MR temperature imaging precision reduction or any other artifacts present due to the coolant flow during MRgFUS sonication. The porcine studies demonstrated skin burns were reduced in size or eliminated when compared to the control condition. CONCLUSION: An open-source design of an MRgFUS active skin cooling system demonstrates device conformability with a reduction of skin burns while ablating superficial tissues.

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.

A Breast-Specific MR Guided Focused Ultrasound Platform and Treatment Protocol: First-in-Human Technical Evaluation.

OBJECTIVE: This paper presents and evaluates a breast-specific magnetic resonance guided focused ultrasound (MRgFUS) system. A first-in-human evaluation demonstrates the novel hardware, a sophisticated tumor targeting algorithm and a volumetric magnetic resonance imaging (MRI) protocol. METHODS: At the time of submission, N = 10 patients with non-palpable T0 stage breast cancer have been treated with the breast MRgFUS system. The described tumor targeting algorithm is evaluated both with a phantom test and in vivo during the breast MRgFUS treatments. Treatments were planned and monitored using volumetric MR-acoustic radiation force imaging (MR-ARFI) and temperature imaging (MRTI). RESULTS: Successful technical treatments were achieved in 80 % of the patients. All patients underwent the treatment with no sedation and 60 % of participants had analgesic support. The total MR treatment time ranged from 73 to 114 minutes. Mean error between desired and achieved targeting in a phantom was 2.9 ±1.8 mm while 6.2 ±1.9 mm was achieved in patient studies, assessed either with MRTI or MR-ARFI measurements. MRTI and MR-ARFI were successful in 60 % and 70 % of patients, respectively. CONCLUSION: The targeting accuracy allows the accurate placement of the focal spot using electronic steering capabilities of the transducer. The use of both volumetric MRTI and MR-ARFI provides complementary treatment planning and monitoring information during the treatment, allowing the treatment of all breast anatomies, including homogeneously fatty breasts.

Learning Multiparametric Biomarkers for Assessing MR-Guided Focused Ultrasound Treatment of Malignant Tumors.

Noninvasive MR-guided focused ultrasound (MRgFUS) treatments are promising alternatives to the surgical removal of malignant tumors. A significant challenge is assessing the viability of treated tissue during and immediately after MRgFUS procedures. Current clinical assessment uses the nonperfused volume (NPV) biomarker immediately after treatment from contrast-enhanced MRI. The NPV has variable accuracy, and the use of contrast agent prevents continuing MRgFUS treatment if tumor coverage is inadequate. This work presents a novel, noncontrast, learned multiparametric MR biomarker that can be used during treatment for intratreatment assessment, validated in a VX2 rabbit tumor model. A deep convolutional neural network was trained on noncontrast multiparametric MR images using the NPV biomarker from follow-up MR imaging (3-5 days after MRgFUS treatment) as the accurate label of nonviable tissue. A novel volume-conserving registration algorithm yielded a voxel-wise correlation between treatment and follow-up NPV, providing a rigorous validation of the biomarker. The learned noncontrast multiparametric MR biomarker predicted the follow-up NPV with an average DICE coefficient of 0.71, substantially outperforming the current clinical standard (DICE coefficient = 0.53). Noncontrast multiparametric MR imaging integrated with a deep convolutional neural network provides a more accurate prediction of MRgFUS treatment outcome than current contrast-based techniques.

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.

A tissue preparation to characterize uterine fibroid tissue properties for thermal therapies.

PURPOSE: Treating uterine fibroids with less invasive therapies such as magnetic resonance-guided focused ultrasound (MRgFUS) is an attractive alternative to surgery. Treatment planning can improve MRgFUS procedures and reduce treatment times, but the tissue properties that currently inform treatment planning tools are not adequate. This study aims to develop an ex vivo uterine fibroid model that can emulate the in vivo environment allowing for characterization of the uterus and fibroid MR, acoustic, and thermal tissue properties while maintaining viability for the necessary postsurgical histopathological assessments. METHODS: Women undergoing a hysterectomy due to fibroid-related symptoms were invited to undergo a preoperative pelvic MRI and to permit postoperative testing of their uterine specimen. Patients that declined or could not be scheduled for a pre-operative MRI were still able to allow post-operative testing of their excised tissue. Following surgical removal of the uterus, nonmorcellated tissues were reperfused with a Krebs-Henseleit buffer solution. An MR-compatible perfusion system was designed to maintain tissue viability inside the MR suite during scanning. MR imaging protocols utilized preoperatively were repeated on whole sample, reperfused ex vivo uterus specimens. Thermal properties including thermal diffusivity and thermal conductivity of the uterus and fibroids were determined using an invasive needle sensor device in 50% of the specimens. Acoustic property measurements (density, speed of sound and attenuation) were obtained for approximately 20% of the tissue samples using both through-transmission and radiation force balance techniques. Differences between fibroid and uterus and in vivo and ex vivo measurements were evaluated with a two-tailed Student t test. RESULTS: Fourteen patients participated in the study and measurements were obtained from 22 unique fibroids. Of the 16 fibroids available for preoperative MRI testing, 69% demonstrated classic hypo-intensity relative to the myometrium, with the remainder presenting with iso- (25%) or hyper-intensity (6%). While thermal diffusivity was not significantly different between fibroid and myometrium tissues (0.217 ± 0.047 and 0.204 ± 0.039 mm2 /s, respectively), the acoustic attenuation in fibroid tissue was significantly higher than myometrium (0.092 ± 0.021 and 0.052 ± 0.023 Np/cm/MHz, respectively). When comparing in vivo with ex vivo MRI T1 and T2 measurements in fibroids and myometrium tissue, the only difference was found in the fibroid T2 property (P < 0.05). Finally, the developed perfusion protocol successfully maintained tissue viability in ex vivo tissues as evaluated through histological analysis. CONCLUSIONS: This study developed an MR-compatible extracorporeal perfusion technique that effectively maintains tissue viability, allowing for the direct measurement of patient-specific MR, thermal, and acoustic property values for both fibroid and myometrium tissues. These measured tissue property values will enable further development and validation of treatment planning models that can be utilized during MRgFUS uterine fibroid treatments.

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.

Experimental assessment of phase aberration correction for breast MRgFUS therapy.

PURPOSE: This study validates that phase aberrations in breast magnetic resonance-guided focussed ultrasound (MRgFUS) therapies can be corrected in a clinically relevant time frame to generate more intense, smaller and more spatially accurate foci. MATERIALS AND METHODS: Hybrid angular spectrum (HAS) ultrasound calculations in an magnetic resonance imaging (MRI)-based tissue model, were used to compute phase aberration corrections for improved experimental MRgFUS heating in four heterogeneous breast-mimicking phantoms (n = 18 total locations). Magnetic resonance(MR) temperature imaging was used to evaluate the maximum temperature rise, focus volume and focus accuracy for uncorrected and phase aberration-corrected sonications. Thermal simulations assessed the effectiveness of the phase aberration correction implementation. RESULTS: In 13 of 18 locations, the maximum temperature rise increased by an average of 30%, focus volume was reduced by 40% and focus accuracy improved from 4.6 to 3.6 mm. Mixed results were observed in five of the 18 locations, with focus accuracy improving from 6.1 to 2.5 mm and the maximum temperature rise decreasing by 8% and focus volume increasing by 10%. Overall, the study demonstrated significant improvements (p < 0.005) in maximum temperature rise, focus volume and focus accuracy. Simulations predicted greater improvements than observed experimentally, suggesting potential for improvement in implementing the technique. The complete phase aberration correction procedure, including model generation, segmentation and phase aberration computations, required less than 45 min per sonication location. CONCLUSION: The significant improvements demonstrated in this study i.e., focus intensity, size and accuracy from phase aberration correction have the potential to improve the efficacy, time-efficiency and safety of breast MRgFUS therapies.

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

PURPOSE: This study investigates the feasibility of non-invasively determining thermal diffusivity (α) and the Pennes perfusion parameter (w) from pre-clinical and clinical magnetic resonance-guided focussed ultrasound (MRgFUS) temperature data. MATERIALS AND METHODS: Pre-clinical MRgFUS experiments were performed in rabbit muscle (N = 3, 28 sonications) using three-dimensional MR thermometry. Eight sonications were made in a clinical QA phantom with two-dimensional thermometry. Retrospective property determination was performed on clinical uterine fibroid (N = 8, 9 sonications) and desmoid tumour (N = 4, 7 sonications) data. The property determination method fits an analytical solution to MRgFUS temperatures in the coronal MR plane, including all temperatures acquired during heating and one cooling image. When possible, additional cooling data were acquired for property determination. RESULTS: Rabbit α and w from Heating Data (α = 0.164 mm2s-1, w = 7.9 kg m-3 s-1) and Heating and Cooling Data (α = 0.146 mm2s-1, w = 3.3 kg m-3 s-1) were within the range of gold-standard invasive measurements, with >50% reduction in variability by including cooling data. QA phantom property determination with cooling data yielded properties within 3% of expected values (α = 0.144 mm2s-1, w = 0.0 kg m-3 s-1), a difference that was not statistically significant (p = 0.053). Uterine fibroid (Heating Data: α = 0.212 mm2s-1, w = 11.0 kg m-3 s-1) and desmoid tumour (Heating & Cooling Data: α = 0.245 mm2s-1, w = 4.7 kg m-3 s-1) properties are feasible but lack independent verification. CONCLUSIONS: Thermal diffusivity and the Pennes perfusion parameter can be obtained from in vivo data and with clinical MRgFUS protocols. Property values are consistently improved by including cooling data. The utility of this property determination method will increase as clinical protocols implement improved temperature imaging.