key: cord-1006910-o90xvg4n
authors: Linssen, Joachim; Ermens, Anthony; Berrevoets, Marvin; Seghezzi, Michela; Previtali, Giulia; van der Sar-van der Brugge, Simone; Russcher, Henk; Verbon, Annelies; Gillis, Judith; Riedl, Jürgen; de Jongh, Eva; Saker, Jarob; Münster, Marion; Munnix, Imke CA; Dofferhof, Anthonius; Scharnhorst, Volkher; Ammerlaan, Heidi; Deiteren, Kathleen; Bakker, Stephan JL; Van Pelt, Lucas Joost; Kluiters-de Hingh, Yvette; Leers, Mathie PG; van der Ven, Andre J
title: A novel haemocytometric COVID-19 prognostic score developed and validated in an observational multicentre European hospital-based study
date: 2020-11-26
journal: eLife
DOI: 10.7554/elife.63195
sha: a490b7449b6af9ae342876c5761c9539e57e8bb4
doc_id: 1006910
cord_uid: o90xvg4n

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.

COVID-19 spans a wide clinical spectrum from asymptomatic to severe pneumonia with multiple organ failure , 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 interpatient 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.

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. Source data 1. Source data for clinical severity and outcome by age of all patients. * Occurrence of comorbidities is shown as a relative frequency expressed as a percentage (patients with a comorbidity divided by total number of patients in whom the presence or absence of comorbidities was recorded) as 2 of the participating study centres did not document the presence or absence of comorbidities. † Absent -comorbidities confirmed to be absent. ‡ Present -one or more comorbidities confirmed to be present. § Body mass index (BMI) measurements were not undertaken as part of this study. It was also up to the discretion of the physician to document obesity as a comorbidity. As this study was undertaken in the early phase of the pandemic, it may not have been common knowledge to all attending physicians that obesity is significant contributor to adverse outcomes in COVID-19. ¶ Other includes gastrointestinal disorders, musculoskeletal disorders, endocrine disorders, lipid disorders, haemochromatosis psoriasis, malnutrition, and Down syndrome.

The online version of this article includes the following source data for Table 1 : Source data 1. Source data for basic demographic characteristics of COVID-19 PCR confirmed patients enrolled for prognostic score development. 

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

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 ). Note: The expanded name of the haemocytometric parameters is shown in brackets in italics. The unit of measure for each parameter is shown in square brackets. A more detailed explanation of the haemocytometric parameters is provided in Table 7 .

Abbreviations: na, not applicable; ns, not significant.

The online version of this article includes the following source data for Table 2 : Source data 1. Source data for demographic and first haemocytometric data available for patients at hospital presentation or within 3 days after admission.

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 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 Source data 1. Source data for number of measurements for each day of hospitalisation that were available per patient group for the trends of lymphocyte-related parameters over 14 days of hospitalisation. Source data 2. Source data for trends of absolute lymphocyte count over 14 days of hospitalisation in critical illness (CI) and non-critical (NC) patients. Source data 3. Source data for trends of neutrophil-to-lymphocyte ratio over 14 days of hospitalisation in critical illness (CI) and non-critical (NC) patients. Source data 4. Source data for trends of AS-LYMPH%/L over 14 days of hospitalisation in critical illness (CI) and non-critical (NC) patients. Source data 5. Source data for trends of RE-Lymph minus AS-LYMPH as a percentage of total lymphocytes over 14 days of hospitalisation in critical illness (CI) and non-critical (NC) patients.

throughout, with values slightly higher in the CI group, but only until day 5, after which the NC and CI groups overlap ( Figure 7B ).

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) . Source data 1. Source data for number of measurements for each day of hospitalisation that were available per patient group for the trends of neutrophil-related parameters over 14 days of hospitalisation. Source data 2. Source data for trends of absolute neutrophil count over 14 days of hospitalisation in critical illness (CI) and non-critical (NC) patients. Source data 3. Source data for trends of NEUT-RI over 14 days of hospitalisation in critical illness (CI) and non-critical (NC) patients. Source data 4. Source data for trends of immature granulocyte count over 14 days of hospitalisation in critical illness (CI) and non-critical (NC) patients. Source data 5. Source data for trends of immature granuclocyte-to-lymphocyte ratio over 14 days of hospitalisation in critical illness (CI) and non-critical (NC) patients.

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 .

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 Figure 5 . Trends of monocyte parameters over 14 days of hospitalisation in critical illness (CI) and non-critical (NC) patients. Note: 14 days of hospitalisation refers to day 0 (day of admission) plus the first 13 days after admission. The normal reference range is depicted by the area between the dotted horizontal lines. Vertical bars indicate standard error of the mean (SEM). (A) Absolute monocyte count (MONO), (B) reactive monocytes as a percentage of monocytes (RE-MONO%/M). The number of sample measurements available per day for the trend analysis for the parameters plotted per patient group are shown in Figure 5 -source data 1. The online version of this article includes the following source data for figure 5:

