- 1School of Biomedical Engineering, Biomedical Instrument Institute, Shanghai Jiao Tong University, Shanghai, China
- 2The Lambe Institute for Translational Medicine and Curam, National University of Ireland Galway, Galway, Ireland
- 3Cardiology Department, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China
- 4University of Michigan-Shanghai Jiao Tong University Joint Institute, Shanghai Jiao Tong University, Shanghai, China
- 5School of Biological Science and Medical Engineering, Southeast University, Nanjing, China
- 6Department of Cardiology, Tongji Hospital, Tongji University School of Medicine, Shanghai, China
- 7Mathematical Sciences Department, Worcester Polytechnic Institute, Worcester, MA, United States
Introduction: Cyclic plaque structural stress has been hypothesized as a mechanism for plaque fatigue and eventually plaque rupture. A novel approach to derive cyclic plaque stress in vivo from optical coherence tomography (OCT) is hereby developed.
Materials and Methods: All intermediate lesions from a previous OCT study were enrolled. OCT cross-sections at representative positions within each lesion were selected for plaque stress analysis. Detailed plaque morphology, including plaque composition, lumen and internal elastic lamina contours, were automatically delineated. OCT-derived vessel and plaque morphology were included in a 2-dimensional finite element analysis, loaded with patient-specific intracoronary pressure tracing data, to calculate the changes in plaque structural stress (ΔPSS) on vessel wall over the cardiac cycle.
Results: A total of 50 lesions from 41 vessels were analyzed. A significant ΔPSS gradient was observed across the plaque, being maximal at the proximal shoulder (45.7 [32.3, 78.6] kPa), intermediate at minimal lumen area (MLA) (39.0 [30.8, 69.1] kPa) and minimal at the distal shoulder (35.1 [28.2, 72.3] kPa; p = 0.046). The presence of lipidic plaques were observed in 82% of the diseased segments. Larger relative lumen deformation and ΔPSS were observed in diseased segments, compared with normal segments (percent diameter change: 8.2 ± 4.2% vs. 6.3 ± 2.3%, p = 0.04; ΔPSS: 59.3 ± 48.2 kPa vs. 27.5 ± 8.2 kPa, p < 0.001). ΔPSS was positively correlated with plaque burden (r = 0.37, p < 0.001) and negatively correlated with fibrous cap thickness (r = −0.25, p = 0.004).
Conclusions: ΔPSS provides a feasible method for assessing plaque biomechanics in vivo from OCT images, consistent with previous biomechanical and clinical studies based on different methodologies. Larger ΔPSS at proximal shoulder and MLA indicates the critical sites for future biomechanical assessment.
Introduction
Spontaneous plaque rupture and subsequent thrombosis are recognized as the leading pathogenic mechanism for acute coronary syndrome (ACS), one of the major causes of mortality worldwide (1–3). Thin cap fibroatheroma (TCFA) has been postulated as the phenotype responsible for plaque rupture (3–8). However, limited specificity was observed for TCFA in predicting future coronary events, urging the need to further define meaningful surrogates for rupture-prone plaque identification (9–12).
From a biomechanical point of view, long-term repetitive superficial stress, generated by the pulsatile coronary pressure wave, might weaken the fibrous cap and ultimately lead to its fatigue and rupture (13, 14). Thus, the evaluation of cyclic plaque structural stress might add prognostic value for future cardiac events and subsequently for ACS prevention. Although direct in vivo measurement of plaque structural stress is not currently feasible, finite element analysis (FEA) might provide a reliable estimation (15). The prerequisites for accurate FEA are precise plaque morphology and composition, known mechanical properties of the different materials and precise model loads. Optical coherence tomography (OCT) provides optimal image resolution, enabling detailed visualization and precise characterization of plaque composition (16). Meanwhile, intracoronary pressure tracing from pressure wire could serve as an accurate load for cyclic plaque stress evaluation using FEA. The aim of this study was to propose a novel method to derive the changes in plaque structural stress during the cardiac cycle in vivo using a combination of OCT images and intracoronary pressure recordings.
Materials and Methods
All patients from a previous prospective optical flow ratio (OFR) study with both OCT and fractional flow reserve (FFR) interrogation were screened for post-hoc analysis (17). Inclusion criteria were intermediate coronary lesions, defined as a diameter stenosis 40–90% by visual estimation. Exclusion criteria for FEA were: 1) bifurcation lesions with a side-branch ≥2mm; 2) diffuse coronary disease in the target vessel; 3) stented lesions; 4) intracoronary thrombus; 5) guidewire artifact in the OCT images; 6) incomplete intracoronary pressure recording, with baseline aortic pressure or distal coronary pressure at rest missing. Detailed description of the OCT acquisition and intracoronary pressure measurement has been previously reported (17).
Representative Position Selection
Five OCT cross-sections were selected at representative positions for each lesion, whenever available (Figures 1C1–C5): 1) proximal and distal references (PR, DR): immediately adjacent cross-sections to the lesion, at the proximal and distal edges, respectively, where neither plaque nor remodeling were observed in OCT and angiography; 2) minimal lumen area (MLA): the cross-section with minimal lumen area in OCT; 3) proximal and distal shoulders (PS, DS): midpoints between PR and MLA or DR and MLA, respectively. In case the PR or DR are not present in certain lesion, the midpoint between proximal edge and MLA or distal edge and MLA will be used for PS or DS, respectively.
Figure 1. Representative example of ΔPSS analysis of an intermediate RCA lesion. Cross-sections (C1-C5) correspond to the five representative positions in the angiography showed in panel (A), with (C1) as the cross-section of proximal reference, (C2) as the cross-section of proximal shoulder, (C3) as the cross-section of minimal lumen area, (C4) as the cross-section of distal shoulder and (C5) as the cross-section of distal reference. The 2-dimensional FEA model is loaded with the position-specific intracoronary pressure derived from intracoronary tracing data, by normalizing the computed OFR pullback curve between the resting aortic pressure tracing data and distal coronary pressure tracing data (B). The three-layer geometric models were reconstructed based on an automatic plaque delineation algorithm (18) (D1–D5), where the lipidic plaque is shown in yellow. The thicknesses of media and adventitia were manually measured from OCT cross-sections (C1–C5). (E1–E5) show the stress distribution with the red triangles pointing the positions with largest ΔPSS. FEA, finite element analysis; RCA, right coronary artery; OCT, optical coherence tomography; OFR, optical flow ratio; ΔPSS, delta plaque structural stress.
Geometric Model Reconstruction
The lumen contours of selected OCT cross-sections were automatically delineated using OctPlus software (Pulse Medical Imaging Technology, Shanghai, China). The contour of the internal elastic lamina (IEL) was then identified or extrapolated from adjacent cross-sections, according to a previously validated method (18). The plaque composition was then automatically analyzed at every cross-section in the region encompassed between IEL and the lumen contour, using artificial intelligence (18). The mechanical-relevant plaque components considered for the current FEA in the intima were lipids, calcium and fibrous tissue. Media and adventitia were also incorporated into the geometric model by measuring their thickness on OCT and offsetting the IEL contour uniformly (Figures 1C1–C5,D1–D5). In case media and adventitia were not visible at the cross-section, their thickness on adjacent OCT cross-section was selected.
Biomechanical Analysis
The simulation was performed with the commercially available FEA software ABAQUS (Version 6.13, Dassault Systemes Simulia Corp., Providence, RI, USA). Intimal lipidic plaque and calcific plaque were modeled as isotropic, hyperelastic materials as described by Mooney-Rivlin strain energy density function:
where I1 and I2 are the first and second strain invariants. C1, C2, D1 and D2 are material constants adapted from previous studies (19, 20).
Intimal fibrous tissue, media and adventitia were modeled as anisotropic, hyperelastic materials, according to Holzapfel model:
where μ, k1 and k2 are material constants adapted from previous studies (21–23).
Position-specific pressure condition derived from intracoronary tracing data was applied to the 2-dimensional (2D) FEA model for each cross-section. From OCT images the OFR pullback was firstly computed using a recently validated software package (OctPlus, Pulse medical imaging technology, Shanghai, China) (17, 24, 25). By normalizing the computed OFR pullback curve between the resting aortic pressure and distal coronary pressure tracing data, the position-specific intracoronary pressure at each OCT cross-section could be precisely estimated, also given the excellent agreement between OFR and FFR (17). The change in coronary pressure, i.e., the relative pressure, was then computed by subtracting the estimated diastolic pressure from the cyclic pressure (Figure 1B), which was used as the mechanical load for simulations. The 2D FEA models were then meshed with three-node or four-node linear, hybrid elements. Large deformation formulation and plane strain assumption were used for simulation. The rotational freedom was restricted to prevent the model from rolling while enabling its radial deformation. By submitting the FEA model to ABAQUS/Explicit Solver, the dynamic change of lumen and the plaque structural stress distribution at each cross-section during the cardiac cycle could then be simulated (Supplementary Video 1).
After the simulation, the maximal superficial von Mises stress on the vessel wall at maximal pressure load moment was denoted as OCT-derived change in plaque structural stress (ΔPSS). The thickness of the superficial layer depends on the size of the meshes in the FEA model which was <50 μm. Lumen diameter change (LDC) and percent lumen diameter change (LDC%) during cardiac contraction were used for presenting lumen deformation and relative lumen deformation. The LDC equals to the maximal lumen diameter minus minimal lumen diameter over the cardiac cycle. The LDC% is computed by dividing LDC with the minimal lumen diameter.
Statistics
Descriptive statistics of continuous variables are reported as mean ± SD or median (quartiles) as appropriate, while those of categorical variables are presented as counts (percentages). The difference between groups was tested using independent sample t-test, Mann-Whitney test or One-way Analysis of Variance (ANOVA), as appropriate. Paired t-test, Wilcoxon signed-rank test or repeated measures ANOVA were used for pair-wise comparison, as appropriate. To evaluate the statistical differences of simulated lumen deformation and plaque structural stress at different locations of the lesion, the generalized estimation equation (GEE) analyses were performed. Statistical assessments were performed with MedCalc version 19.5.6 (MedCalc Software, Ostend, Belgium) and SPSS version 27.0.1.0 (SPSS Inc., Chicago, Illinois). A 2-sided value of p < 0.05 was considered to be statistically significant.
Results
Baseline Clinical and Lesion Characteristics
A total of 83 intermediate lesions from 75 vessels were enrolled. Thirty-three lesions were excluded from FEA due to bifurcation (n = 10), diffuse disease (n = 1), stented lesion (n = 3), intracoronary thrombus (n = 3), guidewire artifact (n = 12) or incomplete intracoronary recording (n = 4), resulting in 50 lesions suitable for biomechanical analysis (Figure 2).
Baseline demographic and lesion characteristics are presented in Tables 1, 2, respectively. The mean age of the patients was 63 ± 11 years. Half of the lesions (50.0%) were located at the proximal segment of the target vessel.
Plaque Morphology and Composition
Quantification of the OCT images was performed using the OctPlus software package (Pulse medical imaging technology, Shanghai, China). Automatic plaque characterization and delineation of the IEL from OCT images by the software was performed using artificial intelligence algorithm and recently validated with high accuracy (18). Average lumen diameter, plaque burden and fibrous cap thickness in the representative positions are presented in Table 3. MLA had the smallest lumen diameter (1.82 [1.60, 1.97] mm) and the largest plaque burden (69.11 ± 9.11%) of all the representative positions. Fibrous cap thickness was numerically smaller in PS than in MLA and DS, but it did not reach statistical significance (PS vs. MLA vs. DS: 154.1 [60.8, 242.2] μm vs. 168.9 [54.2, 264.7] μm vs. 212.4 [92.7, 253.2] μm, p = 0.48). The plaque composition and microfeatures of the lesions included in current study are presented in Table 4. Forty-one lesions (82.0%) have more than one plaque phenotypes. By using the predominant plaque phenotype for each lesion, a total of 28 fibroatheromas, 8 fibrotic plaques and 14 fibrocalcific plaques were included in the present study.
Among all 177 interrogated cross-sections, more than half of them (92 cross-sections) have a good (external elastic lamina circumference ≥270°) or moderate visibility (270° > external elastic lamina circumference ≥ 180°) of the media (54 and 38 cross-sections, respectively). While the media were completely invisible in 14 (7.9%) cross-sections, where the artificial intelligence algorithm used their adjacent cross-sections for IEL delineation. In these cases, the mean number of skipped cross-sections for manual media and adventitia thicknesses measurement was 5.6 ± 6.8 frames.
Lumen Deformation and Plaque Structural Stress
A representative example of ΔPSS analysis is shown in Figure 1 and Supplementary Video 1. The LDC, LDC% and ΔPSS across the lesion are presented in Table 5 and Figure 3. Absolute LDC was smaller at MLA than at both reference cross-sections (PR and DR), but LDC% was larger. Positive correlation was observed between LDC and lumen diameter (r = 0.53, p < 0.01). However, no correlation was observed between LDC% and lumen diameter.
Figure 3. (A–C) Biomechanical results in different locations of the lesion. *p < 0.001. PR, proximal reference; PS, proximal shoulder; MLA, minimal lumen area; DS, distal shoulder; DR, distal reference; LDC, lumen diameter change; LDC%, percent lumen diameter change; ΔPSS, delta plaque structural stress. The blue circles represent the mild outliers, and the red squares represent the extreme outliers.
ΔPSS was significantly smaller at the reference cross-sections than at MLA or at the shoulders (PR vs. PS, MLA and DS: 25.2 [21.2, 32.3] kPa vs. 45.7 [32.3, 78.6] kPa, 39.0 [30.8, 69.1] kPa and 35.1 [28.2, 72.3] kPa; p < 0.001, p < 0.001 and p < 0.001, respectively). In paired analysis per lesion, a significant ΔPSS gradient was observed across the plaque (PS-MLA-DS, p = 0.046; test for linear trend: p = 0.012). ΔPSS at PS was significantly larger than at DS (p = 0.02), while ΔPSS at MLA tended to have intermediate values between PS and DS (Figure 3). Within all 49 lesions with lipidic plaques, the highest ΔPSS occurred at the cross-sections with thinnest cap in 49.0% (24) of these lesions. While for 22 lesions (44.9%), the highest ΔPSS occurred at the cross-sections with lipidic plaques but thicker cap thickness. Three (6.1%) lesions were observed to have the highest ΔPSS at cross-sections without any lipidic plaque.
Normal vs. Diseased Segments
The analyzed OCT cross-sections were divided into two groups: 1) normal segments, comprising PR and DR; 2) diseased segments, comprising PS, MLA and DS. LDC was numerically smaller in diseased segments compared to normal segments, but statistically non-significant (0.17 [0.10, 0.27] mm vs. 0.19 [0.13, 0.25] mm, p = 0.69). While both LDC% and ΔPSS were significantly larger in diseased segments than in normal segments (LDC%: 7.57 [5.38, 10.64]% vs. 6.29 [4.44, 8.17]%, p = 0.002; ΔPSS: 41.8 [29.1, 74.6] kPa vs. 27.7 [21.2, 32.1] kPa, p < 0.001) (Figure 4).
Figure 4. (A–C) Comparison between normal segment and diseased segment. LDC, lumen diameter change; LDC%, percent lumen diameter change; ΔPSS, delta plaque structural stress. The blue circles represent the mild outliers, and the red boxes represent the extreme outliers.
In normal segments, a good correlation was observed between maximal pressure load and ΔPSS and (r = 0.72, p < 0.001). Nonetheless, this correlation was significantly worse in diseased segments (r = 0.42, p < 0.001; difference p = 0.04), wherein relatively low pressure loads often resulted in high ΔPSS.
Correlation of Stress Parameters With Morphological Features of Plaque Vulnerability
ΔPSS was positively correlated with plaque burden (r = 0.37, p < 0.001) while negatively correlated with fibrous cap thickness (r = −0.25, p < 0.001). Good correlation was observed between ΔPSS and LDC% (r = 0.78, p < 0.001).
Moderate correlation was observed between ΔPSS and lipidic plaque area (LPA) (r = 0.44, p < 0.001), while no significant correlation was observed between ΔPSS and calcific plaque area. LPA showed good correlation with both lumen deformation and relative lumen deformation (LPA-LDC: r = 0.55, p < 0.001; LPA-LDC%: r = 0.88, p < 0.001).
Computational Performance of ΔPSS Assessment
Using an off-the-shelf workstation with a quadcore Intel i7-4790 processor (Intel Corporation, Santa Clara, CA; 3.6 GHz) and 8GB of RAM, the average simulation time for analysis of each FEA model was 12.5 ± 10.2 min, and the estimated time for the whole FEA analysis process was <30 min.
Discussion
In this study, we present for the first time a new methodology to calculate the changes in plaque structural stress within the cardiac cycle in vivo using OCT images and FEA simulation. The changes of superficial plaque structural stress on vessel wall are the main focus of our current study since it tends to be closely related to plaque rupture and subsequent acute coronary events (26–29).
The key findings of this study can be summarized as follows: 1) ΔPSS provides a feasible and reasonable approach for OCT-based biomechanical assessment; 2) diseased segments, especially the proximal shoulder and minimal lumen area of the lesion, bear the highest ΔPSS, thus highlighting the critical importance of these sites for future biomechanical studies of plaque vulnerability and prediction of event risk; 3) correlation between ΔPSS and plaque morphology is consistent with previous clinical and imaging studies, thus reassuring the rationale of our method.
This novel method is original in many aspects, including the automatic plaque characterization from OCT images using artificial intelligence, the incorporation of the three-layered structure of the vessel wall into 2D FEA model, the accurate position-specific load derived from intracoronary pressure tracing data. Previous studies performed coronary plaque stress simulation (27), but most of them were based on coronary angiography or intravascular ultrasound (IVUS), with inherent limitations in plaque characterization and inferior image resolution, as compared with OCT. Additionally, rigid-body assumption restricting all deformation was applied to most FEA models hitherto (20, 30, 31).
The first advantage of the current approach is the higher image resolution provided by OCT, enabling the detailed description of coronary lumen, plaque morphology and composition which are essential for FEA analysis. The significant lower correlation between load and ΔPSS in diseased segments than in normal segments (r = 0.42 vs. 0.72, p < 0.001) strongly suggests a crucial role of plaque morphology and composition in determining the actual plaque stress. Chau et al. have also proposed an FEA method to derive plaque stress from OCT images (32). However, a fixed pressure from 0 to 120 mmHg was applied to all models, which might not be realistic since the geometry of the vessel was not imaged under zero-pressure condition. In addition, the same isotropic material properties were applied to the whole artery wall, ignoring the anisotropy of media and adventitia and the mechanical differences between them. Although OCT has limited tissue penetration, as compared with IVUS, and notwithstanding the high attenuation of near-infrared waves in lipids, apparently precluding the imaging of media and adventitia in diseased segments (16), recent studies have proven that the external elastic lamina can be identified for ≥180° in most of the OCT cross-sections (33). Thus, the media contour could be extrapolated considering its circular or elliptical geometry. For OCT cross-sections with invisible IEL, which only accounts for 7.9% of the total cross-sections in this study, our artificial intelligence algorithm will refer to adjacent OCT cross-sections and estimate the contour of IEL. This principle was externally validated for the same specific software used in the current study, providing high diagnostic accuracy for delineation of the media contour and tissue characterization within the plaque, including lipids (90.5%) (18). In this study, we mainly focused on the superficial cyclic plaque structural stress on vessel wall, which is more likely to be affected by the different plaque components in near-lumen regions, as compared to outer layers as media and adventitia. During the analysis, we also observed that the slight deviation, i.e., thickness and contour, in media and adventitia have a relatively small impact on simulation results. While for superficial layers, even a tiny change, e.g., lumen contour and fibrous cap thickness, will cause great difference in ΔPSS. It has also been proven by previous studies that the most crucial impact on stress computation comes from the superficial wall adjacent to lumen, and that the vascular structure determined by OCT provides adequate basis for biomechanical analysis (23, 32). Nevertheless, when focusing on the stress distribution at deeper layer, e.g., behind the plaque or near the outer vessel boundary, more precise simulation results might be achieved by incorporating more accurate delineation of media and adventitia, especially for lesions with large lipidic plaque burden where the outer boundary is invisible due to light attenuation in OCT. In such cases, IVUS with higher tissue penetration could serve as a complemental modality (34). By combining detailed evaluation of plaque morphology provided by OCT and assessment of the entire vessel structure provided by IVUS, a more precise geometric model for FEA should become obtainable, leading to a more accurate biomechanical assessment especially for outer layers (35). In the future, with the development of automatic OCT-IVUS co-registration algorithm and/or hybrid intravascular imaging system with combined OCT and IVUS probes (36), the accuracy of PSS assessment could be further improved with great efficiency.
A complete vessel architecture with three-layer structure was applied for simulation in this study. Very few studies have incorporated the three-layered structure of the vessel into the FEA simulation hitherto. This might be a potential limitation of previous studies because stress distribution heavily depends on the physical properties of the material and the layered arrangement of the vessel wall (21, 37–39). The different composition and stiffness of intima, media and adventitia determine their different roles during loading. The stiffness of intima widely varies depending on the different plaque components. Conversely, media and adventitia use to have more predictable mechanical properties: the structured arrangement of smooth muscle cells in the media confers this layer high resistance to load, while the helically arranged wavy collagen fibrils and elastic fibers in adventitia make it more compliant to pressure than media under normal load (40, 41).
Accurate model load using intracoronary pressure registration has been instrumental to render precise stress estimations. The model, however, was only loaded with the change in coronary pressure, instead of using the whole pressure recording, in order to simplify and expedite the calculation. The stress change over the cardiac cycle in atherosclerotic plaques has been proven to correlate with the incidence of major adverse cardiovascular events (42) and the general relation between plaque morphology and plaque stress remains intact irrespective of the inclusion or exclusion of the initial stress in the model (43). On this rationale, the decision of using only pressure change seems justified. Loading the model with the whole pressure recording might have added marginal accuracy to the estimations, but at the price of exponentially increasing the complexity of the analysis and subsequently the time required for it (43–45). For models to be loaded with the whole pressure recording, accurate computation of zero-pressure state is one of the prerequisites. Several methods have been proposed for the estimation of initial stress. A preshrinkage algorithm by Huang et al. iteratively shrinks the in-vivo plaque geometry before the whole pressure loading to estimate the zero-pressure state and initial stress (44, 45). However, manual adaptation was required in each iteration to compare the computed geometry with the image-derived real geometry and adjust the geometry for the next iteration, which might be too complicated for daily practice. Speelman et al. proposed a Backward Incremental method that requires no manual input (43). However, the vessel geometry under certain fixed intracoronary pressure is one of the prerequisites for the intimal stress estimation, making this method not suitable for in vivo assessment. Theoretically, zero-pressure state and initial stress could be accurately estimated from the geometrical models of the same cross-section at both diastole and systole, where the imaging catheter is required to be fixed at the same position for at least one cardiac cycle. However, this kind of relative stillness between catheter and vessel is technically difficult to achieve, considering the impact of cardiac motion and vessel contraction (46). Thus, the gain in incorporating whole pressure recording might be partly canceled by the relative displacement between the imaging catheter and the analyzed cross-section during cardiac cycle.
The consistency of the plaque stress calculated by this novel method with previous clinical findings and biomechanical studies is reassuring of the validity of our approach. Autopsy studies have linked high plaque stress with plaque rupture (13, 47). Likewise, imaging studies have associated high plaque stress with acute coronary syndromes (14). This evidence supports the biomechanical hypothesis for plaque rupture: the repetitive cyclic stress would end up breaking the plaque, like repetitively bending a paper leads to its weakening and fracture (1). In line with this evidence, normal coronary segments bear significantly lower ΔPSS than diseased segments in our study. Moreover, it is known that strain and plaque stress are higher upstream than downstream in the atherosclerotic plaque (48, 49), i.e., higher in the proximal segments than in the distal segments of the plaque, and that most plaque ruptures occur in these proximal segments of the lesion (30, 48, 50, 51). Our results, finding a plaque stress gradient from the PS to the MLA and ultimately to the DS, are in line with this preceding evidence. In addition, ΔPSS showed positive correlation with LPA and negative correlation with fibrous cap thickness, thus confirming a direct association between ΔPSS and morphologically rupture-prone plaque, i.e., TCFA. Nonetheless, both correlations were relatively weak (r = 0.37 and−0.25, respectively) and only around half (49.0%) of the lesions with lipidic plaques have the highest ΔPSS located at the thinnest cap sites, thus suggesting a mismatch between the histological and biomechanical evaluation of vulnerability. This observation may at least partially account for the limited prognostic value of TCFA alone (9–12).
In this study, the mean LDC% in normal segment is 6.29%, slightly smaller than the range of 7–14% observed by several previous studies from groups of healthy people aged from 8 to 60 (52–54). Since the coronary artery stiffens with age and disease, the relatively smaller LDC% of 6.29% in this study seems reasonable considering that our material property for intimal fibrous tissue was adopted from patients with an average age of 66 (21). The LDC is larger in normal segments compared to the diseased segments, though statistically non-significant. This finding is in line with previous IVUS studies that normal segments are more compliant to deformation (55). Conversely, significantly larger LDC% in diseased segments were observed in the present study (7.57% vs. 6.29%, p = 0.04). A possible explanation for this phenomenon is that 82% of the diseased segments have lipidic plaques, which are softer and thus tend to have larger relative deformation. Good correlation between LPA and %LDC further elucidates this point (r = 0.88, p < 0.001).
Clinical Perspectives
TCFA is currently recognized as a precursor for plaque rupture. However, several issues regarding image-based TCFA detection remain unsolved, including modest interobserver agreement and inconsistent definitions between studies (56). In addition, most TCFAs do not cause symptomatic rupture (9–11), revealing the fact that histological assessment alone is not enough for ACS prevention.
Plaque stress evaluation based on FEA might serve as a supplementary strategy. The location of peak cyclic plaque stress might be helpful to predict the risk of rupture and hence for ACS prevention. However, the length of a focal lesion is 12–30 mm, with 60–150 OCT cross-sections at the highest pullback speed (57, 58). Therefore, plaque stress evaluation would be very time-consuming if the whole lesion were scanned. For future studies, limiting the analysis to high-risk locations, especially for lipid rich plaques which are at higher propensity to rupture, might be instrumental for fast and efficient risk stratification.
Interestingly, LDC% showed good correlation with both ΔPSS (r = 0.78, p < 0.001) and LPA (r = 0.88, p < 0.001) in our study. These findings invite us to explore the possibility of using angiography-derived lumen deformation for the estimation of lipidic burden and for simplified identification of high-risk sites for rupture. Considering the relatively short modeling and simulation time of ΔPSS, the hereby described method might also be integrated into the current OCT-based FFR computation system, where the morphological, histological and biomechanical assessments could be achieved efficiently within one single OCT pullback.
Limitations
The current study is limited by its post-hoc design and its relatively small sample size. Nevertheless, all intermediate lesions with both OCT and intracoronary pressure tracing were enrolled, following strictly predefined inclusion/exclusion criteria, thus minimizing the selection bias. Besides, not all interrogated lesions had PR and DR.
In this study, the FEA models were reconstructed in 2D without incorporating the shear stress and mechanical forces in the axial direction (59–61). The homogenous material property within lipidic and calcific plaque and the uniform assumption for media and adventitia might also introduce error into the simulation. Besides, the effect of residual stress was not considered since it was currently immeasurable. The lumen deformation might also be different from real situation depending on the prevailing diastolic blood pressure levels. In addition, for cross-sections with invisible IEL, the vessel boundaries estimated by the artificial intelligence might be less accurate and lead to imprecise simulation results. However, the impact of this limitation on superficial ΔPSS is negligible. Although the simulation results are in line with clinical findings, due to the complexity of our current approach combining OCT, FFR and FEA model together, there is no single “gold standard” for validation. Besides, only presentative positions were analyzed for each lesion in this feasibility study. We did not perform frame-by-frame analysis considering that the reliability in tissue characterization might be impaired for cross-sections with side branches. Future prospective studies with larger sample sizes are warranted to investigate the prognostic value and the clinical usefulness of ΔPSS and other ΔPSS-related parameters.
The current approach, using intracoronary pressure recordings, is exquisitely accurate, but it makes the method complex and expensive to routine clinical implementation. In addition, only OCT images without guidewire artifact were analyzed in this pilot study. Future simplifications of this approach might facilitate the applicability of this assessment in the cathlab, provided they rendered acceptable accuracy.
Conclusions
Plaque structural stress over the cardiac cycle can be estimated from OCT images, using automatic plaque characterization and FEA, on a feasible fully automated process aided by artificial intelligence. The results of this novel approach are consistent with previous clinical and biomechanical studies, showing higher plaque stress in diseased vs. normal segments. The highest stress at the proximal shoulder and MLA indicates the critical rupture-prone sites for efficient biomechanical assessment in the future.
Data Availability Statement
The dataset can be shared upon reasonable request. Requests to access these datasets should be directed to the corresponding author.
Ethics Statement
The studies involving human participants were reviewed and approved by Local Ethic Committee of Campo de Gibraltar Health Trust. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
JH: analysis of the data and draft of the manuscript. WW and ST: concept design and draft of the manuscript. FY, TX, RL, JW, YL, and XL: analysis of the data. JG-C, PS, YO, LW, and DT: crucial revision of the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported in part by grants from the National Natural Science Foundation of China [grant numbers 82020108015 and 81871460] to ST and [grant number 11972117] to LW, and by the Science Foundation Ireland Research Professorship grant RSF 1413 to WW, and by grants from the Science and Technology Committee of Shanghai Municipality (No. 19411963200) to YL.
Conflict of Interest
WW reports research grants and honoraria from MicroPort, medical advisor of Rede Optimus Research and co-founder of Argonauts, an innovation facilitator. ST reports research grants and consultancy from Pulse medical imaging technology. PS reports personal fees from Sino Medical Science Technology, Philips/Volcano, Xeltis, HeartFlow, and SMT, outside the submitted work.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2021.715995/full#supplementary-material
References
1. Falk E, Shah PK, Fuster V. Coronary plaque disruption. Circulation. (1995) 92:657–71. doi: 10.1161/01.CIR.92.3.657
2. Burke AP, Farb A, Malcom GT, Liang Y, Smialek JE, Virmani R. Coronary risk factors and plaque morphology in men with coronary disease who died suddenly. N Engl J Med. (1997) 336:1276–82. doi: 10.1056/NEJM199705013361802
3. Virmani R, Kolodgie FD, Burke AP, Farb A, Schwartz SM. Lessons from sudden coronary death: a comprehensive morphological classification scheme for atherosclerotic lesions. Arterioscler Thromb Vasc Biol. (2000) 20:1262–75. doi: 10.1161/01.ATV.20.5.1262
4. Cheruvu PK, Finn AV, Gardner C, Caplan J, Goldstein J, Stone GW, et al. Frequency and distribution of thin-cap fibroatheroma and ruptured plaques in human coronary arteries: a pathologic study. J Am Coll Cardiol. (2007) 50:940–9. doi: 10.1016/j.jacc.2007.04.086
5. Finn AV, Nakano M, Narula J, Kolodgie FD, Virmani R. Concept of vulnerable/unstable plaque. Arterioscler Thromb Vasc Biol. (2010) 30:1282–92. doi: 10.1161/ATVBAHA.108.179739
6. Kolodgie FD, Virmani R, Burke AP, Farb A, Weber DK, Kutys R, et al. Pathologic assessment of the vulnerable human coronary plaque. Heart. (2004) 90:1385–91. doi: 10.1136/hrt.2004.041798
7. Virmani R, Burke AP, Farb A, Kolodgie FD. Pathology of the vulnerable plaque. J Am Coll Cardiol. (2006) 47:C13–8. doi: 10.1016/j.jacc.2005.10.065
8. Naghavi M, Libby P, Falk E, Casscells SW, Litovsky S, Rumberger J, et al. From vulnerable plaque to vulnerable patient. Circulation. (2003) 108:1664–72. doi: 10.1161/01.CIR.0000087480.94275.97
9. Stone GW, Maehara A, Lansky AJ, de Bruyne B, Cristea E, Mintz GS, et al. A prospective natural-history study of coronary atherosclerosis. N Engl J Med. (2011) 364:226–35. PubMed PMID: 21247313. doi: 10.1056/NEJMoa1002358
10. Calvert PA, Obaid DR, Osullivan M, Shapiro LM, Mcnab DC, Densem CG, et al. Association between IVUS findings and adverse outcomes in patients with coronary artery disease: the VIVA (VH-IVUS in Vulnerable Atherosclerosis) study. JACC-Cardiovasc Interv. (2011) 4:894–901. doi: 10.1016/j.jcmg.2011.05.005
11. Cheng JM, Garcia-Garcia HM, de Boer SPM, Kardys I, Heo JH, Akkerhuis KM, et al. In vivo detection of high-risk coronary plaques by radiofrequency intravascular ultrasound and cardiovascular outcome: results of the ATHEROREMO-IVUS study. Eur Heart J. (2013) 35:639–47. doi: 10.1093/eurheartj/eht484
12. Prati F, Romagnoli E, Gatto L, La Manna A, Burzotta F, Ozaki Y, et al. Relationship between coronary plaque morphology of the left anterior descending artery and 12 months clinical outcome: the CLIMA study. Eur Heart J. (2019) 41:383–91. doi: 10.1093/eurheartj/ehz520
13. Cheng GC, Loree HM, Kamm RD, Fishbein MC, Lee RT. Distribution of circumferential stress in ruptured and stable atherosclerotic lesions. A structural analysis with histopathological correlation. Circulation. (1993) 87:1179–87. doi: 10.1161/01.cir.87.4.1179
14. Schaar JA, Regar E, Mastik F, McFadden EP, Saia F, Disco C, et al. Incidence of high-strain patterns in human coronary arteries. Circulation. (2004) 109:2716–9. doi: 10.1161/01.CIR.0000131887.65955.3B
15. Gutierrez-Chico JL. Superficial wall stress: the long awaited comprehensive biomechanical parameter to objectify and quantify our intuition. Int J Cardiovasc Imaging. (2018) 34:863–5. doi: 10.1007/s10554-018-1386-1
16. Gutiérrez-Chico JL, Alegría-Barrero E, Teijeiro-Mestre R, Chan PH, Tsujioka H, de Silva R, et al. Optical coherence tomography: from research to practice. Eur Heart J-Cardiovasc Imaging. (2012) 13:370–84. doi: 10.1093/ehjci/jes025
17. Gutierrez-Chico JL, Chen Y, Yu W, Ding D, Huang J, Huang P, et al. Diagnostic accuracy and reproducibility of optical flow ratio for functional evaluation of coronary stenosis in a prospective series. Cardiology J. (2020) 27:350–61. doi: 10.5603/CJ.a2020.0071
18. Chu M, Jia H, Gutiérrez-Chico JL, Maehara A, Ali ZA, Zeng X, et al. Automatic characterisation of human atherosclerotic plaque composition from intravascular optical coherence tomography using artificial intelligence. EuroIntervention. (2021) 17:41–50. doi: 10.4244/EIJ-D-20-01355
19. Huang H, Virmani R, Younis H, Burke AP, Kamm RD, Lee RT. The impact of calcification on the biomechanical stability of atherosclerotic plaques. Circulation. (2001) 103:1051–6. doi: 10.1161/01.CIR.103.8.1051
20. Teng Z, Brown AJ, Calvert PA, Parker RA, Obaid DR, Huang Y, et al. Coronary plaque structural stress is associated with plaque composition and subtype and higher in acute coronary syndrome: the BEACON I (Biomechanical Evaluation of Atheromatous Coronary Arteries) study. Circ-Cardiovasc Imaging. (2014) 7:461–70. doi: 10.1161/CIRCIMAGING.113.001526
21. Holzapfel GA, Sommer G, Gasser CT, Regitnig P. Determination of layer-specific mechanical properties of human coronary arteries with nonatherosclerotic intimal thickening and related constitutive modeling. Am J Physiol-Heart Circul Physiol. (2005) 289:H2048–58. doi: 10.1152/ajpheart.00934.2004
22. Hollander Y, Durban D, Lu X, Kassab GS, Lanir Y. Constitutive modeling of coronary arterial media—comparison of three model classes. J Biomech Eng. (2011) 133:061008. doi: 10.1115/1.4004249
23. Akyildiz AC, Speelman L, van Brummelen H, Gutiérrez MA, Virmani R, van der Lugt A, et al. Effects of intima stiffness and plaque morphology on peak cap stress. BioMed Eng Online. (2011) 10:25. doi: 10.1186/1475-925X-10-25
24. Yu W, Huang J, Jia D, Chen S, Raffel OC, Ding D, et al. Diagnostic accuracy of intracoronary optical coherence tomography-derived fractional flow reserve for assessment of coronary stenosis severity. EuroIntervention. (2019) 15:189–97. doi: 10.4244/EIJ-D-19-00182
25. Huang J, Emori H, Ding D, Kubo T, Yu W, Huang P, et al. Comparison of diagnostic performance of intracoronary optical coherence tomography-based and angiography-based fractional flow reserve for evaluation of coronary stenosis. EuroIntervention. (2020) 16:568–76. doi: 10.4244/EIJ-D-19-01034
26. Kwak BR, Bäck M, Bochaton-Piallat M-L, Caligiuri G, Daemen MJAP, Davies PF, et al. Biomechanical factors in atherosclerosis: mechanisms and clinical implications†. Eur Heart J. (2014) 35:3013–20. doi: 10.1093/eurheartj/ehu353
27. Brown AJ, Teng Z, Calvert PA, Rajani NK, Hennessy O, Nerlekar N, et al. Plaque structural stress estimations improve prediction of future major adverse cardiovascular events after intracoronary imaging. Circ-Cardiovasc Imaging. (2016) 9:e004172. doi: 10.1161/CIRCIMAGING.115.004172
28. Ohayon J, Finet G. Le Floc'h S, Cloutier G, Gharib AM, Heroux J, et al. Biomechanics of atherosclerotic coronary plaque: site, stability and in vivo elasticity modeling. Ann Biomed Eng. (2014) 42:269–79. doi: 10.1007/s10439-013-0888-1
29. Wu X, von Birgelen C, Zhang S, Ding D, Huang J, Tu S. Simultaneous evaluation of plaque stability and ischemic potential of coronary lesions in a fluid-structure interaction analysis. Int J Cardiovasc Imaging. (2019) 35:1563–72. doi: 10.1007/s10554-019-01611-y
30. Costopoulos C, Huang Y, Brown A, Calvert PA, Hoole SP, West N, et al. Plaque rupture in coronary atherosclerosis is associated with increased plaque structural stress. JACC-Cardiovasc Interv. (2017) 10:1472–83. doi: 10.1016/j.jcmg.2017.04.017
31. Akyildiz AC, Speelman L. Velzen Bv, Stevens RRF, Steen AFWvd, Huberts W, et al. Intima heterogeneity in stress assessment of atherosclerotic plaques. Interface Focus. (2018) 8:20170008. doi: 10.1098/rsfs.2017.0008
32. Chau AH, Chan R, Shishkov M, Macneill BD, Iftimia N, Tearney GJ, et al. Mechanical analysis of atherosclerotic plaques based on optical coherence tomography. Ann Biomed Eng. (2004) 32:1494–503. doi: 10.1114/B:ABME.0000049034.75368.4a
33. Ali ZA, Maehara A, Généreux P, Shlofmitz RA, Fabbiocchi F, Nazif TM, et al. Optical coherence tomography compared with intravascular ultrasound and with angiography to guide coronary stent implantation (ILUMIEN III: OPTIMIZE PCI): a randomised controlled trial. Lancet. (2016) 388:2618–28. doi: 10.1016/S0140-6736(16)31922-5
34. Fujii K, Hao H, Shibuya M, Imanaka T, Fukunaga M, Miki K, et al. Accuracy of OCT, grayscale IVUS, and their combination for the diagnosis of coronary TCFA: an ex vivo validation study. JACC-Cardiovasc Imag. (2015) 8:451–60. doi: 10.1016/j.jcmg.2014.10.015
35. Doradla P, Otsuka K, Nadkarni A, Villiger M, Karanasos A. Zandvoort LJv, et al. Biomechanical stress profiling of coronary atherosclerosis: Identifying a multifactorial metric to evaluate plaque rupture risk. JACC-Cardiovasc Imag. (2020) 13:804–16. doi: 10.1016/j.jcmg.2019.01.033
36. Ono M, Kawashima H, Hara H, Gao C, Wang R, Kogame N, et al. Advances in IVUS/OCT and future clinical perspective of novel hybrid catheter system in coronary imaging. Front Cardiovasc Med. (2020) 7:119. doi: 10.3389/fcvm.2020.00119
37. Lu X, Yang J, Zhao JB, Gregersen H, Kassab GS. Shear modulus of porcine coronary artery: contributions of media and adventitia. Am J Physiol-Heart Circul Physiol. (2003) 285:H1966–75. doi: 10.1152/ajpheart.00357
38. Thondapu V, Bourantas CV, Foin N, Jang I-K, Serruys PW, Barlis P. Biomechanical stress in coronary atherosclerosis: emerging insights from computational modelling. Eur Heart J. (2016) 38:81–92. doi: 10.1093/eurheartj/ehv689
39. Holzapfel GA, Sommer G, Regitnig P. Anisotropic mechanical properties of tissue components in human atherosclerotic plaques. J Biomech Eng. (2004) 126:657–65. doi: 10.1115/1.1800557
40. Holzapfel GA, Gasser TC, Ogden RW. A new constitutive framework for arterial wall mechanics and a comparative study of material models. J Elast. (2000) 61:1–48. doi: 10.1023/A:1010835316564
41. Baldewsing RA, de Korte CL, Schaar JA, Mastik F, van der Steen AFW. Finite element modeling and intravascular ultrasound elastography of vulnerable plaques: parameter variation. Ultrasonics. (2004) 42:723–9. doi: 10.1016/j.ultras.2003.11.017
42. Costopoulos C, Maehara A, Huang Y, Brown AJ, Gillard JH, Teng Z, et al. Heterogeneity of plaque structural stress Is increased in plaques leading to MACE: insights from the PROSPECT study. JACC-Cardiovasc Imag. (2020) 13:1206–18. doi: 10.1016/j.jcmg.2019.05.024
43. Speelman L, Akyildiz AC, den Adel B, Wentzel JJ, van der Steen AFW, Virmani R, et al. Initial stress in biomechanical models of atherosclerotic plaques. J Biomech. (2011) 44:2376–82. doi: 10.1016/j.jbiomech.2011.07.004
44. Huang X, Yang C, Yuan C, Liu F, Canton G, Zheng J, et al. Patient-specific artery shrinkage and 3D zero-stress state in multi-component 3D FSI models for carotid atherosclerotic plaques based on in vivo MRI data. Mol Cell Biomech. (2009) 6:121–34. doi: 10.3970/mcb.2009.006.121
45. Huang Y, Teng Z, Sadat U, Hilborne S, Young VE, Graves MJ, et al. Non-uniform shrinkage for obtaining computational start shape for in-vivo MRI-based plaque vulnerability assessment. J Biomech. (2011) 44:2316–9. doi: 10.1016/j.jbiomech.2011.06.014
46. Arbab-Zadeh A, DeMaria AN, Penny WF, Russo RJ, Kimura BJ, Bhargava V. Axial movement of the intravascular ultrasound probe during the cardiac cycle: implications for three-dimensional reconstruction and measurements of coronary dimensions. Am Heart J. (1999) 138:865–72. doi: 10.1016/S0002-8703(99)70011-6
47. Lee RT, Loree HM, Cheng GC, Lieberman EH, Jaramillo N, Schoen FJ. Computational structural analysis based on intravascular ultrasound imaging before in vitro angioplasty: prediction of plaque fracture locations. J Am Coll Cardiol. (1993) 21:777–82. doi: 10.1016/0735-1097(93)90112-E
48. Fukumoto Y, Hiro T, Fujii T, Hashimoto G, Fujimura T, Yamada J, et al. Localized elevation of shear stress is related to coronary plaque rupture. J Am Coll Cardiol. (2008) 51:645–50. doi: 10.1016/j.jacc.2007.10.030
49. Gijsen FJH, Wentzel JJ, Thury A, Mastik F, Schaar JA, Schuurbiers JCH, et al. Strain distribution over plaques in human coronary arteries relates to shear stress. Am J Physiol-Heart Circul Physiol. (2008) 295:H1608–14. PubMed PMID: 18621851. doi: 10.1152/ajpheart.01081.2007
50. Lee JM, Choi G, Hwang D, Park J, Kim HJ, Doh J-H, et al. Impact of longitudinal lesion geometry on location of plaque rupture and clinical presentations. JACC-Cardiovasc Imag. (2017) 10:677–88. doi: 10.1016/j.jcmg.2016.04.012
51. Wentzel JJ, Chatzizisis YS, Gijsen FJH, Giannoglou GD, Feldman CL, Stone PH. Endothelial shear stress in the evolution of coronary atherosclerotic plaque and vascular remodelling: current understanding and remaining questions. Cardiovasc Res. (2012) 96:234–43. doi: 10.1093/cvr/cvs217
52. Buntin CM, Silver FH. Noninvasive assessment of mechanical properties of peripheral arteries. Ann Biomed Eng. (1990) 18:549–66. doi: 10.1007/BF02364617
53. Lee RT, Kamm RD. Vascular mechanics for the cardiologist. J Am Coll Cardiol. (1994) 23:1289–95. doi: 10.1016/0735-1097(94)90369-7
54. Shimazu T, Hori M, Mishima M, Kitabatake A, Kodama K, Nanto S, et al. Clinical assessment of elastic properties of large coronary arteries: pressure-diameter relationship and dynamic incremental elastic modulus. Int J Cardiovasc Imaging. (1986) 13:27–45. doi: 10.1016/0167-5273(86)90077-X
55. Shaw JA, Kingwell BA, Walton AS, Cameron JD, Pillay P, Gatzka CD, et al. Determinants of coronary artery compliance in subjects with and without angiographic coronary artery disease. J Am Coll Cardiol. (2002) 39:1637–43. doi: 10.1016/S0735-1097(02)01842-9
56. Ali ZA, Karimi Galougahi K, Maehara A, Shlofmitz RA, Ben-Yehuda O, Mintz GS, et al. Intracoronary optical coherence tomography 2018: current status and future directions. JACC-Cardiovasc Interv. (2017) 10:2473–87. doi: 10.1016/j.jcin.2017.09.042
57. Smith SC, Dove JT, Jacobs AK, Ward Kennedy J, Kereiakes D, Kern MJ, et al. ACC/AHA guidelines for percutaneous coronary intervention (revision of the 1993 PTCA guidelines)—executive summary. J Am Coll Cardiol. (2001) 37:2215–38. doi: 10.1016/S0735-1097(01)01344-4
58. Antoniadis AP, Papafaklis MI, Takahashi S, Shishido K, Andreou I, Chatzizisis YS, et al. Arterial remodeling and endothelial shear stress exhibit significant longitudinal heterogeneity along the length of coronary plaques. JACC-Cardiovasc Imag. (2016) 9:1007–9. doi: 10.1016/j.jcmg.2016.04.003
59. Dobrin PB. Mechanical properties of arteries. Physiol Rev. (1978) 58:397–460. doi: 10.1152/physrev.1978.58.2.397
60. Pao YC, Lu JT, Ritman EL. Bending and twisting of an in vivo coronary artery at a bifurcation. J Biomech. (1992) 25:287–95. doi: 10.1016/0021-9290(92)90026-W
Keywords: biomechanical assessment, finite element analysis, optical coherence tomography, plaque structural stress, plaque rupture
Citation: Huang J, Yang F, Gutiérrez-Chico JL, Xu T, Wu J, Wang L, Lv R, Lai Y, Liu X, Onuma Y, Tang D, Serruys PW, Wijns W and Tu S (2021) Optical Coherence Tomography-Derived Changes in Plaque Structural Stress Over the Cardiac Cycle: A New Method for Plaque Biomechanical Assessment. Front. Cardiovasc. Med. 8:715995. doi: 10.3389/fcvm.2021.715995
Received: 27 May 2021; Accepted: 04 October 2021;
Published: 04 November 2021.
Edited by:
Rene M. Botnar, King's College London, United KingdomReviewed by:
Kenichiro Otsuka, Massachusetts General Hospital, United StatesYao-Jun Zhang, Xuzhou Cancer Hospital, China
Copyright © 2021 Huang, Yang, Gutiérrez-Chico, Xu, Wu, Wang, Lv, Lai, Liu, Onuma, Tang, Serruys, Wijns and Tu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shengxian Tu, c3h0dSYjeDAwMDQwO3NqdHUuZWR1LmNu