Adapting MRI acoustic radiation force imaging for in vivo human brain focused ultrasound applications

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

OBJECTIVE

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

MATERIALS AND METHODS

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

RESULTS

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

DISCUSSION

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

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Reduced-field of view three-dimensional MR acoustic radiation force imaging with a low-rank reconstruction for targeting transcranial focused ultrasound

PURPOSE

To rapidly image and localize the focus in MR-guided focused ultrasound (FUS) while maintaining a low ultrasound duty cycle to minimize tissue effects.

METHODS

MR-acoustic radiation force imaging (ARFI) is key to targeting FUS procedures such as neuromodulation, and works by encoding ultrasound-induced displacements into the phase of MR images. However, it can require long scan times to cover a volume of tissue, especially when minimizing the FUS dose during targeting is paramount. To simultaneously minimize scan time and the FUS duty cycle, a 2-min three-dimensional (3D) reduced-FOV spin echo ARFI scan with two-dimensional undersampling was implemented at 3T with a FUS duty cycle of 0.85%. The 3D k-space sampling scheme incorporated uniform undersampling in one phase-encoded axis and partial Fourier (PF) sampling in the other. The scan interleaved FUS-on and FUS-off data collection to improve displacement map quality via a joint low-rank image reconstruction. Experiments in agarose and graphite phantoms and living macaque brains for neuromodulation and blood-brain barrier opening studied the effects of the sampling and reconstruction strategy on the acquisition, and evaluated its repeatability and accuracy.

RESULTS

In the phantom, the distances between displacement centroids of 10 prospective reconstructions and a fully sampled reference were below 1 mm. In in vivo brain, the distances between centroids ranged from 1.3 to 2.1 mm. Results in phantom and in vivo brain both showed that the proposed method can recover the FUS focus compared to slower fully sampled scans.

CONCLUSION

The proposed 3D MR-ARFI reduced-FOV method enables rapid imaging of the FUS focus while maintaining a low FUS duty cycle.

Improving in situ acoustic intensity estimates using MR acoustic radiation force imaging in combination with multifrequency MR elastography

PURPOSE

Magnetic resonance acoustic radiation force imaging (MR-ARFI) enables focal spot localization during nonablative transcranial ultrasound therapies. As the acoustic radiation force is proportional to the applied acoustic intensity, measured MR-ARFI displacements could potentially be used to estimate the acoustic intensity at the target. However, variable brain stiffness is an obstacle. The goal of this study was to develop and assess a method to accurately estimate the acoustic intensity at the focus using MR-ARFI displacements in combination with viscoelastic properties obtained with multifrequency MR elastography (MRE).

METHODS

Phantoms with a range of viscoelastic properties were fabricated, and MR-ARFI displacements were acquired within each phantom using multiple acoustic intensities. Voigt model parameters were estimated for each phantom based on storage and loss moduli measured using multifrequency MRE, and these were used to predict the relationship between acoustic intensity and measured displacement.

RESULTS

Using assumed viscoelastic properties, MR-ARFI displacements alone could not accurately estimate acoustic intensity across phantoms. For example, acoustic intensities were underestimated in phantoms stiffer than the assumed stiffness and overestimated in phantoms softer than the assumed stiffness. This error was greatly reduced using individualized viscoelasticity measurements obtained from MRE.

CONCLUSION

We demonstrated that viscoelasticity information from MRE could be used in combination with MR-ARFI displacements to obtain more accurate estimates of acoustic intensity. Additionally, Voigt model viscosity parameters were found to be predictive of the relaxation rate of each phantom's time-varying displacement response, which could be used to optimize patient-specific MR-ARFI pulse sequences.

Transcranial Magnetic Resonance-Guided Histotripsy for Brain Surgery: Pre-clinical Investigation

