Rationale: COVID-19 displays distinct characteristics that suggest a unique pathogenesis. The objective of this study was to compare biomarkers of coagulopathy and outcomes in COVID-19 and non-COVID-19 patients with severe pneumonia.

Methods: Thirty-six non-COVID-19 and 27 COVID-19 non-immunocompromised patients with severe pneumonia were prospectively enrolled, most requiring intensive care. Clinical and biological characteristics (including plasma biomarkers of coagulopathy) were compared.

Results: At similar baseline severity, COVID-19 patients required mechanical ventilation (MV) for significantly longer than non-COVID-19 patients (p = 0.0049) and more frequently developed venous thrombotic complications (p = 0.031). COVID-19 patients had significantly higher plasma concentrations of soluble VCAM1 (sVCAM1) (5,739 ± 3,293 vs. 3,700 ± 2,124 ng/ml; p = 0.009), but lower levels of D-dimers, vWF-A2, sICAM1, sTREM1, VEGF, and P-selectin, compared to non-COVID-19 patients. Principal component analysis identified two main patterns, with a clear distinction between non-COVID-19 and COVID-19 patients. Multivariable regression analysis confirmed that sVCAM1 rising levels were independently associated with a longer duration of MV. Finally, we identified close correlations between sVCAM1 and some features of COVID-19 immune dysregulation (ie. CXCL10, GM-CSF, and IL-10).

Conclusion: We identified specific features of the coagulopathy signature in severe COVID-19 patients, with higher plasma sVCAM1 levels, that were independently associated with the longer duration of mechanical ventilation.

Clinical Trial Registration: ClinicalTrials.gov, identifier: NCT03505281.

Introduction

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for severe pneumonia and acute respiratory distress syndrome (ARDS) in its most severe form. Thrombotic complications have frequently been reported in coronavirus disease 2019 (COVID-19) patients. Coagulopathy was quickly identified as part of the pathogenesis of COVID-19 (1). Elevated levels of D-dimer and fibrin (ogen) degradation products are associated with outcomes in COVID-19 and suggest major coagulation activation (1–4). In addition, endothelial dysfunction and platelet activation are two further hallmarks of the most severe forms of COVID-19 (4–6). On the other hand, coagulopathy is also a classical complication of bacterial pneumonia and ARDS, where D-dimer elevation and platelet decrease are common features of the most severe forms (7, 8). However, COVID-19 displays peculiar characteristics that suggest a unique pathogenesis (9, 10). First, COVID-19-related ARDS is characterized by its protracted course, lasting twice as long as ARDS of other origins (11, 12). Second, our team and others recently identified a unique cytokine response with higher concentrations of CXCL10 and GM-CSF in COVID-19 patients that were independently associated with outcomes (13–16). In addition, COVID-19 ARDS patients experience significantly more thrombotic complications than non-COVID-19 ARDS patients, mainly pulmonary embolism and deep vein thrombosis (2). Autopsy findings have shown thrombosis and microangiopathy in the small vessels and capillaries of the lungs (17, 18) that could contribute to this poor outcome. However, comparisons between COVID-19 and pneumonia of other origins are scarce (2).

Thus, our study addressed this issue by comparing coagulopathy-related biomarkers between non-COVID-19 and COVID-19 patients with severe pneumonia and studying the association between such biomarkers and clinical outcomes.

Materials and Methods

Study Design and Patient Selection

The present work is a prospective, exploratory ancillary study of the ongoing LYMPHONIE project (registered with ClinicalTrials.gov under the number NCT03505281), initiated in November 2018 at the University Hospital of Dijon-Bourgogne (France). Patients were eligible if they had severe community-acquired pneumonia (CAP) associating: (1) pneumonia (≥2 acute signs including cough, purulent sputum, dyspnea, chest pain, temperature <35 or ≥38.5°C, and new radiological pulmonary infiltrate); and (2) severity defined by at least two criteria of the quick-Sequential Organ Failure Assessment (SOFA) score (systolic blood pressure ≤ 100 mm Hg, respiratory rate ≥22, Glasgow score <15) and/or the need for mechanical ventilation (MV) and/or vasopressors; and (3) diagnosed within 48 hours following admission. Non-inclusion criteria were: <18 years, pregnant women, persons under legal protection, decision to limit care, known immune deficiency, chronic disorder known to cause deep lymphopenia (i.e. cirrhosis, lympho- or myeloproliferative syndrome, solid cancer or active systemic lupus), hospitalization for sepsis within the previous 3 months. Non-COVID-19 CAP patients were included until February 20, 2020, one month before the COVID-19 pandemic started in Burgundy, France. COVID-19 patients were eligible if they tested positive for SARS-CoV-2 by reverse transcriptase-polymerase chain reaction (RT-PCR) on one respiratory sample. Oral consent was obtained from the patient or their legal representative. Approval was obtained from the ethics committee (Comité de Protection des Personnes SUD MEDITERRANEE V; 2017-A03404-49). Ethylenediamine tetraacetic acid blood samples were obtained after inclusion of the patient and within 48 hours of hospital admission. After centrifugation of blood samples at 2,000 × g for 10 min at 4°C, plasma was collected and stored at −80°C in the biological resource center (CRB Ferdinand Cabanne; NF S96-900 certification).

Objectives of the Study

