Non-Invasive Identification of Vulnerable Atherosclerotic Plaques Using Texture Analysis in Ultrasound Carotid Elastography: An In Vivo Feasibility Study Validated by Magnetic Resonance Imaging

Novel Imaging-Based Biomarkers for Identifying Carotid Plaque Vulnerability

Carotid artery disease has traditionally been assessed based on the degree of luminal narrowing. However, this approach, which solely relies on carotid stenosis, is currently being questioned with regard to modern risk stratification approaches. Recent guidelines have introduced the concept of the "vulnerable plaque," emphasizing specific features such as thin fibrous caps, large lipid cores, intraplaque hemorrhage, plaque rupture, macrophage infiltration, and neovascularization. In this context, imaging-based biomarkers have emerged as valuable tools for identifying higher-risk patients. Non-invasive imaging modalities and intravascular techniques, including ultrasound, computed tomography, magnetic resonance imaging, intravascular ultrasound, optical coherence tomography, and near-infrared spectroscopy, have played pivotal roles in characterizing and detecting unstable carotid plaques. The aim of this review is to provide an overview of the evolving understanding of carotid artery disease and highlight the significance of imaging techniques in assessing plaque vulnerability and informing clinical decision-making.

Novel ultrasound techniques in the identification of vulnerable plaques-an updated review of the literature

Atherosclerosis is an inflammatory disease partly mediated by lipoproteins. The rupture of vulnerable atherosclerotic plaques and thrombosis are major contributors to the development of acute cardiovascular events. Despite various advances in the treatment of atherosclerosis, there has been no satisfaction in the prevention and assessment of atherosclerotic vascular disease. The identification and classification of vulnerable plaques at an early stage as well as research of new treatments remain a challenge and the ultimate goal in the management of atherosclerosis and cardiovascular disease. The specific morphological features of vulnerable plaques, including intraplaque hemorrhage, large lipid necrotic cores, thin fibrous caps, inflammation, and neovascularisation, make it possible to identify and characterize plaques with a variety of invasive and non-invasive imaging techniques. Notably, the development of novel ultrasound techniques has introduced the traditional assessment of plaque echogenicity and luminal stenosis to a deeper assessment of plaque composition and the molecular field. This review will discuss the advantages and limitations of five currently available ultrasound imaging modalities for assessing plaque vulnerability, based on the biological characteristics of the vulnerable plaque, and their value in terms of clinical diagnosis, prognosis, and treatment efficacy assessment.

Interframe Echo Intensity Variation of Subregions and Whole Plaque in Two-Dimensional Carotid Ultrasonography: Simulations and In Vivo Observations

OBJECTIVES

The risk of cardiovascular disease is associated with the echo intensity of carotid plaques in ultrasound images and their cardiac cycle-induced intensity variations. In this study, we aimed to 1) explore the underlying origin of echo intensity variations by using simulations and 2) evaluate the association between the two-dimensional (2D) spatial distribution of these echo intensity variations and plaque vulnerability.

METHODS

First, we analyzed how out-of-plane motion and compression of simulated scattering spheres of different sizes affect the ultrasound echo intensity. Next, we propose a method to analyze the features of the 2D spatial distribution of interframe plaque echo intensity in carotid ultrasound image sequences and explore their associations with plaque vulnerability in experimental data.

RESULTS

The simulations showed that the magnitude of echo intensity changes was similar for both the out-of-plane motion and compression, but for scattering objects smaller than 1 mm radius, the out-of-plane motion dominated. In experimental data, maps of the 2D spatial distribution of the echo intensity variations had a low correlation with standard B-mode echo intensity distribution, indicating complementary information on plaque tissue composition. In addition, we found the existence of ∼1 mm diameter subregions with pronounced echo intensity variations associated with plaque vulnerability.

CONCLUSIONS

The results indicate that out-of-plane motion contributes to intra-plaque regions of high echo intensity variation. The 2D echo intensity variation maps may provide complementary information for assessing plaque composition and vulnerability. Further studies are needed to verify this method's role in identifying vulnerable plaques and predicting cardiovascular disease risk.

Imaging and Hemodynamic Characteristics of Vulnerable Carotid Plaques and Artificial Intelligence Applications in Plaque Classification and Segmentation

