Compound Ultrasound Strain Imaging for Noninvasive Detection of (Fibro)Atheromatous Plaques: Histopathological Validation in Human Carotid Arteries

Polyvinyl Alcohol Phantoms With Heterogeneous Plaques: Estimation of Pulse Wave Velocity at the Stenotic Region Using Pulse Wave Imaging

OBJECTIVE

Plaque characterization is essential for stroke prevention. In the study reported herein, we describe a heterogeneous phantom manufacturing technique with varying plaque compositions of different stiffness using polyvinyl alcohol (PVA) to emulate stenotic arteries and evaluated the use of pulse wave imaging (PWI) to assess plaque stiffness by comparing derived pulse wave velocities, with the goal of assessing plaque vulnerability and identifying high-risk patients for stroke.

METHODS

Five stenotic phantoms (50% stenosis) were fabricated by pouring PVA solutions into 3-D-printed molds. Two of the phantoms had heterogeneous plaque compositions of soft (E0 = 13 kPa) and intermediate (E0 = 40 kPa) materials and of stiff (E0 = 54 kPa) and intermediate materials. Ultrasound sequences were acquired as the arterial phantoms were connected to a pulsating pump, and PWI was performed on the ultrasound acquisition using normalized cross-correlation to track the pulse-induced phantom wall distension propagations. Pulse wave velocities were estimated by fitting a linear regression line between the arrival time of the peak acceleration of the wall distension waveform and the corresponding location.

RESULTS

Arterial phantoms with heterogeneous plaque stiffness were successfully fabricated. Pulse wave velocities of 2.06, 2.21, 2.49, 2.67 and 3.31 m/s were found in the phantom experiments using PWI for homogeneous soft plaque, the heterogeneous soft and intermediate plaque, homogeneous intermediate plaque, the heterogeneous stiff and intermediate plaque and homogeneous stiff plaque, respectively.

CONCLUSION

A novel arterial phantom building technique was reported with varying heterogenous plaque compositions of different stiffness. The feasibility of using PWI to evaluate plaque stiffness in stenotic arteries was determined and found that PWI can distinguish between plaques of distinct stiffness and composition.

Comparing Focused-Tracked and Plane Wave-Tracked ARFI Log(VoA) In Silico and in Application to Human Carotid Atherosclerotic Plaque, Ex Vivo

A significant risk factor for ischemic stroke is carotid atherosclerotic plaque that is susceptible to rupture, with rupture potential conveyed by plaque morphology. Human carotid plaque composition and structure have been delineated noninvasively and in vivo by evaluating log(VoA), a parameter derived as the decadic log of the second time derivative of displacement induced by an acoustic radiation force impulse (ARFI). In prior work, ARFI-induced displacement was measured using conventional focused tracking; however, this requires a long data acquisition period, thereby reducing framerate. We herein evaluate if ARFI log(VoA) framerate can be increased without a reduction in plaque imaging performance using plane wave tracking instead. In silico, both focused- and plane wave-tracked log(VoA) decreased with increasing echobrightness, quantified as signal-to-noise ratio (SNR), but did not vary with material elasticity for SNRs below 40 dB. For SNRs of 40-60 dB, both focused- and plane wave-tracked log(VoA) varied with SNR and material elasticity. Above 60 dB SNR, both focused- and plane wave-tracked log(VoA) varied with material elasticity alone. This suggests that log(VoA) discriminates features according to a combination of their echobrightness and mechanical property. Further, while both focused- and plane-wave tracked log(VoA) values were artifactually inflated by mechanical reflections at inclusion boundaries, plane wave-tracked log(VoA) was more strongly impacted by off-axis scattering. Applied to three excised human cadaveric carotid plaques with spatially aligned histological validation, both log(VoA) methods detected regions of lipid, collagen, and calcium (CAL) deposits. These findings support that plane wave tracking performs comparably to focused tracking for log(VoA) imaging and that plane wave-tracked log(VoA) is a viable approach to discriminating clinically relevant atherosclerotic plaque features at a 30-fold higher framerate than by focused tracking.

Influence of key parameters on motion artifacts in lateral strain estimation with spatial angular compounding

