Using passive cavitation images to classify high-intensity focused ultrasound lesions

Real-time passive cavitation mapping and B-mode fusion imaging via hybrid adaptive beamformer with modified diagnostic ultrasound platform

The implementation of real-time, convenient and high-resolution passive cavitation imaging (PCM) is crucial for ensuring the safety and effectiveness of ultrasound applications related to cavitation effects. However, the current B-mode ultrasound imaging system cannot achieve these functions. By developing a hybrid adaptive beamforming algorithm, the current work presented a real-time PCM and B-mode fusion imaging technique, using a modified diagnostic ultrasound platform enabling time-division multiplexing external triggering function. The proposed hybrid adaptive beamformer combined the advantages of delay-multiply-and-sum (DMAS) and minimum variance (MV) methods to effectively suppress the side lobe and tail-like artifacts, improving the resolution of PCM images. A high-pass filter was applied to selectively detect cavitation-specific signals while removing the interference from the tissue scatters. The system enabled synchronous visualization of tissue structure and cavitation activity under ultrasound exposure. Both numerical and experimental studies demonstrated that, compared with DAS, MV-DAS and DMAS methods, the proposed MV-DMAS algorithm performed better in both axial and lateral resolutions. This work represented a significant advancement in achieving high-quality real-time B-mode and PCM fusion imaging utilizing commercial medical ultrasound system, providing a powerful tool for synchronous monitoring and manipulating cavitation activity, which would enhance the safety and efficacy of cavitation-based applications.

Effect of pulse parameters on ablation efficiency in dual-frequency HIFU therapy

High-intensity focused ultrasound (HIFU) has now been widely used to ablate various benign and malignant tumors. But it is still critical to increase the ablation efficiency in many clinical applications. Dual-frequency HIFU has been proven to be more efficient in ablation, but the principle on selecting the pulse parameters in this method remains to be explored. In this study, the in vitro lesion areas under different pulse repetition frequencies (PRFs), duty ratios, and frequency differences were compared, cavitation activity was also monitored during HIFU exposure. The results showed that different pulse parameters caused different types of lesions. In HIFU therapy, those pulse parameters that maximize the thermal effect, reduce heat dissipation and generate sufficient cavitation activities should be considered. But the method of evaluating or predicting the damage by using the cavitation dose is only applicable to mechanical damage.

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.

Monitoring cavitation dynamics evolution in tissue mimicking hydrogels for repeated exposures via acoustic cavitation emissions

A 700 kHz histotripsy array is used to generate repeated cavitation events in agarose, gelatin, and polyacrylamide hydrogels. High-speed optical imaging, a broadband hydrophone, and the narrow-band receive elements of the histotripsy array are used to capture bubble dynamics and acoustic cavitation emissions. Bubble radii, lifespan, shockwave amplitudes are noted to be measured in close agreement between the different observation methods. These features also decrease with increasing hydrogel stiffness for all of the tested materials. However, the evolutions of these properties during the repeated irradiations vary significantly across the different material subjects. Bubble maximum radius initially increases, then plateaus, and finally decreases in agarose, but remains constant across exposures in gelatin and polyacrylamide. The bubble lifespan increases monotonically in agarose and gelatin but decreases in polyacrylamide. Collapse shockwave amplitudes were measured to have different-shaped evolutions between all three of the tested materials. Bubble maximum radii, lifespans, and collapse shockwave amplitudes were observed to express evolutions that are dependent on the structure and stiffness of the nucleation medium.

High-spatial-resolution, instantaneous passive cavitation imaging with temporal resolution in histotripsy: a simulation study

Purpose

In histotripsy, a shock wave is transmitted, and the resulting inertial bubble cavitation that disrupts tissue is used for treatment. Therefore, it is necessary to detect when cavitation occurs and track the position of cavitation occurrence using a new passive cavitation (PC) imaging method.

Methods

An integrated PC image, which is constructed by collecting the focused signals at all times, does not provide information on when cavitation occurs and has poor spatial resolution. To solve this problem, we constructed instantaneous PC images by applying delay and sum beamforming at instantaneous time instants. By calculating instantaneous PC images at all data acquisition times, the proposed method can detect cavitation when it occurs by using the property that when signals from the cavitation are focused, their amplitude becomes large, and it can obtain a high-resolution PC image by masking out side lobes in the vicinity of cavitation.