Stroke is a massive public health problem. The rupture of vulnerable carotid atherosclerotic plaques is the most common cause of acute ischemic stroke (AIS) across the world. Currently, vessel wall high-resolution magnetic resonance imaging (VW-HRMRI) is the most appropriate and cost-effective imaging technique to characterize carotid plaque vulnerability and plays an important role in promoting early diagnosis and guiding aggressive clinical therapy to reduce the risk of plaque rupture and AIS. In recent years, great progress has been made in imaging research on vulnerable carotid plaques. This review summarizes developments in the imaging and hemodynamic characteristics of vulnerable carotid plaques on the basis of VW-HRMRI and four-dimensional (4D) flow MRI, and it discusses the relationship between these characteristics and ischemic stroke. In addition, the applications of artificial intelligence in plaque classification and segmentation are reviewed.

Advances in Noninvasive Carotid Wall Imaging with Ultrasound: A Narrative Review

Carotid atherosclerosis is a major cause for stroke, with significant associated disease burden morbidity and mortality in Western societies. Diagnosis, grading and follow-up of carotid atherosclerotic disease relies on imaging, specifically ultrasound (US) as the initial modality of choice. Traditionally, the degree of carotid lumen stenosis was considered the sole risk factor to predict brain ischemia. However, modern research has shown that a variety of other imaging biomarkers, such as plaque echogenicity, surface morphology, intraplaque neovascularization and vasa vasorum contribute to the risk for rupture of carotid atheromas with subsequent cerebrovascular events. Furthermore, the majority of embolic strokes of undetermined origin are probably arteriogenic and are associated with nonstenosing atheromas. Therefore, a state-of-the-art US scan of the carotid arteries should take advantage of recent technical developments and should provide detailed information about potential thrombogenic (/) and emboligenic arterial wall features. This manuscript reviews recent advances in ultrasonographic assessment of vulnerable carotid atherosclerotic plaques and highlights the fields of future development in multiparametric arterial wall imaging, in an attempt to convey the most important take-home messages for clinicians performing carotid ultrasound.

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.

Combining First- and Second-Order Continuity Constraints in Ultrasound Elastography

Ultrasound elastography is a prominent noninvasive medical imaging technique that estimates tissue elastic properties to detect abnormalities in an organ. A common approximation to tissue elastic modulus is tissue strain induced after mechanical stimulation. To compute tissue strain, ultrasound radio frequency (RF) data can be processed using energy-based algorithms. These algorithms suffer from ill-posedness to tackle. A continuity constraint along with the data amplitude similarity is imposed to obtain a unique solution to the time-delay estimation (TDE) problem. Existing energy-based methods exploit the first-order spatial derivative of the displacement field to construct a regularizer. This first-order regularization scheme alone is not fully consistent with the mechanics of tissue deformation while perturbed with an external force. As a consequence, state-of-the-art techniques suffer from two crucial drawbacks. First, the strain map is not sufficiently smooth in uniform tissue regions. Second, the edges of the hard or soft inclusions are not well-defined in the image. Herein, we address these issues by formulating a novel regularizer taking both first- and second-order derivatives of the displacement field into account. The second-order constraint, which is the principal novelty of this work, contributes both to background continuity and edge sharpness by suppressing spurious noisy edges and enhancing strong boundaries. We name the proposed technique: Second-Order Ultrasound eLastography (SOUL). Comparative assessment of qualitative and quantitative results shows that SOUL substantially outperforms three recently developed TDE algorithms called Hybrid, GLUE, and MPWC-Net++. SOUL yields 27.72%, 62.56%, and 81.37% improvements of the signal-to-noise ratio (SNR) and 72.35%, 54.03%, and 65.17% improvements of the contrast-to-noise ratio (CNR) over GLUE with data pertaining to simulation, phantom, and in vivo tissue, respectively. The SOUL code can be downloaded from code.sonography.ai.

Choosing the right therapy for a patient with asymptomatic carotid stenosis

Introduction: Most patients with asymptomatic carotid stenosis (ACS) now have a lower risk with intensive medical therapy than with stenting (CAS) or endarterectomy (CEA); the annual risk of stroke or death with intensive medical therapy is ~ 0.5%, vs. a periprocedural risk with CAS of ~ 2.5-4.1% with CAS, and ~ 1.4-1.8% with CEA. The excess risk of CAS is greater in older patients.Areas covered: Discussed are the need for intensive medical therapy, the nature of intensive medical therapy, approaches to identifying the few patients with ACS who could benefit from CEA or CAS, and which patients would be better suited to CEA vs. CAS.Expert opinion: All patients with ACS are at high risk of cardiovascular events, soshould receive intensive medical therapy including lifestyle modification, intensive lipid-lowering, B vitamins to lower homocysteine (using methylcobalamin rather than cyanocobalamin), and appropriate antithrombotic therapy. High-risk patients who could benefit from intervention can be identified by clinical and imaging features including transcranial Doppler embolus detection, ulceration, intraplaque hemorrhage, reduced cerebrovascular reserve, plaque echolucency, silent infarction on brain imaging, and progression of stenosis. Most patients whose risk of stroke warrants intervention would be better treated with CEA than with CAS.

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.

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 assessment of the intracranial vasculature in an older adult population using iCafe

