Skip to main content

ORIGINAL RESEARCH article

Front. Bioeng. Biotechnol., 23 August 2021
Sec. Biomechanics

Using Optical Coherence Tomography and Intravascular Ultrasound Imaging to Quantify Coronary Plaque Cap Stress/Strain and Progression: A Follow-Up Study Using 3D Thin-Layer Models

Rui LvRui Lv1Akiko MaeharaAkiko Maehara2Mitsuaki MatsumuraMitsuaki Matsumura2Liang WangLiang Wang1Caining ZhangCaining Zhang1Mengde HuangMengde Huang1Xiaoya GuoXiaoya Guo3Habib SamadyHabib Samady4Don. P. Giddens,Don. P. Giddens4,5Jie ZhengJie Zheng6Gary S. MintzGary S. Mintz2Dalin Tang,
Dalin Tang1,7*
  • 1School of Biological Science and Medical Engineering, Southeast University, Nanjing, China
  • 2The Cardiovascular Research Foundation, Columbia University, New York, NY, United States
  • 3School of Science, Nanjing University of Posts and Telecommunications, Nanjing, China
  • 4Department of Medicine, Emory University School of Medicine, Atlanta, GA, United States
  • 5The Wallace H. Coulter Department of Biomedical Engineering, Georgia Institute of Technology, Atlanta, GA, United States
  • 6Mallinckrodt Institute of Radiology, Washington University, St. Louis, MO, United States
  • 7Mathematical Sciences Department, Worcester Polytechnic Institute, Worcester, MA, United States

Accurate plaque cap thickness quantification and cap stress/strain calculations are of fundamental importance for vulnerable plaque research. To overcome uncertainties due to intravascular ultrasound (IVUS) resolution limitation, IVUS and optical coherence tomography (OCT) coronary plaque image data were combined together to obtain accurate and reliable cap thickness data, stress/strain calculations, and reliable plaque progression predictions. IVUS, OCT, and angiography baseline and follow-up data were collected from nine patients (mean age: 69; m: 5) at Cardiovascular Research Foundation with informed consent obtained. IVUS and OCT slices were coregistered and merged to form IVUS + OCT (IO) slices. A total of 114 matched slices (IVUS and OCT, baseline and follow-up) were obtained, and 3D thin-layer models were constructed to obtain stress and strain values. A generalized linear mixed model (GLMM) and least squares support vector machine (LSSVM) method were used to predict cap thickness change using nine morphological and mechanical risk factors. Prediction accuracies by all combinations (511) of those predictors with both IVUS and IO data were compared to identify optimal predictor(s) with their best accuracies. For the nine patients, the average of minimum cap thickness from IVUS was 0.17 mm, which was 26.08% lower than that from IO data (average = 0.23 mm). Patient variations of the individual errors ranged from ‒58.11 to 20.37%. For maximum cap stress between IO and IVUS, patient variations of the individual errors ranged from ‒30.40 to 46.17%. Patient variations of the individual errors of maximum cap strain values ranged from ‒19.90 to 17.65%. For the GLMM method, the optimal combination predictor using IO data had AUC (area under the ROC curve) = 0.926 and highest accuracy = 90.8%, vs. AUC = 0.783 and accuracy = 74.6% using IVUS data. For the LSSVM method, the best combination predictor using IO data had AUC = 0.838 and accuracy = 75.7%, vs. AUC = 0.780 and accuracy = 69.6% using IVUS data. This preliminary study demonstrated improved plaque cap progression prediction accuracy using accurate cap thickness data from IO slices and the differences in cap thickness, stress/strain values, and prediction results between IVUS and IO data. Large-scale studies are needed to verify our findings.

Introduction

Cardiovascular disease is a serious threat to human life and health. Atherosclerotic plaques often rupture without warning, leading to acute cardiovascular syndrome. It is commonly agreed that plaques prone to rupture (also called thin-cap fibroatheroma, TCFA) normally have three main characteristics: large lipid-rich necrotic core, higher prevalence of macrophage infiltration in the fibrous cap, and a fibrous cap with thickness less than 65 μm (Kolodgie et al., 2001). Among the three, cap thickness may be one of the more closely monitored and critical factors for assessing plaque stress, strain, and vulnerability. Quantification of fibrous cap thickness and cap stress/strain conditions plays an important role in plaque progression prediction and vulnerability assessment, which in turn has clinical importance in monitoring disease development and patient management.

Considerable efforts have been made in investigating plaque progression using image-based patient follow-up data (Yang et al., 2010; Samady et al., 2011; Stone et al., 2011; Stone et al., 2012; Hung et al., 2016; Wang et al., 2018). Plaque vessel wall thickness change was often used as a measure for plaque progression. In most available patient follow-up data (follow-up time around 1 year), plaque vessel wall thickness changes were mostly under 100 µm in a year (Yang et al., 2010). With intravascular ultrasound (IVUS) imaging resolution at 150–200 µm and magnetic resonance imaging (MRI) resolution at 200–300 μm, it is extremely challenging to quantify and predict plaque progression with acceptable reliability and accuracy. With its superior resolution at about 10 μm, optical coherence tomography (OCT) is getting acceptance in research community and clinical practices. OCT is able to detect the thin fibrous cap of vulnerable plaques (Kume et al., 2006). However, due to its limited penetration, OCT is not able to “see” the whole vessel. Several groups have combined IVUS and OCT in detecting and analyzing vulnerable plaques, and impressive results have been reported (Sawada et al., 2008; Fujii et al., 2015; Zanchin et al., 2020). As an initial effort, a modelling approach combining IVUS and OCT was introduced by Guo et al. (2019) for cap thickness quantification, more accurate cap stress/strain calculations, and plaque progression prediction. With IVUS and OCT combined, the merged IVUS + OCT (IO) data include the whole vessel from IVUS and accurate cap thickness from OCT. The accurate and reliable image data provide a solid base for further vulnerable plaque investigations. In particular, we should be able to trust cap thickness change measured by OCT and use that to obtain some reliable progression prediction results.

The goal of this study is to extend the initial effort by Guo et al. (2019) to a multipatient study. IVUS, OCT, and angiography data from nine patients with follow-up were used to generate IVUS and IVUS + OCT (IO) slices, where IO slices were made by combining IVUS and OCT segmented contours together. Since cap thickness is one of the most critical and most closely monitored measurable plaque features linked to plaque vulnerability, it was selected as the target for our prediction study. A total of 114 matched (baseline and follow-up) slices were obtained to quantify their cap thickness differences between IVUS and OCT data and their changes between baseline and follow-up. 3D thin-layer models were constructed using both IVUS and IO slices to obtain cap stress and strain data. Two prediction methods (a generalized linear mixed model (GLMM) and least squares support vector machine (LSSVM) method) were used to predict cap thickness changes using nine morphological and mechanical risk factors. Prediction accuracies by all combinations (total 511) of those predictors with both IVUS and IO data were compared to identify the optimal predictor(s) with their best accuracies.

