Numerical evaluation of the skull for human neuromodulation with transcranial focused ultrasound

Model-based navigation of transcranial focused ultrasound neuromodulation in humans: Application to targeting of 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 seconds per solve per CPU yielding pre-computation times under one hour for scalp meshes with up to 4,000 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 yields 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 may become a helpful tool to improve the efficacy of tFUS neuromodulation.

Low-intensity focused ultrasound to the posterior insula reduces temporal summation of pain

BACKGROUND

The insula and dorsal anterior cingulate cortex (dACC) are core brain regions involved in pain processing and central sensitization, a shared mechanism across various chronic pain conditions. Methods to modulate these regions may serve to reduce central sensitization, though it is unclear which target may be most efficacious for different measures of central sensitization.

OBJECTIVE/HYPOTHESIS

Investigate the effect of low-intensity focused ultrasound (LIFU) to the anterior insula (AI), posterior insula (PI), or dACC on conditioned pain modulation (CPM) and temporal summation of pain (TSP).

METHODS

N = 16 volunteers underwent TSP and CPM pain tasks pre/post a 10 min LIFU intervention to either the AI, PI, dACC or Sham stimulation. Pain ratings were collected pre/post LIFU.

RESULTS

Only LIFU to the PI significantly attenuated pain ratings during the TSP protocol. No effects were found for the CPM task for any of the LIFU targets. LIFU pressure modulated group means but did not affect overall group differences.

CONCLUSIONS

LIFU to the PI reduced temporal summation of pain. This may, in part, be due to dosing (pressure) of LIFU. Inhibition of the PI with LIFU may be a future potential therapy in chronic pain populations demonstrating central sensitization. The minimal effective dose of LIFU for efficacious neuromodulation will help to translate LIFU for therapeutic options.

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.

A systematic review of preclinical and clinical transcranial ultrasound neuromodulation and opportunities for functional connectomics

BACKGROUND

Low-intensity transcranial ultrasound has surged forward as a non-invasive and disruptive tool for neuromodulation with applications in basic neuroscience research and the treatment of neurological and psychiatric conditions.

OBJECTIVE

To provide a comprehensive overview and update of preclinical and clinical transcranial low intensity ultrasound for neuromodulation and emphasize the emerging role of functional brain mapping to guide, better understand, and predict responses.

METHODS

A systematic review was conducted by searching the Web of Science and Scopus databases for studies on transcranial ultrasound neuromodulation, both in humans and animals.

RESULTS

187 relevant studies were identified and reviewed, including 116 preclinical and 71 clinical reports with subjects belonging to diverse cohorts. Milestones of ultrasound neuromodulation are described within an overview of the broader landscape. General neural readouts and outcome measures are discussed, potential confounds are noted, and the emerging use of functional magnetic resonance imaging is highlighted.

CONCLUSION

Ultrasound neuromodulation has emerged as a powerful tool to study and treat a range of conditions and its combination with various neural readouts has significantly advanced this platform. In particular, the use of functional magnetic resonance imaging has yielded exciting inferences into ultrasound neuromodulation and has the potential to advance our understanding of brain function, neuromodulatory mechanisms, and ultimately clinical outcomes. It is anticipated that these preclinical and clinical trials are the first of many; that transcranial low intensity focused ultrasound, particularly in combination with functional magnetic resonance imaging, has the potential to enhance treatment for a spectrum of neurological conditions.

Transcranial focused ultrasound to the posterior cingulate cortex modulates default mode network and subjective experience: an fMRI pilot study

Background

Transcranial focused ultrasound (TFUS) is an emerging neuromodulation tool for temporarily altering brain activity and probing network functioning. The effects of TFUS on the default mode network (DMN) are unknown.

Objective

The study examined the effects of transcranial focused ultrasound (TFUS) on the functional connectivity of the default mode network (DMN), specifically by targeting the posterior cingulate cortex (PCC). Additionally, we investigated the subjective effects of TFUS on mood, mindfulness, and self-related processing.

Methods

The study employed a randomized, single-blind design involving 30 healthy subjects. Participants were randomly assigned to either the active TFUS group or the sham TFUS group. Resting-state functional magnetic resonance imaging (rs-fMRI) scans were conducted before and after the TFUS application. To measure subjective effects, the Toronto Mindfulness Scale, the Visual Analog Mood Scale, and the Amsterdam Resting State Questionnaire were administered at baseline and 30 min after sonication. The Self Scale and an unstructured interview were also administered 30 min after sonication.

Results

The active TFUS group exhibited significant reductions in functional connectivity along the midline of the DMN, while the sham TFUS group showed no changes. The active TFUS group demonstrated increased state mindfulness, reduced Global Vigor, and temporary alterations in the sense of ego, sense of time, and recollection of memories. The sham TFUS group showed an increase in state mindfulness, too, with no other subjective effects.

Conclusions

TFUS targeted at the PCC can alter DMN connectivity and cause changes in subjective experience. These findings support the potential of TFUS to serve both as a research tool and as a potential therapeutic intervention.

Effects of the motor cortical theta-burst transcranial-focused ultrasound stimulation on the contralateral motor cortex

