Transcranial focused ultrasound phase correction using the hybrid angular spectrum method

Robust deep learning estimation of cortical bone porosity from MR T1-weighted images for individualized transcranial focused ultrasound planning

Objective

Transcranial focused ultrasound (tFUS) is an emerging neuromodulation approach that has been demonstrated in animals but is difficult to translate to humans because of acoustic attenuation and scattering in the skull. Optimal dose delivery requires subject-specific skull porosity estimates which has traditionally been done using CT. We propose a deep learning (DL) estimation of skull porosity from T1-weighted MRI images which removes the need for radiation-inducing CT scans.

Approach

We evaluate the impact of different DL approaches, including network architecture, input size and dimensionality, multichannel inputs, data augmentation, and loss functions. We also propose back-propagation in the mask (BIM), a method whereby only voxels inside the skull mask contribute to training. We evaluate the robustness of the best model to input image noise and MRI acquisition parameters and propagate porosity estimation errors in thousands of beam propagation scenarios.

Main results

Our best performing model is a cGAN with a ResNet-9 generator with 3D 64×64×64 inputs trained with L1 and L2 losses. The model achieved a mean absolute error of 6.9% in the test set, compared to 9.5% with the pseudo-CT of Izquierdo et al. (38% improvement) and 9.4% with the generic pixel-to-pixel image translation cGAN pix2pix (36% improvement). Acoustic dose distributions in the thalamus were more accurate with our approach than with the pseudo-CT approach of both Burgos et al. and Izquierdo et al, resulting in near-optimal treatment planning and dose estimation at all frequencies compared to CT (reference).

Significance

Our DL approach porosity estimates with ~7% error, is robust to input image noise and MRI acquisition parameters (sequence, coils, field strength) and yields near-optimal treatment planning and dose estimates for both central (thalamus) and lateral brain targets (amygdala) in the 200-1000 kHz frequency range.

Model-based navigation of transcranial focused ultrasound neuromodulation in humans: Application to targeting the amygdala and thalamus

BACKGROUND

Transcranial focused ultrasound (tFUS) neuromodulation has shown promise in animals but is challenging to translate to humans because of the thicker skull that heavily scatters ultrasound waves.

OBJECTIVE

We develop and disseminate a model-based navigation (MBN) tool for acoustic dose delivery in the presence of skull aberrations that is easy to use by non-specialists.

METHODS

We pre-compute acoustic beams for thousands of virtual transducer locations on the scalp of the subject under study. We use the hybrid angular spectrum solver mSOUND, which runs in ∼4 s per solve per CPU yielding pre-computation times under 1 h for scalp meshes with up to 4000 faces and a parallelization factor of 5. We combine this pre-computed set of beam solutions with optical tracking, thus allowing real-time display of the tFUS beam as the operator freely navigates the transducer around the subject' scalp. We assess the impact of MBN versus line-of-sight targeting (LOST) positioning in simulations of 13 subjects.

RESULTS

Our navigation tool has a display refresh rate of ∼10 Hz. In our simulations, MBN increased the acoustic dose in the thalamus and amygdala by 8-67 % compared to LOST and avoided complete target misses that affected 10-20 % of LOST cases. MBN also yielded a lower variability of the deposited dose across subjects than LOST.

CONCLUSIONS

MBN may yield greater and more consistent (less variable) ultrasound dose deposition than transducer placement with line-of-sight targeting, and thus could become a helpful tool to improve the efficacy of tFUS neuromodulation.

A novel transcranial MR Guided focused ultrasound method based on the ultrashort echo time skull acoustic model and phase retrieval techniques