Source data 1. Source data for number of measurements for each day of hospitalisation that were available per patient group for the trends of monocyte parameters over 14 days of hospitalisation. Source data 2. Source data for trends of absolute monocyte count over 14 days of hospitalisation in critical illness (CI) and non-critical (NC) patients. Source data 3. Source data of trends of reactive monocytes as a percentage of total monocyte count over 14 days of hospitalisation in critical illness (CI) and non-critical (NC) patients.

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. The number of sample measurements available per day for the trend analysis for the parameters plotted per patient group are shown in Figure 6 source data 1.

The online version of this article includes the following source data for figure 6:

Source data 1. Source data for number of measurements for each day of hospitalisation that were available per patient group for the trends of red blood cell-related parameters over 14 days of hospitalisation. Source data 2. Source data for trends of corrected haemoglobin values over 14 days of hospitalisation in critical illness (CI) and non-critical (NC) patients. Source data 3. Source data for trends of absolute reticulocyte count over 14 days of hospitalisation in critical illness (CI) and non-critical (NC) patients. Source data 4. Source data for trends of Delta-He values over 14 days of hospitalisation in critical illness (CI) and non-critical (NC) patients. Source data 5. Source data for trends of nucleated red blood cell counts over 14 days of hospitalisation in critical illness (CI) and non-critical (NC) patients.

(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 ).

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 Source data 1. Source data for number of measurements for each day of hospitalisation that were available per patient group for the trends of platelet parameters over 14 days of hospitalisation. Source data 2. Source data for trends of platelet count over 14 days of hospitalisation in critical illness (CI) and non-critical (NC) patients. Source data 3. Source data for trends of platelet-to-lymphocyte count ratio over 14 days of hospitalisation in critical illness (CI) and non-critical (NC) patients. Source data 4. Source data for trends of absolute immature platelet count (IPF#) over 14 days of hospitalisation in critical illness (CI) and non-critical (NC) patients. Source data 5. Source data for trends of immature platelet fraction (IPF%) over 14 days of hospitalisation in critical illness (CI) and non-critical (NC) patients.

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 and platelet-to-lymphocyte ratio (PLR) for development cohort incorporating all measurements over the initial 14-day period of hospitalisation. The number of measurements for each day of hospitalisation that were available per patient group are shown in Figure 8 -source data 1. There were markedly fewer measurements for the second week, notably in the NC group which may have contributed to bias. The online version of this article includes the following source data for figure 8:

Source data 1. Source data for number of measurements for each day of hospitalisation that were available per patient group for prognostic score prediction of clinical severity in the development cohort. Source data 2. Source data for development cohort prognostic score 14-day hospitalisation time horizon (day of admission plus the first 13 days thereafter) comparing non-critical (NC) and critical illness (CI) groups. Source data 3. Source data for development cohort prognostic score 7-day hospitalisation time horizon comparing outcomes (recovered without ICU, recovered with ICU or died). Source data 4. Source data for ROC curves to assess the capability of prediction of critical illness disease progression of the prognostic score, absolute lymphocyte count (LYMPH), neutrophil-to-lymphocyte ratio (NLR), absolute monocyte count (MONO), platelet count (PLT), and platelet-to-lymphocyte ratio (PLR) for development cohort incorporating all measurements over the initial 14-day period of hospitalisation.

lymphoplasmacytoid lymphocytes [Linssen et al., 2007] ) in the prognostic score. Increased lymphocyte activation has been documented in COVID-19 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) . Note: For the primary variables, one point = value above the cut-off value for the best AUC; two points = value above the cut-off value for the best AUC and !80% specificity; three points = value above the cut-off value for the best AUC and >90% specificity; and four points = value above the cut-off value for the best AUC and >95% specificity. The cut-off values for the secondary variables were chosen exclusively based on observed extremes of values in critical disease, with the maximum award of 1 point per variable.

