LOX-1-Expressing Immature Neutrophils Identify Critically-Ill COVID-19 Patients at Risk of Thrombotic Complications

Background Lymphopenia and the neutrophil/lymphocyte ratio may have prognostic value in COVID-19 severity. Objective We investigated neutrophil subsets and functions in blood and bronchoalveolar lavage (BAL) of COVID-19 patients on the basis of patients’ clinical characteristics. Methods We used a multiparametric cytometry profiling based to mature and immature neutrophil markers in 146 critical or severe COVID-19 patients. Results The Discovery study (38 patients, first pandemic wave) showed that 80% of Intensive Care Unit (ICU) patients develop strong myelemia with CD10−CD64+ immature neutrophils (ImNs). Cellular profiling revealed three distinct neutrophil subsets expressing either the lectin‐like oxidized low‐density lipoprotein receptor‐1 (LOX‐1), the interleukin-3 receptor alpha (CD123), or programmed death-ligand 1 (PD-L1) overrepresented in ICU patients compared to non-ICU patients. The proportion of LOX-1- or CD123-expressing ImNs is positively correlated with clinical severity, cytokine storm (IL-1β, IL-6, IL-8, TNFα), acute respiratory distress syndrome (ARDS), and thrombosis. BALs of patients with ARDS were highly enriched in LOX-1-expressing ImN subsets and in antimicrobial neutrophil factors. A validation study (118 patients, second pandemic wave) confirmed and strengthened the association of the proportion of ImN subsets with disease severity, invasive ventilation, and death. Only high proportions of LOX-1-expressing ImNs remained strongly associated with a high risk of severe thrombosis independently of the plasma antimicrobial neutrophil factors, suggesting an independent association of ImN markers with their functions. Conclusion LOX-1-expressing ImNs may help identifying COVID-19 patients at high risk of severity and thrombosis complications.


