The influence of axial image resolution on atherosclerotic plaque stress computations

Fully automated construction of three-dimensional finite element simulations from Optical Coherence Tomography

Despite recent advances in diagnosis and treatment, atherosclerotic coronary artery diseases remain a leading cause of death worldwide. Various imaging modalities and metrics can detect lesions and predict patients at risk; however, identifying unstable lesions is still difficult. Current techniques cannot fully capture the complex morphology-modulated mechanical responses that affect plaque stability, leading to catastrophic failure and mute the benefit of device and drug interventions. Finite Element (FE) simulations utilizing intravascular imaging OCT (Optical Coherence Tomography) are effective in defining physiological stress distributions. However, creating 3D FE simulations of coronary arteries from OCT images is challenging to fully automate given OCT frame sparsity, limited material contrast, and restricted penetration depth. To address such limitations, we developed an algorithmic approach to automatically produce 3D FE-ready digital twins from labeled OCT images. The 3D models are anatomically faithful and recapitulate mechanically relevant tissue lesion components, automatically producing morphologies structurally similar to manually constructed models whilst including more minute details. A mesh convergence study highlighted the ability to reach stress and strain convergence with average errors of just 5.9% and 1.6% respectively in comparison to FE models with approximately twice the number of elements in areas of refinement. Such an automated procedure will enable analysis of large clinical cohorts at a previously unattainable scale and opens the possibility for in-silico methods for patient specific diagnoses and treatment planning for coronary artery disease.

Automated finite element approach to generate anatomical patient-specific biomechanical models of atherosclerotic arteries from virtual histology-intravascular ultrasound

Despite advancements in early detection and treatment, atherosclerosis remains the leading cause of death across all cardiovascular diseases (CVD). Biomechanical analysis of atherosclerotic lesions has the potential to reveal biomechanically instable or rupture-prone regions. Treatment decisions rarely consider the biomechanics of the stenosed lesion due in-part to difficulties in obtaining this information in a clinical setting. Previous 3D FEA approaches have incompletely incorporated the complex curvature of arterial geometry, material heterogeneity, and use of patient-specific data. To address these limitations and clinical need, herein we present a user-friendly fully automated program to reconstruct and simulate the wall mechanics of patient-specific atherosclerotic coronary arteries. The program enables 3D reconstruction from patient-specific data with heterogenous tissue assignment and complex arterial curvature. Eleven arteries with coronary artery disease (CAD) underwent baseline and 6-month follow-up angiographic and virtual histology-intravascular ultrasound (VH-IVUS) imaging. VH-IVUS images were processed to remove background noise, extract VH plaque material data, and luminal and outer contours. Angiography data was used to orient the artery profiles along the 3D centerlines. The resulting surface mesh is then resampled for uniformity and tetrahedralized to generate the volumetric mesh using TetGen. A mesh convergence study revealed edge lengths between 0.04 mm and 0.2 mm produced constituent volumes that were largely unchanged, hence, to save computational resources, a value of 0.2 mm was used throughout. Materials are assigned and finite element analysis (FEA) is then performed to determine stresses and strains across the artery wall. In a representative artery, the highest average effective stress was in calcium elements with 235 kPa while necrotic elements had the lowest average stress, reaching as low as 0.79 kPa. After applying nodal smoothening, the maximum effective stress across 11 arteries remained below 288 kPa, implying biomechanically stable plaques. Indeed, all atherosclerotic plaques remained unruptured at the 6-month longitudinal follow up diagnosis. These results suggest our automated analysis may facilitate assessment of atherosclerotic plaque stability.

A platform for high-fidelity patient-specific structural modelling of atherosclerotic arteries: from intravascular imaging to three-dimensional stress distributions

The pathophysiology of atherosclerotic lesions, including plaque rupture triggered by mechanical failure of the vessel wall, depends directly on the plaque morphology-modulated mechanical response. The complex interplay between lesion morphology and structural behaviour can be studied with high-fidelity computational modelling. However, construction of three-dimensional (3D) and heterogeneous models is challenging, with most previous work focusing on two-dimensional geometries or on single-material lesion compositions. Addressing these limitations, we here present a semi-automatic computational platform, leveraging clinical optical coherence tomography images to effectively reconstruct a 3D patient-specific multi-material model of atherosclerotic plaques, for which the mechanical response is obtained by structural finite-element simulations. To demonstrate the importance of including multi-material plaque components when recovering the mechanical response, a computational case study was conducted in which systematic variation of the intraplaque lipid and calcium was performed. The study demonstrated that the inclusion of various tissue components greatly affected the lesion mechanical response, illustrating the importance of multi-material formulations. This platform accordingly provides a viable foundation for studying how plaque micro-morphology affects plaque mechanical response, allowing for patient-specific assessments and extension into clinically relevant patient cohorts.

