Spatiotemporal Clutter Filtering of Ultrafast Ultrasound Data Highly Increases Doppler and fUltrasound Sensitivity

A transcranial multiple waves suppression method for plane wave imaging based on Radon transform

Transcranial ultrasound imaging presents a significant challenge due to the intricate interplay between ultrasound waves and the heterogeneous human skull. The skull's presence induces distortion, refraction, multiple scattering, and reflection of ultrasound signals, thereby complicating the acquisition of high-quality images. Extracting reflections from the entire waveform is crucial yet exceedingly challenging, as intracranial reflections are often obscured by strong amplitude direct waves and multiple scattering. In this paper, a multiple wave suppression method for ultrasound plane wave imaging is proposed to mitigate the impact of skull interference. Drawing upon prior research, we developed an enhanced high-resolution linear Radon transform using the maximum entropy principle and Bayesian method, facilitating wavefield separation. We detailed the process of wave field separation in the Radon domain through simulation of a model with a high velocity layer. When plane waves emitted at any steering angles, both multiple waves and first arrival waves manifested as distinct energy points. In the brain simulation, we contrasted the characteristic differences between skull reflection and brain-internal signal in Radon domain, and demonstrated that multiples suppression method reduces side and grating lobe levels by approximately 30 dB. Finally, we executed in vitro experiments using a monkey skull to separate weak intracranial reflection signals from strong skull reflections, enhancing the contrast-to-noise ratio by 85 % compared to conventional method using full waveform. This study deeply explores the effect of multiples on effective signal separation, addresses the complexity of wavefield separation, and verifies its efficacy through imaging, thereby significantly advancing ultrasound transcranial imaging techniques.

Dynamic Ultrasound Localization Microscopy Without ECG-Gating

OBJECTIVE

Dynamic Ultrasound Localization Microscopy (DULM) has first been developed for non-invasive Pulsatility measurements in the rodent brain. DULM relies on the localization and tracking of microbubbles (MBs) injected into the bloodstream, to obtain highly resolved velocity and density cine-loops. Previous DULM techniques required ECG-gating, limiting its application to specific datasets, and increasing acquisition time. The objective of this study is to eliminate the need for ECG-gating in DULM experiments by introducing a motion-matching method for time registration.

METHODS

We developed a motion-matching algorithm based on tissue Doppler that leverages the cyclic tissue motion within the brain. Tissue Doppler was estimated for each group of frames in the acquisitions, at multiple locations identified as local maxima in the skin above the skull. Subsequently, each group of frames was time-registered to a reference group by delaying it based on the maximum correlation value between their respective tissue Doppler signals. This synchronization ensured that each group of frames aligned with the brain tissue motion of the reference group, and consequently, with its cardiac cycle. As a result, velocities of MBs could be averaged to retrieve flow velocity variations over time.

RESULTS

Initially validated in ECG-gated acquisitions in a rat model (n = 1), the proposed method was successfully applied in a mice model in 2D (n = 3) and in a feline model in 3D (n = 1). Performing time-registration with the proposed motion-matching method or by using ECG-gating leads to similar results. For the first time, dynamic velocity and density cine-loops were extracted without the need for any information on the animal ECG, and complex dynamic markers such as the Pulsatility index were estimated.

CONCLUSION

Results suggest that DULM can be performed without external gating, enabling the use of DULM on any ULM dataset where enough MBs are detectable. Time registration by motion-matching represents a significant advancement in DULM techniques, making DULM more accessible by simplifying its experimental complexity.

Singular value decomposition with weighting matrix applied for optical-resolution photoacoustic microscopes

The prestige target selectivity and imaging depth of optical-resolution photoacoustic microscope (OR-PAM) have gained attentions to enable advanced intra-cellular visualizations. However, the broad-band nature of photoacoustic signals is prone to noise and artifacts caused by the inefficient light-to-pressure translation, resulting in poor image quality. The present study foresees application of singular value decomposition (SVD) to effectively extract the photoacoustic signals from these noise and artifacts. Although spatiotemporal SVD succeeded in ultrasound flow signal extraction, the conventional multi frame model is not suitable for data acquired with scanning OR-PAM due to the burden of accessing multiple frames. To utilize SVD on the OR-PAM, this study began with exploring SVD applied on multiple A-lines of photoacoustic signal instead of frames. Upon explorations, an obstacle of uncertain presence of unwanted singular vectors was observed. To tackle this, a data-driven weighting matrix was designed to extract relevant singular vectors based on the analyses of temporal-spatial singular vectors. Evaluation on the extraction capability by the SVD with the weighting matrix showed a superior signal quality with efficient computation against past studies. In summary, this study contributes to the field by providing exploration of SVD applied on A-line signals as well as its practical utilization to distinguish and recover photoacoustic signals from noise and artifact components.

What is the place of ultrasound in MSK imaging?

During the past four decades, ultrasound has become popular as an imaging modality applied to the musculoskeletal (MSK) system, particularly outside the USA, due to its low cost, accessibility, and lack of ionizing radiation. A basic requirement in performing these examinations is to have a core group of radiologists and ultrasound technologists with expertise in MSK ultrasound. The extent to which ultrasound will be part of the imaging offered by a particular radiology practice or in an academic institution will vary according to expertise, availability, and reimbursements. A brief discussion of the technical capabilities of the current generation of ultrasound scanners will be followed by a description of some of the more prevalent MSK ultrasound imaging applications. The extent to which training to perform these exams within and outside of Radiology plays a role is discussed. Applications that are unique to ultrasound, such as dynamic evaluation of musculoskeletal anatomy and some, US-guided interventions are an important part of MSK imaging. Ultrasound is increasingly important in the assessment of superficial structures, such as tendons, small joints, and peripheral nerves. These applications help to establish the place of ultrasound as an important part of the Radiologists approach to MSK imaging. Outside of radiology, for a variety of clinical subspecialties, ultrasound already plays an integral role in MSK imaging.

Transcranial ultrafast ultrasound Doppler imaging: A phantom study

Ultrafast ultrasound Doppler imaging facilitates the assessment of cerebral hemodynamics with high spatio-temporal resolution. However, the significant acoustic impedance mismatch between the skull and soft tissue results in phase aberrations, which can compromise the quality of transcranial imaging and introduce biases in velocity and direction quantification of blood flow. This paper proposed an aberration correction method that combines deep learning-based skull sound speed modelling with ray theory to realize transcranial plane-wave imaging and ultrafast Doppler imaging. The method was validated through phantom experiments using a linear array with a center frequency of 6.25 MHz, 128 elements, and a pitch of 0.3 mm. The results demonstrated an improvement in the imaging quality of intracranial targets when using the proposed method. After aberration correction, the average locating deviation decreased from 1.40 mm to 0.27 mm in the axial direction, from 0.50 mm to 0.20 mm in the lateral direction, and the average full-width-at-half-maximum (FWHM) decreased from 1.37 mm to 0.97 mm for point scatterers. For circular inclusions, the average contrast-to-noise ratio (CNR) improved from 8.1 dB to 11.0 dB, and the average eccentricity decreased from 0.36 to 0.26. Furthermore, the proposed method was applied to transcranial ultrafast Doppler flow imaging. The results showed a significant improvement in accuracy and quality for blood Doppler flow imaging. The results in the absence of the skull were considered as the reference, and the average normalized root-mean-square errors of the axial velocity component on the five selected axial profiles were reduced from 17.67% to 8.02% after the correction.

A Pseudo-Spectral Method for Wall Shear Stress Estimation from Doppler Ultrasound Imaging in Coronary Arteries

PURPOSE

The Wall Shear Stress (WSS) is the component tangential to the boundary of the normal stress tensor in an incompressible fluid, and it has been recognized as a quantity of primary importance in predicting possible adverse events in cardiovascular diseases, in general, and in coronary diseases, in particular. The quantification of the WSS in patient-specific settings can be achieved by performing a Computational Fluid Dynamics (CFD) analysis based on patient geometry, or it can be retrieved by a numerical approximation based on blood flow velocity data, e.g., ultrasound (US) Doppler measurements. This paper presents a novel method for WSS quantification from 2D vector Doppler measurements.

METHODS

Images were obtained through unfocused plane waves and transverse oscillation to acquire both in-plane velocity components. These velocity components were processed using pseudo-spectral differentiation techniques based on Fourier approximations of the derivatives to compute the WSS.

RESULTS

Our Pseudo-Spectral Method (PSM) is tested in two vessel phantoms, straight and stenotic, where a steady flow of 15 mL/min is applied. The method is successfully validated against CFD simulations and compared against current techniques based on the assumption of a parabolic velocity profile. The PSM accurately detected Wall Shear Stress (WSS) variations in geometries differing from straight cylinders, and is less sensitive to measurement noise. In particular, when using synthetic data (noise free, e.g., generated by CFD) on cylindrical geometries, the Poiseuille-based methods and PSM have comparable accuracy; on the contrary, when using the data retrieved from US measures, the average error of the WSS obtained with the PSM turned out to be 3 to 9 times smaller than that obtained by state-of-the-art methods.

CONCLUSION

The pseudo-spectral approach allows controlling the approximation errors in the presence of noisy data. This gives a more accurate alternative to the present standard and a less computationally expensive choice compared to CFD, which also requires high-quality data to reconstruct the vessel geometry.