Transcranial ultrasound stimulation (TUS) has been clinically applied as a neuromodulation tool. Particularly, the phase array ultrasound can be applied in TUS to non-invasively focus on the cortex or deep brain. However, the vital phase distortion of the ultrasound induced by the skull limits its clinical application. In the current study, we aimed to develop a hybrid method that combines the ultrashort echo time (UTE) magnetic resonance imaging (MRI) sequences with the prDeep technique to achieve focusing ventral intermediate thalamic nucleus (VIM). The time-reversal (TR) approach of the UTE numerical acoustic model of the skull combined with the prDeep algorithm was used to reduce the number of iterations. The skull acoustic model simulation therapy process was establish to valid this method's prediction and focus performance, and the classical TR method were considered as the gold standard (GS). Our approach could restore 75% of the GS intensity in 25 iteration steps, with a superior the noise immunity. Our findings demonstrate that the phase aberration caused by the skull can be estimated using phase retrieval techniques to achieve a fast and accurate transcranial focus. The method has excellent adaptability and anti-noise capacity for satisfying complex and changeable scenarios.

Ultrasound robotics for precision therapy

In recent years, the application of microrobots in precision therapy has gained significant attention. The small size and maneuverability of these micromachines enable them to potentially access regions that are difficult to reach using traditional methods; thus, reducing off-target toxicities and maximizing treatment effectiveness. Specifically, acoustic actuation has emerged as a promising method to exert control. By harnessing the power of acoustic energy, these small machines potentially navigate the body, assemble at the desired sites, and deliver therapies with enhanced precision and effectiveness. Amidst the enthusiasm surrounding these miniature agents, their translation to clinical environments has proven difficult. The primary objectives of this review are threefold: firstly, to offer an overview of the fundamental acoustic principles employed in the field of microrobots; secondly, to assess their current applications in medical therapies, encompassing tissue targeting, drug delivery or even cell infiltration; and lastly, to delve into the continuous efforts aimed at integrating acoustic microrobots into in vivo applications.

A rapid element pressure field simulation method for transcranial phase correction in focused ultrasound therapy

Transcranial focused ultrasound ablation has emerged as a promising technique for treating neurological disorders. The clinical system exclusively employed the ray tracing method to compute phase aberrations induced by the human skull, taking into account computational time constraints. However, this method compromises slightly on accuracy compared to simulation-based methods. This study evaluates a fast simulation method that simulates the time-harmonic pressure field within the region of interest for effective phase correction. Experimental validation was carried out using a 512-element, 670 kHz hemispherical transducer for fourex vivoskulls. The ray tracing method achieved a restoration ratio of 64.81% ± 4.33% of acoustic intensity normalized to hydrophone measurements. In comparison, the rapid simulation method demonstrated improved results with a restoration ratio of 73.10% ± 7.46%, albeit slightly lower than the full-wave simulation which achieved a restoration ratio of 75.87% ± 5.40%. The rapid simulation methods exhibited computational times that were less than five minutes for parallel computation with 8 threads. The incident angle was calculated, and a maximum difference of 6.8 degrees was found when the fixed position of the skull was changed. Meanwhile, the restoration ratio of acoustic intensity was validated to be above 70% for different target positions away from the geometrical focus of the transducer. The favorable balance between time consumption and correction accuracy makes this method valuable for clinical treatment applications.

Aberration correction in abdominal histotripsy

In transabdominal histotripsy, ultrasound pulses are focused on the body to noninvasively destroy soft tissues via cavitation. However, the ability to focus is limited by phase aberration, or decorrelation of the ultrasound pulses due to spatial variation in the speed of sound throughout heterogeneous tissue. Phase aberration shifts, broadens, and weakens the focus, thereby reducing the safety and efficacy of histotripsy therapy. This paper reviews and discusses aberration effects in histotripsy and in related therapeutic ultrasound techniques (e.g., high intensity focused ultrasound), with an emphasis on aberration by soft tissues. Methods for aberration correction are reviewed and can be classified into two groups: model-based methods, which use segmented images of the tissue as input to an acoustic propagation model to predict and compensate phase differences, and signal-based methods, which use a receive-capable therapy array to detect phase differences by sensing acoustic signals backpropagating from the focus. The relative advantages and disadvantages of both groups of methods are discussed. Importantly, model-based methods can correct focal shift, while signal-based methods can restore substantial focal pressure, suggesting that both methods should be combined in a 2-step approach. Aberration correction will be critical to improving histotripsy treatments and expanding the histotripsy treatment envelope to enable non-invasive, non-thermal histotripsy therapy for more patients.