Histotripsy has been previously applied to target various cranial locations in vitro through an excised human skull. Recently, a transcranial magnetic resonance (MR)-guided histotripsy (tcMRgHt) system was developed, enabling pre-clinical investigations of tcMRgHt for brain surgery. To determine the feasibility of in vivo transcranial histotripsy, tcMRgHt treatment was delivered to eight pigs using a 700-kHz, 128-element, MR-compatible phased-array transducer inside a 3-T magnetic resonance imaging (MRI) scanner. After craniotomy to open an acoustic window to the brain, histotripsy was applied through an excised human calvarium to target the inside of the pig brain based on pre-treatment MRI and fiducial markers. MR images were acquired pre-treatment, immediately post-treatment and 2-4 h post-treatment to evaluate the acute treatment outcome. Successful histotripsy ablation was observed in all pigs. The MR-evident lesions were well confined within the targeted volume, without evidence of excessive brain edema or hemorrhage outside of the target zone. Histology revealed tissue homogenization in the ablation zones with a sharp demarcation between destroyed and unaffected tissue, which correlated well with the radiographic treatment zones on MRI. These results are the first to support the in vivo feasibility of tcMRgHt in the pig brain, enabling further investigation of the use of tcMRgHt for brain surgery.

Techniques, Indications, and Outcomes in Magnetic Resonance-guided Focused Ultrasound Thalamotomy for Tremor

Magnetic resonance (MR)-guided focused ultrasound surgery (MRgFUS) is the latest minimally invasive stereotactic procedure, and thalamotomy using this novel modality has demonstrated its effectiveness and safety, especially for patients with essential tremor (ET) and Parkinson's disease (PD). In Japan, the application of MRgFUS to treat ET and PD has recently been covered by health insurance. Technically, the transducer with 1024 elements emits ultrasound beams, which are then focused on the target with a phase control, resulting in optimal ablation by thermal coagulation. The technical advantages of MRgFUS are continuous intraoperative monitoring of clinical symptoms and MR images and fine adjustment of the target by the steering function. Postoperative tremor control is compatible with other modalities, although long-term follow-up is necessary. The adverse effects are usually transient and acceptable. Prognostic factors for good tremor control include high temperature and large lesion size. A high skull density ratio is a factor to achieve high temperature and large lesioning, but it may not be necessary and sufficient for clinical outcomes. For patients with advanced symptoms such as bilateral tremor or head/neck tremor, deep brain stimulation may be recommended because of the adjustability of stimulation and the possibility of bilateral treatment. Patients have high expectations of MRgFUS because of its non-invasiveness. To perform this treatment safely and effectively, physicians need to understand the technological aspects, the physiological principles. To choose the appropriate modality, physicians also should recognize the clinical advantages and disadvantages of MRgFUS compared to other modalities.

Simultaneous acoustic radiation force imaging and MR thermometry based on a coherent echo-shifted sequence

Background

Simultaneous magnetic resonance (MR) acoustic radiation force imaging (ARFI) and MR thermometry (MRT) (STARFI) based on coherent echo-shifted (cES) sequence was proposed and comprehensively compared to radiofrequency (RF)-spoiled gradient echo (spGRE) STARFI.

Methods

Through use of delicately designed gradients, a collection of echoes was delayed by one repetition time (TR) cycle. The crusher gradient after readout (RO) was used as the displacement encoding gradient (DEG). The sequence was intrinsically sensitive to temperature. High-intensity focused ultrasound (HIFU) pulses were interleaved ON/OFF in successive TRs to separate the phase changes induced by displacement due to acoustic radiation force (ARF) impulses and temperature. Bloch simulation was performed to study the phase sensitivity to displacement of the proposed cES STARFI and spGRE STARFI. The proposed cES sequence was evaluated and compared to spGRE STARFI in ex vivo porcine muscle and ex vivo porcine brain.

Results

