Quantifying effect of intraplaque hemorrhage on critical plaque wall stress in human atherosclerotic plaques using three-dimensional fluid-structure interaction models

Load-sharing characteristics of stenting and post-dilation in heavily calcified coronary artery

In this work, stenting in non-calcified and heavily calcified coronary arteries was quantified in terms of diameter-pressure relationships and load transfer from the balloon to the artery. The efficacy of post-dilation in non-calcified and heavily calcified coronary arteries was also characterized in terms of load sharing and the changes in tissue mechanics. Our results have shown that stent expansion exhibits a cylindrical shape in non-calcified lesions, while it exhibits a dog bone shape in heavily calcified lesions. Load-sharing analysis has shown that only a small portion of the pressure load (1.4 N, 0.8% of total pressure load) was transferred to the non-calcified lesion, while a large amount of the pressure load (19 N, 12%) was transferred to the heavily calcified lesion. In addition, the increasing inflation pressure (from 10 to 20 atm) can effectively increase the minimal lumen diameter (from 1.48 to 2.82 mm) of the heavily calcified lesion, the stress (from 1.5 to 8.4 MPa) and the strain energy in the calcification (1.77 mJ to 26.5 mJ), which are associated with the potential of calcification fracture. Results indicated that increasing inflation pressure can be an effective way to improve the stent expansion if a dog bone shape of the stenting profile is observed. Considering the risk of a balloon burst, our results support the design and application of the high-pressure balloon for post-dilation. This work also sheds some light on the stent design and choice of stent materials for improving the stent expansion at the dog bone region and mitigating stresses on arterial tissues.

High-risk features of basilar artery atherosclerotic plaque

INTRODUCTION

High-resolution magnetic resonance imaging (HR-MRI) is used to characterize atherosclerotic plaque. The present study aimed to determine the high-risk features of the basilar artery (BA) atherosclerotic plaque.

METHODS

Patients with advanced BA stenosis were screened. The features including the ruptured fibrous cap (RFC), lipid core, intraplaque hemorrhage (IPH), plaque enhancement, and calcification were assessed by using high-resolution MRI. The relationship between the features and acute infarction was analyzed.

RESULTS

From 1 June 2014 to 31 December 2018, a total of 143 patients with 76 new strokes were included. RFC was identified in 25% of symptomatic and 10.4% of asymptomatic patients. IPH was identified in 48.7% of symptomatic and 25.4% of asymptomatic patients. RFC (3.157, 95% CI 1.062 to 9.382, p = 0.039) and IPH (2.78, 95% CI 1.127 to 6.505, p = 0.026) were independent risk factors for acute infarction.

CONCLUSION

Our study showed that RFC and IPH of BA plaque were independent risk factors for acute infarction.

Increased Proximal Wall Shear Stress of Basilar Artery Plaques Associated with Ruptured Fibrous Cap

Plaque rupture of the basilar artery is one of the leading causes of posterior circulation stroke. The present study aimed to investigate the role of fluid dynamics in the ruptured fibrous cap of basilar artery plaques. Patients with basilar artery plaques (50−99% stenosis) were screened. Integrity of the fibrous cap was assessed by high-resolution MRI. Computational fluid dynamics models were built based on MR angiography to obtain the wall shear stress and velocity. A total of 176 patients were included. High-resolution MRI identified 35 ruptured fibrous caps of basilar artery plaques. Ruptured fibrous cap was significantly associated with acute infarction (27/35 vs. 96/141, p < 0.05) in the territory of the basilar artery. Proximal wall shear stress of stenosis was positively related with the ruptured fibrous cap (OR 1.564; 95% CI, 1.101−2.222; p = 0.013). The threshold of wall shear stress for the ruptured fibrous cap of basilar artery plaques was 4.84 Pa (Area under ROC 0.732, p = 0.008, 95%CI 0.565−0.899). The present study demonstrated that increased proximal wall shear stress of stenosis was associated with ruptured fibrous caps of basilar artery plaques.