Stress distribution analysis in healthy and stenosed carotid artery models reconstructed from in vivo ultrasonography

Purpose

This study investigated the accuracy of models reconstructed from ultrasound image processing by comparing the radial displacement waveforms of a subject-specific artery model and evaluated stress changes in the proximal shoulder, throat, and distal shoulder of the plaques depending on the degree of carotid artery stenosis.

Methods

Three groups of subjects (healthy and with less than 50% or more carotid stenosis) were evaluated with ultrasonography. Two-dimensional transverse imaging of the common carotid artery was performed to reconstruct the geometry. A longitudinal view of the same region was recorded to extract the Kelvin viscoelastic model parameters. The pulse pressure waveform and the effective pressure of perivascular tissue were loaded onto the internal and external walls of the model. Effective, circumferential, and principal stresses applied to the plaque throat, proximal shoulder, and distal shoulder in the transverse planes were extracted.

Results

The radial displacement waveforms of the model were closely correlated with those of image processing in all three groups. The mean of the effective, circumferential, and principal stresses of the healthy arteries were 15.01±4.93, 12.97±5.07, and 12.39±2.86 kPa, respectively. As stenosis increased from mild to significant, the mean values of the effective, circumferential, and first principal stresses increased significantly (97%, 74%, and 103% at the plaque throat, respectively) (P<0.05). The minimum effective stress was at the lipid pool. The effective stress in calcified areas was higher than in other parts of the artery wall.

Conclusion

This model can discriminate differences in stresses applied to mildly and severely stenotic plaques.

Consideration of stiffness of wall layers is decisive for patient-specific analysis of carotid artery with atheroma

The paper deals with the impact of chosen geometric and material factors on maximal stresses in carotid atherosclerotic plaque calculated using patient-specific finite element models. These stresses are believed to be decisive for the plaque vulnerability but all applied models suffer from inaccuracy of input data, especially when obtained in vivo only. One hundred computational models based on ex vivo MRI are used to investigate the impact of wall thickness, MRI slice thickness, lipid core and fibrous tissue stiffness, and media anisotropy on the calculated peak plaque and peak cap stresses. The investigated factors are taken as continuous in the range based on published experimental results, only the impact of anisotropy is evaluated by comparison with a corresponding isotropic model. Design of Experiment concept is applied to assess the statistical significance of these investigated factors representing uncertainties in the input data of the model. The results show that consideration of realistic properties of arterial wall in the model is decisive for the stress evaluation; assignment of properties of fibrous tissue even to media and adventitia layers as done in some studies may induce up to eightfold overestimation of peak stress. The impact of MRI slice thickness may play a key role when local thin fibrous cap is present. Anisotropy of media layer is insignificant, and the stiffness of fibrous tissue and lipid core may become significant in some combinations.

Automated 3D geometry segmentation of the healthy and diseased carotid artery in free-hand, probe tracked ultrasound images

PURPOSE

Rupture of an arterosclerotic plaque in the carotid artery is a major cause of stroke. Biomechanical analysis of plaques is under development aiming to aid the clinician in the assessment of plaque vulnerability. Patient-specific three-dimensional (3D) geometry assessment of the carotid artery, including the bifurcation, is required as input for these biomechanical models. This requires a high-resolution, 3D, noninvasive imaging modality such as ultrasound (US). In this study, a high-resolution two-dimensional (2D) linear array in combination with a magnetic probe tracking device and automatic segmentation method was used to assess the geometry of the carotid artery. The advantages of using this system over a 3D ultrasound probe are its higher resolution (spatial and temporal) and its larger field of view.

METHODS

A slow sweep (v = ± 5 mm/s) was made over the subject's neck so that the full geometry of the bifurcated geometry of the carotid artery is captured. An automated segmentation pipeline was developed. First, the Star-Kalman method was used to approximate the center and size of the vessels for every frame. Images were filtered with a Gaussian high-pass filter before conversion into the 2D monogenic signals, and multiscale asymmetry features were extracted from these data, enhancing low lateral wall-lumen contrast. These images, in combination with the initial ellipse contours, were used for an active deformable contour model to segment the vessel lumen. To segment the lumen-plaque boundary, Otsu's automatic thresholding method was used. Distension of the wall due to the change in blood pressure was removed using a filter approach. Finally, the contours were converted into a 3D hexahedral mesh for a patient-specific solid mechanics model of the complete arterial wall.