INTRODUCTION
Since the first reports of an outbreak of a severe acute respiratory syndrome caused by coronavirus 2 (SARS-CoV-2) in China in December 2019 (1,2), the coronavirus disease (COVID-19) has grown to be a global public health emergency, with cases of COVID-19 around the world reaching 5.5 million and deaths from the disease standing at more than 90 000 as of May 2021 (for up-to-date data, see https://www.who.int/emergencies/diseases/ novel-coronavirus-2019/situation-reports). The SARS-CoV-2 infection is characterized by a range of symptoms including fever, cough, fatigue, and myalgia in the majority of cases, and occasional headache and diarrhea (1,3). Among reported cases, approximatively 80% present with a mild condition, 13% with a serious condition, and 6% in a critical state requiring intensive care; the latter is associated with a fatality rate of 2-8% of reported cases (4). Some severe cases of COVID-19 progress to acute respiratory distress syndrome (ARDS), which accounts for high mortality related to damage of the alveolar lumen. Numerous patients with ARDS secondary to COVID-19 develop lifethreatening thrombotic complications (5).
Previous coronavirus-related infections have been characterized by the onset of a cytokine storm (6). The inflammatory response measured at both cellular and molecular levels could represent a strong prognostic signature of the disease. The cytokine storm is an uncontrollable inflammatory response leading to viral sepsis, ARDS, respiratory failure, shock, organ failure, or death (7,8). Strong predictive markers are still missing for these complications.
Older age, neutrophilia, and organ and coagulation dysfunction are the major risk factors associated with the development of ARDS and progression to death (9). In addition, serum concentrations of both pro-and antiinflammatory cytokines, including IL-6, TNFa, and IL-10, were increased in the majority of severe cases and were markedly higher than those of moderate cases, suggesting that the cytokine storm might be associated with disease severity (10,11) and leading the way to the development of potential immune-modulatory treatments (3,12). The cytokine storm is associated with a massive influx of innate immune cells, namely neutrophils and monocytes, which could worsen lung injury. The accumulation of innate effectors, in particular neutrophils, in the bronchoalveolar lavage (BAL) of COVID-19 patients correlates with pro-inflammatory cytokines and reflects the outcome of patients (13). However, our knowledge of a particular innate immune actor (i.e., neutrophils and molecular mechanisms) during severe COVID-19 is incomplete.
Increasing clinical data indicate that the neutrophil/lymphocyte ratio (NLR) is a powerful predictive and prognostic indicator of severe COVID-19 (14)(15)(16). Lymphopenia, neutrophilia, and high NLR are associated with a more severe viral infection (14,17). It is well recognized that, in an emergency response to severe infection, the bone marrow releases mature and immature neutrophils (ImNs). ImNs have been shown to be enriched in granule antimicrobial, cytotoxic, and neutrophil extracellular trap (NET)forming proteins (18). More recently, NETs were shown to be associated with COVID-19 severity and microthrombi formation (19)(20)(21).
We previously identified two new CD10−CD64+ neutrophil subsets, expressing either PD-L1 or CD123, that were specific to bacterial sepsis (22). In addition to these markers, previous work showed that LOX-1 is an important mediator of inflammation and neutrophil dysfunction in sepsis and cancers (23,24). To test the hypothesis of a virally driven neutrophil profile, we developed a multiparametric neutrophil profiling strategy based on known neutrophil markers to distinguish COVID-19 phenotypes in patients in critical or severe condition.

Study Participants
Fresh blood samples from 38 (pilot study, Mar-Apr 2020) and 118 (validation study, Sep-Nov 2020) consecutive adult patients with COVID-19 referred to the Department of Internal Medicine 2, Department of Infectious Diseases and Intensive Care Units of Pitie-Salpetrière University Hospital, Paris, were included. The diagnosis of COVID-19 relied on SARS-CoV-2 carriage on the nasopharyngeal swab, as confirmed by real-time reverse transcription PCR analysis. In addition, blood samples were collected from healthy donors obtained from the French blood donation center. When ventilator-associated pneumonia was suspected, patients underwent fiber-optic bronchoscopy and BAL (n = 16) to sample distal respiratory secretions. Leukocyte phenotyping was performed on BAL when possible. Demographic and clinical characteristics are detailed in Tables 1 and 2, respectively.

Study Approval
The studies were conducted in accordance with the Declaration of Helsinki and the International Conference on Harmonization (ICH) Good Clinical Practice (GCP) Guideline and approved by the relevant regulatory and independent ethics committees. In accordance with current French law, informed written consent was obtained from patients or relatives. The studies were registered and approved by the local ethical committee of  Sorbonne Universite/Assistance Publique -Hopitaux de Paris for standard hospitalized patients (N°2020-CER2020-21) and ICU patients (N°CER-2020-31).

Flow Cytometry
One hundred microliters of fresh whole blood collected at admission to the hospital-on anticoagulant citrate-dextrose solution (ACD) for patients in the ICU or on EDTA for patients in standard hospitalization-was stained with a mix of monoclonal antibodies. Samples were diluted in Brilliant  Table S1). One milliliter of BD FACS lysing solution 1X (BD Biosciences) was directly added to the cells to lyse red blood cells, which were incubated for 20 min, centrifuged, and washed with PBS. Leukocytes were resuspended in PBS before analysis with a BD LSRFortessa X-20 (BD Biosciences). FlowJo software 10.0 (FlowJo LLC, Ashland, OR, USA) was used for analysis of marker expression on neutrophils. One hundred microliters of whole blood were stained with a fluorescence minus one (FMO) mix missing antibodies targeting CD123, LOX-1, and PD-L1, in order to determine the threshold of expression of these markers. BAL leukocyte phenotyping was performed similarly after filtration, two wash procedures of BAL cells, and staining with the same antibody mix. Acquired data were normalized and analyzed using the OMIQ platform (https://www.omiq.ai).
To identify neutrophil subsets and visualize all cells in a 2D map where position represents local phenotypic similarity, we used a dimensionality reduction tool: the opt-SNE  Simoa (single molecule array) HD-1 analyzer was used for ultrasensitive multiplex immunodetection of cytokines as described in methods section. Potential association between serum cytokine levels and marker-expressing neutrophil frequencies was evaluated by Spearman's correlation (one-tailed), with significance defined by a p-value < 0.05: *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001; ns, not significant.
implementation of t-SNE. Neutrophils (40 000 events) were randomly taken from the sample for the unsupervised analysis. Cells were also grouped in phenotypically homogeneous clusters using the FlowSOM algorithm.  Table S1).

