key: cord-0886643-dxl2fgmh
authors: Meizlish, Matthew L.; Pine, Alexander B.; Bishai, Jason D.; Goshua, George; Nadelmann, Emily R.; Simonov, Michael; Chang, C-Hong; Zhang, Hanming; Shallow, Marcus; Bahel, Parveen; Owusu, Kent; Yamamoto, Yu; Arora, Tanima; Atri, Deepak S.; Patel, Amisha; Gbyli, Rana; Kwan, Jennifer; Won, Christine H.; Dela Cruz, Charles; Price, Christina; Koff, Jonathan; King, Brett A.; Rinder, Henry M.; Wilson, F. Perry; Hwa, John; Halene, Stephanie; Damsky, William; van Dijk, David; Lee, Alfred I.; Chun, Hyung J.
title: A neutrophil activation signature predicts critical illness and mortality in COVID-19
date: 2021-03-09
journal: Blood Adv
DOI: 10.1182/bloodadvances.2020003568
sha: ecdceff9c7f1e5304997f4bcd0ca10fc89ca4754
doc_id: 886643
cord_uid: dxl2fgmh

Pathologic immune hyperactivation is emerging as a key feature of critical illness in COVID-19, but the mechanisms involved remain poorly understood. We carried out proteomic profiling of plasma from cross-sectional and longitudinal cohorts of hospitalized patients with COVID-19 and analyzed clinical data from our health system database of more than 3300 patients. Using a machine learning algorithm, we identified a prominent signature of neutrophil activation, including resistin, lipocalin-2, hepatocyte growth factor, interleukin-8, and granulocyte colony-stimulating factor, which were the strongest predictors of critical illness. Evidence of neutrophil activation was present on the first day of hospitalization in patients who would only later require transfer to the intensive care unit, thus preceding the onset of critical illness and predicting increased mortality. In the health system database, early elevations in developing and mature neutrophil counts also predicted higher mortality rates. Altogether, these data suggest a central role for neutrophil activation in the pathogenesis of severe COVID-19 and identify molecular markers that distinguish patients at risk of future clinical decompensation.

As the death toll from COVID-19 exceeds 2.1 million worldwide, it remains a pressing concern to understand how the disease causes such a wide spectrum of clinical outcomes. For most patients, COVID-19 manifests as an upper respiratory tract infection that is self-limited. However, the progression of COVID-19 in a large subset of patients to respiratory distress, multiorgan failure, and death has resulted in an enormous global impact. A number of studies have highlighted important roles for monocytes and macrophages in severe COVID-19, 1,2 but our knowledge of the immunologic drivers of critical illness is otherwise limited.

To achieve a deeper understanding of the immunologic phenotypes of COVID-19 across the spectrum of disease severity, we performed proteomic profiling of blood obtained from hospitalized patients with COVID-19 and used a machine learning algorithm to define the biomarkers that best discriminate between critically ill patients and those with milder disease. We discovered a unique neutrophil activation signature composed of neutrophil activators (granulocyte colony-stimulating factor [G-CSF] and interleukin-8 ) and neutrophil-derived effectors (resistin [RETN] , lipocalin-2 [LCN2], and hepatocyte growth factor [HGF]), which had the greatest power of all measured biomarkers to identify critically ill patients. These effector proteins were highly transcriptionally enriched in a developing neutrophil population that was recently identified specifically in critically ill COVID-19 patients and were strongly correlated with absolute neutrophil count (ANC) in our patient cohorts. 3, 4 This signature of neutrophil activation was predictive of in-hospital mortality and, most compellingly, was elevated at the time of hospital admission in patients who only later progressed to critical illness, thus preceding and predicting the onset of critical illness. We identify neutrophil activation as a defining feature of severe COVID-19 that occurs before the onset of critical illness, implicating neutrophils as a central player in the pathogenesis of severe COVID-19 and highlighting opportunities for clinical prediction and therapeutic intervention.