RESULTS

The method was tested on 19 healthy volunteers and on 3 patients. The results were compared to manual segmentation performed by three experienced observers. Results showed an average Hausdorff distance of 0.86 mm and an average similarity index of 0.91 for the common carotid artery (CCA) and 0.88 for the internal and external carotid artery. For the total algorithm, the success rate was 89%, in 4 out of 38 datasets the ICA and ECA were not sufficient visible in the US images. Accurate 3D hexahedral meshes were successfully generated from the segmented images .

CONCLUSIONS

With this method, a subject-specific biomechanical model can be constructed directly from a hand-held 2D US measurement, within 10 min, with a minimal user input. The performance of the proposed segmentation algorithm is comparable to or better than algorithms previously described in literature. Moreover, the algorithm is able to segment the CCA, ICA, and ECA including the carotid bifurcation in transverse B-mode images in both healthy and diseased arteries.

Moderate thickness of lipid core in shoulder region of atherosclerotic plaque determines vulnerable plaque A parametric study

Peak stress in the fibrous cap of atherosclerotic plaque is largely determined by the cap thickness which cannot be accurately estimated in vivo. This parametric study investigates idealized atherosclerotic plaque geometries. Finite element modeling is applied to search for larger morphological features associated with high cap stresses. By varying seven geometrical and two loading parameters, 100 3D model geometries of atherosclerotic plaques in common iliac artery were generated. In each model peak cap stress was calculated, and statistical comparison of the geometries generating the highest and lowest peak cap stresses was performed. The analysis showed that, compared to geometries generating the lowest stresses, those with high peak cap stress had a significantly lower cap thickness, higher stenosis ratio, lower relative lipid core volume, and cap shoulder radius larger than lipid core radius. High cap stress was observed for cap thicknesses up to 0.13 mm. It can be concluded that vulnerable plaques contain thin fibrous cap, large stenosis ratio and only moderate small-radius lipid core which reaches the shoulder region of the fibrous cap.

Effect of distal thickening and stiffening of plaque cap on arterial wall mechanics

To investigate the effect of longitudinal variations of cap thickness and tissue properties on wall stresses and strains along the atherosclerotic stenosis, stenotic plaque models (uniformly thick, distally thickened, homogenous, and distally stiffened) were constructed and subjected to computational stress analyses with due consideration of fluid-structure interactions (FSI). The analysis considered three different cap thicknesses-45, 65, and 200 μm-and tissue properties-soft, fibrous, and hard. The maximum peak cap stress (PCS) and strain were observed in the upstream throat section and demonstrated increases of the order of 345 and 190%, respectively, as the cap thickness was reduced from 200 to 45 μm in uniformly thick models. Distal stiffening increased PCS in the downstream region; however, the overall effect of this increase was rather small. Distal thickening did not affect maximum PCS and strain values for cap thicknesses exceeding 65 μm; however, a noticeable increase in maximum PCS and corresponding longitudinal variation (or spatial gradient) in stress was observed in the very thin (45-μm-thick) cap. It was, therefore, inferred that existence of a rather thin upstream cap demonstrating distal cap thickening indicates an increased risk of plaque progression and rupture. Graphical Abstract ᅟ.

Intima heterogeneity in stress assessment of atherosclerotic plaques