Data Presentation and Statistical Analysis
Statistical analyses of the immunological data and graphic representations were performed with Prism 9 (GraphPad Software Inc.). A nonparametric Mann-Whitney test was used for group comparisons, and one-way ANOVA tests with a Tukey's adjustment was used for multiple group comparison tests. The association between variables was evaluated using Spearman's correlation (one-tailed), with significance defined by a p-value of <0.05. Survival curves were compared with a logrank (Mantel-Cox) test and were considered statistically significant with a p of <0.05. HR values with 95% CI were computed. ROC curves were created using Prism 9. PCAs were performed with R-software 3.3.1.

Increased Proportions of Circulating Immature Neutrophils Expressing Either CD123 or LOX-1 in Critical COVID-19 Patients Are Associated With COVID-19 Severity and Thromboembolic Complications
We designed a first discovery study with 38 individuals and analyzed their neutrophil phenotypes, comparing them with those of patients admitted to the ICU (n = 24) or not (n = 14), within the first day following their admission to the ICU or hospitalization units (see Table 1 and  Table 1). Patients' characteristics confirmed previously published data, with a notably high prevalence of obese COVID-19 patients. The main differences between ICU and non-ICU patients reflected case severity, with a high proportion of patients with lung lesions (observed as ground-glass opacities (GGOs) on chest CT scans) requiring invasive mechanical ventilation (IMV) and resulting in high in-hospital mortality. Within 3 h of drawing blood, we performed whole blood immunostaining using a previously published panel designed to precisely evaluate immature circulating neutrophils (22). Neutrophils were automatically identified and visualized using a dimensionality reduction tool (opt-SNE for optimal implementation of t-distributed stochastic neighbor embedding (t-SNE)) to define an imprint for each sample group ( Figure 1A). Under unsupervised classification, neutrophils from ICU patients were organized in the upper left quadrant of the map, whereas those from the non-ICU patients ended up in the upper right quadrant. This analysis allowed delimitation of two main subsets of neutrophils based on the expression of CD10 and CD64 markers: the ICU-abundant upper left area was composed of neutrophils with mid-to-low expression of CD10 and high expression of CD64, and the non-ICU-abundant upper right area was composed of high-to-mid expression of CD10 and high-to-mid expression of CD64. We next determined whether the identified neutrophil signature would be confirmed by conventional analysis undertaken by experts. Neutrophils were identified with the CD15 neutrophil marker, while excluding prototypical markers of eosinophils and monocytes, respectively CRTH2 and CD14 (see Supplementary Figures S1A and S1B for representative patient samples). An expert-gating strategy confirmed the high abundance of CD10 − CD64 + ImNs among ICU patients compared with non-ICU patients ( Figure 1B). We did not observe any correlation between ImN proportion and age, obesity, hypertension, diabetes, and smoking (data not shown).
We next compared the expression of CD123, LOX-1, and PD-L1 surface molecules, formerly known to be dysregulated in sepsis (22). All three of them were barely co-expressed on neutrophils (Supplementary Figure S1C), which led to the identification of three distinct ImN subpopulations. Subsets of neutrophils expressing either CD123 or LOX-1 were more abundant in ICU than in non-ICU patients ( Figure 1C), unlike those expressing PD-L1, for which overabundance did not reach statistical significance (p = 0.1). Principal component analysis (PCA) revealed that CD123-, LOX-1-, or PD-L1-expressing ImNs contributed independently to the disease severity of patients as evaluated by SAPS II (see Supplementary Figure S2A). In addition, Nonparametric Mann-Whitney test was used to compare differences in cellular abundance of neutrophil subsets between groups, with significance defined by a p-value < 0.05: **p < 0.01; ***p < 0.001 and ns, not significant.
PCA confirmed that the disease severity in ICU patients was associated with the expression of CD123, LOX-1, and PD-L1 on ImNs and distinct patterns of cytokines (see Supplementary Figure  S2B and Table 2). We identified three profiles: (a) patients with high LOX-1-expressing ImN proportions and high IL-1b, IL-6, IL-8, and TNFa serum levels; (b) patients with CD123 and PD-L1 expression and IL-18, IL-22, and IFNg secretion; and (c) patients with a lower severity score associated with high type-1 interferon levels. Thus, ImN subsets expressing either CD123, LOX-1, or PD-L1 may define specific profiles of severity associated with high levels of cytokines. Because the mortality in COVID-19 cases is associated with the virally driven cytokine storm, especially in patients with comorbidities (1), we sought to correlate ImN subsets with severe symptoms such as ARDS (patients receiving intermittent mandatory ventilation (IMV)) or thrombosis ( Figures 1D, E,  respectively). The requirement for IMV or the occurrence of thromboembolic events was associated with higher proportions of ImN expressing CD123 or LOX-1, unlike those expressing PD-L1 (p = 0.1). No differences in the proportions of ImN subsets could be detected between discharged and deceased patients (see Supplementary Figure S2C). These data suggested that increased proportions of circulating ImNs expressing either CD123 or LOX-1 in critical COVID-19 patients are associated with COVID-19 severity and thromboembolic complications.

