Longitudinal and shear mode ultrasound propagation in human skull bone

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.

A review of functional neuromodulation in humans using low-intensity transcranial focused ultrasound

Transcranial ultrasonic neuromodulation is a rapidly burgeoning field where low-intensity transcranial focused ultrasound (tFUS), with exquisite spatial resolution and deep tissue penetration, is used to non-invasively activate or suppress neural activity in specific brain regions. Over the past decade, there has been a rapid increase of tFUS neuromodulation studies in healthy humans and subjects with central nervous system (CNS) disease conditions, including a recent surge of clinical investigations in patients. This narrative review summarized the findings of human neuromodulation studies using either tFUS or unfocused transcranial ultrasound (TUS) reported from 2013 to 2023. The studies were categorized into two separate sections: healthy human research and clinical studies. A total of 42 healthy human investigations were reviewed as grouped by targeted brain regions, including various cortical, subcortical, and deep brain areas including the thalamus. For clinical research, a total of 22 articles were reviewed for each studied CNS disease condition, including chronic pain, disorder of consciousness, Alzheimer's disease, Parkinson's disease, depression, schizophrenia, anxiety disorders, substance use disorder, drug-resistant epilepsy, and stroke. Detailed information on subjects/cohorts, target brain regions, sonication parameters, outcome readouts, and stimulatory efficacies were tabulated for each study. In later sections, considerations for planning tFUS neuromodulation in humans were also concisely discussed. With an excellent safety profile to date, the rapid growth of human tFUS research underscores the increasing interest and recognition of its significant potential in the field of non-invasive brain stimulation (NIBS), offering theranostic potential for neurological and psychiatric disease conditions and neuroscientific tools for functional brain mapping.

Numerical Investigation of Layered Homogeneous Skull Model for Simulations of Transcranial Focused Ultrasound

BACKGROUND AND OBJECTIVES

The influence of the intracranial pressure field must be discussed with the development of a single-element transducer for low-intensity transcranial focused ultrasound because the skull plays a significant role in blocking and dispersing ultrasound wave propagation. Ultrasound propagation is mainly affected by the structure and acoustic properties of the skull; thus, we aimed to investigate the impact of simplifying the acoustic properties of the skull on the simulation of the transcranial pressure field to present guidance for efficient skull modeling in full-wave simulations.

MATERIALS AND METHODS

We constructed a three-dimensional computational model for ultrasound transmission with the same structure but varying acoustic properties of the skull. The structural information and heterogeneous acoustic properties of the skull were acquired from computed tomography images, and we segmented the skull into three layers (3 L), including spongy and compact bones. We then assigned homogeneous acoustic properties to a single layer (1 L) or 3 L of the skull. In addition, we investigated the influence of different types of transducers and different ultrasound frequencies (1.1 MHz, 0.5 MHz, and 0.25 MHz) on the intracranial pressure field to provide a comparison of the heterogenous and homogeneous models.

RESULTS

We indicated the importance of numerical simulations in estimating the intracranial pressure field of the skull owing to beam distortions. When we simplified the skull model, both the 1 L and 3 L models showed contours of the acoustic focus comparable to those of the heterogeneous model. When we evaluated the peak pressure and volume of the acoustic focus, the 1 L model produced a better estimation of peak pressure with a difference <10%, and the 3 L model is suitable to obtain smaller errors in the volume of the acoustic focus.

CONCLUSIONS

In conclusion, we examined the possibility of simplification of skull models using 1 L and 3 L homogeneous properties in the numerical simulation for focused ultrasound. The results show that the layered homogeneous model can provide characteristics comparable to those of the acoustic focus in heterogeneous models.

Ultrasound pulse repetition frequency preferentially activates different neuron populations independent of cell type

Transcranial ultrasound activates mechanosensitive cellular signaling and modulates neural dynamics. Given that intrinsic neuronal activity is limited to a couple hundred hertz and often exhibits frequency preference, we examined whether pulsing ultrasound at physiologic pulse repetition frequencies (PRFs) could selectively influence neuronal activity in the mammalian brain. We performed calcium imaging of individual motor cortex neurons, while delivering 0.35 MHz ultrasound at PRFs of 10, 40, and 140 Hz in awake mice. We found that most neurons were preferentially activated by only one of the three PRFs, highlighting unique cellular effects of physiologic PRFs. Further, ultrasound evoked responses were similar between excitatory neurons and parvalbumin positive interneurons regardless of PRFs, indicating that individual cell sensitivity dominates ultrasound-evoked effects, consistent with the heterogeneous mechanosensitive channel expression we found across single neurons in mice and humans. These results highlight the feasibility of tuning ultrasound neuromodulation effects through varying PRFs.

Device for Multifocal Delivery of Ultrasound Into Deep Brain Regions in Humans

Low-intensity focused ultrasound provides the means to noninvasively stimulate or release drugs in specified deep brain targets. However, successful clinical translations require hardware that maximizes acoustic transmission through the skull, enables flexible electronic steering, and provides accurate and reproducible targeting while minimizing the use of MRI. We have developed a device that addresses these practical requirements. The device delivers ultrasound through the temporal and parietal skull windows, which minimize the attenuation and distortions of the ultrasound by the skull. The device consists of 252 independently controlled elements, which provides the ability to modulate multiple deep brain targets at a high spatiotemporal resolution, without the need to move the device or the subject. And finally, the device uses a mechanical registration method that enables accurate deep brain targeting both inside and outside of the MRI. Using this method, a single MRI scan is necessary for accurate targeting; repeated subsequent treatments can be performed reproducibly in an MRI-free manner. We validated these functions by transiently modulating specific deep brain regions in two patients with treatment-resistant depression.

Intracranial Gene Delivery Mediated by Albumin-Based Nanobubbles and Low-Frequency Ultrasound

Research in the field of high-intensity focused ultrasound (HIFU) for intracranial gene therapy has greatly progressed over the years. However, limitations of conventional HIFU still remain. That is, genes are required to cross the blood-brain barrier (BBB) in order to reach the neurological disordered lesion. In this study, we introduce a novel direct intracranial gene delivery method, bypassing the BBB using human serum albumin-based nanobubbles (NBs) injected through a less invasive intrathecal route via lumbar puncture, followed by intracranial irradiation with low-frequency ultrasound (LoFreqUS). Focusing on both plasmid DNA (pDNA) and messenger RNA (mRNA), our approach utilizes LoFreqUS for deeper tissue acoustic penetration and enhancing gene transfer efficiency. This drug delivery method could be dubbed as the "Spinal Back-Door Approach", an alternative to the "front door" BBB opening method. Experiments showed that NBs effectively responded to LoFreqUS, significantly improving gene transfer in vitro using U-87 MG cell lines. In vivo experiments in mice demonstrated significantly increased gene expression with pDNA; however, we were unable to obtain conclusive results using mRNA. This novel technique, combining albumin-based NBs and LoFreqUS offers a promising, efficient, targeted, and non-invasive solution for central nervous system gene therapy, potentially transforming the treatment landscape for neurological disorders.