Studying the Factors of Human Carotid Atherosclerotic Plaque Rupture, by Calculating Stress/Strain in the Plaque, Based on CEUS Images: A Numerical Study

Carotid plaque neovascularization is one of the major factors for the classification of vulnerable plaque, but the axial force effects of the pulsatile blood flow on the plaque with neovessel and intraplaque hemorrhage was unclear. Together with the severity of stenosis, the fibrous cap thickness, large lipid core, and the neovascularization followed by intraplaque hemorrhage (IPH) have been regarded as high-risk features of plaque rupture. In this work, the effects of these factors were evaluated on the progression and rupture of the carotid atherosclerotic plaques. Five geometries of carotid artery plaque were developed based on contrast-enhanced ultrasound (CEUS) images, which contain two types of neovessel and IPH, and geometry without neovessel and IPH. A one-way fluid-structure interaction model was applied to compute the maximum principal stress and strain in the plaque. For that hyper-elastic and non-linear material, Yeoh 3rd Order strain energy density function was used for components of the plaque. The simulation results indicated that the maximum principal stress of plaque in the carotid artery was higher when the degree of the luminal stenosis increased and the thickness of the fibrous cap decreased. The neovessels within the plaque could introduce a 2.5% increments of deformation in the plaque under the pulsatile blood flow pressure. The IPH also contributed to the increased risk of plaque rupture that a gain of stress was 8.983, 14.526, and 34.47 kPa for the plaque with 50, 65, and 75%, respectively, when comparing stress in the plaque with IPH distributed at the middle to the shoulder of the plaque. In conclusion, neovascularization in the plaque could reduce the stability of the plaque by increasing the stress within the plaque. Also, the risk of plaque rupture increased when large luminal stenosis, thin fibrous cap, and IPH were observed.

Normalized wall index, intraplaque hemorrhage and ulceration of carotid plaques correlate with the severity of ischemic stroke

BACKGROUND AND AIMS

Carotid atherosclerosis is considered an important cause of ischemic stroke. Tthis study aimed to explore the relationship between plaque features and the severity of stroke, and to identify plaque risk factors for the assessment of the severity of ischemic stroke.

METHODS

Symptomatic patients with carotid atherosclerotic plaques were prospectively recruited and underwent high-resolution vessel wall magnetic resonance imaging (VW-MRI). Two trained MRI readers independently identified intraplaque hemorrhage (IPH), calcification (CA), surface CA, deep CA, and ulceration. They measured and calculated the maximum vessel diameter (Max VD), maximum wall thickness (Max WT), total vessel area, lumen area, normalized wall index (NWI), plaque volume, IPH volume, IPH proportion, CA volume, and CA proportion. Patients were divided into two groups according to their National Institutes of Health Stroke Scale (NIHSS) scores (NIHSS ≤1 vs. NIHSS >1). Clinical characteristics and carotid plaque features were compared using the Mann-Whitney U test or Chi-square test as appropriate. Odds ratio (OR) and corresponding 95% confidence interval (CI) of plaque features to distinguish patients with NIHSS >1 were calculated. Spearman's rank correlations or Pearson correlations were determined for plaque features and NIHSS scores.

RESULTS

Of the 97 included patients, 34 (35.05%) with NIHSS >1 had significantly greater NWI (p < 0.05), larger IPH volume (p < 0.01), and greater IPH proportion (p < 0.01), and higher prevalence of IPH (OR, 5.654; 95%CI, 2.272-14.070; p < 0.01) and ulceration (OR, 2.891; 95%CI, 1.090-7.667; p = 0.033) than patients with NIHSS ≤1. Max WT (r = 0.24, p = 0.018), NWI (r = 0.22, p = 0.032), IPH (r = 0.27, p = 0.007), IPH volume (r = 0.35, p < 0.01), IPH proportion (r = 0.28, p = 0.005), and ulceration (r = 0.35, p < 0.01) had positive correlations with NIHSS scores.

