Megahertz rate, volumetric imaging of bubble clouds in sonothrombolysis using a sparse hemispherical receiver array

A Convolutional Neural Network for Beamforming and Image Reconstruction in Passive Cavitation Imaging

Convolutional neural networks (CNNs), initially developed for image processing applications, have recently received significant attention within the field of medical ultrasound imaging. In this study, passive cavitation imaging/mapping (PCI/PAM), which is used to map cavitation sources based on the correlation of signals across an array of receivers, is evaluated. Traditional reconstruction techniques in PCI, such as delay-and-sum, yield high spatial resolution at the cost of a substantial computational time. This results from the resource-intensive process of determining sensor weights for individual pixels in these methodologies. Consequently, the use of conventional algorithms for image reconstruction does not meet the speed requirements that are essential for real-time monitoring. Here, we show that a three-dimensional (3D) convolutional network can learn the image reconstruction algorithm for a 16×16 element matrix probe with a receive frequency ranging from 256 kHz up to 1.0 MHz. The network was trained and evaluated using simulated data representing point sources, resulting in the successful reconstruction of volumetric images with high sensitivity, especially for single isolated sources (100% in the test set). As the number of simultaneous sources increased, the network's ability to detect weaker intensity sources diminished, although it always correctly identified the main lobe. Notably, however, network inference was remarkably fast, completing the task in approximately 178 s for a dataset comprising 650 frames of 413 volume images with signal duration of 20μs. This processing speed is roughly thirty times faster than a parallelized implementation of the traditional time exposure acoustics algorithm on the same GPU device. This would open a new door for PCI application in the real-time monitoring of ultrasound ablation.

Passive Cavitation Imaging Artifact Reduction Using Data-Adaptive Spatial Filtering

Passive cavitation imaging (PCI) with a clinical diagnostic array results in poor axial localization of bubble activity due to the size of the point spread function (PSF). The objective of this study was to determine if data-adaptive spatial filtering improved PCI beamforming performance relative to standard frequency-domain delay, sum, and integrate (DSI) or robust Capon beamforming (RCB). The overall goal was to improve source localization and image quality without sacrificing computation time. Spatial filtering was achieved by applying a pixel-based mask to DSI- or RCB-beamformed images. The masks were derived from DSI, RCB, or phase or amplitude coherence factors (ACFs) using both receiver operating characteristic (ROC) and precision-recall (PR) curve analyses. Spatially filtered passive cavitation images were formed from cavitation emissions based on two simulated sources densities and four source distribution patterns mimicking cavitation emissions induced by an EkoSonic catheter. Beamforming performance was assessed via binary classifier metrics. The difference in sensitivity, specificity, and area under the ROC curve (AUROC) differed by no more than 11% across all algorithms for both source densities and all source patterns. The computational time required for each of the three spatially filtered DSIs was two orders of magnitude less than that required for time-domain RCB and thus this data-adaptive spatial filtering strategy for PCI beamforming is preferable given the similar binary classification performance.

Tri-modality cavitation mapping in shock wave lithotripsy

Shock wave lithotripsy (SWL) has been widely used for non-invasive treatment of kidney stones. Cavitation plays an important role in stone fragmentation, yet it may also contribute to renal injury during SWL. It is therefore crucial to determine the spatiotemporal distributions of cavitation activities to maximize stone fragmentation while minimizing tissue injury. Traditional cavitation detection methods include high-speed optical imaging, active cavitation mapping (ACM), and passive cavitation mapping (PCM). While each of the three methods provides unique information about the dynamics of the bubbles, PCM has most practical applications in biological tissues. To image the dynamics of cavitation bubble collapse, we previously developed a sliding-window PCM (SW-PCM) method to identify each bubble collapse with high temporal and spatial resolution. In this work, to further validate and optimize the SW-PCM method, we have developed tri-modality cavitation imaging that includes three-dimensional high-speed optical imaging, ACM, and PCM seamlessly integrated in a single system. Using the tri-modality system, we imaged and analyzed laser-induced single cavitation bubbles in both free field and constricted space and shock wave-induced cavitation clusters. Collectively, our results have demonstrated the high reliability and spatial-temporal accuracy of the SW-PCM approach, which paves the way for the future in vivo applications on large animals and humans in SWL.

