Broadband Ultrasonic Attenuation Estimation and Compensation With Passive Acoustic Mapping

Histotripsy and Catheter-Directed Lytic: Efficacy in Highly Retracted Porcine Clots In Vitro

OBJECTIVE

Standard treatment for deep vein thrombosis (DVT) involves catheter-directed anticoagulants or thrombolytics, but the chronic thrombi present in many DVT cases are often resistant to this therapy. Histotripsy has been found to be a promising adjuvant treatment, using the mechanical action of cavitating bubble clouds to enhance thrombolytic activity. The objective of this study was to determine if histotripsy enhanced recombinant tissue plasminogen activator (rt-PA) thrombolysis in highly retracted porcine clots in vitro in a flow model of occlusive DVT.

METHODS

Highly retracted porcine whole blood clots were treated for 1 h with either catheter-directed saline (negative control), rt-PA (lytic control), histotripsy, DEFINITY and histotripsy or the combination of rt-PA and histotripsy with or without DEFINITY. Five-cycle, 1.5 MHz histotripsy pulses with a peak negative pressure of 33.2 MPa and pulse repetition frequency of 40 Hz were applied along the clot. B-Mode and passive cavitation images were acquired during histotripsy insonation to monitor bubble activity.

RESULTS

Clots subjected to histotripsy with and without rt-PA exhibited greater thrombolytic efficacy than controls (7.0% flow recovery or lower), and histotripsy with rt-PA was more efficacious than histotripsy with saline (86.1 ± 10.2% compared with 61.7 ± 19.8% flow recovery). The addition of DEFINITY to histotripsy with or without rt-PA did not enhance either thrombolytic efficacy or cavitation dose. Cavitation dose generally did not correlate with thrombolytic efficacy.

CONCLUSION

Enhancement of thrombolytic efficacy was achieved using histotripsy, with and without catheter-directed rt-PA, in the presence of physiologic flow. This suggests these treatments may be effective as therapy for DVT.

Passive Acoustic Mapping for Convex Arrays With the Helical Wave Spectrum Method

Passive acoustic mapping (PAM) has emerged as a valuable imaging modality for monitoring the cavitation activity in focused ultrasound therapies. When it comes to imaging in the human abdomen, convex arrays are preferred due to their large acoustic window. However, existing PAM methods for convex arrays rely on the computationally expensive delay-and-sum (DAS) operation limiting the image reconstruction speed when the field-of-view (FOV) is large. In this work, we propose an efficient and frequency-selective PAM method for convex arrays. This method is based on projecting the helical wave spectrum (HWS) between cylindrical surfaces in the imaging field. Both the in silico and in vitro experiments showed that the HWS method has comparable image quality and similar acoustic cavitation source localization accuracy as the DAS-based methods. Compared to the frequency-domain and time-domain DAS methods, the time-complexity of the HWS method is reduced by one order and two orders of magnitude, respectively. A parallel implementation of the HWS method realized millisecond-level image reconstruction speed. We also show that the HWS method is inherently capable of mapping microbubble (MB) cavitation activity of different status, i.e., no cavitation, stable cavitation, or inertial cavitation. After compensating for the lens effects of the convex array, we further combined PAM formed by the HWS method and B-mode imaging as a real-time dual-mode imaging approach to map the anatomical location where MBs cavitate in a liver phantom experiment. This method may find use in applications where convex arrays are required for cavitation activity monitoring in real time.

Quantifying the Effect of Acoustic Parameters on Temporal and Spatial Cavitation Activity: Gauging Cavitation Dose

OBJECTIVE

Cavitation-enhanced delivery of therapeutic agents is under development for the treatment of cancer and neurodegenerative and cardiovascular diseases, including sonothrombolysis for deep vein thrombosis. The objective of this study was to quantify the spatial and temporal distribution of cavitation activity nucleated by Definity infused through the EKOS catheter over a range of acoustic parameters controlled by the EKOS endovascular system.

METHODS