Materials and Methods

Data Acquisition and Processing

IVUS, OCT, and angiography data including three epicardial coronary vessels were collected at Cardiovascular Research Foundation (CRF) from patients with coronary heart diseases between April 2017 and November 2018 using the protocol approved by the local institute, and informed consent forms were obtained from the patients. Nine lipid-rich plaques were identified from nine patients with good-quality baseline and follow-up data matched for our analysis (average follow-up days: 251). These patients had stable angina pectoris and underwent percutaneous coronary intervention (PCI). Patients with acute coronary syndrome, severe calcified lesion, chronic total occlusion, or chronic kidney disease (Cr > 1.5 mg/ dl) were excluded. Patient demographic data are shown in Table 1. IVUS examination was performed after 0.2 mg of intracoronary nitroglycerin. The IVUS catheter was advanced as far as possible using a commercially available IVUS system: a 40 MHz IVUS catheter (OptiCross, Boston Scientific Corporation, Natick, Massachusetts) with motorized pullback at 0.5 mm/ s. OCT images were acquired with ILUMIEN OPTIS System and Dragonfly or Dragonfly JP Imaging Catheter (St. Jude Medical, Westford, Massachusetts).

TABLE 1
www.frontiersin.org

TABLE 1. Patient demographic and clinical information. F: female; M: male; BL: baseline; FU: follow-up; HT: hypertension; DM: diabetes mellitus; HL: hyperlipidemia.

Coregistration of IVUS and OCT slices was performed by using fiduciary points such as side branches, bifurcations, and calcifications with the assistance of quantitative coronary angiography (QCA) images. IVUS and OCT frames at baseline and follow-up were also matched to quantify changes of plaque morphology and cap thickness using the same matching techniques. Three plaque compositions were considered in segmentation of IVUS and OCT images: lipid-rich necrotic core (short form, lipid) and calcification and other vessel tissue (fibrotic, fibro-fatty, etc.). Segmentation was performed by ImageJ 1.52 v software. For IVUS images, we used manual segmentation. For OCT images, we used the combination of manual segmentation and algorithms. Details of the segmentation methods were previously published (Lv et al., 2020). Small-size plaque components were neglected for simplicity. Matched IVUS and OCT slices were combined together to form IVUS + OCT (IO) slices with OCT providing accurate cap thickness quantifications, while IVUS providing vessel outboundary contours not visible from OCT (Lv et al., 2020). Sample IO and IVUS slices with segmented contours are given in Figure 1 showing the merging process.

FIGURE 1
www.frontiersin.org

FIGURE 1. Sample IVUS and merged IVUS + OCT (IO) slices showing segmentation and merging at baseline and follow-up. (A) IVUS slice and its segmented contours at baseline. (B) Corresponding OCT slice and IVUS + OCT (IO) slice at baseline. (C) IVUS slice and its segmented contour plot at follow-up. (D) Corresponding OCT slice and IVUS + OCT (IO) contour plot at follow-up. Outwall contours of IO slices are taken from IVUS.

IVUS and corresponding IO slices at baseline and follow-up were matched, and 114 matched lipid-slices were selected for analysis. Each slice (IVUS and IO) was divided into 100 parts using 100 evenly spaced lumen points and a Four-Quarter Even-Spacing method so that vessel and cap thickness could be calculated properly and recorded for analysis (Wang et al., 2019). Figure 2 illustrates our Four-Quarter Even-Spacing method with which cap thickness was defined and calculated. It also shows clear cap thickness difference between IVUS and IO slices and between baseline and follow-up. Nine morphological and mechanical risk factors were selected as our predictors to predict cap thickness change from baseline to follow-up: lumen area (LA), plaque area (PA) which is defined as the area between the vessel outwall (see Figure 2) and the lumen contour, plaque burden (PB) defined as the ratio of the plaque area over the area inside the vessel outwall, minimum cap thickness (MinCapT), mean cap thickness (MeanCapT), maximum cap stress (MaxCapS), mean cap stress (MeanCapS), maximum cap strain (MaxCapSn), and mean cap strain (MeanCapSn). Their values were calculated from baseline IVUS and IO slices. Values of cap-related predictors were collected from related cap lumen points (see Figure 2). For better clarity, the distance between the lumen point and its corresponding wall point was defined as the wall thickness. The connecting line from the lumen contour to the edge of the lipid near the lumen was defined as cap thickness (CapT). For a given slice, MinCapT was defined as the minimum of all the fibrous cap point thicknesses for the slice. When analysis was performed at patient level, MinCapT was defined as the minimum of all fibrous cap point thicknesses in all lipid slices of the patient. Slice and patient MeanCapT were defined the same way. Maximum and mean values of stress and strain predictors were calculated similarly. Data for all predictors were collected, calculated, and saved for prediction model use. The work of morphology calculation was accomplished by MATLAB (MATLAB R2018a, MathWorks, United States).

FIGURE 2
www.frontiersin.org

FIGURE 2. The Four-Quarter Even-Spacing method with which cap thickness was defined and calculated. Cap thickness differences between IVUS and IO slices and between baseline and follow-up are shown. (A) IVUS segmented contours with lumen–outwall connection lines at baseline. (B) Corresponding IVUS + OCT (IO) segmented contours with lumen–outwall connection lines at baseline. (C) IVUS slice with connection lines at follow-up. (D) Corresponding IVUS + OCT (IO) slice with connection lines at follow-up. Minimum cap thickness was marked using bold black.

The 3D Thin-Layer Biomechanical Model

For the 114 matched IVUS and IO slices containing a lipid-rich core with fibrous cap, a thin-slice thickness (0.5 mm) was used to make the 3D thin-layer models (see Figure 3). Since IVUS and OCT data were acquired under in vivo condition when the vessel was under stretched and loaded conditions, a 5% axial shrink–stretch combined with our circumferential preshrink process was performed to obtain in vivo slice morphology with proper plaque stress/strain when 5% axial stretch and pressure condition were imposed (Wang et al., 2019). Details of the preshrink procedure are described in our previous studies (Yang et al., 2010; Wang et al., 2015). Vessel tissue was assumed to be hyperelastic, anisotropic, nearly incompressible, and homogeneous (Guo et al., 2018). Plaque components (lipid core and calcification) were assumed to be hyperelasic, isotropic, and nearly incompressible. A modified Mooney–Rivlin material model was used to describe the material properties of vessel tissue with the strain-energy density function given as follows (Wang et al., 2018):

Wiso=c1(I13)+c2(I23)+D1[exp(D2(I13))1],(1)
Waniso=Wiso+K1K2{exp[K2(I41)2]1},(2)