Subsets of Immature Neutrophils Expressing LOX-1 Infiltrate Lung
A pulmonary immune environment during critical COVID-19 infection is one of the major features of disease complications. We thus sought neutrophil subsets in the bronchoalveolar lavages (BALs) when available and compared them to blood samples from patients and healthy donors (Figure 2). With a opt-SNE algorithm, BAL neutrophils were identified in the upper right quadrant of the map; blood neutrophils from COVID-19 patients were more central; and blood neutrophils from healthy donors (HD) were organized in the lower left quadrant of the map (Figure 2A). Automatic clustering using major (CD15, CD10, CD16, CD64) and specific (CD123, LOX-1, PD-L1) neutrophil markers split neutrophil signatures into positive and negative subpopulations for each marker ( Figure 2B and Supplementary Figure S3). This unsupervised analysis allowed the identification of nine clusters, representing three main subsets of neutrophils ( Figure 2C): (a) the mature neutrophils (MatN or subset 8) with high expression of CD15 and CD10 and low expression of CD64; (b) the ImNs with high expression of CD15 and CD64 and low expression of CD10(ImN or subset 6); and (c) the activated neutrophils (ActN or subset 3) with high expression of CD15, CD10, and CD64. Expression of LOX-1 and PD-L1 was spotted on ActN (subset 1 and 2) whereas LOX-1, CD123, and PD-L1 expression was observed on three independent ImN subsets (subsets 4,7 and 5). MatNs were abundant in healthy donor blood, and both ActNs and ImNs were abundant in the blood of COVID-19 patients ( Figure 2D). If ImNs expressing LOX-1, PD-L1, or CD123 represented a few percent of COVID-19 blood neutrophils ( Figure 2D; see also Figure 1C), these subsets were much more present in patient BALs with ImNs expressing LOX-1 being the major subset, representing about 40% of total neutrophils. A profusion of ImNs in COVID-19 BAL was associated with massive production of myeloperoxidase (MPO) and neutrophil elastase (ELA) (Figure 2E), two antimicrobial and cytotoxic proteins known to be highly concentrated in the azurophilic granule of ImN.
These data revealed that ImNs, preferentially those expressing LOX-1, infiltrate bronchoalveolar space in the lungs of COVID-19 patients, where they release their cytotoxic content, suggesting a potential role in disease severity.
Immature Neutrophil Subsets Expressing CD123, LOX-1, or PD-L1 Are Correlated With Clinical Severity, but Only the LOX-1+ Subset Proportion at Entry Is Strongly Associated With Higher Risk of Thrombosis COVID-19 patients from the validation study were segregated into three groups based on severity of disease at the time of admission: 18 in mild, 19 in severe, and 51 in critical condition ( Table 3 and Supplementary Results). The proportion of CD123-, LOX-1-, and PD-L1-expressing ImNs correlated positively with severity ( Figure 3A). Interestingly, abundances of all three ImN subsets were associated with patient death ( Figure 3B) and with patients requiring IMV (see Supplementary Figure S4), but only the LOX-1-expressing ImN subset was associated with thromboembolic events ( Figure 3C), confirming our results in the pilot study (see Figure 1E). We segregated our COVID-19 patients into two groups based on the median proportion of each ImN subset and compared their relative risk of ARDS, thrombosis, and death using a Cox proportional hazards model with other variables (age, gender, hypertension, obesity, and diabetes) ( Figure 3D). Patients with a high abundance of ImN subsets were at higher risk of ARDS requiring IMV and of death. Hypertension and diabetes were risk factors of survival but not for ARDS. Again, only patients with a high abundance of LOX-1-expressing ImNs were at higher risk of thromboembolic complications (HR, 5.99; 95% CI, 2.02-17.81; p = 0.007).
These data from the validation study confirmed that ImN subsets expressing either CD123, PD-L1, or LOX-1 were associated with COVID-19 severity, but only LOX-1 expression remained associated with thromboembolic complications.