Three insonation protocols were compared in an in vitro phantom mimicking venous flow to measure the effect of peak rarefactional pressure, pulse duration and pulse repetition frequency on cavitation activity energy, location and duration. Inertial and stable cavitation activity was quantified using passive cavitation imaging, and a metric of cavitation dose based on energy density was defined.

RESULTS

For all three insonation protocols, cavitation was sustained for the entire 30 min Definity infusion. The evolution of cavitation energy during each pulse duration was similar for all three protocols. For insonation protocols with higher peak rarefactional acoustic pressures, inertial and stable cavitation doses also increased. A complex relationship between the temporal behavior of cavitation energy within each pulse and the pulse repetition frequency affected the cavitation dose for the three insonation protocols. The relative predominance of stable or inertial cavitation dose varied according to insonation schemes. Passive cavitation images revealed the spatial distribution of cavitation activity.

CONCLUSION

Our cavitation dose metric based on energy density enabled the impact of different acoustic parameters on cavitation activity to be measured. Depending on the type of cavitation to be promoted or suppressed, particular pulsing schemes could be employed in future studies, for example, to correlate cavitation dose with sonothrombolytic efficacy.

Passive Cavitation Mapping by Cavitation Source Localization From Aperture-Domain Signals-Part I: Theory and Validation Through Simulations

Passive cavitation mapping (PCM) algorithms for diagnostic ultrasound arrays based on time exposure acoustics (TEA) exhibit poor axial resolution, which is in part due to the diffraction-limited point spread function of the imaging system and poor rejection by the delay-and-sum beamformer. In this article, we adapt a method for speed of sound estimation to be utilized as a cavitation source localization (CSL) approach. This method utilizes a hyperbolic fit to the arrival times of the cavitation signals in the aperture domain, and the coefficients of the fit are related to the position of the cavitation source. Wavefronts exhibiting poor fit to the hyperbolic function are corrected to yield improved source localization. We demonstrate through simulations that this method is capable of accurate estimation of the origin of coherent spherical waves radiating from cavitation/point sources. The average localization error from simulated microbubble sources was 0.12 ± 0.12mm ( 0.15 ± 0.14λ₀ for a 1.78-MHz transmit frequency). In simulations of two simultaneous cavitation sources, the proposed technique had an average localization error of 0.2mm ( 0.23λ₀ ), whereas conventional TEA had an average localization error of 0.81mm ( 0.97λ₀ ). The reconstructed PCM-CSL image showed a significant improvement in resolution compared with the PCM-TEA approach.

Dual-Array Passive Acoustic Mapping for Cavitation Imaging With Enhanced 2-D Resolution

Passive acoustic mapping (PAM) techniques have been developed for the purposes of detecting, localizing, and quantifying cavitation activity during therapeutic ultrasound procedures. Implementation with conventional diagnostic ultrasound arrays has allowed planar mapping of bubble acoustic emissions to be overlaid with B-mode anatomical images, with a variety of beamforming approaches providing enhanced resolution at the cost of extended computation times. However, no passive signal processing techniques implemented to date have overcome the fundamental physical limitation of the conventional diagnostic array aperture that results in point spread functions with axial/lateral beamwidth ratios of nearly an order of magnitude. To mitigate this problem, the use of a pair of orthogonally oriented diagnostic arrays was recently proposed, with potential benefits arising from the substantially expanded range of observation angles. This article presents experiments and simulations intended to demonstrate the performance and limitations of the dual-array system concept. The key finding of this study is that source pair resolution of better than 1 mm is now possible in both dimensions of the imaging plane using a pair of 7.5-MHz center frequency conventional arrays at a distance of 7.6cm. With an eye toward accelerating computations for real-time applications, channel count reductions of up to a factor of eight induce negligible performance losses. Modest sensitivities to sound speed and relative array position uncertainties were identified, but if these can be kept on the order of 1% and 1 mm, respectively, then the proposed methods offer the potential for a step improvement in cavitation monitoring capability.