The minimally achievable TR of cES STARFI was shorter than that of spGRE STARFI, indicating that the cES sequence was more time efficient. It was verified through Bloch simulation and ex vivo experiments that the phase sensitivity to displacement of cES STARFI was higher than that of spGRE STARFI. The optimal trigger delays of cES STARFI and spGRE STARFI in ex vivo porcine muscle were t_offset =-2 and -1 ms, respectively. The displacement-induced phase change to acoustic pressure slopes of cES STARFI were 0.079, 0.079, and 0.047 rad/Mpa across the three muscle samples, while the slopes of spGRE STARFI were only 0.047, 0.052, and 0.027 rad/Mpa. The maximum temperature difference between cES STARFI and spGRE STARFI was 1.1 °C. In ex vivo porcine brain, both the displacement-induced phase-to-noise ratio (PNRd) and the temperature uncertainty of cES STARFI were better than those of spGRE STARFI (P<0.05). The temperature and displacement-induced phase change maps of cES STARFI and spGRE STARFI during HIFU treatment were in good accordance in time and spatial location.

Conclusions

The cES STARFI sequence can provide simultaneous MR-ARFI and temperature measurements during pulsed HIFU applications. Though the exact displacement cannot be quantified directly, the sequence showed increased phase sensitivity compared with the spGRE sequence and provided efficient visualization of the focal spot. cES STARFI could therefore be a desirable alternative to spGRE STARFI in practical applications.

Focused Ultrasound for Noninvasive, Focal Pharmacologic Neurointervention

A long-standing goal of translational neuroscience is the ability to noninvasively deliver therapeutic agents to specific brain regions with high spatiotemporal resolution. Focused ultrasound (FUS) is an emerging technology that can noninvasively deliver energy up the order of 1 kW/cm² with millimeter and millisecond resolution to any point in the human brain with Food and Drug Administration-approved hardware. Although FUS is clinically utilized primarily for focal ablation in conditions such as essential tremor, recent breakthroughs have enabled the use of FUS for drug delivery at lower intensities (i.e., tens of watts per square centimeter) without ablation of the tissue. In this review, we present strategies for image-guided FUS-mediated pharmacologic neurointerventions. First, we discuss blood–brain barrier opening to deliver therapeutic agents of a variety of sizes to the central nervous system. We then describe the use of ultrasound-sensitive nanoparticles to noninvasively deliver small molecules to millimeter-sized structures including superficial cortical regions and deep gray matter regions within the brain without the need for blood–brain barrier opening. We also consider the safety and potential complications of these techniques, with attention to temporal acuity. Finally, we close with a discussion of different methods for mapping the ultrasound field within the brain and describe future avenues of research in ultrasound-targeted drug therapies.

Heterogeneous Angular Spectrum Method for Trans-Skull Imaging and Focusing

Ultrasound, alone or in concert with circulating microbubble contrast agents, has emerged as a promising modality for therapy and imaging of brain diseases. While this has become possible due to advancements in aberration correction methods, a range of applications, including adaptive focusing and tracking of the microbubble dynamics through the human skull, may benefit from even more computationally efficient methods to account for skull aberrations. Here, we derive a general method for the angular spectrum approach (ASA) in a heterogeneous medium, based on a numerical marching scheme to approximate the full implicit solution. We then demonstrate its functionality with simulations for (human) skull-related aberration correction and trans-skull passive acoustic mapping. Our simulations show that the general solution provides accurate trans-skull focusing as compared to the uncorrected case (error in focal point location of 1.0 ± 0.4 mm vs 2.2 ± 0.7 mm) for clinically relevant frequencies (0.25-1.5MHz), apertures (50-100 mm), and targets, with peak focal pressures approximately 30 ± 17% of the free field case, with the effects of skull attenuation and amplitude shading included. In the case of source localization, our method leads to an average of 75% error reduction (from 2.9 ± 1.8 mm to 0.7 ± 0.5 mm) and 40-60% increase in peak intensity, evaluated over the range of frequencies (0.4-1.2 MHz), apertures (50-100 mm), and point source locations (40 mm by 50 mm grid) as compared to the homogeneous medium ASA. Overall, total computation times for both focusing and point source localization of the order milliseconds (166 ± 37 ms, compared with 44 ± 4 ms for the homogeneous ASA formulation) can be attained with this approach. Collectively our findings indicate that the proposed phase correction method based on the ASA could provide a computationally efficient and accurate method for trans-skull transmit focusing and imaging of point scatterers, potentially opening new possibilities for treatment and diagnosis of brain diseases.