Displacement and functional ultrasound (fUS) imaging of displacement-guided focused ultrasound (FUS) neuromodulation in mice

Focused ultrasound (FUS) stimulation is a promising neuromodulation technique with the merits of non-invasiveness, high spatial resolution, and deep penetration depth. However, simultaneous imaging of FUS-induced brain tissue displacement and the subsequent effect of FUS stimulation on brain hemodynamics has proven challenging thus far. In addition, earlier studies lack in situ confirmation of targeting except for the magnetic resonance imaging-guided FUS system-based studies. The purpose of this study is 1) to introduce a fully ultrasonic approach to in situ target, modulate neuronal activity, and monitor the resultant neuromodulation effect by respectively leveraging displacement imaging, FUS, and functional ultrasound (fUS) imaging, and 2) to investigate FUS-evoked cerebral blood volume (CBV) response and the relationship between CBV and displacement. We performed displacement imaging on craniotomized mice to confirm the in situ targeting for neuromodulation site. We recorded hemodynamic responses evoked by FUS while fUS imaging revealed an ipsilateral CBV increase that peaks at 4 s post-FUS. We report a stronger hemodynamic activation in the subcortical region than cortical, showing good agreement with a brain elasticity map that can also be obtained using a similar methodology. We observed dose-dependent CBV responses with peak CBV, activated area, and correlation coefficient increasing with the ultrasonic dose. Furthermore, by mapping displacement and hemodynamic activation, we found that displacement colocalized and linearly correlated with CBV increase. The findings presented herein demonstrated that FUS evokes ipsilateral hemodynamic activation in cortical and subcortical depths while the evoked hemodynamic responses colocalize and correlate with FUS-induced displacement. We anticipate that our findings will help consolidate accurate targeting as well as shedding light on one of the mechanisms behind FUS modulation, i.e., how FUS mechanically displaces brain tissue affecting cerebral hemodynamics and thereby its associated connectivity.

Evaluation of tumor microvasculature with 3D ultrasound localization microscopy based on 2D matrix array

OBJECTIVE

Evaluation of tumor microvascular morphology is of great significance in tumor diagnosis, therapeutic effect prediction, and surgical planning. Recently, two-dimensional ultrasound localization microscopy (2DULM) has demonstrated its superiority in the field of microvascular imaging. However, it suffers from planar dependence and is unintuitive. We propose a novel three-dimensional ultrasound localization microscopy (3DULM) to avoid these limitations.

METHODS

We investigated 3DULM based on a 2D array for tumor microvascular imaging. After intravenous injection of contrast agents, all elements of the 2D array transmit and receive signals to ensure a high and stable frame rate. Microbubble signal extraction, filtering, positioning, tracking, and other processing were used to obtain a 3D vascular map, flow velocity, and flow direction. To verify the effectiveness of 3DULM, it was validated on double helix tubes and rabbit VX2 tumors. Cisplatin was used to verify the ability of 3DULM to detect microvascular changes during tumor treatment.

RESULTS

In vitro, the sizes measured by 3DULM at 3 mm and 13 mm were 178 μ m and 182 μ m , respectively. In the rabbit tumors, we acquired 9000 volumes to reveal vessels about 30 μ m in diameter, which surpasses the diffraction limit of ultrasound in traditional ultrasound imaging, and the results matched with micro-angiography. In addition, there were significant changes in vascular density and curvature between the treatment and control groups.

CONCLUSIONS

The effectiveness of 3DULM was verified in vitro and in vivo. Hence, 3DULM may have potential applications in tumor diagnosis, tumor treatment evaluation, surgical protocol guidance, and cardiovascular disease.

CLINICAL RELEVANCE STATEMENT

3D ultrasound localization microscopy is highly sensitive to microvascular changes; thus, it has clinical potential for tumor diagnosis and treatment evaluation.

KEY POINTS

• 3D ultrasound localization microscopy is demonstrated on double helix tubes and rabbit VX2 tumors. • 3D ultrasound localization microscopy can reveal vessels about 30 μ m in diameter-far smaller than traditional ultrasound. • This form of imaging has potential applications in tumor diagnosis, tumor treatment evaluation, surgical protocol guidance, and cardiovascular disease.

Deleterious effect of sustained neuroinflammation in pediatric traumatic brain injury

INTRODUCTION

Despite improved management of traumatic brain injury (TBI), it still leads to lifelong sequelae and disability, particularly in children. Chronic neuroinflammation (the so-called tertiary phase), in particular, microglia/macrophage and astrocyte reactivity, is among the main mechanisms suspected of playing a role in the generation of lesions associated with TBI. The role of acute neuroinflammation is now well understood, but its persistent effect and impact on the brain, particularly during development, are not. Here, we investigated the long-term effects of pediatric TBI on the brain in a mouse model.

METHODS

Pediatric TBI was induced in mice on postnatal day (P) 7 by weight-drop trauma. The time course of neuroinflammation and myelination was examined in the TBI mice. They were also assessed by magnetic resonance, functional ultrasound, and behavioral tests at P45.

RESULTS

TBI induced robust neuroinflammation, characterized by acute microglia/macrophage and astrocyte reactivity. The long-term consequences of pediatric TBI studied on P45 involved localized scarring astrogliosis, persistent microgliosis associated with a specific transcriptomic signature, and a long-lasting myelination defect consisting of the loss of myelinated axons, a decreased level of myelin binding protein, and severe thinning of the corpus callosum. These results were confirmed by reduced fractional anisotropy, measured by diffusion tensor imaging, and altered inter- and intra-hemispheric connectivity, measured by functional ultrasound imaging. In addition, adolescent mice with pediatric TBI showed persistent social interaction deficits and signs of anxiety and depressive behaviors.

CONCLUSIONS

We show that pediatric TBI induces tertiary neuroinflammatory processes associated with white matter lesions and altered behavior. These results support our model as a model for preclinical studies for tertiary lesions following TBI.

Power-Doppler Ultrasonic Imaging of Peripheral Perfusion in Diabetic Mice

OBJECTIVE

We explored the capabilities of power-Doppler ultrasonic (PD-US) imaging without contrast enhancement for monitoring changes in muscle perfusion over time.

METHODS

Ischemic recovery was observed in healthy and type II diabetic male and female mice with and without exercise. In separate studies, perfusion was measured during and after 5-min ischemic periods and during four-week recovery periods following irreversible femoral ligation. A goal was to assess how well PD-US estimates tracked the diabetic-related changes in endothelial function that influenced perfusion.

RESULTS

The average perfusion recovery time following femoral ligation increased 47% in diabetic males and 74% in diabetic females compared with non-diabetic mice. Flow-mediated dilation in conduit arteries and the reactive hyperemia index in resistive vessels each declined by one half in sedentary diabetic mice compared with sedentary non-diabetic mice. We found that exercise reduced the loss of endothelial function from diabetes in both sexes. The reproducibility of perfusion measurements was limited primarily by our ability to select the same region in muscle and to effectively filter tissue clutter.

CONCLUSIONS/SIGNIFICANCE

PD-US measurements can precisely follow site-specific changes in skeletal muscle perfusion related to diabetes over time, which fills the need for techniques capable of regularly monitoring atherosclerotic changes leading to ischemic vascular pathologies.

Low-rank prior-based Fast-RPCA for clutter filtering and noise suppression in non-contrast ultrasound microvascular imaging

Accurate and real-time separation of blood signal from clutter and noise signals is a critical step in clinical non-contrast ultrasound microvascular imaging. Despite the widespread adoption of singular value decomposition (SVD) and robust principal component analysis (RPCA) for clutter filtering and noise suppression, the SVD's sensitivity to threshold selection, along with the RPCA's limitations in undersampling conditions and heavy computational burden often result in suboptimal performance in complex clinical applications. To address those challenges, this study presents a novel low-rank prior-based fast RPCA (LP-fRPCA) approach to enhance the adaptability and robustness of clutter filtering and noise suppression with reduced computational cost. A low-rank prior constraint is integrated into the non-convex RPCA model to achieve a robust and efficient approximation of clutter subspace, while an accelerated alternating projection iterative algorithm is developed to improve convergence speed and computational efficiency. The performance of the LP-fRPCA method was evaluated against SVD with a tissue/blood threshold (SVD1), SVD with both tissue/blood and blood/noise thresholds (SVD2), and the classical RPCA based on the alternating direction method of multipliers algorithm through phantom and in vivo non-contrast experiments on rabbit kidneys. In the slow flow phantom experiment of 0.2 mm/s, LP-fRPCA achieved an average increase in contrast ratio (CR) of 10.68 dB, 9.37 dB, and 8.66 dB compared to SVD1, SVD2, and RPCA, respectively. In the in vivo rabbit kidney experiment, the power Doppler results demonstrate that the LP-fRPCA method achieved a superior balance in the trade-off between insufficient clutter filtering and excessive suppression of blood flow. Additionally, LP-fRPCA significantly reduced the runtime of RPCA by up to 94-fold. Consequently, the LP-fRPCA method promises to be a potential tool for clinical non-contrast ultrasound microvascular imaging.

Ultrasound Speckle Decorrelation Analysis-Based Velocimetry for 3D-Velocity-Components Measurement Using a 1D Transducer Array