Results

Ultrasound image simulation confirmed that the proposed method has higher resolution than conventional integrated PC imaging and showed that it can determine the position and time of cavitation occurrence as well as the signal strength.

Conclusion

Since the proposed novel PC imaging method can detect each cavitation separately when the incidence of cavitations is low, it can be used to monitor the treatment process of shock wave therapy and histotripsy, in which cavitation is an important mechanism of treatment.

Effect of Overpressure on Acoustic Emissions and Treated Tissue Histology in ex Vivo Bulk Ultrasound Ablation

Bulk ultrasound ablation is a thermal therapy approach in which tissue is heated by unfocused or weakly focused sonication (average intensities on the order of 100 W/cm²) to achieve coagulative necrosis within a few minutes exposure time. Assessing the role of bubble activity, including acoustic cavitation and tissue vaporization, in bulk ultrasound ablation may help in making bulk ultrasound ablation safer and more effective for clinical applications. Here, two series of ex vivo ablation trials were conducted to investigate the role of bubble activity and tissue vaporization in bulk ultrasound ablation. Fresh bovine liver tissue was ablated with unfocused, continuous-wave ultrasound using ultrasound image-ablate arrays sonicating at 31 W/cm² (0.9 MPa amplitude) for either 20 min at a frequency of 3.1 MHz or 10 min at 4.8 MHz. Tissue specimens were maintained at a static overpressure of either 0.52 or 1.2 MPa to suppress bubble activity and tissue vaporization or at atmospheric pressure for control groups. A passive cavitation detector was used to record subharmonic (1.55 or 2.4 MHz), broadband (1.2-1.5 MHz) and low-frequency (5-20 kHz) acoustic emissions. Treated tissue was stained with 2% triphenyl tetrazolium chloride to evaluate thermal lesion dimensions. Subharmonic emissions were significantly reduced in overpressure groups compared with control groups. Correlations observed between acoustic emissions and lesion dimensions were significant and positive for the 3.1-MHz series, but significant and negative for the 4.8-MHz series. The results indicate that for bulk ultrasound ablation, where both acoustic cavitation and tissue vaporization are possible, bubble activity can enhance ablation in the absence of tissue vaporization, but can reduce thermal lesion dimensions in the presence of vaporization.

Use of the Cross-Spectral Density Matrix for Enhanced Passive Ultrasound Imaging of Cavitation

Passive ultrasound imaging is of great interest for cavitation monitoring. Spatiotemporal monitoring of cavitation bubbles in therapeutic applications is possible using an ultrasound imaging probe to passively receive the acoustic signals from the bubbles. Fourier-domain (FD) beamformers have been proposed to process the signals received into maps of the spatial localization of cavitation activity, with reduced computing times with respect to the time-domain approach, and to take advantage of frequency selectivity for cavitation regime characterization. The approaches proposed have been mainly nonadaptive, and these have suffered from low resolution and contrast, due to the many reconstruction artifacts. Inspired by the array-processing literature and in the context of passive ultrasound imaging of cavitation, we propose here a robust estimation of the second-order statistics of data through spatial covariance matrices in the FD or cross-spectral density matrices (CSMs). The benefits of such formalism are illustrated using advanced reconstruction algorithms, such as the robust Capon beamformer, the Pisarenko class beamformer, and the multiple signal classification approach. Through both simulations and experiments in a water tank, we demonstrate that enhanced localization of cavitation activity (i.e., improved resolution and contrast with respect to nonadaptive approaches) is compatible with the rapid and frequency-selective approaches of the FD. Robust estimation of the CSM and the derived adaptive beamformers paves the way to the development of powerful passive ultrasound imaging tools.

Assessment of Collaborative Robot (Cobot)-Assisted Histotripsy for Venous Clot Ablation

OBJECTIVE

