Ex vivo optimisation of a heterogeneous speed of sound model of the human skull for non-invasive transcranial focused ultrasound at 1 MHz

Sustained reduction of essential tremor with low-power non-thermal transcranial focused ultrasound stimulations in humans

BACKGROUND

Transcranial ultrasound stimulation (TUS) is a non-invasive brain stimulation technique; when skull aberrations are compensated for, this technique allows, with millimetric accuracy, circumvention of the invasive surgical procedure associated with deep brain stimulation (DBS) and the limited spatial specificity of transcranial magnetic stimulation.

OBJECTIVE

/hypothesis: We hypothesize that MR-guided low-power TUS can induce a sustained decrease of tremor power in patients suffering from medically refractive essential tremor.

METHODS

The dominant hand only was targeted, and two anatomical sites were sonicated in this exploratory study: the ventral intermediate nucleus of the thalamus (VIM) and the dentato-rubro-thalamic tract (DRT). Patients (N = 9) were equipped with MR-compatible accelerometers attached to their hands to monitor their tremor in real-time during TUS.

RESULTS

VIM neurostimulations followed by a low-duty cycle (5 %) DRT stimulation induced a substantial decrease in the tremor power in four patients, with a minimum of 89.9 % reduction when compared with the baseline power a few minutes after the DRT stimulation. The only patient stimulated in the VIM only and with a low duty cycle (5 %) also experienced a sustained reduction of the tremor (up to 93.4 %). Four patients (N = 4) did not respond. The temperature at target was 37.2 ± 1.4 °C compared to 36.8 ± 1.4 °C for a 3 cm away control point.

CONCLUSIONS

MR-guided low power TUS can induce a substantial and sustained decrease of tremor power. Follow-up studies need to be conducted to reproduce the effect and better to understand the variability of the response amongst patients. MR thermometry during neurostimulations showed no significant thermal rise, supporting a mechanical effect.

Transcranial Blood-Brain Barrier Opening in Alzheimer's Disease Patients Using a Portable Focused Ultrasound System with Real-Time 2-D Cavitation Mapping

Background

Focused ultrasound (FUS) in combination with microbubbles has recently shown great promise in facilitating blood-brain barrier (BBB) opening for drug delivery and immunotherapy in Alzheimer's disease (AD). However, it is currently limited to systems integrated within the MRI suites or requiring post-surgical implants, thus restricting its widespread clinical adoption. In this pilot study, we investigate the clinical safety and feasibility of a portable, non-invasive neuronavigation-guided FUS (NgFUS) system with integrated real-time 2-D microbubble cavitation mapping.

Methods

A phase 1 clinical study with mild to moderate AD patients (N=6) underwent a single session of microbubble-mediated NgFUS to induce transient BBB opening (BBBO). Microbubble activity under FUS was monitored with real-time 2-D cavitation maps and dosing to ensure the efficacy and safety of the NgFUS treatment. Post-operative MRI was used for BBB opening and closure confirmation as well as safety assessment. Changes in AD biomarker levels in both blood serum and extracellular vesicles (EVs) were evaluated, while changes in amyloid-beta (Aβ) load in the brain were assessed through 18F-Florbetapir PET.

Results

BBBO was achieved in 5 out of 6 subjects with an average volume of 983±626 mm3 following FUS at the right frontal lobe both in white and gray matter regions. The outpatient treatment was completed within 34.8±10.7 min. Cavitation dose significantly correlated with the BBBO volume (R2>0.9, N=4), demonstrating the portable NgFUS system's capability of predicting opening volumes. The cavitation maps co-localized closely with the BBBO location, representing the first report of real-time transcranial 2-D cavitation mapping in the human brain. Larger opening volumes correlated with increased levels of AD biomarkers, including Aβ42 (R2=0.74), Tau (R2=0.95), and P-Tau181 (R2=0.86), assayed in serum-derived EVs sampled 3 days after FUS (N=5). From PET scans, subjects showed a lower Aβ load increase in the treated frontal lobe region compared to the contralateral region. Reduction in asymmetry standardized uptake value ratios (SUVR) correlated with the cavitation dose (R2>0.9, N=3). Clinical changes in the mini-mental state examination over 6 months were within the expected range of cognitive decline with no additional changes observed as a result of FUS.

Conclusion

We showed the safety and feasibility of this cost-effective and time-efficient portable NgFUS treatment for BBBO in AD patients with the first demonstration of real-time 2-D cavitation mapping. The cavitation dose correlated with BBBO volume, a slowed increase in pathology, and serum detection of AD proteins. Our study highlights the potential for accessible FUS treatment in AD, with or without drug delivery.

New semi-analytical method for fast transcranial ultrasonic field simulation