where I1=(Cii) and I2=12(I12CijCij), in which  I1 and I2 are the first and second invariants of right Cauchy–Green deformation tensor C=(Cij)=FTF=(Fii)=(xi/aj), (xi) is the current position, (aj)is the original position, and I4=Cij(nc)i(nc)j, where nc is the unit vector in the circumferential direction of the vessel and c1,  c2, D1, D2, K1, and K2 are the material parameters. Material constants of the isotropic Mooney–Rivlin model from the existing literature were used. Lipid: c1 = 0.5 kPa, c2 = 0 kPa, D1 = 0.5 kPa, and D2 = 1.5; calcification: c1 = 92 kPa, c2 = 0 kPa, D1 = 36 kPa, and D2 = 2.0; vessel/fibrous tissue: c1 = ‒262.76 kPa, c2 = 22.9 kPa, D1 = 125.9 kPa, D2 = 2.0, K1 = 7.19 kPa, and K2 = 23.5 (Guo et al., 2018; Guo et al., 2019).

FIGURE 3
www.frontiersin.org

FIGURE 3. The 3D thin-layer model construction process. (A) An in vivo IVUS + OCT (IO) slice. (B) The slice after shrinking. (C) The 3D thin-layer model with a thickness added. (D) The slice with load-added matching in vivo slice morphology. (E) Plaque wall stress under minimum pressure. (F) Plaque wall strain under minimum pressure.

The 3D thin-layer models were solved by finite element software ADINA 9.0 (Adina R & D, Watertown, MA, United States) following our established procedures (Wang et al., 2019). Nonlinear incremental iterative procedures were used to solve the models. Mesh analysis was performed by refining mesh density by 10% until changes of solutions became less than 2%. Figure 4 provides plots of a sample slice contour with cap thickness and stress and strain distributions to show the cap thickness and stress/strain differences between IVUS and IO models. Cap thickness, plaque stress, plaque strain, plaque cap stress, and cap strain were extracted and saved for analysis.

FIGURE 4
www.frontiersin.org

FIGURE 4. An example showing the differences of cap thickness, maximum cap stress, and maximum cap strain in the cap region and maximum stress and strain on the whole lumen for IVUS + OCT (IO) and IVUS. (AC) IO slice—contour, stress, and strain plots. (DF) IVUS slice—contour, stress, and strain plots.

Prediction Models and Strategy

Nine morphological and mechanical risk factors were selected as our predictors to predict cap thickness change from baseline to follow-up (see Data Acquisition and Processing for descriptions): LA, PA, PB, MinCapT, MeanCapT, maximum cap stress (MaxCapS), mean cap stress (MeanCapS), maximum cap strain (MaxCapSn), and mean cap strain (MeanCapSn). As the prediction target variable, Δ(MinCapT) is defined as

Δ(MinCapT)=ΔMinCapT=MinCapTatT2MinCapTatT1,(3)

where T1 = baseline and T2 = follow-up. Δ(MinCapT) of IO was used as the benchmark for model prediction verifications. For a given slice, its binary labels for ΔMinCapT were defined as