Theta-burst transcranial ultrasound stimulation (tbTUS) increases primary motor cortex (M1) excitability for at least 30 min. However, the remote effects of focal M1 tbTUS on the excitability of other cortical areas are unknown. Here, we examined the effects of left M1 tbTUS on right M1 excitability. An 80 s train of active or sham tbTUS was delivered to the left M1 in 20 healthy subjects. Before and after the tbTUS, we measured: (1) corticospinal excitability using motor-evoked potential (MEP) amplitudes from single-pulse transcranial magnetic stimulation (TMS) of left and right M1; (2) interhemispheric inhibition (IHI) from left to right M1 and from right to left M1 using a dual-site paired-pulse TMS paradigm; and (3) intracortical circuits of the right M1 with short-interval intracortical inhibition and intracortical facilitation (ICF) using paired-pulse TMS. Left M1 tbTUS decreased right M1 excitability as shown by decreased MEP amplitudes, increased right M1 ICF and decreased short-interval IHI from left to right hemisphere at interstimulus interval (ISI) of 10 ms but not long-interval IHI at interstimulus interval of 40 ms. The study showed that left M1 tbTUS can change the excitability of remote cortical areas with decreased right M1 excitability and interhemispheric inhibition. The remote effects of tbTUS should be considered when it is used in neuroscience research and as a potential neuromodulation treatment for brain disorders. KEY POINTS: Transcranial ultrasound stimulation (TUS) is a novel non-invasive brain stimulation technique for neuromodulation with the advantages of being able to achieve high spatial resolution and target deep brain structures. A repetitive TUS protocol, with an 80 s train of theta burst patterned TUS (tbTUS), has been shown to increase primary motor cortex (M1) excitability, as well as increase alpha and beta movement-related spectral power in distinct brain regions. In this study, we examined on the effects of the motor cortical tbTUS on the excitability of contralateral M1 measured with MEPs elicited by transcranial magnetic stimulation. We showed that left M1 tbTUS decreased right M1 excitability and left-to-right M1 interhemispheric inhibition, and increased intracortical facilitation of right M1. These results lead to better understand the effects of tbTUS and can help the development of tbTUS for the treatment of neurological and psychiatric disorders and in neuroscience research.

Plasticity-Induced Effects of Theta Burst Transcranial Ultrasound Stimulation in Parkinson's Disease

BACKGROUND

Low-intensity transcranial ultrasound stimulation (TUS) is a noninvasive brain stimulation (NIBS) technique with high spatial specificity. Previous studies showed that TUS delivered in a theta burst pattern (tbTUS) increased motor cortex (MI) excitability up to 30 minutes due to long-term potentiation (LTP)-like plasticity. Studies using other forms of NIBS suggested that cortical plasticity may be impaired in patients with Parkinson's disease (PD).

OBJECTIVE

The aim was to investigate the neurophysiological effects of tbTUS in PD patients off and on dopaminergic medications compared to healthy controls.

METHODS

We studied 20 moderately affected PD patients in on and off dopaminergic medication states (7 with and 13 without dyskinesia) and 17 age-matched healthy controls in a case-controlled study. tbTUS was applied for 80 seconds to the MI. Motor-evoked potentials (MEP), short-interval intracortical inhibition (SICI), and short-interval intracortical facilitation (SICF) were recorded at baseline, and at 5 minutes (T5), T30, and T60 after tbTUS. Motor Unified Parkinson's Disease Rating Scale (mUPDRS) was measured at baseline and T60.

RESULTS

tbTUS significantly increased MEP amplitude at T30 compared to baseline in controls and in PD patients on but not in PD patients off medications. SICI was reduced in PD off medications compared to controls. tbTUS did not change in SICI or SICF. The bradykinesia subscore of mUPDRS was reduced at T60 compared to baseline in PD on but not in the off medication state. The presence of dyskinesia did not affect tbTUS-induced plasticity.

CONCLUSIONS

tbTUS-induced LTP plasticity is impaired in PD patients off medications and is restored by dopaminergic medications. © 2024 The Author(s). Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Noninvasive intervention by transcranial ultrasound stimulation: Modulation of neural circuits and its clinical perspectives

Low-intensity focused transcranial ultrasound stimulation (TUS) is an emerging noninvasive technique capable of stimulating both the cerebral cortex and deep brain structures with high spatial precision. This method is recognized for its potential to comprehensively perturb various brain regions, enabling the modulation of neural circuits, in a manner not achievable through conventional magnetic or electrical brain stimulation techniques. The underlying mechanisms of neuromodulation are based on a phenomenon where mechanical waves of ultrasound kinetically interact with neurons, specifically affecting neuronal membranes and mechanosensitive channels. This interaction induces alterations in the excitability of neurons within the stimulated region. In this review, we briefly present the fundamental principles of ultrasound physics and the physiological mechanisms of TUS neuromodulation. We explain the experimental apparatus and procedures for TUS in humans. Due to the focality, the integration of various methods, including magnetic resonance imaging and magnetic resonance-guided neuronavigation systems, is important to perform TUS experiments for precise targeting. We then review the current state of the literature on TUS neuromodulation, with a particular focus on human subjects, targeting both the cerebral cortex and deep subcortical structures. Finally, we outline future perspectives of TUS in clinical applications in psychiatric and neurological fields.

Translating ultrasound-mediated drug delivery technologies for CNS applications

Ultrasound enhances drug delivery into the central nervous system (CNS) by opening barriers between the blood and CNS and by triggering release of drugs from carriers. A key challenge in translating setups from in vitro to in vivo settings is achieving equivalent acoustic energy delivery. Multiple devices have now been demonstrated to focus ultrasound to the brain, with concepts emerging to also target the spinal cord. Clinical trials to date have used ultrasound to facilitate the opening of the blood-brain barrier. While most have focused on feasibility and safety considerations, therapeutic benefits are beginning to emerge. To advance translation of these technologies for CNS applications, researchers should standardise exposure protocol and fine-tune ultrasound parameters. Computational modelling should be increasingly used as a core component to develop both in vitro and in vivo setups for delivering accurate and reproducible ultrasound to the CNS. This field holds promise for transformative advancements in the management and pharmacological treatment of complex and challenging CNS disorders.

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.

Low-intensity focused ultrasound to the human insular cortex differentially modulates the heartbeat-evoked potential: a proof-of-concept study

Background

The heartbeat evoked potential (HEP) is a brain response time-locked to the heartbeat and a potential marker of interoceptive processing. The insula and dorsal anterior cingulate cortex (dACC) are brain regions that may be involved in generating the HEP. Low-intensity focused ultrasound (LIFU) is a non-invasive neuromodulation technique that can selectively target sub-regions of the insula and dACC to better understand their contributions to the HEP.

Objective

Proof-of-concept study to determine whether LIFU modulation of the anterior insula (AI), posterior insula (PI), and dACC influences the HEP.