Real-Time Passive Acoustic Mapping Using Sparse Matrix Multiplication

Passive acoustic mapping enables the spatiotemporal monitoring of cavitation with circulating microbubbles during focused ultrasound (FUS)-mediated blood-brain barrier opening. However, the computational load for processing large data sets of cavitation maps or more complex algorithms limit the visualization in real-time for treatment monitoring and adjustment. In this study, we implemented a graphical processing unit (GPU)-accelerated sparse matrix-based beamforming and time exposure acoustics in a neuronavigation-guided ultrasound system for real-time spatiotemporal monitoring of cavitation. The system performance was tested in silico through benchmarking, in vitro using nonhuman primate (NHP) and human skull specimens, and demonstrated in vivo in NHPs. We demonstrated the stability of the cavitation map for integration times longer than 62.5 [Formula: see text]. A compromise between real-time displaying and cavitation map quality obtained from beamformed RF data sets with a size of 2000 ×128 ×30 (axial [Formula: see text]) was achieved for an integration time of [Formula: see text], which required a computational time of 0.27 s (frame rate of 3.7 Hz) and could be displayed in real-time between pulses at PRF = 2 Hz. Our benchmarking tests show that the GPU sparse-matrix algorithm processed the RF data set at a computational rate of [Formula: see text]/pixel/sample, which enables adjusting the frame rate and the integration time as needed. The neuronavigation system with real-time implementation of cavitation mapping facilitated the localization of the cavitation activity and helped to identify distortions due to FUS phase aberration. The in vivo test of the method demonstrated the feasibility of GPU-accelerated sparse matrix computing in a close to a clinical condition, where focus distortions exemplify problems during treatment. These experimental conditions show the need for spatiotemporal monitoring of cavitation with real-time capability that enables the operator to correct or halt the sonication in case substantial aberrations are observed.

Doppler Passive Acoustic Mapping

In therapeutic ultrasound using microbubbles, it is essential to drive the microbubbles into the correct type of activity and the correct location to produce the desired biological response. Although passive acoustic mapping (PAM) is capable of locating where microbubble activities are generated, it is well known that microbubbles rapidly move within the ultrasound beam. We propose a technique that can image microbubble movement by estimating their velocities within the focal volume. Microbubbles embedded within a wall-less channel of a tissue-mimicking material were sonicated using 1-MHz focused ultrasound. The acoustic emissions generated by the microbubbles were captured with a linear array (L7-4). PAM with robust Capon beamforming was used to localize the microbubble acoustic emissions. We spectrally analyzed the time trace of each position and isolated the higher harmonics. Microbubble velocity maps were constructed from the position-dependent Doppler shifts at different time points during sonication. Microbubbles moved primarily away from the transducer at velocities on the order of 1 m/s due to primary acoustic radiation forces, producing a time-dependent velocity distribution. We detected microbubble motion both away and toward the receiving array, revealing the influence of acoustic radiation forces and fluid motion due to the ultrasound exposure. High-speed optical images confirmed the acoustically measured microbubble velocities. Doppler PAM enables passive estimation of microbubble motion and may be useful in therapeutic applications, such as drug delivery across the blood-brain barrier, sonoporation, sonothrombolysis, and drug release.

Non-invasive recanalization of deep venous thrombosis by high frequency ultrasound in a swine model with follow-up

AIMS

Pulsed cavitational ultrasound therapy (thombotripsy) allows the accurate fractionation of a distant thrombus. We aimed to evaluate the efficacy and safety of non-invasive thrombotripsy using a robotic assisted and high frequency ultrasound approach to recanalize proximal deep venous thrombosis (DVT) in a swine model.

METHODS

Occlusive thrombosis was obtained with a dual jugular and femoral endoveinous approach. The therapeutic device was composed of a 2.25 MHz focused transducer centered by a linear ultrasound probe, and a robotic arm. The feasibility, security, and efficacy (venous channel patency) assessment after thrombotripsy was performed on 13 pigs with acute occluded DVT. To assess the mid-term efficacy of this technique, 8 pigs were followed up for 14 days after thrombotripsy and compared with 8 control pigs. The primary efficacy endpoint was the venous patency. Safety was assessed by the search for local vessel wall injury and pulmonary embolism.

