Receiver array design for sonothrombolysis treatment monitoring in deep vein thrombosis

A numerical investigation of passive acoustic mapping for monitoring bubble-mediated focused ultrasound treatment of the spinal cord

Focused ultrasound (FUS) combined with intravenous microbubbles (MBs) has been shown to increase drug delivery to the spinal cord in animal models. Eventual clinical translation of such a technique in the sensitive spinal cord requires robust treatment monitoring to ensure efficacy, localization, safety, and provide key intraprocedural feedback. Here, the use of passive acoustic mapping (PAM) of MB emissions with a spine-specific detector array in the context of transvertebral FUS sonications is investigated in silico. Using computed tomography-derived human vertebral geometry, transvertebral detection of MBs is evaluated over varying source locations with and without phase and amplitude corrections (PACs). The impact of prefocal cavitation is studied by simulating concurrent cavitation events in the canal and pre-laminar region. Spatially sensitive application of phase and amplitude is used to balance signal strengths emanating from different axial depths in combination with multiple dynamic ranges to elicit multisource viewing. Collectively, the results of this study encourage the use of PAM in transvertebral FUS applications with PACs to not only localize sources originating in the spinal canal but also multiple sources of innate amplitude mismatches when corrective methods are applied.

High-Pressure Low-Frequency Lateral Mode Phased-Array Transducer System for the Treatment of Deep Vein Thrombosis: An In Vitro Study

Deep vein thrombosis (DVT) can lead to a fatal disease known as pulmonary embolism. Application of high-power ultrasound has been successful in studies to mechanically fragment the clots. Single-element ultrasound transducers were used in most of the studies. Challenges associated with phased arrays, such as high electrical impedance and element breakdown at high voltages, were addressed in the previous study, and a high-power 64-element transducer module was designed and fabricated. In this study, a cylindrical array of 16 modules with the frequency of 260 kHz was modeled and constructed for DVT thrombolysis. The maximum pressure, focal size, and steering ability of the array were examined. In vitro experiments were conducted to assess the performance of the array. The simulated pressure amplitude of 34 MPa at the depth of 55 mm (average femoral vein (FV) distance from the inner surface of the thigh) was in consistent with the experiments and satisfied the purpose of this study. Moreover, the employed module distribution resulted in a focal spot dimension of 2.4×2.8×7.3 mm³ (at the 75% pressure amplitude level) that can be confined in a human FV with the average diameter of 12 mm. In vitro experiments manifested a partial and complete clot breakdown at 11.5- and 15-MPa pressure at the focus. The design and engineering of the array system was succeeded in maintaining the desired pressure and focal size even when steered. The results presented in this study suggest the potential of the designed array system for clinical applications.

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.

High-Power Phased-Array Transducer Module for the Construction of a System for the Treatment of Deep Vein Thrombosis

Blood clot can be disintegrated by high-intensity focused ultrasound alone through inertial cavitation. There are limitations in using single-element ultrasound transducers for this purpose such as lack of steerability and control of the focus in terms of shape and location. Phased-array transducers being able to rapidly scan over the clots can alleviate this problem. A full 3-D control of the ultrasound beam can be achieved by 2-D electronically steerable arrays. However, the required high-pressure amplitude has not been possible with such arrays. In this work, a 2-D 64-element fully populated phased-array transducer module was designed and fabricated for the high-pressure amplitude required for deep vein thrombosis (DVT). Lateral coupling was considered for the transducer design to decrease the electrical impedance and eliminate the need for electrical matching circuit. PZT-4 with a thickness of 0.35 mm, an element surface area of [Formula: see text] mm, and a length of 6 mm showed a mean electrical impedance of 60.4 ± 11.5 measured for each transducer element facilitating effective electric power transfer from the driving electronics. No breakdown was observed when the voltage was increased gradually to 180 ± 3 V_pp. Operation at 180 V_pp was found to be safe over 10,000 repetitions without reduction in the power, resulting in the average pressure amplitude of 1.01 ± 0.09 MPa at 2 mm from the element surface. These pressure amplitude values indicate that an array of eight modules (80 [Formula: see text] mm) is required to reach to the pressure amplitude needed for DVT. Such arrays are practical with the current technology.

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.

Enhanced Detection of Bubble Emissions Through the Intact Spine for Monitoring Ultrasound-Mediated Blood-Spinal Cord Barrier Opening

OBJECTIVE

We previously developed short burst, phase keying (SBPK) focused ultrasound (FUS) to mitigate standing waves in the human vertebral canal. Here, we show microbubble emissions from these pulses can be detected through the human vertebral arch and that these pulses are effective for blood-spinal cord barrier (BSCB) opening.

METHODS

At f₀ = 514 kHz, circulating microbubbles were sonicated through ex vivo human vertebrae (60 kPa-1 MPa) using a dual-aperture approach and SBPK exposures engineered to incorporate pulse inversion (PI). Signals from a 250 kHz receiver were analyzed using PI, short-time Fourier analysis and the maximum projection over the pulse train. In rats (n = 14), SBPK FUS+microbubbles was applied to 3 locations/spinal cord at fixed pressures (∼0.20-0.47 MPa). MRI and histology were used to assess opening and tissue damage.

RESULTS

In human vertebrae between 0.2-0.4 MPa, PI amplified the microbubble/baseline ratio at f₀/2 and 2f₀ by 202 ± 40% (132-291%). This was maximal at 0.4 MPa, coinciding with the onset of broadband emissions. In vivo, opening was achieved at 40/42 locations, with mean MRI enhancement of 46 ± 32%(16%-178%). Using PI, f₀/2 was detected at 14/40 opening locations. At the highest pressures (f₀/2 present) histology showed widespread bleeding throughout the focal region. At the lowest pressures, opening was achieved without bleeding.

CONCLUSION

This study confirmed that PI can increase sensitivity to transvertebral detection of microbubble signals. Preliminary in vivo investigations show that SBPK FUS can increase BSCB permeability without tissue damage.

SIGNIFICANCE

SBPK is a clinically relevant pulse scheme and, in combination with PI, provides a means of mediating and monitoring BSCB opening noninvasively.

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.