The prognostic score values were calculated automatically using a pre-set algorithm, using the above cut-off values to assign points per individual parameter. The aim is to have the formula for the calculation incorporated into the Laboratory Information System software in use in individual laboratories.

Abbreviations: IGLR, immature granulocytes-to-lymphocyte ratio; NLR, neutrophil-to-lymphocyte ratio; AS-LYMPH%/L, antibody synthesising lymphocytes as a proportion of lymphocytes; RE-MONO%/M, reactive monocytes as a proportion of monocytes; DELTA-He, difference in haemoglobinisation of reticulocytes and red blood cells; NRBC, absolute nucleated red blood cell count; HGB, haemoglobin; HYPO-He, percentage of red blood cells that are hypochromic; PLT, platelet count; IPF, absolute immature platelet count; N/L 2 , neutrophil-to-lymphocyte squared ratio.

* Unit of measure for IG (absolute immature granulocyte count) is x 10 3 /mL. † Unit of measure for LYMPH (absolute lymphocyte count) is x 10 3 /mL. ‡ Unit of measure for RE-MONO (absolute reactive monocyte count) is x 10 3 /mL. § Unit of measure for LYMPH% (percentage lymphocyte count). ¶ Unit of measure for RET (absolute reticulocyte count) is x 10 3 /mL. ** Unit of measure for HGB (haemoglobin) is g/dL. † † Unit of measure for Micro-R is %.

Erythropoietic changes have been reported, mostly showing low HGB levels 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 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 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 Note: Recovered refers to patients that survived and were discharged from hospital. By definition, 'Mild' patients are not hospitalised and therefore only day 0 prognostic score values are available.

* Median score refers to the median value obtained for the haemocytometric COVID-19 prognostic score for the patient group for patients with the same length of hospitalisation on the day of measurement. Values in brackets represent the 95% CI for the median. † Sample size was too small to calculate 95% CI.

The online version of this article includes the following source data for Table 4 : Source data 1. Source data for development cohort prognostic score values by day based on clinical severity group, initial management, and outcome. Table 5 continued on next page 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 The online version of this article includes the following source data for 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 Table 6 . Mann-Whitney test for significance of the difference in prognostic score between critical illness (CI) and non-critical (NC) patients. 'md' refers to the Hodges-Lehmann median difference in prognostic score values between the NC and CI groups in each age category.

* Patients that were enrolled into the study from the two hospitals that did not report on comorbidities were excluded from the analysis. † Sample size was too small for analysis in this age group segregated by comorbidities. ‡ Patients older than 84 years and with reported comorbidities were the only group where the difference between CI and NC groups was not significant.

The online version of this article includes the following source data for Table 6 : Source data 1. Mann-Whitney test for significance of the difference in prognostic score between critical illness (CI) and non-critical (NC) patients. Figure 9 . Impact of age and presence of comorbidities on prediction of disease severity. (A) Box and whisker plots of prognostic score values for NC and CI groups segregated by age. The prognostic score can predict severity independent of age, therefore potentially assisting in identifying young patients at risk for severe disease progression as well as older patients not at risk. (B) Box and whisker plots of prognostic score values for NC and CI groups segregated by comorbidities in the 75-84-year-old group, as an illustrative example. The prognostic score is significantly higher in patients with severe disease progression independent of the presence of comorbidities. Please refer to Table 6 for more detailed information on all age groups. The online version of this article includes the following source data for figure 9:

Source data 1. Impact of age on prediction of disease severity.

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 were not met a Figure 10 . Flow chart showing inclusion and exclusion of validation cohort patient. Abbreviations: NC, non-critical patient group; CI, critical illness patient group. (a) The following criteria were used to ensure selection of only those patients for whom the primary presentation at hospital was related to COVID-19: Emergency department location on day 0, provisional diagnosis of pneumonia, if admitted, with admission to a general ward, internal medicine, ICU or anaesthesia (critical care).

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 12 33 2 A) B) Figure 11 . Haemocytometric COVID-19 prognostic score prediction of clinical severity in the validation cohort. (A) Validation cohort prognostic score 14-day hospitalisation time horizon comparing non-critical (NC) and critical illness (CI) groups, (B) ROC curve comparisons of prognostic score and NLR over 14 days. The prognostic score AUC (0.838, 95%CI 0.809-0.864) was significantly higher (P<0.0001) than the NLR AUC (0.673,95% CI 0.637-0.707). The number of measurements for each day of hospitalisation that were available per patient group are shown in Figure 11 -source data 1.