Binary acoustic metasurfaces for dynamic focusing of transcranial ultrasound

Transcranial focused ultrasound (tFUS) is a promising technique for non-invasive and spatially targeted neuromodulation and treatment of brain diseases. Acoustic lenses were designed to correct the skull-induced beam aberration, but these designs could only generate static focused ultrasound beams inside the brain. Here, we designed and 3D printed binary acoustic metasurfaces (BAMs) for skull aberration correction and dynamic ultrasound beam focusing. BAMs were designed by binarizing the phase distribution at the surface of the metasurfaces. The phase distribution was calculated based on time reversal to correct the skull-induced phase aberration. The binarization enabled the ultrasound beam to be dynamically steered along wave propagation direction by adjusting the operation frequency of the incident ultrasound wave. The designed BAMs were manufactured by 3D printing with two coding bits, a polylactic acid unit for bit “1” and a water unit for bit “0.” BAMs for single- and multi-point focusing through the human skull were designed, 3D printed, and validated numerically and experimentally. The proposed BAMs with subwavelength scale in thickness are simple to design, easy to fabric, and capable of correcting skull aberration and achieving dynamic beam steering.

Comparison between MR and CT imaging used to correct for skull-induced phase aberrations during transcranial focused ultrasound

Transcranial focused ultrasound with the InSightec Exablate system uses thermal ablation for the treatment of movement and mood disorders and blood brain barrier disruption for tumor therapy. The system uses computed tomography (CT) images to calculate phase corrections that account for aberrations caused by the human skull. This work investigates whether magnetic resonance (MR) images can be used as an alternative to CT images to calculate phase corrections. Phase corrections were calculated using the gold standard hydrophone method and the standard of care InSightec ray tracing method. MR binary image mask, MR-simulated-CT (MRsimCT), and CT images of three ex vivo human skulls were supplied as inputs to the InSightec ray tracing method. The degassed ex vivo human skulls were sonicated with a 670 kHz hemispherical phased array transducer (InSightec Exablate 4000). 3D raster scans of the beam profiles were acquired using a hydrophone mounted on a 3-axis positioner system. Focal spots were evaluated using six metrics: pressure at the target, peak pressure, intensity at the target, peak intensity, positioning error, and focal spot volume. Targets at the geometric focus and 5 mm lateral to the geometric focus were investigated. There was no statistical difference between any of the metrics at either target using either MRsimCT or CT for phase aberration correction. As opposed to the MRsimCT, the use of CT images for aberration correction requires registration to the treatment day MR images; CT misregistration within a range of ± 2 degrees of rotation error along three dimensions was shown to reduce focal spot intensity by up to 9.4%. MRsimCT images used for phase aberration correction for the skull produce similar results as CT-based correction, while avoiding both CT to MR registration errors and unnecessary patient exposure to ionizing radiation.

Two-step aberration correction: application to transcranial histotripsy