Atherosclerotic plaque rupture is recognized as the primary cause of cardiac and cerebral ischaemic events. High structural plaque stresses have been shown to strongly correlate with plaque rupture. Plaque stresses can be computed with finite-element (FE) models. Current FE models employ homogeneous material properties for the heterogeneous atherosclerotic intima. This study aimed to evaluate the influence of intima heterogeneity on plaque stress computations. Two-dimensional FE models with homogeneous and heterogeneous intima were constructed from histological images of atherosclerotic human coronaries (n = 12). For homogeneous models, a single stiffness value was employed for the entire intima. For heterogeneous models, the intima was subdivided into four clusters based on the histological information and different stiffness values were assigned to the clusters. To cover the reported local intima stiffness range, 100 cluster stiffness combinations were simulated. Peak cap stresses (PCSs) from the homogeneous and heterogeneous models were analysed and compared. By using a global variance-based sensitivity analysis, the influence of the cluster stiffnesses on the PCS variation in the heterogeneous intima models was determined. Per plaque, the median PCS values of the heterogeneous models ranged from 27 to 160 kPa, and the PCS range varied between 43 and 218 kPa. On average, the homogeneous model PCS values differed from the median PCS values of heterogeneous models by 14%. A positive correlation (R² = 0.72) was found between the homogeneous model PCS and the PCS range of the heterogeneous models. Sensitivity analysis showed that the highest main sensitivity index per plaque ranged from 0.26 to 0.83, and the average was 0.47. Intima heterogeneity resulted in substantial changes in PCS, warranting stress analyses with heterogeneous intima properties for plaque-specific, high accuracy stress assessment. Yet, computations with homogeneous intima assumption are still valuable to perform sensitivity analyses or parametric studies for testing the effect of plaque geometry on PCS. Moreover, homogeneous intima models can help identify low PCS, stable type plaques with thick caps. Yet, for thin cap plaques, accurate stiffness measurements of the clusters in the cap and stress analysis with heterogeneous cap properties are required to characterize the plaque stability.

3D Fiber Orientation in Atherosclerotic Carotid Plaques

Atherosclerotic plaque rupture is the primary trigger of fatal cardiovascular events. Fibrillar collagen in atherosclerotic plaques and their directionality are anticipated to play a crucial role in plaque rupture. This study aimed assessing 3D fiber orientations and architecture in atherosclerotic plaques for the first time. Seven carotid plaques were imaged ex-vivo with a state-of-the-art Diffusion Tensor Imaging (DTI) technique, using a high magnetic field (9.4Tesla) MRI scanner. A 3D spin-echo sequence with uni-polar diffusion sensitizing pulsed field gradients was utilized for DTI and fiber directions were assessed from diffusion tensor measurements. The distribution of the 3D fiber orientations in atherosclerotic plaques were quantified and the principal fiber orientations (circumferential, longitudinal or radial) were determined. Overall, 52% of the fiber orientations in the carotid plaque specimens were closest to the circumferential direction, 34% to the longitudinal direction, and 14% to the radial direction. Statistically no significant difference was measured in the amount of the fiber orientations between the concentric and eccentric plaque sites. However, concentric plaque sites showed a distinct structural organization, where the principally longitudinally oriented fibers were closer to the luminal side and the principally circumferentially oriented fibers were located more abluminally. The acquired unique information on 3D plaque fiber direction will help understanding pathobiological mechanisms of atherosclerotic plaque progression and pave the road to more realistic biomechanical plaque modeling for rupture assessment.

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.

MRI-based biomechanical parameters for carotid artery plaque vulnerability assessment

Carotid atherosclerotic plaques are a major cause of ischaemic stroke. The biomechanical environment to which the arterial wall and plaque is subjected to plays an important role in the initiation, progression and rupture of carotid plaques. MRI is frequently used to characterize the morphology of a carotid plaque, but new developments in MRI enable more functional assessment of carotid plaques. In this review, MRI based biomechanical parameters are evaluated on their current status, clinical applicability, and future developments. Blood flow related biomechanical parameters, including endothelial wall shear stress and oscillatory shear index, have been shown to be related to plaque formation. Deriving these parameters directly from MRI flow measurements is feasible and has great potential for future carotid plaque development prediction. Blood pressure induced stresses in a plaque may exceed the tissue strength, potentially leading to plaque rupture. Multi-contrast MRI based stress calculations in combination with tissue strength assessment based on MRI inflammation imaging may provide a plaque stress-strength balance that can be used to assess the plaque rupture risk potential. Direct plaque strain analysis based on dynamic MRI is already able to identify local plaque displacement during the cardiac cycle. However, clinical evidence linking MRI strain to plaque vulnerability is still lacking. MRI based biomechanical parameters may lead to improved assessment of carotid plaque development and rupture risk. However, better MRI systems and faster sequences are required to improve the spatial and temporal resolution, as well as increase the image contrast and signal-to-noise ratio.

The effects of plaque morphology and material properties on peak cap stress in human coronary arteries