Ultrasound Neuromodulation: A Review of Results, Mechanisms and Safety

Ultrasonic neuromodulation is a rapidly growing field, in which low-intensity ultrasound (US) is delivered to nervous system tissue, resulting in transient modulation of neural activity. This review summarizes the findings in the central and peripheral nervous systems from mechanistic studies in cell culture to cognitive behavioral studies in humans. The mechanisms by which US mechanically interacts with neurons and could affect firing are presented. An in-depth safety assessment of current studies shows that parameters for the human studies fall within the safety envelope for US imaging. Challenges associated with accurately targeting US and monitoring the response are described. In conclusion, the literature supports the use of US as a safe, non-invasive brain stimulation modality with improved spatial localization and depth targeting compared with alternative methods. US neurostimulation has the potential to be used both as a scientific instrument to investigate brain function and as a therapeutic modality to modulate brain activity.

Production of acoustic radiation force using ultrasound: methods and applications

INTRODUCTION

Acoustic radiation force (ARF) is used in many biomedical applications. The transfer of momentum in acoustic waves can be used in a multitude of ways to perturb tissue and manipulate cells.

AREAS COVERED

This review will briefly cover the acoustic theory related to ARF, particularly that related to application in tissues. The use of ARF in measurement of mechanical properties will be treated in detail with emphasis on the spatial and temporal modulation of the ARF. Additional topics covered will be the manipulation of particles with ARF, correction of phase aberration with ARF, modulation of cellular behavior with ARF, and bioeffects related to ARF use.

EXPERT COMMENTARY

The use of ARF can be tailored to specific applications for measurements of mechanical properties or correction of focusing for ultrasound beams. Additionally, noncontact manipulation of particles and cells with ARF enables a wide array of applications for tissue engineering and biosensing.

Design and Implementation of a Transmit/Receive Ultrasound Phased Array for Brain Applications

Focused ultrasound phased array systems have attracted increased attention for brain therapy applications. However, such systems currently lack a direct and real-time method to intraoperatively monitor ultrasound pressure distribution for securing treatment. This study proposes a dual-mode ultrasound phased array system design to support transmit/receive operations for concurrent ultrasound exposure and backscattered focal beam reconstruction through a spherically focused ultrasound array. A 256-channel ultrasound transmission system was used to transmit focused ultrasonic energy (full 256 channels), with an extended implementation of multiple-channel receiving function (up to 64 channels) using the same 256-channel ultrasound array. A coherent backscatter-received beam formation algorithm was implemented to map the point spread function (PSF) and focal beam distribution under a free-field/transcranial environment setup, with the backscattering generated from a strong scatterer (a point reflector or a microbubble-perfused tube) or a weakly scattered tissue-mimicking graphite phantom. Our results showed that PSF and focal beam can be successfully reconstructed and visualized in free-field conditions and can also be transcranially reconstructed following skull-induced aberration correction. In vivo experiments were conducted to demonstrate its capability to preoperatively and semiquantitatively map a focal beam to guide blood-brain barrier opening. The proposed system may have potential for real-time guidance of ultrasound brain intervention, and may facilitate the design of a dual-mode ultrasound phased array for brain therapeutic applications.

High sensitivity MR acoustic radiation force imaging using transition band balanced steady-state free precession

PURPOSE

MR acoustic radiation force imaging (MR-ARFI) provides a method to visualize the focal spot of a focused ultrasound (FUS) beam without introducing a significant temperature rise. With conventional spoiled MR-ARFI pulse sequences, the ARFI phase always equals the motion-encoded phase. In this work, MR-ARFI using transition band balanced steady-state free precession (bSSFP) is presented, which improves the sensitivity of MR-ARFI with high acquisition speed.