SliceLabel(ΔMinCapT)={1,ifΔMinCapT<0,1,ifΔMinCapT0,(4)

where Slice Label = 1 means the lipid plaque is progressive (getting worse), while Slice Label = ‒1 means the lipid plaque is nonprogressive.

A 5-fold cross-validation strategy with generalized linear mixed model (GLMM) and least squares support vector machine (LSSVM) methods were adopted to calculate prediction accuracy for each single predictor and all combinations of the nine predictors and quantify prediction accuracy differences using IO and IVUS slice values (Wang et al., 2015). The GLMM model is defined as

yij=E(yij|bj)+εij,(5)
logit(E(yij|bj))=β0+β1x1ij+β2x2ij++bj,(6)

where yij is the binary response of ΔMinCapT and E(yij|bj)=P(yij=1|bj) is the probability on the ith slice of the jth patient. logit(x)=log[x/(1x)] is the binomial link function. x1ij, x2ij, etc. are the risk factors, which were used to predict plaque cap progression. β0, β1, etc. are the fixed-effect coefficients, and bj and εij are the random-effect terms and the random error terms of GLMM. R function glmmPQL was used to estimate the term values by fitting GLMM.

The LSSVM method used Gaussian radial basis function as kernel function, fitted all the samples with least square error, and used the steepest descent method to find the optimal parameters. Radial basis function is given by

K(x,z)=exp(xz2σ2),(7)

where σ2 is the square-coordinated kernel width. MATLAB (MATLAB R2018a, MathWorks, United States) was used to adjust the parameters and complete the prediction work.

All 114 IO and IVUS slices were randomly divided into 5 subgroups, 4 as training subgroups and 1 as test subgroup. The training subgroups were used for model fitting. The test subgroup was used to test the model, and dichotomous prediction is performed. Each subgroup was used once, and the testing results from the five subgroups were combined together as the final result to give prediction accuracy. The binary responses for the measures of cap thickness changes (ΔMinCapT) were adopted to find the best predictor(s). If fibrous cap thickness decreases, it was positive (labeled 1); meanwhile, increase was negative (labeled ‒1). The classified result is judged by accuracy (Acc), sensitivity (Sen), specificity (Spe), and area under the ROC curve (AUC). Accuracy is the probability of correctly predicted slices being detected in all test groups. Sensitivity is the probability of detecting a plaque progression slice in the plaque progression test subgroup. Specificity is the probability of the correct slice being detected in the nonplaque progression test subgroup. The AUC was used to assess the predictive power of each predictor.

Statistical Analysis

MinCapT, MeanCapT, MaxCapS, MeanCapS, MaxCapSn, and MeanCapSn from IO and IVUS models were compared to get differences. Errors by IVUS data and models were calculated using IO results as baseline (gold standard). The Kolmogorov–Smirnov tests were performed to check data normality. Analysis of variance, paired t-tests, and Mann–Whitney U test were used to check if the differences were statistically significant. A p value < 0.05 was considered statistically significant.

Results

IVUS Cap Thickness Data Had Large Errors Compared to OCT Data

Table 2 lists minimum and mean cap thickness values from nine patients based on 114 slices at baseline. Mean cap thickness (MeanCapT) was the average of point-to-point cap thickness for every patient. The average of MinCapT from IVUS was 0.17 mm, which was 26.09% lower than that from IO data (average = 0.23 mm). Patient variations of the individual errors ranged from ‒58.11 to 20.37%, with mean ± SD = ‒20.70 ± 30.32%. Negative error value means IVUS cap thickness values were lower than those from IO values. The average of MeanCapT from IVUS was 0.38 mm, which was 2.56% lower than that from IO data. Patient variations of the individual errors ranged from ‒30.32 to 25.50%, with mean ± SD = 0.16 ± 20.86%. Patient five had the biggest error (in absolute value) for both MinCapT (‒58.11%) and MeanCapT (‒30.32%). Patient six had maximum MinCapT (0.324 mm) based on IO slice, but its IVUS slice had MinCapT 0.141 mm with error ‒56.48%.

TABLE 2
www.frontiersin.org

TABLE 2. Summary of minimum and mean cap thickness from 9 patients based on 114 slices at baseline. Unit for fibrous cap thickness: mm.

Cap Stress and Strain From IVUS and IO Data Have Large Differences and the Differences Have Large Patient Variations

Table 3 summarizes mean and maximum cap stress values from the nine patients based on simulation results from 114 3D thin-layer models. The average of IO MeanCapS of the 114 slices was 67.15 kPa. The average MeanCapS from IVUS slices was 73.44 kPa, which was 9.37% higher than that from IO models. The modest error value could be misleading since positive and negative individual errors canceled against each other in the average. Patient variations of the individual errors for MeanCapS ranged from ‒37.68 to 44.25%, with mean ± SD = 13.93 ± 29.69%. The average of MaxCapS of the 114 IO slices was 84.36 kPa, compared to 88.23 kPa from IVUS slices, which was 4.59% higher than that from IO models. Patient variations of the individual errors ranged from ‒30.40 to 46.17%, with mean ± SD = 7.24 ± 26.47%. Patientwise comparisons of IO and IVUS cap stress indicates that Patient four had the highest error for MeanCapS and MaxCapS (44.25 and 46.17%, respectively).

TABLE 3
www.frontiersin.org

TABLE 3. Summary of mean and maximum cap stress values for the plaques of 9 patients. Unit for stress: kPa.

Table 4 summarizes mean and maximum cap strain values from the nine patients from our 114 3D thin-layer models. The average of IO MeanCapSn of the 114 slices was 0.15 which is the same as that from IVUS slices. However, patient variations of the individual errors ranged from ‒22.92 to 29.63%, with mean ± SD = 4.33 ± 16.10%. The average of MaxCapSn of the 114 IO slices was 0.17, which again was the same as that from the IVUS slices. Patient variations of the individual errors ranged from ‒19.90 to 17.65%, with mean ± SD = 0.11 ± 13.94%. Patientwise comparisons of IO and IVUS mean cap strain indicated that Patient nine had the highest error in 29.63% under MeanCapSn and Patient seven had the lowest error in ‒19.90% under MaxCapSn. Figure 4 gives an example to show the differences of mean cap thickness, maximum cap stress, and maximum cap strain in the cap region and maximum stress and strain on the whole lumen for IO and IVUS.

TABLE 4
www.frontiersin.org

TABLE 4. Summary of mean and maximum cap strain values for 9 patients.

Cap Thickness Progression Prediction Using GLMM: Comparison Between IO and IVUS Results

There were altogether 511 (2^9–1) combinations using nine predictors. Table 5 gives prediction results by two best combinations and nine single predictors for ΔMinCapT using IO and IVUS data based on the GLMM method, respectively. Values for prediction accuracy (Acc), specificity (Spe), and area under the ROC curve (AUC) are given. The optimal combination predictor among all 511 predictors using IO data was the combination of MinCapT, MeanCapS, PB, and PA with the highest AUC = 0.926 and highest accuracy = 90.8% for IO data, which compared favourably with AUC = 0.783 and Acc = 74.6% using IVUS data. The best single predictor was PA with the highest AUC value in IO prediction results with AUC = 0.868 and accuracy = 85%, compared to AUC = 0.637 and accuracy = 64.0% using IVUS data. The best prediction accuracy of IO using combinations was 16.2% higher than that of IVUS (90.8 vs. 74.6%). The AUCs for the best combinations of predictors using IVUS and IO data are given in Figure 5. For the predictor using IO data, the best combination is given by PA + MinCapT + MeanCapS + PB. For IVUS data, the best combination is given by PA + MinCapT + MeanCapT + MeanCapS + MeanCapSn.

TABLE 5
www.frontiersin.org

TABLE 5. ΔMinCapT prediction results using the GLMM method. AUC: area under the ROC curve; Acc: accuracy; Sen: sensitivity; Spe: specificity. Bold indicates the best value.

FIGURE 5
www.frontiersin.org

FIGURE 5. ROC curves of the best combinations of predictors predicting ΔMinCapT. (A) ROC curve using IO data. The best combination is given by PA + MinCapT + MeanCapS + PB. (B) ROC curve using IVUS data. The best combination is given by PA + MinCapT + MeanCapT + MeanCapS + MeanCapSn.

Cap Thickness Progression Prediction Using LSSVM: Comparison Between IO and IVUS Results

Using the LSSVM method to predict ΔMinCapT, Table 6 gives prediction results by two best combinations and nine single predictors using IO and IVUS data, respectively. The optimal combination predictor (best AUC value) among all 511 predictors using IO data was the combination of PA, PB, MinCapT, and MeanCapS with AUC = 0.838 and accuracy = 75.7%. The best predictor for IVUS data was a combination of MeanCapT, MaxCapS, and MeanCapSn with AUC = 0.780 and accuracy = 69.6%. The best single predictor was MinCapT in IO prediction results with AUC = 0.804 and accuracy = 76%, compared to AUC = 0.725 and accuracy = 71.5% using IVUS data. The best prediction accuracy of IO using combinations was 6.1% higher than that of IVUS (75.7 vs. 69.6%).

TABLE 6
www.frontiersin.org

TABLE 6. ΔMinCapT prediction results using the LSSVM method. AUC: area under the ROC curve; Acc: accuracy; Sen: sensitivity; Spe: specificity. Bold indicates the best value.

Combining Morphological and Mechanical Factors Improved Prediction Accuracies

Figure 6 provides bar plots comparing prediction results using only morphological predictors and using combinations of morphological and mechanical predictors. Results from GLMM and LSSVM methods using IO data are shown. For the GLMM method, combining biomechanical and morphological factors helped to improve prediction accuracy from 88.8 to 90.8%. The AUC value improved from 0.911 to 0.926. Sensitivity did not show improvement. Specificity improved from 0.9 to 0.925. The LSSVM method showed similar improvements, while overall performance was not as good as that from the GLMM method.

FIGURE 6
www.frontiersin.org

FIGURE 6. Comparison between the best morphology–only combined predictor and the best biomechanical + morphology combined predictor using IO data based on ΔMinCapT prediction results. (A) IO prediction result based on GLMM. (B) IO prediction result based on LSSVM. AUC: area under the ROC curve; Acc: accuracy; Sen: sensitivity; Spe: specificity.

Discussion

Significance of Accurate Cap Thickness Data for Plaque Progression Assessment

In order to predict plaque progression, one simple truth is that we must be able to acquire the progression data with enough accuracy to predict. With IVUS imaging resolution at 150–200 μm and MRI resolution at 200–300 μm, and with the fact that plaque vessel wall thickness changes were mostly under 100 μm in a year, well below imaging resolutions (Yang et al., 2010), it is extremely challenging to quantify and predict plaque progression with acceptable reliability and accuracy. OCT comes with a superior 10–15 μm resolution, and it gives us hope that we may get accurate data to work with. However, OCT is limited by its penetration and cannot get the whole vessel cross section. Combining IVUS and OCT provides one approach where we have the whole cross section from IVUS and accurate cap thickness from OCT. Since the whole cross section is primarily from IVUS with its resolution limitation, it remains to be questionable to use vessel wall thickness as the measure for plaque progression. In this study, cap thickness was selected as the measure for plaque progression. While it does not give complete information for plaque changes, it does come with the accuracy from OCT and it is one of the most critical plaque features affecting plaque stress/strain conditions. The well-known cap thickness threshold value (65 μm) of vulnerable plaques is usually used as an important criterion for morphological classification (Virmani et al., 2002, Virmani et al., 2006; Cardoso and Weinbaum, 2014; Fenning and Wilensky, 2014; Saba et al., 2014; Arbab-Zadeh and Fuster, 2015; Goncalves et al., 2015; Bala et al., 2018; Yonetsu and Jang, 2018; Zhao et al., 2018; Porcu et al., 2020). For the first time, the progression prediction models were based on image data with resolution good enough to quantify thin cap thickness and the predicted target parameter (cap thickness) has values that are accurate and reliable. With that, prediction accuracies should be more reliable and meaningful.

Cap thickness data from OCT (and IO slices) are considered and used as the gold standard. Therefore, comparison of IVUS and IO data and prediction results is to show the errors (uncertainties) by IVUS images when we have “true” cap thickness obtained by OCT. It should be noted that the average of MeanCapT from IVUS was only 2.56% lower than that from IO data, while patient variations of the individual errors ranged from ‒30.32 to 25.50%, with mean ± SD = 0.16 ± 20.86%. The average of MinCapT from IVUS was 0.17 mm, which was 26.09% lower than that from IO data (average = 0.23 mm). Patient variations of the individual MinCapT errors were even worse: it ranged from ‒58.11 to 20.37%, with mean ± SD = ‒20.70 ± 30.32%. These data indicated that cap thickness data from IVUS have high uncertainties and should be interpreted and used with caution.

Comparison of Prediction Results by GLMM and LSSVM Methods and by IVUS and IO Data

While image data accuracy is the base of any prediction research, selection of prediction methods is also of great interest. The use of prediction models requires their stability, simplicity, extendibility, and good adaptability. The logistic regression model is one of the commonly used models for plaque prediction (Samady et al., 2011; Stone et al., 2012). GLMM is commonly used in clinical medicine, which can be seen as an extended form of GLM without the requirement to meet the normal distribution of dependent variable, and has the characteristics of GLM method at the same time including fixed effect and random effect. Because risk factors were extracted from the slices, the quantitative factors between adjacent slices were neither random nor independent of each other, and GLMM could take this into account. Artificial intelligence (AI) and various machine leaning methods have been introduced to perform predictions with limited success (Han et al., 2020; Kay et al., 2020). For this work, 114 slice data constitute a relatively small data set in general machine learning method applications. SVM is a machine learning method suitable for small samples. LSSVM is a support vector machine optimized by the least squares method. The function maps the space where the samples are mapped to a higher dimensional space, which helps to recognize the targeted group/classification (Guo et al., 2019). Both GLMM and LSSVM methods have the potential to further improve the accuracy of predictions. In this study, the quantitative index based on fibrous cap thickness is more direct and simple. Based on R language, the GLMM method has fast debugging and good stability. The prediction result of GLMM’s optimal combination predictor was better than that of LSSVM (AUC 0.926 vs. 0.830).

Prediction accuracy using IO data (AUC = 0.926 and accuracy = 90.8%, by GLMM) is clearly better that using IVUS data (AUC = 0.783 and accuracy = 74.6%, by GLMM). The comparison seems to be unfair since the prediction target ΔMinCapT was quantified using IO data. Our purpose is not to show that IO data are better since it was used as the gold standard in this study. The purpose of the comparison is to show the impact of IVUS data uncertainty on prediction accuracy. The differences of prediction accuracies using LSSVM were smaller, as shown in Table 6 (AUC: 0.838 vs. 0.780; accuracy: 83 vs. 69.6%).

Cap Stress and Strain Differences Between IVUS and IO Models

Mechanical forces, especially cap stress and strain, play an important role in plaque progression and vulnerability assessment. Cap thickness differences between IVUS and IO data led to large cap stress/strain differences. It is worth noting that overall averaged cap stress and strain values had very modest differences between IO and IVUS data, while individual differences ranged from about −40% to + 40%. Patient variations of the individual errors for mean cap stress (MeanCapS) ranged from ‒37.68 to 44.25%, with mean ± SD = 13.93 ± 29.69%. For maximum cap stress which could be closely related to potential plaque rupture, patient variations of the individual errors ranged from ‒30.40 to 46.17%, with mean ± SD = 7.24 ± 26.47%. For cap strain comparisons, patient variations of individual errors of mean cap strain (MeanCapSn) ranged from ‒22.92 to 29.63%, with mean ± SD = 4.33 ± 16.10%. Patient variations of the individual errors for MinCapSn ranged from ‒19.90 to 17.65%, with mean ± SD = 0.11 ± 13.94%. Errors for strain comparisons were smaller than those from stress comparisons.

Combining Biomechanical and Morphological Factors Could Improve Prediction Accuracy

Combination of morphology and mechanical predictors may give better prediction accuracy than using morphology only or single-factor predictors which was seen again in this work. Using the GLMM method, there is an improvement in accuracy from 88.8 to 90.8%, in Sen + Spe from 0.751 to 1.776 and in AUC value from 0.911 to 0.926 based on IO data (Figure 6A). The corresponding IVUS combined predictors shows the same trend with MaxCapSn, which helps the combination of LA and MeanCapT increase the AUC from 0.756 to 0.778. Figure 6B illustrates a similar result using the LSSVM method. Combining morphological factors and mechanical factors demonstrated a certain ability not to be ignored in plaque progression prediction.

Further Plaque Prediction Strategies

Plaque progression is a long and complex process, and it involves many factors. Data from multichannels should be collected. And, ideally, more time points (i.e., follow-up 1, follow-up 2, etc.) would be desirable and may give better results. However, those patient studies are very expensive, especially for this IVUS + OCT studies. A one-year follow-up is a typical timeframe researchers use for plaque progression investigations (Samady et al., 2011; Stone et al., 2012; D’Ascenzo et al., 2021). OCT has an accuracy of 10–15 µm and can detect relatively small changes in cap thickness. Compared to other publications for plaque progression research, we tend to believe that better accuracy from OCT should lead to more accurate and reliable predictions.

Limitations and Future Directions

Automation of simulation modelling has always been a goal. A 3D thin-layer model helps move forward this process. However, there is still a need for further improvement. Patient-specific vessel material properties were not available for this study, and parameters from our previous studies using VH-IVUS and Cine-IVUS were used (Guo et al., 2019). 3D thin-layer models did not take into account the effects of fluid. The modest improvement of combining morphological and mechanical factors could also be due to the fact that flow shear stress is not included in the predictor list. Biochemical indicators such as hyperlipidemia or diabetes mellitus were not considered. Further refinement of the prediction models is needed.

One major limitation of this study is its small sample size. Clinical significance of our results in predicting plaque progression and regression is limited by the small patient size. Large-scale studies are needed to further verify our findings and bring our research scheme closer to clinical practice. It should be noted that data needed for this study are extremely hard to get because the data set included IVUS and OCT at both baseline and follow-up. Not only IVUS and OCT are very expensive, but also they are invasive procedures. While it is relatively easier to get patients to participate at baseline since they need to do that for their diagnosis anyway, it is much harder to recruit patients to do those invasive imaging procedures (IVUS and OCT) mainly for research purpose. It involves additional risk with only limited benefits to patients. Most insurances do not have coverage for such procedures. So, financial support becomes a limiting factor. These are practical difficulties in the data acquisition front. We will keep trying to get more data, improve our models, and move closer to clinical implementations.

Another major limitation is that predictors in this study only included morphological and biomechanical factors from structure-only models. Atherosclerosis is a pathophysiological process. While plaque rupture may be the result of the mechanical failure of the fibrous cap, there is a cascade of mechanisms involved in the thinning of the cap or the prevention of such thinning that are critical to monitor, so critical that is what is currently monitored in the clinical setting. Additional factors such as plaque erosion, inflammation, blood conditions, genomic factors, and many others should be explored to improve our models and prediction accuracies.

With the same context about vulnerable plaque research in biomechanical setting and current clinical practices, plaque stenosis severity remains to be the commonly quantified, monitored, and used in clinical diagnosis and treatment decision-making (stenting or conservative medication). Researchers have tried in recent years seeking better risk indicators including plaque components (therefore cap thickness) and plaque stress/strain conditions using more advanced imaging modalities and mechanical models. Cap thickness does play an important role in plaque stress/strain calculations. Yet, there are still considerable gaps between research findings and actual clinical implementations.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary material; further inquiries can be directed to the corresponding author.

Ethics Statement

The studies involving human participants were reviewed and approved by the Cardiovascular Research Foundation (CRF) using the protocol approved by the local institute IRB and informed consent forms were obtained from the patients. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author Contributions

RL made major contribution to the writing of the article, model building, data extraction, and prediction analysis; DT was responsible for overall design and execution of the research and a main writer of the manuscript; AM, GM, and MM worked on patient image data preparation; LW, CZ, and XG assisted in and contributed to data segmentation; CZ and MH assisted in and contributed to 3D thin-layer models; HS, DG, JZ, and GM provided clinical guidance and technical support in data processing and analysis.

Funding

This research was supported in part by the National Science Foundation of China under grants 11972117 and 11672001.

Conflict of Interest

The 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.

References

Arbab-Zadeh, A., and Fuster, V. (2015). The Myth of the "Vulnerable Plaque". J. Am. Coll. Cardiol. 65, 846–855. doi:10.1016/j.jacc.2014.11.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Bala, G., Broisat, A., Lahoutte, T., and Hernot, S. (2018). Translating Molecular Imaging of the Vulnerable Plaque-A Vulnerable Project? Mol. Imaging Biol. 20, 337–339. doi:10.1007/s11307-017-1147-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Berry, C., L’Allier, P. L., Grégoire, J., Lespérance, J., Levesque, S., Ibrahim, R., et al. (2007). Comparison of Intravascular Ultrasound and Quantitative Coronary Angiography for the Assessment of Coronary Artery Disease Progression. Circulation 115, 1851–1857. doi:10.1161/CIRCULATIONAHA.106.655654

PubMed Abstract | CrossRef Full Text | Google Scholar

Cardoso, L., and Weinbaum, S. (2014). Changing Views of the Biomechanics of Vulnerable Plaque Rupture: a Review. Ann. Biomed. Eng. 42, 415–431. doi:10.1007/s10439-013-0855-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Chau, A. H., Chan, R. C., Shishkov, M., MacNeill, B., Iftimia, N., Tearney, G. J., et al. (2004). Mechanical Analysis of Atherosclerotic Plaques Based on Optical Coherence Tomography. Ann. Biomed. Eng. 32, 1494–1503. doi:10.1114/b:abme.0000049034.75368.4a

PubMed Abstract | CrossRef Full Text | Google Scholar

Corban, M. T., Eshtehardi, P., Suo, J., McDaniel, M. C., Timmins, L. H., Rassoul-Arzrumly, E., et al. (2014). Combination of Plaque burden, wall Shear Stress, and Plaque Phenotype Has Incremental Value for Prediction of Coronary Atherosclerotic Plaque Progression and Vulnerability. Atherosclerosis 232, 271–276. doi:10.1016/j.atherosclerosis.2013.11.049

PubMed Abstract | CrossRef Full Text | Google Scholar

Costopoulos, C., Timmins, L. H., Huang, Y., Hung, O. Y., Molony, D. S., Brown, A. J., et al. (2019). Impact of Combined Plaque Structural Stress and wall Shear Stress on Coronary Plaque Progression, Regression, and Changes in Composition. Eur. Heart J. 40, 1411–1422. doi:10.1093/eurheartj/ehz132

PubMed Abstract | CrossRef Full Text | Google Scholar

D'Ascenzo, F., Filippo, O. D., Gallone, G., Mittone, G., Deriu, M. A., Iannaccone, M., et al. (2021). Machine Learning-Based Prediction of Adverse Events Following an Acute Coronary Syndrome (PRAISE): a Modelling Study of Pooled Datasets. Lancet 16 (397), 199–207. doi:10.1016/S0140-6736(20)32519-8

CrossRef Full Text | Google Scholar

Diletti, R., Garcia-Garcia, H. M., Gomez-Lara, J., Brugaletta, S., Wykrzykowska, J. J., van Ditzhuijzen, N., et al. (2011). Assessment of Coronary Atherosclerosis Progression and Regression at Bifurcations Using Combined IVUS and OCT. JACC Cardiovasc. Imaging 4, 774–780. doi:10.1016/j.jcmg.2011.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Erlinge, D., Maehara, A., Ben-Yehuda, Ori., Bøtker, H. Erik., Maeng, M., Kjøller-Hansen, L., et al. (2021). Identification of Vulnerable Plaques and Patients by Intracoronary Near-Infrared Spectroscopy and Ultrasound (PROSPECT II): a Prospective Natural History Study. Lancet 397, 985–995. doi:10.1016/S0140-6736(21)00249-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Fenning, R. S., and Wilensky, R. L. (2014). New Insights into the Vulnerable Plaque from Imaging Studies. Curr. Atheroscler. Rep. 16, 397. doi:10.1007/s11883-014-0397-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Fujii, K., Hao, H., Shibuya, M., Imanaka, T., Fukunaga, M., Miki, K., et al. (2015). Accuracy of OCT, Grayscale IVUS, and Their Combination for the Diagnosis of Coronary TCFA. JACC Cardiovasc. Imaging 8, 451–460. doi:10.1016/j.jcmg.2014.10.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Gonçalves, I., den Ruijter, H., Nahrendorf, M., and Pasterkamp, G. (2015). Detecting the Vulnerable Plaque in Patients. J. Intern. Med. 278, 520–530. doi:10.1111/joim.12414

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, X., Giddens, D. P., Molony, D., Yang, C., Samady, H., Zheng, J., et al. (2019). A Multimodality Image-Based Fluid-Structure Interaction Modeling Approach for Prediction of Coronary Plaque Progression Using IVUS and Optical Coherence Tomography Data with Follow-Up. J. Biomech. Eng. 141, 0910031–0910039. doi:10.1115/1.4043866

CrossRef Full Text | Google Scholar

Guo, X., Giddens, D. P., Molony, D., Yang, C., Samady, H., Zheng, J., et al. (2018). Combining IVUS and Optical Coherence Tomography for More Accurate Coronary Cap Thickness Quantification and Stress/Strain Calculations: A Patient-Specific Three-Dimensional Fluid-Structure Interaction Modeling Approach. J. Biomech. Eng. 140, 041005. doi:10.1115/1.4038263

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, D., Kolli, K. K., Al'Aref, S. J., Baskaran, L., van Rosendael, A. R., Gransar, H., et al. (2020). Machine Learning Framework to Identify Individuals at Risk of Rapid Progression of Coronary Atherosclerosis: From the PARADIGM Registry. J. Am. Heart Assoc. 9, e013958. doi:10.1161/JAHA.119.013958

PubMed Abstract | CrossRef Full Text | Google Scholar

Hung, O. Y., Molony, D., Corban, M. T., Rasoul‐Arzrumly, E., Maynard, C., Eshtehardi, P., et al. (2016). Comprehensive Assessment of Coronary Plaque Progression with Advanced Intravascular Imaging, Physiological Measures, and Wall Shear Stress: A Pilot Double‐Blinded Randomized Controlled Clinical Trial of Nebivolol Versus Atenolol in Nonobstructive Coronary Artery Disease. J. Am. Heart Assoc. 5, e002764. doi:10.1161/JAHA.115.002764

PubMed Abstract | CrossRef Full Text | Google Scholar

Kay, F. U., Canan, A., and Abbara, S. (2020). Future Directions in Coronary CT Angiography: CT-Fractional Flow Reserve, Plaque Vulnerability, and Quantitative Plaque Assessment. Korean Circ. J. 50, 185–202. doi:10.4070/kcj.2019.0315

PubMed Abstract | CrossRef Full Text | Google Scholar

Kini, A. S., Vengrenyuk, Y., Yoshimura, T., Matsumura, M., Pena, J., Baber, U., et al. (2017). Fibrous Cap Thickness by Optical Coherence Tomography In Vivo. J. Am. Coll. Cardiol. 69, 644–657. doi:10.1016/j.jacc.2016.10.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Kolodgie, F. D., Burke, A. P., Farb, A., Gold, H. K., Yuan, J., Narula, J., et al. (2001). The Thin-Cap Fibroatheroma: a Type of Vulnerable Plaque: the Major Precursor Lesion to Acute Coronary Syndromes. Curr. Opin. Cardiol. 16, 285–292. doi:10.1097/00001573-200109000-00006

PubMed Abstract | CrossRef Full Text | Google Scholar

Kume, T., Akasaka, T., Kawamoto, T., Okura, H., Watanabe, N., Toyota, E., et al. (2006). Measurement of the Thickness of the Fibrous Cap by Optical Coherence Tomography. Am. Heart J. 152, e1–755. doi:10.1016/j.ahj.2006.06.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Lv, R., Maehara, A., Matsumura, M., Wang, L., Wang, Q., Zhang, C., et al. (2020). Using Optical Coherence Tomography and Intravascular Ultrasound Imaging to Quantify Coronary Plaque Cap Thickness and Vulnerability: a Pilot Study. Biomed. Eng. Online 19, 90. doi:10.1186/s12938-020-00832-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Mintz, G., Nissen, S. E., Anderson, W. D., Bailey, S. R., Erbel, R., Fitzgerald, P. J., et al. (2001). American College of Cardiology Clinical Expert Consensus Document on Standards for Acquisition, Measurement and Reporting of Intravascular Ultrasound Studies (IVUS). A Report of the American College of Cardiology Task Force on Clinical Expert Consensus Documents Developed in Collaboration with the European Society of Cardiology Endorsed by the Society of Cardiac Angiography and Interventions. Eur. J. Echocardio. 2, 299–313. doi:10.1053/euje.2001.0133

PubMed Abstract | CrossRef Full Text | Google Scholar

Ono, M., Kawashima, H., Hara, H., Gao, C., Wang, R., Kogame, N., et al. (2020). Advances in IVUS/OCT and Future Clinical Perspective of Novel Hybrid Catheter System in Coronary Imaging. Front. Cardiovasc. Med. 7, 119. doi:10.3389/fcvm.2020.00119

PubMed Abstract | CrossRef Full Text | Google Scholar

Porcu, M., Anzidei, M., Suri, J. S., Wasserman, B. A., A Wasserman, B., Lucatelli, P., et al. (2020). Carotid Artery Imaging: The Study of Intra-plaque Vascularization and Hemorrhage in the Era of the "vulnerable" Plaque. J. Neuroradiol. 47, 464–472. doi:10.1016/j.neurad.2019.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Richardson, P., Davies, M. J., and Born, G. V. (1989). Influence of Plaque Configuration and Stress Distribution on Fissuring of Coronary Atherosclerotic Plaques. Lancet 334, 941–944. doi:10.1016/s0140-6736(89)90953-7

CrossRef Full Text | Google Scholar

Saba, L., Anzidei, M., Marincola, B. C., Piga, M., Raz, E., Bassareo, P. P., et al. (2014). Imaging of the Carotid Artery Vulnerable Plaque. Cardiovasc. Intervent Radiol. 37, 572–585. doi:10.1007/s00270-013-0711-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Samady, H., Eshtehardi, P., McDaniel, M. C., Suo, J., Dhawan, S. S., Maynard, C., et al. (2011). Coronary Artery wall Shear Stress Is Associated with Progression and Transformation of Atherosclerotic Plaque and Arterial Remodeling in Patients with Coronary Artery Disease. Circulation 124, 779–788. doi:10.1161/CIRCULATIONAHA.111.021824

PubMed Abstract | CrossRef Full Text | Google Scholar

Sawada, T., Shite, J., Garcia-Garcia, H. M., Shinke, T., Watanabe, S., Otake, H., et al. (2008). Feasibility of Combined Use of Intravascular Ultrasound Radiofrequency Data Analysis and Optical Coherence Tomography for Detecting Thin-Cap Fibroatheroma. Eur. Heart J. 29, 1136–1146. doi:10.1093/eurheartj/ehn132

PubMed Abstract | CrossRef Full Text | Google Scholar

Stone, G. W., Maehara, A., Lansky, A. J., de Bruyne, B., Cristea, E., Mintz, G. S., et al. (2011). A Prospective Natural-History Study of Coronary Atherosclerosis. N. Engl. J. Med. 364, 226–235. doi:10.1056/NEJMoa1002358

PubMed Abstract | CrossRef Full Text | Google Scholar

Stone, P. H., Saito, S., Takahashi, S., Makita, Y., Nakamura, S., Kawasaki, T., et al. (2012). Prediction of Progression of Coronary Artery Disease and Clinical Outcomes Using Vascular Profiling of Endothelial Shear Stress and Arterial Plaque Characteristics. Circulation 126, 172–181. doi:10.1161/CIRCULATIONAHA.112.096438

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, D., Kamm, R. D., Yang, C., Zheng, J., Canton, G., Bach, R., et al. (2014). Image-based Modeling for Better Understanding and Assessment of Atherosclerotic Plaque Progression and Vulnerability: Data, Modeling, Validation, Uncertainty and Predictions. J. Biomech. 47, 834–846. doi:10.1016/j.jbiomech.2014.01.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Virmani, R., Burke, A. P., Farb, A., and Kolodgie, F. D. (2006). Pathology of the Vulnerable Plaque. J. Am. Coll. Cardiol. 47, C13–C18. doi:10.1016/j.jacc.2005.10.065

PubMed Abstract | CrossRef Full Text | Google Scholar

Virmani, R., Burke, A. P., Kolodgie, F. D., and Farb, A. (2002). Vulnerable Plaque: the Pathology of Unstable Coronary Lesions. J. Interven Cardiol. 15, 439–446. doi:10.1111/j.1540-8183.2002.tb01087.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Tang, D., Maehara, A., Wu, Z., Yang, C., Muccigrosso, D., et al. (2018). Fluid-structure Interaction Models Based on Patient-specific IVUS at Baseline and Follow-Up for Prediction of Coronary Plaque Progression by Morphological and Biomechanical Factors: A Preliminary Study. J. Biomech. 68, 43–50. doi:10.1016/j.jbiomech.2017.12.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Wu, Z., Yang, C., Zheng, J., Bach, R., Muccigrosso, D., et al. (2015). IVUS-based FSI Models for Human Coronary Plaque Progression Study: Components, Correlation and Predictive Analysis. Ann. Biomed. Eng. 43, 107–121. doi:10.1007/s10439-014-1118-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Q., Tang, D., Wang, L., Canton, G., Wu, Z., Hatsukami, T. S., et al. (2019). Combining Morphological and Biomechanical Factors for Optimal Carotid Plaque Progression Prediction: An MRI-Based Follow-Up Study Using 3D Thin-Layer Models. Int. J. Cardiol. 293, 266–271. doi:10.1016/j.ijcard.2019.07.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, Z., Tian, J., Ma, L., Du, H., Dong, N., Hou, J., et al. (2015). Comparison of Optical Coherence Tomography and Intravascular Ultrasound for Evaluation of Coronary Lipid-Rich Atherosclerotic Plaque Progression and Regression. Eur. Heart J. Cardiovasc. Imaging 16, 1374–1380. doi:10.1093/ehjci/jev104

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, C., Canton, G., Yuan, C., Ferguson, M., Hatsukami, T. S., and Tang, D. (2010). Advanced Human Carotid Plaque Progression Correlates Positively with Flow Shear Stress Using Follow-Up Scan Data: An In Vivo MRI Multi-Patient 3D FSI Study. J. Biomech. 43, 2530–2538. doi:10.1016/j.jbiomech.2010.05.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Yonetsu, T., and Jang, I.-K. (2018). Advances in Intravascular Imaging: New Insights into the Vulnerable Plaque from Imaging Studies. Korean Circ. J. 48, 1–15. doi:10.4070/kcj.2017.0182

PubMed Abstract | CrossRef Full Text | Google Scholar

Zanchin, C., Ueki, Y., Losdat, S., Fahrni, G., Daemen, J., Ondracek, A. S., et al. (2020). In Vivo relationship between Near-Infrared Spectroscopy-Detected Lipid-Rich Plaques and Morphological Plaque Characteristics by Optical Coherence Tomography and Intravascular Ultrasound: a Multimodality Intravascular Imaging Study. Eur. Heart J. Cardiovasc. Imaging 22, 824–834. doi:10.1093/ehjci/jez318

CrossRef Full Text | Google Scholar

Zhao, K. Q., Xie, X., Yin, H. F., Zhao, J. L., Cao, Z. J., Yang, Y., et al. (2018). Clinical Cohort Study of Imaging Evaluation and Postoperative Pathology of Carotid Vulnerable Plaque. Zhonghua Yi Xue Za Zhi 98, 2424–2428. doi:10.3760/cma.j.issn.0376-2491.2018.30.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: coronary plaque, plaque progression prediction, patient-specific coronary plaque models, vulnerable plaque, optical coherence tomography-based coronary models

Citation: Lv R, Maehara A, Matsumura M, Wang L, Zhang C, Huang M, Guo X, Samady H, Giddens DP, Zheng J, Mintz GS and Tang D (2021) Using Optical Coherence Tomography and Intravascular Ultrasound Imaging to Quantify Coronary Plaque Cap Stress/Strain and Progression: A Follow-Up Study Using 3D Thin-Layer Models. Front. Bioeng. Biotechnol. 9:713525. doi: 10.3389/fbioe.2021.713525

Received: 23 May 2021; Accepted: 20 July 2021;
Published: 23 August 2021.

Edited by:

Juan Carlos Del Alamo, University of Washington, United States

Reviewed by:

Hao Gao, University of Glasgow, United Kingdom
Joshua Thomas Morgan, University of California, Riverside, United States
Gador Canton, University of Washington, United States

Copyright © 2021 Lv, Maehara, Matsumura, Wang, Zhang, Huang, Guo, Samady, Giddens, Zheng, Mintz and Tang. 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: Dalin Tang, ZHRhbmdAd3BpLmVkdQ==

Disclaimer: 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.