Objective.To optimize and ensure the safety of ultrasound brain therapy, personalized transcranial ultrasound simulations are very useful. They allow to predict the pressure field, depending on the patient skull and probe position. Most transcranial ultrasound simulations are based on numerical methods which have a long computation time and a high memory usage. The goal of this study is to develop a new semi-analytical field computation method that combines realism and computation speed.Approach.Instead of the classic ray tracing, the ultrasonic paths are computed by time of flight minimization. Then the pressure field is computed using the pencil method. This method requires a smooth and homogeneous skull model. The simulation algorithm, so-called SplineBeam, was numerically validated, by comparison with existing solvers, and experimentally validated by comparison with hydrophone measured pressure fields through anex vivohuman skull.Main results.SplineBeam simulated pressure fields were close to the experimentally measured ones, with a focus position difference of the order of the positioning error and a maximum pressure difference lower than 6.02%. In addition, for those configurations, SplineBeam computation time was lower than another simulation software, k-Wave's, by two orders of magnitude, thanks to its capacity to compute the field only at the focal spot.Significance.These results show the potential of this new method to compute fast and realistic transcranial pressure fields. The combination of this two assets makes it a promising tool for real time transcranial pressure field prediction during ultrasound brain therapy interventions.

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

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

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.

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.

Enhanced Ultrasound Transmission through Skull Using Flexible Matching Layer with Gradual Acoustic Impedance

Transcranial ultrasound imaging and therapy have gained significant attention due to their noninvasive nature, absence of ionizing radiation, and portability. However, the presence of the skull, which has a high acoustic impedance, presents a challenge for the penetration of ultrasound into intracranial tissue. This leads to a low transmission of ultrasound through the skull, hindering energy focusing and imaging quality. To address this challenge, we propose a novel approach that utilizes a flexible matching layer with gradual acoustic impedance to enhance ultrasound transmission through the skull. This matching layer is constructed using Poly(dimethylsiloxane) (PDMS)/tungsten powders as the structural component responsible for the gradual impedance, while agarose serves as the flexible matrix. Our simulation and experimental results demonstrate that the matching layer with an exponential gradual acoustic impedance significantly improves the ultrasound transmission coefficient across a wide frequency range compared to traditional quarter wavelength matching layers. Specifically, at 2 MHz, the maximum transmission coefficient reaches 49.5%, more than four times higher than that of the skull without a matching layer (only 11.7%). Additionally, the good flexibility of our matching layer ensures excellent adhesion to the curved surface of the skull, further enhancing its application potential in transcranial ultrasound imaging and therapy. The improved transmission performance allows for a lower ultrasound transmission power, effectively addressing overheating and safety issues.

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.

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.

Characterization of the Targeting Accuracy of a Neuronavigation-Guided Transcranial FUS System In Vitro, In Vivo, and In Silico

Focused ultrasound (FUS)-enabled liquid biopsy (sonobiopsy) is an emerging technique for the noninvasive and spatiotemporally controlled diagnosis of brain cancer by inducing blood-brain barrier (BBB) disruption to release brain tumor-specific biomarkers into the blood circulation. The feasibility, safety, and efficacy of sonobiopsy were demonstrated in both small and large animal models using magnetic resonance-guided FUS devices. However, the high cost and complex operation of magnetic resonance-guided FUS devices limit the future broad application of sonobiopsy in the clinic. In this study, a neuronavigation-guided sonobiopsy device is developed and its targeting accuracy is characterized in vitro, in vivo, and in silico. The sonobiopsy device integrated a commercially available neuronavigation system (BrainSight) with a nimble, lightweight FUS transducer. Its targeting accuracy was characterized in vitro in a water tank using a hydrophone. The performance of the device in BBB disruption was verified in vivo using a pig model, and the targeting accuracy was quantified by measuring the offset between the target and the actual locations of BBB opening. The feasibility of the FUS device in targeting glioblastoma (GBM) tumors was evaluated in silico using numerical simulation by the k-Wave toolbox in glioblastoma patients. It was found that the targeting accuracy of the neuronavigation-guided sonobiopsy device was 1.7 ± 0.8 mm as measured in the water tank. The neuronavigation-guided FUS device successfully induced BBB disruption in pigs with a targeting accuracy of 3.3 ± 1.4 mm. The targeting accuracy of the FUS transducer at the GBM tumor was 5.5 ± 4.9 mm. Age, sex, and incident locations were found to be not correlated with the targeting accuracy in GBM patients. This study demonstrated that the developed neuronavigation-guided FUS device could target the brain with a high spatial targeting accuracy, paving the foundation for its application in the clinic.

Transcranial 3D ultrasound localization microscopy using a large element matrix array with a multi-lens diffracting layer: anin vitrostudy

