Shear Wave Elastography Quantifies Stiffness in Ex Vivo Porcine Artery with Stiffened Arterial Region

Performance evaluation of commercial and non-commercial shear wave elastography implementations for vascular applications

BACKGROUND

Shear wave elastography (SWE) is mainly used for stiffness estimation of large, homogeneous tissues, such as the liver and breasts. However, little is known about its accuracy and applicability in thin (∼0.5-2 mm) vessel walls. To identify possible performance differences among vendors, we quantified differences in measured wave velocities obtained by commercial SWE implementations of various vendors over different imaging depths in a vessel-mimicking phantom. For reference, we measured SWE values in the cylindrical inclusions and homogeneous background of a commercial SWE phantom. Additionally, we compared the accuracy between a research implementation and the commercially available clinical SWE on an Aixplorer ultrasound system in phantoms and in vivo in patients.

METHODS

SWE measurements were performed over varying depths (0-35 mm) using three ultrasound machines with four ultrasound probes in the homogeneous 20 kPa background and cylindrical targets of 10, 40, and 60 kPa of a multi-purpose phantom (CIRS-040GSE) and in the anterior and posterior wall of a homogeneous polyvinyl alcohol vessel-mimicking phantom. These phantom data, along with in vivo SWE data of carotid arteries in 23 patients with a (prior) head and neck neoplasm, were also acquired in the research and clinical mode of the Aixplorer ultrasound machine. Machine-specific estimated phantom stiffness values (CIRS phantom) or wave velocities (vessel phantom) over all depths were visualized, and the relative error to the reference values and inter-frame variability (interquartile range/median) were calculated. Correlations between SWE values and target/vessel wall depth were explored in phantoms and in vivo using Spearman's correlations. Differences in wave velocities between the anterior and posterior arterial wall were assessed with Wilcoxon signed-rank tests. Intra-class correlation coefficients were calculated for a sample of ten patients as a measure of intra- and interobserver reproducibility of SWE analyses in research and clinical mode.

RESULTS