Immature Neutrophil Subsets and Plasma Levels of MPO and ELA Are Independent Markers of COVID-19 Severity
Because ImNs that recently emigrated from the bone marrow are enriched in granule antimicrobial, cytotoxic, and NET-forming proteins, we measured plasma levels of MPO, ELA, and MPO-DNA complexes representing NET formation in three groups of patients on the basis of severity at the time of admission: 23 in mild, 22 in severe, and 63 in critical condition. MPO and ELA plasma levels were significantly associated with disease severity ( Figure 4A), whereas ELA-DNA complexes were not. In addition, MPO and ELA levels at hospital admission were also significantly increased among COVID-19 patients who later died ( Figure 4B). MPO-DNA complex levels were not associated with survival. There was no association between MPO, ELA, or MPO-DNA complexes and thromboembolic events ( Figure 4C). PCA ( Figure 4D) that combined abundances of ImN subsets and plasma levels of neutrophil microbicidal proteins revealed at least two independent patient profiles at risk of severity: those with high proportions of ImN subsets and those with high plasma levels of MPO, ELA, and NET. To further evaluate the ability of LOX-1 neutrophil markers to segregate patients with thrombosis, we plotted a receiver operating characteristic (ROC) curve ( Figure 4E). The ROC analysis of these abundances indicated the optimal threshold yielding the best separation of the two groups of patients with optimal sensitivity and specificity. The AUC was 0.89 for LOX-1 ImN abundance