Comprehensive quantification of intracranial artery features may help us assess and understand variations of blood supply during brain development and aging. We analyzed vasculature features of 163 participants (age 56-85 years, mean of 71) from a community study to investigate if any of the features varied with age. Three-dimensional time-of-flight magnetic resonance angiography images of these participants were processed in IntraCranial artery feature extraction technique (a recently developed technique to obtain quantitative features of arteries) to divide intracranial vasculatures into anatomical segments and generate 8 morphometry and intensity features for each segment. Overall, increase in age was found negatively associated with number of branches and average order of intracranial arteries while positively associated with tortuosity, which remained after adjusting for cardiovascular risk factors. The associations with number of branches and average order were consistently found between 3 main intracranial artery regions, whereas the association with tortuosity appeared to be present only in middle cerebral artery/distal arteries. The combination of time-of-flight magnetic resonance angiography and IntraCranial artery feature extraction technique may provide an effective way to study vascular conditions and changes in the aging brain.

Diagnostic value of two dimensional shear wave elastography combined with texture analysis in early liver fibrosis

BACKGROUND

Staging diagnosis of liver fibrosis is a prerequisite for timely diagnosis and therapy in patients with chronic hepatitis B. In recent years, ultrasound elastography has become an important method for clinical noninvasive assessment of liver fibrosis stage, but its diagnostic value for early liver fibrosis still needs to be further improved. In this study, the texture analysis was carried out on the basis of two dimensional shear wave elastography (2D-SWE), and the feasibility of 2D-SWE plus texture analysis in the diagnosis of early liver fibrosis was discussed.

AIM

To assess the diagnostic value of 2D-SWE combined with textural analysis in liver fibrosis staging.

METHODS

This study recruited 46 patients with chronic hepatitis B. Patients underwent 2D-SWE and texture analysis; Young's modulus values and textural patterns were obtained, respectively. Textural pattern was analyzed with regard to contrast, correlation, angular second moment (ASM), and homogeneity. Pathological results of biopsy specimens were the gold standard; comparison and assessment of the diagnosis efficiency were conducted for 2D-SWE, texture analysis and their combination.

RESULTS

2D-SWE displayed diagnosis efficiency in early fibrosis, significant fibrosis, severe fibrosis, and early cirrhosis (AUC > 0.7, P < 0.05) with respective AUC values of 0.823 (0.678-0.921), 0.808 (0.662-0.911), 0.920 (0.798-0.980), and 0.855 (0.716-0.943). Contrast and homogeneity displayed independent diagnosis efficiency in liver fibrosis stage (AUC > 0.7, P < 0.05), whereas correlation and ASM showed limited values. AUC of contrast and homogeneity were respectively 0.906 (0.779-0.973), 0.835 (0.693-0.930), 0.807 (0.660-0.910) and 0.925 (0.805-0.983), 0.789 (0.639-0.897), 0.736 (0.582-0.858), 0.705 (0.549-0.883) and 0.798 (0.650-0.904) in four liver fibrosis stages, which exhibited equivalence to 2D-SWE in diagnostic efficiency (P > 0.05). Combined diagnosis (PRE) displayed diagnostic efficiency (AUC > 0.7, P < 0.01) for all fibrosis stages with respective AUC of 0.952 (0.841-0.994), 0.896 (0.766-0.967), 0.978 (0.881-0.999), 0.947 (0.835-0.992). The combined diagnosis showed higher diagnosis efficiency over 2D-SWE in early liver fibrosis (P < 0.05), whereas no significant differences were observed in other comparisons (P > 0.05).

CONCLUSION

Texture analysis was capable of diagnosing liver fibrosis stage, combined diagnosis had obvious advantages in early liver fibrosis, liver fibrosis stage might be related to the hepatic tissue hardness distribution.

In vivo carotid strain imaging using principal strains in longitudinal view