Transcranial Phase Correction Using Pulse-Echo Ultrasound and Deep Learning: A 2-D Numerical Study

Phase aberration caused by human skulls severely degrades the quality of transcranial ultrasound images, posing a major challenge in the practical application of transcranial ultrasound techniques in adults. Aberration can be corrected if the skull profile (i.e., thickness distribution) and speed of sound (SOS) are known. However, accurately estimating the skull profile and SOS using ultrasound with a physics-based approach is challenging due to the complexity of the interaction between ultrasound and the skull. A deep learning approach is proposed herein to estimate the skull profile and SOS using ultrasound radiofrequency (RF) signals backscattered from the skull. A numerical study was performed to test the approach's feasibility. Realistic numerical skull models were constructed from computed tomography (CT) scans of five ex vivo human skulls in this numerical study. Acoustic simulations were performed on 3595 skull segments to generate array-based ultrasound backscattered signals. A deep learning model was developed and trained to estimate skull thickness and SOS from RF channel data. The trained model was shown to be highly accurate. The mean absolute error (MAE) was 0.15 mm (2% error) for thickness estimation and 13 m/s (0.5% error) for SOS estimation. The Pearson correlation coefficient between the estimated and ground-truth values was 0.99 for thickness and 0.95 for SOS. Aberration correction performed using deep-learning-estimated skull thickness and SOS values yielded significantly improved beam focusing (e.g., narrower beams) and transcranial imaging quality (e.g., improved spatial resolution and reduced artifacts) compared with no aberration correction. The results demonstrate the feasibility of the proposed approach for transcranial phase aberration correction.

Quantitative photoacoustic tomography: modeling and inverse problems

Significance

Quantitative photoacoustic tomography (QPAT) exploits the photoacoustic effect with the aim of estimating images of clinically relevant quantities related to the tissue's optical absorption. The technique has two aspects: an acoustic part, where the initial acoustic pressure distribution is estimated from measured photoacoustic time-series, and an optical part, where the distributions of the optical parameters are estimated from the initial pressure.

Aim

Our study is focused on the optical part. In particular, computational modeling of light propagation (forward problem) and numerical solution methodologies of the image reconstruction (inverse problem) are discussed.

Approach

The commonly used mathematical models of how light and sound propagate in biological tissue are reviewed. A short overview of how the acoustic inverse problem is usually treated is given. The optical inverse problem and methods for its solution are reviewed. In addition, some limitations of real-life measurements and their effect on the inverse problems are discussed.

Results

An overview of QPAT with a focus on the optical part was given. Computational modeling and inverse problems of QPAT were addressed, and some key challenges were discussed. Furthermore, the developments for tackling these problems were reviewed. Although modeling of light transport is well-understood and there is a well-developed framework of inverse mathematics for approaching the inverse problem of QPAT, there are still challenges in taking these methodologies to practice.

Conclusions

Modeling and inverse problems of QPAT together were discussed. The scope was limited to the optical part, and the acoustic aspects were discussed only to the extent that they relate to the optical aspect.

Ultrasonic therapies for seizures and drug-resistant epilepsy

Ultrasonic therapy is an increasingly promising approach for the treatment of seizures and drug-resistant epilepsy (DRE). Therapeutic focused ultrasound (FUS) uses thermal or nonthermal energy to either ablate neural tissue or modulate neural activity through high- or low-intensity FUS (HIFU, LIFU), respectively. Both HIFU and LIFU approaches have been investigated for reducing seizure activity in DRE, and additional FUS applications include disrupting the blood-brain barrier in the presence of microbubbles for targeted-drug delivery to the seizure foci. Here, we review the preclinical and clinical studies that have used FUS to treat seizures. Additionally, we review effective FUS parameters and consider limitations and future directions of FUS with respect to the treatment of DRE. While detailed studies to optimize FUS applications are ongoing, FUS has established itself as a potential noninvasive alternative for the treatment of DRE and other neurological disorders.

Transcranial Ultrasonic Focusing by a Phased Array Based on Micro-CT Images

In this paper, we utilize micro-computed tomography (micro-CT) to obtain micro-CT images with a resolution of 60 μm and establish a micro-CT model based on the k-wave toolbox, which can visualize the microstructures in trabecular bone, including pores and bone layers. The transcranial ultrasound phased array focusing field characteristics in the micro-CT model are investigated. The ultrasonic waves are multiply scattered in skull and time delays calculations from the transducer to the focusing point are difficult. For this reason, we adopt the pulse compression method and the linear frequency modulation Barker code to compute the time delay and implement phased array focusing in the micro-CT model. It is shown by the simulation results that ultrasonic loss is mainly caused by scattering from the microstructures of the trabecular bone. The ratio of main and side lobes of the cross-correlation calculation is improved by 5.53 dB using the pulse compression method. The focusing quality and the calculation accuracy of time delay are improved. Meanwhile, the beamwidth at the focal point and the sound pressure amplitude decrease with the increase in the signal frequency. Focusing at different depths indicates that the beamwidth broadens with the increase in the focusing depth, and beam deflection focusing maintains good consistency in the focusing effect at a distance of 9 mm from the focal point. This indicates that the phased-array method has good focusing results and focus tunability in deep cranial brain. In addition, the sound pressure at the focal point can be increased by 8.2% through amplitude regulation, thereby enhancing focusing efficiency. The preliminary experiment verification is conducted with an ex vivo skull. It is shown by the experimental results that the phased array focusing method using pulse compression to calculate the time delay can significantly improve the sound field focusing effect and is a very effective transcranial ultrasound focusing method.

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.

Evaluation of ultrasound sensors for transcranial photoacoustic sensing and imaging

Photoacoustic imaging through skull bone causes strong attenuation and distortion of the acoustic wavefront, which diminishes image contrast and resolution. As a result, transcranial photoacoustic measurements in humans have been challenging to demonstrate. In this study, we investigated the acoustic transmission through the human skull to design an ultrasound sensor suitable for transcranial PA imaging and sensing. We measured the frequency dependent losses of human cranial bones ex vivo, compared the performance of a range of piezoelectric and optical ultrasound sensors, and imaged skull phantoms using a PA tomograph based on a planar Fabry-Perot sensor. All transcranial photoacoustic measurements show the typical effects of frequency and thickness dependent attenuation and aberration associated with acoustic propagation through bone. The performance of plano-concave optical resonator ultrasound sensors was found to be highly suitable for transcranial photoacoustic measurements.

State of the art on microbubble cavitation monitoring and feedback control for blood-brain-barrier opening using focused ultrasound

Focused ultrasound (FUS) is a non-invasive and highly promising method for targeted and reversible blood-brain barrier permeabilization. Numerous preclinical studies aim to optimize the localized delivery of drugs using this method in rodents and non-human primates. Several clinical trials have been initiated to treat various brain diseases in humans using simultaneous BBB permeabilization and drug injection. This review presents the state of the art ofin vitroandin vivocavitation control algorithms for BBB permeabilization using microbubbles (MB) and FUS. Firstly, we describe the different cavitation states, their physical significance in terms of MB behavior and their translation into the spectral composition of the backscattered signal. Next, we report the different indexes calculated and used during the ultrasonic monitoring of cavitation. Finally, the differentin vitroandin vivocavitation control strategies described in the literature are presented and compared.