The main objective of this study was to compare plasma biomarkers of coagulopathy between non-COVID-19 and COVID-19 with severe pneumonia.

The secondary objective was to study the association between the biomarkers of coagulopathy, and clinical outcomes potentially related to coagulopathy.

Variables of Interest, Clinical Outcomes, and Data Collection

Clinical and biological parameters, Charlson comorbidity index and severity scores (SOFA (19), Simplified Acute Physiology Score (SAPS II) (20) and Pneumonia Severity Index (PSI) (21)) were calculated at the time of inclusion. ARDS was defined according to the Berlin definition (22), and septic shock was defined as persistent hypotension requiring vasopressors and a serum lactate level >2 mmol/L despite adequate volume resuscitation. Clinical outcomes were recorded up to 30 days after admission, namely: deep vein thrombosis, pulmonary embolism, 30-day mortality, hospital- and ICU- length of stay and duration of MV. Venous thrombotic events were diagnosed as part of routine care, either with computed tomographic pulmonary angiography (pulmonary embolism) or vascular ultrasonography (deep vein thrombosis). Dedicated clinical research assistants collected all data using a standardized electronic case report form. Automatic checks were generated for missing or incoherent data.

Plasma Biomarkers of Coagulopathy Measurement

Thirteen analytes were quantified in plasma using a Luminex® Human Magnetic assay (R&D Systems, USA) according to the manufacturer's instructions: endothelial and procoagulant response [D-dimer, von Willebrand factor A2 (vWF-A2), u-plasminogen activator (uPA), tissue factor pathway inhibitor (TFPI), soluble thrombomodulin (sTM) receptor from endothelial cells participating in the control of thrombin, soluble triggering receptor expressed on myeloid cells 1 (sTREM-1), soluble intercellular adhesion molecule 1 (sICAM-1) and soluble vascular cell adhesion molecule 1 (sVCAM-1) from endothelial cells, vascular endothelial growth factor (VEGF)], and platelet activation [C-X-C motif chemokine ligand (CXCL) 4, soluble P-selectin, platelet derived growth factor (PDGF)-AA, PDGF-AB/BB]. In addition IL-1β, IL-6, IL-10, CXCL-10, GM-CSF and TNF-α, were also quantified in plasma using a Luminex® Human Magnetic assay (R&D Systems, USA) to study the potential link between immune response and coagulopathy (16). Biomarkers measurements were performed on the same day, by the same person.

Statistical Analysis

Data are described according to COVID-19 status (i.e. non-COVID-19 vs. COVID-19 patients). Continuous variables are expressed as mean ± standard deviation (SD) or median and inter-quartile range (IQR), according to their distribution, and categorical variables as number and percentage. Univariate comparisons were performed using Student's test for means, Wilcoxon Mann-Whitney test for medians and IQRs, and the Chi-square test (or Fisher's exact test when appropriate) for percentages. Coagulopathy biomarkers are presented as boxplots to visualize potential associations with COVID-19 status. The Benjamini–Hochberg procedure was used to control the false discovery rate (FDR) induced by multiple comparisons (23). For each test, a Q-value (adjusted p-value) was calculated, representing the minimum FDR at which the statistical test may be considered as a true positive (i.e. statistically significant given the multiplicity of tests). A result was deemed false positive when the corresponding Q-value was >0.05. To account for potential confounders, multivariable linear regressions considering the relevant coagulopathy biomarkers as an independent variables and COVID-19 status as the main explicative covariate were performed. These models were systematically adjusted for the most clinical relevant variables (i.e. age, sex, Charlson comorbidity index and SOFA score). Absence of serial correlation and heteroscedasticity were assessed using the DW statistic (24) and the White test (25), respectively.

Principal component analysis (PCA) was used to identify potentially significant patterns of 30 variables: the 13 coagulopathy biomarkers and the following 17 clinical and biological variables: PSI, SOFA, CRP, PCT, Creatinine, NT-ProBNP, leucocytes, neutrophils, lymphocytes, monocytes, platelets, prothrombin time (PT), PaO2:FiO2, lactate, MV duration, ICU and hospital length of stay. PCA identifies factors, called principal components, that induce the most variation in the overall data (26). These factors can be expressed as a linear combination of the correlated original variables (OVs). By inversing these formulas, we can express each OV as a linear combination of the factors and coefficients defining these linear combinations may be interpreted as correlation coefficients. Moreover, each factor describes a percent of variation in the OVs. The number of factors to retain was determined using the scree plot (27) and the clinical interpretability of factors (28). Besides, patient OV data can be projected on the planes defined by the factors retained, making it possible to observe patient patterns in a two-dimensional plot.

To strengthen interpretation, Spearman correlations were computed between coagulopathy biomarkers and the most pertinent patient characteristics associated with COVID-19 status in univariate analyses and PCA patterns. Then, multivariable linear regressions were applied to estimate the associations between the selected biomarkers of coagulopathy and MV duration, considered as the independent variable. These models were systematically adjusted for the most pertinent variables (i.e. age, pre-hospital anticoagulation therapy, COVID-19 status, and SOFA score). As above, absence of serial correlation and heteroscedasticity were assessed using the DW statistic (24) and the White test (25), respectively. The interaction between COVID-19 status and the selected biomarkers of coagulopathy was systematically tested. The proportion of variance explained by the models was quantified using the R² coefficient. Measures of association are expressed as mean differences ± standard error (SE). Finally, a heatmap correlation matrix were built between biomarkers of coagulopathy and plasma cytokines concentrations. A p-value < 0.05 was considered statistically significant. Analyses were performed using SAS version 9.4 (SAS Institute Inc., Cary, NC, USA).