Ultrasound velocimetry has been widely used for blood flow imaging. However, the flow measurements are constrained to resolve the in-plane 2D flow components when using a 1D transducer array. In this work, an ultrasound speckle decorrelation analysis-based velocimetry (3C-vUS) is proposed for 3D velocity components measurement using a 1D transducer array. The 3C-vUS theory is first derived and validated with numerical simulations and phantom experiments. The in vivo testing results show that 3C-vUS can accurately measure the blood flow 3D-velocity-components of the human carotid artery at arbitrary probe-to-vessel angles throughout the cardiac cycle. With such capability, the 3C-vUS will alleviate the requirement of operators and promote disease screening for blood flow-related disorders.

ULTRA-SR Challenge: Assessment of Ultrasound Localization and TRacking Algorithms for Super-Resolution Imaging

With the widespread interest and uptake of super-resolution ultrasound (SRUS) through localization and tracking of microbubbles, also known as ultrasound localization microscopy (ULM), many localization and tracking algorithms have been developed. ULM can image many centimeters into tissue in-vivo and track microvascular flow non-invasively with sub-diffraction resolution. In a significant community effort, we organized a challenge, Ultrasound Localization and TRacking Algorithms for Super-Resolution (ULTRA-SR). The aims of this paper are threefold: to describe the challenge organization, data generation, and winning algorithms; to present the metrics and methods for evaluating challenge entrants; and to report results and findings of the evaluation. Realistic ultrasound datasets containing microvascular flow for different clinical ultrasound frequencies were simulated, using vascular flow physics, acoustic field simulation and nonlinear bubble dynamics simulation. Based on these datasets, 38 submissions from 24 research groups were evaluated against ground truth using an evaluation framework with six metrics, three for localization and three for tracking. In-vivo mouse brain and human lymph node data were also provided, and performance assessed by an expert panel. Winning algorithms are described and discussed. The publicly available data with ground truth and the defined metrics for both localization and tracking present a valuable resource for researchers to benchmark algorithms and software, identify optimized methods/software for their data, and provide insight into the current limits of the field. In conclusion, Ultra-SR challenge has provided benchmarking data and tools as well as direct comparison and insights for a number of the state-of-the art localization and tracking algorithms.

Can 3D Multiparametric Ultrasound Imaging Predict Prostate Biopsy Outcome?

OBJECTIVES

To assess the value of 3D multiparametric ultrasound imaging, combining hemodynamic and tissue stiffness quantifications by machine learning, for the prediction of prostate biopsy outcomes.

METHODS

After signing informed consent, 54 biopsy-naïve patients underwent a 3D dynamic contrast-enhanced ultrasound (DCE-US) recording, a multi-plane 2D shear-wave elastography (SWE) scan with manual sweeping from base to apex of the prostate, and received 12-core systematic biopsies (SBx). 3D maps of 18 hemodynamic parameters were extracted from the 3D DCE-US quantification and a 3D SWE elasticity map was reconstructed based on the multi-plane 2D SWE acquisitions. Subsequently, all the 3D maps were segmented and subdivided into 12 regions corresponding to the SBx locations. Per region, the set of 19 computed parameters was further extended by derivation of eight radiomic features per parameter. Based on this feature set, a multiparametric ultrasound approach was implemented using five different classifiers together with a sequential floating forward selection method and hyperparameter tuning. The classification accuracy with respect to the biopsy reference was assessed by a group-k-fold cross-validation procedure, and the performance was evaluated by the Area Under the Receiver Operating Characteristics Curve (AUC).

RESULTS

Of the 54 patients, 20 were found with clinically significant prostate cancer (csPCa) based on SBx. The 18 hemodynamic parameters showed mean AUC values varying from 0.63 to 0.75, and SWE elasticity showed an AUC of 0.66. The multiparametric approach using radiomic features derived from hemodynamic parameters only produced an AUC of 0.81, while the combination of hemodynamic and tissue-stiffness quantifications yielded a significantly improved AUC of 0.85 for csPCa detection (p-value < 0.05) using the Gradient Boosting classifier.

CONCLUSIONS

Our results suggest 3D multiparametric ultrasound imaging combining hemodynamic and tissue-stiffness features to represent a promising diagnostic tool for biopsy outcome prediction, aiding in csPCa localization.

Functional ultrasound neuroimaging reveals mesoscopic organization of saccades in the lateral intraparietal area of posterior parietal cortex

The lateral intraparietal cortex (LIP) located within the posterior parietal cortex (PPC) is an important area for the transformation of spatial information into accurate saccadic eye movements. Despite extensive research, we do not fully understand the functional anatomy of intended movement directions within LIP. This is in part due to technical challenges. Electrophysiology recordings can only record from small regions of the PPC, while fMRI and other whole-brain techniques lack sufficient spatiotemporal resolution. Here, we use functional ultrasound imaging (fUSI), an emerging technique with high sensitivity, large spatial coverage, and good spatial resolution, to determine how movement direction is encoded across PPC. We used fUSI to record local changes in cerebral blood volume in PPC as two monkeys performed memory-guided saccades to targets throughout their visual field. We then analyzed the distribution of preferred directional response fields within each coronal plane of PPC. Many subregions within LIP demonstrated strong directional tuning that was consistent across several months to years. These mesoscopic maps revealed a highly heterogenous organization within LIP with many small patches of neighboring cortex encoding different directions. LIP had a rough topography where anterior LIP represented more contralateral upward movements and posterior LIP represented more contralateral downward movements. These results address two fundamental gaps in our understanding of LIP's functional organization: the neighborhood organization of patches and the broader organization across LIP. These findings were achieved by tracking the same LIP populations across many months to years and developing mesoscopic maps of direction specificity previously unattainable with fMRI or electrophysiology methods.

Detecting heart failure from B-mode ultrasound characterization of arterial pulse waves

This study investigated the sensitivity and specificity of identifying heart failure with reduced ejection fraction (HFrEF) from measurements of the intensity and timing of arterial pulse waves. Previously validated methods combining ultrafast B-mode ultrasound, plane-wave transmission, singular value decomposition (SVD), and speckle tracking were used to characterize the compression and decompression ("S" and "D") waves occurring in early and late systole, respectively, in the carotid arteries of outpatients with left ventricular ejection fraction (LVEF) < 40%, determined by echocardiography, and signs and symptoms of heart failure, or with LVEF ≥ 50% and no signs or symptoms of heart failure. On average, the HFrEF group had significantly reduced S-wave intensity and energy, a greater interval between the R wave of the ECG and the S wave, a reduced interval between the S and D waves, and an increase in the S-wave shift (SWS), a novel metric that characterizes the shift in timing of the S wave away from the R wave of the ECG and toward the D wave (all P < 0.01). Receiver operating characteristics (ROCs) were used to quantify for the first time how well wave metrics classified individual participants. S-wave intensity and energy gave areas under the ROC of 0.76-0.83, the ECG-S-wave interval gave 0.85-0.88, and the S-wave shift gave 0.88-0.92. Hence the methods, which are simple to use and do not require complex interpretation, provide sensitive and specific identification of HFrEF. If similar results were obtained in primary care, they could form the basis of techniques for heart failure screening.NEW & NOTEWORTHY We show that heart failure with reduced ejection fraction can be detected with excellent sensitivity and specificity in individual patients by using B-mode ultrasound to detect altered pulse wave intensity and timing in the carotid artery.

Truly Tiny Acoustic Biomolecules for Ultrasound Imaging and Therapy

Nanotechnology offers significant advantages for medical imaging and therapy, including enhanced contrast and precision targeting. However, integrating these benefits into ultrasonography is challenging due to the size and stability constraints of conventional bubble-based agents. Here bicones, truly tiny acoustic contrast agents based on gas vesicles (GVs), a unique class of air-filled protein nanostructures naturally produced in buoyant microbes, are described. It is shown that these sub-80 nm particles can be effectively detected both in vitro and in vivo, infiltrate tumors via leaky vasculature, deliver potent mechanical effects through ultrasound-induced inertial cavitation, and are easily engineered for molecular targeting, prolonged circulation time, and payload conjugation.

Motion correction of 3D dynamic contrast-enhanced ultrasound imaging without anatomical B-Mode images: Pilot evaluation in eight patients

BACKGROUND

Dynamic contrast-enhanced ultrasound (DCE-US) is highly susceptible to motion artifacts arising from patient movement, respiration, and operator handling and experience. Motion artifacts can be especially problematic in the context of perfusion quantification. In conventional 2D DCE-US, motion correction (MC) algorithms take advantage of accompanying side-by-side anatomical B-Mode images that contain time-stable features. However, current commercial models of 3D DCE-US do not provide side-by-side B-Mode images, which makes MC challenging.

PURPOSE

This work introduces a novel MC algorithm for 3D DCE-US and assesses its efficacy when handling clinical data sets.

METHODS

In brief, the algorithm uses a pyramidal approach whereby short temporal windows consisting of three consecutive frames are created to perform local registrations, which are then registered to a master reference derived from a weighted average of all frames. We applied the algorithm to imaging studies from eight patients with metastatic lesions in the liver and assessed improvements in original versus motion corrected 3D DCE-US cine using: (i) frame-to-frame volumetric overlap of segmented lesions, (ii) normalized correlation coefficient (NCC) between frames (similarity analysis), and (iii) sum of squared errors (SSE), root-mean-squared error (RMSE), and r-squared (R2) quality-of-fit from fitted time-intensity curves (TIC) extracted from a segmented lesion.

