Simulations and measurements of transcranial low-frequency ultrasound therapy: skull-base heating and effective area of treatment

Strategies and safety simulations for ultrasonic cervical spinal cord neuromodulation

Objective. Focused ultrasound spinal cord neuromodulation has been demonstrated in small animals. However, most of the tested neuromodulatory exposures are similar in intensity and exposure duration to the reported small animal threshold for possible spinal cord damage. All efforts must be made to minimize the risk and assure the safety of potential human studies, while maximizing potential treatment efficacy. This requires an understanding of ultrasound propagation and heat deposition within the human spine.Approach. Combined acoustic and thermal modelling was used to assess the pressure and heat distributions produced by a 500 kHz source focused to the C5/C6 level via two approaches (a) the posterior acoustic window between vertebral posterior arches, and (b) the lateral intervertebral foramen from which the C6 spinal nerve exits. Pulse trains of fifty 0.1 s pulses (pulse repetition frequency: 0.33 Hz, free-field spatial peak pulse-averaged intensity: 10 W cm-2) were simulated for four subjects and for ±10 mm translational and ±10∘rotational source positioning errors.Main results.Target pressures ranged between 20%-70% of free-field spatial peak pressures with the posterior approach, and 20%-100% with the lateral approach. When the posterior source was optimally positioned, peak spine heating values were below 1 ∘C, but source mispositioning resulted in bone heating up to 4 ∘C. Heating with the lateral approach did not exceed 2 ∘C within the mispositioning range. There were substantial inter-subject differences in target pressures and peak heating values. Target pressure varied three to four-fold between subjects, depending on approach, while peak heating varied approximately two-fold between subjects. This results in a nearly ten-fold range between subjects in the target pressure achieved per degree of maximum heating.Significance. This study highlights the utility of trans-spine ultrasound simulation software and need for precise source-anatomy positioning to assure the subject-specific safety and efficacy of focused ultrasound spinal cord therapies.

Noninvasive modulation of essential tremor with focused ultrasonic waves

Objective: Transcranial focused low-intensity ultrasound has the potential to noninvasively modulate confined regions deep inside the human brain, which could provide a new tool for causal interrogation of circuit function in humans. However, it has been unclear whether the approach is potent enough to modulate behavior.Approach: To test this, we applied low-intensity ultrasound to a deep brain thalamic target, the ventral intermediate nucleus, in three patients with essential tremor.Main results: Brief, 15 s stimulations of the target at 10% duty cycle with low-intensity ultrasound, repeated less than 30 times over a period of 90 min, nearly abolished tremor (98% and 97% tremor amplitude reduction) in 2 out of 3 patients. The effect was observed within seconds of the stimulation onset and increased with ultrasound exposure time. The effect gradually vanished following the stimulation, suggesting that the stimulation was safe with no harmful long-term consequences detected.Significance: This result demonstrates that low-intensity focused ultrasound can robustly modulate deep brain regions in humans with notable effects on overt motor behavior.

An efficient method for transcranial ultrasound focus correction based on the coupling of boundary integrals and finite elements

Transcranial focused ultrasound is a novel technique for the noninvasive treatment of brain diseases. The success of the treatment greatly depends on achieving precise and efficient intraoperative focus. However, compensating for aberrated ultrasound waves caused by the skull through numerical simulation-based phase corrections is a challenging task due to the significant computational burden involved in solving the acoustic wave equation. In this article, we propose a promising strategy using the coupling of the boundary integral equation method (BIEM) and the finite element method (FEM) to overcome the above limitation. Specifically, we adopt the BIEM to obtain the Robin-to-Dirichlet maps on the boundaries of the skull and then couple the maps to the FEM matrices via a dual interpolation technique, resulting in a computational domain including only the skull. Three simulation experiments were conducted to evaluate the effectiveness of the proposed method, including a convergence test and two skull-induced aberration corrections in 2D and 3D ultrasound. The results show that the method's convergence is guaranteed as the element size decreases, leading to a decrease in pressure error. The computation times for simulating a 500 kHz ultrasound field on a regular desktop computer were found to be 0.47 ± 0.01 s in the 2D case and 43.72 ± 1.49 s in the 3D case, provided that lower-upper decomposition (approximately 13 s in 2D and 2.5 h in 3D) was implemented in advance. We also demonstrated that more accurate transcranial focusing can be achieved by phase correction compared to the noncorrected results (with errors of 1.02 mm vs. 6.45 mm in 2D and 0.28 mm vs. 3.07 mm in 3D). The proposed strategy is valuable for enabling online ultrasound simulations during treatment, facilitating real-time adjustments and interventions.

A phase I clinical trial of sonodynamic therapy combined with temozolomide in the treatment of recurrent glioblastoma

PURPOSE

The prognosis of recurrent glioblastoma (rGBM) is poor, and there is currently no effective treatment strategy. Sonodynamic therapy (SDT) is a new method for cancer treatment that uses a combination of low-frequency ultrasound and sonosensitisers to produce antitumor effects, which have shown good therapeutic effects in preclinical studies. Therefore, we initiated an open, prospective pilot study to evaluate the safety, tolerability, and efficacy of SDT for the treatment of rGBM.

METHODS

Nine patients with rGBM were enrolled who had received multiple treatments, but the nidus continued to progress without additional standard treatments. After MRI localisation, porphyrin drugs were injected, and intermittent low-frequency ultrasound therapy was performed for five days.

RESULTS

None of the nine patients in this clinical trial showed any clinical, neurological, haematological, or skin-targeted adverse effects associated with SDT. After the completion of the trial, one patient maintained stable disease, and eight patients experienced disease progression. Among the eight with progressive disease, the median progression-free survival time was 84 days. Four patients died, and the median overall survival duration after recurrence was 202.5 days.

CONCLUSION

