A total of 413 patients with AIS confirmed by neuroimaging were finally included for the analysis. The mean age of our population was 75.6 years (SD: 0.6 years), 71.4% were hypertensive and 27.1% diabetic. Among this elderly population, around 11.9% of patients had baseline mRS score ≥3. Stroke severity was mild with median NIHSS score of 6 points (IQR: 3-15). ICH incidence was 12.4% and mortality rate was 17.1%. Baseline characteristics of the study patients are presented in Table 1.
Table 1. Baseline characteristics of the study population
Variable
|
n=413
|
Age; years *
|
75.6 (0.6)
|
Female ‡
|
176 (42.6)
|
Hypertension ‡
|
295 (71.4)
|
Type-2 Diabetes ‡
|
112 (27.1)
|
Dyslipidemia ‡
|
201 (48.7)
|
Antiplatelet therapy ‡
|
129 (31.2)
|
Anticoagulant therapy ‡
|
86 (20.8)
|
SBP at admission, mmHg ƚ
|
152 (135-171)
|
DBP at admission, mmHg ƚ
|
82 (73-93)
|
Serum glucose at admission, mg/dl ƚ
|
113 (99.5-142)
|
Neutrophil count at admission, x109/L ƚ
|
5.2 (4.1-6.7)
|
Lymphocyte count at admission, x109/L ƚ
|
1.8 (1.3-2.4)
|
Baseline NIHSS score ƚ
|
6 (3-15)
|
ASPECTS <7 ‡
|
16 (4.01)
|
ASPECTS, points *
|
9.45 (1.24)
|
Baseline mRS score ‡
|
|
mRS 0
|
164 (39.8)
|
mRS 1
|
125 (30.3)
|
mRS 2
|
74 (18.0)
|
mRS 3
|
32 (7.8)
|
mRS 4
|
17 (4.1)
|
Intravenous thrombolysis ‡
|
151 (36.6)
|
Endovascular treatment ‡
|
70 (17.0)
|
Intracranial hemorrhage ‡
|
51 (12.4)
|
3-month mortality ‡
|
67 (17.1)
|
3-month functional independence ‡
|
227 (57.9)
|
SD: standard deviation; IQR: interquartile range (denoted by 25th-75th percentile); SBP: Systolic Blood Pressure; DBP: Diastolic Blood Pressure; NIHSS: National Institute of Health Stroke Scale; ASPECTS: Alberta Stroke Program Early CT Score; mRS: modified Rankin Scale.
* Continuous variables with normal distributions are presented as means (SD).
ƚ Continuous non-normally distributed variables are presented as medians (IQR).
‡ Categorical variables are presented as n (%).
The most common etiology by TOAST was cardioembolic (n=178, 43%), followed by small-vessel occlusion (n=62, 15%), large artery atherosclerosis (n=43, 10.4%) and others causes (n=7, 1.7%) with 30% of strokes with undetermined etiology (n=123).
Plasma calprotectin levels by outcomes
Median plasma calprotectin level was 1.76 µg/mL (IQR, 1.14-2.66). As shown in Figure 2A, patients who died within 3 months showed higher calprotectin levels [median, (IQR) µg/mL; 2.81, (2.06-4.26) vs 1.58, (1.00-2.29), p<0.001]; as did patients who developed ICH after AIS [1.99, (1.36-3.32) vs 1.72, (1.10-2.56), p=0.030, Figure 2B]. Furthermore, lower plasma calprotectin levels after AIS were observed in patients with 3-month FI [1.49, (0.97-2.23) vs 2.2, (1.41-3.38), p<0.001, Figure 2C].
[Insert Figure 2]
Relationship between calprotectin and clinical and analytical variables
Plasma calprotectin levels correlated with neutrophil (Pearson r, 0.28, p<0.001) and lymphocyte count (-0.17, p=0.007); NLR (0.29, p<0.001); CRP (0.50, p<0.001) and fibrinogen (0.27, p=0.001).
Calprotectin levels were significantly associated with clinical stroke severity by NIHSS score in the ANOVA test (p <0.001) and a significant linear trend was observed (p <0.001) (Additional file 1: Supplementary Figure 1A).
When analysing clinical factors, a positive correlation with age was registered (Pearson r, 0.18, p<0.001) with a significant linear trend (p=0.012) in the ANOVA test. In addition, mRS score at baseline was also associated with plasma calprotectin levels (p<0.02,) with a significant positive linear trend (p<0.001).
Patients with confirmed cardioembolic stroke -according to TOAST criteria- showed higher calprotectin levels than patients with non-cardioembolic stroke [1.97, (1.23-2.97) vs 1.66, (0.99-2.39); p<0.004] (Additional file 1: Supplementary Figure 1B).
Predictors of mortality at 90 days
Patients who died within 90 days were older (median age, 86.6 vs 75.7 years ; p<0.001), more frequently females (58.2 vs 41.8% ; p=0.05), had higher DBP (median DBP, 87.5 vs 82.9 mmHg; p=0.04), had worse NIHSS score (median NIHSS score, 17 vs 4; p<0.001),mRS at admission (median mRS, 1.94 vs 0.91; p<0.001) , and a lower ASPECTS (median ASPECT score, 8.5 vs 9.6; p<0.001) than patients who survived at 90 days.
Univariate analysis confirms sex, age, DBP, NIHSS score, baseline mRS and ASPECT score association with mortality (Table 2). In addition, intravenous treatment with tPA was associated with a higher 3-month mortality and a higher incidence of ICH after stroke. Elevated inflammatory markers (neutrophil count, NLR, CRP and calprotectin) were associated with a higher risk of mortality. After adjusting for potential confounders, plasma calprotectin levels remained associated with 3-month mortality [per a log+1 increase OR, (95%CI); 3.06, (1.67-5.61)], in a model adjusted by age, glucose at admission, stroke severity by NIHSS score, ASPECTS, baseline mRS score, and intravenous thrombolysis. Similarly, NLR and CRP also remained associated with 3-month mortality with the same adjustment [per a log+1 increase: OR 1.98, (1.17-3.35); and 1.39, (1.02-1.89) respectively].
As shown in Table 2, after forward stepwise selection of the strongest predictors, only age [per 10-year increase: OR, (95%CI); 1.73, (1.18-2.53)], calprotectin [per a log+1 increase: OR 3.00 (1.71-5.26)], NIHSS score >14 [OR 7.04, (3.08-16.10)], ASPECTS >7 [OR 0.19 (0.05-0.77)] and baseline mRS score [OR 1.57 (1.14-2.15)] were independent predictors of 3-month mortality, whereas NLR and CRP were excluded from the analysis (p>0.05). This predictive model had a good calibration as indicated by the Hosmer-Lemeshow goodness-of-fit test (p value= 0.51). The averaged measure of fitness in k-fold cross-validation with k=5 was 0.30 for this model (range 0.23-0.37).
Table 2. Univariate and multivariate logistic regression models of 3-month mortality with calprotectin and other baseline characteristics.
|
Univariate Logistic Regression
|
Multivariate Logistic Regression (n=393)
|
Variable
|
OR (95%CI)
|
p value
|
OR (95%CI)
|
p value
|
Age (per 10-year increase)
|
2.21(1.65-2.97)
|
<0.001
|
1.73 (1.18-2.53)
|
0.005
|
Sex (female)
|
2.09 (1.23-3.56)
|
0.007
|
ƚ
|
|
Hypertension
|
2.95 (1.41-6.18)
|
0.004
|
ƚ
|
|
Type-2 Diabetes
|
0.64 (0.34-1.21)
|
0.168
|
ƚ
|
|
Dyslipidaemia
|
0.66 (0.39-1.13)
|
0.131
|
ƚ
|
|
SBP at admission *
(per 10mmHg increase)
|
2.35 (0.48-11.48)
|
0.290
|
|
|
DBP at admission (per 10mmHg increase)
|
1.17 (1.00-1.36)
|
0.048
|
ƚ
|
|
Serum Glucose (mg/dL) *
|
1.97 (0.88-4.44)
|
0.101
|
ƚ
|
|
Neutrophil count (xmm3) *
|
3.64 (1.83-7.24)
|
<0.001
|
‡
|
|
Lymphocyte count (xmm3) *
|
0.42 (0.24-0.72)
|
0.002
|
‡
|
|
NLR *
|
2.41 (1.63-3.58)
|
<0.001
|
ƚ
|
|
C-Reactive Protein *
(mg/L)
|
1.02 (1.01-1.03)
|
<0.001
|
ƚ
|
|
Calprotectin µg/mL *
|
4.60 (2.86-7.40)
|
<0.001
|
3.00 (1.71-5.26)
|
<0.001
|
NIHSS score
7-14 (vs 0-7)
>14 (vs 0-7)
|
3.70 (1.55-8.86)
|
0.003
|
2.23 (0.83-6.02)
|
0.112
|
11.86 (5.85-24.05)
|
<0.001
|
7.04 (3.08-16.10)
|
<0.001
|
ASPECTS ≥7
|
0.08 (0.03-0.24)
|
<0.001
|
0.19 (0.05-0.77)
|
0.020
|
Intravenous thrombolysis
|
2.33 (1.37-3.98)
|
0.002
|
ƚ
|
|
Endovascular treatment
|
1.49 (0.77-2.89)
|
0.240
|
|
|
Intracranial Hemorrhage
|
4.04 (2.10-7.77)
|
<0.001
|
ƚ
|
|
baseline mRS score
|
2.10 (1.66-2.65)
|
<0.001
|
1.57 (1.14-2.15)
|
0.005
|
SBP: Systolic Blood Pressure; DBP: Diastolic Blood Pressure; NLR: Neutrophil-to-lymphocyte Ratio; NIHSS: National Institute of Health Stroke Scale; ASPECTS: Alberta Stroke Program Early CT Score; mRS: modified Rankin Scale.
* log transformed.
ƚ Excluded by forward stepwise selection (inclusion p<0.05 and exclusion p>0.10)
‡ Excluded to avoid collinearity between Neutrophils, Lymphocytes and NLR.
ROC curves analysis for significant univariate models showed that the highest AUC was for NIHSS score [area under the curve (AUC), (95%CI): 0.83 (0.77-0.88), Figure 3A] followed by calprotectin as the strongest predictor of mortality across the inflammatory markers tested [AUC (95%CI): calprotectin 0.77 (0.71-0.83)]. The inclusion in the ROC analysis of age, serum glucose, NIHSS score, baseline mRS, intravenous thrombolysis and ASPECTS as covariates, rendered an AUC for 3-month mortality of 0.88 (95%CI: 0.84-0.93; Figure 3B), that was significantly improved when calprotectin was included in the model [AUC (95%CI): 0.90 (0.86-0.94); p=0.04].
ROC analysis for calprotectin rendered a cut-off value of 2.26µg/mL (sensitivity: 72.7%; specificity: 74.3%), which showed a negative predictive value of 92.9% in our population. When stratifying by this optimal cut-off and controlling for age, serum glucose, NIHSS score, baseline mRs, intravenous tPA administration and ASPECTS; patients with calprotectin ≥ 2.26µg/mL were 4 times more likely to die within 90 days than AIS patients with a calprotectin < 2.26µg/mL [OR, (95%CI); 3.98, (1.88-8.41)].
Even though calprotectin levels seemed to better estimate mortality at 90 days, we assessed whether the combination with CRP and NLR might improve mortality prediction. Univariate ROC curves for CRP and NLR rendered cut-off values of 7.1 mg/L and 3.32, respectively. When considering the combination of these three parameters, patients with calprotectin ≥2.26µg/mL, CRP ≥7.1 mg/L and NLR levels ≥3.32 had an estimated probability of dying in the next 90 days of 49%, whereas calculated mortality risk in patients with negative readings for all 3 markers was 3.7%. Further, we stratified our study population according to the presence of none, 1, 2 or 3 of these inflammatory markers above the cut-off values (Figure 3C). At 90 days, 1.6% of patients with negative readings had died. Also, 12.5% of patients with 1 positive marker, 30.7% of patients with 2 positive markers, and 42.3% of patients with 3 positive markers died.
[Insert Figure 3]
Predictors of 3-month functional independence
Variables associated with functional outcome in the univariate analysis were age, admission NIHSS score, intravenous tPA treatment, endovascular treatment, ASPECTS and ICH, as well as serum glucose, NLR, CRP and calprotectin values (Table 3).
In the multivariate analysis, calprotectin and CRP were not associated with 3-month FI, whereas NLR remained associated after adjusting for potential confounders [per a log+1 increase: OR, (95%CI); 0.53, (0.36-0.77); p=0.001]. Similar results were obtained, when the analysis was repeated evaluating the strongest predictors by stepwise selection (Table 3). Likewise, NLR was the only inflammatory marker which remained significantly associated with 3-month FI [0.52, (0.34-0.80); p=0.003]. Univariate ROC analysis for 3-month FI showed a NLR cut-off point for FI of 2.71 [AUC, (95%CI): 0.65, (0.58-0.71)]. When stratifying for this optimal cut-off and controlling for age, serum glucose, NIHSS score, and ICH, patients with NLR < 2.71 were 2.12 times more likely to be functional dependent at 90 days than AIS patients with an NLR ≥ 2.71 [OR, (95%CI); 2.12 (1.19-3.79)].
Table 3. Univariate and multivariate logistic regression models of 3-month functional independence with calprotectin and other baseline characteristics.
|
Univariate Logistic Regression
|
Multivariate Logistic Regression (n=324)
|
Variable
|
OR (95%CI)
|
p value
|
OR (95%CI)
|
p value
|
Age (per 10-year increase)
|
0.57 (0.46-0.71)
|
<0.001
|
0.58 (0.44-0.76)
|
<0.001
|
Sex (female)
|
0.67 (0.43-1.06)
|
0.089
|
|
|
Hypertension
|
0.74 (0.45-1.21)
|
0.232
|
|
|
Type-2 Diabetes
|
0.87 (0.53-1.42)
|
0.572
|
|
|
Dyslipidemia
|
0.96 (0.61-1.50)
|
0.852
|
|
|
SBP at admission *
(per 10mmHg increase)
|
1.40 (0.36-5.36)
|
0.628
|
|
|
DBP at admission (per 10mmHg increase)
|
1.01 (0.89-1.16)
|
0.830
|
|
|
Serum Glucose at admission (mg/dL) *
|
0.34 (0.16-0.71)
|
0.004
|
0.25 (0.10-0.64)
|
0.004
|
Neutrophil count (xmm3) *
|
0.37 (0.20-0.66)
|
0.001
|
‡
|
|
Lymphocyte count (xmm3) *
|
2.40 (1.50-3.86)
|
<0.001
|
‡
|
|
NLR *
|
0.44 (0.31-0.62)
|
<0.001
|
0.52 (0.34-0.80)
|
0.003
|
C-Reactive Protein (mg/L) *
|
0.64 (0.53-0.78)
|
<0.001
|
ƚ
|
|
Calprotectin µg/mL *
|
0.49 (0.34-0.70)
|
<0.001
|
ƚ
|
|
NIHSS score at admission
7-14 (vs 0-7)
>14 (vs 0-7)
|
0.16 (0.08-0.31)
|
<0.001
|
0.14 (0.07-0.31)
|
<0.001
|
0.12 (0.07-0.21)
|
<0.001
|
0.16 (0.08-0.32)
|
<0.001
|
ASPECTS at baseline ≥7
|
18.95 (2.37-151.57)
|
0.006
|
ƚ
|
|
Intravenous thrombolysis
|
0.47 (0.30-0.74)
|
0.001
|
ƚ
|
|
Endovascular treatment
|
0.51 (0.29-0.89)
|
0.018
|
ƚ
|
|
Intracranial Hemorrhage
|
0.13 (0.06-0.28)
|
<0.001
|
0.32 (0.13-0.83)
|
0.019
|
mRS score at baseline
|
0.66 (0.51-0.86)
|
0.002
|
|
|
SBP: Systolic Blood Pressure; DBP: Diastolic Blood Pressure; NLR: Neutrophil-to-lymphocyte Ratio; NIHSS: National Institute of Health Stroke Scale; ASPECTS: Alberta Stroke Program Early CT Score; mRS: modified Rankin Scale.
* log transformed.
ƚ Excluded by forward stepwise selection (inclusion p<0.05 and exclusion p>0.10)
‡ Excluded by authors to avoid collinearity between Neutrophils, Lymphocytes and NLR.
Predictors of intracranial hemorrhage after ischemic stroke
Development of an ICH after AIS was associated with admission NIHSS score, intravenous tPA treatment and endovascular treatment, along with calprotectin and CRP in univariate analyses. Detailed results are presented in Additional file 2: Supplementary Table 1. Furthermore, ICH was more frequent in females. However, in multivariate analyses only sex, NIHSS score ≥14 and intravenous tPA treatment remains statistically significant, and none inflammatory markers were independently associated with ICH after AIS.
Presence of S100A9 in thrombi retrieved from stroke patients
S100A9 protein was present in all thrombi analysed and quantification of thrombus constituents revealed that stroke thrombi contained on average 12.8% (IQR:1.34-28.89) red blood cells, 15.05% (4.84-26.65) platelets, 0.62% (0.35-1.17) leukocytes, and 3.52% (1.17-6.68) S100A9. Distribution pattern of S100A9 through the thrombus seemed to be related with leukocyte distribution and was primarily found at the interface between red blood cell-rich and platelet-rich areas. Besides the described distribution, S100A9 was also present within platelets islets (Figure 4A-B).
[Insert Figure 4]
Interestingly, a positive correlation was observed between S100A9 and platelets (Pearson r 0.46, p<0.002); leukocytes (0.45, p<0.01) and neutrophil elastase (0.70, p<0.001) thrombi content. When evaluating the amount of S100A9 in thrombi and its association with calprotectin circulating levels, we did not see a correlation. Furthermore, no correlation between thrombus S100A9 content and age, sex, functional outcome and stroke severity was observed (not shown).
Finally, we observed a tendency to higher thrombi S100A9 amount in cardioembolic thrombi and in those who had died (Figure 4C-D). The small sample size of atherothrombotic thrombi and the low number of deaths did not allow us to reach statistical significance (Additional file 2: Supplementary Table 2).