CONCLUSIONS

NWI, IPH, and ulceration of carotid atherosclerotic plaque based on high-resolution VW-MRI may be useful indicators for assessing the severity of ischemic stroke in patients with atherosclerosis. NIHSS score is related to max WT, NWI, IPH, IPH volume, IPH proportion, and ulceration.

Combining morphological and biomechanical factors for optimal carotid plaque progression prediction: An MRI-based follow-up study using 3D thin-layer models

Plaque progression prediction is of fundamental significance to cardiovascular research and disease diagnosis, prevention, and treatment. Magnetic resonance image (MRI) data of carotid atherosclerotic plaques were acquired from 20 patients with consent obtained. 3D thin-layer models were constructed to calculate plaque stress and strain. Data for ten morphological and biomechanical risk factors were extracted for analysis. Wall thickness increase (WTI), plaque burden increase (PBI) and plaque area increase (PAI) were chosen as three measures for plaque progression. Generalized linear mixed models (GLMM) with 5-fold cross-validation strategy were used to calculate prediction accuracy and identify optimal predictor. The optimal predictor for PBI was the combination of lumen area (LA), plaque area (PA), lipid percent (LP), wall thickness (WT), maximum plaque wall stress (MPWS) and maximum plaque wall strain (MPWSn) with prediction accuracy = 1.4146 (area under the receiver operating characteristic curve (AUC) value is 0.7158), while PA, plaque burden (PB), WT, LP, minimum cap thickness, MPWS and MPWSn was the best for WTI (accuracy = 1.3140, AUC = 0.6552), and a combination of PA, PB, WT, MPWS, MPWSn and average plaque wall strain (APWSn) was the best for PAI with prediction accuracy = 1.3025 (AUC = 0.6657). The combinational predictors improved prediction accuracy by 9.95%, 4.01% and 1.96% over the best single predictors for PAI, PBI and WTI (AUC values improved by 9.78%, 9.45%, and 2.14%), respectively. This suggests that combining both morphological and biomechanical risk factors could lead to better patient screening strategies.

Association Between Carotid Atherosclerotic Plaque Calcification and Intraplaque Hemorrhage

Objective—

Carotid intraplaque hemorrhage (IPH) is associated with cardiovascular events. Calcification, which frequently accompanies IPH, may play a role in IPH occurrence. In this study, we aimed to investigate the associations between calcification characteristics and IPH in carotid plaques.

Approach and results—

One hundred seventeen patients with cerebrovascular symptoms and carotid plaques detected by ultrasound were recruited and underwent multicontrast magnetic resonance imaging. Advanced carotid plaques with composition measured by magnetic resonance imaging were included in the analysis. Carotid calcifications were divided into the following categories: surface, mixed, and deep calcification. They were also classified into single and multiple calcifications according to quantity. Logistic regression models utilizing generalized estimating equations were performed to evaluate the relationship between calcification and IPH. Of 117 subjects, 85 with 142 plaques were included in the final analysis, whereas 32 were excluded because of lack of plaque compositions. Of the 142 plaques, 40 (28.2%) had IPH. Plaques with IPH showed greater prevalence of calcification than those without (87.5% versus 55.9%; P =0.005). After adjusting for age, low-density lipoprotein, maximum wall thickness, and maximum soft plaque thickness, multiple calcifications (odd ratio, 10.1; 95% confidence interval, 3.3–30.4), surface calcification (odd ratio, 29.4; 95% confidence interval, 4.1–210.8), and mixed calcifications (odd ratio, 27.9; 95% confidence interval, 7.3–107.1) were found to be strongly associated with the presence of IPH (all P <0.05).

Conclusions—