We first conducted a study of 49 adult patients admitted to Yale-New Haven Hospital between 13 and 24 April 2020 with a confirmed diagnosis of COVID-19 via polymerase chain reaction (cross-sectional cohort). The protocol was approved by the Yale University Institutional Review Board (IRB; IRB 2000027792). Plasma was analyzed from 40 patients treated in a medical intensive care unit (ICU) and 9 patients treated on standard medical floors (non-ICU) in our hospital. Admission to the ICU was based on established ICU admission guidelines (supplemental Table 1 ). For this cohort, blood was collected at a single time point during each patient's hospitalization (supplemental Figure 1 ). Blood samples from 13 additional asymptomatic, nonhospitalized controls were also analyzed after signed consents were obtained for a separate IRB-approved protocol (IRB 1401013259).

We also analyzed blood samples obtained longitudinally on days 1 (within 24 hours), 4, and 7 of hospitalization from a separate cohort of 23 consecutive adult patients admitted for treatment of laboratory-confirmed COVID-19 between 23 and 28 May 2020 who remained hospitalized until at least day 4 (longitudinal cohort). These patients were identified prospectively based on real-time electronic chart review of hospital admissions, and analysis was performed for consecutive patients who met the following criteria: (1) COVID-19-positive polymerase chain reaction test; (2) admission for the purpose of COVID-19 treatment; (3) patient identification in a time frame allowing for sample collection within 24 hours of admission; and (4) collection of at least 2 longitudinal samples (resulting in the exclusion of 8 patients discharged before day 4 sample collection). 

Blood was collected in 3.2% sodium citrate tubes and centrifuged at 2000g for 20 minutes at room temperature; the resulting plasma supernatant was frozen at 280°C and used for further testing. Biomarker profiling analyses were conducted at Eve Technologies (Calgary, Alberta, Canada). For the cross-sectional cohort, the following assays were performed: Human Cytokine 71-Plex, Human Complement Panels 1 and 2, Human SAA & ADAMTS13, and Human Adipokine 5-Plex. For the longitudinal cohort, the following assays were performed: Human Cytokine 48-Plex, Human Complement Panel 1, Human Adipokine 5-Plex, and Human MMP 9-Plex and TIMP 4-Plex. Five control samples were evaluated concurrently with the longitudinal cohort samples. Heatmaps were generated using concentrations of circulating biomarkers obtained from biomarker profiling analyses using Heatmapper as described. 5 For confirmation, RETN levels were also measured by enzyme-linked immunosorbent assay (ELISA) (category no. DRSN00; R&D Systems) following a fivefold sample dilution per the manufacturer's protocol. Linear regression analysis was performed to assess the relationship between ELISA-based concentrations and values reported by Eve Technologies based on multiplex biomarker profiling.

Values for biomarker profiles in the cross-sectional cohort were logtransformed. A pseudocount of half the minimum observed nonzero value per biomarker was added to each observed zero value before log 10 transformation. IL-6 was excluded from modeling because 38 of 40 patients in the ICU subset of the cross-sectional cohort had received IL-6 receptor blockade before blood collection. The remaining biomarker values were used in a principal component analysis (PCA) and a random forest classifier using the scikit-learn Python package. 6 Data were partitioned into 66-33 train-validation split (training cohort: 26 ICU patients, 7 non-ICU patients, 3 controls; validation cohort: 14 ICU patients, 2 non-ICU patients, 5 controls), and all biomarkers were then used to predict ICU status. A maximal tree depth of 10 was used with the minimal costcomplexity set to 0.02; all other parameters were set to their defaults. Feature importance was assessed using mean decrease in impurity. Data restricted to 5 biomarkers of interest with high feature importance were then used in a separate random forest classifier, using the same methods.