A simulation study on the sensitivity of transcranial ray-tracing ultrasound modeling to skull properties

In transcranial focused ultrasound therapies, such as treating essential tremor via thermal ablation in the thalamus, acoustic energy is focused through the skull using a phased-array transducer. Ray tracing is a computationally efficient method that can correct skull-induced phase aberrations via per-element phase delay calculations using patient-specific computed tomography (CT) data. However, recent studies show that variations in CT-derived Hounsfield unit may account for only 50% of the speed of sound variability in human skull specimens, potentially limiting clinical transcranial ultrasound applications. Therefore, understanding the sensitivity of treatment planning methods to material parameter variations is essential. The present work uses a ray-tracing simulation model to explore how imprecision in model inputs, arising from clinically significant uncertainties in skull properties or considerations of acoustic phenomena, affects acoustic focusing quality through the skull. We propose and validate new methods to optimize ray-tracing skull simulations for clinical treatment planning, relevant for predicting intracranial target's thermal rise, using experimental data from ex-vivo human skulls.

Current and Emerging Systems for Focused Ultrasound-Mediated Blood-Brain Barrier Opening

With an ever-growing list of neurological applications of focused ultrasound (FUS), there has been a consequent increase in the variety of systems for delivering ultrasound energy to the brain. Specifically, recent successful pilot clinical trials of blood-brain barrier (BBB) opening with FUS have generated substantial interest in the future applications of this relatively novel therapy, with divergent, purpose-built technologies emerging. With many of these technologies at various stages of pre-clinical and clinical investigation, this article seeks to provide an overview and analysis of the numerous medical devices in active use and under development for FUS-mediated BBB opening.

Photoacoustic digital brain and deep-learning-assisted image reconstruction

Photoacoustic tomography (PAT) is a newly developed medical imaging modality, which combines the advantages of pure optical imaging and ultrasound imaging, owning both high optical contrast and deep penetration depth. Very recently, PAT is studied in human brain imaging. Nevertheless, while ultrasound waves are passing through the human skull tissues, the strong acoustic attenuation and aberration will happen, which causes photoacoustic signals' distortion. In this work, we use 180 T1 weighted magnetic resonance imaging (MRI) human brain volumes along with the corresponding magnetic resonance angiography (MRA) brain volumes, and segment them to generate the 2D human brain numerical phantoms for PAT. The numerical phantoms contain six kinds of tissues, which are scalp, skull, white matter, gray matter, blood vessel and cerebrospinal fluid. For every numerical phantom, Monte-Carlo based optical simulation is deployed to obtain the photoacoustic initial pressure based on optical properties of human brain. Then, two different k-wave models are used for the skull-involved acoustic simulation, which are fluid media model and viscoelastic media model. The former one only considers longitudinal wave propagation, and the latter model takes shear wave into consideration. Then, the PA sinograms with skull-induced aberration is taken as the input of U-net, and the skull-stripped ones are regarded as the supervision of U-net to train the network. Experimental result shows that the skull's acoustic aberration can be effectively alleviated after U-net correction, achieving conspicuous improvement in quality of PAT human brain images reconstructed from the corrected PA signals, which can clearly show the cerebral artery distribution inside the human skull.

Physics-informed neural networks for transcranial ultrasound wave propagation

Transcranial ultrasound imaging has been playing an increasingly important role in the non-invasive treatment of brain disorders. However, the conventional mesh-based numerical wave solvers, which are an integral part of imaging algorithms, suffer from limitations such as high computational cost and discretization error in predicting the wavefield passing through the skull. In this paper, we explore the use of physics-informed neural networks (PINNs) for predicting the transcranial ultrasound wave propagation. The wave equation, two sets of time snapshots data and a boundary condition (BC) are embedded as physical constraints in the loss function during training. The proposed approach has been validated by solving the two-dimensional (2D) acoustic wave equation under three increasingly complex spatially varying velocity models. Our cases demonstrate that due to the meshless nature of PINNs, they can be flexibly applied to different wave equations and types of BCs. By adding physics constraints to the loss function, PINNs can predict wavefields far outside the training data, providing ideas for improving the generalization capability of existing deep learning methods. The proposed approach offers exciting perspectives because of the powerful framework and simple implementation. We conclude with a summary of the strengths, limitations and further research directions of this work.

Ultrasound as a versatile tool for short- and long-term improvement and monitoring of brain function

Treating the brain with focused ultrasound (FUS) at low intensities elicits diverse responses in neurons, astroglia, and the extracellular matrix. In combination with intravenously injected microbubbles, FUS also opens the blood-brain barrier (BBB) and facilitates focal drug delivery. However, an incompletely understood cellular specificity and a wide parameter space currently limit the optimal application of FUS in preclinical and human studies. In this perspective, we discuss how different FUS modalities can be utilized to achieve short- and long-term improvements, thereby potentially treating brain disorders. We review the ongoing efforts to determine which parameters induce neuronal inhibition versus activation and how mechanoreceptors and signaling cascades are activated to induce long-term changes, including memory improvements. We suggest that optimal FUS treatments may require different FUS modalities and devices, depending on the targeted brain area or local pathology, and will be greatly enhanced by new techniques for monitoring FUS efficacy.

Effect of skull porosity on ultrasound transmission and wave mode conversion at large incidence angles

BACKGROUND

Transcranial ultrasound imaging and therapy depend on the efficient transmission of acoustic energy through the skull. Multiple previous studies have concluded that a large incidence angle should be avoided during transcranial-focused ultrasound therapy to ensure transmission through the skull. Alternatively, some other studies have shown that longitudinal-to-shear wave mode conversion might improve transmission through the skull when the incidence angle is increased above the critical angle (i.e., 25° to 30°).

PURPOSE

The effect of skull porosity on the transmission of ultrasound through the skull at varying incidence angles was investigated for the first time to elucidate why transmission through the skull at large angles of incidence is decreased in some cases but improved in other cases.

METHODS

Transcranial ultrasound transmission at varying incidence angles (0°-50°) was investigated in phantoms and ex vivo skull samples with varying bone porosity (0% to 28.54% ± 3.36%) using both numerical and experimental methods. First, the elastic acoustic wave transmission through the skull was simulated using micro-computed tomography data of ex vivo skull samples. The trans-skull pressure was compared between skull segments having three levels of porosity, that is, low porosity (2.65% ± 0.03%), medium porosity (13.41% ± 0.12%), and high porosity (26.9%). Next, transmission through two 3D-printed resin skull phantoms (compact vs. porous phantoms) was experimentally measured to test the effect of porous microstructure alone on ultrasound transmission through flat plates. Finally, the effect of skull porosity on ultrasound transmission was investigated experimentally by comparing transmission through two ex vivo human skull segments having similar thicknesses but different porosities (13.78% ± 2.05% vs. 28.54% ± 3.36%).