Objective: Phase aberration correction is essential in transcranial histotripsy to compensate for focal distortion caused by the heterogeneity of the intact skull bone. This paper improves the 2-step aberration correction (AC) method that has been previously presented and develops an AC workflow that fits in the clinical environment, in which the computed tomography (CT)-based analytical approach was first implemented, followed by a cavitation-based approach using the shockwaves from the acoustic cavitation emission (ACE).Approach:A 700 kHz, 360-element hemispherical transducer array capable of transmit-and-receive on all channels was used to transcranially generate histotripsy-induced cavitation and acquire ACE shockwaves. For CT-AC, two ray-tracing models were investigated: a forward ray-tracing model (transducer-to-focus) in the open-source software Kranion, and an in-house backward ray-tracing model (focus-to-transducer) accounting for refraction and the sound speed variation in skulls. Co-registration was achieved by aligning the skull CT data to the skull surface map reconstructed using the acoustic pulse-echo method. For ACE-AC, the ACE signals from the collapses of generated bubbles were aligned by cross-correlation to estimate the corresponding time delays.Main results:The performance of the 2-step method was tested with 3 excised human calvariums placed at 2 different locations in the transducer array. Results showed that the 2-step AC achieved 90 ± 7% peak focal pressure compared to the gold standard hydrophone correction. It also reduced the focal shift from 0.84 to 0.30 mm and the focal volume from 10.6 to 2.0 mm3on average compared to the no AC cases.Significance:The 2-step AC yielded better refocusing compared to either CT-AC or ACE-AC alone and can be implemented in real-time for transcranial histotripsy brain therapy.

Differential evolution method to find optimal location of a single-element transducer for transcranial focused ultrasound therapy

BACKGROUND AND OBJECTIVE

Focused ultrasound (FUS) has been receiving growing attention as a noninvasive brain stimulation tool because of its superior spatial specificity and depth penetrability. However, the large mismatch of acoustic properties between the skull and water can disrupt and shift the acoustic focus in the brain. In this paper, we present a numerical method to find the optimal location of a single-element FUS transducer, which creates focus on the target region.

METHODS

The score function, representing the superposition of acoustic waves according to the relative phase difference and transmissibility, was defined based on time-reversal invariance of acoustic waves and depending on the spatial location of the transducer. The optimal location of the transducer was then determined using a differential evolution algorithm. To assess the proposed method, we conducted a forward simulation and compared the resulting focal location to the desired target point. We also performed experimental validation by measuring the acoustic pressure field through an ex vivo human skull in a water tank.

RESULTS

The numerical results indicated that the score function had a positive proportional relationship with the acoustic pressure at the target. Moreover, for the optimized transducer location, both the numerical and experimental results showed that the normalized acoustic pressure at the target was higher than 0.9.

CONCLUSIONS

In this study, we developed an optimization method to place a single-element transducer that effectively transmits acoustic energy to the targeted region in the brain. Our numerical and experimental results demonstrate that the proposed method can provide an optimal transducer location for safe and efficient FUS treatment.

Spatio-Spectral Ultrasound Characterization of Reflection and Transmission Through Bone With Temperature Dependence

Transcranial focused ultrasound (tFUS) is a promising approach for the treatment of neurological disorders. It has proven useful in several clinical applications, with promising outcomes reported in the recent literature. Furthermore, it is currently being investigated in a range of neuromodulation (NM) and ablative applications, including epilepsy. In this application, tFUS access through the temporal window is the key to optimizing the treatment safety and efficacy. Traditional approaches have utilized transducers with low operating frequencies for tFUS applications. Modern array transducers and driving systems allow for more intelligent use of the temporal window by exploiting the spatio-spectral transmission bandwidth to a specified target or targets within the brain. To demonstrate the feasibility of this approach, we have investigated the ultrasound reflection and transmission characteristics for different access points within the temporal window of human skull samples ex vivo. Different transmit-receive (Rx) configurations are used for characterization of the spatio-spectral variability in reflection and transmission through the temporal window. In this article, we show results from a dual-piston transducer set up in the frequency range of 2-7 MHz. Broadband pulses as well as synthesized orthogonal frequency division multiplexed (OFDM) waveforms were used. The latter was used to improve the magnitude and phase measurements in 100-kHz subbands within the 2-7 MHz spectral window. A temperature-controlled water bath was used to characterize the change in reflection and transmission characteristics with temperature in the 25°C-43°C range. The measured values of the complex reflection and transmission coefficients exhibited significant variations with space, frequency, and temperature. On the other hand, the measured transmission phase varied more with location and frequency, with smaller sensitivity to temperature. A measurement-based hybrid angular spectrum (HAS) simulation through the human temporal bone was used to demonstrate the dependence of focusing gain on the skull profile and spatial distribution of change of speed of sound (SOS) at different skull temperatures.