There was a high variability in obtained SWE values among ultrasound machines, probes, and, in some cases, with depth. Compared to the homogeneous CIRS-background, this variation was more pronounced for the inclusions and the vessel-mimicking phantom. Furthermore, higher stiffnesses were generally underestimated. In the vessel-mimicking phantom, anterior wave velocities were (incorrectly) higher than posterior wave velocities (3.4-5.6 m/s versus 2.9-5.9 m/s, p ≤ 0.005 for 3/4 probes) and remarkably correlated with measurement depth for most machines (Spearman's ρ = -0.873-0.969, p < 0.001 for 3/4 probes). In the Aixplorer's research mode, this difference was smaller (3.3-3.9 m/s versus 3.2-3.6 m/s, p = 0.005) and values did not correlate with measurement depth (Spearman's ρ = 0.039-0.659, p ≥ 0.002). In vivo, wave velocities were higher in the posterior than the anterior vessel wall in research (left p = 0.001, right p < 0.001) but not in clinical mode (left: p = 0.114, right: p = 0.483). Yet, wave velocities correlated with vessel wall depth in clinical (Spearman's ρ = 0.574-0.698, p < 0.001) but not in research mode (Spearman's ρ = -0.080-0.466, p ≥ 0.003).

CONCLUSIONS

We observed more variation in SWE values among ultrasound machines and probes in tissue with high stiffness and thin-walled geometry than in low stiffness, homogeneous tissue. Together with a depth-correlation in some machines, where carotid arteries have a fixed location, this calls for caution in interpreting SWE results in clinical practice for vascular applications.

On the Role of Coherent Plane Wave Compounding in Shear Wave Elasticity Imaging: The Convolution Effect and Its Implications

OBJECTIVE

The clinical applicability of shear wave elasticity imaging (SWEI) has been confounded by its appreciable inter-system variability and unsatisfactory sensitivity. While SWEI relies on plane wave imaging (PWI) to achieve real-time rendering, it has been rarely noticed that PWI can affect SWEI's performance. This work is aimed at demonstrating that the use of coherent plane wave compounding (CPWC) can be a factor causing SWEI's underperformance.

METHODS

We presented a model to formally describe the slow-time behavior of CPWC in motion tracking. This model reveals that CPWC introduces temporal convolution on the observed motion, making the motion sampling process a low-pass filter (LPF). For validation, shear waves were produced in a phantom in the same way but sampled via PWI using different compounding numbers (CN) and pulse repetition frequencies (PRF), with the obtained signals compared with the inferences drawn from our model. Similar experiments were performed to reconstruct two small targets in the phantom in order to appraise the impact of CPWC on SWEI's sensitivity.

DISCUSSION

The validation experiment shows that the measurements match well with the model inferences, which verifies the LPF nature of CPWC. The phantom study also shows that either increasing CN or decreasing PRF can cause the loss of high-frequency motion information, leading to blurred target delineation by SWEI.

CONCLUSION

The convolution effect can help understand the variability of SWEI. Researchers should beware this effect when working on SWEI standardization. Clinicians using SWEI should also be cautious because this effect makes it harder to identify small lesions.

Extending arterial stiffness assessment along the circumference using beam-steered ARFI and wave-tracking: A proof-of-principle study in phantoms and ex vivo

Background

To fully quantify arterial wall and plaque stiffness, acoustic radiation force impulse (ARFI)-induced wave-tracking along the entire vessel circumference is desired. However, attenuation and guided wave behavior in thin vessel walls limits wave-tracking to short trajectories. This study investigated the potential of beam-steered ARFI and wave-tracking to extend group velocity estimation over a larger proportion of the circumference compared to conventional 0° ARFI-induced wave-tracking.

Methods

Seven vessel-mimicking polyvinyl alcohol cryogel phantoms with various dimensions and compositions and an ex vivo human carotid artery were imaged in a dedicated setup. For every 20⁰ phantom rotation, transverse group wave velocity measurements were performed with an Aixplorer Ultimate system and SL18-5 transducer using 0⁰/20⁰/-20⁰-angled ultrasound pushes. Transmural angular wave velocities were derived along 60⁰-trajectories. A 360⁰-angular velocity map was composed from the top-wall 60⁰-trajectories 0°-data, averaged over all physical phantom rotations (reference). For each phantom rotation, 360⁰-angular velocity maps were composed using 0°-data (0⁰-approach) or data from all angles (beam-steered approach). Percentages of rotations with visible waves and relative angular velocity errors compared to the reference map as function of the circumferential angle were determined for both approaches.

Results

Reference 360°-angular velocity maps could be derived for all samples, representing their stiffness. Beam-steering decreased the proportion of the circumference where waves were untraceable by 20% in phantoms and 10% ex vivo, mainly at 0° push locations. Relative errors were similar for both approaches (phantoms: 10-15%, ex vivo: 15-35%).

Conclusion

Beam-steering enables wave-tracking along a higher proportion of the wall circumference than 0⁰ ARFI-induced wave-tracking.

Effects of acute cervical stretching on arterial wall elastic properties

Purpose: Acute (immediate) or regular (mid- or long-term) stretching increases arterial compliance and reduces arterial stiffness. Stretching is widely known to induce arterial functional factor changes, but it is unclear whether stretching alters arterial structural factors. Ultrasound shear wave elastography can quantify the distribution of tissue elastic properties as an index of arterial structural factors. This study thus aimed to examine the effects of acute cervical stretching on arterial wall tissue elastic properties. Methods: Seventeen healthy young adults participated in two different trials for 15 min in random order on separate days: a resting and sitting trial (CON) and a supervised cervical stretching trial (CS). In CS, subjects performed 10 different stretches. At each site, the stretch was held for 30 s followed by a 10-s relaxation period. In CON, subjects rested on a chair for 15 min. Results: After the experiment, carotid arterial compliance, assessed by combined ultrasound imaging and applanation tonometry, was significantly increased in CS, but not in CON. However, there was no significant change in tissue elasticity properties of the arterial wall in either trial, as assessed by ultrasound shear wave elastography. Conclusion: Acute cervical stretching significantly increased carotid artery compliance in young participants, but did not reduce elastic tissue properties (i.e., arterial structural factors) of the carotid artery wall. These results strongly suggest that changes in structural factors have little relation to stretching-induced acute increases in arterial compliance.

Vascular wall motion detection models based on long short-term memory in plane-wave-based ultrasound imaging

Objective.Vascular wall motion can be used to diagnose cardiovascular diseases. In this study, long short-term memory (LSTM) neural networks were used to track vascular wall motion in plane-wave-based ultrasound imaging.Approach.The proposed LSTM and convolutional LSTM (ConvLSTM) models were trained using ultrasound data from simulations and tested experimentally using a tissue-mimicking vascular phantom and anin vivostudy using a carotid artery. The performance of the models in the simulation was evaluated using the mean square error from axial and lateral motions and compared with the cross-correlation (XCorr) method. Statistical analysis was performed using the Bland-Altman plot, Pearson correlation coefficient, and linear regression in comparison with the manually annotated ground truth.Main results.For thein vivodata, the median error and 95% limit of agreement from the Bland-Altman analysis were (0.01, 0.13), (0.02, 0.19), and (0.03, 0.18), the Pearson correlation coefficients were 0.97, 0.94, and 0.94, respectively, and the linear equations were 0.89x+ 0.02, 0.84x+ 0.03, and 0.88x+ 0.03 from linear regression for the ConvLSTM model, LSTM model, and XCorr method, respectively. In the longitudinal and transverse views of the carotid artery, the LSTM-based models outperformed the XCorr method. Overall, the ConvLSTM model was superior to the LSTM model and XCorr method.Significance.This study demonstrated that vascular wall motion can be tracked accurately and precisely using plane-wave-based ultrasound imaging and the proposed LSTM-based models.

Full waveform inversion for arterial viscoelasticity

Objective. Arterial viscosity is emerging as an important biomarker, in addition to the widely used arterial elasticity. This paper presents an approach to estimate arterial viscoelasticity using shear wave elastography (SWE).Approach. While dispersion characteristics are often used to estimate elasticity from SWE data, they are not sufficiently sensitive to viscosity. Driven by this, we develop a full waveform inversion (FWI) methodology, based on directly matching predicted and measured wall velocity in space and time, to simultaneously estimate both elasticity and viscosity. Specifically, we propose to minimize an objective function capturing the correlation between measured and predicted responses of the anterior wall of the artery.Results. The objective function is shown to be well-behaving (generally convex), leading us to effectively use gradient optimization to invert for both elasticity and viscosity. The resulting methodology is verified with synthetic data polluted with noise, leading to the conclusion that the proposed FWI is effective in estimating arterial viscoelasticity.Significance. Accurate estimation of arterial viscoelasticity, not just elasticity, provides a more precise characterization of arterial mechanical properties, potentially leading to a better indicator of arterial health.

Safety of arterial shear wave elastography-ex-vivoassessment of induced strain and strain rates

Shear wave elastography (SWE) is a promising technique for characterizing carotid plaques and assessing local arterial stiffness. The mechanical stress to which the tissue is subjected during SWE using acoustic radiation force (ARF), leading to strain at a certain strain rate, is still relatively unknown. Because SWE is increasingly used for arterial applications where the mechanical stress could potentially lead to significant consequences, it is important to understand the risks of SWE- induced strain and strain rate. The aim of this study was to investigate the safety of SWE in terms of induced arterial strain and strain rateex-vivoand in a human carotid arteryin-vivo. SWE was performed on six porcine aortae as a model of the human carotid artery using different combinations of ARF push parameters (push voltage: 60/90 V, aperture width: f/1.0/1.5, push length: 100/150/200 μs) and distance to push position. The largest induced strain and strain rate were 1.46 % and 54 s^-1(90 V, f/1.0, 200 μs), respectively. Moreover, the SWE-induced strains and strain rates increased with increasing push voltage, aperture, push length, and decreasing distance between the region of interest and the push. In the human carotid artery, the SWE-induced maximum strain was 0.06 % and the maximum strain rate was 1.58 s^-1, compared with the maximum absolute strain and strain rate of 12.61 % and 5.12 s^-1, respectively, induced by blood pressure variations in the cardiac cycle. Our results indicate thatex-vivoarterial SWE does not expose the artery to higher strain rate than normal blood pressure variations, and to strain one order of magnitude higher than normal blood pressure variations, at the push settings and distances from the region of interest used in this study.

The combined importance of finite dimensions, anisotropy, and pre-stress in acoustoelastography

Dynamic elastography, whether based on magnetic resonance, ultrasound, or optical modalities, attempts to reconstruct quantitative maps of the viscoelastic properties of biological tissue, properties that are altered by disease and injury, by noninvasively measuring mechanical wave motion in the tissue. Most reconstruction strategies that have been developed neglect boundary conditions, including quasistatic tensile or compressive loading resulting in a nonzero prestress. Significant prestress is inherent to the functional role of some biological tissues currently being studied using elastography, such as skeletal and cardiac muscle, arterial walls, and the cornea. In the present article, we review how prestress alters both bulk mechanical wave motion and wave motion in one- and two-dimensional waveguides. Key findings are linked to studies on skeletal muscle and the human cornea, as one- and two-dimensional waveguide examples. This study highlights the underappreciated combined acoustoelastic and waveguide challenge to elastography. Can elastography truly determine viscoelastic properties of a material when what it is measuring is affected by both these material properties and unknown prestress and other boundary conditions?

Microscopic multifrequency magnetic resonance elastography of ex vivo abdominal aortic aneurysms for extracellular matrix imaging in a mouse model

An abdominal aortic aneurysm (AAA) is a permanent dilatation of the abdominal aorta, usually accompanied by thrombus formation. The current clinical imaging modalities cannot reliably visualize the thrombus composition. Remodeling of the extracellular matrix (ECM) during AAA development leads to stiffness changes, providing a potential imaging marker. 14 apolipoprotein E-deficient mice underwent surgery for angiotensin II-loaded osmotic minipump implantation. 4 weeks post-op, 5 animals developed an AAA. The aneurysm was imaged ex vivo by microscopic multifrequency magnetic resonance elastography (µMMRE) with an in-plane resolution of 40 microns. Experiments were performed on a 7-Tesla preclinical magnetic resonance imaging scanner with drive frequencies between 1000 Hz and 1400 Hz. Shear wave speed (SWS) maps indicating stiffness were computed based on tomoelastography multifrequency inversion. As control, the aortas of 5 C57BL/6J mice were examined with the same imaging protocol. The regional variation of SWS in the thrombus ranging from 0.44 ± 0.07 to 1.20 ± 0.31 m/s was correlated fairly strong with regional histology-quantified ECM accumulation (R² = 0.79). Our results suggest that stiffness changes in aneurysmal thrombus reflect ECM remodeling, which is critical for AAA risk assessment. In the future, µMMRE could be used for a mechanics-based clinical characterization of AAAs in patients. STATEMENT OF SIGNIFICANCE: To our knowledge, this is the first study mapping the stiffness of abdominal aortic aneurysms with microscopic resolution of 40 µm. Our work revealed that stiffness critically changes due to extracellular matrix (ECM) remodeling in the aneurysmal thrombus. We were able to image various levels of ECM remodeling in the aneurysm reflected in distinct shear wave speed patterns with a strong correlation to regional histology-quantified ECM accumulation. The generated results are significant for the application of microscopic multifrequency magnetic resonance elastography for quantification of pathological remodeling of the ECM and may be of great interest for detailed characterization of AAAs in patients.

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.

Cylindrical Transducer Array for Intravascular Shear Wave Elasticity Imaging: Preliminary Development

We present an intravascular ultrasound (IVUS) transducer array designed to enable shear wave elasticity imaging (SWEI) of arteries for the detection and characterization of atherosclerotic soft plaques. Using a custom dicing fixture, we have fabricated single-element and axially-segmented array transducer prototypes from 4.6-Fr to 7.6-Fr piezoceramic tubes, respectively. Focused excitation of the array prototype at 4 MHz yielded a focal gain of 5× in intensity, for an estimated 60 mW/cm² [Formula: see text] and 1.6-MPa negative peak pressure at 4.5-mm range in water. The single-element transducer generated a peak radial displacement of [Formula: see text] in a uniform elasticity phantom, with axial shear waves detectable by an external linear array probe up to 5 mm away from the excitation plane. In a vessel phantom with a soft inclusion, the array prototype generated peak displacements of 2.2 and [Formula: see text] in the soft inclusion and vessel wall regions, respectively. A SWEI image of the vessel phantom was reconstructed, with measured shear wave speed (SWS) of 1.66 ± 0.91 m/s and 0.97 ± 0.59 m/s for the soft inclusion and vessel wall regions, respectively. The array prototype was also used to obtain a SWEI image of an ex vivo porcine artery, with a mean SWS of 3.97 ± 1.12 m/s. These results suggest that a cylindrical intravascular ultrasound (IVUS) transducer array could be made capable of SWEI for atherosclerotic plaque detection in coronary arteries.

Evaluation of carotid artery stiffness in patients with coronary artery disease using acoustic radiation force impulse elastography

OBJECTIVES

We evaluated carotid artery stiffness in patients with coronary artery disease (CAD) using acoustic radiation force impulse (ARFI) elastography and investigated the relationship between stiffness and CAD.

METHODS

This study examined 76 CAD patients (aged 60.7 ± 11) and 70 healthy individuals (aged 59.6 ± 9). The left common carotid artery mean shear wave velocity (LCmeanSWV) and the right common carotid artery mean shear wave velocity (RCmeanSWV) of the anterior walls were measured using ARFI elastography, and the results of the patient group and the healthy group were compared. The common carotid intima-media thickness (CIMT) was measured in both groups and compared with mean SWV.

RESULTS

The RCmeanSWVs in the patient and healthy groups were 3.47 ± 1 m/s and 2.69 ± 0.90 m/s, respectively (p < 0.001). The LCmeanSWVs in the patient and healthy groups were 3.60 ± 0.9 m/s and 2.90 ± 0.80 m/s, respectively (p < 0.001). A significant correlation was found between the RCmeanSWV and the right CIMT values and between the LCmeanSWV and the left CIMT values (r = 0.231, p = 0.03 and r = 0.331, p = 0.002, respectively).

CONCLUSIONS

The mean SWV values of carotid arteries of CAD patients measured with ARFI elastography were significantly higher than the mean SWV values of the carotid arteries of healthy individuals. Thus, measurement of the carotid artery SWV could be a potential tool in the risk evaluation of cardiovascular disease. Nonetheless, new studies are required to determine whether this method serves as a useful additional tool.

Ultrasound shear wave elastography imaging of common carotid arteries in patients with Spontaneous Coronary Artery Dissection (SCAD)

BACKGROUND

Shear wave elastography (SWE) is emerging as a valuable clinical tool for a variety of conditions. The aim of this pilot study was to assess the potential of SWE imaging of the common carotid arteries (CCA) in patients with spontaneous coronary artery dissection (SCAD), a rare but potentially life-threatening condition, hypothesized to be linked to changes in vessel wall elasticity.

METHODS

Ultrasound shear wave elastography (SWE) estimates of artery wall elasticity were obtained from the left and right CCAs of 89 confirmed SCAD patients and 38 non-dissection controls. SWE images obtained over multiple cardiac cycles were analysed by a blinded observer to estimate elasticity in the form of a Young's Modulus (YM) value, across regions of interest (ROI) located within the anterior and posterior CCA walls.

RESULTS

YM estimates ranged from 17 to 133 kPa in SCAD patients compared to 34 to 87 kPa in non-dissection controls. The mean YM of 55 [standard deviation (SD): 21] kPa in SCAD patients was not significantly different to the mean of 57 [SD: 12] kPa in controls, p = 0.32. The difference between groups was 2 kPa [95% Confidence Interval - 11, 4].

CONCLUSIONS

SWE imaging of CCAs in SCAD patients is feasible although the clinical benefit is limited by relatively high variability of YM values which may have contributed to our finding of no significant difference between SCAD patients and non-dissection controls.

Wave-based optical coherence elastography: The 10-year perspective

After 10 years of progress and innovation, optical coherence elastography (OCE) based on the propagation of mechanical waves has become one of the major and the most studied OCE branches, producing a fundamental impact in the quantitative and nondestructive biomechanical characterization of tissues. Preceding previous progress made in ultrasound and magnetic resonance elastography; wave-based OCE has pushed to the limit the advance of three major pillars: (1) implementation of novel wave excitation methods in tissues, (2) understanding new types of mechanical waves in complex boundary conditions by proposing advance analytical and numerical models, and (3) the development of novel estimators capable of retrieving quantitative 2D/3D biomechanical information of tissues. This remarkable progress promoted a major advance in answering basic science questions and the improvement of medical disease diagnosis and treatment monitoring in several types of tissues leading, ultimately, to the first attempts of clinical trials and translational research aiming to have wave-based OCE working in clinical environments. This paper summarizes the fundamental up-to-date principles and categories of wave-based OCE, revises the timeline and the state-of-the-art techniques and applications lying in those categories, and concludes with a discussion on the current challenges and future directions, including clinical translation research.

Toward improved accuracy in shear wave elastography of arteries through controlling the arterial response to ultrasound perturbation in-silico and in phantoms

Dispersion-based inversion has been proposed as a viable direction for materials characterization of arteries, allowing clinicians to better study cardiovascular conditions using shear wave elastography. However, these methods rely ona prioriknowledge of the vibrational modes dominating the propagating waves induced by acoustic radiation force excitation: differences between anticipated and real modal content are known to yield errors in the inversion. We seek to improve the accuracy of this process by modeling the artery as a fluid-immersed cylindrical waveguide and building an analytical framework to prescribe radiation force excitations that will selectively excite certain waveguide modes using ultrasound acoustic radiation force. We show that all even-numbered waveguide modes can be eliminated from the arterial response to perturbation, and confirm the efficacy of this approach within silicotests that show that odd modes are preferentially excited. Finally, by analyzing data from phantom tests, we find a set of ultrasound focal parameters that demonstrate the viability of inducing the desired odd-mode response in experiments.

Multimodal guided wave inversion for arterial stiffness: methodology and validation in phantoms

Arterial stiffness is an important biomarker for many cardiovascular diseases. Shear wave elastography is a recent technique aimed at estimating local arterial stiffness using guided wave inversion (GWI), i.e. matching the computed and measured wave dispersion. This paper develops and validates a new GWI approach by synthesizing various recent observations and algorithms: (a) refinements to signal processing to obtain more accurate experimental dispersion curves; (b) an efficient forward model to compute theoretical dispersion curves for immersed, incompressible cylindrical waveguides; (c) an optimization framework based on the recent observation that the measured dispersion curve is multimodal, i.e. it matches for not one but two different wave modes in two different frequency ranges. The resulting inversion approach is validated using extensive experimental data from rubber tube phantoms, not only for modulus estimation but also to simultaneously estimate modulus and wall thickness. The observations indicate that the modulus estimates are best performed with the information on wall thickness. The approach, which takes less than half a minute to run, is shown to be accurate, with the modulus estimated with less than 4% error for 70% of the experiments.

Hydrophone Spatial Averaging Correction for Acoustic Exposure Measurements From Arrays-Part I: Theory and Impact on Diagnostic Safety Indexes

This article reports underestimation of mechanical index (MI) and nonscanned thermal index for bone near focus (TIB) due to hydrophone spatial averaging effects that occur during acoustic output measurements for clinical linear and phased arrays. TIB is the appropriate version of thermal index (TI) for fetal imaging after ten weeks from the last menstrual period according to the American Institute of Ultrasound in Medicine (AIUM). Spatial averaging is particularly troublesome for highly focused beams and nonlinear, nonscanned modes such as acoustic radiation force impulse (ARFI) and pulsed Doppler. MI and variants of TI (e.g., TIB), which are displayed in real-time during imaging, are often not corrected for hydrophone spatial averaging because a standardized method for doing so does not exist for linear and phased arrays. A novel analytic inverse-filter method to correct for spatial averaging for pressure waves from linear and phased arrays is derived in this article (Part I) and experimentally validated in a companion article (Part II). A simulation was developed to estimate potential spatial-averaging errors for typical clinical ultrasound imaging systems based on the theoretical inverse filter and specifications for 124 scanner/transducer combinations from the U.S. Food and Drug Administration (FDA) 510(k) database from 2015 to 2019. Specifications included center frequency, aperture size, acoustic output parameters, hydrophone geometrical sensitive element diameter, etc. Correction for hydrophone spatial averaging using the inverse filter suggests that maximally achievable values for MI, TIB, thermal dose ( t ₄₃ ), and spatial-peak-temporal-average intensity ( [Formula: see text]) for typical clinical systems are potentially higher than uncorrected values by (means ± standard deviations) 9% ± 4% (ARFI MI), 19% ± 15% (ARFI TIB), 50% ± 41% (ARFI t ₄₃ ), 43% ± 39% (ARFI [Formula: see text]), 7% ± 5% (pulsed Doppler MI), 15% ± 11% (pulsed Doppler TIB), 42% ± 31% (pulsed Doppler t ₄₃ ), and 33% ± 27% (pulsed Doppler [Formula: see text]). These values correspond to frequencies of 3.2 ± 1.3 (ARFI) and 4.1 ± 1.4 MHz (pulsed Doppler), and the model predicts that they would increase with frequency. Inverse filtering for hydrophone spatial averaging significantly improves the accuracy of estimates of MI, TIB, t ₄₃ , and [Formula: see text] for ARFI and pulsed Doppler signals.

High Frequency Ultrasound Elastography for Estimating the Viscoelastic Properties of the Cornea Using Lamb Wave Model

OBJECTIVE

Estimating the elasticity distribution in the cornea is important because corneal elasticity is usually influenced by corneal pathologies and surgical treatments, especially for early corneal sclerosis. Because the thickness of the cornea is typically less than 1 mm, high-resolution ultrasound elastography as well as the Lamb wave model is required for viscoelastic property estimation. In the present study, an array high-frequency ultrasound (HFUS) elastography method based on ultrafast ultrasound imaging was proposed for estimating the viscoelastic properties of porcine cornea.

METHODS

The elastic wave was generated by an external vibrator, after which the wave propagation image was obtained using a 40-MHz array transducer. Viscoelasticity estimation was performed by fitting the phase velocity curve using the Lamb wave model. The performance of the proposed HFUS elastography system was verified using 2-mm-thick thin-layer gelatin phantoms with gelatin concentrations of 7% and 12%. Ex vivo experiments were carried out using fresh porcine cornea with artificial sclerosing.

RESULTS

Experimental results showed that the estimated elasticity was close to the standard value obtained in the phantom study when the Lamb wave model was used for elasticity measurement. However, the error between the standard elasticity values and the elasticity values estimated using group shear wave velocity was large. In the ex vivo eyeball experiments, the estimated elasticities and viscosities were respectively 9.1 ± 1.3 kPa and 0.5 ± 0.10 Pa·s for a healthy cornea and respectively 15.9 ± 2.1 kPa and 1.1 ± 0.12 Pa·s for a cornea with artificial sclerosis. A 3D HFUS elastography was also obtained for distinguishing the region of sclerosis in the cornea.

CONCLUSION

The experimental results demonstrated that the proposed HFUS elastography method has high potential for the clinical diagnosis of corneal diseases compared with other HFUS single-element transducer elastography systems.

Vascular Shear Wave Elastography in Atherosclerotic Arteries: A Systematic Review

Ischemic stroke is a leading cause of death and disability worldwide, so adequate prevention strategies are crucial. However, current stroke risk stratification is based on epidemiologic studies and is still suboptimal for individual patients. The aim of this systematic review was to provide a literature overview on the feasibility and diagnostic value of vascular shear wave elastography (SWE) using ultrasound (US) in (mimicked) human and non-human arteries affected by different stages of atherosclerotic diseases or diseases related to atherosclerosis. An online search was conducted on Pubmed, Embase, Web of Science and IEEE databases to identify studies using US SWE for the assessment of vascular elasticity. A quality assessment was performed using Quality Assessment of Diagnostic Accuracy Studies 2 (QUADAS-2) checklist, and relevant data were extracted. A total of 19 studies were included: 10 with human patients and 9 with non-human subjects (i.e., [excised] animal arteries and polyvinyl alcohol phantoms). All studies revealed the feasibility of using US SWE to assess individually stiffness of the arterial wall and plaques. Quantitative elasticity values were highly variable between studies. However, within studies, SWE could detect statistically significant elasticity differences in patient/subject characteristics and could distinguish different plaque types with good reproducibility. US SWE, with its unique ability to assess the elasticity of the vessel wall and plaque throughout the cardiac cycle, might be a good candidate to improve stroke risk stratification. However, more clinical studies have to be performed to assess this technique's exact clinical value.

Combined spatiotemporal and frequency-dependent shear wave elastography enables detection of vulnerable carotid plaques as validated by MRI

Fatal cerebrovascular events are often caused by rupture of atherosclerotic plaques. However, rupture-prone plaques are often distinguished by their internal composition rather than degree of luminal narrowing, and conventional imaging techniques might thus fail to detect such culprit lesions. In this feasibility study, we investigate the potential of ultrasound shear wave elastography (SWE) to detect vulnerable carotid plaques, evaluating group velocity and frequency-dependent phase velocities as novel biomarkers for plaque vulnerability. In total, 27 carotid plaques from 20 patients were scanned by ultrasound SWE and magnetic resonance imaging (MRI). SWE output was quantified as group velocity and frequency-dependent phase velocities, respectively, with results correlated to intraplaque constituents identified by MRI. Overall, vulnerable lesions graded as American Heart Association (AHA) type VI showed significantly higher group and phase velocity compared to any other AHA type. A selection of correlations with intraplaque components could also be identified with group and phase velocity (lipid-rich necrotic core content, fibrous cap structure, intraplaque hemorrhage), complementing the clinical lesion classification. In conclusion, we demonstrate the ability to detect vulnerable carotid plaques using combined SWE, with group velocity and frequency-dependent phase velocity providing potentially complementary information on plaque characteristics. With such, the method represents a promising non-invasive approach for refined atherosclerotic risk prediction.

In vivo evaluation of posterior eye elasticity using shaker-based optical coherence elastography

Age-related macular degeneration (AMD) is a progressive retinal disease and becomes the leading cause of blindness. It is well established that early detection is the key to preservation of functional vision. However, it is very difficult to diagnose AMD in very early stages, before structural changes are evident. Consequently, investigating the biomechanical properties of the retina maybe essential for understanding its physiological function. In this study, we present a shear wave-based quantitative method for estimating the elasticity of the posterior eye using shaker-based optical coherence elastography. This technique has been developed and validated on both a homogeneous phantom and a healthy rabbit in vivo. The shear wave speed from the ganglion side to the photoreceptor side of the rabbit eye is 4.1 m/s, 4.9 m/s, and 6.7 m/s, respectively. In addition, the most stiff sclera region has an average shear wave speed of 9.1 m/s. The results demonstrate the feasibility of using this technique to quantify biomechanical properties of the posterior eye and its potential translation to the clinical study.

Impact statement

Herein, we propose a potentially clinical applicable shaker-based optical coherence elastography (OCE) technique to characterize the biomechanical properties of the posterior eye, including different layers of the retina. Compared with either acoustic radiation force OCE or air-puff OCE, the newly developed method can induce sufficient shear wave propagation at the posterior eye with high resolution and large field of view.

Lagrangian carotid strain imaging indices normalized to blood pressure for vulnerable plaque

OBJECTIVE

Ultrasound Lagrangian carotid strain imaging (LCSI) utilizes physiological deformation caused by arterial pressure variations to generate strain tensor maps of the vessel walls and plaques. LCSI has been criticized for the lack of normalization of magnitude-based strain indices to physiological stimuli, namely blood pressure. We evaluated the impact of normalization of magnitude-based strain indices to blood pressure measured immediately after the acquisition of radiofrequency (RF) data loops for LCSI.

MATERIALS AND METHODS

A complete clinical ultrasound examination along with RF data loops for LCSI was performed on 50 patients (30 males and 20 females) who presented with >60% carotid stenosis and were scheduled for carotid endarterectomy. Cognition was assessed using the 60-minute neuropsychological test protocol.

RESULTS

For axial strains correlation of maximum accumulated strain indices (MASI), cognition scores were -0.46 for non-normalized and -0.45, -0.49, -0.37, and -0.48 for systolic, diastolic, pulse pressure, and mean arterial pressure normalized data, respectively. The corresponding area under the curve (AUC) values for classifiers designed using maximum likelihood estimation of a binormal distribution with a median-split of the executive function cognition scores were 0.73, 0.70, 0.71, 0.70, and 0.71, respectively.

CONCLUSIONS

No significant differences in the AUC estimates were obtained between normalized and non-normalized magnitude-based strain indices.

Quantitative Time-Harmonic Ultrasound Elastography of the Abdominal Aorta and Inferior Vena Cava

The purpose of this study was to evaluate the sensitivity of quantitative time-harmonic ultrasound elastography (THE) of the inferior vena cava (IVC) and abdominal aorta (AA) to changes in central volume status. THE of the IVC and AA was performed in 20 healthy volunteers before and after oral intake of 1 L of water and before or during passive leg raising to augment venous filling. Compound maps of shear wave speed (SWS) as surrogate measures of vessel wall stiffness were generated within the full field of view from multifrequency harmonic wave fields. SWS was measured in regions of the IVC and AA. Blood pressure, stroke volume, cardiac output and pulse wave velocity were recorded. Statistical significance of SWS changes was tested using one-way repeated-measures analysis of variance. SWS measured in the IVC increased from 1.71 ± 0.1 m/s before water intake to 1.82 ± 0.1 m/s during passive leg raising and, further, to 1.87 ± 0.1 m/s after hydration and to 1.95 ± 0.1 m/s with hydration plus passive leg raising (p < 0.001). SWS in the AA did not change significantly after hydration (2.14 ± 0.13 m/s vs. 2.15 ± 0.16 m/s; p = 0.792). SWS was significantly higher in the AA than in the IVC across all experiments (p < 0.001). Water drinking did not significantly influence blood pressure, pulse wave velocity and cardiac output (all p values >0.1), whereas stroke volume increased significantly (p = 0.031). Time-harmonic ultrasound elastography enables quantification of the wall stiffness of the large abdominal vessels and is sensitive to different volume and pressure states in the IVC.

Bidirectional Ultrasound Elastographic Imaging Framework for Non-invasive Assessment of the Non-linear Behavior of a Physiologically Pressurized Artery

Studies of non-destructive bidirectional ultrasound assessment of non-linear mechanical behavior of the artery are scarce in the literature. We hereby propose derivation of a strain-shear modulus relationship as a new graphical diagnostic index using an ultrasound elastographic imaging framework, which encompasses our in-house bidirectional vascular guided wave imaging (VGWI) and ultrasound strain imaging (USI). This framework is used to assess arterial non-linearity in two orthogonal (i.e., longitudinal and circumferential) directions in the absence of non-invasive pressure measurement. Bidirectional VGWI estimates longitudinal (μ_L) and transverse (μ_T) shear moduli, whereas USI estimates radial strain (ɛ_r). Vessel-mimicking phantoms (with and without longitudinal pre-stretch) and in vitro porcine aortas under static and/or dynamic physiologic intraluminal pressure loads were examined. ɛ_r was found to be a suitable alternative to intraluminal pressure for representation of cyclic loading on the artery wall. Results revealed that μ_T values of all samples examined increased non-linearly with ε_r magnitude and more drastically than μ_L, whereas μ_L values of only the pre-stretched phantoms and aortas increased with ɛ_r magnitude. As a new graphical representation of arterial non-linearity and function, strain-shear modulus loops derived by the proposed framework over two consecutive dynamic loading cycles differentiated sample pre-conditions and corroborated direction-dependent non-linear mechanical behaviors of the aorta with high estimation repeatability.

Arterial Stiffness Assessment by Shear Wave Elastography and Ultrafast Pulse Wave Imaging: Comparison with Reference Techniques in Normotensives and Hypertensives

Shear wave elastography and ultrafast imaging of the carotid artery pulse wave were performed in 27 normotensive participants and 29 age- and sex-matched patients with essential hypertension, and compared with reference techniques: carotid-femoral pulse wave velocity (cfPWV) determined via arterial tonometry and carotid stiffness (carPWV) determined via echotracking. Shear wave speed in the carotid anterior (a-SWS) and posterior (p-SWS) walls were assessed throughout the cardiac cycle. Ultrafast PWV was measured in early systole (ufPWV-FW) and in end-systole (dicrotic notch, ufPWV-DN). Shear wave speed in the carotid anterior appeared to be the best candidate to evaluate arterial stiffness from ultrafast imaging. In univariate analysis, a-SWS was associated with carPWV (r = 0.56, p = 0.003) and carotid-to-femoral PWV (r = 0.66, p < 0.001). In multivariate analysis, a-SWS was independently associated with age (R² = 0.14, p = 0.02) and blood pressure (R² = 0.21, p = 0.004). Moreover, a-SWS increased with blood pressure throughout the cardiac cycle and did not differ between normotensive participants and patients with essential hypertension when compared at similar blood pressures.

Experimental Investigation of Guided Wave Imaging in Thin Soft Media under Various Coupling Conditions

Guided wave imaging for the artery remains in its infancy in clinical practice mainly because of complex arterial microstructure, hemodynamics and boundary conditions. Despite the theoretically known potential effect of the surrounding medium on guided wave propagation in thin media in non-destructive testing, experimental evidence pertaining to thin soft materials, such as the artery, is relatively scarce in the relevant literature. Therefore, this study first evaluated the propagating guided wave generated by acoustic radiation force in polyvinyl alcohol-based hydrogel plates differing in thickness and stiffness under various material coupling conditions (water and polyvinyl alcohol bulk). A thin-walled polyvinyl alcohol hollow cylindrical phantom coupled by softer gelatin-agar phantoms and an excised porcine aorta surrounded by water and pork belly were further examined. Guided waves in the thin structure and shear waves in the bulk media were captured by ultrafast ultrasound imaging, and guided wave dispersion as a function of the frequency-thickness product was analyzed using the zero-order anti-symmetric Lamb wave model to estimate the shear modulus of each thin medium studied. Results confirmed the deviated shear modulus estimates from the ground truth for thin plates, the thin-walled hollow cylindrical phantom and the porcine aorta bounded by stiffness-unmatched bulk medium. The findings indicated the need for (i) careful interpretation of estimated shear moduli of thin structure bounded by bulk media and (ii) a generalized guided wave model that takes into account the effect of coupling medium.

Quantitative Assessment of Thin-Layer Tissue Viscoelastic Properties Using Ultrasonic Micro-Elastography With Lamb Wave Model

Characterizing the viscoelastic properties of thin-layer tissues with micro-level thickness has long remained challenging. Recently, several micro-elastography techniques have been developed to improve the spatial resolution. However, most of these techniques have not considered the medium boundary conditions when evaluating the viscoelastic properties of thin-layer tissues such as arteries and corneas; this might lead to estimation bias or errors. This paper aims to integrate the Lamb wave model with our previously developed ultrasonic micro-elastography imaging system for obtaining accurate viscoelastic properties in thin-layer tissues. A 4.5-MHz ring transducer was used to generate an acoustic radiation force for inducing tissue displacements to produce guided wave, and the wave propagation was detected using a confocally aligned 40-MHz needle transducer. The phase velocity and attenuation were obtained from k-space by both the impulse and the harmonic methods. The measured phase velocity was fit using the Lamb wave model with the Kelvin-Voigt model. Phantom experiments were conducted using 7% and 12% gelatin and 1.5% agar phantoms with different thicknesses (2, 3, and 4 mm). Biological experiments were performed on porcine cornea and rabbit carotid artery ex vivo. Thin-layer phantoms with different thicknesses were confirmed to have the same elasticity; this was consistent with the estimates of bulk phantoms from mechanical tests and the shear wave rheological model. The trend of the measured attenuations was also confirmed with the viscosity results obtained using the Lamb wave model. Through the impulse and harmonic methods, the shear viscoelasticity values were estimated to be 8.2 kPa for $0.9~\text {Pa}{\cdot} \text {s}$ and 9.6 kPa for $0.8~\text {Pa}{\cdot} \text {s}$ in the cornea and 27.9 kPa for $0.1~\text {Pa}\cdot \text {s}$ and 26.5 kPa for $0.1~\text {Pa}\cdot \text {s}$ in the artery.

Ultrasonic Microelastography to Assess Biomechanical Properties of the Cornea

OBJECTIVE

To both qualitatively and quantitatively investigate corneal biomechanical properties through an ultrasonic microelastography imaging system, which is potentially useful in the diagnosis of diseases, such as keratoconus, postrefractive keratectasia, and tracking treatment such as cross-linking surgery.

METHODS

Our imaging system has a dual-frequency configuration, including a 4.5 MHz ring transducer to push the tissue and a confocally aligned 40 MHz needle transducer to track micron-level displacement. Two-dimensional/three-dimensional acoustic radiation force impulse (ARFI) imaging and Young's modulus in the region of interest were performed on ex vivo porcine corneas that were either cross-linked using formalin solution or preloaded with intraocular pressure (IOPs) from 5 to 30 mmHg.

RESULTS

The increase of corneal stiffness and the change in cross-linked volume following formalin crosslinking could be precisely observed in the ARFI images and reflected by the reconstructed Young's modulus while the B-mode structural images remained almost unchanged. In addition, the relationship between the stiffness of the cornea and IOPs was investigated among 12 porcine corneas. The corneal stiffness is significantly different at various IOPs and has a tendency to become stiffer with increasing IOP.

CONCLUSION

Our results demonstrate the principle of using ultrasonic microelastography techniques to image the biomechanical properties of the cornea. Integrating high-resolution ARFI imaging labeled with reconstructed Young's modulus and structural imaging of the cornea can potentially lead to a routinely performed imaging modality in the field of ophthalmology.

Fourier-Domain Shift Matching: A Robust Time-of-Flight Approach for Shear Wave Speed Estimation

Our primary objective of this work was to design and test a new time-of-flight (TOF) method that allows measurements of shear wave speed (SWS) following impulsive excitation in soft tissues. Particularly, under the assumption of the local plane shear wave, this work named the Fourier-domain shift matching (FDSM) method, estimates SWS by aligning a series of shear waveforms either temporally or spatially using a solution space deduced by characteristic curves of the well-known 1-D wave equation. The proposed SWS estimation method was tested using computer-simulated data, and tissue-mimicking phantom and ex vivo tissue experiments. Its performance was then compared with three other known TOF methods: lateral time-to-peak (TTP) method with robust random sampling consensus (RANSAC) fitting method, Radon sum transformation method, and a modified cross correlation method. Hereafter, these three TOF methods are referred to as the TTP-RANSAC, Radon sum, and X-corr methods, respectively. In addition to an adapted form of the 2-D Fourier transform (2-D FT)-based method in which the (group) SWS was approximated by averaging phase SWS values was considered for comparison. Based on data evaluated, we found that the overall performance of the above-mentioned temporal implementation of the proposed FDSM method was most similar to the established Radon sum method (correlation = 0.99, scale factor = 1.03, and mean difference = 0.07 m/s), and the 2-D FT (correlation = 0.98, scale factor = 1.00, and mean difference = 0.10 m/s) at high signal quality. However, results obtained from the 2-D FT method diverged (correlation = 0.201) from these of the proposed temporal implementation in the presence of diminished signal quality, whereas the agreement between the Radon sum approach and the proposed temporal implementation largely remained the same (correlation = 0.98).

Multidirectional Estimation of Arterial Stiffness Using Vascular Guided Wave Imaging with Geometry Correction

We previously found that vascular guided wave imaging (VGWI) could non-invasively quantify transmural wall stiffness in both the longitudinal (r-z plane, 0°) and circumferential (r-θ plane, 90°) directions of soft hollow cylinders. Arterial stiffness estimation in multiple directions warrants further comprehensive characterization of arterial health, especially in the presence of asymmetric plaques, but is currently lacking. This study therefore investigated the multidirectional estimation of the arterial Young's modulus in a finite-element model, in vitro artery-mimicking phantoms and an excised porcine aorta. A longitudinal pre-stretch of 20% and/or lumen pressure (15 or 70 mm Hg) was additionally introduced to pre-condition the phantoms for emulating the intrinsic mechanical anisotropy of the real artery. The guided wave propagation was approximated by a zero-order antisymmetric Lamb wave model. Shape factor, which was defined as the ratio of inner radius to thickness, was calculated over the entire segment of each planar cross section of the hollow cylindrical structure at a full rotation (0°-360° at 10° increments) about the radial axis. The view-dependent geometry of the cross segment was found to affect the guided wave propagation, causing Young's modulus overestimation in four angular intervals along the propagation pathway, all of which corresponded to wall regions with low shape factors (<1.5). As validated by mechanical tensile testing, the results indicate not only that excluding the propagation pathway with low shape factors could correct the overestimation of Young's modulus, but also that VGWI could portray the anisotropy of hollow cylindrical structures and the porcine aorta based on the derived fractional anisotropy values from multidirectional modulus estimates. This study may serve as an important step toward 3-D assessment of the mechanical properties of the artery.

Quantification of carotid artery plaque stability with multiple region of interest based ultrasound strain indices and relationship with cognition

Vulnerability and instability in carotid artery plaque has been assessed based on strain variations using noninvasive ultrasound imaging. We previously demonstrated that carotid plaques with higher strain indices in a region of interest (ROI) correlated to patients with lower cognition, probably due to cerebrovascular emboli arising from these unstable plaques. This work attempts to characterize the strain distribution throughout the entire plaque region instead of being restricted to a single localized ROI. Multiple ROIs are selected within the entire plaque region, based on thresholds determined by the maximum and average strains in the entire plaque, enabling generation of additional relevant strain indices. Ultrasound strain imaging of carotid plaques, was performed on 60 human patients using an 18L6 transducer coupled to a Siemens Acuson S2000 system to acquire radiofrequency data over several cardiac cycles. Patients also underwent a battery of neuropsychological tests under a protocol based on National Institute of Neurological Disorders and Stroke and Canadian Stroke Network guidelines. Correlation of strain indices with composite cognitive index of executive function revealed a negative association relating high strain to poor cognition. Patients grouped into high and low cognition groups were then classified using these additional strain indices. One of our newer indices, namely the average L  -  1 norm with plaque (AL1NWP) presented with significantly improved correlation with executive function when compared to our previously reported maximum accumulated strain indices. An optimal combination of three of the new indices generated classifiers of patient cognition with an area under the curve (AUC) of 0.880, 0.921 and 0.905 for all (n  =  60), symptomatic (n  =  33) and asymptomatic patients (n  =  27) whereas classifiers using maximum accumulated strain indices alone provided AUC values of 0.817, 0.815 and 0.813 respectively.