Overall, there were relatively few measurements per day for the NC group which has contributed to greater variance per time point. The online version of this article includes the following source data for figure 11:

Source data 1. Source data for number of measurements for each day of hospitalisation that were available per patient group for prognostic score prediction of clinical severity in the validation cohort. Source data 2. Source data for haemocytometric COVID-19 prognostic score prediction of clinical severity in the validation cohort over 14 day hospitalisation time horizon comparing non-critical (NC) and critical illness (CI) groups. Source data 3. Source data for ROC curve comparisons of prognostic score and NLR over 14 days.

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.

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.

Study design, sample size, and participants

Whilst 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).

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

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

Step 5:

Step 6:

Step 7:

Step 4: Parameters identified in step 1. All samples from days 0-3 Determination of prognostic score cut-off value to best discriminate between NC and CI during first 3 days of hospitalisation (using all full profile (CBC-DIFF-RET) samples from days 0-3)

Validation of prognostic score and cut-offs using an independent patient cohort Prognostic score 14-day time horizon plot (using all full profile samples from days 0-13) Identification of additional parameters to improve discrimination between CI and NC groups (using all samples from days 0-3) (Maximum 1 point per additional parameter identified)

All parameters All samples from days 0-3 Prognostic score parameters

Step 1:

Step 2:

Univariate analysis with regard to critical or fatal outcome 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.0001); 2) >20% outside reference range (CI group); 3) only dominant parameter selected if ≥2 eligible parameters are interdependent.

Step 3:

ROC curve analysis of individual parameters. AUC values used to select best cut-off value to discriminate between NC and CI groups.

Points per parameter assigned based on fulfilment of 1) precondition 2) optimal cut-off value 3) with extra points for increasing specificity and likelihood ratio (Range 0-4 points per variable) Prognostic Score Figure 12 . Flow chart illustrating steps involved in prognostic score development. Abbreviations: NC, non-critical patient group; CI, critical illness patient group; ROC, receiver operating characteristics (curve); AUC, area under the curve.

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

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

draft, Writing -review and editing; Marion Mü nster, Formal analysis, Validation, Visualization, Writing -original draft, Writing -review and editing; Andre J van der Ven, Conceptualization, Formal analysis, Validation, Visualization, Methodology, Writing -original draft, Project administration, Writing -review and editing

Simone van der Sar-van der Brugge http://orcid.org/0000-0001-6462-075X Lucas Joost Van Pelt http://orcid.org/0000-0001-5538-1806

Andre J van der Ven https://orcid.org/0000-0003-1833-3391

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. . Supplementary file 1. Basic demographic characteristics of COVID-19 PCR confirmed patients enrolled for prognostic score development demographic data by hospital.

. Supplementary file 2. Novel parameters of different manufacturers in relation to possible adaptability of the haemocytometric COVID-19 prognostic score.

. Transparent reporting form . Reporting standard 1. TRIPOD checklist.

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. 

Immature granulocytes index as early marker of Sepsis

Improved flagging rates on the sysmex XE-5000 compared with the XE-2100 reduce the number of manual film reviews and increase laboratory productivity

Laboratory diagnosis of Anemia: are the old and new red cell parameters useful in classification and treatment, how?

Estimation of risk factors for COVID-19 mortality -preliminary results

White blood cell counts: reference methodology

Risk factors of fatal outcome in hospitalized subjects with coronavirus disease 2019 From a Nationwide Analysis in China

Reactive lymphocytes in patients with COVID-19

Evaluation of nucleated red blood cells in the peripheral blood of hematological diseases

Hematologic parameters in patients with COVID-19 infection

A clinical risk score to identify patients with COVID-19 at high risk of critical care admission or death: an observational cohort study

Clinical characteristics of coronavirus disease 2019 in China

Hematologic, biochemical and immune biomarker abnormalities associated with severe illness and mortality in coronavirus disease 2019 (COVID-19): a meta-analysis

Reticulated platelets: analytical aspects and clinical utility

Prediction for progression risk in patients with COVID-19 pneumonia: the CALL score

The association between severe COVID-19 and low platelet count: evidence from 31 observational studies involving 7613 participants