RESULTS

We succeeded in treating all pigs except two with no accessible femoral vein. After median treatment duration of 23 minutes of cavitation, all treated DVT were fully recanalized acutely. At 14 days, in the treated group, six of the eight pigs had a persistent patent vein and two pigs had a venous reocclusion. In the control group all pigs had a persistent venous occlusion. At sacrifice, no local vein nor arterial wall damage were observed as well as no evidence of pulmonary embolism in all pigs.

CONCLUSION

High frequency thrombotripsy seems to be effective and safe for non-invasive venous recanalization of DVT.

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

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

Ultrasound-Responsive Cavitation Nuclei for Therapy and Drug Delivery

Therapeutic ultrasound strategies that harness the mechanical activity of cavitation nuclei for beneficial tissue bio-effects are actively under development. The mechanical oscillations of circulating microbubbles, the most widely investigated cavitation nuclei, which may also encapsulate or shield a therapeutic agent in the bloodstream, trigger and promote localized uptake. Oscillating microbubbles can create stresses either on nearby tissue or in surrounding fluid to enhance drug penetration and efficacy in the brain, spinal cord, vasculature, immune system, biofilm or tumors. This review summarizes recent investigations that have elucidated interactions of ultrasound and cavitation nuclei with cells, the treatment of tumors, immunotherapy, the blood-brain and blood-spinal cord barriers, sonothrombolysis, cardiovascular drug delivery and sonobactericide. In particular, an overview of salient ultrasound features, drug delivery vehicles, therapeutic transport routes and pre-clinical and clinical studies is provided. Successful implementation of ultrasound and cavitation nuclei-mediated drug delivery has the potential to change the way drugs are administered systemically, resulting in more effective therapeutics and less-invasive treatments.

Correlation Between Brain Tissue Damage and Inertial Cavitation Dose Quantified Using Passive Cavitation Imaging

Focused ultrasound (FUS)-induced cavitation-mediated brain therapies have become emerging therapeutic modalities for neurologic diseases. Cavitation monitoring is essential to ensure the safety of all cavitation-mediated therapeutic techniques as inertial cavitation can be associated with tissue damage. The objective of this study was to reveal the correlation between the inertial cavitation dose, quantified by passive cavitation imaging (PCI), and brain tissue histologic-level damage induced by FUS in combination with microbubbles. An ultrasound image-guided FUS system consisting of a single-element FUS transducer (1.5 MHz) and a co-axially aligned 128-element linear ultrasound imaging array was used to perform FUS treatment of mice. Mice were sonicated by FUS with different peak negative pressures (0.5 MPa, 1.1 MPa, 4.0 MPa and 6.5 MPa) in the presence of systemically injected microbubbles. The acoustic emissions from the FUS-activated microbubbles were passively detected by the imaging array. The pre-beamformed channel data were acquired and processed offline using the frequency-domain delay, sum and integration algorithm to generate inertial cavitation maps. All the mice were sacrificed after the FUS treatment, and their brains were harvested and processed for hematoxylin and eosin staining. The obtained inertial cavitation maps revealed the dynamic changes of microbubble behaviors during FUS treatment at different pressure levels. It was found that the inertial cavitation dose quantified based on PCI had a linear correlation with the scale of histologic-level tissue damage. Findings from this study suggested that PCI can be used to predict histologic-level tissue damage associated with the FUS-induced cavitation.

Comparison study of passive acoustic mapping and high-speed photography for monitoring in situ cavitation bubbles

The spatiotemporal accuracy of passive acoustic mapping (PAM) for monitoring in situ cavitation bubbles has not been assessed directly via optical means. Here, the cavitation bubbles are monitored from two image sequences obtained simultaneously with PAM and high-speed photography (HSP). The temporal accuracy of PAM for detecting cavitation nucleation and the spatial resolution for cavitation localization are compared with those measured from HSP. The results show that PAM has a temporal accuracy of 20 μs. Mean differences in the spatial locations of PAM and HSP are as small as 10.0 and 30.5 μm along the lateral and axial directions, respectively.

Advances in acoustic monitoring and control of focused ultrasound-mediated increases in blood-brain barrier permeability

