Wall Shear Rate Measurement: Validation of a New Method Through Multiphysics Simulations

Adaptive Wall Shear Stress Imaging in Phantoms, Simulations and In Vivo

WSS measurement is challenging since it requires sensitive flow measurements at a distance close to the wall. The aim of this study is to develop an ultrasound imaging technique which combines vector flow imaging with an unsupervised data clustering approach that automatically detects the region close to the wall with optimally linear flow profile, to provide direct and robust WSS estimation. The proposed technique was evaluated in phantoms, mimicking normal and atherosclerotic vessels, and spatially registered Fluid Structure Interaction (FSI) simulations. A relative error of 6.7% and 19.8% was obtained for peak systolic (WSSPS) and end diastolic (WSSED) WSS in the straight phantom, while in the stenotic phantom, a good similarity was found between measured and simulated WSS distribution, with a correlation coefficient, R, of 0.89 and 0.85 for WSSPS and WSSED, respectively. Moreover, the feasibility of the technique to detect pre-clinical atherosclerosis was tested in an atherosclerotic swine model. Six swines were fed atherogenic diet, while their left carotid artery was ligated in order to disturb flow patterns. Ligated arterial segments that were exposed to low WSSPS and WSS characterized by high frequency oscillations at baseline, developed either moderately or highly stenotic plaques (p < 0.05). Finally, feasibility of the technique was demonstrated in normal and atherosclerotic human subjects. Atherosclerotic carotid arteries with low stenosis had lower WSSPS as compared to control subjects (p < 0.01), while in one subject with high stenosis, elevated WSS was found on an arterial segment, which coincided with plaque rupture site, as determined through histological examination.

Time-Resolved Wall Shear Rate Mapping Using High-Frame-Rate Ultrasound Imaging

In atherosclerosis, low wall shear stress (WSS) is known to favor plaque development, while high WSS increases plaque rupture risk. To improve plaque diagnostics, WSS monitoring is crucial. Here, we propose wall shear imaging (WASHI), a noninvasive contrast-free framework that leverages high-frame-rate ultrasound (HiFRUS) to map the wall shear rate (WSR) that relates to WSS by the blood viscosity coefficient. Our method measures WSR as the tangential flow velocity gradient along the arterial wall from the flow vector field derived using a multi-angle vector Doppler technique. To improve the WSR estimation performance, WASHI semiautomatically tracks the wall position throughout the cardiac cycle. WASHI was first evaluated with an in vitro linear WSR gradient model; the estimated WSR was consistent with theoretical values (an average error of 4.6% ± 12.4 %). The framework was then tested on healthy and diseased carotid bifurcation models. In both scenarios, key spatiotemporal dynamics of WSR were noted: 1) oscillating shear patterns were present in the carotid bulb and downstream to the internal carotid artery (ICA) where retrograde flow occurs; and 2) high WSR was observed particularly in the diseased model where the measured WSR peaked at 810 [Formula: see text] due to flow jetting. We also showed that WASHI could consistently track arterial wall motion to map its WSR. Overall, WASHI enables high temporal resolution mapping of WSR that could facilitate investigations on causal effects between WSS and atherosclerosis.

Ultrasound deep learning for monitoring of flow-vessel dynamics in murine carotid artery