Strain imaging can reveal the changes in tissue mechanical properties related to pathological alterations by estimating tissue strains in the lateral and axial directions of ultrasound imaging. The estimation performance in the lateral direction is usually worse than that in the axial direction. Spatial angular compounding (SAC) has been demonstrated to improve the quality of lateral estimation by deriving the lateral displacements using axial displacements obtained from multi-angle transmissions. However, motion and deformation of tissues during multiple transmissions may cause motion artifacts, and thus deteriorate the quality of strain estimation. These artifacts can be reduced by choosing appropriate imaging parameters. However, few studies have been conducted to evaluate the influences of key parameters in strain estimation, such as the pulse repetition frequency (PRF), the number of steering angles (NSA), and the maximum steering angles (MSA), in terms of performance optimization. Therefore, this study aims to investigate the effects of these parameters through simulations and phantom experiments. The performance of strain estimation is evaluated by measuring the root-mean-square error (RMSE) and the standard deviation (SD) in the simulations and phantom experiments, respectively. The contrast-to-noise ratio (CNR) of strain images is calculated in both the simulations and phantom experiments. The results show that motion artifacts in strain estimation can be reduced by increasing the PRF to 1 kHz. When the PRF reaches 1 kHz, further increase of the PRF shows little obvious improvement in strain estimation. An increase in the NSA can cause larger motion artifacts and deteriorate the quality of strain images, and the improvement of strain estimation is limited when the NSA is increased from 3 to 7. An NSA of 3 is thus recommended to balance the influences of motion artifacts and the improvement for strain estimation. The MSA has little influence on the motion artifacts, while increased MSA can achieve improved lateral estimation performance at the cost of a smaller imaging region. In light of the lateral strain estimation performance and imaging region, an MSA of 15° is recommended. The influences of these key parameters obtained from this study may provide insights for parameter optimization in strain estimation with SAC to minimize the effects of motion artifacts.

Pulse Wave Imaging Coupled With Vector Flow Mapping: A Phantom, Simulation, and In Vivo Study

Pulse wave imaging (PWI) is an ultrasound imaging modality that estimates the wall stiffness of an imaged arterial segment by tracking the pulse wave propagation. The aim of the present study is to integrate PWI with vector flow imaging, enabling simultaneous and co-localized mapping of vessel wall mechanical properties and 2-D flow patterns. Two vector flow imaging techniques were implemented using the PWI acquisition sequence: 1) multiangle vector Doppler and 2) a cross-correlation-based vector flow imaging (CC VFI) method. The two vector flow imaging techniques were evaluated in vitro using a vessel phantom with an embedded plaque, along with spatially registered fluid structure interaction (FSI) simulations with the same geometry and inlet flow as the phantom setup. The flow magnitude and vector direction obtained through simulations and phantom experiments were compared in a prestenotic and stenotic segment of the phantom and at five different time frames. In most comparisons, CC VFI provided significantly lower bias or precision than the vector Doppler method ( ) indicating better performance. In addition, the proposed technique was applied to the carotid arteries of nonatherosclerotic subjects of different ages to investigate the relationship between PWI-derived compliance of the arterial wall and flow velocity in vivo. Spearman's rank-order test revealed positive correlation between compliance and peak flow velocity magnitude ( rs = 0.90 and ), while significantly lower compliance ( ) and lower peak flow velocity magnitude ( ) were determined in older (54-73 y.o.) compared with young (24-32 y.o.) subjects. Finally, initial feasibility was shown in an atherosclerotic common carotid artery in vivo. The proposed imaging modality successfully provided information on blood flow patterns and arterial wall stiffness and is expected to provide additional insight in studying carotid artery biomechanics, as well as aid in carotid artery disease diagnosis and monitoring.

Carotid Plaque Fibrous Cap Thickness Measurement by ARFI Variance of Acceleration: In Vivo Human Results

This study evaluates the performance of an acoustic radiation force impulse (ARFI)-based outcome parameter, the decadic logarithm of the variance of acceleration, or log(VoA), for measuring carotid fibrous cap thickness. Carotid plaque fibrous cap thickness measurement by log(VoA) was compared to that by ARFI peak displacement (PD) in patients undergoing clinically indicated carotid endarterectomy using a spatially-matched histological validation standard. Fibrous caps in parametric log(VoA) and PD images were automatically segmented using a custom clustering algorithm, and a pathologist with expertise in atherosclerosis hand-delineated fibrous caps in histology. Over 10 fibrous caps, log(VoA)-derived thickness was more strongly correlated to histological thickness than PD-derived thickness, with Pearson correlation values of 0.98 for log(VoA) compared to 0.89 for PD. The log(VoA)-derived cap thickness also had better agreement with histology-measured thickness, as assessed by the concordance correlation coefficient (0.95 versus 0.62), and, by Bland-Altman analysis, was more consistent than PD-derived fibrous cap thickness. These results suggest that ARFI log(VoA) enables improved discrimination of fibrous cap thickness relative to ARFI PD and further contributes to the growing body of evidence demonstrating ARFI's overall relevance to delineating the structure and composition of carotid atherosclerotic plaque for stroke risk prediction.