Results

Characteristics of the Study Population

Sixty-three patients were enrolled (36 in the non-COVID-19 group, and 27 in the COVID-19 group). Bacterial, viral or mixed etiologies were proven in 10 (28%), 10 (28%) and 3 (8%) patients from the non-COVID-19 group, respectively (Supplementary Table 1). Mean age was marginally lower in the COVID-19 group as compared to the non-COVID-19 group (62.5 ± 10.9 vs. 68.0 ± 13.0; p = 0.07). Other demographic and comorbidity data were not statistically different between the two groups, with the exception of the prehospital anticoagulation therapy, which was more frequent in non-COVID-19 patients (p = 0.033; Table 1). Although the PaO₂:FiO₂ ratio and SOFA score were not different between groups (p = 0.35 and p = 0.52, respectively), COVID-19 patients less often had septic shock (0 vs. 11 (31%); p = 0.0015). In COVID-19 patients, we observed significantly lower arterial lactates (p = 0.01), serum creatinine (p = 0.02), NT-proBNP (p = 0.0003), procalcitonin (p = 0.007) and CRP (p = 0.004) levels (Table 1). White blood cell and platelet counts did not differ significantly between groups, except for monocytes, which were lower (p = 0.05) in COVID-19 patients (Table 1). All patients were treated with antibiotics.

TABLE 1

Table 1. Baseline characteristics and outcomes of the study population (LYMPHONIE study, 2018–2020).

At Similar Baseline Respiratory Severity, COVID-19 Patients Were More Likely to Have Venous Thrombotic Complications and Had a Significantly Longer Duration of MV

Venous thrombotic events (i.e. deep vein thrombosis or pulmonary embolism) were more frequently diagnosed in COVID-19 patients than in non-COVID-19 patients [7 (26%) vs. 2 (6%); p = 0.0312)], with a median time from admission to diagnosis of 8 (IQR = 6–11) and 11 (IQR = 9–14) days respectively (p = 0.86). COVID-19 patients had a significantly longer duration of MV (15 days, IQR = 7–22 vs. 4 days, IQR = 0–14.5; p = 0.0049), as well as longer ICU stay (p = 0.0274) (Table 1). The 30-day mortality rate was 6% (n = 2) in the non-COVID-19 group and 4% (n = 1) in the COVID-19 group (p = 1.00).

COVID-19 Patients Displayed a Unique Plasma Signature of Coagulopathy-Associated Biomarkers

COVID-19 patients had higher plasma sVCAM1 (Figure 1). However, in COVID-19 patients, we observed lower plasma levels of D-dimers, vWF-A2, sICAM1, sTREM1, VEGF, and sP-selectin as compared to non-COVID-19 patients, and a non-significant difference for plasma levels of sTM, TFPI, uPA, CXCL4, PDGF-AA and PDGF-AB/BB (Figure 1, Supplementary Table 2).

FIGURE 1

Figure 1. Box plot showing plasma concentrations of biomarkers of coagulopathy in non-COVID-19 and COVID-19 patients (Lymphonie study 2018–2020). Plasma concentration of biomarkers of coagulation and/or endothelium and platelet activation was measured within 48 hours of hospitalization in 36 non-COVID-19 and 27 COVID-19 patients with severe pneumonia. COVID-19 patients had significantly higher plasma concentrations of sVCAM1 (D), lower plasma concentration of D-dimers, vWF-A2, sICAM1, sTREM1, VEGF, sP-selectin (A,B,E,F,I,K) and a non-significant difference for plasma levels of sTM, TFPI, uPA, CXCL4, PDGF-AA and PDGF-AB/BB. For each cytokine, P-values from Student t test and Q-value of the false discovery rate (FDR) are indicated. COVID-19, coronavirus disease 2019; CXCL4, platelet C-X-C motif chemokine ligand 4; FDR, false discovery rate; PDGF, platelet derived growth factor; sICAM-1, soluble intercellular adhesion molecule 1; sTREM-1, soluble triggering receptor expressed on myeloid cells 1; sVCAM-1, soluble vascular cell adhesion molecule 1; sTM, soluble thrombomodulin; TFPI, tissue factor pathway inhibitor; uPA, u-plasminogen activator; VEGF, vascular endothelial growth factor; vWF-A2, von Willebrand factor A2 (LYMPHONIE study, 2018–2020).

PCA identified Two Patterns Linking Coagulopathy-Associated Biomarkers and Outcomes, With a Clear Distinction Between Non-COVID-19 and COVID-19 Patients