Surface calcification and multiple calcifications in carotid atherosclerotic plaques are independently associated with the presence of IPH, suggesting that both quantity and location of calcification may play important roles in the occurrence of IPH. These findings may provide novel insights for understanding mechanisms of IPH.

3D MRI-based multicomponent thin layer structure only plaque models for atherosclerotic plaques

MRI-based fluid-structure interactions (FSI) models for atherosclerotic plaques have been developed to perform mechanical analysis to investigate the association of plaque wall stress (PWS) with cardiovascular disease. However, the time consuming 3D FSI model construction process is a great hinder for its clinical implementations. In this study, a 3D thin-layer structure only (TLS) plaque model was proposed as an approximation with much less computational cost to 3D FSI models for better clinical implementation potential. 192 TLS models were constructed based on 192 ex vivo MRI Images of 12 human coronary atherosclerotic plaques. Plaque stresses were extracted from all lumen nodal points. The maximum value of Plaque wall stress (MPWS) and average value of plaque wall stress (APWS) of each slice were used to compare with those from corresponding FSI models. The relative errors for MPWS and APWS were 9.76% and 9.89%, respectively. Both MPWS and APWS values obtained from TLS models showed very good correlation with those from 3D FSI models. Correlation results from TLS models were in consistent with FSI models. Our results indicated that the proposed 3D TLS plaque models may be used as a good approximation to 3D FSI models with much less computational cost. With further validation, 3D TLS models may be possibly used to replace FSI models to save time and perform mechanical analysis for atherosclerotic plaques for clinical implementation.

Quest for the Vulnerable Atheroma: Carotid Stenosis and Diametric Strain--A Feasibility Study

The Bernoulli effect may result in eruption of a vulnerable carotid atheroma, causing a stroke. We measured electrocardiography (ECG)-registered QRS intra-stenotic blood velocity and atheroma strain dynamics in carotid artery walls using ultrasonic tissue Doppler methods, providing displacement and time resolutions of 0.1 μm and 3.7 ms. Of 22 arteries, 1 had a peak systolic velocity (PSV) >280 cm/s, 4 had PSVs between 165 and 280 cm/s and 17 had PSVs <165 cm/s. Eight arteries with PSVs <65 cm/s and 4 of 9 with PSVs between 65 and 165 cm/s had normal systolic diametric expansion (0% and 7%) and corresponding systolic wall thinning. The remaining 10 arteries had abnormal systolic strain dynamics, 2 with diametric reduction (>-0.05 mm), 2 with extreme wall expansion (>0.1 mm), 2 with extreme wall thinning (>-0.1 mm) and 4 with combinations. Decreases in systolic diameter and/or extreme systolic arterial wall thickening may indicate imminent atheroma rupture.

Piecewise Pulse Wave Imaging (pPWI) for Detection and Monitoring of Focal Vascular Disease in Murine Aortas and Carotids In Vivo

Atherosclerosis and Abdominal Aortic Aneurysms (AAAs) are two common vascular diseases associated with mechanical changes in the arterial wall. Pulse Wave Imaging (PWI), a technique developed by our group to assess and quantify the mechanical properties of the aortic wall in vivo, may provide valuable diagnostic information. This work implements piecewise PWI (pPWI), an enhanced version of PWI designed for focal vascular diseases. Localized, sub-regional PWVs and PWI moduli ( EPWI ) were estimated within 2-4 mm wall segments of murine normal, atherosclerotic and aneurysmal arteries. Overall, stiffness was found to increase in the atherosclerotic cases. The mean sub-regional PWV was found to be 2.57±0.18 m/s for the normal aortas (n = 7) with a corresponding mean EPWI of 43.82±5.86 kPa. A significant increase ( (p ≤ 0.001)) in the group means of the sub-regional PWVs was found between the normal aortas and the aortas of mice on high-fat diet for 20 ( 3.30±0.36 m/s) and 30 weeks ( 3.56±0.29 m/s). The mean of the sub-regional PWVs ( 1.57±0.78 m/s) and EPWI values ( 19.23±15.47 kPa) decreased significantly in the aneurysmal aortas (p ≤ 0.05) . Furthermore, the mean coefficient of determination (r(2)) of the normal aortas was significantly higher (p ≤ 0.05) than those of the aneurysmal and atherosclerotic cases. These findings demonstrated that pPWI may be able to provide useful biomarkers for monitoring focal vascular diseases.

