A novel haemocytometric COVID-19 prognostic score developed and validated in an observational multicentre European hospital-based study
Abstract
COVID-19 induces haemocytometric changes. Complete blood count changes, including new cell activation parameters, from 982 confirmed COVID-19 adult patients from 11 European hospitals were retrospectively analysed for distinctive patterns based on age, gender, clinical severity, symptom duration, and hospital days. The observed haemocytometric patterns formed the basis to develop a multi-haemocytometric-parameter prognostic score to predict, during the first three days after presentation, which patients will recover without ventilation or deteriorate within a two-week timeframe, needing intensive care or with fatal outcome. The prognostic score, with ROC curve AUC at baseline of 0.753 (95% CI 0.723–0.781) increasing to 0.875 (95% CI 0.806–0.926) on day 3, was superior to any individual parameter at distinguishing between clinical severity. Findings were confirmed in a validation cohort. Aim is that the score and haemocytometry results are simultaneously provided by analyser software, enabling wide applicability of the score as haemocytometry is commonly requested in COVID-19 patients.
Introduction
COVID-19 spans a wide clinical spectrum from asymptomatic to severe pneumonia with multiple organ failure (Guan et al., 2020), majorly threatening global health, including that of Europe (Sun et al., 2020a). Early identification of critical patients may reduce mortality by timely interventions (Sun et al., 2020a). Many studies have explored the diagnostic or prognostic value of various factors including age, sex, CT scan, biochemical, and haematological parameters (Galloway et al., 2020; Ji et al., 2020; Lu et al., 2020; Shi et al., 2020a; Shi et al., 2020b; Sun et al., 2020b; Wang et al., 2020a; Yan et al., 2020; Caramelo et al., 2020; Yuan et al., 2020; Chen et al., 2020). Most studies were however geographically limited, had high risk for bias, and had no validation cohort (Wynants et al., 2020). C-reactive protein, ferritin, D-dimer, albumin, urea nitrogen, bilirubin, and lactate dehydrogenase (LDH) levels are cited as indirect indicators of the presence and severity of COVID-19 (Henry et al., 2020; Ji et al., 2020; Liu et al., 2020a; Luo et al., 2020; Sun et al., 2020b; Wu et al., 2020; Yang et al., 2020a; Zhou et al., 2020), as are complete blood count (CBC) and differential count (DIFF) changes, specifically lymphopenia, neutrophilia, high neutrophil-to-lymphocyte ratio (NLR) and thrombocytopenia (Henry et al., 2020; Jiang et al., 2020; Kermali et al., 2020; Khartabil et al., 2020; Lippi et al., 2020a; Lippi and Plebani, 2020b). All aforementioned parameters are widely available, but their value is constrained by significant inter-patient variability and limited specificity.
Newer haematology analysers are capable of functional characterisation of blood cells (Hoffmann, 2014; Buttarello, 2016; Chabot-Richards and George, 2015), including measurement of immune cell activation (Kono et al., 2018; Linssen et al., 2007), which has shown promise in screening for infectious diseases (Prodjosoewojo et al., 2019). The aim of this study is to develop and validate a prognostic score using only haemocytometric data for COVID-19 patients presenting at hospitals, to predict within three days of hospital admission, who will deteriorate and require intensive care unit (ICU) transfer within 14 days of admission. Importantly, our intended purpose of this score is to assist with objective risk stratification to support patient management decision making early on, and thus facilitate timely interventions, such as need for ICU or not, before symptoms of severe illness become clinically overt, with the intention to improve patient outcomes, and not to predict mortality. A secondary aim of this study is to document trends of haematology parameters over time, specifically for the newer parameters, as most of the published data in COVID-19 patients focussed on traditional parameters such as lymphocyte counts, platelet counts, and NLR.
Results
Patient characteristics
In total 999 patients were enrolled in the development cohort (Figure 1). Seventeen patients with underlying haematological malignancies or currently undergoing chemotherapy, were subsequently excluded. Nine hundred eighty-two patients with 2587 haematology measurements (day 0–13), were included to analyse and document temporal haemocytometric data trends. Median age was 71 years (range 18-96) and 68% of the patients were male. Patient distribution by sex, clinical severity, and comorbidities is shown in Table 1 with a breakdown by hospital shown as Supplementary file 1. After excluding 59 patients with missing day 0–3 CBC-Diff data, the remaining 923 patients with 1587 measurements were used for the prognostic score development.
Of the 923 patients, 64 (6.9%) were not hospitalised, 687 (74.4%) were admitted to general wards and 172 (18.6%) went directly to ICU, with a further 9 ICU transfers after day 3. The mortality rate for ICU patients (29.3%; 53/181) and general ward patients (27.4%, 186/678) was comparable (Figure 1).
Patients who died or were critically ill, were significantly older than those less severely ill (median age 74 vs 65 years, p<0.001). Although males outnumbered females (631; 292) in patients that had a severe disease progression, mortality rates were independent of sex. Characteristics of the 923 patients are presented in Table 2. Distribution of clinical severity by age is shown in Figure 2.
Haemocytometry data trends over 14 days of hospitalisation in critical illness (CI) and non-critical (NC) patients
The haematological changes derived from the 2587 measurements, taken from 982 patients at various time points (at the discretion of the attending physician) from day of admission (day 0) up to 13 days, are shown in Figure 3 (lymphocyte-related parameters), Figure 4 (neutrophil-related parameters), Figure 5 (monocyte parameters), Figure 6 (red blood cell-related parameters), and Figure 7 (platelet parameters) and further described below with findings grouped along haemopoietic cell lines.
-
Figure 3—source data 1
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig3-data1-v2.xlsx
-
Figure 3—source data 2
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig3-data2-v2.csv
-
Figure 3—source data 3
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig3-data3-v2.csv
-
Figure 3—source data 4
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig3-data4-v2.csv
-
Figure 3—source data 5
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig3-data5-v2.csv
-
Figure 4—source data 1
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig4-data1-v2.xlsx
-
Figure 4—source data 2
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig4-data2-v2.csv
-
Figure 4—source data 3
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig4-data3-v2.csv
-
Figure 4—source data 4
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig4-data4-v2.csv
-
Figure 4—source data 5
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig4-data5-v2.csv
-
Figure 5—source data 1
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig5-data1-v2.xlsx
-
Figure 5—source data 2
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig5-data2-v2.csv
-
Figure 5—source data 3
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig5-data3-v2.csv
-
Figure 6—source data 1
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig6-data1-v2.xlsx
-
Figure 6—source data 2
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig6-data2-v2.csv
-
Figure 6—source data 3
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig6-data3-v2.csv
-
Figure 6—source data 4
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig6-data4-v2.csv
-
Figure 6—source data 5
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig6-data5-v2.csv
-
Figure 7—source data 1
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig7-data1-v2.xlsx
-
Figure 7—source data 2
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig7-data2-v2.csv
-
Figure 7—source data 3
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig7-data3-v2.csv
-
Figure 7—source data 4
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig7-data4-v2.csv
-
Figure 7—source data 5
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig7-data5-v2.csv
WBC findings
In NC and CI groups, lymphopenia is present for 7 and 10 days respectively and normalises thereafter, as lymphocyte numbers tend to increase after 5 days in both groups. (Figure 3A) The NLR increases in the CI group compared to the NC group, and then gradually decreases again. The differences between the groups remain significant over time (p<0.001) (Figure 3B). Neutrophil counts are normal and remain stable in the NC group, whereas values are mildly elevated and continue to rise over time in the CI group (Figure 4A). This increase in neutrophils is accompanied by a mildly elevated NEUT-RI level in the CI group (Figure 4B). Absolute differences in baseline IG, although statistically significant (p<0.001), are small. After day 2, there is a marked rise in IG in the CI but not in the NC group (Figure 4C). The immature granulocyte-to-lymphocyte ratio (IGLR) trend mirrors that of IG (Figure 4D). AS-LYMPH%/L remains abnormally elevated in both groups from baseline throughout the first two weeks of hospitalisation (Figure 3C). The RE-LYMPH, minus AS-LYMPH (a subset of RE-LYMPH, which is absent in healthy individuals) indicates that the reactive lymphocytes observed in both NC and CI groups are largely comprised of AS-LYMPH (Figure 3D), in keeping with a predominant B-cell rather than T-cell response. Absolute monocyte counts are within normal limits for the NC group and remain stable throughout. Likewise, values are largely normal within the CI group, but show an upward trend with mild monocytosis evident from about day 10 onwards (Figure 5A). RE-MONO%/M remains stable over time, and largely within normal limits in the NC group. In contrast, in the CI group monocyte activation increases up to days 3 to 4, with values returning to within the reference range after a week (Figure 5B).
RBC findings
In all patients, there is a gradual drop in HGB, also after adjusting for age and sex (Figure 6A), with differences between the groups becoming increasingly wider from day 5 onwards. After day 7, HGB continues to decline only in the CI group. RET# remain low in both groups despite dropping HGB in the first week but RET# shows a consistent rise thereafter in the CI group, towards the upper limit of the normal reference range (Figure 6B).
The Delta-He is negative and remains relatively stable in the NC group (Figure 6C). In contrast, Delta-He drops progressively in the CI group, reaching its nadir at about day 7, and then rises towards zero, primarily due to an improvement in the RET-He values (data not shown). NRBC# is almost zero in the NC group (within normal range) but rises sharply and progressively at about day 5 in the CI group (Figure 6D).
PLT findings
PLT are largely within the normal range but show a progressive upward trend over time for both groups, with patients in the CI group manifesting with mild thrombocytosis from about day 10 onwards (Figure 7A). The IPF# initially is within the normal reference range for both groups but over time the CI group shows a gradual increase, exceeding the upper limit of the reference range in parallel to PLT (Figure 7C), whereas IPF(%) remains within normal limits throughout for both groups (Figure 7D). The platelet-to-lymphocyte ratio (PLR) is abnormally elevated for both groups throughout, with values slightly higher in the CI group, but only until day 5, after which the NC and CI groups overlap (Figure 7B).
Formulation and performance of the haemocytometric COVID-19 prognostic score in the development cohort
In the analysis of 1587 samples, 923 patients (days 0-3), six variables (parameters or ratios thereof) were identified that fulfilled the prespecified selection criteria, namely sufficient discriminatory power between NC and CI patient groups (p≤0.001) and at least 20% of all CI results outside of the normal range. These were NLR (p<0.0001, 77.1%), IGLR (p<0.0001, 49,4%), RE-MONO%/M (p<0.0001, 37.3%), AS-LYMPH%/L (p<0.0001, 90,4%), Delta-He (p<0.0001, 85,7%), and NRBC (p<0.0001, 46.0%). Other parameters also fulfilled these criteria, but these were not selected as they are interdependent on others already included. These were WBC, NEUT#, LYMPH# (represented by NLR, IGLR, AS-LYMPH), MONO# (represented by RE-MONO), RE-LYMPH (represented by AS-LYMPH), RBC, HGB, and HYPO-He (represented by DELTA-He) as well as NEUT-RI which is commonly elevated in bacterial infection (Prodjosoewojo et al., 2019).
These six variables (NLR, IGLR, RE-MONO%/M, AS-LYMPH%/L, Delta-He, NRBC) were assigned a score from 0 to maximum of 4 points each. Applying this score to the CI and NC patient groups, no overlap (SEM) was noticed in the mean values except on day 13, where the available data points were few (Figure 8A). These combined six variables had a sensitivity of 68% (249/366) in identifying CI patients in the first 3 days. However, clear haematological abnormalities, especially HGB, HYPO-He, PLT, and IPF# were observed in 8 of the undetected CI patients. These parameters were therefore added to the score with a maximum of 1 point each. The sum of 10 variables makes up the final prognostic score value, with a minimum of 0 and a theoretical maximum 28 points. The prognostic score cut-off values and point allocation matrix is shown in Table 3.
-
Figure 8—source data 1
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig8-data1-v2.xlsx
-
Figure 8—source data 2
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig8-data2-v2.csv
-
Figure 8—source data 3
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig8-data3-v2.csv
-
Figure 8—source data 4
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig8-data4-v2.csv
Prognostic score performance in the 923 patients using 3 as the cut-off, correctly identified 70.5% (95% CI 66–75) of patients finally classified as critical illness (CI group) on days 0–3 with sensitivities of 62.3% (95% CI 55–69), 74.5% (95% CI 61–85), 75.6% (95% CI 64–85), and 87.0% (95% CI 77–94) on days 0, 1, 2, and 3 respectively, and 93% (95% CI 81-99) on day 6. Moreover, all patients classified with critical illness, whether they recovered or died, had an increasing score value, compared to those that recovered without ICU intervention (Figure 8B). Specifically, the scores for non-ICU patients who subsequently deteriorated and died (n=186, Figure 1), were already notably higher during the initial phase of hospitalisation with a median value of 3 (95% CI 2-4), 4 (95% CI 1.5-6.7), 5 (95% CI 4.0-6.0), and 5 (95% CI 4.4-7.6) on days 0, 1, 2, and 3 respectively, (Table 4), compared with 1 (95% CI 1.0-2.0) to 2 (95% CI 1.0-3.0) for patients who recovered without ICU intervention.
In patients classified as having ‘critical illness’, irrespective of whether the final outcome was death (with or without ICU intervention) or recovery, the prognostic score values increased progressively over time with no statistically significant difference (p>0.05) in prognostic score values between survivors and non-survivors, confirming that our score predicts a critical clinical course, but not the outcome, that is recovery or death.
For the initial period of hospitalisation (<4 days) and for the 14-day prediction time horizon, the prognostic score was better than NLR at differentiating clinical severity, with a higher AUC at all time points (Table 5). Notably, the cut-off value that determined the best AUC was consistent for the prognostic score (≥3 or >4) over time, whereas it ranged from >4.7 to >11.6 for NLR. AUC values for platelet-to-lymphocyte ratio (PLR), LYMPH, MONO, and PLT were all low (range 0.501-0.647) with highly variable cut-off values at the different time points. The AUC comparisons for the prognostic score and these individual variables is shown in Figure 8C.
In investigating if the score can predict severity independent of the classical risk factors such as age and presence of comorbidities, using a Mann-Whitney test, it was found that the prognostic score was significantly higher in the CI group than the NC group across all age groups and for all age groups segregated by the presence or absence of comorbidities, with the exception of patients 84 years and older with reported comorbidities (Table 6 and Figure 9). The median difference in prognostic score values between the NC and CI groups ranged from 2 to 7 points.
Prognostic score validation
For 202 patients (Figure 10), 65.8% male and median age 70 years (range 22-93), with 217 (165 CI; 52 NC) CBC-DIFF-RET day 0–3 measurements available, the prognostic score gave an AUC of 0.797 (95% CI 0.724-0.829), which is comparable to the prognostic score performance in the development cohort (AUC 0.753, 95% CI 0.723-0.781). Mortality rate was 20.8% (42/202) and outcome was correctly identified in 72% (91/127) of CI patients (development cohort 70.5%). Except for day 1, score values of the NC and CI groups over the 14-day period (718 measurements) did not overlap (SEM) (Figure 11A), with AUC of prognostic score superior to NLR (Figure 11B).
-
Figure 11—source data 1
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig11-data1-v2.xlsx
-
Figure 11—source data 2
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig11-data2-v2.csv
-
Figure 11—source data 3
- https://cdn.elifesciences.org/articles/63195/elife-63195-fig11-data3-v2.csv
Discussion
We showed that SARS-CoV-2 infection is accompanied by haemocytometric changes over time and that distinct haemocytometric parameters, combined in a COVID-19 prognostic score, can be used early on to identify those patients likely to deteriorate thereafter and thus may benefit from ICU admission. Moreover, our data suggest that parameters reflecting the activation or functional status of blood cells are better disease severity indicators than traditional parameters, such as lymphocyte or platelet counts.
In COVID-19, lymphopenia has been assigned a key role based on a higher incidence and greater suppression observed in ICU patients (Terpos et al., 2020). Our data indicate that lymphocyte count had no significant prognostic value (AUC 0.550; 95% CI 0.523-0.576), although relative presence of lymphocytes did, hence the incorporation of NLR, IGLR, and AS-LYMPH%/L (reflecting lymphoplasmacytoid lymphocytes [Linssen et al., 2007]) in the prognostic score. Increased lymphocyte activation has been documented in COVID-19 (Fan et al., 2020; Chong et al., 2020). Our findings concur with a previous study (Yip et al., 2020) that reported AS-LYMPH as a strong predictor of clinical severity in COVID-19 patients. Furthermore, a specific “hourglass” appearance on the WBC scattergram on Sysmex analysers, representative of lymphoplasmacytoid lymphocytes, was reported to have a high positive predictive value to detect COVID-19 (Osman et al., 2020).
In contrast to lymphocytes, neutrophils, including precursors, have a tendency to increase in COVID-19 (Mitra et al., 2020). IGs, representing metamyelocytes, myelocytes,, and promyelocytes (Briggs et al., 2011), were commonly present in our patient population, especially in those more severely ill. The importance of IGs in management of sepsis has been reported (Ayres et al., 2019; Nierhaus et al., 2013). Of note, increases in neutrophil counts and neutrophil activity (NEUT-RI) were dissimilar, unlike observations in bacterial infections (Park et al., 2015). RE-MONO%/M is abnormal in critical illness cases only, in line with a previous report attributing a key role for monocytes and macrophages in severe COVID-19 (Merad and Martin, 2020).
Erythropoietic changes have been reported, mostly showing low HGB levels (Fan et al., 2020; Lippi and Plebani, 2020b; Sun et al., 2020c). We found that whilst HGB levels decrease in COVID-19 patients, the erythropoietic response to anaemia, indicated by RET and reticulocyte production index (data not shown), were mostly normal. Haemoglobinisation of reticulocytes, as indicated by negative DELTA-He levels, is however significantly compromised, specifically in more severe cases, possibly due to ongoing inflammation (Weimann et al., 2016). NRBCs are absent in peripheral blood of healthy adults. Their presence, without reticulocytosis, in severe COVID-19 cases indicates haematopoietic stress, probably due to prolonged hypoxia or inflammation (Danise et al., 2011). Furthermore, NRBCs, were reported as a marker of disease severity in ARDS patients, indicating a higher risk of death (Menk et al., 2018).
Contrary to other studies (Henry et al., 2020; Jiang et al., 2020; Wang et al., 2020b; Yang et al., 2020b), PLT at presentation were similar between CI and NC cases, mostly within normal limits with no sign of increased platelet consumption as IPF also remained normal. PLT, and IPF#, tended to increase with disease severity. Higher PLT have been previously reported in severe COVID-19 (Qu et al., 2020). So thrombocytosis, more than thrombocytopenia may be linked to severe COVID-19 which is in contrast to guidelines to identify severe pneumonia (Metlay et al., 2019).
A recent meta-analysis (Zeng et al., 2020) concluded that severe COVID-19 patients had higher neutrophil counts and NLR, and lower lymphocyte counts than those with non-severe COVID-19, and that these basic parameters might help clinicians to predict the severity and prognosis of COVID-19. Although our findings concur with their observed clinical severity-based WBC differences, the discriminating power, early on during hospitalisation and thus value to determine prognosis, was insufficient. A previous report about the prognostic value of NLR (Liu et al., 2020b) is also not supported by our data. Altogether, our findings indicate that new parameters, reflecting functional status of blood cells, are more frequently outside reference ranges in COVID-19 than classical parameters such as lymphocytes, neutrophils or platelets. However, none of the measured parameters, traditional or novel, alone could discriminate patients based on disease severity. The prognostic score we developed used multiple parameters, representing the three haemopoietic cell lines, with the aim to distinguish CI from NC COVID-19 patients. Different cell lines may not all be equally affected by COVID-19 at the same time. Thus, it is not unexpected that multiple parameters together are better at predicting disease severity, as duration of symptoms at presentation differs widely in individual patients.
The patient group is inherently heterogeneous as health-seeking behaviour may vary widely between patients: some patients may present early whilst others present late at health facilities. Furthermore, recall when symptoms started may be unclear at time of presentation to the hospital. Therefore, a patient on day 3 of hospitalisation could be in an earlier phase of infection than a patient on day 0. In this regard, as the intended use of the prognostic score was to identify at an early stage (once medical care had been sought) who is likely to deteriorate, to keep it as simple and practical as possible, we analysed the data from day 0–3 as a collective dataset for the purpose of identifying which parameters, and at what cut-off values, should be incorporated in the score. Likewise, we did not compare score values from the early phase of hospitalisation (days 0–3) with the later phase (days 4–13) in keeping with the intention to predict outcome early on before signs of critical illness may become clinically overt. We did however plot the data over time to document the patterns as this may be useful in the future to assess response to specific COVID-19 therapies.
The prognostic score correctly identified 70.5% of these patients during first 3 days after hospital admission, with similar performance confirmed in the validation cohort (72.0%). Prognostic score trends over 14 days confirmed a stable clinical course in NC patients and disease progression in CI patients, peaking on day 6 (sensitivity 93%). In the development cohort there was a distinct progressive upward trend of prognostic score values from day 0, peaking on day 6 in the CI group. In the validation cohort, the score also peaked on day 6, but the day 0 score value started relatively high, dropping on day 1, giving the appearance of convergence with the NC group on day 1. The validation cohort was comprised of patients from only two hospitals. Study enrolment at these two hospitals took place at different times of the pandemic, during which time ICU capacity was ramped up significantly. In this regard we speculate that ICU access differed at the two hospitals, hence patients with a similar degree of disease severity may have been admitted to ICU at one hospital and remained on the general ward at the other. Furthermore, all day 0 patients in the CI group were from a single hospital, and numbers were small relative to the development cohort. Our ethics approval did not permit us to review the individual patient clinical records of six patients in the CI group that already had very high (10–16 points) prognostic score values on day 0, four of whom also had no further follow-up samples which would have brought the score value averages on the days following admission, upwards. Any such inter-hospital variability would have been masked in the development cohort, because of the much larger sample size and representation from nine hospitals, compared to the validation cohort where such differences may have been more readily exposed. We therefore speculate that the high average day 0 score value patients in the validation cohort are outliers, and that the day 1 data is more representative of the expected prognostic score values in the validation cohort data set.
As has been widely reported in the literature, our data too shows that males were predominantly affected, and that disease severity was associated with increasing age and presence of comorbidities in general (Table 1, Figure 2). However, not all young patients had a mild course, and not all old patients with comorbidities were critical. Systemic inflammation is an important factor driving disease severity. Our prognostic score, incorporating the activation status of immune cells, may therefore have additional value, especially on an individual patient level, over classical risk factors such as age, gender and comorbidities in discriminating between NC and CI patients, (Table 6, Figure 9). As such, our prognostic score may assist in identifying any patients at risk for severe disease progression, being young or old, male or female, with or without comorbidity and by doing so, support individualised treatment decisions with objective data.
Notably, the mortality rate was relatively high in patients on the general ward, comparable to that of patients admitted to ICU, possibly due to ICU bed shortages or unfamiliarity with COVID-19 at that time. Whatever reason, we assumed that the need for more intensive treatment should have been considered for all patients that died in the general ward. Once a patient is overtly critically ill, clinical judgement will suffice to prioritise intensive care for such a patient. Importantly, our score strives to identify the deranged immune response as a harbinger of criticality (organ failure) before such organ failure is clinically evident. Our data (Figure 8B) shows that prognostic score values of patients with a critical clinical course (CI group) were indistinguishable based on outcome, namely whether they recovered or died, confirming that the clinical applicability of our score is to predict on an individual patient level who is likely to have a critical clinical course, but not to predict mortality.
In our study, the prognostic score was calculated retrospectively using the haemocytometric data from each individual sample measurement exported into Microsoft Excel (Microsoft Corp., USA). For future clinical practice, we envisage that the prognostic score will be automatically calculated by the laboratory information system and that the result is given on request, together with complete blood count analysis. Besides serving as a risk stratification tool for clinical decision making early on, we postulate that the prognostic score, which provides a snapshot in time of the phase of an individual’s immune response, may be promising as an aid to patient selection for future clinical trials exploring therapeutic options for COVID-19.
Most previous risk score development studies were limited by small sample size, single centre, univariate analysis, and lack of a validation cohort. In contrast, the 4C mortality score, developed and validated by the International Severe Acute Respiratory and Emerging Infections Consortium (ISARIC), incorporating patient demographic information (age, sex, comorbidities), clinical observations (respiratory rate, peripheral oxygen saturation, level of consciousness) and blood parameters (C-reactive protein (CRP), urea level), incorporated more than 35,000 and 22,000 patients into the development and validation cohorts respectively, and showed better performance at predicting mortality in COVID-19 patients than all previously published scores (Knight et al., 2020). The ISARIC investigators conclude that their score is easy to use as commonly available at hospital admission. Our score has several advantages: firstly, only objective measurements of a haematology analyser are used whilst the ISARIC 4C mortality score uses clinical observation parameters which may be subject to interpretation; secondly, our score is aimed to be automatically generated whilst 4C mortality score needs to be calculated; thirdly, different laboratory measurements (oxygen saturation, CRP, urea) are required to calculate the 4C mortality score, whilst our score uses the globally most commonly requested laboratory examination for patients attending health facilities; fourthly, haematology analysers are widely available globally more so than CRP measurements; and fifthly, the 4C mortality score is aimed at patients admitted to the hospital, whilst our score was developed for patients presenting at the hospital (of whom some were never admitted).
The strength of our study is the inclusion of a relatively large group of confirmed COVID-19 cases from multiple centres and countries, including a validation cohort, Furthermore, we believe that an advantage of our prognostic score is that all input data required to calculate the score value are generated from a single haematology profile test, which is the most common routinely requested baseline blood test in all patients globally.
A limitation of our study is its retrospective nature, as data retrieved from hospital records were sometimes incomplete. Our study was performed at a time that COVID-19 was a new disease entity that constrained the health care in many of the participating centres which may have affected management decisions and therefore study outcome parameters. Data from out-patient settings including more mild cases, and from nursing homes that usually accommodate high risk patients, are needed. Furthermore, clinical data collection was limited, including comorbidity affecting COVID-19 susceptibility and ICU admission decisions, notably as the demand for ICU beds was greater than the availability at the time of our study. Importantly, the conditions of our fast-tracked ethics clearance to facilitate rapid study initiation did not permit data collection about bacterial superinfections and medication (antibiotics, corticosteroids), whilst these factors may affect outcome and haemocytometric parameters.
Finally, our prognostic score includes Sysmex unique parameters. This is a limitation as the score is not universally applicable to all haematology analysers, although the concept is transferable (see Supplementary file 2 for parameters available on other manufacturer haematology platforms). However, it is the very ability to quantify blood cell activation, a reflection of the general immune response status of an individual, that has rendered our prognostic score (which incorporates cell activation parameters such as reactive monocytes and antibody-synthesising lymphocytes amongst others) better than using only standard parameters, such as neutrophil-to-lymphocyte ratio, which are universally available on all systems, at least in our patient dataset.
Conclusions
Our finding of potential usefulness of extended haemocytometry may be impactful as Sysmex haematology analysers are widely available. Haematology blood profile requests are common, inexpensive, quick, highly standardised, quality-controlled baseline tests. Furthermore, this investigation is requested in febrile patients and those with non-febrile conditions. As the latter patients are at higher risk for serious COVID-19, early recognition is important to provide supportive care.
Materials and methods
Study design, sample size, and participants
Request a detailed protocolWhilst it was not possible to calculate an appropriate sample size due to the rapid escalation of the COVID-19 pandemic and concomitant resource constraints experienced by the study centres during the time of planned data collection, the study team set a minimum target of 500 patients, of which at least 250 were admitted to intensive care, and that the study would remain open for enrolment of patients until 6 April 2020 to increase the patient numbers as much as possible.
In this explorative multicentre study patients were prospectively enrolled into a prognostic score development cohort from 21 February to 6 April 2020, with follow-up to document clinical outcome until 9 June 2020, from seven hospitals in the Netherlands and one each in Italy and Belgium.
Data analysis and prognostic score development was performed retrospectively from 9 to 29 June 2020. Thereafter, eligible patients presenting between 7 April to 15 May 2020 and 3 March and 17 May 2020 respectively at two hospitals in the Netherlands were retrospectively enrolled into a validation cohort.
Inclusion criteria for both cohorts were: primary presentation at participating hospitals; RT-PCR confirmed COVID-19; age ≥18 years;≥1 CBC-DIFF (±reticulocyte measurement (RET)) analysed on a Sysmex XN series haematology analyser (Kobe, Japan) as part of routine care; initial management decision record (self-isolation, general ward admission, ICU admission). Any documented comorbidities, duration of symptoms, subsequent ICU transfer, and outcome (recovered, died, unknown) were recorded.
Patients younger than 18 years were excluded as ethical approval was limited to adult subjects. Patients with an unknown outcome were excluded from analysis. Also, patients with a known underlying condition, such as a haematological malignancy or disorder, or concurrent chemotherapy, were excluded as such conditions are associated with abnormal blood cell counts which may either mask the severity of disease or falsely amplify it.
All CBC-DIFFs (+/- RET) done at presentation and throughout hospitalisation, of eligible patients, were included for retrospective analysis. All patients that met the inclusion criteria were included in the longitudinal haemocytometric data trend analysis (secondary aim) whereas only those patients that had at least one haematology profile measurement available on one or more days, on day 0, 1, 2, or 3 were included in the prognostic score development (primary aim).
Due to the rapidly changing dynamics of the clinical course of COVID-19, and the timing of blood testing being entirely at the discretion of the clinician, if a patient had more than one measurement on a single day all measurements were included in the analysis as our objective is to have a prognostic score that is universally applicable to measurements taken at any time point between days 0-3.
Reference data were obtained from a random sample of 12,782 healthy individuals from a large population-based cohort in the Netherlands (http://www.lifelines.nl).
Haemocytometry
View detailed protocolParameters included in this study are shown in Table 7.
In brief, standard CBC-DIFF parameters plus nucleated red blood cells (NRBC), immature granulocytes (IG), neutrophil reactivity index (NEUT-RI), and neutrophil granularity index (NEUT-GI) were measured. Neutrophil-to-lymphocyte ratio (NLR) and immature granulocyte-to-lymphocyte ratio (IGLR) were calculated.
Where available, we assessed reticulocyte count (RET), reticulocyte haemoglobin content (RET-He, the optical measurement of reticulocyte haemoglobinisation), the difference in reticulocyte and RBC haemoglobin content (Delta-He) (which is also a real time marker for iron bioavailability, with a negative value being an indirect marker of monocyte activation), the percentage of hypochromic (HYPO-He), and microcytic RBC (MICRO-R).
Anonymous XN-analyser raw data files were provided to Sysmex Europe collaborators for analysis using virtual analyser software, specifically to obtain reactive lymphocyte (RE-LYMPH), antibody-synthesising lymphocyte (AS-LYMPH), and reactive monocyte (RE-MONO) counts. If the analyser was equipped with a PLT-F (optical platelet count) measurement channel, and the initial CBC platelet count result triggered a reflex PLT-F measurement, or this was included in the initial request as a default profile, then immature platelet fraction (IPF#, IPF%) values measured in this channel were included in the data analysis. However, as most analysers in this study were not equipped with the PLT-F channel, the IPF# and IPF(%) values were partly derived from the RET channel measurement using virtual analyser software.
Definitions
Disease severity was scored as follows: Mild: no hospitalisation, recovered; Moderate: ≤5 days hospitalisation without ICU/ventilation, recovered; Severe: >5 days hospitalisation without ICU/ventilation, recovered; Critical: ICU/ventilation at any stage of hospitalisation, recovered; Fatal: death.
Outcomes were classified as either recovered, died or unknown. ‘Recovered’ refers to all patients that had shown significant improvement and/or resolution of their clinical symptoms and thus discharged from hospital at the discretion of the attending physician. All patients that were never hospitalised (mild disease classification) were assumed to have recovered.
‘Died’ refers to all patients that died in hospital. Patients who were still hospitalised at the time of study termination or those that had been transferred to another healthcare facility and who’s outcome could not be retrieved, were classified as having an ‘unknown’ outcome. ICU admission, whether direct or a subsequent transfer from the general ward, was not considered an outcome, but rather a criterion to define critical illness.
Two disease severity groups were defined, assuming that all patients that died would have needed ICU admission, and would have been admitted to ICU, had an ICU bed been available: (1) non-critical illness (NC) group comprising patients classified as mild, moderate or severe, and (2) the critical illness (CI) group comprised of critical (ICU survivors) and all fatal outcome patients, irrespective of general ward or ICU admission.
Day 0 refers to the day of first presentation at hospital, and the day of admission for those patients requiring hospitalisation. Day 1 refers to one day after the day of admission, or alternatively, one hospital bed night. Day 2 refers to 2 days, and so on.
Trend analysis of single parameters and development of a prognostic score to predict disease severity and progression in COVID-19 patients
Request a detailed protocolAll available haemocytometric measurements from days 0-13 were included in the individual haematology parameter longitudinal analysis. Haemocytometric data were grouped according to clinical severity, symptom duration and days of hospitalisation, and analysed up to day 14 and compared with healthy controls, to identify specific patterns and trends. Haemoglobin values were adjusted for age and gender.
The prognostic score development process is outlined in Figure 12. The day 0 to 3 data points from the trend analysis were utilised for this purpose. In brief, each parameter was analysed univariately with regard to non-critical or critical illness classification and backward-selection multiple logistic regression analysis was conducted to select parameters for inclusion in the prognostic score, using the following criteria: (1) significant difference on days 0-3 (p≤0.001); (2)>20% abnormal in the CI group; (3) only dominant parameter selected if ≥2 eligible parameters are inter-dependent.
Receiver operating characteristics (ROC) curve analysis of selected parameters was applied for calculation of optimal cut-off values associated with sensitivity and specificity and used as the baseline for points assignment, with further increments for increasing specificity and likelihood ratio in discriminating between NC and CI groups (maximum four points). The score represents the sum of individual parameter points.
Score values were calculated for day 0-3 measurements to determine the cut-off value that best discriminate between the NC and CI groups. To further enhance score sensitivity without compromising specificity, additional parameters were included (maximum one point each) if they were significantly abnormal from day 4 onwards with >95% specificity.
The 14-day prognostic score time horizon was plotted using all available measurements and its performance in predicting disease severity validated in an independent patient cohort.
Only haemocytometry parameters were used as predictors in the prognostic score development. As these data are generated from automated haematology analysers, and do not rely on interpretation, predictors for the prognostic score were automatically blinded.
Also, the disease severity groups were defined in advance of patient enrolment, and patients were classified according to objective data (length of hospital stay, general ward or ICU, recovered and discharged from hospital or died) providing by the enrolling study centre prior to the commencement of data analysis and score development. In this regard, the authors involved in score developed had no influence over assessment of prediction of disease severity of individual patients, and hence assessment of outcome was deemed to have been blinded.
Statistics
MedCalc Statistical Software version 19.2.1 (MedCalc Software Ltd, Ostend, Belgium; https://www.medcalc.org; 2020) was used. Differences between patient groups were assessed using Student`s t-test for normality and if rejected (D'Agostino-Pearson test for normal distribution), the Mann-Whitney U-test for non-normally distributed variables was used. The Hodges-Lehmann median difference (md) was used to describe the differences in values of non-standard distributed parameters. To evaluate predictive values, we calculated the AUC of ROC curves, and 95% confidence intervals for each of the first 14 days of hospitalisation. This was done for both the development and validation cohort. The influence of two categorical variables was visualised using the clustered multiple comparison graph with standard error of mean (SEM). Confoundance of prognostic score, age, gender, and presence of comorbidities was tested by logistic regression. Box and whisker plots were used to visualise comparisons of multiple groups.
Missing data
Request a detailed protocolThose samples that did not have a RET channel measurement were excluded from the trend analysis of all parameters measured or derived from this measurement channel (Table 7 ). Only samples with full profile (CBC-DIFF-RET) measurements were included for determination of cut-off values for prognostic score calculation. Only patients with full profile measurements were included in the validation cohort. No imputation was done for missing data. All patients that were enrolled from those centres that did not provide any information on the presence or absence of comorbidities, were excluded from the statistical analysis of prognostic score performance versus presence of comorbidities (classical risk factors) in predicting disease severity outcome.
Data availability
All data analysed during this study are included in the manuscript and supporting files. Source data files have been provided for figures 2, 3, 4, 5, 6, 7, 8, 9 and 11, and tables 1, 2, 4, 5, and 6.
References
-
Immature granulocytes index as early marker of SepsisInternational Journal of Laboratory Hematology 41:392–396.https://doi.org/10.1111/ijlh.12990
-
Improved flagging rates on the sysmex XE-5000 compared with the XE-2100 reduce the number of manual film reviews and increase laboratory productivityAmerican Journal of Clinical Pathology 136:309–316.https://doi.org/10.1309/AJCPDLR4KGKAFW4W
-
Laboratory diagnosis of Anemia: are the old and new red cell parameters useful in classification and treatment, how?International Journal of Laboratory Hematology 38:123–132.https://doi.org/10.1111/ijlh.12500
-
White blood cell counts: reference methodologyClinics in Laboratory Medicine 35:11–24.https://doi.org/10.1016/j.cll.2014.10.007
-
Reactive lymphocytes in patients with COVID-19British Journal of Haematology 189:844.https://doi.org/10.1111/bjh.16690
-
Evaluation of nucleated red blood cells in the peripheral blood of hematological diseasesClinical Chemistry and Laboratory Medicine 50:357–360.https://doi.org/10.1515/CCLM.2011.766
-
Hematologic parameters in patients with COVID-19 infectionAmerican Journal of Hematology 95:E131.https://doi.org/10.1002/ajh.25774
-
Clinical characteristics of coronavirus disease 2019 in ChinaNew England Journal of Medicine 382:1708–1720.https://doi.org/10.1056/NEJMoa2002032
-
Hematologic, biochemical and immune biomarker abnormalities associated with severe illness and mortality in coronavirus disease 2019 (COVID-19): a meta-analysisClinical Chemistry and Laboratory Medicine 58:1021–1028.https://doi.org/10.1515/cclm-2020-0369
-
Reticulated platelets: analytical aspects and clinical utilityClinical Chemistry and Laboratory Medicine 52:1107–1117.https://doi.org/10.1515/cclm-2014-0165
-
Prediction for progression risk in patients with COVID-19 pneumonia: the CALL scoreClinical Infectious Diseases 71:1393–1399.https://doi.org/10.1093/cid/ciaa414
-
The association between severe COVID-19 and low platelet count: evidence from 31 observational studies involving 7613 participantsBritish Journal of Haematology 190:e29–e33.https://doi.org/10.1111/bjh.16817
-
A summary of the diagnostic and prognostic value of hemocytometry markers in COVID-19 patientsCritical Reviews in Clinical Laboratory Sciences 57:415–431.https://doi.org/10.1080/10408363.2020.1774736
-
Detection of activated neutrophils by reactive oxygen species production using a hematology analyzerJournal of Immunological Methods 463:122–126.https://doi.org/10.1016/j.jim.2018.10.004
-
Laboratory abnormalities in patients with COVID-2019 infectionClinical Chemistry and Laboratory Medicine 58:1131–1134.https://doi.org/10.1515/cclm-2020-0198
-
Prognostic value of interleukin-6, C-reactive protein, and procalcitonin in patients with COVID-19Journal of Clinical Virology 127:104370.https://doi.org/10.1016/j.jcv.2020.104370
-
Prognostic value of C-reactive protein in patients with COVID-19Clinical Infectious Diseases : An Official Publication of the Infectious Diseases Society of America 71:2174–2179.https://doi.org/10.1093/cid/ciaa641
-
Pathological inflammation in patients with COVID-19: a key role for monocytes and macrophagesNature Reviews Immunology 20:355–362.https://doi.org/10.1038/s41577-020-0331-4
-
Diagnosis and treatment of adults with Community-acquired pneumonia an official clinical practice guideline of the american thoracic society and infectious diseases society of americaAmerican Journal of Respiratory and Critical Care Medicine 200:e45–e67.https://doi.org/10.1164/rccm.201908-1581ST
-
Leukoerythroblastic reaction in a patient with COVID-19 infectionAmerican Journal of Hematology 95:999–1000.https://doi.org/10.1002/ajh.25793
-
Rapid screening of COVID-19 patients using white blood cell scattergrams, a study on 381 patientsBritish Journal of Haematology 190:718–722.https://doi.org/10.1111/bjh.16943
-
Sepsis affects most routine and cell population data (CPD) obtained using the sysmex XN-2000 blood cell analyzer: neutrophil-related CPD NE-SFL and NE-WY provide useful information for detecting SepsisInternational Journal of Laboratory Hematology 37:190–198.https://doi.org/10.1111/ijlh.12261
-
Platelet-to-lymphocyte ratio is associated with prognosis in patients with coronavirus disease-19Journal of Medical Virology 92:1533–1541.https://doi.org/10.1002/jmv.25767
-
Radiological findings from 81 patients with COVID-19 pneumonia in Wuhan, China: a descriptive studyThe Lancet Infectious Diseases 20:425–434.https://doi.org/10.1016/S1473-3099(20)30086-4
-
Abnormalities of peripheral blood system in patients with COVID-19 in Wenzhou, ChinaClinica Chimica Acta 507:174–180.https://doi.org/10.1016/j.cca.2020.04.024
-
Hematological findings and complications of COVID-19American Journal of Hematology 95:834–847.https://doi.org/10.1002/ajh.25829
-
The diagnostic and predictive role of NLR, d-NLR and PLR in COVID-19 patientsInternational Immunopharmacology 84:106504.https://doi.org/10.1016/j.intimp.2020.106504
-
Thrombocytopenia and its association with mortality in patients with COVID-19Journal of Thrombosis and Haemostasis 18:1469–1472.https://doi.org/10.1111/jth.14848
-
Temporal changes in immune blood cell parameters in COVID-19 infection and recovery from severe infectionBritish Journal of Haematology 190:33–36.https://doi.org/10.1111/bjh.16847
-
Can we predict the severity of coronavirus disease 2019 with a routine blood test?Polish Archives of Internal Medicine 130:400–406.https://doi.org/10.20452/pamw.15331
Article and author information
Author details
Funding
Sysmex Europe GMBH provided free of charge reagents for the study. No monetary payments were made to any of the investigators. Joachim Linssen, Jarob Saker and Marion Münster are full-time employees of Sysmex Europe GMBH and Andre van der Ven has an ad hoc consultancy agreement with Sysmex Europe GMBH.
Acknowledgements
We would like to thank all patients involved in this study, as well as doctors, nurses, and researchers working together to fight against COVID-19. We also acknowledge Claudia Wienefoet of Sysmex Netherlands/Belgium for her support to the Dutch and Belgian participating centres and Sysmex Europe throughout the duration of this study. Sysmex Europe provided free of charge reagents to the study sites. No other funders were involved in this study.
Ethics
Human subjects: The study was reviewed by all participating centre ethics committees with approval granted in Italy (Registration Number 54/20) and Belgium (Registration Number 3002020000105) and exemption in the Netherlands, with need for informed consent waived by all.
Copyright
© 2020, Linssen et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 34,856
- views
-
- 1,905
- downloads
-
- 42
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Medicine
- Microbiology and Infectious Disease
- Epidemiology and Global Health
- Immunology and Inflammation
eLife has published articles on a wide range of infectious diseases, including COVID-19, influenza, tuberculosis, HIV/AIDS, malaria and typhoid fever.
-
- Cell Biology
- Medicine
The IncRNA Malat1 was initially believed to be dispensable for physiology due to the lack of observable phenotypes in Malat1 knockout (KO) mice. However, our study challenges this conclusion. We found that both Malat1 KO and conditional KO mice in the osteoblast lineage exhibit significant osteoporosis. Mechanistically, Malat1 acts as an intrinsic regulator in osteoblasts to promote osteogenesis. Interestingly, Malat1 does not directly affect osteoclastogenesis but inhibits osteoclastogenesis in a non-autonomous manner in vivo via integrating crosstalk between multiple cell types, including osteoblasts, osteoclasts, and chondrocytes. Our findings substantiate the existence of a novel remodeling network in which Malat1 serves as a central regulator by binding to β-catenin and functioning through the β-catenin-OPG/Jagged1 pathway in osteoblasts and chondrocytes. In pathological conditions, Malat1 significantly promotes bone regeneration in fracture healing. Bone homeostasis and regeneration are crucial to well-being. Our discoveries establish a previous unrecognized paradigm model of Malat1 function in the skeletal system, providing novel mechanistic insights into how a lncRNA integrates cellular crosstalk and molecular networks to fine tune tissue homeostasis, remodeling and repair.