Methods

In a within-subject, repeated-measures design, healthy human participants (n=16) received 10 minutes of stereotaxically targeted LIFU to the AI, PI, dACC or Sham at rest during continuous electroencephalography (EEG) and electrocardiography (ECG) recording on separate days. Primary outcome was change in HEP amplitudes. Relationships between LIFU pressure and HEP changes were examined using linear mixed modelling. Peripheral indices of visceromotor output including heart rate and heart rate variability (HRV) were explored between conditions.

Results

Relative to sham, LIFU to the PI, but not AI or dACC, decreased HEP amplitudes; this was partially explained by increased LIFU pressure. LIFU did not affect time or frequency dependent measures of HRV.

Conclusions

These results demonstrate the ability to modulate HEP amplitudes via non-invasive targeting of key interoceptive brain regions. Our findings have implications for the causal role of these areas in bottom-up heart-brain communication that could guide future work investigating the HEP as a marker of interoceptive processing in healthy and clinical populations.

Transcranial pulse stimulation in Alzheimer's disease

BACKGROUND

Transcranial pulse stimulation (TPS) is a novel noninvasive ultrasonic brain stimulation that can increase cortical and corticospinal excitability, induce neuroplasticity, and increase functional connectivity within the brain. Several trials have confirmed its potential in treating Alzheimer's disease (AD).

OBJECTIVE

To investigate the effect and safety of TPS on AD.

DESIGN

A systematic review.

METHODS

PubMed, Embase via Ovid, Web of Science, Cochrane Library, CNKI (China National Knowledge Infrastructure), VIP (China Science and Technology Journal Database), and WanFang were searched from inception to April 1, 2023. Study selection, data extraction, and quality evaluation of the studies were conducted by two reviewers independently, with any controversy resolved by consensus. The Methodological Index for Nonrandomized Studies was used to assess the risk of bias.

RESULTS