Several arterial diseases are closely related with mechanical properties of the blood vessel and interactions of flow-vessel dynamics such as mean flow velocity, wall shear stress (WSS) and vascular strain. However, there is an opportunity to improve the measurement accuracy of vascular properties and hemodynamics by adopting deep learning-based ultrasound imaging for flow-vessel dynamics (DL-UFV). In this study, the DL-UFV is proposed by devising an integrated neural network for super-resolved localization and vessel wall segmentation, and it is also combined with tissue motion estimation and flow measurement techniques such as speckle image velocimetry and speckle tracking velocimetry for measuring velocity field information of blood flow. Performance of the DL-UFV is verified by comparing with other conventional techniques in tissue-mimicking phantoms. After the performance verification, in vivo feasibility is demonstrated in the murine carotid artery with different pathologies: aging and diabetes mellitus (DM). The mutual comparison of flow-vessel dynamics and histological analyses shows correlations between the immunoreactive region and abnormal flow-vessel dynamics interactions. The DL-UFV improves biases in measurements of velocity, WSS, and strain with up to 4.6-fold, 15.1-fold, and 22.2-fold in the tissue-mimicking phantom, respectively. Mean flow velocities and WSS values of the DM group decrease by 30% and 20% of those of the control group, respectively. Mean flow velocities and WSS values of the aging group (34.11 cm/s and 13.17 dyne/cm²) are slightly smaller than those of the control group (36.22 cm/s and 14.25 dyne/cm²). However, the strain values of the aging and DM groups are much smaller than those of the control group (p < 0.05). This study shows that the DL-UFV performs better than the conventional ultrasound-based flow and strain measurement techniques for measuring vascular stiffness and complicated flow-vessel dynamics. Furthermore, the DL-UFV demonstrates its excellent performance in the analysis of the hemodynamic and hemorheological effects of DM and aging on the flow and vascular characteristics. This work provides useful hemodynamic information, including mean flow velocity, WSS and strain with high-resolution for diagnosing the pathogenesis of arterial diseases. This information can be used for monitoring progression and regression of atherosclerotic diseases in clinical practice.

Tapered Vector Doppler for Improved Quantification of Low Velocity Blood Flow

A new vector velocity estimation scheme is developed, termed tapered vector Doppler (TVD), aiming to improve the accuracy of low velocity flow estimation. This is done by assessing the effects of singular value decomposition (SVD) and finite impulse response (FIR) filters and designing an estimator which accounts for signal loss due to filtering. Synthetic data created using a combination of in vivo recordings and flow simulations were used to investigate scenarios with low blood flow, in combination with true clutter motion. Using this approach, the accuracy and precision of TVD was investigated for a range of clutter-to-blood and signal-to-noise ratios. The results indicated that for the investigated carotid application and setup, the SVD filter performed as a frequency-based filter. For both SVD and FIR filters, suppression of the clutter signal resulted in large bias and variance in the estimated blood velocity magnitude and direction close to the vessel walls. Application of the proposed tapering technique yielded significant improvement in the accuracy and precision of near-wall vector velocity measurements, compared to non-TVD and weighted least squares approaches. In synthetic data, for a blood SNR of 5 dB, and in a near-wall region where the average blood velocity was 9 cm/s, the use of tapering reduced the average velocity magnitude bias from 26.3 to 1.4 cm/s. Complex flow in a carotid bifurcation was used to demonstrate the in vivo performance of TVD, and it was shown that tapering enables vector velocity estimation less affected by clutter and clutter filtering than what could be obtained by adaptive filter design only.

Determining Haemodynamic Wall Shear Stress in the Rabbit Aorta In Vivo Using Contrast-Enhanced Ultrasound Image Velocimetry

Abnormal blood flow and wall shear stress (WSS) can cause and be caused by cardiovascular disease. To date, however, no standard method has been established for mapping WSS in vivo. Here we demonstrate wide-field assessment of WSS in the rabbit abdominal aorta using contrast-enhanced ultrasound image velocimetry (UIV). Flow and WSS measurements were made independent of beam angle, curvature or branching. Measurements were validated in an in silico model of the rabbit thoracic aorta with moving walls and pulsatile flow. Mean errors over a cardiac cycle for velocity and WSS were 0.34 and 1.69%, respectively. In vivo time average WSS in a straight segment of the suprarenal aorta correlated highly with simulations (PC = 0.99) with a mean deviation of 0.29 Pa or 5.16%. To assess fundamental plausibility of the measurement, UIV WSS was compared to an analytic approximation derived from the Poiseuille equation; the discrepancy was 17%. Mapping of WSS was also demonstrated in regions of arterial branching. High time average WSS (TAWSS_xz = 3.4 Pa) and oscillatory flow (OSI_xz = 0.3) were observed near the origin of conduit arteries. In conclusion, we have demonstrated that contrast-enhanced UIV is capable of measuring spatiotemporal variation in flow velocity, arterial wall location and hence WSS in vivo with high accuracy over a large field of view.

Biomechanical Forces and Atherosclerosis: From Mechanism to Diagnosis and Treatment