THEORY AND METHODS

Motion-encoding gradients (MEG) are inserted into bSSFP sequences for MR-ARFI. By applying an ultrasound pulse during the MEG, motion-encoded phase is generated, which leads to an amplified change in the image phase when operating in the bSSFP transition band. MR-ARFI was performed on a homemade gel phantom using both the proposed technique and a spoiled gradient echo ARFI sequence with identical MEG and FUS, and ARFI images were compared.

RESULTS

The bSSFP-ARFI sequence generated an ARFI image phase that is more than 5 times larger than the motion-encoded phase in a few seconds with 2DFT readout. By keeping FUS pulses as short as 1.45 ms, temperature rise was insignificant during the measurement.

CONCLUSION

bSSFP-ARFI has enhanced sensitivity compared with conventional MR-ARFI pulse sequences and could provide an efficient way to visualize the focal spot. Magn Reson Med 79:1532-1537, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Optimization of acoustic radiation force imaging: Influence of timing parameters on sensitivity

PURPOSE

Optimization of timing parameters for MR-guided ARFI to achieve the highest displacement signal-to-noise ratio (SNR_d ).

THEORY AND METHODS

In MR-guided ARFI the phase change induced by motion encoding gradients (MEGs) is measured to assess tissue displacement. The sensitivity of this encoding procedure depends on several timing parameters, such as the MEG duration and the offset time between ultrasound (US) and MEG. Furthermore, mechanical and MR tissue constants and MEG schemes (bipolar or three-lobed) influence SNR_d . Optimal timing parameters were determined in simulations for bipolar and three-lobed MEGs, and the results were compared with measurements. To provide clinically usable timing parameters, physiologically relevant ranges of tissue constants were considered.

RESULTS

For the considered ranges of tissue constants, optimal timing parameters provide only 6% higher SNR_d for bipolar than for three-lobed MEG. Three-lobed MEG is less sensitive to motion as confirmed in phantom experiments. Bipolar MEG can use approximately 1.5-fold shorter MEG durations.

CONCLUSION

Both bipolar and three-lobed MEGs can yield approximately the same SNR_d if the optimal timing parameters are chosen. Bipolar MEG allows for shorter durations, which is preferable if deposition of US energy needs to be minimized, and three-lobed MEG is more suitable when residual motion compensation is necessary. Magn Reson Med 79:981-986, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

MR techniques for guiding high-intensity focused ultrasound (HIFU) treatments

To make full use of the ability of magnetic resonance (MR) to guide high-intensity focused ultrasound (HIFU) treatment, effort has been made to improve techniques for thermometry, motion tracking, and sound beam visualization. For monitoring rapid temperature elevation with proton resonance frequency (PRF) shift, data acquisition and processing can be accelerated with parallel imaging and/or sparse sampling in conjunction with appropriate signal processing methods. Thermometry should be robust against tissue motion, motion-induced magnetic field variation, and susceptibility change. Thus, multibaseline, referenceless, or hybrid techniques have become important. In cases with adipose or bony tissues, for which PRF shift cannot be used, thermometry with relaxation times or signal intensity may be utilized. Motion tracking is crucial not only for thermometry but also for targeting the focus of an ultrasound in moving organs such as the liver, kidney, or heart. Various techniques for motion tracking, such as those based on an anatomical image atlas with optical-flow displacement detection, a navigator echo to seize the diaphragm position, and/or rapid imaging to track vessel positions, have been proposed. Techniques for avoiding the ribcage and near-field heating have also been examined. MR acoustic radiation force imaging (MR-ARFI) is an alternative to thermometry that can identify the location and shape of the focal spot and sound beam path. This technique could be useful for treating heterogeneous tissue regions or performing transcranial therapy. All of these developments, which will be discussed further in this review, expand the applicability of HIFU treatments to a variety of clinical targets while maintaining safety and precision.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 4 J. Magn. Reson. Imaging 2018;47:316-331.