RESULTS

Numerical simulations indicated that an increase in transmission pressure occurs at large incidence angles for skull segments having low porosities but not for those with high porosity. In experimental studies, a similar phenomenon was observed. Specifically, for the low porosity skull sample (13.78% ± 2.05%), the normalized pressure was 0.25 when the incidence angle increased to 35°. However, for the high porosity sample (28.54% ± 3.36%), the pressure was no more than 0.1 at large incidence angles.

CONCLUSIONS

These results indicate that the skull porosity has an evident effect on the transmission of ultrasound at large incidence angles. The wave mode conversion at large, oblique incidence angles could enhance the transmission of ultrasound through parts of the skull having lower porosity in the trabecular layer. However, for transcranial ultrasound therapy in the presence of highly porous trabecular bone, transmission at a normal incidence angle is preferable relative to oblique incidence angles due to the higher transmission efficiency.

Optimizing transcranial ultrasound delivery at large incident angles by leveraging cranial leaky guided wave dispersion

We investigate the role of leaky guided waves in transcranial ultrasound transmission in temporal and parietal bones at large incidence angles. Our numerical and experimental results show that the dispersion characteristics of the fundamental leaky guided wave mode with longitudinal polarization can be leveraged to estimate the critical angle above which efficient shear mode conversion takes place, and below which major transmission drops can be expected. Simulations that employ a numerical propagator matrix and a Semi-Analytical approach establish the transcranial dispersion characteristics and transmission coefficients at different incident angles. Experimental transmission tests conducted at 500 kHz and radiation tests performed in the 200-800 kHz range confirm the numerical findings in terms of transmitted peak pressure and frequency-radiation angle spectra, based on which the connection between critical angles, dispersion and transmission is demonstrated. Our results support the identification of transcranial ultrasound strategies that leverage shear mode conversion, which is less sensitive to phase aberrations compared to normal incidence ultrasound. These findings can also enable higher transmission rates in cranial bones with low porosity by leveraging dispersion information extracted through signal processing, without requiring measurement of geometric and mechanical properties of the cranial bone.

Transcranial ultrasound simulations: A review

Transcranial ultrasound is more and more used for therapy and imaging of the brain. However, the skull is a highly attenuating and aberrating medium, with different structures and acoustic properties among samples and even within a sample. Thus, case-specific simulations are needed to perform transcranial focused ultrasound interventions safely. In this article, we provide a review of the different methods used to model the skull and to simulate ultrasound propagation through it.

Numerical and experimental evaluation of low-intensity transcranial focused ultrasound wave propagation using human skulls for brain neuromodulation

BACKGROUND

Low-intensity transcranial focused ultrasound (tFUS) has gained considerable attention as a promising noninvasive neuromodulatory technique for human brains. However, the complex morphology of the skull hinders scholars from precisely predicting the acoustic energy transmitted and the region of the brain impacted during the sonication. This is due to the fact that different ultrasound frequencies and skull morphology variations greatly affect wave propagation through the skull.

PURPOSE

Although the acoustic properties of human skull have been studied for tFUS applications, such as tumor ablation using a multielement phased array, there is no consensus about how to choose a single-element focused ultrasound (FUS) transducer with a suitable frequency for neuromodulation. There are interests in exploring the magnitude and dimension of tFUS beam through human parietal bone for modulating specific brain lobes. Herein, we aim to investigate the wave propagation of tFUS on human skulls to understand and address the concerns above.

METHODS

Both experimental measurements and numerical modeling were conducted to investigate the transmission efficiency and beam pattern of tFUS on five human skulls (C3 and C4 regions) using single-element FUS transducers with six different frequencies (150-1500 kHz). The degassed skull was placed in a water tank, and a calibrated hydrophone was utilized to measure acoustic pressure past it. The cranial computed tomography scan data of each skull were obtained to derive a high-resolution acoustic model (grid point spacing: 0.25 mm) in simulations. Meanwhile, we modified the power-law exponent of acoustic attenuation coefficient to validate numerical modeling and enabled it to be served as a prediction tool, based on the experimental measurements.

RESULTS

The transmission efficiency and -6 dB beamwidth were evaluated and compared for various frequencies. An exponential decrease in transmission efficiency and a logarithmic decrease of -6 dB beamwidth with an increase in ultrasound frequency were observed. It is found that a >750 kHz ultrasound leads to a relatively lower tFUS transmission efficiency (<5%), whereas a <350 kHz ultrasound contributes to a relatively broader beamwidth (>5 mm). Based on these observations, we further analyzed the dependence of tFUS wave propagation on FUS transducer aperture size.

CONCLUSIONS

We successfully studied tFUS wave propagation through human skulls at different frequencies experimentally and numerically. The findings have important implications to predict tFUS wave propagation for ultrasound neuromodulation in clinical applications, and guide researchers to develop advanced ultrasound transducers as neural interfaces.

Three-layer model with absorption for conservative estimation of the maximum acoustic transmission coefficient through the human skull for transcranial ultrasound stimulation

Transcranial ultrasound stimulation (TUS) has been shown to be a safe and effective technique for non-invasive superficial and deep brain stimulation. Safe and efficient translation to humans requires estimating the acoustic attenuation of the human skull. Nevertheless, there are no international guidelines for estimating the impact of the skull bone. A tissue independent, arbitrary derating was developed by the U.S. Food and Drug Administration to take into account tissue absorption (0.3 dB/cm-MHz) for diagnostic ultrasound. However, for the case of transcranial ultrasound imaging, the FDA model does not take into account the insertion loss induced by the skull bone, nor the absorption by brain tissue. Therefore, the estimated absorption is overly conservative which could potentially limit TUS applications if the same guidelines were to be adopted. Here we propose a three-layer model including bone absorption to calculate the maximum pressure transmission through the human skull for frequencies ranging between 100 kHz and 1.5 MHz. The calculated pressure transmission decreases with the frequency and the thickness of the bone, with peaks for each thickness corresponding to a multiple of half the wavelength. The 95th percentile maximum transmission was calculated over the accessible surface of 20 human skulls for 12 typical diameters of the ultrasound beam on the skull surface, and varies between 40% and 78%. To facilitate the safe adjustment of the acoustic pressure for short ultrasound pulses, such as transcranial imaging or transcranial ultrasound stimulation, a table summarizes the maximum pressure transmission for each ultrasound beam diameter and each frequency.

Sound out the impaired perfusion: Photoacoustic imaging in preclinical ischemic stroke

Acoustically detecting the optical absorption contrast, photoacoustic imaging (PAI) is a highly versatile imaging modality that can provide anatomical, functional, molecular, and metabolic information of biological tissues. PAI is highly scalable and can probe the same biological process at various length scales ranging from single cells (microscopic) to the whole organ (macroscopic). Using hemoglobin as the endogenous contrast, PAI is capable of label-free imaging of blood vessels in the brain and mapping hemodynamic functions such as blood oxygenation and blood flow. These imaging merits make PAI a great tool for studying ischemic stroke, particularly for probing into hemodynamic changes and impaired cerebral blood perfusion as a consequence of stroke. In this narrative review, we aim to summarize the scientific progresses in the past decade by using PAI to monitor cerebral blood vessel impairment and restoration after ischemic stroke, mostly in the preclinical setting. We also outline and discuss the major technological barriers and challenges that need to be overcome so that PAI can play a more significant role in preclinical stroke research, and more importantly, accelerate its translation to be a useful clinical diagnosis and management tool for human strokes.