The article provides an overview of current views on the role of biomechanical forces in the pathogenesis of atherosclerosis. The importance of biomechanical forces in maintaining vascular homeostasis is considered. We provide descriptions of mechanosensing and mechanotransduction. The roles of wall shear stress and circumferential wall stress in the initiation, progression and destabilization of atherosclerotic plaque are described. The data on the possibilities of assessing biomechanical factors in clinical practice and the clinical significance of this approach are presented. The article concludes with a discussion on current therapeutic approaches based on the modulation of biomechanical forces.

Quantitative assessment of the conjunctival microcirculation using a smartphone and slit-lamp biomicroscope

PURPOSE

The conjunctival microcirculation is a readily-accessible vascular bed for quantitative haemodynamic assessment and has been studied previously using a digital charge-coupled device (CCD). Smartphone video imaging of the conjunctiva, and haemodynamic parameter quantification, represents a novel approach. We report the feasibility of smartphone video acquisition and subsequent haemodynamic measure quantification via semi-automated means.

METHODS

Using an Apple iPhone 6 s and a Topcon SL-D4 slit-lamp biomicroscope, we obtained videos of the conjunctival microcirculation in 4 fields of view per patient, for 17 low cardiovascular risk patients. After image registration and processing, we quantified the diameter, mean axial velocity, mean blood volume flow, and wall shear rate for each vessel studied. Vessels were grouped into quartiles based on their diameter i.e. group 1 (<11 μm), 2 (11-16 μm), 3 (16-22 μm) and 4 (>22 μm).

RESULTS

From the 17 healthy controls (mean QRISK3 6.6%), we obtained quantifiable haemodynamics from 626 vessel segments. The mean diameter of microvessels, across all sites, was 21.1μm (range 5.8-58 μm). Mean axial velocity was 0.50mm/s (range 0.11-1mm/s) and there was a modestly positive correlation (r 0.322) seen with increasing diameter, best appreciated when comparing group 4 to the remaining groups (p < .0001). Blood volume flow (mean 145.61pl/s, range 7.05-1178.81pl/s) was strongly correlated with increasing diameter (r 0.943, p < .0001) and wall shear rate (mean 157.31 s^-1, range 37.37-841.66 s^-1) negatively correlated with increasing diameter (r - 0.703, p < .0001).

CONCLUSIONS

We, for the first time, report the successful assessment and quantification of the conjunctival microcirculatory haemodynamics using a smartphone-based system.

Continuous Simultaneous Recording of Brachial Artery Distension and Wall Shear Rate: A New Boost for Flow-Mediated Vasodilation

Vascular ultrasound has been extensively applied in the clinical setting to noninvasively assess the endothelial function by means of the so-called brachial artery flow mediated dilation (FMD). Despite the usefulness in large-scale epidemiological studies, this approach has revealed some pitfalls for assessing vascular physiology and health in individual subjects. Mainly, a reliable FMD examination should be based on the simultaneous and reliable measurement of both the stimulus, i.e., the wall shear rate (WSR), and the response, i.e., the diameter change. However, multiple technical, practical, and methodological challenges must be faced to meet this goal. In this work, we present the technical developments needed to implement a system to enable the extensive and reliable clinical ultrasound FMD examination. It integrates both a hardware part, i.e., an upgraded version of the ultrasound advanced open platform (ULA-OP), and a software part, i.e., a signal processing and data analysis platform. The system was applied for a two-center pilot clinical study on 35 young and healthy volunteers. Therefore, we present here the results of a statistical analysis on magnitude, time-course, and kinetic parameters of WSR and diameter trends that allowed us to accurately explore the vasodilatory response to the dynamic WSR changes. Our observations demonstrate that a direct and accurate estimation of WSR stimulus by multigate spectral Doppler allows understanding brachial artery vasodilatory response to reactive hyperemia. Drawing inferences on WSR stimulus from the diameter response along with an inaccurate estimation of WSR may cause further uncertainties for the accurate interpretation of the FMD response.

Data-Adaptive Coherent Demodulator for High Dynamics Pulse-Wave Ultrasound Applications