The role of biomarkers in diagnosis of COVID-19 -A systematic review

A summary of the diagnostic and prognostic value of hemocytometry markers in COVID-19 patients

Risk stratification of patients admitted to hospital with covid-19 using the ISARIC WHO clinical characterisation protocol: development and validation of the 4C mortality score

Detection of activated neutrophils by reactive oxygen species production using a hematology analyzer

Identification and quantification of high fluorescence-stained lymphocytes as antibody synthesizing/secreting cells using the automated routine hematology analyzer XE-2100

Thrombocytopenia is associated with severe coronavirus disease 2019 (COVID-19) infections: a meta-analysis

Laboratory abnormalities in patients with COVID-2019 infection

Prognostic value of interleukin-6, C-reactive protein, and procalcitonin in patients with COVID-19

Longitudinal characteristics of lymphocyte responses and cytokine profiles in the peripheral blood of SARS-CoV-2 infected patients

ACP risk grade: a simple mortality index for patients with confirmed or suspected severe acute respiratory syndrome coronavirus 2 disease (COVID-19) during the early stage of outbreak in Wuhan

Prognostic value of C-reactive protein in patients with COVID-19. Clinical Infectious Diseases : An Official Publication of the Infectious Diseases Society of

Nucleated red blood cells as predictors of mortality in patients with acute respiratory distress syndrome (ARDS): an observational study

Pathological inflammation in patients with COVID-19: a key role for monocytes and macrophages

Diagnosis and treatment of adults with Community-acquired pneumonia an official clinical practice guideline of the american thoracic society and infectious diseases society of america

Leukoerythroblastic reaction in a patient with COVID-19 infection

Revisiting the white blood cell count: immature granulocytes count as a diagnostic marker to discriminate between SIRS and Sepsis-a prospective, observational study

Rapid screening of COVID-19 patients using white blood cell scattergrams, a study on 381 patients

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 Sepsis

A novel diagnostic algorithm equipped on an automated hematology analyzer to differentiate between common causes of febrile illness in southeast asia

Platelet-tolymphocyte ratio is associated with prognosis in patients with coronavirus disease-19

Host susceptibility to severe COVID-19 and establishment of a host risk score: findings of 487 cases outside Wuhan

Radiological findings from 81 patients with COVID-19 pneumonia in Wuhan, China: a descriptive study

Lower mortality of COVID-19 by early recognition and intervention: experience from jiangsu province

Characteristics and prognostic factors of disease severity in patients with COVID-19: the Beijing experience

Abnormalities of peripheral blood system in patients with COVID-19 in Wenzhou, China

Hematological findings and complications of COVID-19

Coronavirus disease 2019 in elderly patients: characteristics and prognostic factors based on 4-week follow-up

Clinical course and outcome of 107 patients infected with the novel coronavirus, SARS-CoV-2, discharged from two hospitals in Wuhan

Delta-He, Ret-He and a new diagnostic plot for differential diagnosis and therapy monitoring of patients suffering from various Disease-Specific types of Anemia

Risk factors associated with acute respiratory distress syndrome and death in patients with coronavirus disease

Prediction models for diagnosis and prognosis of covid-19 infection: systematic review and critical appraisal

Prediction of criticality in patients with severe Covid-19 infection using three clinical features: a machine learning-based prognostic model with clinical data in Wuhan

The diagnostic and predictive role of NLR, d-NLR and PLR in COVID-19 patients

Thrombocytopenia and its association with mortality in patients with COVID-19

Temporal changes in immune blood cell parameters in COVID-19 infection and recovery from severe infection

Association of radiologic findings with mortality of patients infected with 2019 novel coronavirus in Wuhan

Can we predict the severity of coronavirus disease 2019 with a routine blood test?

Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study

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. 

* Standard parameters available in the complete blood count (CBC), differential (DIFF) or reticulocyte measurement (RET) channels across multiple manufacturer haematology analyser platforms. † Depending on region, the immature platelet count (IPF# and IPF%) are obtained either from the RET or PLT-F channel. PLT-F refers to fluorescent optical platelet measurement.

Joachim Linssen, Jarob Saker, Marion Mü nster: is a permanent employee of Sysmex Europe GMBH who provided free of charge study reagents to the study centres. Andre J van der Ven: has an ad hoc consultancy agreement with Sysmex Europe GMBH who provided free of charge study reagents to the study centres. The other authors declare that no competing interests exist.

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.