3.1 Clinical characteristics
Of the 909 patients enrolled, 559 were Gensini ‘progressors’ (61.50%) and 283 were Crouse ‘progressors’ (31.13%). The age of the patients ranged from 26 to 90 years, with a mean age of 64.51 years. We divided the cohorts into four groups based on the Crouse and Gensini scores (Figure 1). 240 patients (26.40%) with no progression in both sites were defined as the C-G- group, 386 patients (42.46%) undergoing coronary plaque progression but no carotid plaque progression were assigned to the C-G+ group, 110 subjects (12.10%) with carotid plaque progression alone belonged to the C+G- group, and the remaining 173 patients (19.03%) with progression of both sites were classified as the C+G+ group. 909 patients were randomly assigned to the training (n=636) and validation cohort (n=273) in a 7:3 ratio. Mean baseline levels of low-density lipoprotein cholesterol (LDL-C) were 2.45±0.93 and 2.39±0.85 mmol/L for the training and validation cohort, respectively, and LDL-C levels at 1-year follow-up were 1.87±0.67 and 1.86±0.6 mmol/L, respectively. Optimal LDL-C levels were achieved in 122 (16.58%) and 30 (17.34%) patients in the training and validation groups, as recommended by the 2019 ACC/AHA guidelines,which state that very high-risk patients have a ≥50% reduction in LDL-C from baseline and an LDL-C target of <1.4 mmol/L (<55mg/dL) [17]. The baseline characteristics of patients in the training and validation groups are presented in Table 1. There were no significant differences between the two datasets for each parameter.
3.2 Relationship between carotid and coronary atherosclerosis progression
Among the 909 participants enrolled, 173 patients had progression at both sites. The McNemar’s test was conducted to determine whether carotid plaque progression was associated with coronary artery lesions. Chi-squared statistic was 152.47, P-value was 5.00×10−35 and the Kappa value was equal to -0.004 (<0.75). Additionally, based on the observed probability differences, we estimated that a sample size of approximately 36.4 would be required to achieve a target power of 0.80. The significantly smaller estimated sample size (36.4) compared to the actual sample size (909) suggests that plaque progression was not consistent at both sites.
Table1 Characteristics of training and validation dataset
Characteristic
|
Training(n=636)
|
Validation(n=273)
|
P value
|
Male (n (%))
|
506 (79.56)
|
203 (74.36)
|
0.0827
|
Age(mean±SD)
|
63.64 (8.95)
|
64.51 (9.19)
|
0.1817
|
BMI (mean±SD)
|
25.06 (3.16)
|
24.8 (3.37)
|
0.2674
|
Smoking (n (%))
|
177 (27.83)
|
73 (26.74)
|
0.7358
|
Physical patameters
|
|
|
|
SBP (mean±SD)
|
137.39 (19.57)
|
136.68 (19.67)
|
0.6148
|
DBP (mean±SD)
|
76.02 (11.57)
|
75.18 (11.03)
|
0.3084
|
HR (mean±SD)
|
76.98 (11.05)
|
76.38 (11.45)
|
0.458
|
Medical history
|
|
|
|
Hypertension (n (%))
|
443 (69.65)
|
194 (71.06)
|
0.6708
|
Hyperlipidemia (n (%))
|
202(31.76)
|
91(33.33)
|
0.642
|
Diabetes (n (%))
|
222 (34.91)
|
100 (36.63)
|
0.6183
|
Cerebrovascular disease (n (%))
|
77 (12.11)
|
30 (10.99)
|
0.6316
|
Arrhythmia (n (%))
|
65 (10.22)
|
36 (13.19)
|
0.192
|
Coronary heart disease (n (%))
|
47 (7.39)
|
13 (4.76)
|
0.1435
|
Kidney disease (n (%))
|
39 (6.13)
|
19 (6.96)
|
0.6398
|
Carotid ultrasound
|
|
|
|
Plaque thickness of right common carotid artery (mean±SD)
|
1.68 (1.06)
|
1.75 (1.12)
|
0.3576
|
Plaque thickness of right internal carotid (mean±SD)
|
0.52 (0.99)
|
0.55 (1.03)
|
0.6872
|
Plaque thickness of left common carotid artery (mean±SD)
|
1.9 (1.03)
|
1.93 (0.95)
|
0.6586
|
Plaque thickness of left internal carotid (mean±SD)
|
0.53 (1.01)
|
0.55 (0.96)
|
0.7332
|
Crouse (mean±SD)
|
4.63 (2.71)
|
4.79 (2.77)
|
0.4234
|
Delta of Crouse (mean±SD)
|
0.36 (1.85)
|
0.46 (1.81)
|
0.464
|
Carotid plaque progressors(n (%))
|
191 (30.03)
|
92 (33.7)
|
0.2736
|
Coronary angiography
|
|
|
|
Branches of coronary artery disease(mean±SD)
|
2.3 (0.82)
|
2.22 (0.86)
|
0.1834
|
Gensini (mean±SD)
|
7.34 (9.65)
|
6.44 (7.87)
|
0.1756
|
Delta of Gensini (mean±SD)
|
6.69 (13.14)
|
5.84 (12.84)
|
0.3729
|
Coronary plaque progressors (n (%))
|
402 (63.21)
|
157 (57.51)
|
0.1056
|
Laboratory examinations
|
|
|
|
Hb (mean±SD)
|
138.28 (13.57)
|
138.03 (14.8)
|
0.7976
|
RBC (mean±SD)
|
4.42 (0.45)
|
4.54 (2.06)
|
0.1509
|
WBC (mean±SD)
|
8.02 (34.82)
|
6.41 (2.14)
|
0.4445
|
PLT (mean±SD)
|
188.42 (71.37)
|
179.96 (49.59)
|
0.0753
|
Lymphocyte (mean±SD)
|
28.41 (8.44)
|
29.18 (8.7)
|
0.2085
|
Neotrophil (mean±SD)
|
59.95 (9.33)
|
58.73 (9.24)
|
0.0691
|
Monocyte (mean±SD)
|
8.2 (1.91)
|
8.46 (2.08)
|
0.0676
|
Basophilic granulocyte (mean±SD)
|
0.58 (0.32)
|
0.58 (0.31)
|
0.7822
|
Eosinophilic granulocyte (mean±SD)
|
2.83 (2.16)
|
3.05 (2.2)
|
0.1627
|
ALT (mean±SD)
|
27.75 (27.53)
|
25.48 (18.64)
|
0.2124
|
AST (mean±SD)
|
33.14 (58.67)
|
30.48 (56.61)
|
0.526
|
TBIL (mean±SD)
|
13.2 (5.35)
|
13.33 (6.03)
|
0.7507
|
DBIL (mean±SD)
|
2.42 (1.04)
|
2.48 (1.01)
|
0.4555
|
Cr (mean±SD)
|
82.67 (29.89)
|
82.71 (55.36)
|
0.9894
|
UN (mean±SD)
|
5.82 (1.7)
|
5.86 (1.82)
|
0.746
|
eGFR (mean±SD)
|
84.47 (17.98)
|
85.6 (17.69)
|
0.3841
|
UA (mean±SD)
|
329.38 (86.25)
|
328.9 (86.42)
|
0.052
|
TC (mean±SD)
|
4.11 (1.14)
|
4.04 (1.08)
|
0.4125
|
TG (mean±SD)
|
1.66 (1.06)
|
1.58 (0.97)
|
0.3077
|
HDLC (mean±SD)
|
1.09 (0.27)
|
1.10 (0.28)
|
0.0608
|
LDLC -baseline(mean±SD)
|
2.45 (0.93)
|
2.41 (0.9)
|
0.5463
|
LDLC-follow-up
|
1.87(0.67)
|
1.86(0.6)
|
0.9034
|
LDLC-1450(n(%))
|
122(16.58%)
|
30(17.34%)
|
0.8083
|
Non-HDLC (mean±SD)
|
3.03 (1.11)
|
2.93 (1.06)
|
0.2382
|
RC (mean±SD)
|
0.57 (0.43)
|
0.52 (0.35)
|
0.0693
|
ApoAI (mean±SD)
|
1.2 (0.2)
|
1.23 (0.22)
|
0.0517
|
ApoB (mean±SD)
|
0.82 (0.39)
|
0.8 (0.26)
|
0.441
|
Lp(a) (mean±SD)
|
0.25 (0.24)
|
0.24 (0.22)
|
0.7196
|
ApoE (mean±SD)
|
3.63 (1.06)
|
3.61 (0.96)
|
0.8637
|
Fibrinogen (mean±SD)
|
2.9 (0.71)
|
2.89 (0.63)
|
0.8909
|
INR (mean±SD)
|
0.96 (0.11)
|
0.96 (0.08)
|
0.3962
|
DDI (mean±SD)
|
0.42 (0.66)
|
0.36 (0.55)
|
0.2007
|
PT (mean±SD)
|
11.44 (1.88)
|
11.34 (1)
|
0.3845
|
TT (mean±SD)
|
18.29 (2.48)
|
18.16 (1.54)
|
0.426
|
APTT (mean±SD)
|
30.16 (5.83)
|
30.42 (7.86)
|
0.5676
|
CK (mean±SD)
|
225.36 (827.89)
|
163.75 (511.36)
|
0.2547
|
LDH (mean±SD)
|
185.74 (143.19)
|
167 (113.31)
|
0.0552
|
MYO (mean±SD)
|
55.38 (238.22)
|
38.4 (98.5)
|
0.2559
|
TNI (mean±SD)
|
3.21 (17.86)
|
2.47 (17.64)
|
0.5656
|
CKMB (mean±SD)
|
9.6 (42.76)
|
6.78 (31.07)
|
0.3253
|
CRP (mean±SD)
|
5.38 (14.96)
|
4.28 (12.41)
|
0.289
|
INS (mean±SD)
|
13.45 (26.86)
|
13.62 (28.21)
|
0.9345
|
HbAlc (mean±SD)
|
6.38 (1.18)
|
6.46 (1.31)
|
0.3724
|
ProBNP (mean±SD)
|
409.52 (1555.69)
|
326.96 (976.64)
|
0.4177
|
BMI, body mass index; SBP , systolic blood pressure; DBP, diastolic blood pressure; HR, heart rate; Hb, hemoglobin; RBC, red blood cell; WBC, white blood cell; PLT, platelet; ALT, alanine transaminase; AST, aspartic transaminase; TBIL, total bilirubin; DBIL, direct bilirubin; Cr, creatinine; UN, urea nitrogen; eGFR, estimated glomerular filtration rate; UA, uric acid; TC, total cholesterol; TG, triglyceride; HDLC, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol; LDLC-1450, an LDL-C reduction of ≥50% from baseline and LDL-C goal of <1.4 mmol/L(<55mg/dL); Non-HDLC, non-high-density lipoprotein cholesterol; RC, remnant cholesterol; ApoA1 , apolipoprotein A-I; ApoB , apolipoprotein B; Lp(a), lipoprotein a; ApoE , apolipoprotein E; INR, international; DDI, d-dimer ; TT, thrombin time ; APTT, activated partial thromboplastin time; CKMB, creatine kinase isoenzyme; LDH, lactate dehydrogenase; MYO, myoglobin; CRP, C-reactive protein; INS, insulin;HbA1c , hemoglobin A1c; ProBNP , pro–brain natriuretic peptide; ISR, in-stent restenosis. Compared with Training cohort: ∗P < 0.05; ∗∗P < 0.01; ∗∗∗P < 0.001.
3.3 Defining outcomes and identifying covariates for nomograms
Although the progression of carotid and coronary plaques was not consistent during the 1-year follow-up, 173 subjects experienced concurrent plaque progression at both sites. Hence, we wanted to determine the distinguishing characteristics of these 173 co-progressors in C+G+ group. The goal of our nomogram was to predict the probability of carotid and coronary plaque co-progression. Therefore, C+G+ group was set as case (co-progressors) and the other patients as control (non-progressors). LASSO regression, which can adjust variable selection and regularization and avoid overfitting when fitting the generalized linear model, was used to screen for risk factors (Figure 2). The 57 variables including age, gender, body mass index (BMI), smoking, heart rate (HR), hypertension, dyslipidemia, diabetes mellitus, cerebrovascular disease, arrhythmia, kidney disease, lipid profiles, and other laboratory parameters listed in the Table 1 were taken into LASSO regression and the optimal lambda with 1se criteria resulted in five variables containing the number of involved vessels, BMI, dyslipidemia, diabetes mellitus, and uric acid (Figure 3). Multivariable odds ratios indicated that the number of involved vessels, BMI, dyslipidemia, diabetes mellitus, and uric acid were risk factors for co-progression of carotid and coronary plaque in both training and validation cohorts (Table 2).
3.4 Development and explanation of an ASPRF nomogram for predicting the probability of carotid and coronary plaque co-progression
Based on the regression results, a nomogram incorporating 5 significant risk factors was created (Figure 3). The value of each variable corresponds to a score on the point scale axis. The total score was calculated by adding each individual score, and the probability of carotid and coronary plaque co-progression was estimated by projecting the total score to the lower total point scale. For example, the patient was confirmed to have two lesioned vessels (score=10), a BMI of 32.1 (score=67.5), had dyslipidemia (score=22.5) and diabetes mellitus (score=17.5), and had a uric acid level of 520μmol/L (score=45). Summing these individual scores gave a total point of 162.5, implying an 80% probability of co-progression of carotid and coronary plaque.
3.5 Validation of the ASPRF Nomogram
The C-index of the nomogram for prediction of carotid and coronary plaque co-progression were 0.837 (95% CI= 0.779-0.895) and 0.802 (95% CI=0.762-0.842) in the training and validation cohorts, respectively. The AUCs were 0.823 (95% CI= 0.769-0.887) and 0.803 (95% CI= 0.762-0.842) respectively, indicating satisfactory discrimination of the ASPRF nomogram (Figure 4 A, B). The calibration results showed favorable consistencies between predicted and observed plaque co-progression probability. Brier scores (range = 0-1), a measurement of the difference between the algorithm’s predicted probabilities and actual outcomes, with lower Brier scores representing more accurate predictions of actual events, were 0.106 and 0.124 for the training and validation cohorts, respectively (Figure 4 C, D).
3.6 Decision curve analysis (DCA) of the ASPRF Nomogram
DCA is always used to evaluate whether a model is beneficial for clinical decision making, that is, who should receive treatment or therapeutic interventions[18]. The decision curves for the ASPRF nomogram in the training and validation cohorts are presented in Figure 5. We are concerned with the net benefit (NB) in the DCA, namely income minus expenditure. In our study, the “income” represents true positives - cases of carotid and coronary plaque co-progression; the “expenditure” represents false positives - unnecessary invasive operations, such as CAG rechecking. In the DCA plot, the Y-axis represents the NB, and the X-axis is the threshold probability. As shown in Figure 5, the green and blue lines correspond to predictive models in the training and validation cohorts, respectively, exhibiting high NB over a wide range of risk thresholds from 0 to 0.8. It is conducive to understand DCA to know that the benefit is good, and the referred risk threshold varies depending on the specific clinical scenario.
The clinical impact curve visualized the estimated number of people considered high risk and true positives in the range of 0 to 0.6 by using our prediction model. For example, if 1000 participants were screened using a risk threshold of 0.2, 350 patients would be classified as high-risk for co-progression of carotid and coronary plaques, whereas 150 participants would be true positives. Clinical prediction of coronary plaque progression in an additional 650 participants through our nomogram model may prevent them from undergoing unnecessary invasive CAG manipulations.
Table 2 LASSO and Multivariable logistic regression analysis of progression of coronary and Carotid plaque in the training (n=636) and validation(n=273) dataset
|
|
|
Training
|
|
Validation
|
Variable
|
Coefficient
|
OR
|
95% CI
|
p value
|
Coefficient
|
OR
|
95% CI
|
p value
|
Number of vessels disease
|
0.022
|
1.841
|
1.352-2.559
|
0.003**
|
0.021
|
2.065
|
1.558-4.478
|
0.006**
|
Body mass index
|
0.020
|
1.250
|
1.161-1.349
|
<0.001***
|
0.017
|
1.258
|
1.558-4.478
|
<0.001***
|
dyslipidemia
|
0.081
|
2.771
|
1.769-4.355
|
<0.001***
|
0.067
|
2.826
|
1.354-5.911
|
<0.006**
|
Diabetes
|
0.077
|
2.544
|
1.815-4.488
|
<0.001***
|
0.070
|
1.637
|
0.787-3.406
|
0.185
|
Uric acid
|
0.000
|
1.004
|
1.002-1.007
|
0.001**
|
0.000
|
1.004
|
1.002-1.008
|
0.043*
|
Coefficient in LASSO regression and odd ratio (OR) in multivariable logistic regression of the five risk factors were summarized. LASSO, least absolute shrinkage and selection operator; CI, confident interval. ∗P < 0.05; ∗∗P < 0.01; ∗∗∗P < 0.001.