In PCA, three factors were retained, which together accounted for 51.6% of all the information contained in the 30 correlated original variables (Supplementary Table 3, Supplementary Figure 1). Two of the three factors in the final pattern were associated with outcomes. Factor 1 associated (1) baseline severity and extra-respiratory organ dysfunction (SOFA and PSI scores, lactate, creatinine and NT-ProBNP levels) and (2) “classical” markers of coagulation and endothelium activation (i.e. D-dimer, vWF-A2, sTM, sICAM1, TREM1, TFPI, uPA, VEGF). Factor 3 associated: (1) outcomes (duration of MV, ICU and hospital length of stay) and (2) sVCAM1 (Supplementary Table 3). The projection of patient data onto the directions defined by Factors 1 and 3 showed a clear separation between COVID-19 and non-COVID-19 patients. By plotting patients according to different origins of pneumonia, we observed that most of the non-COVID-19 patients with proven bacterial (or mixed) origin had high Factor 1 values, while COVID-19 patients had high Factor 3 and low Factor 1 values (Figure 2). Factor 2 associated: (1) white blood cell and platelets levels and (2) markers of platelet activation and growth factors (ie. CXCL4, sP-selectin, PDGF-AA, PDGF-AB/BB) (Table 2). Planes defined by factor 2 did not enable any discrimination between patients according to COVID-19 status.

FIGURE 2

Figure 2. Two-dimensional score plot of principal component analysis according to pneumonia etiology (Lymphonie study 2018–2020). Principal component analysis (PCA) was used to identify potentially significant patterns of 30 variables: 13 plasma biomarkers of coagulopathy and 17 clinical and biological variables from 63 patients with severe pneumonia [non-COVID-19 (n = 36), COVID-19 (n = 27)]. Factors 1 and 3 were used to build a two-dimensional score plot of PCA. There is a clear separation between COVID-19 patients (red circles) and non-COVID-19 patients: bacterial (blue triangles), mixed (blue diamonds), viral (gray crosses) and other non-documented severe community acquired pneumonia (gray Xs). Bacterial and mixed severe pneumonia were represented by a high factor 1, while COVID-19 patients had a high factor 3 (LYMPHONIE study, 2018-2020).

TABLE 2

Table 2. Multivariable linear regression factors associated with duration of mechanical ventilation (variation in minutes) in the first 30 days of hospitalization among 63 patients with severe pneumonia (R² = 41.3%) (LYMPHONIE study, 2018–2020).

Plasma sVCAM1 Was Independently Associated With the Duration of Mechanical Ventilation

As the particular severity of COVID-19 patients was represented by a significantly longer duration of MV, we investigated whether coagulopathy-associated biomarkers was associated with this poor outcome. First, we observed a significant correlation between the duration of MV and sVCAM1 (p < 0.0001) and uPA (p = 0.001) (Supplementary Table 4). In multivariable linear regression, sVCAM1 was independently associated with a longer duration of MV (p = 0.027), after adjustment for age, COVID-19 status, SOFA score and prehospital anticoagulation therapy (Table 2). We estimated an excess of 2.8 ± 1.2 days of MV per increase of 3,000 ng/mL of sVCAM1 (p = 0.027).

In addition, multivariable analysis showed that COVID-19 status was independently associated with sVCAM1. In severe COVID-19 patients, we estimated a proper excess of 2,247 ± 695 ng/ml for sVCAM1 (p = 0.002), independently of age, sex, Charlson and SOFA scores (Supplementary Table 5). Autocorrelation and heteroscedasticity were non-significant in the models.

sVCAM1 Was Strongly Correlated With CXCL10, GMCSF and IL-10 Concentrations

Since sVCAM1 is upregulated on immune-mediated activation, we investigated whether sVCAM1 was correlated with the production of cytokines. We observed strong correlations between sVCAM1 and each of CXCL10, GM-CSF and IL-10 (respectively r = 0.71, p < 0.0001; r = 0.66, p < 0.0001 and r = 0.70, p < 0.0001), but not with IL-1β, IL-6 and TNF-α (Supplementary Figure 2).

Discussion

Our study comparing non-COVID-19 with COVID-19 patients with severe pneumonia identified two main results: (1) COVID-19 patients displayed a distinct coagulopathy signature, with significantly higher concentrations of sVCAM1, but lower concentrations of some classical biomarkers of coagulation, endothelial or platelet activation; (2) sVCAM1 was independently associated with the duration of MV and strongly correlated with some recently reported features of the dysregulated immune response in severe COVID-19 patients.

Patients with severe COVID-19 display particular characteristics as compared to pneumonia of other origins, notably with a higher risk of thrombotic complications in COVID-19 (2), and a longer duration of MV, as observed in our study. There is an urgent need to understand the peculiar pathogenesis of severe COVID-19, with a view to finding new therapeutic strategies likely to improve outcomes in the most severe forms (e.g. ARDS). Ackermann et al. showed that patients who died from COVID-19 had a distinctive vascular histologic pattern with severe endothelial injury associated with the presence of intracellular SARS-CoV-2 and peripheral T-cell infiltration (18). Different studies have highlighted complement activation (e.g. the C5aR1 axis (29)), platelet activation (4), immune-mediated thrombotic mechanisms (6), endothelial dysfunction, antiphospholipid antibody, and renin-angiotensin system dysregulation as important actors of COVID-19-associated coagulopathy (30). However, the pathophysiological process driving the poorer outcomes in severe COVID-19 pneumonia remains unclear, and comparisons with pneumonia of other origins are scarce.

