Study population
A total of 2,186 T2DM patients were enrolled in the study, 496 patients were excluded from the analysis as survival status (n=460) and diabetes duration (n=36) was not available, thus a total of 1,690 T2DM patients were analyzed. Detailed description of the baseline characteristics is displayed in Table 1.
Table 1. Baseline characteristics of the overall study cohort.
Characteristics
|
Overall cohort (n=1690)
|
Demographics
|
|
Age, years (IQR)
|
63 [54-69]
|
Female, n (%)
|
783 (46)
|
Diabetes duration, years (IQR)
|
10 [5-19]
|
Hypertension, n (%)
|
1135 (67)
|
Dyslipidaemia, n (%)
|
1152 (68)
|
Smoking, n (%)
|
339 (20)
|
BMI, kg/m² (IQR)
|
28.7 [25.4-32.7]
|
Cardiovascular disease
|
|
PCI, n (%)
|
66 (4)
|
PAD, n (%)
|
173 (10)
|
CeVD, n (%)
|
99 (6)
|
CABG, n (%)
|
62 (4)
|
Medications
|
|
Statins, n (%)
|
764 (45)
|
Acetylsalicylic acid, n (%)
|
634 (38)
|
Insulin, n (%)
|
888 (53)
|
Oral antidiabetics, n (%)
|
984 (58)
|
Laboratory parameters
|
|
NT-proBNP, pg/ml (IQR)
|
122 [59-266]
|
Albumin/creatinine ratio, mg/mmol (IQR)
|
0.87 [0.35-2.94]
|
eGFR, ml/min (IQR)
|
72.7 [60.3-85.3]
|
LDL cholesterol, mg/dl (IQR)
|
102 [82-123]
|
HbA1c, % (IQR)
|
7.2 [6.5-8.1]
|
BMI, body mass index; PCI, percutaneous coronary intervention; PAD. peripheral artery disease; CeVD, cerebrovascular disease; CABG, coronary artery bypass graft; NT-proBNP, N-terminal pro-B-type natriuretic peptide; eGFR, estimated glomerular filtration rate; LDL, low-density lipoprotein
Median age of the total study population was 63 years (IQR 54-69), 783 (46%) of the patients were female. CVD was present in 311 (18.4%) patients. According to the ESC/EASD risk model criteria, 25 (1.5%) were classified as moderate, 252 (14.9%) as high and 1125 (66.6%) as very high risk. A total of 288 patients (17.0%) were not classifiable based on the stated ESC/EASD criteria, as 280 patients had diabetes duration less than 10 years with 1 or 2 (but < 3) established CV risk factors and 8 patients presented with a diabetes duration longer than 10 years without any risk factors. Detailed characteristics for the ESC/EASD risk strata are presented in Additional file 1: Table S3. In the overall cohort, 654 patients (39%) had a calculated SCORE risk estimate below 5%, 525 patients (31%) between 5 and 10% and 511 patients (30%) above 10%. The calculated SCORE risk estimates for 10-year fatal increased with ESC/EASD risk category (0% [IQR 0-0] vs 6% [IQR 2-10] vs 8% [IQR 4-12], p<0.001 for the moderate, high and very-high risk category).
Distribution of patients in the ESC/EASD and SCORE risk strata as well as proportion of patients with normal (n=871) and elevated (n=819) NT-proBNP levels at a cut-off 125 pg/ml within these groups are illustrated in Figure 1 (p<0.001 for both models).
Association of NT-proBNP, the ESC/EASD risk strata and SCORE with the primary endpoint CVD death and the secondary outcome all-cause death at 10 years
During 10 years of follow-up, 448 (26.5%) patients died, CVD death accounted to 44.9% (n=201) of all deaths. The cumulative incidences for all risk models stratified into three groups with regards to both endpoints are shown in Figure 2. SCORE and NT-proBNP were both significantly associated with the primary outcome CVD death and the secondary outcome all-cause death (p<0.001 for both) while the ESC/EASD risk model was only associated with all-cause death (moderate risk vs. high risk: p=0.046 and vs. very high risk: p=0.031). Table 2 shows the results of the univariate Cox regression analysis.
Table 2. Association of NT-proBNP, the ESC/EASD and SCORE risk model with outcome in unselected patients with T2DM (n=1690).
|
10-y Cardiovascular death
|
10-y All-cause death
|
|
HR [95% CI]
|
P
|
C-index [95% CI]
|
HR [95%CI]
|
P
|
C-index [95% CI]
|
ESC/EASD risk model
|
-
|
-
|
0.53 [0.50 to 0.56]
|
-
|
|
0.52 [0.50 to 0.54]
|
moderate
|
reference
|
-
|
-
|
reference
|
-
|
-
|
high
|
23.68 [0.08 to 7013.72]
|
0.276
|
-
|
7.49 [1.04 to 53.94]
|
0.046
|
-
|
very high
|
20.88 [0.25 to 1747.39]
|
0.178
|
-
|
8.63 [1.21 to 61.42]
|
0.031
|
-
|
SCORE, %*
|
-
|
-
|
0.63 [0.60 to 0.67]
|
-
|
-
|
0.62 [0.60 to 0.64]
|
<5%
|
reference
|
-
|
-
|
reference
|
-
|
-
|
5-10%
|
3.10 (2.08 to 4.63)
|
<0.001
|
-
|
2.90 (2.24 to 3.76)
|
<0.001
|
-
|
<10%
|
3.79 (2.56 to 5.60)
|
<0.001
|
-
|
3.25 (2.52 to 4.21)
|
<0.001
|
-
|
SCORE, %
|
1.05 [1.04 to 1.07]
|
<0.001
|
0.64 [0.60 to 0.67]
|
1.06 [1.04 to 1.07]
|
<0.001
|
0.66 [0.63 to 0.68]
|
NT-proBNP, tertiles
|
-
|
-
|
0.75 [0.73 to 0.78]
|
-
|
-
|
0.70 [0.67 to 0.72]
|
tertile 1
|
reference
|
-
|
-
|
reference
|
-
|
-
|
tertile 2
|
4.24 (2.19 to 8.21)
|
<0.001
|
-
|
1.76 (1.31 to 2.38)
|
<0.001
|
-
|
tertile 3
|
17.21 (9.3 to 31.78)
|
<0.001
|
-
|
4.90 (3.76 to 6.39)
|
<0.001
|
-
|
NT-proBNP, pg/ml§
|
5.72 [4.68 to 7.00]
|
<0.001
|
0.80 [0.77 to 0.83]
|
3.32 [2.90 to 3.79]
|
<0.001
|
0.73 [0.70 to 0.76]
|
NT-proBNP, > 125 pg/ml
|
7.15 [4.85 to 10.53]
|
<0.001
|
0.71 [0.68 to 0.74]
|
3.14 [2.56 to 3.86]
|
<0.001
|
0.66 [0.63 to 0.68]
|
*refers to SCORE treated as categorical variable (cut-off: <5%, 5-10%, >10%)
§refers to ln-transformed NT-proBNP per 1-IQR increase
Patients with NT-proBNP >125pg/ml had a 7.2-fold and 3.1-fold risk to suffer from death from CVD causes or all-causes at 10 years, respectively, compared to individuals with NT-proBNP within the range considered as normal (p<0.001). NT-proBNP remained a strong predictor of risk irrespective of traditional confounders as age, eGFR, sex, hypertension, smoking and albumin/creatinine ratio for both CVD death (ln[NT-proBNP per IQR increase] adjusted HR: 3.92 [2.92-5.27], p<0.001) and all-cause death (ln[NT-proBNP per IQR increase] adjusted HR: 2.29 [1.89-2.78], p<0.001). Patients with NT-proBNP >125pg/ml had a 7.2-fold and 3.1-fold risk to suffer from death from CVD causes or all-causes at 10 years, respectively, compared to individuals with NT-proBNP within the range considered as normal (p<0.001).
Discriminatory performance of NT-proBNP, the ESC/EASD risk strata and SCORE in in unselected T2DM patients
ROC curves for all risk models and the endpoints 10-year CVD and all-cause death are shown in Figure 3. Additional file 1: Table S4 provides information on sensitivity, specifity, negative predictive value and positive predictive value of the respective risk assessments. The best discriminatory power to predict 10-year CVD and all-cause death was observed for NT-proBNP indicated by highest C-indices compared to both the ESC/EASD and SCORE risk model (C-index: CVD death: 0.80 vs 0.53 vs 0.64; all-cause death: 0.73 vs 0.52 vs 0.66, p<0.001 for all comparisons).
When NT-proBNP was entered as a categorical variable based on tertiles (1st tertile: 59 pg/mL [IQR 59 - 59], 2nd tertile: 122 pg/mL [IQR 90 - 156], 3rd tertile: 376 pg/mL [IQR 267 - 648]) the results remained virtually unchanged (C-index: CVD death: 0.75 vs 0.53 vs 0.64, p<0.001 for all comparisons; all-cause death: 0.70 vs 0.52 vs 0.66, p<0.001 for comparison NT-proBNP vs ESC/EASD, p=0.014 for NT-proBNP vs SCORE).
Net individual risk prediction was significantly improved when assessing NT-proBNP with the ESC/EASD model with a continuous overall NRI of 0.6149 (SE 0.591, p<0.001) for all-cause mortality and of 0.7418 (SE 0.0793, p<0.001) for CVD mortality. Similar results were observed for NT-proBNP and the SCORE model with a continuous overall NRI of 0.6235 (SE 0.0551, p<0.001) for all-cause mortality and of 0.7982 (SE 0.0751, p<0.001) for CVD mortality. Furthermore, a detailed reclassification table comparing NT-proBNP tertiles with the ESC/EASD model as well as the SCORE model for all presented outcomes is given in Additional file 1: Table S5.
Predictive performance of NT-proBNP, the ESC/EASD risk strata and SCORE in distinct subgroups
The following subgroups were investigated: i) patients without CVD and ii) patients without CVD and aged 40-64 years according to SCORE derivation cohort.
Cox regression analysis for the subgroups regarding 10-year fatal outcome are shown in Additional file 1: Table S6. NT-proBNP and SCORE were equally associated with 10-year CVD and all-cause death in subgroups of patients without established CVD and in patients without CVD and aged between 40 and 64 years (p<0.001 for all). The ESC/EASD risk model was only associated with 10-year all-cause death in patients without CVD when comparing the moderate risk with the high risk category (p=0.046).
Figure 4 displays C-statistics of the risk assessments according to the T2DM population studied, the respective ROC graphics are shown in Additional file 1: Figure S1. NT-proBNP was characterized by robustly highest C-indices across both subgroups comparable to the unselected cohort for both endpoints. The ESC/EASD model was characterized by poor C-indices in all groups. The performance of the SCORE risk prediction model improved with progressing exclusivity of patient criteria, performing best in the cohort closest to its derivation population, i.e. in patients without CVD and aged 40-64 years. NT-proBNP outperformed SCORE with regards to CVD death and all-cause death in the subgroup without CVD and was statistically comparable, but showing numerically higher C-indices, in the subgroup without CVD and aged 40-64 years.
Observed vs predicted risk estimate for 10-year fatal CVD by SCORE
The SCORE risk algorithm underestimated the actual risk of CVD death in unselected T2DM patients (observed vs predicted CVD fatal risk: 13% vs 8% [IQR 4-12]; X2 =70.1, p <0.001) and in T2DM patients without CVD (observed vs predicted CVD fatal risk: 10% vs 8% [IQR 4-12]; X2=19.8, p <0.001). When investigating only individuals without CVD and aged 40-64 years, the observed risk was 6% compared with a median predicted risk of 4% [IQR 2-8]. Here, goodness-of-fit for SCORE risk estimate was good with a X2 of 4.3 (p = 0.234).
Predictive performance of NT-proBNP, the ESC/EASD risk strata and SCORE for the secondary outcomes 5-year CVD and all-cause hospitalization
Over a follow-up of 5 years, 1053 (62.3%) patients were hospitalized due to any causes and 367 (21.7%) patients due to unplanned CVD events. Risk for all-cause hospitalization increased by 7% and for CVD hospitalization by 12% per 100pg/ml increase in NT-proBNP (p<0.001 for both). Similarly, SCORE was associated with increased risk for CVD and all-cause hospitalizations (p<0.001 for both) while the ESC/EASD risk stratification model was only associated with all-cause hospitalizations (moderate risk vs. high risk: p=0.033 and vs. very high risk: p=0.002). Cox regression analysis is presented in Additional file 1: Table S7. Figure 5 presents cumulative incidence of 5-year hospitalizations for NT-proBNP. In terms of discriminatory accuracy NT-proBNP was superior to the ESC/EASD risk model for both outcomes (C-index: CVD hospitalization: 0.74 vs 0.54; all-cause hospitalization: 0.62 vs 0.55; p<0.001 for all comparisons). When comparing NT-proBNP with SCORE, NT-proBNP showed superior discrimination for CVD hospitalization (C-index: 0.74 vs 0.62, p<0.001) but not for all-cause hospitalization (C-index: 0.62 vs 0.59, p=0.09). The NRI significantly improved when assessing NT-proBNP with the ESC/EASD model with a continuous overall NRI of 0.3223 (SE 0.0557, p<0.001) for all-cause hospitalizations and of 0.6961 (0.0631, p<0.001) for CVD hospitalizations. Similar results were obtained for NT-proBNP and the SCORE model with a continuous overall NRI of 0.3255 (SE 0.0502, p<0.001) for all-cause hospitalizations and of 0.694 (0.0590, p<0.001) for CVD hospitalizations.
Outcome specifity of NT-proBNP, the ESC/EASD risk strata and SCORE
As indicated by C-statistics, NT-proBNP yielded better discrimination for CVD than for all-cause death (0.80 vs 0.73, p≤0.05) and for CVD than for all-cause hospitalization (0.74 vs 0.62, p≤0.05). No difference for the ESC/EASD risk strata (CVD vs all-cause death: 0.53 vs 0.52, p>0.05; CVD vs all-cause hospitalization: 0.54 vs 0.55, p>0.05) or the SCORE risk estimation (CVD vs all-cause death: 0.64 vs 0.66, p>0.05; CVD vs all-cause hospitalization: 0.62 vs 0.59, p>0.05) could be observed.