Towards controlled drug delivery in brain tumors with microbubble-enhanced focused ultrasound

Brain tumors are particularly challenging malignancies, due to their location in a structurally and functionally distinct part of the human body - the central nervous system (CNS). The CNS is separated and protected by a unique system of brain and blood vessel cells which together prevent most bloodborne therapeutics from entering the brain tumor microenvironment (TME). Recently, great strides have been made through microbubble (MB) ultrasound contrast agents in conjunction with ultrasound energy to locally increase the permeability of brain vessels and modulate the brain TME. As we elaborate in this review, this physical method can effectively deliver a wide range of anticancer agents, including chemotherapeutics, antibodies, and nanoparticle drug conjugates across a range of preclinical brain tumors, including high grade glioma (glioblastoma), diffuse intrinsic pontine gliomas, and brain metastasis. Moreover, recent evidence suggests that this technology can promote the effective delivery of novel immunotherapeutic agents, including immune check-point inhibitors and chimeric antigen receptor T cells, among others. With early clinical studies demonstrating safety, and several Phase I/II trials testing the preclinical findings underway, this technology is making firm steps towards shaping the future treatments of primary and metastatic brain cancer. By elaborating on its key components, including ultrasound systems and MB technology, along with methods for closed-loop spatial and temporal control of MB activity, we highlight how this technology can be tuned to enable new, personalized treatment strategies for primary brain malignancies and brain metastases.

Precise Ultrasound Neuromodulation in a Deep Brain Region Using Nano Gas Vesicles as Actuators

Ultrasound is a promising new modality for non-invasive neuromodulation. Applied transcranially, it can be focused down to the millimeter or centimeter range. The ability to improve the treatment's spatial resolution to a targeted brain region could help to improve its effectiveness, depending upon the application. The present paper details a neurostimulation scheme using gas-filled nanostructures, gas vesicles (GVs), as actuators for improving the efficacy and precision of ultrasound stimuli. Sonicated primary neurons display dose-dependent, repeatable Ca²⁺ responses, closely synced to stimuli, and increased nuclear expression of the activation marker c-Fos in the presence of GVs. GV-mediated ultrasound triggered rapid and reversible Ca²⁺ responses in vivo and could selectively evoke neuronal activation in a deep-seated brain region. Further investigation indicate that mechanosensitive ion channels are important mediators of this effect. GVs themselves and the treatment scheme are also found not to induce significant cytotoxicity, apoptosis, or membrane poration in treated cells. Altogether, this study demonstrates a simple and effective method to achieve enhanced and better-targeted neurostimulation with non-invasive low-intensity ultrasound.

Modeling Focused-Ultrasound Response for Non-Invasive Treatment Using Machine Learning

The interactions between body tissues and a focused ultrasound beam can be evaluated using various numerical models. Among these, the Rayleigh-Sommerfeld and angular spectrum methods are considered to be the most effective in terms of accuracy. However, they are computationally expensive, which is one of the underlying issues of most computational models. Typically, evaluations using these models require a significant amount of time (hours to days) if realistic scenarios such as tissue inhomogeneity or non-linearity are considered. This study aims to address this issue by developing a rapid estimation model for ultrasound therapy using a machine learning algorithm. Several machine learning models were trained on a very-large dataset (19,227 simulations), and the performance of these models were evaluated with metrics such as Root Mean Squared Error (RMSE), R-squared (R²), Akaike Information Criterion (AIC), and Bayesian Information Criterion (BIC). The resulted random forest provides superior accuracy with an R² value of 0.997, an RMSE of 0.0123, an AIC of -82.56, and a BIC of -81.65 on an external test dataset. The results indicate the efficacy of the random forest-based model for the focused ultrasound response, and practical adoption of this approach will improve the therapeutic planning process by minimizing simulation time.