RESULTS

We noted improvements in frame-to-frame lesion overlap across all patients, from 68% ± 13% without correction to 83% ± 3% with MC (p = 0.023). Frame-to-frame similarity as assessed by NCC also improved on two different sets of time points from 0.694 ± 0.057 (original cine) to 0.862 ± 0.049 (corresponding MC cine) and 0.723 ± 0.066 to 0.886 ± 0.036 (p ≤ 0.001 for both). TIC analysis displayed a significant decrease in RMSE (p = 0.018) and a significant increase in R2 goodness-of-fit (p = 0.029) for the patient cohort.

CONCLUSIONS

Overall, results suggest decreases in 3D DCE-US motion after applying the proposed algorithm.

Ultrafast 3-D Transcutaneous Super Resolution Ultrasound Using Row-Column Array Specific Coherence-Based Beamforming and Rolling Acoustic Sub-aperture Processing: In Vitro, in Rabbit and in Human Study

OBJECTIVE

This study aimed to realise 3-D super-resolution ultrasound imaging transcutaneously with a row-column array which has far fewer independent electronic channels and a wider field of view than typical fully addressed 2-D matrix arrays. The in vivo image quality of the row-column array is generally poor, particularly when imaging non-invasively. This study aimed to develop a suite of image formation and post-processing methods to improve image quality and demonstrate the feasibility of ultrasound localisation microscopy using a row-column array, transcutaneously on a rabbit model and in a human.

METHODS

To achieve this, a processing pipeline was developed which included a new type of rolling window image reconstruction, which integrated a row-column array specific coherence-based beamforming technique with acoustic sub-aperture processing. This and other processing steps reduced the 'secondary' lobe artefacts, and noise and increased the effective frame rate, thereby enabling ultrasound localisation images to be produced.

RESULTS

Using an in vitro cross tube, it was found that the procedure reduced the percentage of 'false' locations from ∼26% to ∼15% compared to orthogonal plane wave compounding. Additionally, it was found that the noise could be reduced by ∼7 dB and the effective frame rate was increased to over 4000 fps. In vivo, ultrasound localisation microscopy was used to produce images non-invasively of a rabbit kidney and a human thyroid.

CONCLUSION

It has been demonstrated that the proposed methods using a row-column array can produce large field of view super-resolution microvascular images in vivo and in a human non-invasively.

Vasa vasorum interna in the carotid wall of active forms of Takayasu arteritis evidenced by ultrasound localization microscopy

Introduction: Takayasu arteritis (TA) is associated with microvascularization of the wall of large arteries and is related to inflammation. Ultrasound localization microscopy (ULM), combining ultrafast ultrasound imaging with microbubble (MB) injection, can track the path of MBs within the arterial wall and thus provide imaging of the vasa vasorum. From the analysis of MB tracks in the common carotid arteries of patients with active TA, we report the presence of microvessels in connection with the carotid lumen (i.e., vasa vasorum interna [VVI]). Methods: ULM maps were obtained on five patients with active disease in the observational single-center series of the TAK-UF study. MB tracks connected to the carotid lumen were automatically identified, allowing the reconstruction of VVI. Results: MB tracking allows us to observe a microvascular network on the inner part of the wall, with some vessels in communication with the carotid lumen. This type of vessel was identified in all patients with active TA (n = 5) with a median of 2.2 [1.1-3.0] vessels per acquisition (2D longitudinal view of 3 cm of the common carotid artery). The blood flow within these vessels is mainly centrifugal; that is, toward the adventitia (88% [54-100] of MB tracks with flow directed to the outer part of the wall). Conclusion: VVI are present in humans in the case of active TA and emphasize the involvement of the intima in the pathological process. ClinicalTrials.gov Identifier: NCT03956394.

Transthoracic ultrasound localization microscopy of myocardial vasculature in patients

Myocardial microvasculature and haemodynamics are indicative of potential microvascular diseases for patients with symptoms of coronary heart disease in the absence of obstructive coronary arteries. However, imaging microvascular structure and flow within the myocardium is challenging owing to the small size of the vessels and the constant movement of the patient's heart. Here we show the feasibility of transthoracic ultrasound localization microscopy for imaging myocardial microvasculature and haemodynamics in explanted pig hearts and in patients in vivo. Through a customized data-acquisition and processing pipeline with a cardiac phased-array probe, we leveraged motion correction and tracking to reconstruct the dynamics of microcirculation. For four patients, two of whom had impaired myocardial function, we obtained super-resolution images of myocardial vascular structure and flow using data acquired within a breath hold. Myocardial ultrasound localization microscopy may facilitate the understanding of myocardial microcirculation and the management of patients with cardiac microvascular diseases.

Functional Cerebral Neurovascular Mapping During Focused Ultrasound Peripheral Neuromodulation of Neuropathic Pain

BACKGROUND

Nociceptive pain is required for healthy function, yet, neuropathic pain (disease or injury) can be severely debilitating. Though a wide-array of treatment options are available, they are often systemic and/or invasive. As a promising neuromodulation treatment, Focused ultrasound (FUS) is a noninvasive and highly spatially-targeted technique shown to stimulate neural activity, yet, effects on pain signaling are currently unknown.

OBJECTIVE

Develop and validate a method for studying FUS nerve stimulation modulation of pain-evoked neural responses in vivo.

METHODS

We developed a high-resolution functional ultrasound (fUS) method capable of mapping cortical responses in healthy and neuropathic pain mice in response to FUS neuromodulation treatment.

RESULTS

FUS-evoked hemodynamic responses are correlated with the intensity of peripheral neuromodulation. We confirm functional connectivity is altered in neuropathic mice and demonstrate that FUS can modulate neuropathic pain-evoked hemodynamics.

CONCLUSIONS

The findings presented herein provides evidence for an FUS-based nerve pain method and validates the fUS technique developed for monitoring pain-evoked hemodynamics.

SIGNIFICANCE

We anticipate that the findings presented herein describe a noninvasive and flexible nerve modulation technique for pain mitigation, furthering evidence for clinical translation.

TOBE-Costas Arrays for Fast High-Resolution 3-D Power Doppler Imaging

Two-dimensional sparse arrays and row-column arrays are both alternatives to 2-D fully addressed arrays with lower channel counts. Row-column arrays have recently demonstrated fast 3-D structural and flow imaging but commonly suffer from high grating lobes or require multiplexing to achieve better quality. Two-dimensional sparse arrays enable full-volume acquisitions for each transmit event, but plane-wave transmissions with them usually lack quality in terms of uniformity of wavefronts. Here, we propose a novel architecture that combines both types of these arrays in one aperture, enabling imaging using row-column or sparse arrays alone or a hybrid imaging scheme where the row-column array is used in transmission and a 2-D sparse array in reception. This hybrid imaging scheme can potentially solve the shortcomings of each of these approaches. The sparse array layout chosen is a Costas array, characterized by having only one element per row and column, facilitating its integration with row-column arrays. We simulate images acquired with TOBE-Costas arrays using the hybrid imaging scheme and compare them to row-column and sparse spiral arrays of equivalent aperture size (128λ × 128λ at 7.5 MHz) in ultrafast plane-wave imaging of point targets and 3-D power Doppler imaging of synthetic flow phantoms. Our simulation results show that TOBE-Costas arrays exhibit superior resolution and lower sidelobe levels compared with plane-wave compounding with row-column arrays. Compared with density-tapered spiral arrays, they provide a larger field of view and finer resolution.

Photoacoustic vector tomography for deep haemodynamic imaging

Imaging deep haemodynamics non-invasively remains a quest. Although optical imaging techniques can be used to measure blood flow, they are generally limited to imaging within ∼1 mm below the skin's surface. Here we show that such optical diffusion limit can be broken through by leveraging the spatial heterogeneity of blood and its photoacoustic contrast. Specifically, successive single-shot wide-field photoacoustic images of blood vessels can be used to visualize the frame-to-frame propagation of blood and to estimate blood flow speed and direction pixel-wise. The method, which we named photoacoustic vector tomography (PAVT), allows for the quantification of haemodynamics in veins more than 5 mm deep, as we show for regions in the hands and arms of healthy volunteers. PAVT may offer advantages for the diagnosis and monitoring of vascular diseases and for the mapping of the function of the circulatory system.

Functional ultrasound imaging of the human spinal cord

Utilizing the first in-human functional ultrasound imaging (fUSI) of the spinal cord, we demonstrate the integration of spinal functional responses to electrical stimulation. We record and characterize the hemodynamic responses of the spinal cord to a neuromodulatory intervention commonly used for treating pain and increasingly used for the restoration of sensorimotor and autonomic function. We found that the hemodynamic response to stimulation reflects a spatiotemporal modulation of the spinal cord circuitry not previously recognized. Our analytical capability offers a mechanism to assess blood flow changes with a new level of spatial and temporal precision in vivo and demonstrates that fUSI can decode the functional state of spinal networks in a single trial, which is of fundamental importance for developing real-time closed-loop neuromodulation systems. This work is a critical step toward developing a vital technique to study spinal cord function and effects of clinical neuromodulation.