Heart attacks are often caused by rupture of caps of atherosclerotic plaques in coronary arteries. Cap rupture occurs when cap stress exceeds cap strength. We investigated the effects of plaque morphology and material properties on cap stress. Histological data from 77 coronary lesions were obtained and segmented. In these patient-specific cross sections, peak cap stresses were computed by using finite element analyses. The finite element analyses were 2D, assumed isotropic material behavior, and ignored residual stresses. To represent the wide spread in material properties, we applied soft and stiff material models for the intima. Measures of geometric plaque features for all lesions were determined and their relations to peak cap stress were examined using regression analyses. Patient-specific geometrical plaque features greatly influence peak cap stresses. Especially, local irregularities in lumen and necrotic core shape as well as a thin intima layer near the shoulder of the plaque induce local stress maxima. For stiff models, cap stress increased with decreasing cap thickness and increasing lumen radius (R = 0.79). For soft models, this relationship changed: increasing lumen radius and increasing lumen curvature were associated with increased cap stress (R = 0.66). The results of this study imply that not only accurate assessment of plaque geometry, but also of intima properties is essential for cap stress analyses in atherosclerotic plaques in human coronary arteries.

A computer-simulation study on the effects of MRI voxel dimensions on carotid plaque lipid-core and fibrous cap segmentation and stress modeling

Background

The benefits of a decreased slice thickness and/or in-plane voxel size in carotid MRI for atherosclerotic plaque component quantification accuracy and biomechanical peak cap stress analysis have not yet been investigated in detail because of practical limitations.

Methods

In order to provide a methodology that allows such an investigation in detail, numerical simulations of a T1-weighted, contrast-enhanced, 2D MRI sequence were employed. Both the slice thickness (2 mm, 1 mm, and 0.5 mm) and the in plane acquired voxel size (0.62x0.62 mm² and 0.31x0.31 mm²) were varied. This virtual MRI approach was applied to 8 histology-based 3D patient carotid atherosclerotic plaque models.

Results

A decreased slice thickness did not result in major improvements in lumen, vessel wall, and lipid-rich necrotic core size measurements. At 0.62x0.62 mm² in-plane, only a 0.5 mm slice thickness resulted in improved minimum fibrous cap thickness measurements (a 2–3 fold reduction in measurement error) and only marginally improved peak cap stress computations. Acquiring voxels of 0.31x0.31 mm² in-plane, however, led to either similar or significantly larger improvements in plaque component quantification and computed peak cap stress.

Conclusions

This study provides evidence that for currently-used 2D carotid MRI protocols, a decreased slice thickness might not be more beneficial for plaque measurement accuracy than a decreased in-plane voxel size. The MRI simulations performed indicate that not a reduced slice thickness (i.e. more isotropic imaging), but the acquisition of anisotropic voxels with a relatively smaller in-plane voxel size could improve carotid plaque quantification and computed peak cap stress accuracy.

The influence of computational strategy on prediction of mechanical stress in carotid atherosclerotic plaques: comparison of 2D structure-only, 3D structure-only, one-way and fully coupled fluid-structure interaction analyses

BACKGROUND

Compositional and morphological features of carotid atherosclerotic plaques provide complementary information to luminal stenosis in predicting clinical presentations. However, they alone cannot predict cerebrovascular risk. Mechanical stress within the plaque induced by cyclical changes in blood pressure has potential to assess plaque vulnerability. Various modeling strategies have been employed to predict stress, including 2D and 3D structure-only, 3D one-way and fully coupled fluid-structure interaction (FSI) simulations. However, differences in stress predictions using different strategies have not been assessed.

METHODS

Maximum principal stress (Stress-P1) within 8 human carotid atherosclerotic plaques was calculated based on geometry reconstructed from in vivo computerized tomography and high resolution, multi-sequence magnetic resonance images. Stress-P1 within the diseased region predicted by 2D and 3D structure-only, and 3D one-way FSI simulations were compared to 3D fully coupled FSI analysis.

RESULTS

Compared to 3D fully coupled FSI, 2D structure-only simulation significantly overestimated stress level (94.1 kPa [65.2, 117.3] vs. 85.5 kPa [64.4, 113.6]; median [inter-quartile range], p=0.0004). However, when slices around the bifurcation region were excluded, stresses predicted by 2D structure-only simulations showed a good correlation (R(2)=0.69) with values obtained from 3D fully coupled FSI analysis. 3D structure-only model produced a small yet statistically significant stress overestimation compared to 3D fully coupled FSI (86.8 kPa [66.3, 115.8] vs. 85.5 kPa [64.4, 113.6]; p<0.0001). In contrast, one-way FSI underestimated stress compared to 3D fully coupled FSI (78.8 kPa [61.1, 100.4] vs. 85.5 kPa [64.4, 113.7]; p<0.0001).