Objective. Early diagnosis and acute knowledge of cerebral disease require to map the microflows of the whole brain. Recently, ultrasound localization microscopy (ULM) was applied to map and quantify blood microflows in 2D in the brain of adult patients down to the micron scale. Whole brain 3D clinical ULM remains challenging due to the transcranial energy loss which reduces significantly the imaging sensitivity.Approach. Large aperture probes with a large surface can increase both the field of view and sensitivity. However, a large active surface implies thousands of acoustic elements, which limits clinical translation. In a previous simulation study, we developed a new probe concept combining a limited number of elements and a large aperture. It is based on large elements, to increase sensitivity, and a multi-lens diffracting layer to improve the focusing quality. In this study, a 16 elements prototype, driven at 1 MHz frequency, was made andin vitroexperiments were performed to validate the imaging capabilities of this new probe concept.Main results. First, pressure fields emitted from a large single transducer element without and with diverging lens were compared. Low directivity was measured for the large element with the diverging lens while maintaining high transmit pressure. The focusing quality of 4 × 3cm matrix arrays of 16 elements without/with lenses were compared.In vitroexperiments in a water tank and through a human skull were achieved to localize and track microbubbles in tubes.Significance.ULM was achieved demonstrating the strong potential of multi-lens diffracting layer to enable microcirculation assessment over a large field of view through the bones.

Coregistered transcranial optoacoustic and magnetic resonance angiography of the human brain

Imaging modalities capable of visualizing the human brain have led to major advances in neurology and brain research. Multi-spectral optoacoustic tomography (MSOT) has gained importance for studying cerebral function in rodent models due to its unique capability to map changes in multiple hemodynamic parameters and to directly visualize neural activity within the brain. The technique further provides molecular imaging capabilities that can facilitate early disease diagnosis and treatment monitoring. However, transcranial imaging of the human brain is hampered by acoustic attenuation and other distortions introduced by the skull. Here, we demonstrate non-invasive transcranial MSOT angiography of pial veins through the temporal bone of an adult healthy volunteer. Time-of-flight (TOF) magnetic resonance angiography (MRA) and T₁-weighted structural magnetic resonance imaging (MRI) were further acquired to facilitate anatomical registration and interpretation. The superior middle cerebral vein in the temporal cortex was identified in the MSOT images, matching its location observed in the TOF-MRA images. These initial results pave the way toward the application of MSOT in clinical brain imaging.

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.

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.

Binary acoustic metasurfaces for dynamic focusing of transcranial ultrasound

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

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

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

Two-step aberration correction: application to transcranial histotripsy

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

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.

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.

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.

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

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

Transcranial focused ultrasound phase correction using the hybrid angular spectrum method

The InSightec Exablate system is the standard of care used for transcranial focused ultrasound ablation treatments in the United States. The system calculates phase corrections that account for aberrations caused by the human skull. This work investigates whether skull aberration correction can be improved by comparing the standard of care InSightec ray tracing method with the hybrid angular spectrum (HAS) method and the gold standard hydrophone method. Three degassed ex vivo human skulls were sonicated with a 670 kHz hemispherical phased array transducer (InSightec Exablate 4000). Phase corrections were calculated using four different methods (straight ray tracing, InSightec ray tracing, HAS, and hydrophone) and were used to drive the transducer. 3D raster scans of the beam profiles were acquired using a hydrophone mounted on a 3-axis positioner system. Focal spots were evaluated using six metrics: pressure at the target, peak pressure, intensity at the target, peak intensity, positioning error, and focal spot volume. For three skulls, the InSightec ray tracing method achieved 52 ± 21% normalized target intensity (normalized to hydrophone), 76 ± 17% normalized peak intensity, and 0.72 ± 0.47 mm positioning error. The HAS method achieved 74 ± 9% normalized target intensity, 81 ± 9% normalized peak intensity, and 0.35 ± 0.09 mm positioning error. The InSightec-to-HAS improvement in focal spot targeting provides promise in improving treatment outcomes. These improvements to skull aberration correction are also highly relevant for the applications of focused ultrasound neuromodulation and blood brain barrier opening, which are currently being translated for human use.

An efficient and accurate parallel hybrid acoustic signal correction method for transcranial ultrasound