Adaptive Pulse Wave Imaging: Automated Spatial Vessel Wall Inhomogeneity Detection in Phantoms and in-Vivo

Imaging arterial mechanical properties may improve vascular disease diagnosis. Pulse wave velocity (PWV) is a marker of arterial stiffness linked to cardio-vascular mortality. Pulse wave imaging (PWI) is a technique for imaging the pulse wave propagation at high spatial and temporal resolution. In this paper, we introduce adaptive PWI, a technique for the automated partition of heterogeneous arteries into individual segments characterized by most homogeneous pulse wave propagation, allowing for more robust PWV estimation. This technique was validated in a silicone phantom with a soft-stiff interface. The mean detection error of the interface was 4.67 ± 0.73 mm and 3.64 ± 0.14 mm in the stiff-to-soft and soft-to-stiff pulse wave transmission direction, respectively. This technique was tested in monitoring the progression of atherosclerosis in mouse aortas in vivo ( n = 11 ). The PWV was found to already increase at the early stage of 10 weeks of high-fat diet (3.17 ± 0.67 m/sec compared to baseline 2.55 ± 0.47 m/sec, ) and further increase after 20 weeks of high-fat diet (3.76±1.20 m/sec). The number of detected segments of the imaged aortas monotonically increased with the duration of high-fat diet indicating an increase in arterial wall property inhomogeneity. The performance of adaptive PWI was also tested in aneurysmal mouse aortas in vivo. Aneurysmal boundaries were detected with a mean error of 0.68±0.44 mm. Finally, initial feasibility was shown in the carotid arteries of healthy and atherosclerotic human subjects in vivo ( n = 3 each). Consequently, adaptive PWI was successful in detecting stiffness inhomogeneity at its early onset and monitoring atherosclerosis progression in vivo.

Delineation of Human Carotid Plaque Features In Vivo by Exploiting Displacement Variance

While in vivo acoustic radiation force impulse (ARFI)-induced peak displacement (PD) has been demonstrated to have high sensitivity and specificity for differentiating soft from stiff plaque components in patients with carotid plaque, the parameter exhibits poorer performance for distinguishing between plaque features with similar stiffness. To improve discrimination of carotid plaque features relative to PD, we hypothesize that signal correlation and signal-to-noise ratio (SNR) can be combined, outright or via displacement variance. Plaque feature detection by displacement variance, evaluated as the decadic logarithm of the variance of acceleration and termed "log(VoA)," was compared to that achieved by exploiting SNR, cross correlation coefficient, and ARFI-induced PD outcome metrics. Parametric images were rendered for 25 patients undergoing carotid endarterectomy, with spatially matched histology confirming plaque composition and structure. On average, across all plaques, log(VoA) was the only outcome metric with values that statistically differed between regions of lipid-rich necrotic core (LRNC), intraplaque hemorrhage (IPH), collagen (COL), and calcium (CAL). Further, log(VoA) achieved the highest contrast-to-noise ratio (CNR) for discriminating between LRNC and IPH, COL and CAL, and grouped soft (LRNC and IPH) and stiff (COL and CAL) plaque components. More specifically, relative to the previously demonstrated ARFI PD parameter, log(VoA) achieved 73% higher CNR between LRNC and IPH and 59% higher CNR between COL and CAL. These results suggest that log(VoA) enhances the differentiation of LRNC, IPH, COL, and CAL in human carotid plaques, in vivo, which is clinically relevant to improving stroke risk prediction and medical management.

Non-contrast agent based small vessel imaging of human thyroid using motion corrected power Doppler imaging