Carotid plaque rupture can result in stroke or transient ischemic attack that can be devastating for patients. Ultrasound strain imaging provides a noninvasive method to identify unstable plaque likely to rupture. Axial, lateral and shear strains in carotid plaque have been shown to be linked to carotid plaque instability. Recently, there has been interest in using principal strains, which do not depend on angle of insonification of the carotid artery for quantifying instability in plaque along the longitudinal view. In this work relationships between angle dependent axial, lateral and shear strain along with axis independent principal strains are compared. Three strain indices were defined, 1) Average Mean Strain (AMS), 2) Maximum Mean Strain (MMS) and 3) Mean Standard Deviation (MSD) to identify relationships between these five strain image types in a group of 76 in vivo patients. The maximum principal strain demonstrated the highest strain values when compared to axial strain for all patients with a linear regression slope of 1.6 and a y intercept of 2.4 percent strain for AMS. The maximum shear strain when compared to shear strain had a slope of 1.15 and a y intercept of 0.21 percent for AMS. Next, the effect of insonification angle, which is the angle subtended by the artery at the location of plaque was studied. Patients were divided into three sub groups, i.e. less than 5 degrees (n = 31), between 5 and 10 degrees (n = 24) and above 10 degrees (n = 21). The angle of insonification did not make a significant difference between the three angle groups when comparing the relationship between the angle dependent and independent strain values.

Assessment of Carotid Artery Plaque Components With Machine Learning Classification Using Homodyned-K Parametric Maps and Elastograms

Quantitative ultrasound (QUS) imaging methods, including elastography, echogenicity analysis, and speckle statistical modeling, are available from a single ultrasound (US) radio-frequency data acquisition. Since these US imaging methods provide complementary quantitative tissue information, characterization of carotid artery plaques may gain from their combination. Sixty-six patients with symptomatic ( n = 26 ) and asymptomatic ( n = 40 ) carotid atherosclerotic plaques were included in the study. Of these, 31 underwent magnetic resonance imaging (MRI) to characterize plaque vulnerability and quantify plaque components. US radio-frequency data sequence acquisitions were performed on all patients and were used to compute noninvasive vascular US elastography and other QUS features. Additional QUS features were computed from three types of images: homodyned-K (HK) parametric maps, Nakagami parametric maps, and log-compressed B-mode images. The following six classification tasks were performed: detection of 1) a small area of lipid; 2) a large area of lipid; 3) a large area of calcification; 4) the presence of a ruptured fibrous cap; 5) differentiation of MRI-based classification of nonvulnerable carotid plaques from neovascularized or vulnerable ones; and 6) confirmation of symptomatic versus asymptomatic patients. Feature selection was first applied to reduce the number of QUS parameters to a maximum of three per classification task. A random forest machine learning algorithm was then used to perform classifications. Areas under receiver-operating curves (AUCs) were computed with a bootstrap method. For all tasks, statistically significant higher AUCs were achieved with features based on elastography, HK parametric maps, and B-mode gray levels, when compared to elastography alone or other QUS alone ( ). For detection of a large area of lipid, the combination yielding the highest AUC (0.90, 95% CI 0.80-0.92, ) was based on elastography, HK, and B-mode gray-level features. To detect a large area of calcification, the highest AUC (0.95, 95% CI 0.94-0.96, ) was based on HK and B-mode gray level features. For other tasks, AUCs varied between 0.79 and 0.97. None of the best combinations contained Nakagami features. This study shows the added value of combining different features computed from a single US acquisition with machine learning to characterize carotid artery plaques.

Interoperator Reproducibility of Carotid Elastography for Identification of Vulnerable Atherosclerotic Plaques

Ultrasound-based carotid elastography has been developed to evaluate the vulnerability of carotid atherosclerotic plaques. The aim of this study was to investigate the in vivo interoperator reproducibility of carotid elastography for the identification of vulnerable plaques, with high-resolution magnetic resonance imaging (MRI) as reference. Ultrasound radio-frequency data of 45 carotid arteries (including 53 plaques) from 32 volunteers were acquired separately by two experienced operators in the longitudinal view and then were used to estimate the interframe axial strain rate (ASR) with a two-step optical flow method. The maximum 99th percentile of absolute ASR of all plaques in a carotid artery was used as the elastographic index. MRI scanning was also performed on each volunteer to identify the vulnerable plaque. The results showed no systematic bias in the Bland-Altman plot and an intraclass correlation coefficient of 0.66 between the two operators. In addition, no statistical significance was found between the receiver operating characteristic (ROC) curves from the two operators ( ), and their areas under the ROC curves were 0.83 and 0.77, respectively. Using the mean measurements of the two operators as the classification criterion, a sensitivity of 71.4%, a specificity of 87.1%, and an accuracy of 82.2% were obtained with a cutoff value of 1.37 [Formula: see text]. This study validates the interoperator reproducibility of ultrasound-based carotid elastography for identifying vulnerable carotid plaques.