Time Efficient Ultrasound Localization Microscopy Based on A Novel Radial Basis Function 2D Interpolation

Ultrasound localization microscopy (ULM) allows for the generation of super-resolved (SR) images of the vasculature by precisely localizing intravenously injected microbubbles. Although SR images may be useful for diagnosing and treating patients, their use in the clinical context is limited by the need for prolonged acquisition times and high frame rates. The primary goal of our study is to relax the requirement of high frame rates to obtain SR images. To this end, we propose a new time-efficient ULM (TEULM) pipeline built on a cutting-edge interpolation method. More specifically, we suggest employing Radial Basis Functions (RBFs) as interpolators to estimate the missing values in the 2-dimensional (2D) spatio-temporal structures. To evaluate this strategy, we first mimic the data acquisition at a reduced frame rate by applying a down-sampling (DS = 2, 4, 8, and 10) factor to high frame rate ULM data. Then, we up-sample the data to the original frame rate using the suggested interpolation to reconstruct the missing frames. Finally, using both the original high frame rate data and the interpolated one, we reconstruct SR images using the ULM framework steps. We evaluate the proposed TEULM using four in vivo datasets, a Rat brain (dataset A), a Rat kidney (dataset B), a Rat tumor (dataset C) and a Rat brain bolus (dataset D), interpolating at the in-phase and quadrature (IQ) level. Results demonstrate the effectiveness of TEULM in recovering vascular structures, even at a DS rate of 10 (corresponding to a frame rate of sub-100Hz). In conclusion, the proposed technique is successful in reconstructing accurate SR images while requiring frame rates of one order of magnitude lower than standard ULM.

A New Method of Plane-Wave Ultrasound Imaging Based on Reverse Time Migration

Coherent plane-wave compounding technique enables rapid ultrasound imaging with comparable image quality to traditional B-mode imaging that relies on focused beam transmission. However, existing methods assume homogeneity in the imaged medium, neglecting the heterogeneity in sound velocities and densities present in real tissues, resulting in noise reverberation. This study introduces the Reverse Time Migration (RTM) method for ultrasound plane-wave imaging to overcome this limitation, which is combined with a method for estimating the speed of sound in layered media. Simulation results in a homogeneous background demonstrate that RTM reduces side lobes and grating lobes by approximately 30 dB, enhancing the contrast-to-noise ratio by 20% compared to conventional delay and sum (DAS) beamforming. Moreover, RTM achieves superior imaging outcomes with fewer compounding angles. The lateral resolution of the RTM with 5-9 angle compounding is able to achieve the effectiveness of the DAS method with 15-19 angle compounding, and the CNR of the RTM with 11-angle compounding is almost the same as that of the DAS with 21-angle compounding. In a heterogeneous background, experimental simulations and in vitro wire phantom experiments confirm RTM's capability to correct depth imaging, focusing reflected waves on point targets. In vitro porcine tissue experiments enable accurate imaging of layer interfaces by estimating the velocities of multiple layers containing muscle and fat. The proposed imaging procedure optimizes velocity estimation in complex media, compensates for the impact of velocity differences, provides more reliable imaging results.

Adaptive higher-order singular value decomposition clutter filter for ultrafast Doppler imaging of coronary flow under non-negligible tissue motion

BACKGROUND AND OBJECTIVE

With the development of advanced clutter-filtering techniques by singular value decomposition (SVD) and leveraging favorable acquisition settings such as open-chest imaging by a linear high-frequency probe and plane waves, several studies have shown the feasibility of cardiac flow measurements during the entire cardiac cycle, ranging from coronary flow to myocardial perfusion. When applying these techniques in a routine clinical setting, using transthoracic ultrasound imaging, new challenges emerge. Firstly, a smaller aperture is needed that can fit between ribs. Consequently, diverging waves are employed instead of plane waves to achieve an adequate field of view. Secondly, to ensure imaging at a larger depth, the maximum pulse repetition frequency has to be reduced. Lastly, in comparison to the open-chest scenario, tissue motion induced by the heartbeat is significantly stronger. The latter complicates substantially the distinction between clutter and blood signals.

METHODS

This study investigates a strategy to overcome these challenges by diverging wave imaging with an optimal number of tilt angles, in combination with dedicated clutter-filtering techniques. In particular, a novel, adaptive, higher-order SVD (HOSVD) clutter filter, which utilizes spatial, temporal, and angular information of the received ultrasound signals, is proposed to enhance clutter and blood separation.

RESULTS

When non-negligible tissue motion is present, using fewer tilt angles not only reduces the decorrelation between the received waveforms but also allows for collecting more temporal samples at a given ensemble duration, contributing to improved Doppler performance. The addition of a third angular dimension enables the application of HOSVD, providing greater flexibility in selecting blood separation thresholds from a 3-D tensor. This differs from the conventional threshold selection method in a 2-D spatiotemporal space using SVD. Exhaustive threshold search has shown a significant improvement in Contrast and Contrast-to-Noise ratio for Power Doppler images filtered with HOSVD compared to the SVD-based clutter filter.

CONCLUSION

With the improved settings, the obtained Power Doppler images show the feasibility of measuring coronary flow under the influence of non-negligible tissue motion in both in vitro and ex vivo.

Transcranial volumetric imaging using a conformal ultrasound patch

Accurate and continuous monitoring of cerebral blood flow is valuable for clinical neurocritical care and fundamental neurovascular research. Transcranial Doppler (TCD) ultrasonography is a widely used non-invasive method for evaluating cerebral blood flow1, but the conventional rigid design severely limits the measurement accuracy of the complex three-dimensional (3D) vascular networks and the practicality for prolonged recording2. Here we report a conformal ultrasound patch for hands-free volumetric imaging and continuous monitoring of cerebral blood flow. The 2 MHz ultrasound waves reduce the attenuation and phase aberration caused by the skull, and the copper mesh shielding layer provides conformal contact to the skin while improving the signal-to-noise ratio by 5 dB. Ultrafast ultrasound imaging based on diverging waves can accurately render the circle of Willis in 3D and minimize human errors during examinations. Focused ultrasound waves allow the recording of blood flow spectra at selected locations continuously. The high accuracy of the conformal ultrasound patch was confirmed in comparison with a conventional TCD probe on 36 participants, showing a mean difference and standard deviation of difference as -1.51 ± 4.34 cm s-1, -0.84 ± 3.06 cm s-1 and -0.50 ± 2.55 cm s-1 for peak systolic velocity, mean flow velocity, and end diastolic velocity, respectively. The measurement success rate was 70.6%, compared with 75.3% for a conventional TCD probe. Furthermore, we demonstrate continuous blood flow spectra during different interventions and identify cascades of intracranial B waves during drowsiness within 4 h of recording.

Context-aware deep learning enables high-efficacy localization of high concentration microbubbles for super-resolution ultrasound localization microscopy

Ultrasound localization microscopy (ULM) enables deep tissue microvascular imaging by localizing and tracking intravenously injected microbubbles circulating in the bloodstream. However, conventional localization techniques require spatially isolated microbubbles, resulting in prolonged imaging time to obtain detailed microvascular maps. Here, we introduce LOcalization with Context Awareness (LOCA)-ULM, a deep learning-based microbubble simulation and localization pipeline designed to enhance localization performance in high microbubble concentrations. In silico, LOCA-ULM enhanced microbubble detection accuracy to 97.8% and reduced the missing rate to 23.8%, outperforming conventional and deep learning-based localization methods up to 17.4% in accuracy and 37.6% in missing rate reduction. In in vivo rat brain imaging, LOCA-ULM revealed dense cerebrovascular networks and spatially adjacent microvessels undetected by conventional ULM. We further demonstrate the superior localization performance of LOCA-ULM in functional ULM (fULM) where LOCA-ULM significantly increased the functional imaging sensitivity of fULM to hemodynamic responses invoked by whisker stimulations in the rat brain.

Functional ultrasound (fUS) imaging of displacement-guided focused ultrasound (FUS) neuromodulation in mice

Focused ultrasound (FUS) stimulation is a promising neuromodulation technique with the merits of non-invasiveness, high spatial resolution, and deep penetration depth. However, simultaneous imaging of FUS-induced brain tissue displacement and the subsequent effect of FUS stimulation on brain hemodynamics has proven challenging thus far. In addition, earlier studies lack in situ confirmation of targeting except for the magnetic resonance imaging-guided FUS system-based studies. The purpose of this study is 1) to introduce a fully ultrasonic approach to in situ target, modulate neuronal activity, and monitor the resultant neuromodulation effect by respectively leveraging displacement imaging, FUS, and functional ultrasound (fUS) imaging, and 2) to investigate FUS-evoked cerebral blood volume (CBV) response and the relationship between CBV and displacement. We performed displacement imaging on craniotomized mice to confirm the in targeting for neuromodulation site. We recorded hemodynamic responses evoked by FUS and fUS revealed an ipsilateral CBV increase that peaks at 4 s post-FUS. We saw a stronger hemodynamic activation in the subcortical region than cortical, showing good agreement with the brain elasticity map that can also be obtained using a similar methodology. We observed dose-dependent CBV response with peak CBV, activated area, and correlation coefficient increasing with ultrasonic dose. Furthermore, by mapping displacement and hemodynamic activation, we found that displacement colocalizes and linearly correlates with CBV increase. The findings presented herein demonstrated that FUS evokes ipsilateral hemodynamic activation in cortical and subcortical depths and the evoked hemodynamic responses colocalized and correlate with FUS-induced displacement. We anticipate that our findings will help consolidate accurate targeting as well as an understanding of how FUS displaces brain tissue and affects cerebral hemodynamics.