A tool for monitoring cell type-specific focused ultrasound neuromodulation and control of chronic epilepsy

Focused ultrasound (FUS) is a powerful tool for noninvasive modulation of deep brain activity with promising therapeutic potential for refractory epilepsy; however, tools for examining FUS effects on specific cell types within the deep brain do not yet exist. Consequently, how cell types within heterogeneous networks can be modulated and whether parameters can be identified to bias these networks in the context of complex behaviors remains unknown. To address this, we developed a fiber Photometry Coupled focused Ultrasound System (PhoCUS) for simultaneously monitoring FUS effects on neural activity of subcortical genetically targeted cell types in freely behaving animals. We identified a parameter set that selectively increases activity of parvalbumin interneurons while suppressing excitatory neurons in the hippocampus. A net inhibitory effect localized to the hippocampus was further confirmed through whole brain metabolic imaging. Finally, these inhibitory selective parameters achieved significant spike suppression in the kainate model of chronic temporal lobe epilepsy, opening the door for future noninvasive therapies.

S.M.A.R.T. F.U.S: Surrogate Model of Attenuation and Refraction in Transcranial Focused Ultrasound

Low-intensity focused ultrasound (LIFU) is an increasingly applied method for achieving non-invasive brain stimulation. However, transmission of ultrasound through the human skull can substantially affect focal point characteristics of LIFU, including dramatic attenuation in intensity and refraction of focal point location. These effects depend on a high-dimensional parameter space, making these effects difficult to estimate from previous work. Instead, focal point properties of LIFU experiments are often estimated using numerical simulation of LIFU sonication through skull. However, this procedure presents many entry barriers to even computationally savvy investigators and often requires expensive computational hardware, impeding LIFU research. We present a novel MATLAB toolbox (data: doi:10.5068/D1QD60; Matlab Scripts: https://doi.org/10.5281/zenodo.5811122) for rapidly estimating beam properties of LIFU transmitted through bone. Users provide specific values for frequency of LIFU, bone thickness, angle at which LIFU is applied, depth of the LIFU focal point, and diameter of the transducer used and receive an estimation of the degree of refraction/attenuation expected for the given parameters.

Design and Construction of a Low-Frequency Ultrasound Acquisition Device for 2-D Brain Imaging Using Full-Waveform Inversion

The main techniques used to image the brain and obtain structural data are magnetic resonance imaging and X-ray computed tomography. These techniques produce images with high spatial resolution, but with the disadvantage of requiring very large equipment with special installation needs. In addition, X-ray tomography uses ionizing radiation, which limits their use. Ultrasound imaging is a safe technology that is delivered using compact and mobile devices. However, conventional ultrasound reconstruction techniques have failed to obtain images of the brain because of, fundamentally, the presence of the skull and the distortion that it produces on ultrasound. Recent studies have indicated that full-waveform inversion, a computational technique originally from Earth science, has the potential to generate accurate 3-D images of the brain. This technology can overcome the limitations of conventional ultrasound imaging, but a prototype for transcranial applications does not yet exist. Here, we investigate different designs of an annular array of ultrasound transducers to optimize the number of elements and rotations needed to conduct transcranial imaging with full-waveform inversion. This device uses small-diameter, low-frequency transducers that readily propagate ultrasound through the skull with good signal-to-noise ratios. It also incorporates the use of rotations to produce a high-density coverage of the target and acquire redundant traces that are beneficial for full-waveform inversion. We have built a ring of 40 transducers to illustrate that this design is capable of reconstructing images of the brain, retrieving its anatomy and acoustic properties with millimeter resolution. Laboratory results reveal the ability of this device to successfully image a 2.5-D brain- and skull-mimicking phantom using full-waveform inversion. To our knowledge, this is the first prototype ever used for transcranial-like imaging. The importance of these findings and their implications for the design of a 3-D reconstruction system with possible clinical applications are discussed.

Low-intensity focused ultrasound attenuates early traumatic brain injury by OX-A/NF-κB/NLRP3 signaling pathway

Background: Traumatic brain injury (TBI) is a serious hazard to human health and is characterized by high rates of disability and mortality. It is necessary to explore new effective treatment methods to reduce the impact of TBI on individuals and society. As an emerging neuromodulation technique, ultrasound is used to treat some neurological diseases, but the neuroprotective mechanism of low-intensity focused ultrasound (LIFUS) in TBI remains unclear. We aimed to investigate the protective effects and potential mechanisms of LIFUS in TBI.

Methods: A rat model of TBI was established using the free-fall method. After establishing the TBI model, the hypothalamus region was covered with LIFUS radiation, and an orexin receptor 1 (OXR1) antagonist (SB334867) was injected intraperitoneally. Neurobehavioral examination, Nissl staining, hematoxylin and eosin staining of the brain tissue, and brain water content, were performed 3 days later. Western blotting, quantitative real-time polymerase chain reaction, immunofluorescence staining, and immunohistochemical staining, were used to evaluate the neuroprotective mechanisms of LIFUS.

Results: LIFUS improved tissue damage, neurological deficits, and brain edema. LIFUS can increase the expression of orexin-A (OX-A) and OXR1, significantly inhibit the activation of nuclear factor-κB (NF-κB) protein and nucleotide-binding domain-like receptor protein 3 (NLRP3) inflammasome after TBI, and reduce the release of pro-inflammatory factors after TBI; however, SB334867 can reverse this effect.

Conclusions: This study suggests that LIFUS may play a neuroprotective role by promoting the release of OX-A from the hypothalamus and inhibiting the inflammatory response after TBI through the OX-A /NF-κB/NLRP3 pathway.

Benchmark problems for transcranial ultrasound simulation: Intercomparison of compressional wave models

Computational models of acoustic wave propagation are frequently used in transcranial ultrasound therapy, for example, to calculate the intracranial pressure field or to calculate phase delays to correct for skull distortions. To allow intercomparison between the different modeling tools and techniques used by the community, an international working group was convened to formulate a set of numerical benchmarks. Here, these benchmarks are presented, along with intercomparison results. Nine different benchmarks of increasing geometric complexity are defined. These include a single-layer planar bone immersed in water, a multi-layer bone, and a whole skull. Two transducer configurations are considered (a focused bowl and a plane piston operating at 500 kHz), giving a total of 18 permutations of the benchmarks. Eleven different modeling tools are used to compute the benchmark results. The models span a wide range of numerical techniques, including the finite-difference time-domain method, angular spectrum method, pseudospectral method, boundary-element method, and spectral-element method. Good agreement is found between the models, particularly for the position, size, and magnitude of the acoustic focus within the skull. When comparing results for each model with every other model in a cross-comparison, the median values for each benchmark for the difference in focal pressure and position are less than 10% and 1 mm, respectively. The benchmark definitions, model results, and intercomparison codes are freely available to facilitate further comparisons.

Benchmark problems for transcranial ultrasound simulation: Intercomparison of compressional wave models

Computational models of acoustic wave propagation are frequently used in transcranial ultrasound therapy, for example, to calculate the intracranial pressure field or to calculate phase delays to correct for skull distortions. To allow intercomparison between the different modeling tools and techniques used by the community, an international working group was convened to formulate a set of numerical benchmarks. Here, these benchmarks are presented, along with intercomparison results. Nine different benchmarks of increasing geometric complexity are defined. These include a single-layer planar bone immersed in water, a multi-layer bone, and a whole skull. Two transducer configurations are considered (a focused bowl and a plane piston operating at 500 kHz), giving a total of 18 permutations of the benchmarks. Eleven different modeling tools are used to compute the benchmark results. The models span a wide range of numerical techniques, including the finite-difference time-domain method, angular spectrum method, pseudospectral method, boundary-element method, and spectral-element method. Good agreement is found between the models, particularly for the position, size, and magnitude of the acoustic focus within the skull. When comparing results for each model with every other model in a cross-comparison, the median values for each benchmark for the difference in focal pressure and position are less than 10% and 1 mm, respectively. The benchmark definitions, model results, and intercomparison codes are freely available to facilitate further comparisons.

Clinical Intervention Using Focused Ultrasound (FUS) Stimulation of the Brain in Diverse Neurological Disorders

Various surgical techniques and pharmaceutical treatments have been developed to improve the current technologies of treating brain diseases. Focused ultrasound (FUS) is a new brain stimulation modality that can exert a therapeutic effect on diseased brain cells, with this effect ranging from permanent ablation of the pathological neural circuit to transient excitatory/inhibitory modulation of the neural activity depending on the acoustic energy of choice. With the development of intraoperative imaging technology, FUS has become a clinically available noninvasive neurosurgical option with visual feedback. Over the past 10 years, FUS has shown enormous potential. It can deliver acoustic energy through the physical barrier of the brain and eliminate abnormal brain cells to treat patients with Parkinson's disease and essential tremor. In addition, FUS can help introduce potentially beneficial therapeutics at the exact brain region where they need to be, bypassing the brain's function barrier, which can be applied for a wide range of central nervous system disorders. In this review, we introduce the current FDA-approved clinical applications of FUS, ranging from thermal ablation to blood barrier opening, as well as the emerging applications of FUS in the context of pain control, epilepsy, and neuromodulation. We also discuss the expansion of future applications and challenges. Broadening FUS technologies requires a deep understanding of the effect of ultrasound when targeting various brain structures in diverse disease conditions in the context of skull interface, anatomical structure inside the brain, and pathology.

An Ultrasound-Guided Hemispherical Phased Array for Microbubble-Mediated Ultrasound Therapy

GOAL

To develop a low-cost magnetic resonance imaging (MRI)-free transcranial focused ultrasound (FUS) system for microbubble-mediated therapy.

METHODS

A 128-element 11 MHz array for skull localization was integrated within a 256-module multi-frequency (306/612/1224 kHz) dual-mode phased array. The system's transcranial transmit and receive performance was evaluated with ex-vivo human skullcaps using phase aberration corrections calculated from computed tomography (CT)-based simulations via ultrasound-based (USCT) and landmark-based (LMCT) registrations, and a gold-standard fixed source emitter (FSE)-based method.

RESULTS

Displacement and rotation registration errors of 1.4 ± 0.4 mm and 2.1 ± 0.2 ^° were obtained using USCT, resulting in sub-millimeter transmit targeting errors driven at 306 kHz (0.9 ± 0.2 mm) and 612 kHz (0.9 ± 0.3 mm), and source localization errors of 1.0 ± 0.3 mm and 0.6 ± 0.2 mm at receive frequencies of 306 kHz and 612 kHz, respectively (mean ± SD). Similar errors were obtained using LMCT and no significant differences between these two approaches were found on either transmit (p = 0.64/0.99) or receive (p = 0.45/0.36) at 306 kHz/612kHz. During volumetric multi-point exposures, approximately 70% and 60% of the transmit frames in which microbubble activity was detected via FSE were recovered using USCT when imaging at the second-harmonic and half-harmonic, respectively, compared to 60% and 69% using LMCT.

CONCLUSION

This low-cost ultrasound-guided transcranial FUS system affords USCT skull registration with accuracy comparable to LMCT methods.

SIGNIFICANCE

Such systems have great potential to advance the adoption of microbubble-mediated FUS brain therapy by improving access to the technology.

Refraction-Corrected Transcranial Ultrasound Imaging Through the Human Temporal Window Using a Single Probe

Transcranial ultrasound imaging (TUI) is a diagnostic modality with numerous applications, but unfortunately, it is hindered by phase aberration caused by the skull. In this article, we propose to reconstruct a transcranial B-mode image with a refraction-corrected synthetic aperture imaging (SAI) scheme. First, the compressional sound velocity of the aberrator (i.e., the skull) is estimated using the bidirectional headwave technique. The medium is described with four layers (i.e., lens, water, skull, and water), and a fast marching method calculates the travel times between individual array elements and image pixels. Finally, a delay-and-sum algorithm is used for image reconstruction with coherent compounding. The point spread function (PSF) in a wire phantom image and reconstructed with the conventional technique (using a constant sound speed throughout the medium), and the proposed method was quantified with numerical synthetic data and experiments with a bone-mimicking plate and a human skull, compared with the PSF achieved in a ground truth image of the medium without the aberrator (i.e., the bone plate or skull). A phased-array transducer (P4-1, ATL/Philips, 2.5 MHz, 96 elements, pitch = 0.295 mm) was used for the experiments. The results with the synthetic signals, the bone-mimicking plate, and the skull indicated that the proposed method reconstructs the scatterers with an average lateral/axial localization error of 0.06/0.14 mm, 0.11/0.13 mm, and 1.0/0.32 mm, respectively. With the human skull, an average contrast ratio (CR) and full-width-half-maximum (FWHM) of 37.1 dB and 1.75 mm were obtained with the proposed approach, respectively. This corresponds to an improvement of CR and FWHM by 7.1 dB and 36% compared with the conventional method, respectively. These numbers were 12.7 dB and 41% with the bone-mimicking plate.

Increasing the transmission efficiency of transcranial ultrasound using a dual-mode conversion technique based on Lamb waves

Transcranial focused ultrasound (FUS) is a noninvasive treatment for brain tumors and neuromodulation. Based on normal incidence, conventional FUS techniques use a focused or an array of ultrasonic transducers to overcome the attenuation and absorption of ultrasound in the skull; however, this remains the main limitation of using FUS. A dual-mode conversion technique based on Lamb waves is proposed to achieve high transmission efficiency. This concept was validated using the finite element analysis (FEA) and experiments based on changes in the incident angle. Aluminum, plexiglass, and a human skull were used as materials with different attenuations. The transmission loss was calculated for each material, and the results were compared with the reflectance function of the Lamb waves. Oblique incidence based on dual-mode conversion exhibited a better transmission efficiency than that of a normal incidence for all of the specimens. The total transmission losses for the materials were 13.7, 15.46, and 3.91 dB less than those associated with the normal incidence. A wedge transducer was designed and fabricated to implement the proposed method. The results demonstrated the potential applicability of the dual-mode conversion technique for the human skull.

A Spatial Multitarget Ultrasound Neuromodulation System Using High-Powered 2-D Array Transducer

Transcranial focused ultrasound (tFUS) is increasingly used in experimental neuroscience due to its neuromodulatory effectiveness in animal studies. However, achieving multitarget tFUS in small animals is typically limited by transducer size, energy transfer efficiency, and brain volume. The objective of this work was to construct an ultrasound system for multitarget neuromodulation in small animals. First, a miniaturized high-powered 2-D array transducer was developed. The phase delay of each array element was calculated based on the multifocal time-reversal method, generating multiple foci simultaneously in a 3-D field. The effects of the axial focal length, interfocus spacing (lateral distance between the two focal centers), and the number of foci on the focal properties of the pressure field were examined through numerical simulations. In-vitro ultrasonic measurements and transcranial simulations on a rat skull were conducted. The minimum interfocus spacing separating two -6-dB foci and the peak full-width at half-maximum were positively correlated with axial focal length; the relative relationship between the interfocus spacing and pressure field properties was similar for each axial focal length. The maximum acoustic pressure and spatial average intensity at focus in deionized water were 2.21 MPa and 133 W/cm², respectively. The simulated and experimental results were compared, demonstrating agreement in both peak position and focus shape. The ultrasound system can provide a neuroscientific platform for evaluating the feasibility of multitarget ultrasound stimulation treatment protocols, thus improving the understanding of functional neuroanatomy.

Establishing density-dependent longitudinal sound speed in the vertebral lamina

Focused ultrasound treatments of the spinal cord may be facilitated using a phased array transducer and beamforming to correct spine-induced focal aberrations. Simulations can non-invasively calculate aberration corrections using x-ray computed tomography (CT) data that are correlated to density (ρ) and longitudinal sound speed (c_L). We aimed to optimize vertebral lamina-specific c_L(ρ) functions at a physiological temperature (37 °C) to maximize time domain simulation accuracy. Odd-numbered ex vivo human thoracic vertebrae were imaged with a clinical CT-scanner (0.511 × 0.511 × 0.5 mm), then sonicated with a transducer (514 kHz) focused on the canal via the vertebral lamina. Vertebra-induced signal time shifts were extracted from pressure waveforms recorded within the canals. Measurements were repeated 5× per vertebra, with 2.5 mm vertical vertebra shifts between measurements. Linear functions relating c_L with CT-derived density were optimized. The optimized function was c_L(ρ)=0.35(ρ-ρ_w)+ c_L,w m/s, where w denotes water, giving the tested laminae a mean bulk density of 1600 ± 30 kg/m³ and a mean bulk c_L of 1670 ± 60 m/s. The optimized lamina c_L(ρ) function was accurate to λ/16 when implemented in a multi-layered ray acoustics model. This modelling accuracy will improve trans-spine ultrasound beamforming.

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.

Acoustic properties across the human skull

Transcranial ultrasound is emerging as a noninvasive tool for targeted treatments of brain disorders. Transcranial ultrasound has been used for remotely mediated surgeries, transient opening of the blood-brain barrier, local drug delivery, and neuromodulation. However, all applications have been limited by the severe attenuation and phase distortion of ultrasound by the skull. Here, we characterized the dependence of the aberrations on specific anatomical segments of the skull. In particular, we measured ultrasound propagation properties throughout the perimeter of intact human skulls at 500 kHz. We found that the parietal bone provides substantially higher transmission (average pressure transmission 31 ± 7%) and smaller phase distortion (242 ± 44 degrees) than frontal (13 ± 2%, 425 ± 47 degrees) and occipital bone regions (16 ± 4%, 416 ± 35 degrees). In addition, we found that across skull regions, transmission strongly anti-correlated (R=-0.79) and phase distortion correlated (R=0.85) with skull thickness. This information guides the design, positioning, and skull correction functionality of next-generation devices for effective, safe, and reproducible transcranial focused ultrasound therapies.

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.

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.

3-D Transcranial Microbubble Cavitation Localization by Four Sensors

Cavitation is the fundamental physical mechanism of various focused ultrasound (FUS)-mediated therapies in the brain. Accurately knowing the three-dimensional (3-D) location of cavitation in real-time can improve the targeting accuracy and avoid off-target tissue damage. Existing techniques for 3-D passive transcranial cavitation detection require the use of expensive and complicated hemispherical phased arrays with 128 or 256 elements. The objective of this study was to investigate the feasibility of using four sensors for transcranial 3-D localization of cavitation. Differential microbubble cavitation detection combined with the time difference of arrival algorithm was developed for the localization using the four sensors. Numerical simulation using k-Wave toolbox was performed to validate the proposed method for transcranial cavitation source localization. The sensors with a center frequency of 2.25 MHz and a 6 dB bandwidth of 1.39 MHz were used to locate cavitation generated by FUS (500 kHz) sonication of microbubbles that were injected into a tube positioned inside an ex vivo human skullcap. Cavitation emissions from the microbubbles were detected transcranially using the four sensors. Both simulation and experimental studies found that the proposed method achieved accurate 3-D cavitation localization. When the cavitation source was located within 30 mm from the geometric center of the sensor network, the accuracy of the localization method with the skull was measured to be 1.9±1.0 mm, which was not significantly different from that without the skull (1.7 ± 0.5 mm). The accuracy decreased as the cavitation source was away from the geometric center of the sensor network. It also decreased as the pulse length increased. Its accuracy was not significantly affected by the sensor position relative to the skull. In summary, four sensors combined with the proposed localization algorithm offer a simple approach for 3-D transcranial cavitation localization.

Verification of Blood-Brain Barrier Disruption Based on the Clinical Validation Platform Using a Rat Model with Human Skull

Methods to improve drug delivery efficiency through blood-brain barrier disruption (BBBD) based on microbubbles and focused ultrasound (FUS) are continuously being studied. However, most studies are being conducted in preclinical trial environments using small animals. The use of the human skull shows differences between the clinical and preclinical trials. BBBD results from preclinical trials are difficult to represent in clinical trials because various distortions of ultrasound by the human skull are excluded in the former. Therefore, in our study, a clinical validation platform based on a preclinical trial environment, using a human skull fragment and a rat model, was developed to induce BBBD under conditions similar to clinical trials. For this, a human skull fragment was inserted between the rat head and a 250 kHz FUS transducer, and optimal ultrasound parameters for the free field (without human skull fragment) and human skull (with human skull fragment) were derived by 300 mV_pp and 700 mV_pp, respectively. BBBD was analyzed according to each case using magnetic resonance images, Evans blue dye, cavitation, and histology. Although it was confirmed using magnetic resonance images and Evans blue dye that a BBB opening was induced in each case, multiple BBB openings were observed in the brain tissues. This phenomenon was analyzed by numerical simulation, and it was confirmed to be due to standing waves owing to the small skull size of the rat model. The stable cavitation doses (SCD_h and SCD_u) in the human skull decreased by 13.6- and 5.3-fold, respectively, compared to those in the free field. Additionally, the inertial cavitation dose in the human skull decreased by 1.05-fold compared to that of the free field. For the histological analysis, although some extravasated red blood cells were observed in each case, it was evaluated as recoverable based on our previous study results. Therefore, our proposed platform can help deduct optimal ultrasound parameters and BBBD results for clinical trials in the preclinical trials with small animals because it considers variables relevant to the human skull.

Transcranial focused ultrasound modulates cortical and thalamic motor activity in awake sheep

Transcranial application of pulsed low-intensity focused ultrasound (FUS) modulates the excitability of region-specific brain areas, and anesthetic confounders on brain activity warrant the evaluation of the technique in awake animals. We examined the neuromodulatory effects of FUS in unanesthetized sheep by developing a custom-fit headgear capable of reproducibly placing an acoustic focus on the unilateral motor cortex (M1) and corresponding thalamic area. The efferent responses to sonication, based on the acoustic parameters previously identified in anesthetized sheep, were measured using electromyography (EMG) from both hind limbs across three experimental conditions: on-target sonication, off-target sonication, and without sonication. Excitatory sonication yielded greater amplitude of EMG signals obtained from the hind limb contralateral to sonication than that from the ipsilateral limb. Spurious appearance of motion-related EMG signals limited the amount of analyzed data (~ 10% selection of acquired data) during excitatory sonication, and the averaged EMG response rates elicited by the M1 and thalamic stimulations were 7.5 ± 1.4% and 6.7 ± 1.5%, respectively. Suppressive sonication, while sheep walked on the treadmill, temporarily reduced the EMG amplitude from the limb contralateral to sonication. No significant change was found in the EMG amplitudes during the off-target sonication. Behavioral observation throughout the study and histological analysis showed no sign of brain tissue damage caused by the acoustic stimulation. Marginal response rates observed during excitatory sonication call for technical refinement to reduce motion artifacts during EMG acquisitions as well as acoustic aberration correction schemes to improve spatial accuracy of sonication. Yet, our results indicate that low-intensity FUS modulated the excitability of regional brain tissues reversibly and safely in awake sheep, supporting its potential in theragnostic applications.

Region-specific effects of ultrasound on individual neurons in the awake mammalian brain

Ultrasound modulates brain activity. However, it remains unclear how ultrasound affects individual neurons in the brain, where neural circuit architecture is intact and different brain regions exhibit distinct tissue properties. Using a high-resolution calcium imaging technique, we characterized the effect of ultrasound stimulation on thousands of individual neurons in the hippocampus and the motor cortex of awake mice. We found that brief 100-ms-long ultrasound pulses increase intracellular calcium in a large fraction of individual neurons in both brain regions. Ultrasound-evoked calcium response in hippocampal neurons exhibits a rapid onset with a latency shorter than 50 ms. The evoked response in the hippocampus is shorter in duration and smaller in magnitude than that in the motor cortex. These results demonstrate that noninvasive ultrasound stimulation transiently increases intracellular calcium in individual neurons in awake mice, and the evoked response profiles are brain region specific.

Effect of Skull Porous Trabecular Structure on Transcranial Ultrasound Imaging in the Presence of Elastic Wave Mode Conversion at Varying Incidence Angle

With the advancement of aberration correction techniques, transcranial ultrasound imaging has exhibited great potential in applications such as imaging neurological function and guiding therapeutic ultrasound. However, the feasibility of transcranial imaging varies among individuals because of the differences in skull acoustic properties. To better understand the fundamental mechanisms underlying the variation in imaging performance, the effect of the structure of the porous trabecular bone on transcranial imaging performance (i.e., target localization errors and resolution) was investigated for the first time through the use of elastic wave simulations and experiments. Simulation studies using high-resolution computed tomography data from ex vivo skull samples revealed that imaging at large incidence angles reduced the target localization error for skulls having low porosity; however, as skull porosity increased, large angles of incidence resulted in degradation of resolution and increased target localization errors. Experimental results indicate that imaging at normal incidence introduced a localization error of 1.85 ± 0.10 mm, while imaging at a large incidence angle (40°) resulted in an increased localization error of 6.54 ± 1.33 mm and caused a single point target to no longer appear as a single, coherent target in the resulting image, which is consistent with simulation results. This first investigation of the effects of skull microstructure on transcranial ultrasound imaging indicates that imaging performance is highly dependent on the porosity of the skull, particularly at non-normal angles of incidence.

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.

Photoacoustic computed tomography for functional human brain imaging [Invited]

The successes of magnetic resonance imaging and modern optical imaging of human brain function have stimulated the development of complementary modalities that offer molecular specificity, fine spatiotemporal resolution, and sufficient penetration simultaneously. By virtue of its rich optical contrast, acoustic resolution, and imaging depth far beyond the optical transport mean free path (∼1 mm in biological tissues), photoacoustic computed tomography (PACT) offers a promising complementary modality. In this article, PACT for functional human brain imaging is reviewed in its hardware, reconstruction algorithms, in vivo demonstration, and potential roadmap.

Comparison of compressional and elastic wave simulations for patient-specific planning prior to transcranial photoacoustic-guided neurosurgery

SIGNIFICANCE

Simulations have the potential to be a powerful tool when planning the placement of photoacoustic imaging system components for surgical guidance. While elastic simulations (which include both compressional and shear waves) are expected to more accurately represent the physical transcranial acoustic wave propagation process, these simulations are more time-consuming and memory-intensive than the compressional-wave-only simulations that our group previously used to identify optimal acoustic windows for transcranial photoacoustic imaging.

AIM

We present qualitative and quantitative comparisons of compressional and elastic wave simulations to determine which option is more suitable for preoperative surgical planning.

APPROACH

Compressional and elastic photoacoustic k-Wave simulations were performed based on a computed tomography volume of a human cadaver head. Photoacoustic sources were placed in the locations of the internal carotid arteries and likely positions of neurosurgical instrument tips. Transducers received signals from three previously identified optimal acoustic windows (i.e., the ocular, nasal, and temporal regions). Target detectability, image-based target size estimates, and target-to-instrument distances were measured using the generalized contrast-to-noise ratio (gCNR), resolution, and relative source distances, respectively, for each simulation method.

RESULTS

The gCNR was equivalent between compressional and elastic simulations. The areas of the -6  dB contours of point spread functions utilized to measure resolution differed by 0.33 to 3.35  mm2. Target-to-instrument distance measurements were within 1.24 mm of the true distances.

CONCLUSIONS

These results indicate that it is likely sufficient to utilize the less time-consuming, less memory-intensive compressional wave simulations for presurgical planning.