The application of bubble-based ablation with the focus ultrasound therapy histotripsy is gaining traction for the treatment of venous thrombosis, among other pathologies. For extensive clot burden, the histotripsy source must be translated to ensure uniform bubble activity throughout the vascular obstruction. The purpose of this study was to evaluate the targeting accuracy of a histotripsy system comprised of a focused source, ultrasound image guidance, and a collaborative robot (cobot) positioner. The system was designed with a primary emphasis for treating deep vein thrombosis.

METHODS

Studies to test treatment planning and targeting bubble activity with the histotripsy-cobot system were conducted in an in vitro clot model. A tissue-mimicking phantom was also targeted with the system, and the predicted and actual areas of liquefaction were compared to gauge the spatial accuracy of ablation.

RESULTS

The system provided submillimeter accuracy for both tracking along an intended path (within 0.6 mm of a model vessel) and targeting bubble activity within the venous clot model (0.7 mm from the center of the clot). Good correlation was observed between the planned and actual liquefaction locations in the tissue phantom, with an average Dice similarity coefficient of 77.8%, and average Hausdorff distance of 1.6 mm.

CONCLUSION

Cobots provide an effective means to apply histotripsy pulses over a treatment volume, with the ablation precision contingent on the quality of image guidance.

SIGNIFICANCE

Overall, these results demonstrate cobots can be used to guide histotripsy ablation for targets that extend beyond the natural focus of the transducer.

Cavitation Emissions Nucleated by Definity Infused through an EkoSonic Catheter in a Flow Phantom

The EkoSonic endovascular system has been cleared by the U.S. Food and Drug Administration for the controlled and selective infusion of physician specified fluids, including thrombolytics, into the peripheral vasculature and the pulmonary arteries. The objective of this study was to explore whether this catheter technology could sustain cavitation nucleated by infused Definity, to support subsequent studies of ultrasound-mediated drug delivery to diseased arteries. The concentration and attenuation spectroscopy of Definity were assayed before and after infusion at 0.3, 2.0 and 4.0 mL/min through the EkoSonic catheter. PCI was used to map and quantify stable and inertial cavitation as a function of Definity concentration in a flow phantom mimicking the porcine femoral artery. The 2.0 mL/min infusion rate yielded the highest surviving Definity concentration and acoustic attenuation. Cavitation was sustained throughout each 15 ms ultrasound pulse, as well as throughout the 3 min infusion. These results demonstrate a potential pathway to use cavitation nucleation to promote drug delivery with the EkoSonic endovascular system.

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.

Evaluation of capacitive micromachined ultrasonic transducers for passive monitoring of microbubble-assisted ultrasound therapies

Passive cavitation detection can be performed to monitor microbubble activity during brain therapy. Microbubbles under ultrasound exposure generate a response characterized by multiple nonlinear emissions. Here, the wide bandwidth of capacitive micromachined ultrasonic transducers (CMUTs) was exploited to monitor the microbubble signature through a rat skull and a macaque skull. The intrinsic nonlinearity of the CMUTs was characterized in receive mode. Indeed, undesirable nonlinear components generated by the CMUTs must be minimized as they can mask the microbubble harmonic response. The microbubble signature at harmonic and ultra-harmonic components (0.5-6 MHz) was successfully extracted through a rat skull using moderate bias voltage.

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.

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.

In vitro assessment of stiffness-dependent histotripsy bubble cloud activity in gel phantoms and blood clots