In our study, we showed firstly that, regardless of baseline clinical severity, plasma concentrations of numerous biomarkers characterizing coagulation/endothelial activation (i.e. D-dimers, vWF-A2, sICAM1, sTREM1, VEGF) or platelet activation (sP-selectin) were elevated in both groups, but lower in COVID-19 compared to non-COVID-19 patients. Other markers (e.g. sTM, TFPI, uPA, CXCL4, PDGF-AA and PDGF-AB/BB) did not appear to differ between the two groups. Coagulopathy is a hallmark of ARDS regardless of its origin (8), and could contribute to the pathogenesis of lung injury and subsequent respiratory failure (31). Indeed, we observed that most of the biomarkers characterizing coagulation/endothelial and platelet activation were correlated with the baseline severity, in particular in proven bacterial pneumonia (higher Factor 1 in the PCA analysis), but not with the outcomes. In addition, the higher concentrations of D-dimers in non-COVID-19 patients could be explained by a greater activation of fibrinolysis, related to the higher rate of septic shock, as compared to COVID-19 patients. In contrast, even if these later frequently required vasopressors, they never met the criteria for septic shock in our study, considering the low concentrations of lactate.

However, plasma sVCAM1 was the only biomarker found to be higher in COVID-19 patients. Correlation and PCA analyses identified a distinct pattern linking sVCAM1 to the duration of MV, which represents the characteristic outcome associated with COVID-19 and a major concern for clinicians and health systems. Finally, multivariable analysis confirmed that sVCAM1 was independently associated with the duration of MV. Based on these analyses, we hypothesized that sVCAM1 could drive the dysregulated coagulopathy leading in turn to severe lung injury and prolonged need for MV in COVID-19 ARDS patients.

VCAM-1 is an integrin expressed in endothelial cells, but also on leukocytes, and is upregulated on immune-mediated activation (pro-inflammatory cytokines, including TNF-α and IFNs) (32). VCAM1 interacts with α4β1 integrins on leukocytes, and in turn, mediates their adhesion and activation to vascular endothelium and engages their migration (33). VCAM1 has been shown to regulate leukocyte trafficking and immune response against viral infections (34). Our results are consistent with vascular endothelial dysfunctions described in COVID-19 patients, and also mediated by SARS-CoV-2 adhesion on endothelial cells (35). However, in parallel, ICAM-1, a transmembrane integrin, mediates the migration of leukocytes into the inflamed tissues (36). Both integrins were found to be elevated in the serum of COVID-19 patients and related to severity of disease (37). Since circulating levels of these isoforms reflect enhanced expression of the respective adhesion molecules (38), we postulate that overexpression of VCAM1 in COVID-19 patients could participate in massive recruitment of immune cells in the lung, contributing to lung injury and prolonged MV. However, massive endothelial adhesion (high expression of VCAM1) with less extensive diapedesis (lower level of sICAM1 in COVID-19 compared to non-COVID-19 patients) could also contribute to the massive sequestration of immune cells in the endothelial compartment, consequently contributing to a higher risk of thrombosis (36).

Accordingly, anti-VCAM1 humanized monoclonal antibody (36) could represent a promising therapeutic strategy, in order to dampen coagulopathy and improve outcomes in severe forms of COVID-19. However, since VCAM1 expression is induced by inflammatory mediators (e.g. TNF-α, IL-1β, IFN-α and γ) (32, 39), targeting the immune response is also another possible way to attenuate the deleterious effect of VCAM1-mediated endothelial and lung injury. Interestingly, we observed that sVCAM-1 was highly correlated with GM-CSF, the T_H1 chemokine CXCL10 and IL-10, which our team and others recently identified as the peculiar “cytokine/chemokine signature” associated with poorer outcome in COVID-19 patients (13, 14, 16). A more robust T_H1 immune response, and subsequent myeloid cell activation in COVID-19 could explain this peculiar cytokine signature, in which VCAM1 could be the consequence, but also a protagonist of the auto-inflammatory loop mediating endothelial and lung immune cell infiltration and tissue damage. Mitigating this immune response could also be a promising therapeutic avenue, since targeting upstream mediators would impact on both immune dysregulation and coagulopathy.

Our study may suffer from a lack of power given the large number of variables studied. However, it was an exploratory study in the context of a pandemic emergency and the multiplicity of tests was accounted for with the appropriate statistical techniques. In addition, we used several statistical methods to assess the association between COVID-19 status and biomarkers, for which data are scarce in the available literature. Second, we were unable to measure VCAM-1 expression in the tissues, even though no study to date dealing with COVID-19 was able to assess such data. This was a single center study with a small sample size and a short follow-up. A control group with non-severe COVID-19, and a follow-up measurement of clinical and biological data would have been more accurate to better assess the prognostic value of plasma biomarkers of coagulopathy in COVID-19. Since every patient did not receive systematically CT-pulmonary angiography and ultrasonography, venous thrombotic events may have been underestimated. Finally, we did not measure all the paths involved in coagulopathy (e.g. complement, or other markers of platelet activation), although we did observe lower levels of sP-selectin in COVID-19 patients, suggesting lower levels of platelet activation (4, 29). In addition, other biomarkers of endotheliopathy, such as angiopoietin 2, would have been worth evaluating in our study since it was found as a predictive biomarker of ICU admission in COVID-19 patients (40).

In conclusion, COVID-19 patients displayed a distinct coagulopathy signature with significantly higher concentrations of sVCAM1, but lower concentrations of some classical biomarkers of coagulation or platelet activation, compared to non-COVID-19 patients. This could result into a higher level of immune-mediated endotheliopathy. sVCAM1 was independently associated with the duration of MV. Further research is needed to investigate the role of VCAM-1 (or the immune response driving its expression), in the pathogenesis of severe COVID-19.

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 author/s.