Phased-control focused ultrasound transducers provide a new and noninvasive treatment method for brain disease. However, improving the accuracy of phase correction and reducing the calculation time during treatment have always been contradictory constraints. In this paper, a hybrid acoustic signal correction (HASC) method combined with k-Wave stage and holography stage was introduced for phase correction and simulation of transcranial focused ultrasound. The k-Wave stage is mainly used to calculate the sound field in a heterogeneous medium (skull), which divides the sound field calculation process into paths that can be calculated in parallel, and the transcranial correction phase can also be obtained during the calculation. The holography stage is sufficient to simulate the acoustic field in the homogenous intracranial medium after ultrasound transmitting through the skull. The agreement of the k-space corrected pseudospectral time domain method and HASC method was assessed by statistical methods: linear regression between the two methods provided a slope of 0.9735, intercept of 0.0078, and R ² of 0.9982. The Bland-Altman method provided a bias of 0.0015 and 95% limits of agreement 0.065 apart. We demonstrated that the difference in sound intensity at the focal point corrected by HASC and time reversal phase correction method was 0.2% and 0.5% in the results of simulation and experiment, respectively. Not only that, the phase calculation time by the HASC phase correction method can be reduced to 11 min on a multi GPU array, which has clinical potential for ultrasound treatment of brain therapy.

Computationally Efficient Transcranial Ultrasonic Focusing: Taking Advantage of the High Correlation Length of the Human Skull

The phase correction necessary for transcranial ultrasound therapy requires numerical simulation to noninvasively assess the phase shift induced by the skull bone. Ideally, the numerical simulations need to be fast enough for clinical implementation in a brain therapy protocol and to provide accurate estimation of the phase shift to optimize the refocusing through the skull. In this article, we experimentally performed transcranial ultrasound focusing at 900 kHz on N = 5 human skulls. To reduce the computation time, we propose here to perform the numerical simulation at 450 kHz and use the corresponding phase shifts experimentally at 900 kHz. We demonstrate that a 450-kHz simulation restores 94.2% of the pressure when compared with a simulation performed at 900 kHz and 85.0% of the gold standard pressure obtained by an invasive time reversal procedure based on the signal recorded by a hydrophone placed at the target. From a 900- to 450-kHz simulation, the grid size is divided by 8, and the computation time is divided by 10.

Attenuation of the human skull at broadband frequencies by using a carbon nanotube composite photoacoustic transducer

The shockwave generated from a focused carbon nanotube (CNT) composite photoacoustic transducer has a wide frequency band that reaches several MHz in a single pulse. The objective of this study was to measure the transmission characteristics of a shockwave generated by a CNT composite photoacoustic transducer through Asian skulls and compare the results with numerical simulation ones. Three Korean cadaver skulls were used, and five sites were measured for each skull. The average densities and sound speeds of the three skulls were calculated from computed tomography images. The sound pressure after skull penetration was about 11% of the one before skull penetration. High-frequency energy was mostly attenuated. The average attenuation coefficients measured at the five sites of the three skulls were 3.59 ± 0.29, 5.99 ± 1.07, and 3.90 ± 0.86 np/cm/MHz. These values were higher than those previously measured at 270, 836, and 1402 kHz from other groups. The attenuation coefficients simulated by Sim4life were slightly smaller than the experimental values, with similar trends at most sites. The attenuation coefficients varied with measurement sites, skull shape, and thickness. These results may provide important data for future applications of shockwaves in noninvasive neurological treatments.

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

BACKGROUND

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

OBJECTIVES

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

An open-source phase correction toolkit for transcranial focused ultrasound

BACKGROUND

The phase correction on transcranial focused ultrasound is essential to regulate unwanted focal point shift caused by skull bone aberration. The aim of the current study was to design and investigate the feasibility of a ray-based phase correction toolkit for transcranial focused ultrasound.

RESULTS

The peak pressure at focal area was improved by 140.5 ± 7.0% on target I and 134.8 ± 19.1% on target II using proposed phase correction toolkit, respectively. A total computation time of 402.1 ± 24.5 milliseconds was achieved for each sonication.

CONCLUSION

The designed ray-based phase correction software can be used as a lightweight toolkit to compensate aberrated phase within clinical environment.

Virtual clinical trials in medical imaging: a review

The accelerating complexity and variety of medical imaging devices and methods have outpaced the ability to evaluate and optimize their design and clinical use. This is a significant and increasing challenge for both scientific investigations and clinical applications. Evaluations would ideally be done using clinical imaging trials. These experiments, however, are often not practical due to ethical limitations, expense, time requirements, or lack of ground truth. Virtual clinical trials (VCTs) (also known as in silico imaging trials or virtual imaging trials) offer an alternative means to efficiently evaluate medical imaging technologies virtually. They do so by simulating the patients, imaging systems, and interpreters. The field of VCTs has been constantly advanced over the past decades in multiple areas. We summarize the major developments and current status of the field of VCTs in medical imaging. We review the core components of a VCT: computational phantoms, simulators of different imaging modalities, and interpretation models. We also highlight some of the applications of VCTs across various imaging modalities.

A Noninvasive Ultrasound Resonance Method for Detecting Skull Induced Phase Shifts May Provide a Signal for Adaptive Focusing

OBJECTIVE