Ultrafast Doppler Imaging for Early Detection of Synovitis in Rheumatoid Arthritis

OBJECTIVE

Current medical ultrasound systems possess limited sensitivity in detecting slow and weak blood flow during the early stages of rheumatoid arthritis (RA), leading to potential misdiagnosis. Ultrafast Doppler is capable of detecting slow and weak flow. This study was aimed at evaluating the diagnostic value of ultrafast Doppler for RA.

METHODS

Thirty-three RA patients (19 established, 14 early stage) and 15 healthy participants were enrolled. A programmable imaging platform with ultrafast Doppler capability was used. The benchmark was a clinical system with conventional Doppler imaging. Standardized dorsal long-axis scanning of both wrists was performed. Both ultrafast and conventional power Doppler (PD) images were quantitatively analyzed with computer assistance and semiquantitatively scored with the Outcome Measures in Rheumatology (OMERACT) scoring system.

RESULTS

Ultrafast PD revealed more blood area than conventional PD in both RA wrists and healthy wrists. Ultrafast PD OMERACT was positive in 65 of 66 RA wrists and 26 of 30 healthy wrists (sensitivity [SEN] = 0.985, accuracy [ACC] = 0.719), while conventional PD OMERACT was positive in 28 of 66 RA wrists and 0 of 30 healthy wrists (SEN = 0.424, ACC = 0.604). Ultrafast PD revealed a higher synovial PD area, dilated vessels and PD brightness in RA wrists. Peak synovial PD brightness had the best diagnostic value for RA (area under the receiver operating characteristic curve = 0.802, SEN = 0.909, ACC = 0.813). For early-stage RA patients, ultrafast peak synovial PD brightness had higher sensitivity and accuracy than conventional PD indexes.

CONCLUSION

Ultrafast PD had an increase of 0.561 in sensitivity and 0.209 in accuracy when compared with conventional PD. With its high sensitivity, ultrafast PD can detect early synovitis and identify RA patients during the early phase.

Concurrent functional ultrasound imaging with graphene-based DC-coupled electrophysiology as a platform to study slow brain signals and cerebral blood flow under control and pathophysiological brain states

Current methodology used to investigate how shifts in brain states associated with regional cerebral blood volume (CBV) change in deep brain areas, are limited by either the spatiotemporal resolution of the CBV techniques, and/or compatibility with electrophysiological recordings; particularly in relation to spontaneous brain activity and the study of individual events. Additionally, infraslow brain signals (<0.1 Hz), including spreading depolarisations, DC-shifts and infraslow oscillations (ISO), are poorly captured by traditional AC-coupled electrographic recordings; yet these very slow brain signals can profoundly change CBV. To gain an improved understanding of how infraslow brain signals couple to CBV we present a new method for concurrent CBV with wide bandwidth electrophysiological mapping using simultaneous functional ultrasound imaging (fUS) and graphene-based field effect transistor (gFET) DC-coupled electrophysiological acquisitions. To validate the feasibility of this methodology visually-evoked neurovascular coupling (NVC) responses were examined. gFET recordings are not affected by concurrent fUS imaging, and epidural placement of gFET arrays within the imaging window did not deteriorate fUS signal quality. To examine directly the impact of infra-slow potential shifts on CBV, cortical spreading depolarisations (CSDs) were induced. A biphasic pattern of decreased, followed by increased CBV, propagating throughout the ipsilateral cortex, and a delayed decrease in deeper subcortical brain regions was observed. In a model of acute seizures, CBV oscillations were observed prior to seizure initiation. Individual seizures occurred on the rising phase of both infraslow brain signal and CBV oscillations. When seizures co-occurred with CSDs, CBV responses were larger in amplitude, with delayed CBV decreases in subcortical structures. Overall, our data demonstrate that gFETs are highly compatible with fUS and allow concurrent examination of wide bandwidth electrophysiology and CBV. This graphene-enabled technological advance has the potential to improve our understanding of how infraslow brain signals relate to CBV changes in control and pathological brain states.

Transcranial Functional Ultrasound Imaging Detects Focused Ultrasound Neuromodulation Induced Hemodynamic Changes in Mouse and Nonhuman Primate Brains In Vivo

Focused ultrasound (FUS) is an emerging noinvasive technique for neuromodulation in the central nervous system (CNS). To evaluate the effects of FUS-induced neuromodulation, many studies used behavioral changes, functional magnetic resonance imaging (fMRI) or electroencephalography (EEG). However, behavioral readouts are often not easily mapped to specific brain activity, EEG has low spatial resolution limited to the surface of the brain and fMRI requires a large importable scanner that limits additional readouts and manipulations. In this context, functional ultrasound imaging (fUSI) holds promise to directly monitor the effects of FUS neuromodulation with high spatiotemporal resolution in a large field of view, with a comparatively simple and flexible setup. fUSI uses ultrafast Power Doppler Imaging (PDI) to measure changes in cerebral blood volume, which correlates well with neuronal activity and local field potentials. We designed a setup that aligns a FUS transducer with a linear array to allow immediate subsequent monitoring of the hemodynamic response with fUSI during and after FUS neuromodulation. We established a positive correlation between FUS pressure and the size of the activated area, as well as changes in cerebral blood volume (CBV) and found that unilateral sonications produce bilateral hemodynamic changes with ipsilateral accentuation in mice. We further demonstrated the ability to perform fully noninvasive, transcranial FUS-fUSI in nonhuman primates for the first time by using a lower-frequency transducer configuration.

Automatic Max-Likelihood Envelope Detection Algorithm for Quantitative High-Frame-Rate Ultrasound for Neonatal Brain Monitoring

OBJECTIVE

Post-operative brain injury in neonates may result from disturbed cerebral perfusion, but accurate peri-operative monitoring is lacking. High-frame-rate (HFR) cerebral ultrasound could visualize and quantify flow in all detectable vessels using spectral Doppler; however, automated quantification in small vessels is challenging because of low signal amplitude. We have developed an automatic envelope detection algorithm for HFR pulsed wave spectral Doppler signals, enabling neonatal brain quantitative parameter maps during and after surgery.

METHODS

HFR ultrasound data from high-risk neonatal surgeries were recorded with a custom HFR mode (frame rate = 1000 Hz) on a Zonare ZS3 system. A pulsed wave Doppler spectrogram was calculated for each pixel containing blood flow in the image, and spectral peak velocity was tracked using a max-likelihood estimation algorithm of signal and noise regions in the spectrogram, where the most likely cross-over point marks the blood flow velocity. The resulting peak systolic velocity (PSV), end-diastolic velocity (EDV) and resistivity index (RI) were compared with other detection schemes, manual tracking and RIs from regular pulsed wave Doppler measurements in 10 neonates.

RESULTS

Envelope detection was successful in both high- and low-quality arterial and venous flow spectrograms. Our technique had the lowest root mean square error for EDV, PSV and RI (0.46 cm/s, 0.53 cm/s and 0.15, respectively) when compared with manual tracking. There was good agreement between the clinical pulsed wave Doppler RI and HFR measurement with a mean difference of 0.07.

CONCLUSION

The max-likelihood algorithm is a promising approach to accurate, automated cerebral blood flow monitoring with HFR imaging in neonates.

Artificial intelligence-based speckle featurization and localization for ultrasound speckle tracking velocimetry

Deep learning-based super-resolution ultrasound (DL-SRU) framework has been successful in improving spatial resolution and measuring the velocity field information of a blood flows by localizing and tracking speckle signals of red blood cells (RBCs) without using any contrast agents. However, DL-SRU can localize only a small part of the speckle signals of blood flow owing to ambiguity problems encountered in the classification of blood flow signals from ultrasound B-mode images and the building up of suitable datasets required for training artificial neural networks, as well as the structural limitations of the neural network itself. An artificial intelligence-based speckle featurization and localization (AI-SFL) framework is proposed in this study. It includes a machine learning-based algorithm for classifying blood flow signals from ultrasound B-mode images, dimensionality reduction for featurizing speckle patterns of the classified blood flow signals by approximating them with quantitative values. A novel and robust neural network (ResSU-net) is trained using the online data generation (ODG) method and the extracted speckle features. The super-resolution performance of the proposed AI-SFL and ODG method is evaluated and compared with the results of previous U-net and conventional data augmentation methods under in silico conditions. The predicted locations of RBCs by the AI-SFL and DL-SRU for speckle patterns of blood flow are applied to a PTV algorithm to measure quantitative velocity fields of the flow. Finally, the feasibility of the proposed AI-SFL framework for measuring real blood flows is verified under in vivo conditions.

Equivalent-time-active-cavitation-imaging enables vascular-resolution blood-brain-barrier-opening-therapy planning