Ethics Statement

The studies involving human participants were reviewed and approved by (Comité de Protection des Personnes SUD MEDITERRANEE V; 2017-A03404-49). Oral consent was obtained from the patient or their legal representatives.

Author Contributions

MB, AB, CB, and LP: Concept and design. MB, J-PQ, MN, BB, and P-EC: Recruitment of patients. MB, EM, AB, CB, and LP: Acquisition, analysis, or interpretation of data. MB, EM, and AB: Drafting of the manuscript. MB, EM, AB, J-PQ, MN, BB, P-EC, CB, and LP: Critical revision. CB and LP: Supervision. All authors contributed to the article and approved the submitted version.

LYMPHONIE STUDY GROUP CONTRIBUTORS

Pascal Andreu, François Aptel, Auguste Dargent, Marie Labruyère, Audrey Large, Sébastien Prin, (Department of Intensive Care, Dijon Bourgogne University Hospital, Dijon, France); Philippe Bonniaud, Guillaume Beltramo, Marjolaine Georges (Department of Pneumology, Dijon Bourgogne University Hospital, Dijon, France); Philip Bielefeld, Hervé Devilliers, Suzanne Mouries-Martin (Department of Internal Medicine and Systemic Diseases, Dijon Bourgogne University Hospital, Dijon, France); Bernard Bonnotte, Alexandre Guilhem (Department of Internal Medicine and Clinical Immunology, Dijon Bourgogne University Hospital, Dijon, France); Pascal Chavanet, Marielle Buisson (Infectious Diseases Department, Dijon Bourgogne University Hospital, Dijon, France); Alain Putot, Jeremy Barben (Geriatrics Internal Medicine Department, Dijon Bourgogne University Hospital, Dijon, France), Serge Monier (cytometry platform, Faculty of Medicine, Dijon, France).

Funding