Five studies were included in this review, with a total of 99 patients with AD. For cognitive performance, TPS significantly improved the scores of the CERAD (Consortium to Establish a Registry for Alzheimer's Disease) test battery, Alzheimer's Disease Assessment Scale (cognitive), Montreal Cognitive Assessment, and Mini-Mental Status Examination. For depressive symptoms, TPS significantly reduced the scores of the Alzheimer's Disease Assessment Scale (affective), Geriatric Depression Score, and Beck Depression Inventory. By functional magnetic resonance imaging, studies have shown that TPS improved cognitive performance in AD patients by increasing functional connectivity in the hippocampus, parahippocampal cortex, precuneus, and parietal cortex, and activating cortical activity in the bilateral hippocampus. TPS alleviated depressive symptoms in AD patients by decreasing functional connectivity between the ventromedial network (left frontal orbital cortex) and the salience network (right anterior insula). Adverse events in this review, including headache, worsening mood, jaw pain, nausea, and drowsiness, were reversible and lasted no longer than 1 day. No serious adverse events or complications were observed.

CONCLUSIONS

TPS is promising in improving cognitive performance and reducing depressive symptoms in patients with AD. TPS may be a safe adjunct therapy in the treatment of AD. However, these findings lacked a sham control and were limited by the small sample size of the included studies. Further research may be needed to better explore the potential of TPS.

PATIENT AND PUBLIC INVOLVEMENT

Patients and the public were not involved in this study.

Improving Sonication Efficiency in Transcranial MR-Guided Focused Ultrasound Treatment: A Patient-Data Simulation Study

The potential improvement in sonication efficiency achieved by tilting the focused ultrasound (FUS) transducer of the transcranial MR-guided FUS system is presented. A total of 56 cases of patient treatment data were used. The relative position of the clinical FUS transducer to the patient's head was reconstructed, and region-specific skull density and porosity were calculated based on the patient's CT volume image. The total transmission coefficient of acoustic waves emitted from each channel was calculated. Then, the total energy penetrating the human skull-which represents the sonication efficiency-was estimated. As a result, improved sonication efficiency was by titling the FUS transducer to a more appropriate angle achieved in all 56 treatment cases. This simulation result suggests the potential improvement in transcranial-focused ultrasound treatment by simply adjusting the transducer angle.

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.

Real-Time Acoustic Simulation Framework for tFUS: A Feasibility Study Using Navigation System

Transcranial focused ultrasound (tFUS), in which acoustic energy is focused on a small region in the brain through the skull, is a non-invasive therapeutic method with high spatial resolution and depth penetration. Image-guided navigation has been widely utilized to visualize the location of acoustic focus in the cranial cavity. However, this system is often inaccurate because of the significant aberrations caused by the skull. Therefore, acoustic simulations using a numerical solver have been widely adopted to compensate for this inaccuracy. Although the simulation can predict the intracranial acoustic pressure field, real-time application during tFUS treatment is almost impossible due to the high computational cost. In this study, we propose a neural network-based real-time acoustic simulation framework and test its feasibility by implementing a simulation-guided navigation (SGN) system. Real-time acoustic simulation is performed using a 3D conditional generative adversarial network (3D-cGAN) model featuring residual blocks and multiple loss functions. This network was trained by the conventional numerical acoustic simulation program (i.e., k-Wave). The SGN system is then implemented by integrating real-time acoustic simulation with a conventional image-guided navigation system. The proposed system can provide simulation results with a frame rate of 5 Hz (i.e., about 0.2 s), including all processing times. In numerical validation (3D-cGAN vs. k-Wave), the average peak intracranial pressure error was 6.8 ± 5.5%, and the average acoustic focus position error was 5.3 ± 7.7 mm. In experimental validation using a skull phantom (3D-cGAN vs. actual measurement), the average peak intracranial pressure error was 4.5%, and the average acoustic focus position error was 6.6 mm. These results demonstrate that the SGN system can predict the intracranial acoustic field according to transducer placement in real-time.

A transducer positioning method for transcranial focused ultrasound treatment of brain tumors

Purpose

As a non-invasive method for brain diseases, transcranial focused ultrasound (tFUS) offers higher spatial precision and regulation depth. Due to the altered path and intensity of sonication penetrating the skull, the focus and intensity in the skull are difficult to determine, making the use of ultrasound therapy for cancer treatment experimental and not widely available. The deficiency can be effectively addressed by numerical simulation methods, which enable the optimization of sonication modulation parameters and the determination of precise transducer positioning.

Methods

A 3D skull model was established using binarized brain CT images. The selection of the transducer matrix was performed using the radius positioning (RP) method after identifying the intracranial target region. Simulations were performed, encompassing acoustic pressure (AP), acoustic field, and temperature field, in order to provide compelling evidence of the safety of tFUS in sonication-induced thermal effects.

Results

It was found that the angle of sonication path to the coronal plane obtained at all precision and frequency models did not exceed 10° and 15° to the transverse plane. The results of thermal effects illustrated that the peak temperatures of tFUS were 43.73°C, which did not reach the point of tissue degeneration. Once positioned, tFUS effectively delivers a Full Width at Half Maximum (FWHM) stimulation that targets tumors with diameters of up to 3.72 mm in a one-off. The original precision model showed an attenuation of 24.47 ± 6.13 mm in length and 2.40 ± 1.42 mm in width for the FWHM of sonication after penetrating the skull.

Conclusion

The vector angles of the sonication path in each direction were determined based on the transducer positioning results. It has been suggested that when time is limited for precise transducer positioning, fixing the transducer on the horizontal surface of the target region can also yield positive results for stimulation. This framework used a new transducer localization method to offer a reliable basis for further research and offered new methods for the use of tFUS in brain tumor-related research.

Generating Patient-Specific Acoustic Simulations for Transcranial Focused Ultrasound Procedures Based on Optical Tracking Information

Optical tracking is a real-time transducer positioning method for transcranial focused ultrasound (tFUS) procedures, but the predicted focus from optical tracking typically does not incorporate subject-specific skull information. Acoustic simulations can estimate the pressure field when propagating through the cranium but rely on accurately replicating the positioning of the transducer and skull in a simulated space. Here, we develop and characterize the accuracy of a workflow that creates simulation grids based on optical tracking information in a neuronavigated phantom with and without transmission through an ex vivo skull cap. The software pipeline could replicate the geometry of the tFUS procedure within the limits of the optical tracking system (transcranial target registration error (TRE): 3.9 ± 0.7 mm). The simulated focus and the free-field focus predicted by optical tracking had low Euclidean distance errors of 0.5±0.1 and 1.2±0.4 mm for phantom and skull cap, respectively, and some skull-specific effects were captured by the simulation. However, the TRE of simulation informed by optical tracking was 4.6±0.2, which is as large or greater than the focal spot size used by many tFUS systems. By updating the position of the transducer using the original TRE offset, we reduced the simulated TRE to 1.1 ± 0.4 mm. Our study describes a software pipeline for treatment planning, evaluates its accuracy, and demonstrates an approach using MR-acoustic radiation force imaging as a method to improve dosimetry. Overall, our software pipeline helps estimate acoustic exposure, and our study highlights the need for image feedback to increase the accuracy of tFUS dosimetry.

Transcranial focused ultrasound-mediated neurochemical and functional connectivity changes in deep cortical regions in humans

Low-intensity transcranial ultrasound stimulation (TUS) is an emerging non-invasive technique for focally modulating human brain function. The mechanisms and neurochemical substrates underlying TUS neuromodulation in humans and how these relate to excitation and inhibition are still poorly understood. In 24 healthy controls, we separately stimulated two deep cortical regions and investigated the effects of theta-burst TUS, a protocol shown to increase corticospinal excitability, on the inhibitory neurotransmitter gamma-aminobutyric acid (GABA) and functional connectivity. We show that theta-burst TUS in humans selectively reduces GABA levels in the posterior cingulate, but not the dorsal anterior cingulate cortex. Functional connectivity increased following TUS in both regions. Our findings suggest that TUS changes overall excitability by reducing GABAergic inhibition and that changes in TUS-mediated neuroplasticity last at least 50 mins after stimulation. The difference in TUS effects on the posterior and anterior cingulate could suggest state- or location-dependency of the TUS effect-both mechanisms increasingly recognized to influence the brain's response to neuromodulation.

A head template for computational dose modelling for transcranial focused ultrasound stimulation

Transcranial focused Ultrasound Stimulation (TUS) at low intensities is emerging as a novel non-invasive brain stimulation method with higher spatial resolution than established transcranial stimulation methods and the ability to selectively stimulate also deep brain areas. Accurate control of the focus position and strength of the TUS acoustic waves is important to enable a beneficial use of the high spatial resolution and to ensure safety. As the human skull causes strong attenuation and distortion of the waves, simulations of the transmitted waves are needed to accurately determine the TUS dose distribution inside the cranial cavity. The simulations require information of the skull morphology and its acoustic properties. Ideally, they are informed by computed tomography (CT) images of the individual head. However, suited individual imaging data is often not readily available. For this reason, we here introduce and validate a head template that can be used to estimate the average effects of the skull on the TUS acoustic wave in the population. The template was created from CT images of the heads of 29 individuals of different ages (between 20-50 years), gender and ethnicity using an iterative non-linear co-registration procedure. For validation, we compared acoustic and thermal simulations based on the template to the average of the simulation results of all 29 individual datasets. Acoustic simulations were performed for a model of a focused transducer driven at 500 kHz, placed at 24 standardized positions by means of the EEG 10-10 system. Additional simulations at 250 kHz and 750 kHz at 16 of the positions were used for further confirmation. The amount of ultrasound-induced heating at 500 kHz was estimated for the same 16 transducer positions. Our results show that the template represents the median of the acoustic pressure and temperature maps from the individuals reasonably well in most cases. This underpins the usefulness of the template for the planning and optimization of TUS interventions in studies of healthy young adults. Our results further indicate that the amount of variability between the individual simulation results depends on the position. Specifically, the simulated ultrasound-induced heating inside the skull exhibited strong interindividual variability for three posterior positions close to the midline, caused by a high variability of the local skull shape and composition. This should be taken into account when interpreting simulation results based on the template.

Mechanisms of theta burst transcranial ultrasound induced plasticity in the human motor cortex

BACKGROUND

Transcranial ultrasound stimulation (TUS) is a novel non-invasive brain stimulation technique with high depth penetrance and spatial resolution. Theta-burst TUS (tbTUS) is a plasticity-inducing protocol which increases motor cortical excitability for up to 30 min following 80s of sonication. While this protocol may have therapeutic potential for the treatment of psychiatric and neurological disorders, the mechanisms of action of TUS remain unclear.

OBJECTIVE

We conducted the first pharmacological study to examine the mechanisms of TUS in human primary motor cortex. By administering brain-active drugs with known mechanisms of action, we aimed to elucidate the mechanisms of tbTUS.

METHODS

Fourteen healthy subjects participated in a within-subjects randomized, double-blind, cross-over study with five visits. At each visit, one of four study drugs (carbamazepine - Na+ channel blocker, nimodipine - L-type Ca2+ channel blocker, lorazepam - positive allosteric modulator of gamma-aminobutyric acid (GABA) type A receptor, dextromethorphan - N-methyl-d-aspartate receptor antagonist) or placebo was administered in random order, followed by tbTUS.

RESULTS

The plasticity effects of tbTUS on motor cortex excitability measured by motor-evoked potential amplitudes elicited by transcranial magnetic stimulation were reduced by all study drugs compared to placebo.

CONCLUSION

tbTUS may induce NMDA-dependent synaptic plasticity since the effects are blocked by increased GABAA receptor activities and voltage-gated Na+ and Ca2+ channels blockers. These results are consistent with the hypotheses that tbTUS induced long-term potentiation-like mechanisms and that TUS involves activation of mechanosensitive Na+ and Ca2+ channels. Alternatively, non-specific pharmacologically induced changes in excitatory/inhibitory balance might have interfered with the effects of tbTUS.

Global sonication of the human intracranial space via a jumbo planar transducer

Contrary to conditioning a Focused Ultrasound (FUS) beam to sonicate a localized region of the human brain, the goal of this investigation was to explore the prospect of distributing homogeneous ultrasound energy over the entire brain space with a large cranium-wide ultrasound beam. Recent ultrasound preclincal studies utilizing large or whole brain stimulation regions create a demand for expanding the treatment envelope of transcranial pulsed-low intensity ultrasound towards Global Brain Sonication (GBS) for potential human investigation. Here, we conduct ultrasound field characterizations when transmitting pulsed ultrasound through human skull specimens using a 1-3 piezocomposite planar transducer operating at 464 kHz with an active single-element surface of 30 × 30 cm. Through computational simulation and hydrophone scanning methodology, ultrasound wave behavior and dose homogeneity in the brain space were evaluated under various trajectories of sonication using the planar transducer. Clinically relevant pulse parameters used for transcranial therapeutic ultrasound applications were used in the experiments. Simulations and empirical testing revealed that dose homogeneity and acoustic intensity over the brain space are influenced by sonication trajectory, skull lens effects, and acoustic wave reflections. The transducer can emit a spatial peak pulse average intensity of 4.03 W/cm² (0.24 MPa) measured in the free-field at 464 kHz with electrical power of 1 kW. The simulation showed that approximately 99 % of the cranial volume was exposed with <30 % of the maximum external acoustic intensity being transmitted into the skull. The transmission loss across all sonication trajectories is similar to previously reported FUS studies. A marker for GBS dose homogeneity is introduced to score the ultrasound pressure field uniformity in the intracranial space. Results of this study identify the initial challenges of exposing the entire human brain space with ultrasound using a large cranium-wide sonication beam intended for global brain therapeutic modulation.

Acoustic-responsive carbon dioxide-loaded liposomes for efficient drug release

The role of liposomes as drug carriers has been investigated. Ultrasound-based drug release methods have been developed for on-demand drug delivery. However, the acoustic responses of current liposome carriers result in low drug release efficiency. In this study, CO₂-loaded liposomes were synthesized under high pressure from supercritical CO₂ and irradiated with ultrasound at 237 kHz to demonstrate their superior acoustic responsiveness. When liposomes containing fluorescent drug models were irradiated with ultrasound under acoustic pressure conditions that are safe for the human body, CO₂-loaded liposomes synthesized using supercritical CO₂ had 17.1 times higher release efficiency than liposomes synthesized using the conventional Bangham method. In particular, the release efficiency of CO₂-loaded liposomes synthesized using supercritical CO₂ and monoethanolamine was 19.8 times higher than liposomes synthesized using the conventional Bangham method. These findings on the release efficiency of acoustic-responsive liposomes suggest an alternative liposome synthesis strategy for on-demand release of drugs by ultrasound irradiation in future therapies.

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.

Weak Ultrasound Contributes to Neuromodulatory Effects in the Rat Motor Cortex

Transcranial focused ultrasound (tFUS) is a novel neuromodulating technique. It has been demonstrated that the neuromodulatory effects can be induced by weak ultrasound exposure levels (spatial-peak temporal average intensity, I_SPTA < 10 mW/cm²) in vitro. However, fewer studies have examined the use of weak tFUS to potentially induce long-lasting neuromodulatory responses in vivo. The purpose of this study was to determine the lower-bound threshold of tFUS stimulation for inducing neuromodulation in the motor cortex of rats. A total of 94 Sprague-Dawley rats were used. The sonication region aimed at the motor cortex under weak tFUS exposure (I_SPTA of 0.338-12.15 mW/cm²). The neuromodulatory effects of tFUS on the motor cortex were evaluated by the changes in motor-evoked potentials (MEPs) elicited by transcranial magnetic stimulation (TMS). In addition to histology analysis, the in vitro cell culture was used to confirm the neuromodulatory mechanisms following tFUS stimulation. In the results, the dose-dependent inhibitory effects of tFUS were found, showing increased intensities of tFUS suppressed MEPs and lasted for 30 min. Weak tFUS significantly decreased the expression of excitatory neurons and increased the expression of inhibitory GABAergic neurons. The PIEZO-1 proteins of GABAergic neurons were found to involve in the inhibitory neuromodulation. In conclusion, we show the use of weak ultrasound to induce long-lasting neuromodulatory effects and explore the potential use of weak ultrasound for future clinical neuromodulatory applications.

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.

Evaluation of the effects of focused ultrasound stimulation on the central nervous system through a multiscale simulation approach

The lack of sensory feedback represents one of the main drawbacks of commercial upper limb prosthesis. Transcranial Focused Ultrasound Stimulation (tFUS) seems to be a valid non-invasive technique for restoring sensory feedback allowing to deliver acoustic energy to cortical sensory areas with high spatial resolution and depth penetration. This paper aims at studying in simulation the use of tFUS on cortical sensory areas to evaluate its effects in terms of latency ad firing rate of the cells response, for understanding if these parameters influence the safety and the efficacy of the stimulation. In this paper, in order to study the propagation of the ultrasound wave from the transducer to the cortical cells, a multiscale approach was implemented by building a macroscopic model, which estimates the pressure profile in a simplified 2D human head geometry, and coupling it with the SONIC microscale model, that describes the electrical behaviour of a cortical neuron. The influence of the stimulation parameters and of the skull thickness on the latency and the firing rate are evaluated and the obtained behaviour is linked to the sensory response obtained on human subjects. Results have shown that slight changes in the transducer position should not affect the efficacy of the stimulation; however, high skull thickness leads to lower cells activation. These results will be useful for evaluating safety and effectiveness of tFUS for sensory feedback in closed-loop prosthetic systems.

Application of subject-specific helmets for the study of human visuomotor behavior using transcranial focused ultrasound: a pilot study

BACKGROUND AND OBJECTIVE

As a novel non-invasive human brain stimulation method, transcranial focused ultrasound (tFUS) is receiving growing attention due to its superior spatial specificity and depth penetrability. Since the focal point of tFUS needs to be fixated precisely to the target brain region during stimulation, a critical issue is to identify and maintain the accurate position and orientation of the tFUS transducer relative to the subject's head. This study aims to propose the entire framework of tFUS stimulation integrating the methods previously proposed by the authors for tFUS transducer configuration optimization and a subject-specific 3D-printed helmet, and to validate this complete setup in a human behavioral neuromodulation study.

METHODS

To find the optimal configuration of the tFUS transducer, a numerical method based on subject-specific tFUS beamlines simulation was used. Then, the subject-specific 3D-printed helmet has been applied to effectively secure the transducer at the estimated optimal configuration. To validate this tFUS framework, a common behavioral neuromodulation paradigm was chosen; the effect of the dorsolateral prefrontal cortex (DLPFC) stimulation on anti-saccade (AS) behavior. While human participants (n=2) were performing AS tasks, tFUS stimulations were randomly applied to the left DLPFC right after the fixation target disappeared.

RESULTS

The neuromodulation result strongly suggests that the cortical stimulation using the proposed tFUS setup is effective in significantly reducing the error rates of anti-saccades (about -10 %p for S1 and -16 %p for S2), whereas no significant effect was observed on their latencies. These observed behavioral effects are consistent with the previous results based on conventional brain stimulation or lesion studies.

CONCLUSIONS

The proposed subject-specific tFUS framework has been effectively used in human neuromodulation study. The result suggests that the tFUS stimulation targeted to the DLPFC can generate a neuromodulatory effect on AS behavior.

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.

Transcranial ultrasound stimulation to human middle temporal complex improves visual motion detection and modulates electrophysiological responses

BACKGROUND

Transcranial ultrasound stimulation (TUS) holds promise as a novel technology for non-invasive neuromodulation, with greater spatial precision than other available methods and the ability to target deep brain structures. However, its safety and efficacy for behavioural and electrophysiological modulation remains controversial and it is not yet clear whether it can be used to manipulate the neural mechanisms supporting higher cognitive function in humans. Moreover, concerns have been raised about a potential TUS-induced auditory confound.

OBJECTIVES

We aimed to investigate whether TUS can be used to modulate higher-order visual function in humans in an anatomically-specific way whilst controlling for auditory confounds.

METHODS

We used participant-specific skull maps, functional localisation of brain targets, acoustic modelling and neuronavigation to guide TUS delivery to human visual motion processing cortex (hMT+) whilst participants performed a visual motion detection task. We compared the effects of hMT+ stimulation with sham and control site stimulation and examined EEG data for modulation of task-specific event-related potentials. An auditory mask was applied which prevented participants from distinguishing between stimulation and sham trials.

RESULTS

Compared with sham and control site stimulation, TUS to hMT+ improved accuracy and reduced response times of visual motion detection. TUS also led to modulation of the task-specific event-related EEG potential. The amplitude of this modulation correlated with the performance benefit induced by TUS. No pathological changes were observed comparing structural MRI obtained before and after stimulation.

CONCLUSIONS

The results demonstrate for the first time the precision, efficacy and safety of TUS for stimulation of higher-order cortex and cognitive function in humans whilst controlling for auditory confounds.

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

OBJECTIVE

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

METHODS

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

RESULTS

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

SIGNIFICANCE

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

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.

Modelling transcranial ultrasound neuromodulation: an energy-based multiscale framework

Several animal and human studies have now established the potential of low intensity, low frequency transcranial ultrasound (TUS) for non-invasive neuromodulation. Paradoxically, the underlying mechanisms through which TUS neuromodulation operates are still unclear, and a consensus on the identification of optimal sonication parameters still remains elusive. One emerging hypothesis based on thermodynamical considerations attributes the acoustic-induced nerve activity alterations to the mechanical energy and/or entropy conversions occurring during TUS action. Here, we propose a multiscale modelling framework to examine the energy states of neuromodulation under TUS. First, macroscopic tissue-level acoustic simulations of the sonication of a whole monkey brain are conducted under different sonication protocols. For each one of them, mechanical loading conditions of the received waves in the anterior cingulate cortex region are recorded and exported into a microscopic cell-level 3D viscoelastic finite element model of neuronal axon embedded extracellular medium. Pulse-averaged elastically stored and viscously dissipated energy rate densities during axon deformation are finally computed under different sonication incident angles and are mapped against distinct combinations of sonication parameters of the TUS. The proposed multiscale framework allows for the analysis of vibrational patterns of the axons and its comparison against the spectrograms of stimulating ultrasound. The results are in agreement with literature data on neuromodulation, demonstrating the potential of this framework to identify optimised acoustic parameters in TUS neuromodulation. The proposed approach is finally discussed in the context of multiphysics energetic considerations, argued here to be a promising avenue towards a scalable framework for TUS in silico predictions. STATEMENT OF SIGNIFICANCE: Low-intensity transcranial ultrasound (TUS) is poised to become a leading neuromodulation technique for the treatment of neurological disorders. Paradoxically, how it operates at the cellular scale remains unknown, hampering progress in personalised treatment. To this end, models of the multiphysics of neurons able to upscale results to the organ scale are required. We propose here to achieve this by considering an axon submitted to an ultrasound wave extracted from a simulation at the organ scale. Doing so, information pertaining to both stored and dissipated axonal energies can be extracted for a given head/brain morphology. This two-scale multiphysics energetic approach is a promising scalable framework for in silico predictions in the context of personalised TUS treatment.

A Comparative Feasibility Study for Transcranial Extracorporeal Shock Wave Therapy

The potential beneficial regenerative and stimulatory extracorporeal shock wave therapy (ESWT) applications to the central nervous system have garnered interest in recent years. Treatment zones for these indications are acoustically shielded by bones, which heavily impact generated sound fields. We present the results of high-resolution tissue-realistic simulations, comparing the viability of different ESWT applicators in their use for transcranial applications. The performances of electrohydraulic, electromagnetic, and piezoelectric transducers for key reflector geometries are compared. Based on density information obtained from CT imaging of the head, we utilized the non-linear wave propagation toolset Matlab k-Wave to obtain spatial therapeutic sound field geometries and waveforms. In order to understand the reliability of results on the appropriate modeling of the skull, three different bone attenuation models were compared. We find that all currently clinically ESWT applicator technologies show significant retention of peak pressures and energies past the bone barrier. Electromagnetic transducers maintain a significantly higher energy flux density compared to other technologies while low focusing strength piezoelectric applicators have the weakest transmissions. Attenuation estimates provide insights into sound field degradation and energy losses, indicating that effective transcranial therapies can readily be attained with current applicators. Furthermore, the presented approach will allow for future targeted in silico development and the design of applicators and therapy plans to ultimately improve therapeutic outcomes.

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

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

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.

Mechanistic insights into ultrasonic neurostimulation of disconnected neurons using single short pulses

Ultrasonic neurostimulation is a potentially potent noninvasive therapy, whose mechanism has yet to be elucidated. We designed a system capable of applying ultrasound with minimal reflections to neuronal cultures. Synaptic transmission was pharmacologically controlled, eliminating network effects, enabling examination of single-cell processes. Short single pulses of low-intensity ultrasound were applied, and time-locked responses were examined using calcium imaging. Low-pressure (0.35MPa) ultrasound directly stimulated ∼20% of pharmacologically disconnected neurons, regardless of membrane poration. Stimulation was resistant to the blockade of several purinergic receptor and mechanosensitive ion channel types. Stimulation was blocked, however, by suppression of action potentials. Surprisingly, even extremely short (4μs) pulses were effective, stimulating ∼8% of the neurons. Lower-pressure pulses (0.35MPa) were less effective than higher-pressure ones (0.65MPa). Attrition effects dominated, with no indication of compromised viability. Our results detract from theories implicating cavitation, heating, non-transient membrane pores >1.5nm, pre-synaptic release, or gradual effects. They implicate a post-synaptic mechanism upstream of the action potential, and narrow down the list of possible targets involved.

Current State of Potential Mechanisms Supporting Low Intensity Focused Ultrasound for Neuromodulation

Low intensity focused ultrasound (LIFU) has been gaining traction as a non-invasive neuromodulation technology due to its superior spatial specificity relative to transcranial electrical/magnetic stimulation. Despite a growing literature of LIFU-induced behavioral modifications, the mechanisms of action supporting LIFU's parameter-dependent excitatory and suppressive effects are not fully understood. This review provides a comprehensive introduction to the underlying mechanics of both acoustic energy and neuronal membranes, defining the primary variables for a subsequent review of the field's proposed mechanisms supporting LIFU's neuromodulatory effects. An exhaustive review of the empirical literature was also conducted and studies were grouped based on the sonication parameters used and behavioral effects observed, with the goal of linking empirical findings to the proposed theoretical mechanisms and evaluating which model best fits the existing data. A neuronal intramembrane cavitation excitation model, which accounts for differential effects as a function of cell-type, emerged as a possible explanation for the range of excitatory effects found in the literature. The suppressive and other findings need additional theoretical mechanisms and these theoretical mechanisms need to have established relationships to sonication parameters.

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.

Mathematical Model of Ultrasound Attenuation With Skull Thickness for Transcranial-Focused Ultrasound

Transcranial-focused ultrasound (tFUS) has potential for both neuromodulation and neuroimaging. Due to the influence of head tissue, especially the skull, its attenuation is a key issue affecting precise focusing. The objective of the present study was to construct a mathematical model of ultrasound attenuation inclusive of skull thickness. First, combined with real skull phantom experiments and simulation experiments, tFUS attenuation of different head tissues was investigated. Furthermore, based on the system identification method, a mathematical model of ultrasound attenuation was constructed taking skull thickness into account. Finally, the performance of the mathematical model was tested, and its potential applications were investigated. For different head tissues, including scalp, skull, and brain tissue, the skull was found to be the biggest influencing factor for ultrasound attenuation, the attenuation caused by it being 4.70 times and 7.06 times that of attenuation caused by the brain and scalp, respectively. Consistent with the results of both the simulation and phantom experiments, the attenuation of the mathematical model increased as the skull thickness increased. The average error of the mathematical model was 1.87% in the phantom experiment. In addition, the experimental results show that the devised mathematical model is suitable for different initial pressures and different skulls with correlation coefficients higher than 0.99. Both simulation and phantom experiments validated the effectiveness of the proposed mathematical model. It can be concluded from this experiment that the proposed mathematical model can accurately calculate the tFUS attenuation and can significantly contribute to further research and application of tFUS.

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.

First evidence of long-term effects of transcranial pulse stimulation (TPS) on the human brain

Background

With the high spatial resolution and the potential to reach deep brain structures, ultrasound-based brain stimulation techniques offer new opportunities to non-invasively treat neurological and psychiatric disorders. However, little is known about long-term effects of ultrasound-based brain stimulation. Applying a longitudinal design, we comprehensively investigated neuromodulation induced by ultrasound brain stimulation to provide first sham-controlled evidence of long-term effects on the human brain and behavior.

Methods

Twelve healthy participants received three sham and three verum sessions with transcranial pulse stimulation (TPS) focused on the cortical somatosensory representation of the right hand. One week before and after the sham and verum TPS applications, comprehensive structural and functional resting state MRI investigations and behavioral tests targeting tactile spatial discrimination and sensorimotor dexterity were performed.

Results

Compared to sham, global efficiency significantly increased within the cortical sensorimotor network after verum TPS, indicating an upregulation of the stimulated functional brain network. Axial diffusivity in left sensorimotor areas decreased after verum TPS, demonstrating an improved axonal status in the stimulated area.

Conclusions

TPS increased the functional and structural coupling within the stimulated left primary somatosensory cortex and adjacent sensorimotor areas up to one week after the last stimulation. These findings suggest that TPS induces neuroplastic changes that go beyond the spatial and temporal stimulation settings encouraging further clinical applications.

Low Intensity Focused Ultrasound for Non-invasive and Reversible Deep Brain Neuromodulation-A Paradigm Shift in Psychiatric Research

This article describes an emerging non-invasive neuromodulatory technology, called low intensity focused ultrasound (LIFU). This technology is potentially paradigm shifting as it can deliver non-invasive and reversible deep brain neuromodulation through acoustic sonication, at millimeter precision. Low intensity focused ultrasound's spatial precision, yet non-invasive nature sets it apart from current technologies, such as transcranial magnetic or electrical stimulation and deep brain stimulation. Additionally, its reversible effects allow for the causal study of deep brain regions implicated in psychiatric illness. Studies to date have demonstrated that LIFU can safely modulate human brain activity at cortical and subcortical levels. Due to its novelty, most researchers and clinicians are not aware of the potential applications and promise of this technique, underscoring the need for foundational papers to introduce the community to LIFU. This mini-review and synthesis of recent advances examines several key papers on LIFU administered to humans, describes the population under study, parameters used, and relevant findings that may guide future research. We conclude with a concise overview of some of the more pressing questions to date, considerations when interpreting new data from an emerging field, and highlight the opportunities and challenges in this exciting new area of study.

Induction of Human Motor Cortex Plasticity by Theta Burst Transcranial Ultrasound Stimulation

OBJECTIVE

Transcranial ultrasound stimulation (TUS) is a promising non-invasive brain stimulation technique with advantages of high spatial precision and ability to target deep brain regions. This study aimed to develop a TUS protocol to effectively induce brain plasticity in human subjects.

METHODS

An 80 s train of theta burst patterned TUS (tbTUS), regularly patterned TUS (rTUS) with the same sonication duration and sham tbTUS were delievered to the motor cortex in healthy subjects. Transcranial magnetic stimulation (TMS) was used to examine changes in corticospinal excitability, intracortical inhibition and facilitation, and the site of plasticity induction. The effects of motor cortical tbTUS on a visuo-motor task and the effects of occipital cortex tbTUS on motor cortical excitability were also tested.

RESULTS

The tbTUS produced consistent increase in corticospinal excitability for at least 30 minutes while rTUS and sham tbTUS produced no significant change. tbTUS decreased short-interval intracortical inhibiton and increased intracortical facilitation. The effects of TMS in different current directions suggested that the site of the plastic changes was within the motor cortex. tbTUS to the occipital cortex did not change motor cortical excitability. Motor cortical tbTUS shortened movement time in a visuo-motor task.

INTERPRETATION

tbTUS is a novel and efficient paradigm to induce cortical plasticity in human. It has the potential to be developed for neuromodulation treatment for neurological and psychiatric disorders, and to advance neuroscience research. This article is protected by copyright. All rights reserved.

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.

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.

Acoustic simulation for transcranial focused ultrasound using GAN-based synthetic CT

Transcranial focused ultrasound (tFUS) is a promising non-invasive technique for treating neurological and psychiatric disorders. One of the challenges for tFUS is the disruption of wave propagation through the skull. Consequently, despite the risks associated with exposure to ionizing radiation, computed tomography (CT) is required to estimate the acoustic transmission through the skull. This study aims to generate synthetic CT (sCT) from T1-weighted magnetic resonance imaging (MRI) and investigate its applicability to tFUS acoustic simulation. We trained a 3D conditional generative adversarial network (3D-cGAN) with 15 subjects. We then assessed image quality with 15 test subjects: mean absolute error (MAE) = 85.72± 9.50 HU (head) and 280.25±24.02 HU (skull), dice coefficient similarity (DSC) = 0.88±0.02 (skull). In terms of skull density ratio (SDR) and skull thickness (ST), no significant difference was found between sCT and real CT (rCT). When the acoustic simulation results of rCT and sCT were compared, the intracranial peak acoustic pressure ratio was found to be less than 4%, and the distance between focal points less than 1 mm.