MR thermometry for focused ultrasound monitoring utilizing model predictive filtering and ultrasound beam modeling

Background

A major challenge in using magnetic resonance temperature imaging (MRTI) to monitor focused ultrasound (FUS) applications is achieving high spatio-temporal resolution over a large field of view (FOV). This is important to accurately monitor all ultrasound (US) power depositions. Magnetic resonance (MR) subsampling in conjunction with thermal model-based reconstruction of the MRTI utilizing Pennes bioheat transfer equation (PBTE) is one promising approach. The thermal properties used in the thermal model are often estimated from a pre-treatment, low-power sonication.

Methods

In this proof-of-concept study we investigate the use of US simulations computed using the hybrid angular spectrum (HAS) method to estimate the US power deposition density Q, thereby avoiding the pre-treatment sonication and any potential tissue damage. MRTI reconstructions are performed using a thermal model-based reconstruction method called model predictive filtering (MPF). Experiments are performed in a homogeneous gelatin phantom and in a gelatin phantom with embedded plastic skull. MPF reconstructions are compared to separate sonications imaged with fully sampled data over a smaller FOV. Temperature root-mean-square errors (RMSE) and focal spot positions and shapes are evaluated.

Results

HAS simulations accurately predict the location of the focal spot (to within 1 mm) in both phantoms. Accurate temperature maps (RMSE below 1 °C), where the location of the focal spot agrees well with fully sampled “truth” (to within 1 mm), are also achieved in both phantoms.

Conclusions

HAS simulations can be used to accurately predict the focal spot location in homogeneous media and when focusing through an aberrating plastic skull. The HAS simulated power deposition (Q) patterns can be used in the MPF thermal model-based reconstruction to obtain accurate temperature maps with high spatio-temporal resolution over large FOVs.

Variable ultrasound trigger delay for improved magnetic resonance acoustic radiation force imaging

Magnetic resonance acoustic radiation force imaging (MR-ARFI) allows the quantification of microscopic displacements induced by ultrasound pulses, which are proportional to the local acoustic intensity. This study describes a new method to acquire MR-ARFI maps, which reduces the measurement noise in the quantification of displacement as well as improving its robustness in the presence of motion. Two MR-ARFI sequences were compared in this study. The first sequence 'variable MSG' involves switching the polarity of the motion sensitive gradient (MSG) between odd and even image frames. The second sequence named 'static MSG' involves a variable ultrasound trigger delay to sonicate during the first or second MSG for odd and even image frames, respectively. As previously published, the data acquired with a variable MSG required the use of reference data acquired prior to any sonication to process displacement maps. In contrary, data acquired with a static MSG were converted to displacement maps without using reference data acquired prior to the sonication. Displacement maps acquired with both sequences were compared by performing sonications for three different conditions: in a polyacrylamide phantom, in the leg muscle of a freely breathing pig and in the leg muscle of pig under apnea. The comparison of images acquired at even image frames and odd image frames indicates that the sequence with a static MSG provides a significantly better steady state (p < 0.001 based on a Student's t-test) than the images acquired with a variable MSG. In addition no reference data prior to sonication were required to process displacement maps for data acquired with a static MSG. The absence of reference data prior to sonication provided a 41% reduction of the spatial distribution of noise (p < 0.001 based on a Student's t-test) and reduced the sensitivity to motion for displacements acquired with a static MSG. No significant differences were expected and observed for thermal maps acquired with a variable MSG and a static MSG. The use of a static MSG with a variable ultrasound trigger delay improves the ARFI displacement map quality without additional acquisition time and remains compatible with the simultaneous acquisition of MR thermal maps.

Evaluation of a three-dimensional MR acoustic radiation force imaging pulse sequence using a novel unbalanced bipolar motion encoding gradient

PURPOSE