This work was supported by grants (to MB) from the AOIc2020 (Appel d'Offre Interne COVID-19), from the INSERM (Centre de Recherche UMR 1231), the national research agency (ANR) Investissements d'Avenir Grant (ANR-11 LABX-0021-01, Labex Lipstic), the FEDER, the Regional Council of Bourgogne Franche-Comte and by crowdfunding (https://thellie.org/covid-19).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We thank all the investigators of the Lymphonie study group, the CIC-EC1432 (Alexandra Lamotte-Felin, Lydie Rossie and Delphine Pecqueur), the Intensive Care Medicine Research team (Solenne Villot, Mathilde Audry and Samantha Lima), the Cytometry core facility, Berengere Grammatico from the laboratory of virology and Shaliha Bechoua from the CRB (Centre de Ressource Biologique). We also thank all the donors (crowdfunding, https://thellie.org/covid-19), and the patients and healthcare team. We thank Fiona Ecarnot, PhD (EA3920, University of Burgundy Franche-Comté, Besançon, France) for proofreading and editing the manuscript.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2021.675191/full#supplementary-material

Abbreviations

ARDS, acute respiratory distress syndrome; ASAT, aspartate aminotransferase; COVID-19, coronavirus disease 2019; CXCL, C-X-C motif chemokine ligand; FDR, false discovery rate; GM-CSF, Granulocyte-macrophage colony-stimulating factor; ICU, intensive care unit; IQR, interquartile range; IL, interleukin; MV, mechanical ventilation; NT-ProBNP, N-Terminal Fragment of the Prohormone Brain-Type Natriuretic Peptide; PaO₂:FiO₂, arterial pressure of oxygen / oxygen inspiratory fraction; PDGF, platelet derived growth factor; PT, prothrombin time; RT-PCR, reverse transcriptase polymerase chain reaction; SAPS II, simplified acute physiology score II; SARS-CoV-2, severe acute respiratory syndrome coronavirus 2; SOFA, Sequential Organ Failure Assessment; vWF-A2, von Willebrand factor A2; uPA, u-plasminogen activator; TFPI, tissue factor pathway inhibitor; sTM, soluble thrombomodulin; sTREM-1, soluble triggering receptor expressed on myeloid cells 1; sICAM-1, soluble intercellular adhesion molecule 1; sVCAM-1, soluble vascular cell adhesion molecule 1; VEGF, vascular endothelial growth factor.

References

1. Liao D, Zhou F, Luo L, Xu M, Wang H, Xia J, et al. Haematological characteristics and risk factors in the classification and prognosis evaluation of COVID-19: a retrospective cohort study. Lancet Haematol. (2020) 7:e671–8. doi: 10.1016/S2352-3026(20)30217-9

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Helms J, Tacquard C, Severac F, Leonard-Lorant I, Ohana M, Delabranche X, et al. High risk of thrombosis in patients with severe SARS-CoV-2 infection: a multicenter prospective cohort study. Intensive Care Med. (2020) 1:67–72 doi: 10.1007/s00134-020-06062-x

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Tang N, Li D, Wang X, Sun Z. Abnormal coagulation parameters are associated with poor prognosis in patients with novel coronavirus pneumonia. J Thromb Haemost. (2020) 18:844–7. doi: 10.1111/jth.14768

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Hottz ED, Azevedo-Quintanilha IG, Palhinha L, Teixeira L, Barreto EA, Pão CRR, et al. Platelet activation and platelet-monocyte aggregates formation trigger tissue factor expression in severe COVID-19 patients. Blood. (2020) 136:1330–41 doi: 10.1182/blood.2020007252

CrossRef Full Text | Google Scholar

5. Goshua G, Pine AB, Meizlish ML, Chang C-H, Zhang H, Bahel P, et al. Endotheliopathy in COVID-19-associated coagulopathy: evidence from a single-centre, cross-sectional study. Lancet Haematol. (2020) 7:e575–82. doi: 10.1016/S2352-3026(20)30216-7

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Nicolai L, Leunig A, Brambs S, Kaiser R, Weinberger T, Weigand M, et al. Immunothrombotic dysregulation in COVID-19 pneumonia is associated with respiratory failure and coagulopathy. Circulation. (2020) 142:142–89 doi: 10.1161/CIRCULATIONAHA.120.048488

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Brands X, Haak BW, Klarenbeek AM, Otto NA, Faber DR, Lutter R, et al. Concurrent immune suppression and hyperinflammation in patients with community-acquired pneumonia. Front Immunol. (2020) 11:796. doi: 10.3389/fimmu.2020.00796

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Glas GJ, Van Der Sluijs KF, Schultz MJ, Hofstra J-JH, Van Der Poll T, Levi M. Bronchoalveolar hemostasis in lung injury and acute respiratory distress syndrome. J Thromb Haemost. (2013) 11:17–25. doi: 10.1111/jth.12047

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Grasselli G, Zangrillo A, Zanella A, Antonelli M, Cabrini L, Castelli A, et al. Baseline characteristics and outcomes of 1591 patients infected with SARS-CoV-2 admitted to ICUs of the Lombardy region, Italy. JAMA. (2020) 323:1574–81. doi: 10.1001/jama.2020.5394

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Iba T, Levy JH, Connors JM, Warkentin TE, Thachil J, Levi M. The unique characteristics of COVID-19 coagulopathy. Crit Care. (2020) 24:360. doi: 10.1186/s13054-020-03077-0

CrossRef Full Text | Google Scholar

11. Bellani G, Laffey JG, Pham T, Fan E, Brochard L, Esteban A, et al. Epidemiology, patterns of care, and mortality for patients with acute respiratory distress syndrome in intensive care units in 50 countries. JAMA. (2016) 315:788–800. doi: 10.1001/jama.2016.0291

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Cummings MJ, Baldwin MR, Abrams D, Jacobson SD, Meyer BJ, Balough EM, et al. Epidemiology, clinical course, and outcomes of critically ill adults with COVID-19 in New York City: a prospective cohort study. Lancet. (2020) 395:1763–70. doi: 10.1016/S0140-6736(20)31189-2

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Hue S, Beldi-Ferchiou A, Bendib I, Surenaud M, Fourati S, Frapard T, et al. Uncontrolled innate and impaired adaptive immune responses in patients with COVID-19 ARDS. Am J Respir Crit Care Med. (2020) 202:202–19 doi: 10.1164/rccm.202005-1885OC

CrossRef Full Text | Google Scholar

14. Blot M, Jacquier M, Aho Glele L-S, Beltramo G, Nguyen M, Bonniaud P, et al. CXCL10 could drive longer duration of mechanical ventilation during COVID-19 ARDS. Crit Care. (2020) 24:632. doi: 10.1186/s13054-020-03328-0

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Blot M, Bourredjem A, Binquet C, Piroth L, LYMPHONIE Study Group. Is interleukin 6 the right target in COVID-19 severe pneumonia? Am J Respir Crit Care Med. (2021) 203:139–140. doi: 10.1164/rccm.202007-2924LE

CrossRef Full Text | Google Scholar

16. Blot M, Bour JB, Quenot JP, Bourredjem A, Nguyen M, Guy J et al. The dysregulated innate immune response in severe COVID-19 pneumonia that could drive poorer outcome. J Transl Med. (2020) 18:1–4. doi: 10.1186/s12967-020-02646-9

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Fox SE, Akmatbekov A, Harbert JL, Li G, Quincy Brown J, Vander Heide RS. Pulmonary and cardiac pathology in African American patients with COVID-19: an autopsy series from New Orleans. Lancet Respir Med. (2020) 8:681–6. doi: 10.1016/S2213-2600(20)30243-5

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Ackermann M, Verleden SE, Kuehnel M, Haverich A, Welte T, Laenger F, et al. Pulmonary vascular endothelialitis, thrombosis, and angiogenesis in Covid-19. N Engl J Med. (2020) 383:383–8 doi: 10.1056/NEJMoa2015432

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Vincent JL, Moreno R, Takala J, Willatts S, De Mendonça A, Bruining H, et al. The SOFA (Sepsis-related Organ Failure Assessment) score to describe organ dysfunction/failure. On behalf of the Working Group on Sepsis-Related Problems of the European Society of Intensive Care Medicine. Intensive Care Med. (1996) 22:707–10. doi: 10.1007/BF01709751

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Le Gall JR, Lemeshow S, Saulnier F. A new Simplified Acute Physiology Score (SAPS II) based on a European/North American multicenter study. JAMA. (1993) 270:2957–63. doi: 10.1001/jama.270.24.2957

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Fine MJ, Auble TE, Yealy DM, Hanusa BH, Weissfeld LA, Singer DE, et al. A prediction rule to identify low-risk patients with community-acquired pneumonia. N Engl J Med. (1997) 336:243–50. doi: 10.1056/NEJM199701233360402

PubMed Abstract | CrossRef Full Text | Google Scholar

22. ARDS Definition Task Force, Ranieri VM, Rubenfeld GD, Thompson BT, Ferguson ND, Caldwell E, et al. Acute respiratory distress syndrome: the Berlin Definition. JAMA. (2012) 307:2526–33. doi: 10.1001/jama.2012.5669

CrossRef Full Text

23. Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A practical and powerful approach to multiple testing. J R Stat Soc B. (1995) 57:289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x

CrossRef Full Text | Google Scholar

24. Durbin J, Watson GS. Testing for serial correlation in least squares regression. II. Biometrika. (1951) 38:159–78. doi: 10.1093/biomet/38.1-2.159

PubMed Abstract | CrossRef Full Text | Google Scholar

25. White H. A heteroskedasticity-consistent covariance matrix estimator and a direct test for heteroskedasticity. Econometrica. (1980) 48:817–38. doi: 10.2307/1912934

CrossRef Full Text | Google Scholar

26. Sainani KL. Introduction to principal components analysis. PM R. (2014) 6:275–8. doi: 10.1016/j.pmrj.2014.02.001

CrossRef Full Text | Google Scholar

27. Cattell RB. The scree test for the number of factors. Multivariate Behav Res. (1966) 1:245–76. doi: 10.1207/s15327906mbr0102_10

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Cattell RB, Vogelmann S. A comprehensive trial of the scree and kg criteria for determining the number of factors. Multivariate Behav Res. (1977) 12:289–325. doi: 10.1207/s15327906mbr1203_2

CrossRef Full Text | Google Scholar

29. Carvelli J, Demaria O, Vély F, Batista L, Benmansour NC, Fares J, et al. Association of COVID-19 inflammation with activation of the C5a-C5aR1 axis. Nature. (2020) 588:146–50. doi: 10.1038/s41586-020-2600-6

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Manne BK, Denorme F, Middleton EA, Portier I, Rowley JW, Stubben CJ, et al. Platelet gene expression and function in COVID-19 patients. Blood. (2020) 136:1317–29. doi: 10.1182/blood.2020007214

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Welty-Wolf KE, Carraway MS, Ortel TL, Piantadosi CA. Coagulation and inflammation in acute lung injury. Thromb Haemost. (2002) 88:17–25. doi: 10.1055/s-0037-1613147

CrossRef Full Text | Google Scholar

32. Lechleitner S, Gille J, Johnson DR, Petzelbauer P. Interferon enhances tumor necrosis factor–induced vascular cell adhesion molecule 1 (CD106) expression in human endothelial cells by an interferon-related factor 1–dependent pathway. J Exp Med. (1998) 187:2023–30. doi: 10.1084/jem.187.12.2023

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Bazan-Socha S, Bukiej A, Marcinkiewicz C, Musial J. Integrins in pulmonary inflammatory diseases. Curr Pharm Des. (2005) 11:893–901. doi: 10.2174/1381612053381710

CrossRef Full Text | Google Scholar

34. Ou R, Zhang M, Huang L, Flavell RA, Koni PA, Moskophidis D. Regulation of immune response and inflammatory reactions against viral infection by VCAM-1. J Virol. (2008) 82:2952–65. doi: 10.1128/JVI.02191-07

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Huertas A, Montani D, Savale L, Pichon J, Tu L, Parent F, et al. Endothelial cell dysfunction: a major player in SARS-CoV-2 infection (COVID-19)? Eur Respir J. (2020) 56:2001634. doi: 10.1183/13993003.01634-2020

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Cox D, Brennan M, Moran N. Integrins as therapeutic targets: lessons and opportunities. Nat Rev Drug Discov. (2010) 9:804–20. doi: 10.1038/nrd3266

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Tong M, Jiang Y, Xia D, Xiong Y, Zheng Q, Chen F, et al. Elevated expression of serum endothelial cell adhesion molecules in COVID-19 patients. J Infect Dis. (2020) 222:894–8. doi: 10.1093/infdis/jiaa349

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Spronk PE, Bootsma H, Huitema MG, Limburg PC, Kallenberg CG. Levels of soluble VCAM-1, soluble ICAM-1, and soluble E-selectin during disease exacerbations in patients with systemic lupus erythematosus (SLE); a long term prospective study. Clin Exp Immunol. (1994) 97:439–44. doi: 10.1111/j.1365-2249.1994.tb06107.x

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Nagata M, Sedgwick JB, Kita H, Busse WW. Granulocyte macrophage colony-stimulating factor augments ICAM-1 and VCAM-1 activation of eosinophil function. Am J Respir Cell Mol Biol. (1998) 19:158–66. doi: 10.1165/ajrcmb.19.1.3001

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Smadja DM, Guerin CL, Chocron R, Yatim N, Boussier J, Gendron N, et al. Angiopoietin-2 as a marker of endothelial activation is a good predictor factor for intensive care unit admission of COVID-19 patients. Angiogenesis. (2020) 23:611–20. doi: 10.1007/s10456-020-09730-0

PubMed Abstract | CrossRef Full Text | Google Scholar