The number of patients in this study was small; therefore, a long-term survival benefit was not demonstrated. However, this study suggests that SDT has potential as a treatment for rGBM and warrants further exploration. Trial information: Chinese Clinical Trial Registry ( http://www.chictr.org.cn/ ): ChiCTR2200065992. November 2, 2022, retrospectively registered.

Large-Volume Focused-Ultrasound Mild Hyperthermia for Improving Blood-Brain Tumor Barrier Permeability Application

Mild hyperthermia can locally enhance permeability of the blood-tumor barrier in brain tumors, improving delivery of antitumor nanodrugs. However, a clinical transcranial focused ultrasound (FUS) system does not provide this modality yet. The study aimed at the development of the transcranial FUS technique dedicated for large-volume mild hyperthermia in the brain. Acoustic pressure, multiple-foci, temperature and thermal dose induced by FUS were simulated in the brain through the skull. A 1-MHz, 114-element, spherical helmet transducer was fabricated to verify large-volume hyperthermia in the phantom. The simulated results showed that two foci were simultaneously formed at (2, 0, 0) and (−2, 0, 0) and at (0, 2, 0) and (0, −2, 0), using the phases of focusing pattern 1 and the phases of focusing pattern 2, respectively. Switching two focusing patterns at 5 Hz produced a hyperthermic zone with an ellipsoid of 7 mm × 6 mm × 11 mm in the brain and the temperature was 41–45 °C in the ellipsoid as the maximum intensity was 150 W/cm² and sonication time was 3 min. The phased array driven by switching two mode phases generated a 41 °C-contour region of 10 ± 1 mm × 8 ± 2 mm × 13 ± 2 mm in the phantom after 3-min sonication. Therefore, we have demonstrated our developed FUS technique for large-volume mild hyperthermia.

Numerical investigation of the energy distribution of Low-intensity transcranial focused ultrasound neuromodulation for hippocampus

OBJECTIVE

Ultrasonic neuromodulation as a safe and non-invasive brain stimulation method that delivers a low-intensity, focused ultrasound to nervous system tissue in a targeted area of the brain. The objective of this study is to numerically investigate the ultrasound wave propagation and the energy distribution within the brain tissues using customized single element focused ultrasound transducers (SEFT), targeting the hippocampus.

METHODS

A high resolution detailed human head model with seven tissue types was constructed from magnetic resonance imaging (MRI). A full-wave finite-difference time-domain simulation platform, Sim4life, was then used to simulate a 3D non-linear ultrasound wave equation to the specific region of interest, the hippocampus. Three customized SEFT were used to test the effect of transducer positions, and another customized transducer was used to compare the sensitivity effect on heterogeneous and homogeneous brain models. Finally, the sensitivity and performance of low intensity focusing ultrasound stimulation were evaluated.

RESULTS

An optimized application of SEFT was customized to deliver 100 W/m² intensity of energy deposition at the hippocampus region. About 85.65% of the generated volume beam was delivered to the targeted hippocampus region and the beam overlap parameter was affected by different transducer positions. Deflection angle changes of SEFT at the range of ± 5% did not have a significant effect on energy delivery and position displacement. Only 0.5% of peak pressure change was observed between heterogeneous and homogeneous brain models. The sensitivity analysis also showed that the sound speed is the most influential acoustic parameter.

SIGNIFICANCE

This study demonstrated that ultrasound neuromodulation targeting the depth brain tissue of the hippocampus could be a potential and promising alternative method to some non-acoustic brain stimulation modalities. In the numerical study of ultrasound brain stimulations, ultrasound parameters and the brain model need to be properly determined to simulate the ultrasonic neuromodulations.

Feasibility of Upper Cranial Nerve Sonication in Human Application via Neuronavigated Single-Element Pulsed Focused Ultrasound

Sonicating deep brain regions with pulsed focused ultrasound using magnetic resonance imaging-guided neuronavigation single-element piezoelectric transducers is a new area of exploration for neuromodulation. Upper cranial nerves such as the trigeminal nerve and other nerves responsible for sensory/motor functions in the head may be potential targets for ultrasound pain therapy. The location of upper cranial nerves close to the skull base poses additional challenges when compared with conventional cortical or middle brain targets. In the work described here, a series of computational and empirical testing methods using human skull specimens were conducted to assess the feasibility of sonicating the trigeminal pathway near the sphenoid bone region. The results indicate a transducer with a focal length of 120 mm and diameter of 85 mm (350 kHz) can deliver sonication to upper cranial nerve regions with spatial accuracy comparable to that of focused ultrasound brain targets used in previous human studies. Temperature measurements in cortical bone and in the skull base with embedded thermocouples yield evidence of minimal bone heating. Conventional pulse parameters were found to cause reverberation interference patterns near the cranial floor; therefore, changes in pulse cycles and pulse repetition frequency were examined for reducing standing waves. Limitations and considerations for conducting ultradeep focal targeting in human applications are discussed.

Combined Therapy Planning, Real-Time Monitoring, and Low Intensity Focused Ultrasound Treatment Using a Diagnostic Imaging Array

Low intensity focused ultrasound (FUS) therapies use low intensity focused ultrasound waves, typically in combination with microbubbles, to non-invasively induce a variety of therapeutic effects. FUS therapies require pre-therapy planning and real-time monitoring during treatment to ensure the FUS beam is correctly targeted to the desired tissue region. To facilitate more streamlined FUS treatments, we present a system for pre-therapy planning, real-time FUS beam visualization, and low intensity FUS treatment using a single diagnostic imaging array. Therapy planning was accomplished by manually segmenting a B-mode image captured by the imaging array and calculating a sonication pattern for the treatment based on the user-input region of interest. For real-time monitoring, the imaging array transmitted a visualization pulse which was focused to the same location as the FUS therapy beam and ultrasonic backscatter from this pulse was used to reconstruct the intensity field of the FUS beam. The therapy planning and beam monitoring techniques were demonstrated in a tissue-mimicking phantom and in a rat tumor in vivo while a mock FUS treatment was carried out. The FUS pulse from the imaging array was excited with an MI of 0.78, which suggests that the array could be used to administer select low intensity FUS treatments involving microbubble activation.

Numerical Evaluation of the Effects of Transducer Displacement on Transcranial Focused Ultrasound in the Rat Brain

Focused ultrasound is a promising therapeutic technique, as it involves the focusing of an ultrasonic beam with sufficient acoustic energy into a target brain region with high precision. Low-intensity ultrasound transmission by a single-element transducer is mostly established for neuromodulation applications and blood–brain barrier disruption for drug delivery. However, transducer positioning errors can occur without fine control over the sonication, which can affect repeatability and lead to reliability problems. The objective of this study was to determine whether the target brain region would be stable under small displacement (0.5 mm) of the transducer based on numerical simulations. Computed-tomography-derived three-dimensional models of a rat head were constructed to investigate the effects of transducer displacement in the caudate putamen (CP) and thalamus (TH). Using three different frequencies (1.1, 0.69, and 0.25 MHz), the transducer was displaced by 0.5 mm in each of the following six directions: superior, interior, anterior, posterior, left, and right. The maximum value of the intracranial pressure field was calculated, and the targeting errors were determined by the full-width-at-half-maximum (FWHM) overlap between the free water space (FWHM_water) and transcranial transmission (FWHM_base). When the transducer was positioned directly above the target region, a clear distinction between the target regions was observed, resulting in 88.3%, 81.5%, and 84.5% FWHM_water for the CP and 65.6%, 76.3%, and 64.4% FWHM_water for the TH at 1.1, 0.69, and 0.25 MHz, respectively. Small transducer displacements induced both enhancement and reduction of the peak pressure and targeting errors, compared with when the transducer was displaced in water. Small transducer displacement to the left resulted in the lowest stability, with 34.8% and 55.0% targeting accuracy (FWHM_water) at 1.1 and 0.69 MHz in the TH, respectively. In addition, the maximum pressure was reduced by up to 11% by the transducer displacement. This work provides the targeting errors induced by transducer displacements through a preclinical study and recommends that attention be paid to determining the initial sonication foci in the transverse plane in the cases of small animals.

Technical Comparison of Treatment Efficiency of Magnetic Resonance-Guided Focused Ultrasound Thalamotomy and Pallidotomy in Skull Density Ratio-Matched Patient Cohorts

Objective

MR-guided focused ultrasound (MRgFUS) is increasingly being used to treat patients with essential tremor (ET) and Parkinson's disease (PD) with thalamotomy and pallidotomy, respectively. Pallidotomy is performed off-center within the cranium compared to thalamotomy and may present challenges to therapeutic lesioning due to this location. However, the impact of target location on treatment efficiency and ability to create therapeutic lesions has not been studied. This study aimed to compare the physical efficiency of MRgFUS thalamotomy and pallidotomy.

Methods

Treatment characteristics were compared between patients treated with thalamotomy (n = 20) or pallidotomy (n = 20), matched by skull density ratios (SDR). Aspects of treatment efficiency were compared between these groups. Demographic and comparative statistics were conducted to assess these differences. Acoustic field simulations were performed to compare and validate the simulated temperature profile for VIM and GPi ablation.

Results

Lower SDR values were associated with greater energy requirement for thalamotomy (R² = 0.197, p = 0.049) and pallidotomy (R² = 0.342, p = 0.007). The impact of low SDR on efficiency reduction was greater for pallidotomy, approaching significance (p = 0.061). A nearly two-fold increase in energy was needed to reach 50°C in pallidotomy (10.9kJ) than in thalamotomy (5.7kJ), (p = 0.002). Despite lower energy requirement, the maximum average temperature reached was higher in thalamotomy (56.7°C) than in pallidotomy (55.0°C), (p = 0.017). Mean incident angle of acoustic beams was lesser in thalamotomy (12.7°) than in pallidotomy (18.6°), (p < 0.001). For all patients, a lesser mean incident angle correlated with a higher maximum average temperature reached (R² = 0.124, p = 0.026), and less energy needed to reach 50°C (R²=0.134, p = 0.020). Greater skull thickness was associated with a higher maximum energy for a single sonication for thalamotomy (R² = 0.206, p = 0.045) and pallidotomy (R² = 0.403, p = 0.003). An acoustic and temperature field simulation validated similar findings for thalamotomy and pallidotomy in a single patient.

Conclusion

The centrally located VIM offers a more efficient location for therapeutic lesioning compared to GPi pallidotomy in SDR matched cohort of patients. The impact on therapeutic lesioning with lower SDR may be greater for pallidotomy patients. As newer off-center targets are investigated, these findings can inform patient selection and treatment requirements for lesion production.

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.

Transcranial MR-Guided Histotripsy System

Histotripsy has been previously shown to treat a wide range of locations through excised human skulls in vitro. In this article, a transcranial magnetic resonance (MR)-guided histotripsy (tcMRgHt) system was developed, characterized, and tested in the in vivo pig brain through an excised human skull. A 700-kHz, 128-element MR-compatible phased-array ultrasound transducer with a focal depth of 15 cm was designed and fabricated in-house. Support structures were also constructed to facilitate transcranial treatment. The tcMRgHt array was acoustically characterized with a peak negative pressure up to 137 MPa in free field, 72 MPa through an excised human skull with aberration correction, and 48.4 MPa without aberration correction. The electronic focal steering range through the skull was 33.5 mm laterally and 50 mm axially, where a peak negative pressure above the 26-MPa cavitation intrinsic threshold can be achieved. The MR compatibility of the tcMRgHt system was assessed quantitatively using SNR, B0 field map, and B1 field map in a clinical 3T magnetic resonance imaging (MRI) scanner. Transcranial treatment using electronic focal steering was validated in red blood cell phantoms and in vivo pig brain through an excised human skull. In two pigs, targeted cerebral tissue was successfully treated through the human skull as confirmed by MRI. Excessive bleeding or edema was not observed in the peri-target zones by the time of pig euthanasia. These results demonstrated the feasibility of using this preclinical tcMRgHt system for in vivo transcranial treatment in a swine model.

A Generalized Split-Step Angular Spectrum Method for Efficient Simulation of Wave Propagation in Heterogeneous Media

Angular spectrum (AS) methods enable efficient calculation of wave propagation from one plane to another inside homogeneous media. For wave propagation in heterogeneous media such as biological tissues, AS methods cannot be applied directly. Split-stepping techniques decompose the heterogeneous domain into homogeneous and perturbation parts, and provide a solution for forward wave propagation by propagating the incident wave in both frequency-space and frequency-wavenumber domains. Recently, a split-step hybrid angular spectrum (HAS) method was proposed for plane wave propagation of focused ultrasound beams. In this study, we extend these methods to enable simulation of acoustic pressure field for an arbitrary source distribution, by decomposing the source and reflection spectra into orthogonal propagation direction components, propagating each component separately, and summing all components to get the total field. We show that our method can efficiently simulate the pressure field of arbitrary sources in heterogeneous media. The accuracy of the method was analyzed comparing the resultant pressure field with pseudospectral time domain (PSTD) solution for breast tomography and hemispherical transcranial-focused ultrasound simulation models. Eighty times acceleration was achieved for a 3-D breast simulation model compared to PSTD solution with 0.005 normalized root mean-squared difference (NRMSD) between two solutions. For the hemispherical phased array, aberrations due to skull were accurately calculated in a single simulation run as evidenced by the resultant-focused ultrasound beam simulations, which had 0.001 NRMSD with 40 times acceleration factor compared to the PSTD method.

Comparison Between Ray-Tracing and Full-Wave Simulation for Transcranial Ultrasound Focusing on a Clinical System Using the Transfer Matrix Formalism

Only one high-intensity focused ultrasound device has been clinically approved for transcranial brain surgery at the time of writing. The device operates within 650 and 720 kHz and corrects the phase distortions induced by the skull of each patient using a multielement phased array. Phase correction is estimated adaptively using a proprietary algorithm based on computed-tomography (CT) images of the patient's skull. In this article, we assess the performance of the phase correction computed by the clinical device and compare it to: 1) the correction obtained with a previously validated full-wave simulation algorithm using an open-source pseudo-spectral toolbox and 2) a hydrophone-based correction performed invasively to measure the aberrations induced by the skull at 650 kHz. For the full-wave simulation, three different mappings between CT Hounsfield units and the longitudinal speed of sound inside the skull were tested. All methods are compared with the exact same setup due to transfer matrices acquired with the clinical system for N = 5 skulls and T = 2 different targets for each skull. We show that the clinical ray-tracing software and the full-wave simulation restore, respectively, 84% ± 5% and 86% ± 5% of the pressure obtained with hydrophone-based correction for targets located in central brain regions. On the second target (off-center), we also report that the performance of both algorithms degrades when the average incident angles of the acoustic beam at the skull surface increase. When incident angles are higher than 20°, the restored pressure drops below 75% of the pressure restored with hydrophone-based correction.

Ultrasound for the Brain: A Review of Physical and Engineering Principles, and Clinical Applications

The emergence of new ultrasound technologies has improved our understanding of the brain functions and offered new opportunities for the treatment of brain diseases. Ultrasound has become a valuable tool in preclinical animal and clinical studies as it not only provides information about the structure and function of brain tissues but can also be used as a therapy alternative for brain diseases. High-resolution cerebral flow images with high sensitivity can be acquired using novel functional ultrasound and super-resolution ultrasound imaging techniques. The noninvasive treatment of essential tremors has been clinically approved and it has been demonstrated that the ultrasound technology can revolutionize the currently existing treatment methods. Microbubble-mediated ultrasound can remotely open the blood-brain barrier enabling targeted drug delivery in the brain. More recently, ultrasound neuromodulation received a great amount of attention due to its noninvasive and deep penetration features and potential therapeutic benefits. This review provides a thorough introduction to the current state-of-the-art research on brain ultrasound and also introduces basic knowledge of brain ultrasound including the acoustic properties of the brain/skull and engineering techniques for ultrasound. Ultrasound is expected to play an increasingly important role in the diagnosis and therapy of brain diseases.

Echo-Focusing in Transcranial Focused Ultrasound Thalamotomy for Essential Tremor: A Feasibility Study

BACKGROUND

Transcranial magnetic resonance-guided focused ultrasound (TcMRgFUS) systems currently employ computed tomography (CT)-based aberration corrections, which may provide suboptimal trans-skull focusing.

OBJECTIVES

The objective of this study was to evaluate a contrast agent microbubble imaging-based transcranial focusing method, echo-focusing (EF), during TcMRgFUS for essential tremor.

METHODS

A clinical trial of TcMRgFUS thalamotomy using EF for the treatment of essential tremor was conducted (NCT03935581; funded by InSightec [Tirat Carmel, Israel]). Patients (n = 12) were injected with Definity (Lantheus Medical Imaging, North Billerica, MA) microbubbles, and EF was performed using a research feature add-on to a commercial TcMRgFUS system (ExAblate Neuro, InSightec). Subablative thermal sonications carried out using (1) EF and (2) CT-based aberration corrections were compared via magnetic resonance thermometry, and the optimal focusing method for each patient was employed for TcMRgFUS thalamotomy.

RESULTS

EF aberration corrections provided increased sonication efficiency, decreased focal size, and equivalent targeting accuracy relative to CT-based focusing. EF aberration corrections were employed successfully for lesion formation in all 12 patients, 3 of whom had previously undergone unsuccessful TcMRgFUS thalamotomy via CT-based focusing. There were no adverse events related directly to the EF procedure.

CONCLUSIONS

EF is feasible and appears safe during TcMRgFUS thalamotomy for essential tremor and improves on the trans-skull focal quality provided by existing CT-based focusing methods. © 2020 International Parkinson and Movement Disorder Society.

Neuronavigated Repetitive Transcranial Ultrasound Stimulation Induces Long-Lasting and Reversible Effects on Oculomotor Performance in Non-human Primates

Since the late 2010s, Transcranial Ultrasound Stimulation (TUS) has been used experimentally to carryout safe, non-invasive stimulation of the brain with better spatial resolution than Transcranial Magnetic Stimulation (TMS). This innovative stimulation method has emerged as a novel and valuable device for studying brain function in humans and animals. In particular, single pulses of TUS directed to oculomotor regions have been shown to modulate visuomotor behavior of non-human primates during 100 ms ultrasound pulses. In the present study, a sustained effect was induced by applying 20-s trains of neuronavigated repetitive Transcranial Ultrasound Stimulation (rTUS) to oculomotor regions of the frontal cortex in three non-human primates performing an antisaccade task. With the help of MRI imaging and a frame-less stereotactic neuronavigation system (SNS), we were able to demonstrate that neuronavigated TUS (outside of the MRI scanner) is an efficient tool to carry out neuromodulation procedures in non-human primates. We found that, following neuronavigated rTUS, saccades were significantly modified, resulting in shorter latencies compared to no-rTUS trials. This behavioral modulation was maintained for up to 20 min. Oculomotor behavior returned to baseline after 18-31 min and could not be significantly distinguished from the no-rTUS condition. This study is the first to show that neuronavigated rTUS can have a persistent effect on monkey behavior with a quantified return-time to baseline. The specificity of the effects could not be explained by auditory confounds.

Ultrafast three-dimensional microbubble imaging in vivo predicts tissue damage volume distributions during nonthermal brain ablation

Transcranial magnetic resonance imaging (MRI)-guided focused ultrasound (FUS) thermal ablation is under clinical investigation for non-invasive neurosurgery, though its use is restricted to central brain targets due primarily to skull heating effects. The combination of FUS and contrast agent microbubbles greatly reduces the ultrasound exposure levels needed to ablate brain tissue and may help facilitate the use of transcranial FUS ablation throughout the brain. However, sources of variability exist during microbubble-mediated FUS procedures that necessitate the continued development of systems and methods for online treatment monitoring and control, to ensure that excessive and/or off-target bioeffects are not induced from the exposures. Methods: Megahertz-rate three-dimensional (3D) microbubble imaging in vivo was performed during nonthermal ablation in rabbit brain using a clinical-scale prototype transmit/receive hemispherical phased array system. Results:In-vivo volumetric acoustic imaging over microsecond timescales uncovered spatiotemporal microbubble dynamics hidden by conventional whole-burst temporal averaging. Sonication-aggregate ultrafast 3D source field intensity data were predictive of microbubble-mediated tissue damage volume distributions measured post-treatment using MRI and confirmed via histopathology. Temporal under-sampling of acoustic emissions, which is common practice in the field, was found to impede performance and highlighted the importance of capturing adequate data for treatment monitoring and control purposes. Conclusion: The predictive capability of ultrafast 3D microbubble imaging, reported here for the first time, will enable future microbubble-mediated FUS treatments with unparalleled precision and accuracy, and will accelerate the clinical translation of nonthermal tissue ablation procedures both in the brain and throughout the body.

A Basal Forebrain-Cingulate Circuit in Macaques Decides It Is Time to Act

The medial frontal cortex has been linked to voluntary action, but an explanation of why decisions to act emerge at particular points in time has been lacking. We show that, in macaques, decisions about whether and when to act are predicted by a set of features defining the animal’s current and past context; for example, respectively, cues indicating the current average rate of reward and recent previous voluntary action decisions. We show that activity in two brain areas—the anterior cingulate cortex and basal forebrain—tracks these contextual factors and mediates their effects on behavior in distinct ways. We use focused transcranial ultrasound to selectively and effectively stimulate deep in the brain, even as deep as the basal forebrain, and demonstrate that alteration of activity in the two areas changes decisions about when to act.

•

Likelihood and timing of voluntary action in macaques can be partially predicted

•

Recent experience and present context influence when voluntary action occurs

•

A basal forebrain-cingulate circuit mediated effects of these factors on behavior

•

Stimulation of this circuit by ultrasound changed decisions about when to act

Method to optimize the placement of a single-element transducer for transcranial focused ultrasound

BACKGROUND AND OBJECTIVE

Transcranial focused ultrasound (tFUS) is a promising neuromodulation technique because of its non-invasiveness and high spatial resolution (within millimeter scale). However, the presence of the skull can lead to disrupting and shifting the acoustic focus in the brain. In this study, we propose a computationally efficient way to determine the optimal position of a single-element focused ultrasound transducer which can effectively deliver acoustic energy to the brain target. We hypothesized that the placement of a single element transducer with the lowest average reflection coefficient would be the optimal position.

METHODS

The reflection coefficient is defined by the ratio of the amplitude of the reflected wave to the incident wave. To calculate the reflection coefficient, we assumed ultrasound waves as straight lines (beam lines). At each beam line, the reflection coefficient was calculated from the incidence angle at the skull interface (outer/inner skull surfaces). The average reflection coefficient (ARC) was calculated at each possible placement of the transducer using a custom-built software. For comparison purposes, acoustic simulations (k-Wave MATLAB toolbox) which numerically solved the linear wave equation were performed with the same transducer positions used in the ARC calculation. In addition, the experimental validation of our proposed method was also performed by measuring acoustic wave propagation through the calvaria skull phantom in water. The accuracy of our method was defined as the distance between the two optimal transducer placements which were determined from the acoustic simulations and from the ARC method.

RESULT

Simulated acoustic pressure distribution corresponding to each ARC showed an inverse relationship with peak acoustic pressures produced in the brain. In comparison to the acoustic simulations, the accuracy of our method was 5.07 ± 4.27 mm when targeting the cortical region in the brain. The computing time of ARC calculations were 0.08% of the time required for acoustic pressure simulations.

CONCLUSION

We calculated the ARC to find the optimal position of the tFUS transducer used in the present study. The optimal placement of the transducer was found when the ARC was the lowest. Our numerical and experimental results showed that the proposed ARC method can effectively be used to find the optimal position of a single-element tFUS transducer for targeting the cortex region of the brain in a computationally inexpensive way.

Offline impact of transcranial focused ultrasound on cortical activation in primates

To understand brain circuits it is necessary both to record and manipulate their activity. Transcranial ultrasound stimulation (TUS) is a promising non-invasive brain stimulation technique. To date, investigations report short-lived neuromodulatory effects, but to deliver on its full potential for research and therapy, ultrasound protocols are required that induce longer-lasting 'offline' changes. Here, we present a TUS protocol that modulates brain activation in macaques for more than one hour after 40 s of stimulation, while circumventing auditory confounds. Normally activity in brain areas reflects activity in interconnected regions but TUS caused stimulated areas to interact more selectively with the rest of the brain. In a within-subject design, we observe regionally specific TUS effects for two medial frontal brain regions - supplementary motor area and frontal polar cortex. Independently of these site-specific effects, TUS also induced signal changes in the meningeal compartment. TUS effects were temporary and not associated with microstructural changes.

Offline impact of transcranial focused ultrasound on cortical activation in primates

To understand brain circuits it is necessary both to record and manipulate their activity. Transcranial ultrasound stimulation (TUS) is a promising non-invasive brain stimulation technique. To date, investigations report short-lived neuromodulatory effects, but to deliver on its full potential for research and therapy, ultrasound protocols are required that induce longer-lasting 'offline' changes. Here, we present a TUS protocol that modulates brain activation in macaques for more than one hour after 40 s of stimulation, while circumventing auditory confounds. Normally activity in brain areas reflects activity in interconnected regions but TUS caused stimulated areas to interact more selectively with the rest of the brain. In a within-subject design, we observe regionally specific TUS effects for two medial frontal brain regions - supplementary motor area and frontal polar cortex. Independently of these site-specific effects, TUS also induced signal changes in the meningeal compartment. TUS effects were temporary and not associated with microstructural changes.

Spatial filters suppress ripple artifacts in the computation of acoustic fields with the angular spectrum method

The angular spectrum method (ASM) is an effective tool for propagating wave fields between parallel planes through decomposition of the field into a series of independent plane waves. One source of error is interference from mirror sources introduced through the inherent periodicity of the fast Fourier transform (FFT) used to implement this method numerically. Here, spatial filters attenuate waves propagating at large angles, which are sensitive to mirror sources. Simulations show that this suppresses the ripple artifact whilst preserving the accuracy of the ASM-computed fields. To achieve comparable performance without filtering requires up to a 13.5-fold increase in computation time.

Simulating transvertebral ultrasound propagation with a multi-layered ray acoustics model

The simulation accuracy of transvertebral ultrasound propagation using a multi-layered ray acoustics model based on CT-derived vertebral geometry was investigated through comparison with experimental measurements of pressure fields in ex vivo human vertebral foramen. A spherically focused transducer (5 cm diameter, f-number 1.2, 514 kHz) was geometrically focused to the centre of individual thoracic vertebral foramen, through the posterior bony elements. Transducer propagation paths through the laminae and the spinous processes were tested. Simulation transducer-vertebra configurations were registered to experiment transducer-vertebra configurations, and simulation accuracy of the simulation model was evaluated for predicting maximum transmitted pressure to the canal, voxel pressure in the canal, and focal distortion. Accuracy in predicting maximum transmitted pressure was calculated by vertebra, and it is shown that simulation predicts maximum pressure with a greater degree of accuracy than a vertebra-specific insertion loss. Simulation error in voxel pressure was evaluated using root-mean-square error and cross-correlation, and found to be similar to the water-only case. Simulation accuracy in predicting focal distortion was evaluated by comparing experiment and simulation maximum pressure location and weighted  >50% focal volume location. Average simulation error across all measurements and simulations in maximum pressure location and weighted  >50% focal volume location were 2.3 mm and 1.5 mm, respectively. These errors are small relative to the dimensions of the transducer focus (4.9 mm full width half maximum), the spinal cord (10 mm diameter), and vertebral canal diameter (15-20 mm diameter). These results suggest that ray acoustics can be applied to simulating transvertebral ultrasound propagation.

Inside/outside the brain binary cavitation localization based on the lowpass filter effect of the skull on the harmonic content: a proof of concept study

Cavitation activity induced by ultrasound may occur during high intensity focused ultrasound (HIFU) treatment, due to bubble nucleation under high peak negative pressure, and during blood-brain-barrier (BBB) disruption, due to injected ultrasound contrast agents (UCAs). Such microbubble activity has to be monitored to assess the safety and efficiency of ultrasonic brain treatments. In this study, we aim at assessing whether cavitation occurs within cerebral tissue by binary discriminating cavitation activity originating from the inside or the outside of the skull. The results were obtained from both in vitro experiments mimicking BBB opening, by using UCA flow, and in vitro thermal necrosis in calf brain samples. The sonication was applied using a 1 MHz focused transducer and the acoustic response of the microbubbles was recorded with a wideband passive cavitation detector. The spectral content of the recorded signal was used to localize microbubble activity. Since the skull acts as a low pass filter, the ratio of high harmonics to low harmonics is lower for cavitation events located inside the skull compared to events outside the skull. Experiments showed that the ratio of the 5/2 ultraharmonic to the 1/2 subharmonic for binary localization cavitation activity achieves 100% sensitivity and specificity for both monkey and human skulls. The harmonic ratio of the fourth to the second harmonic provided 100% sensitivity and 96% and 46% specificity on a non-human primate for thermal necrosis and BBB opening, respectively. Nonetheless, the harmonic ratio remains promising for human applications, as the experiments showed 100% sensitivity and 100% specificity for both thermal necrosis and BBB opening through the human skull. The study requires further validation on a larger number of skull samples.

A 200-1380-kHz Quadrifrequency Focused Ultrasound Transducer for Neurostimulation in Rodents and Primates: Transcranial In Vitro Calibration and Numerical Study of the Influence of Skull Cavity

Low intensity transcranial focused ultrasound has been demonstrated to produce neuromodulation in both animals and humans. Primarily for technical reasons, frequency is one of the most poorly investigated critical wave parameters. We propose the use of a quadri-band transducer capable of operating at 200, 320, 850, and 1380 kHz for further investigation of the frequency dependence of neuromodulation efficacy while keeping the position of the transducer fixed with respect to the subject's head. This paper presents the results of the transducer calibration in water, in vitro transmission measurements through a monkey skull flap, 3-D simulations based on both a μ -computed tomography ( μ CT)-scan of a rat and on CT-scans of two macaques. A maximum peak pressure greater than 0.52 MPa is expected at each frequency in rat and macaque heads. According to the literature, our transducer can achieve neuromodulation in rodents and primates at each four frequencies. The impact of standing waves is shown to be most prominent at the lowest frequencies.

A Tikhonov Regularization Scheme for Focus Rotations With Focused Ultrasound-Phased Arrays

Phased arrays have a wide range of applications in focused ultrasound therapy. By using an array of individually driven transducer elements, it is possible to steer a focus through space electronically and compensate for acoustically heterogeneous media with phase delays. In this paper, the concept of focusing an ultrasound-phased array is expanded to include a method to control the orientation of the focus using a Tikhonov regularization scheme. It is then shown that the Tikhonov regularization parameter used to solve the ill-posed focus rotation problem plays an important role in the balance between quality focusing and array efficiency. Finally, the technique is applied to the synthesis of multiple foci, showing that this method allows for multiple independent spatial rotations.

A multi-frequency sparse hemispherical ultrasound phased array for microbubble-mediated transcranial therapy and simultaneous cavitation mapping

Focused ultrasound (FUS) phased arrays show promise for non-invasive brain therapy. However, the majority of them are limited to a single transmit/receive frequency and therefore lack the versatility to expose and monitor the treatment volume. Multi-frequency arrays could offer variable transmit focal sizes under a fixed aperture, and detect different spectral content on receive for imaging purposes. Here, a three-frequency (306, 612, and 1224 kHz) sparse hemispherical ultrasound phased array (31.8 cm aperture; 128 transducer modules) was constructed and evaluated for microbubble-mediated transcranial therapy and simultaneous cavitation mapping. The array is able to perform effective electronic beam steering over a volume spanning (-40, 40) and (-30, 50) mm in the lateral and axial directions, respectively. The focal size at the geometric center is approximately 0.9 (2.1) mm, 1.7 (3.9) mm, and 3.1 (6.5) mm in lateral (axial) pressure full width at half maximum (FWHM) at 1224, 612, and 306 kHz, respectively. The array was also found capable of dual-frequency excitation and simultaneous multi-foci sonication, which enables the future exploration of more complex exposure strategies. Passive acoustic mapping of dilute microbubble clouds demonstrated that the point spread function of the receive array has a lateral (axial) intensity FWHM between 0.8-3.5 mm (1.7-11.7 mm) over a volume spanning (-25, 25) mm in both the lateral and axial directions, depending on the transmit/receive frequency combination and the imaging location. The device enabled both half and second harmonic imaging through the intact skull, which may be useful for improving the contrast-to-tissue ratio or imaging resolution, respectively. Preliminary in vivo experiments demonstrated the system's ability to induce blood-brain barrier opening and simultaneously spatially map microbubble cavitation activity in a rat model. This work presents a tool to investigate optimal strategies for non-thermal FUS brain therapy and concurrent microbubble cavitation monitoring through the availability of multiple frequencies.

Nonthermal ablation of deep brain targets: A simulation study on a large animal model

PURPOSE

Thermal ablation with transcranial MRI-guided focused ultrasound (FUS) is currently limited to central brain targets because of heating and other beam effects caused by the presence of the skull. Recently, it was shown that it is possible to ablate tissues without depositing thermal energy by driving intravenously administered microbubbles to inertial cavitation using low-duty-cycle burst sonications. A recent study demonstrated that this ablation method could ablate tissue volumes near the skull base in nonhuman primates without thermally damaging the nearby bone. However, blood-brain disruption was observed in the prefocal region, and in some cases, this region contained small areas of tissue damage. The objective of this study was to analyze the experimental model with simulations and to interpret the cause of these effects.

METHODS

The authors simulated prior experiments where nonthermal ablation was performed in the brain in anesthetized rhesus macaques using a 220 kHz clinical prototype transcranial MRI-guided FUS system. Low-duty-cycle sonications were applied at deep brain targets with the ultrasound contrast agent Definity. For simulations, a 3D pseudospectral finite difference time domain tool was used. The effects of shear mode conversion, focal steering, skull aberrations, nonlinear propagation, and the presence of skull base on the pressure field were investigated using acoustic and elastic wave propagation models.

RESULTS

The simulation results were in agreement with the experimental findings in the prefocal region. In the postfocal region, however, side lobes were predicted by the simulations, but no effects were evident in the experiments. The main beam was not affected by the different simulated scenarios except for a shift of about 1 mm in peak position due to skull aberrations. However, the authors observed differences in the volume, amplitude, and distribution of the side lobes. In the experiments, a single element passive cavitation detector was used to measure the inertial cavitation threshold and to determine the pressure amplitude to use for ablation. Simulations of the detector's acoustic field suggest that its maximum sensitivity was in the lower part of the main beam, which may have led to excessive exposure levels in the experiments that may have contributed to damage in the prefocal area.

CONCLUSIONS

Overall, these results suggest that case-specific full wave simulations before the procedure can be useful to predict the focal and the prefocal side lobes and the extent of the resulting bioeffects produced by nonthermal ablation. Such simulations can also be used to optimally position passive cavitation detectors. The disagreement between the simulations and the experiments in the postfocal region may have been due to shielding of the ultrasound field due to microbubble activity in the focal region. Future efforts should include the effects of microbubble activity and vascularization on the pressure field.

Image-guided ultrasound phased arrays are a disruptive technology for non-invasive therapy

Focused ultrasound offers a non-invasive way of depositing acoustic energy deep into the body, which can be harnessed for a broad spectrum of therapeutic purposes, including tissue ablation, the targeting of therapeutic agents, and stem cell delivery. Phased array transducers enable electronic control over the beam geometry and direction, and can be tailored to provide optimal energy deposition patterns for a given therapeutic application. Their use in combination with modern medical imaging for therapy guidance allows precise targeting, online monitoring, and post-treatment evaluation of the ultrasound-mediated bioeffects. In the past there have been some technical obstacles hindering the construction of large aperture, high-power, densely-populated phased arrays and, as a result, they have not been fully exploited for therapy delivery to date. However, recent research has made the construction of such arrays feasible, and it is expected that their continued development will both greatly improve the safety and efficacy of existing ultrasound therapies as well as enable treatments that are not currently possible with existing technology. This review will summarize the basic principles, current statures, and future potential of image-guided ultrasound phased arrays for therapy.

Computational exploration of wave propagation and heating from transcranial focused ultrasound for neuromodulation

OBJECTIVE

While ultrasound is largely established for use in diagnostic imaging, its application for neuromodulation is relatively new and crudely understood. The objective of the present study was to investigate the effects of tissue properties and geometry on the wave propagation and heating in the context of transcranial neuromodulation.

APPROACH

A computational model of transcranial-focused ultrasound was constructed and validated against empirical data. The models were then incrementally extended to investigate a number of issues related to the use of ultrasound for neuromodulation, including the effect on wave propagation of variations in geometry of skull and gyral anatomy as well as the effect of multiple tissue and media layers, including scalp, skull, CSF, and gray/white matter. In addition, a sensitivity analysis was run to characterize the influence of acoustic properties of intracranial tissues. Finally, the heating associated with ultrasonic stimulation waveforms designed for neuromodulation was modeled.

MAIN RESULTS

The wave propagation of a transcranially focused ultrasound beam is significantly influenced by the cranial domain. The half maximum acoustic beam intensity profiles are insensitive overall to small changes in material properties, though the inclusion of sulci in models results in greater peak intensity values compared to a model without sulci (1%-30% greater). Finally, heating using currently employed stimulation parameters in humans is highest in bone (0.16 °C) and is negligible in brain (4.27 × 10(-3) °C) for a 0.5 s exposure.

SIGNIFICANCE

Ultrasound for noninvasive neuromodulation holds great promise and appeal for its non-invasiveness, high spatial resolution and deep focal lengths. Here we show gross brain anatomy and biological material properties to have limited effect on ultrasound wave propagation and to result in safe heating levels in the skull and brain.

Ultrasonic neuromodulation

Ultrasonic waves can be non-invasively steered and focused into mm-scale regions across the human body and brain, and their application in generating controlled artificial modulation of neuronal activity could therefore potentially have profound implications for neural science and engineering. Ultrasonic neuro-modulation phenomena were experimentally observed and studied for nearly a century, with recent discoveries on direct neural excitation and suppression sparking a new wave of investigations in models ranging from rodents to humans. In this paper we review the physics, engineering and scientific aspects of ultrasonic fields, their control in both space and time, and their effect on neuronal activity, including a survey of both the field's foundational history and of recent findings. We describe key constraints encountered in this field, as well as key engineering systems developed to surmount them. In closing, the state of the art is discussed, with an emphasis on emerging research and clinical directions.

Experimental demonstration of passive acoustic imaging in the human skull cavity using CT-based aberration corrections

PURPOSE

Experimentally verify a previously described technique for performing passive acoustic imaging through an intact human skull using noninvasive, computed tomography (CT)-based aberration corrections Jones et al. [Phys. Med. Biol. 58, 4981-5005 (2013)].

METHODS

A sparse hemispherical receiver array (30 cm diameter) consisting of 128 piezoceramic discs (2.5 mm diameter, 612 kHz center frequency) was used to passively listen through ex vivo human skullcaps (n = 4) to acoustic emissions from a narrow-band fixed source (1 mm diameter, 516 kHz center frequency) and from ultrasound-stimulated (5 cycle bursts, 1 Hz pulse repetition frequency, estimated in situ peak negative pressure 0.11-0.33 MPa, 306 kHz driving frequency) Definity™ microbubbles flowing through a thin-walled tube phantom. Initial in vivo feasibility testing of the method was performed. The performance of the method was assessed through comparisons to images generated without skull corrections, with invasive source-based corrections, and with water-path control images.

RESULTS

For source locations at least 25 mm from the inner skull surface, the modified reconstruction algorithm successfully restored a single focus within the skull cavity at a location within 1.25 mm from the true position of the narrow-band source. The results obtained from imaging single bubbles are in good agreement with numerical simulations of point source emitters and the authors' previous experimental measurements using source-based skull corrections O'Reilly et al. [IEEE Trans. Biomed. Eng. 61, 1285-1294 (2014)]. In a rat model, microbubble activity was mapped through an intact human skull at pressure levels below and above the threshold for focused ultrasound-induced blood-brain barrier opening. During bursts that led to coherent bubble activity, the location of maximum intensity in images generated with CT-based skull corrections was found to deviate by less than 1 mm, on average, from the position obtained using source-based corrections.

CONCLUSIONS

Taken together, these results demonstrate the feasibility of using the method to guide bubble-mediated ultrasound therapies in the brain. The technique may also have application in ultrasound-based cerebral angiography.

Minimizing eddy currents induced in the ground plane of a large phased-array ultrasound applicator for echo-planar imaging-based MR thermometry

BACKGROUND

The study aims to investigate different ground plane segmentation designs of an ultrasound transducer to reduce gradient field induced eddy currents and the associated geometric distortion and temperature map errors in echo-planar imaging (EPI)-based MR thermometry in transcranial magnetic resonance (MR)-guided focused ultrasound (tcMRgFUS).

METHODS

Six different ground plane segmentations were considered and the efficacy of each in suppressing eddy currents was investigated in silico and in operando. For the latter case, the segmented ground planes were implemented in a transducer mockup model for validation. Robust spoiled gradient (SPGR) echo sequences and multi-shot EPI sequences were acquired. For each sequence and pattern, geometric distortions were quantified in the magnitude images and expressed in millimeters. Phase images were used for extracting the temperature maps on the basis of the temperature-dependent proton resonance frequency shift phenomenon. The means, standard deviations, and signal-to-noise ratios (SNRs) were extracted and contrasted with the geometric distortions of all patterns.

RESULTS

The geometric distortion analysis and temperature map evaluations showed that more than one pattern could be considered the best-performing transducer. In the sagittal plane, the star (d) (3.46 ± 2.33 mm) and star-ring patterns (f) (2.72 ± 2.8 mm) showed smaller geometric distortions than the currently available seven-segment sheet (c) (5.54 ± 4.21 mm) and were both comparable to the reference scenario (a) (2.77 ± 2.24 mm). Contrasting these results with the temperature maps revealed that (d) performs as well as (a) in SPGR and EPI.

CONCLUSIONS

We demonstrated that segmenting the transducer ground plane into a star pattern reduces eddy currents to a level wherein multi-plane EPI for accurate MR thermometry in tcMRgFUS is feasible.

Comparison of analytical and numerical approaches for CT-based aberration correction in transcranial passive acoustic imaging

Computed tomography (CT)-based aberration corrections are employed in transcranial ultrasound both for therapy and imaging. In this study, analytical and numerical approaches for calculating aberration corrections based on CT data were compared, with a particular focus on their application to transcranial passive imaging. Two models were investigated: a three-dimensional full-wave numerical model (Connor and Hynynen 2004 IEEE Trans. Biomed. Eng. 51 1693-706) based on the Westervelt equation, and an analytical method (Clement and Hynynen 2002 Ultrasound Med. Biol. 28 617-24) similar to that currently employed by commercial brain therapy systems. Trans-skull time delay corrections calculated from each model were applied to data acquired by a sparse hemispherical (30 cm diameter) receiver array (128 piezoceramic discs: 2.5 mm diameter, 612 kHz center frequency) passively listening through ex vivo human skullcaps (n  =  4) to emissions from a narrow-band, fixed source emitter (1 mm diameter, 516 kHz center frequency). Measurements were taken at various locations within the cranial cavity by moving the source around the field using a three-axis positioning system. Images generated through passive beamforming using CT-based skull corrections were compared with those obtained through an invasive source-based approach, as well as images formed without skull corrections, using the main lobe volume, positional shift, peak sidelobe ratio, and image signal-to-noise ratio as metrics for image quality. For each CT-based model, corrections achieved by allowing for heterogeneous skull acoustical parameters in simulation outperformed the corresponding case where homogeneous parameters were assumed. Of the CT-based methods investigated, the full-wave model provided the best imaging results at the cost of computational complexity. These results highlight the importance of accurately modeling trans-skull propagation when calculating CT-based aberration corrections. Although presented in an imaging context, our results may also be applicable to the problem of transmit focusing through the skull.

MR-Guided Transcranial Focused Ultrasound

Previous chapters introduced the ability of using focused ultrasound to ablate tissues. It has led to various clinical applications in the treatment of uterine fibroid, prostate or liver cancers. Nevertheless, treating the brain non-invasively with focused ultrasound has been considered beyond reach for almost a century: The skull bone protects the brain from mechanical injuries, but it also reflects and refracts ultrasound, making it difficult to target the brain with focused ultrasound. Fortunately, aberration correction techniques have been developed recently and thermal lesioning in the thalamus has been achieved clinically. This chapter introduces the aberration effect of the skull bone and how it can be corrected non-invasively. It also presents the latest clinical results obtained with thermal ablation and introduces novel non-thermal approaches that could revolutionize brain therapy in the future.

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.

Temperature-dependent MR signals in cortical bone: potential for monitoring temperature changes during high-intensity focused ultrasound treatment in bone

PURPOSE

Because existing magnetic resonance thermometry techniques do not provide temperature information within bone, high-intensity focused ultrasound (HIFU) exposures in bone are monitored using temperature changes in adjacent soft tissues. In this study, the potential to monitor temperature changes in cortical bone using a short TE gradient echo sequence is evaluated.

METHODS

The feasibility of this proposed method was initially evaluated by measuring the temperature dependence of the gradient echo signal during cooling of cortical bone samples implanted with fiber-optic temperature sensors. A subsequent experiment involved heating a cortical bone sample using a clinical MR-HIFU system.

RESULTS

A consistent relationship between temperature change and the change in magnitude signal was observed within and between cortical bone samples. For the two-dimensional gradient echo sequence implemented in this study, a least-squares linear fit determined the percentage change in signal to be (0.90 ± 0.01)%/°C. This relationship was used to estimate temperature changes observed in the HIFU experiment and these temperatures agreed well with those measured from an implanted fiber-optic sensor.

CONCLUSION

This method appears capable of displaying changes related to temperature in cortical bone and could improve the safety of MR-HIFU treatments. Further investigations into the sensitivity of the technique in vivo are warranted.

In silico study of low-frequency transcranial ultrasound fields in acute ischemic stroke patients

Ultrasound in the sub-megahertz range enhances thrombolysis and may be applied transcranially to ischemic stroke patients. The consistency of transcranial insonification needs to be evaluated. Acoustic and thermal simulations based on computed-tomography (CT) scans of 20 patients were performed. An unfocused 120-kHz transducer allowed homogeneous insonification of the thrombus, and positioning based on external landmarks performed similarly to an optimized placement based on CT data. With a weakly focused 500-kHz transducer, the landmark-based positioning underperformed. The predicted inter-patient variation of in situ acoustic pressure was similar with both the 120 and 500-kHz transducers for the optimized placement (18.0-26.4% relative standard deviation). The simulated maximum acoustic pressure in intervening tissues was 2.6 ± 0.6 and 2.0 ± 0.7 times the pressure in the thrombus for the 120-kHz and 500-kHz transducers, respectively. A 1 W/cm(2) insonification of the thrombus caused a 3.8 ± 2.2 °C increase in the bone for the 120-kHz transducer, and a 13.4 ± 3.3 °C increase for the 500-kHz transducer. Contralateral local maxima up to 1.1 times the pressure amplitude in the targeted zone were predicted for the 120-kHz transducer. We established two transducer placement approaches, one based on analysis of a head CT and the other using simple external, visible landmarks. Both approaches allowed consistent insonification of the thrombus.

Numerical simulations of clinical focused ultrasound functional neurosurgery

A computational model utilizing grid and finite difference methods were developed to simulate focused ultrasound functional neurosurgery interventions. The model couples the propagation of ultrasound in fluids (soft tissues) and solids (skull) with acoustic and visco-elastic wave equations. The computational model was applied to simulate clinical focused ultrasound functional neurosurgery treatments performed in patients suffering from therapy resistant chronic neuropathic pain. Datasets of five patients were used to derive the treatment geometry. Eight sonications performed in the treatments were then simulated with the developed model. Computations were performed by driving the simulated phased array ultrasound transducer with the acoustic parameters used in the treatments. Resulting focal temperatures and size of the thermal foci were compared quantitatively, in addition to qualitative inspection of the simulated pressure and temperature fields. This study found that the computational model and the simulation parameters predicted an average of 24 ± 13% lower focal temperature elevations than observed in the treatments. The size of the simulated thermal focus was found to be 40 ± 13% smaller in the anterior-posterior direction and 22 ± 14% smaller in the inferior-superior direction than in the treatments. The location of the simulated thermal focus was off from the prescribed target by 0.3 ± 0.1 mm, while the peak focal temperature elevation observed in the measurements was off by 1.6 ± 0.6 mm. Although the results of the simulations suggest that there could be some inaccuracies in either the tissue parameters used, or in the simulation methods, the simulations were able to predict the focal spot locations and temperature elevations adequately for initial treatment planning performed to assess, for example, the feasibility of sonication. The accuracy of the simulations could be improved if more precise ultrasound tissue properties (especially of the skull bone) could be obtained.

Cavitation-based third ventriculostomy using MRI-guided focused ultrasound

OBJECT

Transcranial focused ultrasound is increasingly being investigated as a minimally invasive treatment for a range of intracranial pathologies. At higher peak rarefaction pressures than those used for thermal ablation, focused ultrasound can initiate inertial cavitation and create holes in the brain by fractionation of the tissue elements. The authors investigated the technical feasibility of using MRI-guided focused ultrasound to perform a third ventriculostomy as a possible noninvasive alternative to endoscopic third ventriculostomy for hydrocephalus.

METHODS

A craniectomy was performed in male pigs weighing 13-19 kg to expose the supratentorial brain, leaving the dura mater intact. Seven pigs were treated through the craniectomy, while 2 pigs were treated through ex vivo human skulls placed in the beam path. Registration and targeting was done using T2-weighted MRI sequences. For transcranial treatments a CT scan was used to correct the beam from aberrations due to the skull and maintain a small, high-intensity focus. Sonications were performed at both 650 kHz and 230 kHz at a range of intensities, and the in situ pressures were estimated both from simulations and experimental data to establish a threshold for tissue fractionation in the brain.

RESULTS

In craniectomized animals at 650 kHz, a peak pressure ≥ 22.7 MPa for 1 second was needed to reliably create a ventriculostomy. Transcranially at this frequency the ExAblate 4000 was unable to generate the required intensity to fractionate tissue, although cavitation was initiated. At 230 kHz, ventriculostomy was successful through the skull with a peak pressure of 8.8 MPa.

CONCLUSIONS

This is the first study to suggest that it is possible to perform a completely noninvasive third ventriculostomy using ultrasound. This may pave the way for future studies and eventually provide an alternative means for the creation of CSF communications in the brain, including perforation of the septum pellucidum or intraventricular membranes.

Surface vibration and nearby cavitation of an ex vivo bovine femur exposed to high intensity focused ultrasound

The acoustic pressure distribution, thermal ablation, and sonochemiluminescence (SCL) generated by cavitation near the surface of an ex vivo bovine femur were investigated at normal and oblique incidences of high intensity focused ultrasound (HIFU), as were the characteristics of bone surface vibrations. The acoustic pressure at the HIFU focus, the width of thermal ablation, and the SCL intensity in the pre-focal region were 1.3 MPa, 7 mm, and 454 electrons, respectively, in the control group at normal incidence, and they respectively increased to 1.5 MPa, 12 mm and 968 electrons in the presence of the bone. At oblique incidence from the left, the acoustic pressure at 3 mm to the right of the HIFU focus was 0.6 MPa and decreased to 0.4 MPa at 3 mm to the left of the focus. The thermal ablation was 20 mm in width and extended along the front surface of the bone to the right of the HIFU focus. The SCL intensity on the right of the HIFU focus was 394 electrons and was 362 electrons on the left. The presence of bone would directionally change the spatial distribution of acoustic pressure, thermal and cavitation effects for oblique incidence of HIFU.