Pulse-Wave Doppler (PWD) ultrasound has been applied to the detection of blood flow for a long time; recently the same method was also proven effective in the monitoring of industrial fluids and suspensions flowing in pipes. In a PWD investigation, bursts of ultrasounds at 0.5–10 MHz are periodically transmitted in the medium under test. The received signal is amplified, sampled at tens of MHz, and digitally processed in a Field Programmable Gate Array (FPGA). First processing step is a coherent demodulation. Unfortunately, the weak echoes reflected from the fluid particles are received together with the echoes from the high-reflective pipe walls, whose amplitude can be 30–40 dB higher. This represents a challenge for the input dynamics of the system and the demodulator, which should clearly detect the weak fluid signal while not saturating at the pipe wall components. In this paper, a numerical demodulator architecture is presented capable of auto-tuning its internal dynamics to adapt to the feature of the actual input signal. The proposed demodulator is integrated into a system for the detection of the velocity profile of fluids flowing in pipes. Simulations and experiments with the system connected to a flow-rig show that the data-adaptive demodulator produces a noise reduction of at least of 20 dB with respect to different approaches, and recovers a correct velocity profile even when the input data are sampled at 8 bits only instead of the typical 12–16 bits.

Real-Time Blood Velocity Vector Measurement Over a 2-D Region

Quantitative blood velocity measurements, as currently implemented in commercial ultrasound scanners, are based on pulsed-wave (PW) spectral Doppler and are limited to detect the axial component of the velocity in a single sample volume. On the other hand, vector Doppler methods produce angle-independent estimates by, e.g., combining the frequency shifts measured from different directions. Moreover, thanks to the transmission of plane waves, the investigation of a 2-D region is possible with high temporal resolution, but, unfortunately, the clinical use of these methods is hampered by the massive calculation power required for their real-time execution. In this paper, we present a novel approach based on the transmission of plane waves and the simultaneous reception of echoes from 16 distinct subapertures of a linear array probe, which produces eight lines distributed over a 2-D region. The method was implemented on the ULAO-OP 256 research scanner and tested both in phantom and in vivo. A continuous real-time refresh rate of 36 Hz was achieved in duplex combination with a standard B-mode at pulse repetition frequency of 8 kHz. Accuracies of -11% on velocity and of 2°on angle measurements have been obtained in phantom experiments. Accompanying movies show how the method improves the quantitative measurements of blood velocities and details the flow configurations in the carotid artery of a volunteer.

Brachial artery vasodilatory response and wall shear rate determined by multigate Doppler in a healthy young cohort

Wall shear rate (WSR) is an important stimulus for the brachial artery flow-mediated dilation (FMD) response. However, WSR estimation near the arterial wall by conventional Doppler is inherently difficult. To overcome this limitation, we utilized multigate Doppler to accurately determine the WSR stimulus near the vessel wall simultaneously with the FMD response using an integrated FMD system [Ultrasound Advanced Open Platform (ULA-OP)]. Using the system, we aimed to perform a detailed analysis of WSR-FMD response and establish novel WSR parameters in a healthy young population. Data from 33 young healthy individuals (27.5 ± 4.9 yr, 19 females) were analyzed. FMD was assessed with reactive hyperemia using ULA-OP. All acquired raw data were postprocessed using custom-designed software to obtain WSR and diameter parameters. The acquired velocity data revealed that nonparabolic flow profiles within the cardiac cycle and under different flow states, with heterogeneity between participants. We also identified seven WSR magnitude and four WSR time-course parameters. Among them, WSR area under the curve until its return to baseline was the strongest predictor of the absolute ( R² = 0.25) and percent ( R² = 0.31) diameter changes in response to reactive hyperemia. For the first time, we identified mono- and biphasic WSR stimulus patterns within our cohort that produced different magnitudes of FMD response [absolute diameter change: 0.24 ± 0.10 mm (monophasic) vs. 0.17 ± 0.09 mm (biphasic), P < 0.05]. We concluded that accurate and detailed measurement of the WSR stimulus is important to comprehensively understand the FMD response and that this advance in current FMD technology could be important to better understand vascular physiology and pathology. NEW & NOTEWORTHY An estimation of wall shear rate (WSR) near the arterial wall by conventional Doppler ultrasound is inherently difficult. Using a recently developed integrated flow-mediated dilation ultrasound system, we were able to accurately estimate WSR near the wall and identified a number of novel WSR variables that may prove to be useful in the measurement of endothelial function, an important biomarker of vascular physiology and disease.