There may be a need to perform dynamic skull aberration corrections during the non-invasive high-intensity transcranial treatment with magnetic resonance imaging (MRI) -guided focused ultrasound in order to accurately and rapidly restore the focus in the brain.

METHODS

This could possibly be accomplished by using an ultrasound-based correction method based on the skulls' thickness resonance frequencies. The focus of a 500 kHz transducer was centered in the ex vivo human skull caps at different temperatures. The pulse-echoed signals reflected from the skulls were analyzed in the frequency domain to reveal the resonance frequencies for the phase shift calculation. The accuracy was compared to both hydrophone and computed tomography (CT) based analytical methods.

RESULTS

Around 73% of the measurements (n = 784) were in the optimal constructive interference region, with a 15° decrease in the average phase error compared to the previous study. In the best implementation, it performed approximately the same or better than the CT based analytical method currently in clinical use. Linear correlation was found between the resonance frequencies or skull induced phase shifts and the skull temperature with an average rate of -0.4 kHz/°C and 2.6 deg/°C, respectively.

CONCLUSION

The ultrasound based resonance method has shown the feasibility of detecting heating-induced changes of skull phase shift non-invasively and accurately.

SIGNIFICANCE

Since the technique can be made MRI compatible and integrated in the therapy arrays, it may enable temperature tracking and adaptive focusing during high-intensity transcranial ultrasound treatments, to prevent skull overheating and preserve the transcranial focusing integrity.

Super-Resolution Imaging Through the Human Skull

High-resolution transcranial ultrasound imaging in humans has been a persistent challenge for ultrasound due to the imaging degradation effects from aberration and reverberation. These mechanisms depend strongly on skull morphology and have high variability across individuals. Here, we demonstrate the feasibility of human transcranial super-resolution imaging using a geometrical focusing approach to efficiently concentrate energy at the region of interest, and a phase correction focusing approach that takes the skull morphology into account. It is shown that using the proposed focused super-resolution method, we can image a 208- [Formula: see text] microtube behind a human skull phantom in both an out-of-plane and an in-plane configuration. Individual phase correction profiles for the temporal region of the human skull were calculated and subsequently applied to transmit-receive a custom focused super-resolution imaging sequence through a human skull phantom, targeting the 208- [Formula: see text] diameter microtube at 68.5 mm in depth and at 2.5 MHz. Microbubble contrast agents were diluted to a concentration of 1.6×10⁶ bubbles/mL and perfused through the microtube. It is shown that by correcting for the skull aberration, the RF signal amplitude from the tube improved by a factor of 1.6 in the out-of-plane focused emission case. The lateral registration error of the tube's position, which in the uncorrected case was 990 [Formula: see text], was reduced to as low as 50 [Formula: see text] in the corrected case as measured in the B-mode images. Sensitivity in microbubble detection for the phase-corrected case increased by a factor of 1.48 in the out-of-plane imaging case, while, in the in-plane target case, it improved by a factor of 1.31 while achieving an axial registration correction from an initial 1885- [Formula: see text] error for the uncorrected emission, to a 284- [Formula: see text] error for the corrected counterpart. These findings suggest that super-resolution imaging may be used far more generally as a clinical imaging modality in the brain.

On the Evaluation of the Suitability of the Materials Used to 3D Print Holographic Acoustic Lenses to Correct Transcranial Focused Ultrasound Aberrations

The correction of transcranial focused ultrasound aberrations is a relevant topic for enhancing various non-invasive medical treatments. Presently, the most widely accepted method to improve focusing is the emission through multi-element phased arrays; however, a new disruptive technology, based on 3D printed holographic acoustic lenses, has recently been proposed, overcoming the spatial limitations of phased arrays due to the submillimetric precision of the latest generation of 3D printers. This work aims to optimize this recent solution. Particularly, the preferred acoustic properties of the polymers used for printing the lenses are systematically analyzed, paying special attention to the effect of p-wave speed and its relationship to the achievable voxel size of 3D printers. Results from simulations and experiments clearly show that, given a particular voxel size, there are optimal ranges for lens thickness and p-wave speed, fairly independent of the emitted frequency, the transducer aperture, or the transducer-target distance.

Simulation of nonlinear trans-skull focusing and formation of shocks in brain using a fully populated ultrasound array with aberration correction