Subject-Specific Fully-Coupled and One-Way Fluid-Structure Interaction Models for Modeling of Carotid Atherosclerotic Plaques in Humans

Background

Hemodynamics play an important role in the development and progression of carotid atherosclerosis, and may be important in the assessment of plaque vulnerability. The aim of this study was to develop a system to assess the hemodynamics of carotid atherosclerotic plaques using subject-specific fluid-structure interaction (FSI) models based on magnetic resonance imaging (MRI).

Material/Methods

Models of carotid bifurcations (n=86 with plaques from 52 patients, n=14 normal carotids from 12 participants) were obtained at the Department of Radiology, Beijing Tian Tan Hospital between 2010 and 2013. The maximum von Mises stress, minimum pressure, and flow velocity values were assessed at the most stenotic site in patients, or at the carotid bifurcations in healthy volunteers. Results of one-way FSI were compared with fully-coupled FSI for the plaques of 19 randomly selected models.

Results

The maximum von Mises stress and the minimum pressure and velocity were significantly increased in the stenosis group compared with controls based on one-way FSI (all P<0.05). The maximum von Mises stress and the minimum pressure were significantly higher and the velocity was significantly lower based on fully coupled FSI compared with on-way FSI (all P<0.05). Although there were differences in numerical values, both methods were equivalent. The maximum von Mises stress of vulnerable plaques was significantly higher than stable plaques (P<0.001). The maximum von Mises stress of the group with fibrous cap defect was significantly higher than the group without fibrous cap defect (P=0.001).

Conclusions

The hemodynamics of atherosclerotic plaques can be assessed noninvasively using subject-specific models of FSI based on MRI.

Haemodynamical stress in mouse aortic arch with atherosclerotic plaques: Preliminary study of plaque progression

Atherosclerotic plaques develop at particular sites in the arterial tree, and this regional localisation depends largely on haemodynamic parameters (such as wall shear stress; WSS) as described in the literature. Plaque rupture can result in heart attack or stroke and hence understanding the development and vulnerability of atherosclerotic plaques is critically important. The purpose of this study is to characterise the haemodynamics of blood flow in the mouse aortic arch using numerical modelling. The geometries are digitalised from synchrotron imaging and realistic pulsatile blood flow is considered under rigid wall assumptions. Two cases are considered; arteries with and without plaque. Mice that are fed under fat diet present plaques in the aortic arch whose size is dependent on the number of weeks under the diet. The plaque distribution in the region is however relatively constant through the different samples. This result underlines the influence of the geometry and consequently of the wall shear stresses for plaque formation with plaques growing in region of relative low shear stresses. A discussion of the flow field in real geometry in the presence and absence of plaques is conducted. The presence of plaques was shown to alter the blood flow and hence WSS distribution, with regions of localised high WSS, mainly on the wall of the brachiocephalic artery where luminal narrowing is most pronounced. In addition, arch plaques are shown to induce recirculation in the blood flow, a phenomenon with potential influence on the progression of the plaques. The oscillatory shear index and the relative residence time have been calculated on the geometry with plaques to show the presence of this recirculation in the arch, an approach that may be useful for future studies on plaque progression.

Image-based modeling for better understanding and assessment of atherosclerotic plaque progression and vulnerability: data, modeling, validation, uncertainty and predictions