MR guided focused ultrasound procedures require accurate focal spot localization in three dimensions. This study presents a three-dimensional (3D) pulse sequence for acoustic radiation force imaging (ARFI) that efficiently localizes the focal spot by means of ultrasound induced tissue displacement over a large field-of-view.

METHODS

A novel unbalanced bipolar motion encoding gradient was implemented to maximize time available for motion encoding, reduce echo times, and allow for longer echo train lengths. Two advanced features, kz reduction factor (KZRF) and kz -level interleaving, were implemented to reduce tissue heating. Studies in gelatin phantoms compared the location of peak displacement and temperature measured by 3D MR thermometry. MR-ARFI induced tissue heating was evaluated through a parametric study of sequence parameters and MR thermometry measurements during repeated application of ARFI sonication patterns. Sequence performance was characterized in the presence of respiration and tissue inhomogeneity.

RESULTS

The location of peak displacement and temperature rise agreed within 0.2 ± 0.1 mm and 0.5 ± 0.3 mm in the transverse and longitudinal direction, respectively. The 3D displacement maps were acquired safely, and the KZRF and kz -level interleaving features reduced tissue heating by 51%. High quality displacement maps were obtained despite respiration and tissue inhomogeneities.

CONCLUSION

This sequence provides a safe, accurate, and simple approach to localizing the focal spot in three dimensions with a single scan. Magn Reson Med 76:803-813, 2016. © 2015 Wiley Periodicals, Inc.

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.

Transcranial MRI-Guided Focused Ultrasound: A Review of the Technologic and Neurologic Applications

OBJECTIVE

This article reviews the physical principles of MRI-guided focused ultra-sound and discusses current and potential applications of this exciting technology.

CONCLUSION

MRI-guided focused ultrasound is a new minimally invasive method of targeted tissue thermal ablation that may be of use to treat central neuropathic pain, essential tremor, Parkinson tremor, and brain tumors. The system has also been used to temporarily disrupt the blood-brain barrier to allow targeted drug delivery to brain tumors.

Transcranial phase aberration correction using beam simulations and MR-ARFI

PURPOSE

Transcranial magnetic resonance-guided focused ultrasound surgery is a noninvasive technique for causing selective tissue necrosis. Variations in density, thickness, and shape of the skull cause aberrations in the location and shape of the focal zone. In this paper, the authors propose a hybrid simulation-MR-ARFI technique to achieve aberration correction for transcranial MR-guided focused ultrasound surgery. The technique uses ultrasound beam propagation simulations with MR Acoustic Radiation Force Imaging (MR-ARFI) to correct skull-caused phase aberrations.

METHODS

Skull-based numerical aberrations were obtained from a MR-guided focused ultrasound patient treatment and were added to all elements of the InSightec conformal bone focused ultrasound surgery transducer during transmission. In the first experiment, the 1024 aberrations derived from a human skull were condensed into 16 aberrations by averaging over the transducer area of 64 elements. In the second experiment, all 1024 aberrations were applied to the transducer. The aberrated MR-ARFI images were used in the hybrid simulation-MR-ARFI technique to find 16 estimated aberrations. These estimated aberrations were subtracted from the original aberrations to result in the corrected images. Each aberration experiment (16-aberration and 1024-aberration) was repeated three times.

RESULTS

The corrected MR-ARFI image was compared to the aberrated image and the ideal image (image with zero aberrations) for each experiment. The hybrid simulation-MR-ARFI technique resulted in an average increase in focal MR-ARFI phase of 44% for the 16-aberration case and 52% for the 1024-aberration case, and recovered 83% and 39% of the ideal MR-ARFI phase for the 16-aberrations and 1024-aberration case, respectively.

CONCLUSIONS

Using one MR-ARFI image and noa priori information about the applied phase aberrations, the hybrid simulation-MR-ARFI technique improved the maximum MR-ARFI phase of the beam's focus.

MR-acoustic radiation force imaging (MR-ARFI) and susceptibility weighted imaging (SWI) to visualize calcifications in ex vivo swine brain

PURPOSE