Multi-element high-intensity focused ultrasound phased arrays in the shape of hemispheres are currently used in clinics for thermal lesioning in deep brain structures. Certain side effects of overheating non-targeted tissues and skull bones have been revealed. Here, an approach is developed to mitigate these effects. A specific design of a fully populated 256-element 1-MHz array shaped as a spherical segment (F-number, F_# = 1) and filled by randomly distributed equal-area polygonal elements is proposed. Capability of the array to generate high-amplitude shock fronts at the focus is tested in simulations by combining three numerical algorithms for linear and nonlinear field modeling and aberration correction. The algorithms are based on the combination of the Rayleigh integral, a linear pseudo-spectral time domain Kelvin-Voigt model, and nonlinear Westervelt model to account for the effects of inhomogeneities, aberrations, reflections, absorption, nonlinearity, and shear waves in the skull. It is shown that the proposed array can generate nonlinear waveforms with shock amplitudes >60 MPa at the focus deep inside the brain without exceeding the existing technical limitation on the intensity of 40 W/cm² at the array elements. Such shock amplitudes are sufficient for mechanical ablation of brain tissues using the boiling histotripsy approach and implementation of other shock-based therapies.

Virtual Brain Projection for Evaluating Trans-skull Beam Behavior of Transcranial Ultrasound Devices

Focused ultrasound single-element piezoelectric transducers constitute a promising method to deliver ultrasound to the brain in low-intensity applications, but are subject to defocusing and high attenuation because of transmission through the skull. Here, a novel virtual brain projection method is used to superimpose a magnetic resonance image of the brain in ex vivo human skulls to provide targets during trans-skull focused ultrasound single-element piezoelectric transducer pressure field mapping. Positions of the transducer, skull and hydrophone are tracked in real time using a stereoscopic navigation camera and 3-D Slicer software. Virtual targets of the left dorsolateral prefrontal cortex, left hippocampus and cerebellar vermis were chosen to illustrate the method's flexibility in evaluating focal-zone beam distortion and attenuation. The regions are of interest as non-invasive brain stimulation targets in the treatment of neuropsychiatric disorders via repeated ultrasound exposure. The technical approach can facilitate the assessment of transcranial ultrasound device operator positioning reliability, intracranial beam behavior and computational model validation.

A rapid beam simulation framework for transcranial focused ultrasound

Transcranial focused ultrasound is a non-invasive therapeutic modality that can be used to treat essential tremor. Beams of energy are focused into a small spot in the thalamus, resulting in tissue heating and ablation. Here, we report on a rapid 3D numeric simulation framework that can be used to predict focal spot characteristics prior to the application of ultrasound. By comparing with magnetic resonance proton resonance frequency shift thermometry (MR thermometry) data acquired during treatments of essential tremor, we verified that our simulation framework can be used to predict focal spot position, and with patient-specific calibration, predict focal spot temperature rise. Preliminary data suggests that lateral smearing of the focal spot can be simulated. The framework may also be relevant for other therapeutic ultrasound applications such as blood brain barrier opening and neuromodulation.

Accumulated thermal dose in MRI-guided focused ultrasound for essential tremor: repeated sonications with low focal temperatures

OBJECTIVE

The object of this study was to correlate lesion size with accumulated thermal dose (ATD) in transcranial MRI-guided focused ultrasound (MRgFUS) treatments of essential tremor with focal temperatures limited to 50°C-54°C.

METHODS

Seventy-five patients with medically refractory essential tremor underwent MRgFUS thalamotomy at the authors' institution. Intraoperative MR thermometry was performed to measure the induced temperature and thermal dose distributions (proton resonance frequency shift coefficient = -0.00909 ppm/°C). In 19 patients, it was not possible to raise the focal temperature above 54°C because of unfavorable skull characteristics and/or the pain associated with cranial heating. In this patient subset, sonications with focal temperatures between 50°C and 54°C were repeated (5.1 ± 1.5, mean ± standard deviation) to accumulate a sufficient thermal dose for lesion formation. The ATD profile sizes (17, 40, 100, 200, and 240 cumulative equivalent minutes at 43°C [CEM43]) calculated by combining axial MR thermometry data from individual sonications were correlated with the corresponding lesion sizes measured on axial T1-weighted (T1w) and T2-weighted (T2w) MR images acquired 1 day posttreatment. Manual corrections were applied to the MR thermometry data prior to thermal dose accumulation to compensate for off-resonance-induced spatial-shifting artifacts.

RESULTS

Mean lesion sizes measured on T2w MRI (5.0 ± 1.4 mm) were, on average, 28% larger than those measured on T1w MRI (3.9 ± 1.4 mm). The ATD thresholds found to provide the best correlation with lesion sizes measured on T2w and T1w MRI were 100 CEM43 (regression slope = 0.97, R2 = 0.66) and 200 CEM43 (regression slope = 0.98, R2 = 0.89), respectively, consistent with data from a previous study of MRgFUS thalamotomy via repeated sonications at higher focal temperatures (≥ 55°C). Two-way linear mixed-effects analysis revealed that dominant tremor subscores on the Fahn-Tolosa-Marin Clinical Rating Scale for Tremor (CRST) were statistically different from baseline at 3 months and 1 year posttreatment in both low-temperature (50°C-54°C) and high-temperature (≥ 55°C) patient cohorts. No significant fixed effect on the dominant tremor scores was found for the temperature cohort factor.

CONCLUSIONS

In transcranial MRgFUS thalamotomy for essential tremor, repeated sonications with focal temperatures between 50°C and 54°C can accumulate a sufficient thermal dose to generate lesions for clinically relevant tremor suppression up to 1 year posttreatment, and the ATD can be used to predict the size of the resulting ablation zones measured on MRI. These data will serve to guide future clinical MRgFUS brain procedures, particularly those in which focal temperatures are limited to below 55°C.

Steering Capabilities of an Acoustic Lens for Transcranial Therapy: Numerical and Experimental Studies

For successful brain therapy, transcranial focused ultrasound must compensate for the time shifts induced locally by the skull. The patient-specific phase profile is currently generated by multi-element arrays which, over time, have tended toward increasing element count. We recently introduced a new approach, consisting of a single-element transducer coupled to an acoustic lens of controlled thickness. By adjusting the local thickness of the lens, we were able to induce phase differences which compensated those induced by the skull. Nevertheless, such an approach suffers from an apparent limitation: the lens is a priori designed for one specific target. In this paper, we demonstrate the possibility of taking advantage of the isoplanatic angle of the aberrating skull in order to steer the focus by mechanically moving the transducer/acoustic lens pair around its initial focusing position. This study, conducted on three human skull samples, demonstrates that tilting of the transducer with the lens restores a single -3 dB focal volume at 914 kHz for a steering up to ±11 mm in the transverse direction, and ±10 mm in the longitudinal direction, around the initial focal region.

Enhanced Numerical Method for the Design of 3-D-Printed Holographic Acoustic Lenses for Aberration Correction of Single-Element Transcranial Focused Ultrasound

The correction of transcranial focused ultrasound aberrations is a relevant issue for enhancing various non-invasive medical treatments. The emission through multi-element phased arrays has been the most widely accepted method to improve focusing in recent years; however, the number and size of transducers represent a bottleneck that limits the focusing accuracy of the technique. To overcome this limitation, a new disruptive technology, based on 3-D-printed acoustic lenses, has recently been proposed. As the submillimeter precision of the latest generation of 3-D printers has been proven to overcome the spatial limitations of phased arrays, a new challenge is to improve the accuracy of the numerical simulations required to design this type of ultrasound lens. In the study described here, we evaluated two improvements in the numerical model applied in previous works for the design of 3-D-printed lenses: (i) allowing the propagation of shear waves in the skull by means of its simulation as an isotropic solid and (ii) introduction of absorption into the set of equations that describes the dynamics of the wave in both fluid and solid media. The results obtained in the numerical simulations are evidence that the inclusion of both s-waves and absorption significantly improves focusing.

Modeling ultrasound propagation through material of increasing geometrical complexity

Ultrasound is increasingly being recognized as a neuromodulatory and therapeutic tool, inducing a broad range of bio-effects in the tissue of experimental animals and humans. To achieve these effects in a predictable manner in the human brain, the thick cancellous skull presents a problem, causing attenuation. In order to overcome this challenge, as a first step, the acoustic properties of a set of simple bone-modeling resin samples that displayed an increasing geometrical complexity (increasing step sizes) were analyzed. Using two Non-Destructive Testing (NDT) transducers, we found that Wiener deconvolution predicted the Ultrasound Acoustic Response (UAR) and attenuation caused by the samples. However, whereas the UAR of samples with step sizes larger than the wavelength could be accurately estimated, the prediction was not accurate when the sample had a smaller step size. Furthermore, a Finite Element Analysis (FEA) performed in ANSYS determined that the scattering and refraction of sound waves was significantly higher in complex samples with smaller step sizes compared to simple samples with a larger step size. Together, this reveals an interaction of frequency and geometrical complexity in predicting the UAR and attenuation. These findings could in future be applied to poro-visco-elastic materials that better model the human skull.

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

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

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

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

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

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

Measurements of the Relationship Between CT Hounsfield Units and Acoustic Velocity and How It Changes With Photon Energy and Reconstruction Method

Transcranial magnetic resonance-guided focused ultrasound continues to gain traction as a noninvasive treatment option for a variety of pathologies. Focusing ultrasound through the skull can be accomplished by adding a phase correction to each element of a hemispherical transducer array. The phase corrections are determined with acoustic simulations that rely on speed of sound estimates derived from CT scans. While several studies have investigated the relationship between acoustic velocity and CT Hounsfield units (HUs), these studies have largely ignored the impact of X-ray energy, reconstruction method, and reconstruction kernel on the measured HU, and therefore the estimated velocity, and none have measured the relationship directly. In this paper, 91 ex vivo human skull fragments from two skulls are imaged by 80 CT scans with a variety of energies and reconstruction methods. The average HU from each fragment is found for each scan and correlated with the speed of sound measured using a through transmission technique in that fragment. As measured by the -squared value, the results show that CT is able to account for 23%-53% of the variation in velocity in the human skull. Both the X-ray energy and the reconstruction technique significantly alter the -squared value and the linear relationship between HU and speed of sound in bone. Accounting for these variations will lead to more accurate phase corrections and more efficient transmission of acoustic energy through the skull.

The reduction in treatment efficiency at high acoustic powers during MR-guided transcranial focused ultrasound thalamotomy for Essential Tremor

PURPOSE

To analyze clinical data indicating a reduction in the induced energy-temperature efficiency relationship during transcranial focused ultrasound (FUS) Essential Tremor (ET) thalamotomy treatments at higher acoustic powers, establish its relationship with the spatial distribution of the focal temperature elevation, and explore its cause.

METHODS

A retrospective observational study of patients (n = 19) treated between July 2015 and August 2016 for (ET) by FUS thalamotomy was performed. These data were analyzed to compare the relationships between the applied power, the applied energy, the resultant peak temperature achieved in the brain, and the dispersion of the focal volume. Full ethics approval was received and all patients provided signed informed consent forms before the initiation of the study. Computer simulations, animal experiments, and clinical system tests were performed to determine the effects of skull heating, changes in brain properties and transducer acoustic output, respectively. All animal procedures were approved by the Animal Care and Use Committee and conformed to the guidelines set out by the Canadian Council on Animal Care. MATLAB was used to perform statistical analysis.

RESULTS

The reduction in the energy efficiency relationship during treatment correlates with the increase in size of the focal volume at higher sonication powers. A linear relationship exists showing that a decrease in treatment efficiency correlates positively with an increase in the focal size over the course of treatment (P < 0.01), supporting the hypothesis of transient skull and tissue heating causing acoustic aberrations leading to a decrease in efficiency. Changes in thermal conductivity, perfusion, absorption rates in the brain, as well as ultrasound transducer acoustic output levels were found to have minimal effects on the observed reduction in efficiency.

CONCLUSIONS

The reduction in energy-temperature efficiency during high-power FUS treatments correlated with observed increases in the size of the focal volume and is likely caused by transient changes in the tissue and skull during heating.

3D-printed adaptive acoustic lens as a disruptive technology for transcranial ultrasound therapy using single-element transducers

The development of multi-element arrays for better control of the shape of ultrasonic beams has opened the way for focusing through highly aberrating media, such as the human skull. As a result, the use of brain therapy with transcranial-focused ultrasound has rapidly grown. Although effective, such technology is expensive. We propose a disruptive, low-cost approach that consists of focusing a 1 MHz ultrasound beam through a human skull with a single-element transducer coupled with a tailored silicone acoustic lens cast in a 3D-printed mold and designed using computed tomography-based numerical acoustic simulation. We demonstrate on N  =  3 human skulls that adding lens-based aberration correction to a single-element transducer increases the deposited energy on the target 10 fold.

Catheter Hydrophone Aberration Correction for Transcranial Histotripsy Treatment of Intracerebral Hemorrhage: Proof-of-Concept

Histotripsy is a minimally invasive ultrasound therapy that has shown rapid liquefaction of blood clots through human skullcaps in an in vitro intracerebral hemorrhage model. However, the efficiency of these treatments can be compromised if the skull-induced aberrations are uncorrected. We have developed a catheter hydrophone which can perform aberration correction (AC) and drain the liquefied clot following histotripsy treatment. Histotripsy pulses were delivered through an excised human skullcap using a 256-element, 500-kHz hemisphere array transducer with a 15-cm focal distance. A custom hydrophone was fabricated using a mm PZT-5h crystal interfaced to a coaxial cable and integrated into a drainage catheter. An AC algorithm was developed to correct the aberrations introduced between histotripsy pulses from each array element. An increase in focal pressure of up to 60% was achieved at the geometric focus and 27%-62% across a range of electronic steering locations. The sagittal and axial -6-dB beam widths decreased from 4.6 to 2.2 mm in the sagittal direction and 8 to 4.4 mm in the axial direction, compared to 1.5 and 3 mm in the absence of aberration. After performing AC, lesions with diameters ranging from 0.24 to 1.35 mm were generated using electronic steering over a mm grid in a tissue-mimicking phantom. An average volume of 4.07 ± 0.91 mL was liquefied and drained after using electronic steering to treat a 4.2-mL spherical volume in in vitro bovine clots through the skullcap.