DISCUSSION
Stratification of patients by using biomarkers remains an unmet need in COVID-19 patient care. Here, we sought to identify innate immune cellular signatures that may help predict the outcome of COVID-19 patients with severe symptoms. Conventional whole blood flow cytometry identified classical hallmarks of severe infection, such as neutrophilia and myelemia, but also revealed three novel neutrophil subsets able to distinguish between patients requiring ICU and those who do not. All three immature CD10 − CD64 + neutrophil subsets expressing either CD123, LOX-1, or CD123 strongly correlated with severity scores commonly used in clinical practice: SAPS II and the Sequential Organ Failure Assessment (SOFA). They were also associated with distinct cytokine profiles. We found that high blood proportions of LOX-1+ neutrophils could be an important predictive signature of thromboembolic events, while D-dimer levels, an important negative predictor of thrombosis, were found to be high in a vast majority of patients. The strength of this study is that the main results were observed in the pilot study and confirmed in the validation study, representing samples from more than 150 patients collected over the course of two independent pandemic waves in France. Previously, we reported that the proportion of CD123expressing ImNs was correlated with bacterial sepsis severity (22). Here, we showed that the expression of CD123 on CD10 − CD64 + neutrophils was associated with COVID-19 severity (higher SOFA or SAPS II scores, longer length of stay, and higher risk of death). Our results also suggested that both CD123 and its cognate ligand, the IL-3 cytokine, play an important role in sepsis. High IL-3 plasma levels have been associated with lung inflammation, lung injury, and high mortality rates in an animal model, as well as in humans (25,26). IL-3 neutralization and anti-CD123 treatment improved mice outcome by decreasing inflammation and mortality rates (25). IL-3 promotes emergency myelopoiesis, exacerbating proinflammatory cytokine secretion and, consequently, systemic inflammation, organ dysfunction, and death. However, we did not observe similar results in COVID-19, as we reported an inverse correlation between IL-3 levels and SOFA scores. The association between CD123 expression on ImN and high serum levels of IL-17, IL-22, and IFNg, to the best of our knowledge, has never been reported and may reveal a yet unidentified link between innate and adaptive immune responses. These findings open the way to new therapeutic opportunities aiming to control the excessive inflammation induced by SARS-CoV-2 infection.
Single-cell technology approaches to assessing dysregulation in the myeloid compartment have identified subsets of dysfunctional neutrophils in COVID-19 patients (27,28). PD-L1-expressing neutrophils were detected only in patients with severe COVID-19. The authors revealed that this subset displayed surface markers and gene expression profiles reminiscent of myeloid-derived suppressor cells (MDSCs), suggesting that they may be involved in immune regulation by suppressing other immune cell activity. PD-L1-expressing MDSCs were previously described in HIV (29) and in cancer patients (30); they were also suggested to be involved in the immunosuppressive activity of T cell functions. Whether these PD-L1 neutrophils are dysfunctional or ImNs remains unclear. Emergency granulopoiesis in response to infections is well-documented (31). It mobilizes ImNs equipped with stocks of innate defense armory that are packaged into different granule subsets (32). ProNETosis gene signatures were observed in both pro-(MPO, ELA, and PRTN3) and pre-(PADI4) neutrophils and were associated with COVID-19 severity. Other reports have shown that NETosis increased in COVID-19 patients with ARDS and thrombotic complications (33,34). We measured high levels of MPO and ELA in patients with critical conditions, but we did not observe any association between NETosis and patients' clinical characteristics in the validation study. A potential bias of our study is that patients' samples were obtained on the very first day of hospital admission and NETosis may have only been at its early stages. Active anticoagulant therapy administered to every patient may have prevented it (35).  , or critical (n = 51) clinical condition observed in patients from the validation study. One-way ANOVA with a Tukey's adjustment was used to compare the three groups (critical vs severe or critical vs mild), with significance defined as follows: *p < 0.05; and ***p < 0.001. (B) Box plots (min to max distribution) of the abundance of CD10−CD64+ neutrophil subsets among discharged (n = 74) and deceased (n = 14) patients. Nonparametric Mann-Whitney test was used to compare differences between groups, with significance defined by a p-value < 0.05: * **p < 0.01; ***p < 0.001; ns, not significant. (C) Box plots (min to max distribution) of the abundance of CD10−CD64+ neutrophil subsets among patients without (n = 75) or with (n = 12) thromboembolic complications. Nonparametric Mann-Whitney test was used to compare differences in cellular abundance of neutrophil subsets between groups, with significance defined by a p-value < 0.05: *p < 0.05; **p < 0.01. (D) Forest plots comparing hazard ratio (log-rank test) for survival, invasive mechanical ventilation (IMV), and thrombotic events (thrombosis) in 118 patients according to median frequency of ImN subsets (reference group with ImN subset frequencies below the median) and comorbidities (reference group with no comorbidities or age below 61). Log-rank (Mantel-Cox) test was used to compare HR between groups, with significance defined by a p-value < 0.05: *p < 0.05; **p < 0.01; and ***p < 0.001. LOX-1 was originally discovered as a receptor for oxLDL (36). oxLDL binding to LOX-1 results in oxLDL internalization and cellular accumulation, which triggers intracellular signaling processes leading to pro-apoptotic, pro-oxidant, and proinflammatory pathways, causing cell dysfunction associated with atherosclerosis and increased cardiovascular disease (CVD) risk (37). Binding of endothelial LOX-1 by ligands induces superoxide generation, inhibits nitric oxide production, enhances endothelial adhesiveness for leukocytes, and induces expression of chemokines (38).
In addition, LOX-1 binds other ligands with poor or no structural similarities to lipoproteins (for review (39) including phospholipids, proteins, bacteria and cells. It may have worked as host-defense molecule and still maintains some amino acid sequence homology with C-type lectin receptor.
LOX-1 expression and functions on neutrophils remain elusive. LOX-1 is barely detected on neutrophils during homeostasis, while its expression increases on neutrophils from human cancer patients (23) and in murine sepsis (24,40). LOX-1 deletion in a murine model of polymicrobial sepsis reduced IL-6 and TNFa levels in the blood and lungs, which enhanced bacterial clearance and prevented neutrophil activation (24). LOX-1 was identified as a marker on granulocytic MDSCs able to suppress T cell activity (23) and produced high amount ROS. Recently, the work of Cabrera et al. revealed that ImN of severe COVID-19 patients had significantly more LOX-1+ cells as compared to mild cases and healthy controls (41) but T cell suppressive activity could not be evaluated. However, LOX-1 is mostly acknowledged for its role in atherosclerosis. LOX-1 is a class E scavenger receptor contributing to the formation of atherosclerotic plaques by promoting endothelial cell activation, macrophage foam cell formation, and smooth muscle cell migration and proliferation. LOX-1 activation induces NF-kB activation, leading to pro-inflammatory cytokine release, endoplasmic reticulum stress, and ROS production, which could damage the microenvironment (38,39).
In this study, LOX-1 expression on ImNs seems to be detrimental for patients, as it was associated with the secretion of several pro-inflammatory cytokines, such as IL-6, IL-1b, and TNFa, with severity and a higher risk of thrombosis. In severe cases of COVID-19, the integrity of the lung is compromised by an exaggerated immune response leading to ARDS (17,42). Mechanisms contributing to microcirculation disorders in sepsis are capillary leakage, leukocyte adhesion and infiltration, and intravascular coagulation, leading to thrombus formation. Over the course of systemic inflammatory diseases such as sepsis, the microenvironment is highly oxidative, leading to an increase of oxidized low-density lipoprotein (oxLDL) in plasma, which triggers LOX-1 overexpression through a positive feedback loop. In physiological conditions, the increase of LOX-1 expression, especially by endothelial cells, leads to an increase of LDL uptake by vessel walls, which activates the specific Oct-1/SIRT1 thrombosis protective pathway (43,44). In this study, we observed an increase in vascular thrombotic events among individuals displaying a high frequency of immature LOX-1+ neutrophils. Whether thrombosis in COVID-19 patients results from functionally diverted neutrophils expressing LOX-1 or from its expression on endothelial and smooth muscle cells remains to be elucidated. We also observed a slight correlation between LOX-1-expressing ImNs and D-dimer levels (Spearman's correlation, r = 0.42 p = 0.023). Ddimer level is a marker of activation of coagulation and fibrinolysis. However, in our study, the predictive score of D-dimers for thromboembolic events was not significant (ROC test, AUC = 0.79, p = 0.003, n = 116) compared to LOX-1-expressing neutrophil abundance in the blood (AUC = 0.89, p < 0.0001, n = 118). These results suggest that the high LOX-1 expression by ImNs might be a useful tool for predicting thromboembolic events among critically ill COVID-19 patients. The overexpression of LOX-1 may also be found in other cell types that might trigger the prothrombotic ERK1/2 pathway. The polymorphic LOX-1 gene is also intensively associated with increased susceptibility to myocardial diseases. LOX-1 should thus be considered a potential target for therapeutic interventions.
The first step in viral entry of SARS-CoV-2 into host cells is binding of a viral envelope protein, known as the SPIKE (S) protein, to a cell host cell membrane receptor-enzyme protein, namely the angiotensin-converting enzyme 2 (ACE2) protein. ACE2 converts Angiotensin II (ANGII) the main effector of the renin-angiotensin system (RAS) to the vasodilator ANG1-7, which binds to MAS, a receptor that signals to reduce the activity of the ANGII receptor AT1R. Because ACE2 functions as both a SARS-CoV-2 receptor and a RAS regulator, the potential interference between of SARS-CoV-2 infection and RAS-controlled pathophysiology such as in heart, kidney and liver diseases was suspected. It has been suggested that the activation of the RAS could pre-dispose the patients with comorbidities to severe COVID-19 (45). Recently, patients with pre-existing cardiovascular diseases have a higher mortality rate and are five-times more likely to experience severe COVID-19 symptoms. LOX-1 is one of the major receptors of ox-LDL and thus an additional player in RAS (46). Dyslipidemia up-regulates and activates AT1R, and RAS activation in turn up-regulates and activates LOX-1, and then facilitates uptake of ox-LDL into endothelial cells. Activation of both receptors may predispose patients to severe COVID-19 infection. In this context, it is worth noting that plasma ox-LDL and LOX-1 messengers were increased in children with Kawasaki Disease (KD). Of note, recent clinical data suggest that SARS-CoV-2 may cause a pediatric multisystem inflammatory syndrome reminiscent of Kawasaki Disease (KD) in children (47).
In conclusion, we outline new measurable potential biomarkers of COVID-19 severity among immature circulating neutrophils: CD123, LOX-1, and PD-L1 surface markers. These markers are significantly correlated with disease severity, and even more so with thromboembolic events in the case of LOX-1.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

ETHICS STATEMENT
The studies involving human participants were reviewed and approved by local ethical committee of Sorbonne Universite/ Assistance Publique -Hôpitaux de Paris for standard hospitalized patients (N°2020-CER2020-21) and ICU patients (N°CER-2020-31). The patients/participants provided their written informed consent to participate in this study.

AUTHOR CONTRIBUTIONS
BC and CC designed the study and performed data analysis. LA, PR, NG, CP, OB, and AB performed experimental work and compiled and analyzed the data. PQ was responsible for clinical data mining and analysis. KD and DS performed cytokine dosage and analysis. MM, ZA, JM, C-EL, GG, and AG provided patient sample access. GG, AG, BC, and CC provided financial support. BC, LA, PQ, and CC wrote the manuscript. All authors contributed to the article and approved the submitted version.

FUNDING
The study was supported by the Fondation de France, « Tous unis contre le virus » framework Alliance (Fondation de France, AP-HP, Institut Pasteur) in collaboration with the Agence nationale de la Recherche (ANR Flash COVID19 program ICOVID and ANR COV7 program Neutrosets), by the Programme Hospitalier de Recherche Clinique PHRC-20-0375 COVID-19, and by the SARS-CoV-2 Program of the Faculty of Medicine at Sorbonne University (I-COVID programs). The program was supported by AG2R-Reǵion Ile de France. LA and PR are recipients of postdoctoral fellowships from the European Union's Horizon 2020 research and innovation program funding, under grant agreement No. 681137.