CONCLUSIONS

A 3D structure-only model seems to be a computationally inexpensive yet reasonably accurate approximation for stress within carotid atherosclerotic plaques with mild to moderate luminal stenosis as compared to fully coupled FSI analysis.

The Influence of Inaccuracies in Carotid MRI Segmentation on Atherosclerotic Plaque Stress Computations

Biomechanical finite element analysis (FEA) based on in vivo carotid magnetic resonance imaging (MRI) can be used to assess carotid plaque vulnerability noninvasively by computing peak cap stress. However, the accuracy of MRI plaque segmentation and the influence this has on FEA has remained unreported due to the lack of a reliable submillimeter ground truth. In this study, we quantify this influence using novel numerical simulations of carotid MRI. Histological sections from carotid plaques from 12 patients were used to create 33 ground truth plaque models. These models were subjected to numerical computer simulations of a currently used clinically applied 3.0 T T1-weighted black-blood carotid MRI protocol (in-plane acquisition voxel size of 0.62 × 0.62 mm2) to generate simulated in vivo MR images from a known underlying ground truth. The simulated images were manually segmented by three MRI readers. FEA models based on the MRI segmentations were compared with the FEA models based on the ground truth. MRI-based FEA model peak cap stress was consistently underestimated, but still correlated (R) moderately with the ground truth stress: R = 0.71, R = 0.47, and R = 0.76 for the three MRI readers respectively (p < 0.01). Peak plaque stretch was underestimated as well. The peak cap stress in thick-cap, low stress plaques was substantially more accurately and precisely predicted (error of −12 ± 44 kPa) than the peak cap stress in plaques with caps thinner than the acquisition voxel size (error of −177 ± 168 kPa). For reliable MRI-based FEA to compute the peak cap stress of carotid plaques with thin caps, the current clinically used in-plane acquisition voxel size (∼0.6 mm) is inadequate. FEA plaque stress computations would be considerably more reliable if they would be used to identify thick-cap carotid plaques with low stresses instead.

Computational approaches for analyzing the mechanics of atherosclerotic plaques: a review

Vulnerable and stable atherosclerotic plaques are heterogeneous living materials with peculiar mechanical behaviors depending on geometry, composition, loading and boundary conditions. Computational approaches have the potential to characterize the three-dimensional stress/strain distributions in patient-specific diseased arteries of different types and sclerotic morphologies and to estimate the risk of plaque rupture which is the main trigger of acute cardiovascular events. This review article attempts to summarize a few finite element (FE) studies for different vessel types, and how these studies were performed focusing on the used stress measure, inclusion of residual stress, used imaging modality and material model. In addition to histology the most used imaging modalities are described, the most common nonlinear material models and the limited number of models for plaque rupture used for such studies are provided in more detail. A critical discussion on stress measures and threshold stress values for plaque rupture used within the FE studies emphasizes the need to develop a more location and tissue-specific threshold value, and a more appropriate failure criterion. With this addition future FE studies should also consider more advanced strain-energy functions which then fit better to location and tissue-specific experimental data.

Numerical simulations of carotid MRI quantify the accuracy in measuring atherosclerotic plaque components in vivo

PURPOSE

Atherosclerotic carotid plaques can be quantified in vivo by MRI. However, the accuracy in segmentation and quantification of components such as the thin fibrous cap (FC) and lipid-rich necrotic core (LRNC) remains unknown due to the lack of a submillimeter scale ground truth.

METHODS

A novel approach was taken by numerically simulating in vivo carotid MRI providing a ground truth comparison. Upon evaluation of a simulated clinical protocol, MR readers segmented simulated images of cross-sectional plaque geometries derived from histological data of 12 patients.

RESULTS

MR readers showed high correlation (R) and intraclass correlation (ICC) in measuring the luminal area (R = 0.996, ICC = 0.99), vessel wall area (R = 0.96, ICC = 0.94) and LRNC area (R = 0.95, ICC = 0.94). LRNC area was underestimated (mean error, -24%). Minimum FC thickness showed a mediocre correlation and intraclass correlation (R = 0.71, ICC = 0.69).

CONCLUSION

Current clinical MRI can quantify carotid plaques but shows limitations for thin FC thickness quantification. These limitations could influence the reliability of carotid MRI for assessing plaque rupture risk associated with FC thickness. Overall, MRI simulations provide a feasible methodology for assessing segmentation and quantification accuracy, as well as for improving scan protocol design.