Transcranial focused ultrasound (FUS) combined with intravenously circulating microbubbles can transiently and selectively increase blood-brain barrier permeability to enable targeted drug delivery to the central nervous system, and is a technique that has the potential to revolutionize the way neurological diseases are managed in medical practice. Clinical testing of this approach is currently underway in patients with brain tumors, early Alzheimer's disease, and amyotrophic lateral sclerosis. A major challenge that needs to be addressed in order for widespread clinical adoption of FUS-mediated blood-brain barrier permeabilization to occur is the development of systems and methods for real-time treatment monitoring and control, to ensure that safe and effective acoustic exposure levels are maintained throughout the procedures. This review gives a basic overview of the oscillation dynamics, acoustic emissions, and biological effects associated with ultrasound-stimulated microbubbles in vivo, and provides a summary of recent advances in acoustic-based strategies for detecting, controlling, and mapping microbubble activity in the brain. Further development of next-generation clinical FUS brain devices tailored towards microbubble-mediated applications is warranted and required for translation of this potentially disruptive technology into routine clinical practice.

Cavitation dose painting for focused ultrasound-induced blood-brain barrier disruption

Focused ultrasound combined with microbubble for blood-brain barrier disruption (FUS-BBBD) is a promising technique for noninvasive and localized brain drug delivery. This study demonstrates that passive cavitation imaging (PCI) is capable of predicting the location and concentration of nanoclusters delivered by FUS-BBBD. During FUS-BBBD treatment of mice, the acoustic emissions from FUS-activated microbubbles were passively detected by an ultrasound imaging system and processed offline using a frequency-domain PCI algorithm. After the FUS treatment, radiolabeled gold nanoclusters, ⁶⁴Cu-AuNCs, were intravenously injected into the mice and imaged by positron emission tomography/computed tomography (PET/CT). The centers of the stable cavitation dose (SCD) maps obtained by PCI and the corresponding centers of the ⁶⁴Cu-AuNCs concentration maps obtained by PET coincided within 0.3 ± 0.4 mm and 1.6 ± 1.1 mm in the transverse and axial directions of the FUS beam, respectively. The SCD maps were found to be linearly correlated with the ⁶⁴Cu-AuNCs concentration maps on a pixel-by-pixel level. These findings suggest that SCD maps can spatially "paint" the delivered nanocluster concentration, a technique that we named as cavitation dose painting. This PCI-based cavitation dose painting technique in combination with FUS-BBBD opens new horizons in spatially targeted and modulated brain drug delivery.

Receiver array design for sonothrombolysis treatment monitoring in deep vein thrombosis

High intensity focused ultrasound (HIFU) can disintegrate blood clots through the generation and stimulation of bubble clouds within thrombi. This work examined the design of a device to image bubble clouds for monitoring cavitation-based HIFU treatments of deep vein thrombosis (DVT). Acoustic propagation simulations were carried out on multi-layered models of the human thigh using two patient data sets from the Visible Human Project. The design considerations included the number of receivers (32, 64, 128, 256, and 512), their spatial positioning, and the effective angular array aperture (100° and 180° about geometric focus). Imaging array performance was evaluated for source frequencies of 250, 750, and 1500 kHz. Receiver sizes were fixed relative to the wavelength (pistons, diameter  =  λ/2) and noise was added at levels that scaled with receiver area. With a 100° angular aperture the long axis size of the  -3 dB main lobe was ~1.2λ-i.e. on the order of the vessel diameter at 250 kHz (~7 mm). Increasing the array aperture to span 180° about the geometric focus reduced the long axis by a factor of ~2. The smaller main lobe sizes achieved by imaging at higher frequencies came at the cost of increased levels of sensitivity to phase aberrations induced during acoustic propagation through the intervening soft tissue layers. With noise added to receiver signals, images could be reconstructed with peak sidelobe ratios  <  -3 dB using single-cycle integration times for source frequencies of 250 and 750 kHz (NRx  ⩾  128). At 1500 kHz, longer integration times and/or higher element counts were required to achieve similar peak sidelobe ratios. Our results suggest that a modest number of receivers(i.e. NRx  =  128) arranged on a semi-cylindrical shell may be sufficient to enable passive acoustic imaging with single-cycle integration times (i.e. volumetric rates up to 0.75 MHz) for monitoring cavitation-based HIFU treatments of DVT.