As a bubble-based ablative therapy, the efficacy of histotripsy has been demonstrated in healthy or acutely diseased models. Chronic conditions associated with stiff tissues may require additional bubble activity prior to histotripsy liquefaction. In this study, histotripsy pulses were generated in agarose phantoms of Young's moduli ranging from 12.3 to 142 kPa, and in vitro clot models with mild and strong platelet-activated retraction. Bubble cloud emissions were tracked with passive cavitation imaging, and the threshold acoustic power associated with phantom liquefaction was extracted with receiver operator characteristic analysis. The power of histotripsy-generated emissions and the degree of liquefaction were tabulated for both clot models. For the agarose phantoms, the acoustic power associated with liquefaction increased with Young's modulus. When grouped based on agarose concentration, only two arms displayed a significant difference in the liquefaction threshold acoustic power (22.1 kPa versus 142 kPa Young's modulus). The bubble cloud dynamics tracked with passive cavitation imaging indicated no strong changes in the bubble dynamics based on the phantom stiffness. For identical histotripsy exposure, the power of acoustic emissions and degree of clot lysis did not vary based on the clot model. Overall, these results indicate that a fixed threshold acoustic power mapped with passive cavitation imaging can be utilized for predicting histotripsy liquefaction over a wide range of tissue stiffness.

Passive cavitation mapping using dual apodization with cross-correlation in ultrasound therapy monitoring

Recently, passive acoustic mapping (PAM) has been successfully applied for dynamic monitoring of ultrasound therapy by beamforming acoustic emissions of cavitation activity during ultrasound exposure. The most widely used PAM algorithm in the literature is time exposure acoustics (TEA), which is a standard delay, sum, and integrate algorithm. However, it results in large point spread function (PSF) and serious imaging artifacts for the case where a narrow-aperture receiving array such as a standard B-mode linear array is used, therefore degrading the quality of cavitation image. To address these challenges, in this paper, we proposed a novel PAM algorithm namely dual apodization with cross-correlation (DAX)-based TEA, in which DAX was originally used as a reconstruction algorithm in medical ultrasound imaging. In the proposed algorithm, two sets of signals were beamformed by two receive apodization functions with alternating elements enabled, and the cross-correlation coefficient of the two signals served as a weighting factor that would be multiplied to the sum of the two signals. The performance of the proposed algorithm was tested on simulated channel data obtained using a multi-bubble model, and experiments were also performed in an in vitro vessel phantom with flowing microbubbles as cavitation nuclei. The reconstructed cavitation images were evaluated quantitatively using established quality metrics including full width at half maximum (FWHM), A_-6dB area, and signal-to-noise ratio (SNR). The results suggested that the proposed algorithm significantly outperformed the conventionally used TEA algorithm. This work may have the potential of providing a useful tool for highly accurate localization of cavitation activity during ultrasound therapy.

Observation and modulation of the dissolution of histotripsy-induced bubble clouds with high-frame rate plane wave imaging

Focused ultrasound therapies are a noninvasive means to ablate tissue. Histotripsy utilizes short ultrasound pulses with sufficient tension to nucleate bubble clouds that impart lethal strain to the surrounding tissues. Tracking bubble cloud dissolution between the application of histotripsy pulses is critical to ensure treatment efficacy. In this study, plane wave B-mode imaging was employed to monitor bubble cloud motion and grayscale at frame rates up to 11.25 kHz. Minimal changes in the area or position of the bubble clouds were observed 50 ms post excitation. The bubble cloud grayscale was observed to decrease with the square root of time, indicating a diffusion-driven process. These results were qualitatively consistent with an analytic model of gas diffusion during the histotripsy process. Finally, the rate of bubble cloud dissolution was found to be dependent on the output of the imaging pulse, indicating an interaction between the bubble cloud and imaging parameters. Overall, these results highlight the utility of plane wave B-mode imaging for monitoring histotripsy bubble clouds.

For Whom the Bubble Grows: Physical Principles of Bubble Nucleation and Dynamics in Histotripsy Ultrasound Therapy

Histotripsy is a focused ultrasound therapy for non-invasive tissue ablation. Unlike thermally ablative forms of therapeutic ultrasound, histotripsy relies on the mechanical action of bubble clouds for tissue destruction. Although acoustic bubble activity is often characterized as chaotic, the short-duration histotripsy pulses produce a unique and consistent type of cavitation for tissue destruction. In this review, the action of histotripsy-induced bubbles is discussed. Sources of bubble nuclei are reviewed, and bubble activity over the course of single and multiple pulses is outlined. Recent innovations in terms of novel acoustic excitations, exogenous nuclei for targeted ablation and histotripsy-enhanced drug delivery and image guidance metrics are discussed. Finally, gaps in knowledge of the histotripsy process are highlighted, along with suggested means to expedite widespread clinical utilization of histotripsy.

Simultaneous Passive Acoustic Mapping and Magnetic Resonance Thermometry for Monitoring of Cavitation-Enhanced Tumor Ablation in Rabbits Using Focused Ultrasound and Phase-Shift Nanoemulsions

Thermal ablation of solid tumors via focused ultrasound (FUS) is a non-invasive image-guided alternative to conventional surgical resection. However, the usefulness of the technique is limited in vascularized organs because of convection of heat, resulting in long sonication times and unpredictable thermal lesion formation. Acoustic cavitation has been found to enhance heating but requires use of exogenous nuclei and sufficient acoustic monitoring. In this study, we employed phase-shift nanoemulsions (PSNEs) to promote cavitation and incorporated passive acoustic mapping (PAM) alongside conventional magnetic resonance imaging (MRI) thermometry within the bore of a clinical MRI scanner. Simultaneous PAM and MRI thermometry were performed in an in vivo rabbit tumor model, with and without PSNE to promote cavitation. Vaporization and cavitation of the nanoemulsion could be detected using PAM, which led to accelerated heating, monitored with MRI thermometry. The maximum heating assessed from MRI was well correlated with the integrated acoustic emissions, illustrating cavitation-enhanced heating. Examination of tissue revealed thermal lesions that were larger in the presence of PSNE, in agreement with the thermometry data. Using fixed exposure conditions over 94 sonications in multiple animals revealed an increase in the mean amplitude of acoustic emissions and resulting temperature rise, but with significant variability between sonications, further illustrating the need for real-time monitoring. The results indicate the utility of combined PAM and MRI for monitoring of tumor ablation and provide further evidence for the ability of PSNEs to promote cavitation-enhanced lesioning.

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.

Pulse inversion enhances the passive mapping of microbubble-based ultrasound therapy

Therapeutic ultrasound combined with preformed circulating microbubbles has enabled non-invasive and targeted drug delivery into the brain, tumors, and blood clots. Monitoring the microbubble activity is essential for the success of such therapies; however, skull and tissues limit our ability to detect low acoustic signals. Here, we show that by emitting consecutive therapeutic pulses of inverse polarity, the sensitivity in the detection of weak bubble acoustic signals during blood-brain barrier opening is enhanced compared to therapeutic pulses of the same polarity. Synchronous passive mapping of the cavitation activity was conducted using delay-and-sum beamforming with absolute time delays, which offers superior spatial resolution compared to the existing asynchronous passive imaging techniques. Sonication with pulse inversion allowed filter-free suppression of the tissue signals by up to 8 dB in a tissue-mimicking phantom and by 7 dB in vivo, compared to exposure without pulse inversion, enabling enhanced passive mapping of microbubble activity. Both therapeutic schemes resulted in similar free-field microbubble activation in vitro and efficient blood-brain barrier opening in vivo.

Ultrasonic Cavitation-Enabled Treatment for Therapy of Hypertrophic Cardiomyopathy: Proof of Principle

Ultrasound myocardial cavitation-enabled treatment was applied to the SS-16^BN rat model of hypertrophic cardiomyopathy for proof of the principle underlying myocardial reduction therapy. A focused ultrasound transducer was targeted using 10-MHz imaging (10 S, GE Vivid 7) to the left ventricular wall of anesthetized rats in a warmed water bath. Pulse bursts of 4-MPa peak rarefactional pressure amplitude were intermittently triggered 1:8 heartbeats during a 10-min infusion of a microbubble suspension. Methylprednisolone was given to reduce initial inflammation, and Losartan was given to reduce fibrosis in the healing tissue. At 28 d post therapy, myocardial cavitation-enabled treatment significantly reduced the targeted wall thickness by 16.2% (p <0.01) relative to shams, with myocardial strain rate and endocardial displacement reduced by 34% and 29%, respectively, which are sufficient for therapeutic treatment. Premature electrocardiogram complexes and plasma troponin measurements were found to identify optimal and suboptimal treatment cohorts and would aid in achieving the desired impact. With clinical translation, myocardial cavitation-enabled treatment should fill the need for a new non-invasive hypertrophic cardiomyopathy therapy option.

Frequency-sum beamforming for passive cavitation imaging

Beamforming includes a variety of spatial filtering techniques that may be used for determining sound source locations from near-field sensor array recordings. For this scenario, beamforming resolution depends on the acoustic frequency, array geometry, and target location. Random scattering in the medium between the source and the array may degrade beamforming resolution with higher frequencies being more susceptible to degradation. The performance of frequency-sum (FS) beamforming for reducing such sensitivity to mild scattering while increasing resolution is reported here. FS beamforming was used with a data-dependent [minimum variance (MV)] or data-independent (delay-and-sum, DAS) weight vector to produce higher frequency information from lower frequency signal components via a quadratic product of complex signal amplitudes. The current findings and comparisons are based on simulations and passive cavitation imaging experiments using 3 MHz and 6 MHz emissions recorded by a 128-element linear array. FS beamforming results are compared to conventional DAS and MV beamforming using four metrics: point spread function (PSF) size, axial and lateral contrast, and computation time. FS beamforming produces a smaller PSF than conventional DAS beamforming with less computation time than MV beamforming in free space and mild scattering environments. However, it may fail when multiple unknown sound sources are present.

Evaluation of a Three Hydrophones Method for 2-Dimensional Cavitation Localization

Cavitation is a critical parameter in various therapeutic applications involving ultrasound (US) such as histotrispy, lithothripsy, drug delivery, and cavitation-enhanced hyperthermia. A cavitation exposure outside the region of interest may lead to suboptimal treatment efficacy or in a worse case, to safety issues. Current methods of localizing cavitation are based on imaging approaches, such as beamforming the cavitation signals received passively by a US imager. These methods, although efficient, require expensive equipment, which may discourage potential future developments. We propose a threehydrophone method to localize the cavitation cloud source. Firstly, the delays between the three receptors are measured by detecting the maximum of their inter-correlations. Then, the position of the source is calculated by either minimizing a cost function or solving hyperbolic equations. After a numerical validation, the method was assessed experimentally. This method was able to track a source displacement with accuracy similar to the size of the cavitation cloud (2-4 millimeters). This light and versatile method provides interesting perspectives since localization can be executed in real time and the extension to three-dimensional localization seems straightforward.

Power cavitation-guided blood-brain barrier opening with focused ultrasound and microbubbles

Image-guided monitoring of microbubble-based focused ultrasound (FUS) therapies relies on the accurate localization of FUS-stimulated microbubble activity (i.e. acoustic cavitation). Passive cavitation imaging with ultrasound arrays can achieve this, but with insufficient spatial resolution. In this study, we address this limitation and perform high-resolution monitoring of acoustic cavitation-mediated blood-brain barrier (BBB) opening with a new technique called power cavitation imaging. By synchronizing the FUS transmit and passive receive acquisition, high-resolution passive cavitation imaging was achieved by using delay and sum beamforming with absolute time delays. Since the axial image resolution is now dependent on the duration of the received acoustic cavitation emission, short pulses of FUS were used to limit its duration. Image sets were acquired at high-frame rates for calculation of power cavitation images analogous to power Doppler imaging. Power cavitation imaging displays the mean intensity of acoustic cavitation over time and was correlated with areas of acoustic cavitation-induced BBB opening. Power cavitation-guided BBB opening with FUS could constitute a standalone system that may not require MRI guidance during the procedure. The same technique can be used for other acoustic cavitation-based FUS therapies, for both safety and guidance.

Diffraction Effects and Compensation in Passive Acoustic Mapping

Over the last decade, a variety of noninvasive techniques have been developed to monitor therapeutic ultrasound procedures in support of safety or efficacy assessments. One class of methods employs diagnostic ultrasound arrays to sense acoustic emissions, thereby providing a means to passively detect, localize, and quantify the strength of nonlinear sources, including cavitation. Real array element diffraction patterns may differ substantially from those presumed in existing beamforming algorithms. However, diffraction compensation has received limited treatment in passive and active imaging, and measured diffraction data have yet to be used for array response correction. The objectives of this paper were to identify differences between ideal and real element diffraction patterns, and to quantify the impact of diffraction correction on cavitation mapping beamformer performance. These objectives were addressed by performing calibration measurements on a diagnostic linear array, using the results to calculate diffraction correction terms, and applying the corrections to cavitation emission data collected from soft tissue phantom experiments. Measured diffraction patterns were found to differ significantly from those of ideal element forms, particularly at higher frequencies and shorter distances from the array. Diffraction compensation of array data resulted in cavitation energy estimates elevated by as much as a factor of 5, accompanied by the elimination of a substantial bias between two established beamforming algorithms. These results illustrate the importance of using measured array responses to validate analytical field models and to minimize observation biases in imaging applications where quantitative analyses are critical for assessment of therapeutic safety and efficacy.

Post Hoc Analysis of Passive Cavitation Imaging for Classification of Histotripsy-Induced Liquefaction in Vitro

Histotripsy utilizes focused ultrasound to generate bubble clouds for transcutaneous tissue liquefaction. Bubble activity maps are under development to provide image guidance and monitor treatment progress. The aim of this paper was to investigate the feasibility of using plane wave B-mode and passive cavitation images to be used as binary classifiers of histotripsy-induced liquefaction. Prostate tissue phantoms were exposed to histotripsy pulses over a range of pulse durations (5- ) and peak negative pressures (12-23 MPa). Acoustic emissions were recorded during the insonation and beamformed to form passive cavitation images. Plane wave B-mode images were acquired following the insonation to detect the hyperechoic bubble cloud. Phantom samples were sectioned and stained to delineate the liquefaction zone. Correlation between passive cavitation and plane wave B-mode images and the liquefaction zone was assessed using receiver operating characteristic (ROC) curve analysis. Liquefaction of the phantom was observed for all the insonation conditions. The area under the ROC (0.94 versus 0.82), accuracy (0.90 versus 0.83), and sensitivity (0.81 versus 0.49) was greater for passive cavitation images relative to B-mode images ( ) along the azimuth of the liquefaction zone. The specificity was greater than 0.9 for both imaging modalities. These results demonstrate a stronger correlation between histotripsy-induced liquefaction and passive cavitation imaging compared with the plane wave B-mode imaging, albeit with limited passive cavitation image range resolution.

Real-time monitoring of controllable cavitation erosion in a vessel phantom with passive acoustic mapping

Cavitation erosion in blood vessel plays an important role in ultrasound thrombolysis, drug delivery, and other clinical applications. The controllable superficial vessel erosion based on ultrasonic standing wave (USW) has been used to effectively prevent vessel ruptures and haemorrhages, and optical method is used to observe the experiments. But optical method can only work in transparent media. Compared with standard B-mode imaging, passive acoustic mapping (PAM) can monitor erosion in real time and has better sensitivity of cavitation detection. However, the conventionally used PAM has limitations in imaging resolution and artifacts. In this study, a unique PAM method that combined the robust Capon beamformer (RCB) with the sign coherence factor (SCF) was proposed to monitor the superficial vessel erosion in real time. The performance of the proposed method was validated by simulations. In vitro experiments showed that the lateral (axial) resolution of the proposed PAM was 2.31±0.51 (3.19±0.38) times higher than time exposure acoustics (TEA)-based PAM and 1.73±0.38 (1.76±0.48) times higher than RCB-based PAM, and the cavitation-to-artifact ratio (CAR) of the proposed PAM could be improved by 22.5±3.2dB and 7.1±1.2dB compared with TEA and RCB-based PAM. These results showed that the proposed PAM can precisely monitor the superficial vessel erosion and the erosion shift after USW modulation. This work may have the potential of developing a useful tool for precise spatial control and real-time monitoring of the superficial vessel erosion.

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

It is well established that high intensity focused ultrasound can be used to disintegrate clots. This approach has the potential to rapidly and noninvasively resolve clot causing occlusions in cardiovascular diseases such as deep vein thrombosis (DVT). However, lack of an appropriate treatment monitoring tool is currently a limiting factor in its widespread adoption. Here we conduct cavitation imaging with a large aperture, sparse hemispherical receiver array during sonothrombolysis with multi-cycle burst exposures (0.1 or 1 ms burst lengths) at 1.51 MHz. It was found that bubble cloud generation on imaging correlated with the locations of clot degradation, as identified with high frequency (30 MHz) ultrasound following exposures. 3D images could be formed at integration times as short as 1 µs, revealing the initiation and rapid development of cavitation clouds. Equating to megahertz frame rates, this is an order of magnitude faster than any other imaging technique available for in vivo application. Collectively, these results suggest that the development of a device to perform DVT therapy procedures would benefit greatly from the integration of receivers tailored to bubble activity imaging.

Multi-focal HIFU reduces cavitation in mild-hyperthermia

Background

Mild-hyperthermia therapy (40–45 °C) with high-intensity focused ultrasound (HIFU) is a technique being considered in a number of different treatments such as thermally activated drug delivery, immune-stimulation, and as a chemotherapy adjuvant. Mechanical damage and loss of cell viability associated with HIFU-induced acoustic cavitation may pose a risk during these treatments or may hinder their success. Here we present a method that achieves mild heating and reduces cavitation by using a multi-focused HIFU beam. We quantify cavitation level and temperature rise in multi-focal sonications and compare it to single-focus sonications at the transducer geometric focus.

Methods

Continuous wave sonications were performed with the Sonalleve V2 transducer in gel phantoms and pork at 5, 10, 20, 40, 60, 80 acoustic watts for 30 s. Cavitation activity was measured with two ultrasound (US) imaging probes, both by computing the raw channel variance and using passive acoustic mapping (PAM). Temperature rise was measured with MR thermometry at 3 T. Cavitation and heating were compared for single- and multi-focal sonication geometries. Multi-focal sonications used four points equally spaced on a ring of either 4 mm or 8 mm diameter. Single-focus sonications were not steered.

Results

Multi-focal sonication generated distinct foci that were visible in MRI thermal maps in both phantoms and pork, and visible in PAM images in phantoms only. Cavitation activity (measured by channel variance) and mean PAM image value were highly correlated (r > 0.9). In phantoms, cavitation exponentially decreased over the 30-second sonication, consistent with depletion of cavitation nuclei. In pork, sporadic spikes signaling cavitation were observed with single focusing only. In both materials, the widest beam reduced average and peak cavitation level by a factor of two or more at each power tested when compared to a single focus. The widest beam reduced peak temperature by at least 10 °C at powers above 5 W, and created heating that was more spatially diffuse than single focus, resulting in more voxels in the mild heating (3–8 °C) range.

Conclusions

Multi-focal HIFU can be used to achieve mild temperature elevation and reduce cavitation activity.

Quantitative Frequency-Domain Passive Cavitation Imaging

Passive cavitation detection has been an instrumental technique for measuring cavitation dynamics, elucidating concomitant bioeffects, and guiding ultrasound therapies. Recently, techniques have been developed to create images of cavitation activity to provide investigators with a more complete set of information. These techniques use arrays to record and subsequently beamform received cavitation emissions, rather than processing emissions received on a single-element transducer. In this paper, the methods for performing frequency-domain delay, sum, and integrate passive imaging are outlined. The method can be applied to any passively acquired acoustic scattering or emissions, including cavitation emissions. To compare data across different systems, techniques for normalizing Fourier transformed data and converting the data to the acoustic energy received by the array are described. A discussion of hardware requirements and alternative imaging approaches is additionally outlined. Examples are provided in MATLAB.

High-contrast active cavitation imaging technique based on multiple bubble wavelet transform

In this study, a unique method that combines the ultrafast active cavitation imaging technique with multiple bubble wavelet transform (MBWT) for improving cavitation detection contrast was presented. The bubble wavelet was constructed by the modified Keller-Miksis equation that considered the mutual effect among bubbles. A three-dimensional spatial model was applied to simulate the spatial distribution of multiple bubbles. The effects of four parameters on the signal-to-noise ratio (SNR) of cavitation images were evaluated, including the following: initial radii of bubbles, scale factor in the wavelet transform, number of bubbles, and the minimum inter-bubble distance. And the other two spatial models and cavitation bubble size distributions were introduced in the MBWT method. The results suggested that in the free-field experiments, the averaged SNR of images acquired by the MBWT method was improved by 7.16 ± 0.09 dB and 3.14 ± 0.14 dB compared with the values of images acquired by the B-mode and single bubble wavelet transform (SBWT) methods. In addition, in the tissue experiments, the averaged cavitation-to-tissue ratio of cavitation images acquired by the MBWT method was improved by 4.69 ± 0.25 dB and 1.74± 0.29 dB compared with that of images acquired by B-mode and SBWT methods.