Objective. Linking cavitation and anatomy was found to be important for predictable outcomes in focused-ultrasound blood-brain-barrier-opening and requires high resolution cavitation mapping. However, cavitation mapping techniques for planning and monitoring of therapeutic procedures either (1) do not leverage the full resolution capabilities of ultrasound imaging or (2) place constraints on the length of the therapeutic pulse. This study aimed to develop a high-resolution technique that could resolve vascular anatomy in the cavitation map.Approach. Herein, we develop BandPass-sampled-equivalent-time-active-cavitation-imaging (BP-ETACI), derived from bandpass sampling and dual-frequency contrast imaging at 12.5 MHz to produce cavitation maps prior and during blood-brain barrier opening with long therapeutic bursts using a 1.5 MHz focused transducer in the brain of C57BL/6 mice.Main results. The BP-ETACI cavitation maps were found to correlate with the vascular anatomy in ultrasound localization microscopy vascular maps and in histological sections. Cavitation maps produced from non-blood-brain-barrier disrupting doses showed the same cavitation-bearing vasculature as maps produced over entire blood-brain-barrier opening procedures, allowing use for (1) monitoring focused-ultrasound blood-brain-barrier-opening (FUS-BBBO), but also for (2) therapy planning and target verification.Significance. BP-ETACI is versatile, created high resolution cavitation maps in the mouse brain and is easily translatable to existing FUS-BBBO experiments. As such, it provides a means to further study cavitation phenomena in FUS-BBBO.

Low-Cost Scalable PCB-Based 2-D Transducer Arrays for Volumetric Photoacoustic Imaging

Photoacoustic (PA) imaging provides deep tissue molecular imaging of chromophores with optical absorption contrast and ultrasonic resolution. Present PA imaging techniques are predominantly limited to one 2D plane per acquisition. 2D ultrasound transducers, required for real-time 3D PA imaging, are high-cost, complex to fabricate and have limited scalability in design. We present novel PCB-based 2D matrix ultrasound transducer arrays that are capable of being bulk manufactured at low-cost without using laborious ultrasound fabrication tools. The 2D ultrasound array specifications are easily scalable with respect to widely available PCB design and fabrication tools at low cost. To demonstrate scalability, we fabricated low (11 MHz) frequency 8x8 matrix array and high (40 MHz) frequency 4x4 matrix array by directly bonding an undiced polyvinylidene fluoride (PVDF) piezoelectric material of desired thickness to the custom designed PCB substrate. Characterization results demonstrate wideband PA receive sensitivity for both low (87%) and high (188%) frequency arrays. Volumetric PA imaging results of light absorbing targets inside optical scattering medium demonstrate improved spatial resolution and field of view with increase in aperture size.

Ultrasound Flow Imaging Study on Rat Brain with Ultrasound and Light Stimulations

Functional ultrasound (fUS) flow imaging provides a non-invasive method for the in vivo study of cerebral blood flow and neural activity. This study used functional flow imaging to investigate rat brain's response to ultrasound and colored-light stimuli. Male Long-Evan rats were exposed to direct full-field strobe flashes light and ultrasound stimulation to their retinas, while brain activity was measured using high-frequency ultrasound imaging. Our study found that light stimuli, particularly blue light, elicited strong responses in the visual cortex and lateral geniculate nucleus (LGN), as evidenced by changes in cerebral blood volume (CBV). In contrast, ultrasound stimulation elicited responses undetectable with fUS flow imaging, although these were observable when directly measuring the brain's electrical signals. These findings suggest that fUS flow imaging can effectively differentiate neural responses to visual stimuli, with potential applications in understanding visual processing and developing new diagnostic tools.

Quantification of brain-wide vascular resistivity via ultrafast Doppler in human neonates helps early detection of white matter injury

Preterm birth is associated with cerebrovascular development disruption and can induce white matter injuries (WMI). Transfontanellar ultrasound Doppler is the most widely used clinical imaging technique to monitor neonatal cerebral vascularisation and haemodynamics based on vascular indexes such as the resistivity index (RI); however, it has poor predictive value for brain damage. Indeed, these RI measurements are currently limited to large vessels, leading to a very limited probing of the brain's vascularisation, which may hinder prognosis. Here we show that ultrafast Doppler imaging (UfD) enables simultaneous quantification, in the whole field of view, of the local RI and vessel diameter, even in small vessels. Combining both pieces of information, we defined two new comprehensive resistivity parameters of the vascular trees. First, we showed that our technique is more sensitive in the early characterisation of the RI modifications between term and preterm neonates and for the first time we could show that the RI depends both on the vessel diameter and vascular territory. We then showed that our parameters can be used for early prediction of WMI. Our results demonstrate the potential of UfD to provide new biomarkers and pave the way for continuous monitoring of neonatal brain resistivity.

Dynamic Coronary Blood Flow Velocity and Wall Shear Stress Estimation Using Ultrasound in an Ex Vivo Porcine Heart

PURPOSE

Wall shear stress (WSS) is a critically important physical factor contributing to atherosclerosis. Mapping the spatial distribution of local, oscillatory WSS can identify important mechanisms underlying the progression of coronary artery disease.

METHODS

In this study, blood flow velocity and time-varying WSS were estimated in the left anterior descending (LAD) coronary artery of an ex vivo beating porcine heart using ultrasound with an 18 MHz linear array transducer aligned with the LAD in a forward-viewing orientation. A pulsatile heart loop with physiologically-accurate flow was created using a pulsatile pump. The coronary artery wall motion was compensated using a local block matching technique. Next, 2D and 3D velocity magnitude and WSS maps in the LAD coronary artery were estimated at different time points in the cardiac cycle using an ultrafast Doppler approach. The blood flow velocity estimated using the presented approach was compared with a commercially-available, calibrated single element blood flow velocity measurement system.

RESULTS

The resulting root mean square error (RMSE) of 2D velocity magnitude acquired from a high frequency, linear array transducer was less than 8% of the maximum velocity estimated by the commercial system.

CONCLUSION

When implemented in a forward-viewing intravascular ultrasound device, the presented approach will enable dynamic estimation of WSS, an indicator of plaque vulnerability in coronary arteries.

Time-Lagged Functional Ultrasound for Multi-Parametric Cerebral Hemodynamic Imaging

We introduce an ultrasound speckle decorrelation-based time-lagged functional ultrasound technique (tl-fUS) for the quantification of the relative changes in cerebral blood flow speed (rCBF [Formula: see text]), cerebral blood volume (rCBV) and cerebral blood flow (rCBF) during functional stimulations. Numerical simulations, phantom validations, and in vivo mouse brain experiments were performed to test the capability of tl-fUS to parse out and quantify the ratio change of these hemodynamic parameters. The blood volume change was found to be more prominent in arterioles compared to venules and the peak blood flow changes were around 2.5 times the peak blood volume change during brain activation, agreeing with previous observations in the literature. The tl-fUS shows the ability of distinguishing the relative changes of rCBFspeed, rCBV, and rCBF, which can inform specific physiological interpretations of the fUS measurements.

Four-dimensional computational ultrasound imaging of brain hemodynamics

Four-dimensional ultrasound imaging of complex biological systems such as the brain is technically challenging because of the spatiotemporal sampling requirements. We present computational ultrasound imaging (cUSi), an imaging method that uses complex ultrasound fields that can be generated with simple hardware and a physical wave prediction model to alleviate the sampling constraints. cUSi allows for high-resolution four-dimensional imaging of brain hemodynamics in awake and anesthetized mice.

Quantitative Biomarkers Derived from a Novel, Contrast-Free Ultrasound, High-Definition Microvessel Imaging for Differentiating Choroidal Tumors

Angiogenesis has an essential role in the de novo evolution of choroidal melanoma as well as choroidal nevus transformation into melanoma. Differentiating early-stage melanoma from nevus is of high clinical importance; thus, imaging techniques that provide objective information regarding tumor microvasculature structures could aid accurate early detection. Herein, we investigated the feasibility of quantitative high-definition microvessel imaging (qHDMI) for differentiation of choroidal tumors in humans. This new ultrasound-based technique encompasses a series of morphological filtering and vessel enhancement techniques, enabling the visualization of tumor microvessels as small as 150 microns and extracting vessel morphological features as new tumor biomarkers. Distributional differences between the malignant melanomas and benign nevi were tested on 37 patients with choroidal tumors using a non-parametric Wilcoxon rank-sum test, and statistical significance was declared for biomarkers with p-values < 0.05. The ocular oncology diagnosis was choroidal melanoma (malignant) in 21 and choroidal nevus (benign) in 15 patients. The mean thickness of benign and malignant masses was 1.70 ± 0.40 mm and 3.81 ± 2.63 mm, respectively. Six HDMI biomarkers, including number of vessel segments (p = 0.003), number of branch points (p = 0.003), vessel density (p = 0.03), maximum tortuosity (p = 0.001), microvessel fractal dimension (p = 0.002), and maximum diameter (p = 0.003) exhibited significant distributional differences between the two groups. Contrast-free HDMI provided noninvasive imaging and quantification of microvessels of choroidal tumors. The results of this pilot study indicate the potential use of qHDMI as a complementary tool for characterization of small ocular tumors and early detection of choroidal melanoma.

A Quantitative Method for the Evaluation of Deep Vein Thrombosis in a Murine Model Using Three-Dimensional Ultrasound Imaging