Processed count matrices with deidentified metadata and embeddings from a single-cell RNA sequencing (scRNAseq) dataset published by Wilk et al 3 were downloaded from the COVID-19 Cell Atlas (https://www.covid19cellatlas.org/#2iki20) hosted by the Wellcome Sanger Institute. This dataset was further analyzed using Seurat (version 3.0, https://satijalab.org/seurat). Clusters in this object were renamed based on cell type, using UMAPs and provided heatmaps. A new metadata column, ARDSstatus, was added to the object, which grouped samples into Control, COVID-19-non-acute respiratory distress syndrome (ARDS), and COVID-19-ARDS categories. Violin and feature plots were generated to examine specific genes of interest. Cells identified as developing neutrophil were made a subset from the overall object to examine this population more closely.

To evaluate differences in mean age across groups, we used 1-way analysis of variance with post hoc Tukey's multiple comparisons tests in the cross-sectional cohort and Kruskal-Wallis test in the longitudinal cohort, because the samples in the former cohort but not the latter conformed to the normal distribution. For comparisons of sex proportions, we used x 2 tests in the cross-sectional cohort and Fisher's exact test in the longitudinal cohort. In both crosssectional and longitudinal cohorts, we examined differences in proportions of patients with comorbidities using Fisher's exact test.

To compare biomarker values, we used 2-sample tests. Because most sample distributions did not satisfy Anderson-Darling and D'Agostino-Pearson tests of normality, we used the unpaired 2tailed Mann-Whitney U test. For samples that conformed to the normal distribution, we used unpaired 2-sided Student t tests with Welch correction for unequal variances, where applicable; where appropriate, P values were corrected for multiple comparisons using the false discovery rate procedure with the discovery rate Q of 5%. P , .05 was considered significant.

Correlation coefficients between ANC and each biomarker were computed using Spearman's rank method. For Kaplan-Meier survival analyses, we used median values of evaluated variables as cutpoints to classify patients into high and low groups, based on whether values for a particular biomarker or cell count were above or below the cutpoint. We then used log-rank test to compare the 2 groups for each evaluated variable. To determine P values for cell type-specific enrichment of genes in scRNAseq data, we used the Wilcoxon rank sum test. All statistical analyses were carried out using GraphPad Prism (v8.4.3; GraphPad Software, San Diego, CA), Stata (v16; StataCorp. College Station, TX), and R (v4; R Core Team, 2020).

Forty-nine adult COVID-19 patients (40 in the medical ICU and 9 in non-ICU units) and 13 nonhospitalized controls were included in our initial cross-sectional cohort (Table 1) . We performed multiplexed biomarker profiling of plasma to measure circulating concentrations of 78 proteins with immunologic functions ( Figure 1A ; supplemental Figure 2 ). PCA showed separation of the control, non-ICU COVID-19, and ICU COVID-19 samples, indicating that these circulating markers capture a spectrum of illness severity ( Figure 1B ).

To assess whether this profile of plasma biomarkers could distinguish between critically ill and noncritically ill patients, we applied a random forest machine learning prediction model. The model was trained on data from two-thirds of the patients, and its ability to predict ICU status was tested on the remaining one-third. The model accurately identified 14 of 14 patients in the ICU and 7 of 7 subjects not in the ICU (Figure 2A ). To avoid introducing bias based on treatment, the random forest model and PCA analyses excluded IL-6 because most (38 of 40) ICU patients received tocilizumab, which is known to increase IL-6 levels. 2 These results demonstrate that our proteomic profile provides a reliable circulating signature of critical illness in COVID-19.

To gain insight into the specific plasma proteins that may affect disease severity, we next examined the importance of each feature in the random forest prediction model. Five of the top 6 features contributing to the model were proteins related to neutrophil activation ( Figure 2B -C). Three of the top 4 features were RETN, LCN2, and HGF, each produced by neutrophils, stored in neutrophil secondary granules, and released on neutrophil activation. 7-12 The next 2 highest ranking features were IL-8 and G-CSF, which stimulate neutrophil chemotaxis and development, respectively. 13, 14 Next, we assessed whether these 5 proteins alone could discriminate between critically ill (ICU) and noncritically ill (non-ICU) patients. We repeated random forest modeling using only the markers related to neutrophil activation (RETN, LCN2, HGF, IL-8, and G-CSF) and found that circulating levels of this selective panel also accurately classified all test subjects into ICU and non-ICU categories ( Figure 2A) . Interestingly, each of the neutrophil activation markers showed greater discriminatory power for critical illness than previously described monocyte/macrophage markers (with the exception of MIG/CXCL9) that were also included in our proteomic panel ( Figure 1A ; Figure 2B ; supplemental Figure 2 ), 2,15-17 suggesting that neutrophils may play a central role in critical illness associated with COVID-19.

Although the neutrophil granule proteins RETN, LCN2, and HGF are produced in large quantities by neutrophils, they can also be secreted by other cell types. To determine whether increases in these proteins were associated with other indicators of neutrophil activity, we analyzed the correlation of the ANC at the time of blood draw with all circulating markers measured. We found that RETN, HGF, and LCN2 were the 3 proteins that most strongly correlated with ANC ( Figure 3A ), further suggesting that neutrophils are the likely common source of these proteins. We then analyzed data from a scRNAseq study of peripheral blood mononuclear cells (PBMCs) from patients with COVID-19. 3 We found that RETN and LCN2 transcripts were highly enriched in a circulating cell population designated as developing neutrophils (P , 10e 2300 for each), which was detected almost exclusively in severely ill patients with ARDS ( Figure 3B ; supplemental Figure 3 ). 3 Another highly enriched marker (P , 10e 2300 ) for this population was matrix metallopeptidase 8 (MMP8) or neutrophil collagenase ( Figure 3B ). Notably, constituents of neutrophil granules, like RETN, LCN2, and MMP8, are transcribed at earlier stages of neutrophil development, packaged into granules, and later released from mature neutrophils on degranulation. 18 This set of observations strongly supports the conclusion that neutrophils are the primary source of these circulating markers of critical illness in COVID-19.

In our initial cross-sectional cohort, blood was collected during the course of hospitalization, limiting our capacity to determine whether the observed neutrophil signature was a consequence of critical illness or whether it preceded its onset. To address this, we established a second longitudinal cohort of patients. Proteomic plasma profiling was conducted on blood samples collected serially, starting on day 1 (within 24 hours of hospital admission), from 23 consecutive patients admitted for treatment of confirmed COVID-19 who remained hospitalized for $4 days ( Figure 4A ; supplemental Figures 4 and 5) . To further test the role of neutrophil activation, we included a panel of matrix metalloproteinases (MMPs) and tissue inhibitor of metalloproteinases including MMP-8, a marker of the developing neutrophil population described above ( Figure 3B ). 19 Consistent with our findings in the cross-sectional cohort, RETN, HGF, and LCN2 in the day 1 samples of the longitudinal cohort showed strong correlations with ANC, and MMP-8 was the protein most highly correlated with ANC, further supporting a neutrophilic source of these proteins in COVID-19 (supplemental Figure 6 ).

Sixteen patients in the longitudinal cohort were initially admitted to lower-acuity non-ICU units, and 7 were admitted directly to the ICU (ICU-Admit) ( Table 2) . Of the patients admitted to non-ICU units, 9 remained in those units until discharge (non-ICU), whereas 7 were later transferred to the ICU (ICU-Transfer). Consistent with our observations in the cross-sectional cohort, ICU-Admit patients had significantly higher levels of all neutrophil activation markers than non-ICU patients on day 1 of hospitalization ( Figure 4B ; supplemental Figure 4b ). As predicted, the same pattern was observed for MMP-8, which, like RETN, HGF, and LCN2, is known to be stored and released from neutrophil secondary granules and is transcriptionally enriched in developing neutrophils in COVID-19 patients ( Figure 3B ). 20, 21 Remarkably, the levels of the neutrophil activation markers on day 1 of hospitalization were also significantly elevated in ICU-Transfer patients (who were at that time in non-ICU units), at levels comparable to ICU-Admit patients and well above levels in non-ICU patients ( Figure 4B ). In both ICU-Transfer and ICU-Admit patients, the neutrophil activation markers remained elevated from the day of admission to day 7 and did not change appreciably over time, whereas these levels in non-ICU patients remained stably low (supplemental Figure 4b) . Thus, despite the fact that patients in the ICU-Transfer group were not critically ill at the time of the day 1 blood draw, evidence of neutrophil activation was already present. The neutrophil activation signature identified patients who were primed for eventual transfer to the ICU before the onset of critical illness.

Previously reported markers of macrophage activation, including IL-6, IL-10, tumor necrosis factor-a (TNF-a), and MIG/CXCL9, were also elevated in ICU-Admit and ICU-Transfer patients compared with non-ICU patients on day 1 of hospitalization (supplemental Figure 5) . Notably, most of these markers (as well as IL-8 and G-CSF) were also significantly elevated in non-ICU patients compared with controls, whereas the neutrophil granule proteins (RETN, LCN2, HGF, MMP-8) that may more directly reflect neutrophil activation were not elevated in non-ICU patients compared with controls and were only significantly different in patients who later became critically ill ( Figure 4B ; supplemental Figure 5 ). This suggests that proteins released by neutrophil degranulation may be more specific to severe COVID-19 than elevations in traditional cytokines. Altogether, these longitudinal data offer the first evidence, to our knowledge, that neutrophil activation and monocyte/macrophage activation precede the onset of critical illness.

To further validate our findings, we conducted additional analyses using RETN, the factor most important in distinguishing critical illness in our random forest prediction model. Among patients initially admitted to non-ICU units (non-ICU, ICU-Transfer), those with day 1 RETN levels above the median value were much more likely to later require ICU transfer ( Figure 4C ). Moreover, patients with day 1 RETN levels above the median were significantly less likely to survive ( Figure 4D ). To confirm these results, we also measured RETN levels in the longitudinal cohort by a separate ELISA assay and found that concentrations measured by ELISA showed a strong linear relationship with those measured by multiplex profiling (R 2 5 0.9863). In addition to RETN, we found that higher levels on hospital day 1 of LCN2, G-CSF, IL-8, and MMP-8, as well as IL-6, IL-10, TNF-a, IL-1RA, and M-CSF, were significantly associated with mortality (supplemental Figure 7) . Meanwhile, we found that the neutrophil to lymphocyte ratio (NLR), previously described as a prognostic indicator in COVID-19, 16 was not significantly different between non-ICU and ICU-Transfer patients on day 1 of hospitalization (supplemental Figure 8a ), suggesting that NLR does not distinguish critical illness as clearly as the molecular markers, which may more directly reflect cellular activation. Among established clinical markers of severe COVID-19, C-reactive protein and D-dimer were elevated in ICU-Transfer and ICU-Admit patients compared with non-ICU patients (supplemental Figure 8a ), but none of these markers were significantly different between patients who did and did not survive their hospitalization (supplemental Figure 8b) , unlike the neutrophil and macrophage activation markers (supplemental Figure 7) .

Last, we explored whether these mechanistic insights regarding the role of neutrophils in severe COVID-19 could be validated in a large patient population hospitalized with COVID-19. Using a database of 3325 patients admitted to the Yale-New Haven Health System who tested positive for SARS-CoV-2 (DOM-CovX cohort; Table 3 ), we examined the first recorded values of immature granulocyte (IG) and neutrophil counts during the course of hospitalization. We used these early values to avoid potential confounding effects of hospitalization, such as therapeutic interventions and secondary infections. We found that in-hospital mortality was significantly higher among patients with elevated initial IG absolute count ( Figure 5A ), IG percent ( Figure 5B ), and ANC ( Figure 5C ). Interestingly, we did not find that initial absolute monocyte count was associated with increased mortality ( Figure 5D ). Although our molecular markers that would allow us to assess neutrophil activation are not available in this large cohort, these findings confirm that early signs of neutrophil development are associated with future mortality in one of the largest COVID-19 patient cohorts assessed to date. Transfer indicates patients who were admitted to a non-ICU unit and were transferred to an ICU unit during hospitalization; ICU-Admit indicates patients who were admitted directly to an ICU unit. Asterisks denote statistically significant differences between groups (*P , .05, **P , .01, ***P , .001, ****P , 0.0001). (C) Kaplan-Meier curve depicting the likelihood of ICU admission depending on resistin (RETN) levels on day 1, using the median of the group (non-ICU and ICU-Transfer) as the cutoff value (pg/mL).

(D) Kaplan-Meier curve depicting the likelihood of survival depending on resistin levels on day 1, using the median of the group (entire longitudinal cohort) as the cutoff value (pg/mL). ns, not significant.

Our analyses of multiple hospitalized patient cohorts with COVID-19 reveal a prominent signature of neutrophil development and activation in critically ill patients with COVID-19. Moreover, we demonstrate for the first time that elevations in circulating markers of neutrophil and macrophage activation precede the onset of critical illness, identifying noncritically ill patients who are at risk of becoming critically ill.

The host immune response is essential for pathogen resistance, but it can also cause costly tissue injury. 22 Neutrophils are often the first responders to infection and injury yet can also cause significant collateral damage. 23 Although important roles for monocytes and macrophages in severe COVID-19 have been well described, 15 a potential role for neutrophils has not received the same attention. However, suggestive evidence has begun to emerge, including the prognostic value of NLR, 24 the identification of a population of immature neutrophils in the blood of critically ill COVID-19 patients, and the observation of increased neutrophils in bronchoalveolar lavage. 1, 3, 4, 25 It has remained unclear, however, whether there is heightened neutrophil activation in severe COVID-19 and whether the observed innate immune response is a consequence or potential cause of critical illness. Using proteomic profiling and a machine learning prediction algorithm, we determined that several markers of neutrophil activation had the greatest discriminatory power for detecting critical illness, outperforming many other cytokines previously associated with severe COVID-19. 16, 17 This neutrophil signature includes G-CSF and IL-8, which can drive neutrophil activation and have been detected by other groups, with a recent study reporting an association of IL-8 levels with survival. 2, 26 Unlike in prior studies, we identify a group of proteins that are released from neutrophils themselves. RETN, LCN2, and MMP8 are all novel markers in COVID-19 (whereas HGF has been reported), 26 and the fact that they are among the strongest predictors of critical illness offers perhaps the most direct evidence to date that neutrophil activation is a hallmark of severe disease. These signs of conspicuous neutrophil activation may not have been detected in other analyses because neutrophil granule proteins are not included in standard cytokine panels.

One limitation of this study is that we did not directly demonstrate that neutrophils are the source of these proteins in COVID-19. However, multiple lines of evidence support this conclusion. First, RETN, LCN2, HGF, and MMP8 are all well-established products of neutrophils and are each known to be stored in the same subcellular compartment, from which they would be released together on neutrophil activation. [8] [9] [10] [11] 27, 28 Second, consistent with this, we found that each of these proteins was highly correlated with neutrophil counts and with one another. Third, by interrogating published scRNAseq data from blood samples of COVID-19 patients, where a developing neutrophil population was recently identified specifically in critically ill patients, 3 we found that RETN, LCN2, and MMP-8 were highly enriched in these neutrophil precursors and were not expressed at high levels in any other circulating cell population. Fourth, our findings are supported by other studies that have observed shifts in neutrophil populations in both the blood and lungs of patients with severe COVID-19. 1, 3, 4, 25 Importantly, we have not relied on the specificity of any single biomarker but rather the combination of all of these proteins (including MMP-8, which we measured as an additional test of our hypothesis) to draw the inference that their high levels are likely attributable to neutrophil activation. Nonetheless, we cannot exclude contributions from other cellular sources.

The temporal relationship between markers of immune activation and the onset of critical illness in COVID-19 has not been previously established. We reasoned that if neutrophils play a causal role in the development of critical illness, then evidence of neutrophil activation should be present before its onset. Remarkably, in our longitudinal cohort of patients, we found that the neutrophil activation signature was already elevated on day 1 of hospitalization in patients who were not critically ill at that time but later developed critical illness and required transfer to the ICU. Among all patients, day 1 elevations in neutrophil activation markers also predicted increased mortality. These observations additionally held for a few cytokines associated with monocyte/ macrophage activation (IL-6, IL-10, TNF-a), although these proteins were also markedly increased in non-ICU patients compared with controls, indicating that they are less specific to critical illness than the neutrophil signature. Altogether, these findings indicate that neutrophil activation specifically, and innate immune activation more broadly, precede the onset of critical illness and may play a key role in its pathogenesis.

Our findings suggest a model of neutrophil development and activation in the pathogenesis of severe COVID-19 ( Figure 6 ). We hypothesize that high levels of G-CSF stimulate emergency granulopoiesis to increase neutrophil production and that IL-8 (CXCL8) drives neutrophil migration into the lungs and perhaps other tissues. 13, 14, 29 There, neutrophils are activated and release RETN, LCN2, HGF, MMP8, and other proteins with antimicrobial and other inflammatory functions that also cause significant collateral damage to the lungs, vasculature, and other organs.

The occurrence of emergency granulopoiesis in COVID-19, which we propose in our model, is suggested by the presence of a developing neutrophil population specifically in severely ill patients. 3, 4 We find that these developing neutrophils transcribe the genes encoding the circulating proteins (RETN, LCN2, MMP8) that we detect in the blood of critically ill patients. This is consistent with the sequence of events described in the neutrophil literature: neutrophil secondary granule contents are transcribed by the transcription factor C/EBPe during the myelocyte stage of development, are packaged into granules, and are later released from mature neutrophils. 9, 18, 28, 30, 31 Typically, this neutrophil development takes place in the bone marrow, but during emergency granulopoiesis (driven by G-CSF, which we detect at elevated levels in patients with severe disease), immature granulocytes can be detected in the circulation ( Figure 6 ). 13 Our model accounts for the detection of developing neutrophils in the blood of only severely ill COVID-19 patients and is supported by our analysis of a hospitalwide dataset, where we find that the earliest counts of immature granulocytes predict mortality.

The neutrophil granule proteins that we detected may simply be markers of neutrophil activation or may themselves be detrimental to patients with COVID-19. RETN regulates production of multiple cytokines that are elevated in severely ill COVID-19 patients, including IL-6, IL-8, and TNF-a, 32 and is associated with increased endothelial permeability and thrombosis risk, another feature of severe COVID-19. 33 LCN2 is upregulated in multiple autoimmune diseases. 34, 35 Elevated levels of MMP-8 are associated with worsening clinical outcomes in pediatric patients with ARDS. 35 Neutrophils are known to contribute to acute lung injury in other contexts, including viral lung infections, and their role in COVID-19 may represent an exaggerated version of that pathophysiology. 36, 37 Interestingly, a recent study demonstrated a higher degree of neutrophil activation in severe COVID-19 than in influenza pneumonia. 38 Our findings raise essential questions. What are the initial triggers for neutrophil development, recruitment, and activation in COVID-19? Is there a threshold of injury to the lungs or vascular beds that stimulates neutrophil activation, 39, 40 setting in motion the cascade of events that propagate critical illness in COVID-19? Directly testing the causal role of neutrophil activation in COVID-19 will require either animal models that recapitulate neutrophil activation or clinical trials of targeted therapeutic interventions. Therapeutic strategies such as inhibition of G-CSF or IL-8 should be considered, but the potential benefits of these approaches will need to be weighed carefully against the risks of immunosuppression. Although multiple trials are underway for granulocyte-macrophage colony-stimulating factor (GM-CSF) inhibitors, this growth factor has a distinct functional role from G-CSF, and we observe an association of only G-CSF levels with disease severity in this study. 41, 42 To our knowledge, there are no ongoing clinical trials of G-CSF inhibitors in COVID-19. Overall, our study highlights a central role for neutrophil activation in the pathogenesis of severe COVID-19, which may help guide the development of new therapeutic strategies and more accurate predictive markers of severe disease. (955) 285 (1300) 109 (1452) 43 (1508) 15 (1532) 7 (1540) 1364 (0) 570 (647) 245 (909) 103 (1017) 49 (1058) 21 (1080) 7 (1092) IG ≤ median IG > median Number at risk (number censored) (908) 254 (1215) 109 (1338) 44 (1393) 14 (1421) 7 (1427) 1497 (0) 648 (694) 276 (994) 103 (1131) 48 (1173) 22 (1191) com.

Single-cell landscape of bronchoalveolar immune cells in patients with COVID-19

An inflammatory cytokine signature predicts COVID-19 severity and survival

A single-cell atlas of the peripheral immune response in patients with severe COVID-19

Severe COVID-19 is marked by a dysregulated myeloid cell compartment

Heatmapper: web-enabled heat mapping for all

BindsNET: a machine learning-oriented spiking neural networks library in Python

Human resistin promotes neutrophil proinflammatory activation and neutrophil extracellular trap formation and increases severity of acute lung injury

Resistin is stored in neutrophil granules being released upon challenge with inflammatory stimuli

Identification of murine and human XCP1 genes as C/EBP-epsilon-dependent members of FIZZ/ Resistin gene family

Differential role of neutrophils and alveolar macrophages in hepatocyte growth factor production in pulmonary fibrosis

Presence of a mobilizable intracellular pool of hepatocyte growth factor in human polymorphonuclear neutrophils

Molecular cloning and expression of a cDNA encoding NGAL: a lipocalin expressed in human neutrophils

Emergency granulopoiesis

Neutrophil recruitment and function in health and inflammation

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

Clinical features of patients infected with 2019 novel coronavirus in Wuhan

Sinai Immunology Review Project. Immunology of COVID-19: current state of the science

The ontogeny of a neutrophil: mechanisms of granulopoiesis and homeostasis

Human neutrophil collagenase. A distinct gene product with homology to other matrix metalloproteinases

Subendothelial infiltration of neutrophil granulocytes and liberation of matrix-destabilizing enzymes in an experimental model of human neo-intima

Human neutrophil collagenase expression is C/EBP-dependent during myeloid development

Disease tolerance as a defense strategy

The balancing act of neutrophils

Neutrophil-to-lymphocyte ratio as an independent risk factor for mortality in hospitalized patients with COVID-19

heightened innate immune responses in the respiratory tract of COVID-19 patients

Elevated plasma levels of selective cytokines in COVID-19 patients reflect viral load and lung injury

Lipocalin-2 from both myeloid cells and the epithelium combats Klebsiella pneumoniae lung infection in mice

CEBPE-mutant specific granule deficiency correlates with aberrant granule organization and substantial proteome alterations in neutrophils

Mechanisms of degranulation in neutrophils

Granulopoiesis and granules of human neutrophils

Regulation of neutrophil and eosinophil secondary granule gene expression by transcription factors C/EBP e and PU

Resistin induces expression of proinflammatory cytokines and chemokines in human articular chondrocytes via transcription and messenger RNA stabilization

Resistin: functional roles and therapeutic considerations for cardiovascular disease

The novel role of neutrophil gelatinase-B associated lipocalin (NGAL)/Lipocalin-2 as a biomarker for lupus nephritis

Early elevation of matrix metalloproteinase-8 and -9 in pediatric ARDS is associated with an increased risk of prolonged mechanical ventilation

Excessive neutrophils and neutrophil extracellular traps contribute to acute lung injury of influenza pneumonitis

Contribution of neutrophils to acute lung injury

Vascular neutrophilic inflammation and immunothrombosis distinguish severe COVID-19 from influenza pneumonia

Endotheliopathy in COVID-19-associated coagulopathy: evidence from a single-centre, cross-sectional study

Endothelial cell infection and endotheliitis in COVID-19

GM-CSF-based treatments in COVID-19: reconciling opposing therapeutic approaches

Granulocyte/macrophage colony-stimulating factor-deficient mice show no major perturbation of hematopoiesis but develop a characteristic pulmonary pathology

The authors thank all front-line staff taking care of our community of patients affected by COVID-19 and are grateful for the multidisciplinary collaboration at the Yale School of Medicine and Yale-New Haven Health System in the care of our community. This work was supported by a gift donation from Jack Levin to the Benign Hematology program at Yale, the National Institutes of Health, National Heart, Lung, and Blood Institute (HL142818