Medical imaging and image-based modeling have made considerable progress in recent years in identifying atherosclerotic plaque morphological and mechanical risk factors which may be used in developing improved patient screening strategies. However, a clear understanding is needed about what we have achieved and what is really needed to translate research to actual clinical practices and bring benefits to public health. Lack of in vivo data and clinical events to serve as gold standard to validate model predictions is a severe limitation. While this perspective paper provides a review of the key steps and findings of our group in image-based models for human carotid and coronary plaques and a limited review of related work by other groups, we also focus on grand challenges and uncertainties facing the researchers in the field to develop more accurate and predictive patient screening tools.

Plaque hemorrhage in carotid artery disease: pathogenesis, clinical and biomechanical considerations

Stroke remains the most prevalent disabling illness today, with internal carotid artery luminal stenosis due to atheroma formation responsible for the majority of ischemic cerebrovascular events. Severity of luminal stenosis continues to dictate both patient risk stratification and the likelihood of surgical intervention. But there is growing evidence to suggest that plaque morphology may help improve pre-existing risk stratification criteria. Plaque components such a fibrous tissue, lipid rich necrotic core and calcium have been well investigated but plaque hemorrhage (PH) has been somewhat overlooked. In this review we discuss the pathogenesis of PH, its role in dictating plaque vulnerability, PH imaging techniques, marterial properties of atherosclerotic tissues, in particular, those obtained based on in vivo measurements and effect of PH in modulating local biomechanics.

Higher critical plaque wall stress in patients who died of coronary artery disease compared with those who died of other causes: a 3D FSI study based on ex vivo MRI of coronary plaques

Mechanical forces play an important role in the rupture of vulnerable plaques. This process is often associated with cardiovascular syndromes, such as heart attack and stroke. In this study, magnetic resonance imaging (MRI)-based models were used to investigate the association between plaque wall stress (PWS) and coronary artery disease (CAD). Ex vivo MRI data of coronary plaques from 12 patients were used to construct 12 three-dimensional (3D) fluid-structure interaction (FSI) computational models. Six of the patients had died from CAD and six had died from non-CAD causes. PWS was assessed using all nodal points on the lumen surface of each plaque. The maximum PWS from all possible vulnerable sites of each plaque was defined as the 3D critical plaque wall stress (CPWS). Mean 3D CPWS in the CAD group was 94.3% higher than that in the non-CAD group (265.6 vs. 136.7 kPa, P=0.0029). There was no statistically significant difference in global maximum plaque wall stress (GMPWS) between the two groups (P=0.347). There was also no statistically significant difference in plaque burden between the CAD group (84.4±5%) and the non-CAD group (82.0±8%, P=0.552). The results indicate that plaques from patients who died from CAD were associated with higher CPWS compared with those from patients who died from non-CAD causes. With further validation, analysis of CPWS may prove to be an important component in assessment of plaque vulnerability.

Fatigue Crack Propagation Analysis of Plaque Rupture

Rupture of atheromatous plaque is the major cause of stroke or heart attack. Considering that the cardiovascular system is a classic fatigue environment, plaque rupture was treated as a chronic fatigue crack growth process in this study. Fracture mechanics theory was introduced to describe the stress status at the crack tip and Paris' law was used to calculate the crack growth rate. The effect of anatomical variation of an idealized plaque cross-section model was investigated. The crack initiation was considered to be either at the maximum circumferential stress location or at any other possible locations around the lumen. Although the crack automatically initialized at the maximum circumferential stress location usually propagated faster than others, it was not necessarily the most critical location where the fatigue life reached its minimum. We found that the fatigue life was minimum for cracks initialized in the following three regions: the midcap zone, the shoulder zone, and the backside zone. The anatomical variation has a significant influence on the fatigue life. Either a decrease in cap thickness or an increase in lipid pool size resulted in a significant decrease in fatigue life. Comparing to the previously used stress analysis, this fatigue model provides some possible explanations of plaque rupture at a low stress level in a pulsatile cardiovascular environment, and the method proposed here may be useful for further investigation of the mechanism of plaque rupture based on in vivo patient data.