Deep vein thrombosis (DVT) is a life-threatening condition that can lead to its sequelae pulmonary embolism (PE) or post-thrombotic syndrome (PTS). Murine models of DVT are frequently used in early-stage disease research and to assess potential therapies. This creates the need for the reliable and easy quantification of blood clots. In this paper, we present a novel high-frequency 3D ultrasound approach for the quantitative evaluation of the volume of DVT in an in vitro model and an in vivo murine model. The proposed method involves the use of a high-resolution ultrasound acquisition system and semiautomatic segmentation of the clot. The measured 3D volume of blood clots was validated to be correlated with in vitro blood clot weights with an R2 of 0.89. Additionally, the method was confirmed with an R2 of 0.91 in the in vivo mouse model with a cylindrical volume from macroscopic measurement. We anticipate that the proposed method will be useful in pharmacological or therapeutic studies in murine models of DVT.

Superficial Bifurcated Microflow Phantom for High-Frequency Ultrasound Applications

OBJECTIVE

To evaluate and optimize high-frequency ultrasound (HFUS) imaging techniques that visualize the morphology of microscale vasculatures, many studies have used flow phantoms with straight channels. However, the previous phantoms lack the complexity of microvessels to simulate a realistic vascular environment in a shallow depth. This study was aimed at devising a new protocol for fabrication of a microflow phantom with bifurcated geometry at a superficial region.

METHODS

The proposed protocol involved the following features: (i) a bifurcated flow tract model 300 µm in diameter was debossed on the surface of a tissue slab made of polyvinyl alcohol cryogel, and (ii) a wall-less lumen was created via bonding tissue slabs to put a lid on the debossed flow tract. The structure of the created microflow phantom was evaluated using 2-D and 3-D power Doppler imaging with a 30 MHz HFUS modality.

RESULTS

Ultrasound imaging revealed that the desired flow tract with bifurcation was successfully created in the phantom at a depth of 2-5 mm from the ultrasound probe. The diameters of the flow tract measured in the axial direction were 307 ± 3.7 µm in the parent branch and 232 ± 18.2 and 256 ± 23.3 µm in the two daughter branches, respectively.

CONCLUSION

The experiments revealed that the proposed protocol for creating a microscale intricate flow tract with desired dimensions and depth is valid. This new phantom will facilitate further improvement in the ultrasound technologies for the precise visualization of superficial complex vasculatures such as those in skin layers.

Decoding motor plans using a closed-loop ultrasonic brain-machine interface

Brain-machine interfaces (BMIs) enable people living with chronic paralysis to control computers, robots and more with nothing but thought. Existing BMIs have trade-offs across invasiveness, performance, spatial coverage and spatiotemporal resolution. Functional ultrasound (fUS) neuroimaging is an emerging technology that balances these attributes and may complement existing BMI recording technologies. In this study, we use fUS to demonstrate a successful implementation of a closed-loop ultrasonic BMI. We streamed fUS data from the posterior parietal cortex of two rhesus macaque monkeys while they performed eye and hand movements. After training, the monkeys controlled up to eight movement directions using the BMI. We also developed a method for pretraining the BMI using data from previous sessions. This enabled immediate control on subsequent days, even those that occurred months apart, without requiring extensive recalibration. These findings establish the feasibility of ultrasonic BMIs, paving the way for a new class of less-invasive (epidural) interfaces that generalize across extended time periods and promise to restore function to people with neurological impairments.

Ultrasound Localization Microscopy for Breast Cancer Imaging in Patients: Protocol Optimization and Comparison with Shear Wave Elastography

OBJECTIVE

Ultrasound localization microscopy (ULM) has gained increasing attention in recent years because of its ability to visualize blood vessels at super-resolution. The field of oncology, in particular, could benefit from detailed vascular characterization, for example, for diagnosis and therapy monitoring. This study was aimed at refining ULM for breast cancer patients by optimizing the measurement protocol, identifying translational challenges and combining ULM and shear wave elastography.

METHODS

We computed ULM images of 11 patients with breast cancer by recording contrast-enhanced ultrasound (CEUS) sequences and post-processing them in an offline pipeline. For CEUS, two different doses and injection speeds of SonoVue were applied. The best injection protocol was determined based on quantitative parameters derived from so-called occurrence maps. In addition, a suitable measurement time window was determined, also considering the occurrence of motion. ULM results were compared with shear wave elastography and histological vessel density.

RESULTS

At the higher dose and injection speed, the highest number of microbubbles, number of tracks and vessel coverage were achieved, leading to the most detailed representation of tumor vasculature. Even at the highest concentration, no significant overlay of microbubble signals occurred. Motion significantly reduced the number of usable frames, thus limiting the measurement window to 3.5 min. ULM vessel coverage was comparable to the histological vessel fraction and correlated significantly with mean tumor elasticity.

CONCLUSION

The settings for microbubble injection strongly influence ULM images, thus requiring optimized protocols for different indications. Patient and examiner motion was identified as the main translational challenge for ULM.

Ultrasound Particle Image Velocimetry to Investigate Potential Hemodynamic Causes of Limb Thrombosis After Endovascular Aneurysm Repair With the Anaconda Device

PURPOSE

To identify potential hemodynamic predictors for limb thrombosis (LT) following endovascular aneurysm repair with the Anaconda endograft in a patient-specific phantom.

MATERIALS AND METHODS

A thin-walled flow phantom, based on a patient's aortic anatomy and treated with an Anaconda endograft, that presented with a left-sided LT was fabricated. Contrast-enhanced ultrasound particle image velocimetry was performed to quantify time-resolved velocity fields. Measurements were performed in the same phantom with and without the Anaconda endograft, to investigate the impact of the endograft on the local flow fields. Hemodynamic parameters, namely vector complexity (VC) and residence time (RT), were calculated for both iliac arteries.

RESULTS

In both limbs, the vector fields were mostly unidirectional during the peak systolic and end-systolic velocity phases before and after endograft placement. Local vortical structures and complex flow fields were observed at the diastolic and transitional flow phases. The average VC was higher (0.11) in the phantom with endograft, compared to the phantom without endograft (0.05). Notably, in both left and right iliac arteries, the anterior wall regions corresponded to a 2- and 4-fold increase in VC in the phantom with endograft, respectively. RT simulations showed values of 1.3 to 6 seconds in the phantom without endograft. A higher RT (up to 25 seconds) was observed in the phantom with endograft, in which the left iliac artery, with LT in follow-up, showed 2 fluid stasis regions.

CONCLUSION

This in vitro study shows that unfavorable hemodynamics were present mostly in the limb that thrombosed during follow-up, with the highest VC and longest RT. These parameters might be valuable in predicting the occurrence of LT in the future.

CLINICAL IMPACT

This in-vitro study aimed to identify potential hemodynamic predictors for limb thrombosis following EVAR using ultrasound particle image velocimetry (echoPIV) technique. It was shown that unfavorable hemodynamic norms were present mostly in the thrombosed limb. Owing to the in-vivo feasibility of the echoPIV, future efforts should focus on the evaluation of these hemodynamic norms in clinical trials. Thereafter, using echoPIV as a bedside technique in hospitals becomes more promising. Performing echoPIV in pre-op phase may provide valuable insights for surgeons to enhance treatment planning. EchoPIV is also applicable for follow-up sessions to evaluate treatment progress and avoid/predict complications.

Clutter filtering of angular domain data for contrast-free ultrafast microvascular imaging

Objective. Contrast-free microvascular imaging is clinically valuable for the assessment of physiological status and the early diagnosis of diseases. Effective clutter filtering is essential for microvascular visualization without contrast enhancement. Singular value decomposition (SVD)-based spatiotemporal filter has been widely used to suppress clutter. However, clinical real-time imaging relies on short ensembles (dozens of frames), which limits the implementation of SVD filtering due to the large error of eigen-correlated estimations and high dependence on optimal threshold when used in such ensembles.Approach. To address the above challenges of imaging in short ensembles, two optimized filters of angular domain data are proposed in this paper: grouped angle SVD (GA-SVD) and angular-coherence-based higher-order SVD (AC-HOSVD). GA-SVD applies SVD to the concatenation of all angular data to improve clutter rejection performance in short ensembles, while AC-HOSVD applies HOSVD to the angular data tensor and utilizes angular coherence in addition to spatial and temporal features for filtering. Feasible threshold selection strategies in each feature space are provided. The clutter rejection performance of the proposed filters and SVD was evaluated with Doppler phantom andin vivostudies at different cases. Moreover, the robustness of the filters was explored under wrong singular value threshold estimation, and their computational complexity was studied.Main results. Qualitative and quantitative results indicated that GA-SVD and AC-HOSVD can effectively improve clutter rejection performance in short ensembles, especially AC-HOSVD. Notably, the proposed methods using 20 frames had similar image quality to SVD using 100 frames.In vivostudies showed that compared to SVD, GA-SVD increased the signal-to-noise-ratio (SNR) by 6.03 dB on average, and AC-HOSVD increased the SNR by 8.93 dB on average. Furthermore, AC-HOSVD remained better power Doppler image quality under non-optimal thresholds, followed by GA-SVD.Significance. The proposed filters can greatly enhance contrast-free microvascular visualization in short ensembles and have potential for different clinical translations due to the performance differences.