Singular value based spatiotemporal clutter filtering (SVD-STF) can significantly improve the sensitivity of blood flow imaging in small vessels without using contrast agents. However, despite effective clutter filtering, large physiological motion in thyroid imaging can impact coherent integration of the Doppler signal and degrade the visualization of the underlying vasculature. In this study, we hypothesize that motion correction of the clutter filtered Doppler ensemble, prior to the power Doppler estimation, can considerably improve the visualization of smalls vessels in suspicious thyroid nodules. We corroborated this hypothesis by conducting in vivo experiments on 10 female patients in the age group 44-82 yrs, with at least one thyroid nodule suspicious of malignancy, with recommendation for fine needle aspiration biopsy. Ultrasound images were acquired using a clinical ultrasound scanner, implemented with compounded plane wave imaging. Axial and lateral displacements associated with the thyroid nodules were estimated using 2D normalized cross-correlation. Subsequently, the tissue clutter associated with the Doppler ensemble was suppressed using SVD-STF. Motion correction of the clutter-filtered Doppler ensemble was achieved using a spline based sub-pixel interpolation. The results demonstrated that power Doppler images of thyroid nodules were noticeably degraded due to large physiological motion of the pulsating carotid artery in the proximity. The resultant power Doppler images were corrupted with signal distortion, motion blurring and occurrence of artificial shadow vessels and displayed visibly low signal-to-background contrast. In contrast, the power Doppler images obtained from the motion corrected ultrasound data addressed the issue and considerabley improved the visualization of blood flow. The signal-to-noise ratio and the contrast-to-noise ratio increased by up to 15.2 dB and 12.1 dB, respectively. Across the ten subjects, the highest improvement was observed for the nodule with the largest motion. These preliminary results show the ability of using motion correction to improve the visualization of small vessel blood flow in thyroid, without using any contrast agents. The results of this feasibility study were encouraging, and warrant further development and more in vivo validation in moving tissues and organs.

Simultaneous Vascular Strain and Blood Vector Velocity Imaging Using High-Frequency Versus Conventional-Frequency Plane Wave Ultrasound: A Phantom Study

Plaque strain and blood vector velocity imaging of stenosed arteries are expected to aid in diagnosis and prevention of cerebrovascular disease. Ultrafast plane wave imaging enables simultaneous strain and velocity estimation. Multiple ultrasound vendors are introducing high-frequency ultrasound probes and systems. This paper investigates whether the use of high-frequency ultrafast ultrasound is beneficial for assessing blood velocities and strain in arteries. The performance of strain and blood flow velocity estimation was compared between a high-frequency transducer (MS250, fc = 21 MHz) and a clinically utilized transducer (L12-5, fc = 9 MHz). Quantitative analysis based on straight tube phantom experiments revealed that the MS250 outperformed the L12-5 in the superficial region: low velocities near the wall were more accurately estimated and wall strains were better resolved. At greater than 2-cm echo depth, the L12-5 performed better due to the high attenuation of the MS250 probe. Qualitative comparison using a perfused patient-specific carotid bifurcation phantom confirmed these findings. Thus, in conclusion, for strain and blood velocity estimation for depths up to ~2 cm, a high-frequency probe is recommended.

Visualizing Angle-Independent Principal Strains in the Longitudinal View of the Carotid Artery: Phantom and In Vivo Evaluation

Non-invasive vascular elastography can evaluate the stiffness of the carotid artery by visualizing the vascular strain distribution. Axial strain estimates of the longitudinal cross section of the carotid artery are sensitive to the angle between the artery and the transducer. Anatomical variations in branching and arching of the carotid artery can affect the assessment of arterial stiffness. In this study, we hypothesized that principal strain elastograms computed using compounded plane wave imaging can reliably visualize the strain distribution in the carotid artery, independent of the transducer angle. We corroborated this hypothesis by conducting phantom and in vivo studies using a commercial ultrasound scanner (Sonix RP, Ultrasonix Medical Corp., Richmond, BC, Canada). The phantom studies were conducted using a homogeneous cryogel vessel phantom. The goal of the phantom study was to assess the feasibility of visualizing the radial deformation in the longitudinal plane of the vessel phantom, independent of the transducer angle (±30°, ±20°, ±10° and 0°). The in vivo studies were conducted on 20 healthy human volunteers in the age group 50-60 y. All echo imaging was performed at a transmit frequency of 5 MHz and sampling frequency of 40 MHz. The elastograms obtained from the phantom study revealed that for straight vessels, which had their lumen parallel to the transducer, principal strains were similar to axial strains. At non-parallel configurations (angles ±30°, ±20° and ±10°), the magnitudes of the mean principal strains were within 2.5% of the parallel configuration (0° angle) estimates and, thus, were observed to be relatively unaffected by change in angle. However, in comparison, the magnitude of the axial strain decreased with increase in angle because of coordinate dependency. Further, the pilot in vivo study indicated that the principal and axial strain elastograms were similar for subjects with relatively straight arteries. However, for arteries with arched geometry, axial strains were significantly lower (p <0.01) than the corresponding principal vascular strains, which was consistent with the results obtained from the phantom study. In conclusion, the results of the phantom and in vivo studies revealed that principal strain elastograms computed using CPW imaging could reliably visualize angle-independent vascular strains in the longitudinal plane of the carotid artery.