Spatiotemporal Assessment of the Cellular Safety of Cavitation-Based Therapies by Passive Acoustic Mapping

Many useful therapeutic bio-effects can be generated using ultrasound-induced cavitation. However, cavitation is also capable of causing unwanted cellular and vascular damage, which should be monitored to ensure treatment safety. In this work, the unique opportunity provided by passive acoustic mapping (PAM) to quantify cavitation dose across an entire volume of interest during therapy is utilised to provide setup-independent measures of spatially localised cavitation dose. This spatiotemporally quantifiable cavitation dose is then related to the level of cellular damage generated. The cavitation-mediated destruction of equine red blood cells mixed with one of two types of cavitation nuclei at a variety of concentrations is investigated. The blood is placed within a 0.5-MHz ultrasound field and exposed to a range of peak rarefactional pressures up to 2 MPa, with 50 to 50,000 cycle pulses maintaining a 5% duty cycle. Two co-planar linear arrays at 90° to each other are used to generate 400-µm-resolution frequency domain robust capon beamforming PAM maps, which are then used to generate estimates of cavitation dose. A relationship between this cavitation dose and the levels of haemolysis generated was found which was comparable regardless of the applied acoustic pressure, pulse length, cavitation agent type or concentration used. PAM was then used to monitor cellular damage in multiple locations within a tissue phantom simultaneously, with the damage-cavitation dose relationship being similar for the two experimental models tested. These results lay the groundwork for this method to be applied to other measures of safety, allowing for improved ultrasound monitoring of cavitation-based therapies.

Assessment of corneal viscoelasticity using elastic wave optical coherence elastography

The corneal viscoelasticity have great clinical significance, such as the early diagnosis of keratoconus. In this work, an analysis method which utilized the elastic wave velocity, frequency and energy attenuation to assess the corneal viscoelasticity is presented. Using phase-resolved optical coherence tomography, the spatial-temporal displacement map is derived. The phase velocity dispersion curve and center frequency are obtained by transforming the displacement map into the wavenumber-frequency domain through the 2D fast Fourier transform (FFT). The shear modulus is calculated through Rayleigh wave equation using the phase velocity in the high frequency. The normalized energy distribution is plotted by transforming the displacement map into the spatial-frequency domain through the 1D FFT. The energy attenuation coefficient is derived by exponential fitting to calculate the viscous modulus. Different concentrations of tissue-mimicking phantoms and porcine corneas are imaged to validate this method, which demonstrates that the method has the capability to assess the corneal viscoelasticity.

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.

Compensation of array lens effects for improved co-registration of passive acoustic mapping and B-mode images for cavitation monitoring

Passive acoustic mapping (PAM) techniques offer a simple means of spatio-temporal cavitation monitoring during therapeutic ultrasound procedures. Implementation with a conventional diagnostic ultrasound system allows natural integration of PAM with B-mode imaging. However, the refracting properties of diagnostic array lenses may introduce PAM image registration errors that could lead to inaccuracies in treatment monitoring and guidance. To address these concerns, this paper presents lens characterization of two different array designs, analytical estimation of lens-induced source mapping errors in simple media, and experimental demonstration and correction of lens effects, reducing the depth-averaged image co-registration errors to no more than 0.52 mm.

Compensation of array lens effects for improved co-registration of passive acoustic mapping and B-mode images for cavitation monitoring

Passive acoustic mapping (PAM) techniques offer a simple means of spatio-temporal cavitation monitoring during therapeutic ultrasound procedures. Implementation with a conventional diagnostic ultrasound system allows natural integration of PAM with B-mode imaging. However, the refracting properties of diagnostic array lenses may introduce PAM image registration errors that could lead to inaccuracies in treatment monitoring and guidance. To address these concerns, this paper presents lens characterization of two different array designs, analytical estimation of lens-induced source mapping errors in simple media, and experimental demonstration and correction of lens effects, reducing the depth-averaged image co-registration errors to no more than 0.52 mm.