To present the use of MR-acoustic radiation force imaging (MR-ARFI) and susceptibility weighted imaging (SWI) to visualize calcifications in ex vivo brain tissue as a planning indicator for MR-guided focused ultrasound (MRgFUS).

MATERIALS AND METHODS

Calcifications were implanted in ex vivo swine brain and imaged using SWI, MR-ARFI, and computed tomography (CT). SWI-filtered phase images used 3D gradient recalled echo (GRE) images with a Fourier-based unwrapping algorithm. The MR-ARFI pulse sequence used a 2DFT spin-echo with repeated bipolar encoding gradients in the direction of the longitudinal ultrasound beam. MR-ARFI interrogations scanned a subregion (14 × 10 × 12 mm) of the brain surrounding the calcification. They were combined into a single displacement weighted map, using the sum of squares method. Calcification size estimates were based on image profiles plotted along the ±x and ±z direction, at the full-width half-maximum.

RESULTS

Both MR-ARFI and SWI were able to visualize the calcifications. The contrast ratio was 150 for CT, 12 for SWI, and 12 for MR-ARFI. Profile measures were 1.35 × 1.28 mm on CT, 1.24 × 1.73 mm on SWI, and 2.45 × 3.02 mm on MR-ARFI. MR-ARFI displacement showed a linear increase with acoustic power (20-80 W), and also increased with calcification size.

CONCLUSION

The use of SWI-filtered phase and MR-ARFI have the potential to provide a clinical indicator of calcification relevance in the planning of a transcranial MRgFUS treatment.

Application of Zernike polynomials towards accelerated adaptive focusing of transcranial high intensity focused ultrasound

PURPOSE

To study the phase aberrations produced by human skulls during transcranial magnetic resonance imaging guided focused ultrasound surgery (MRgFUS), to demonstrate the potential of Zernike polynomials (ZPs) to accelerate the adaptive focusing process, and to investigate the benefits of using phase corrections obtained in previous studies to provide the initial guess for correction of a new data set.

METHODS

The five phase aberration data sets, analyzed here, were calculated based on preoperative computerized tomography (CT) images of the head obtained during previous transcranial MRgFUS treatments performed using a clinical prototype hemispherical transducer. The noniterative adaptive focusing algorithm [Larrat et al., "MR-guided adaptive focusing of ultrasound," IEEE Trans. Ultrason. Ferroelectr. Freq. Control 57(8), 1734-1747 (2010)] was modified by replacing Hadamard encoding with Zernike encoding. The algorithm was tested in simulations to correct the patients' phase aberrations. MR acoustic radiation force imaging (MR-ARFI) was used to visualize the effect of the phase aberration correction on the focusing of a hemispherical transducer. In addition, two methods for constructing initial phase correction estimate based on previous patient's data were investigated. The benefits of the initial estimates in the Zernike-based algorithm were analyzed by measuring their effect on the ultrasound intensity at the focus and on the number of ZP modes necessary to achieve 90% of the intensity of the nonaberrated case.

RESULTS

Covariance of the pairs of the phase aberrations data sets showed high correlation between aberration data of several patients and suggested that subgroups can be based on level of correlation. Simulation of the Zernike-based algorithm demonstrated the overall greater correction effectiveness of the low modes of ZPs. The focal intensity achieves 90% of nonaberrated intensity using fewer than 170 modes of ZPs. The initial estimates based on using the average of the phase aberration data from the individual subgroups of subjects was shown to increase the intensity at the focal spot for the five subjects.

CONCLUSIONS

The application of ZPs to phase aberration correction was shown to be beneficial for adaptive focusing of transcranial ultrasound. The skull-based phase aberrations were found to be well approximated by the number of ZP modes representing only a fraction of the number of elements in the hemispherical transducer. Implementing the initial phase aberration estimate together with Zernike-based algorithm can be used to improve the robustness and can potentially greatly increase the viability of MR-ARFI-based focusing for a clinical transcranial MRgFUS therapy.