Endothelial dysfunction determines severe COVID-19 in combination with dysregulated lymphocyte responses and cytokine networks

Dynamic changes in the entire immune cell composition in severe COVID-19 patients with memory T cell development and plasmablast expansion

A total of N=42 patients with n=86 samples from MHH with confirmed SARS-CoV-2 infection were analyzed at various time points after symptom onset. N=25 patients (n=58 samples) have been admitted to ICU and N=17 patients (n=28 samples) were convalescent former ICU patients (CONV). N=29 unexposed individuals (UE, n=36) served as control group (Suppl. Tab. 1). To analyze the leukocyte dynamics during COVID-19, we quantified their absolute numbers and proportions in ICU, CONV and UE. Unsupervised cluster analyses of our immune phenotyping data sets revealed a clear separation of ICU from CONV and UE, indicating significant dynamic changes in the lymphocyte compartment during acute COVID-19 (Fig. 1 A-B, S1A). As known from other studies (14,15), both absolute numbers as well as frequencies of granulocytes were highly elevated in blood of ICU patients while monocytes were stable and only proportionally reduced (Fig. 1C, S1B). In contrast, absolute numbers as well as frequencies of lymphocytes, especially T cells, were significantly decreased in ICU compared to CONV and UE samples, with some patients showing strong lymphopenia (<1000 lymphocytes/µl blood). The number of NK cells was also decreased, whereas B cell numbers were unchanged in ICU compared to CONV and UE. Leukocyte numbers and frequencies of CONV displayed no significant differences compared to UE (Fig. C, S1B). Interestingly, four out of five deceased ICU patients clustered close to each other and exhibited decreased numbers of T cells and an activated immune phenotype (Fig. 1B, S1A). As reported earlier, ICU patients show a strong decline in lymphocyte numbers, including T and NK cells, leading to lymphopenia along with an unusual expansion of granulocytes. Furthermore, comparable leukocyte numbers and frequencies of CONV and UE indicate recovery at the cellular level in the periphery of CONV.

• Download figure
• Open in new tab

Fig. 1: Dynamic changes in the entire immune cell composition in severe COVID-19 patients with memory T cell development and plasmablast expansion.

Immune cell distribution reflected by absolute numbers in patient blood was analyzed using TruCount analyses. (A) Principle component analysis. Multigroup comparison with a p-value cut off of 0.041 was used to identify significant differences among UE (n=36), ICU (n=58) and CONV (n=28). (B) Heatmap analysis. Multigroup comparison with a p-value cut off of 0.041 was used to identify significant differences among UE (n=36), ICU (n=58) and CONV (n=28). (D-F) Numbers of different immune cells in patient blood. T cells: naïve (CCR7+CD45RO-), central memory (CM, CCR7+CD45RO+), effector memory (EM, CCR7-CD45RO+) and TEMRA (CCR7-CD45RO-); B cells: naïve (IgD+CD27-) memory (mem, CD27+IgD-), switch precursor (switch pre, CD27+IgD+), effector memory (eff mem, IgD-CD27-) and plasmablasts (CD19+CD20-CD27+CD38+) Black triangles represent last samples from deceased patients.

UE: Unexposed donors, ICU: Intensive care unit patients, CONV: convalescent patients, CM: central memory, EM: effector memory, mem eff: memory effector

Statistical analysis: ANOVA test with Turkey multiple comparison test or Kruskal-Wallis with test with Dunn’s multiple comparison test were performed. *p < 0.05, **p<0.01, ***p<0.001, ****p<0.0001.

To further investigate the T cell phenotypes in severe COVID-19, the distribution of naïve vs. memory CD4⁺ and CD8⁺ T cells was determined for ICU, CONV and UE based on expression of the chemokine receptor CCR7 and the CD45RO memory marker (Fig. S2A). Especially ICU showed highly variable numbers and proportions of all T cell subsets. CD8⁺ T cell numbers, particularly CCR7⁺CD45RO^- naïve and CCR7⁺CD45RO⁺ central memory (CM) T cells were decreased in ICU compared to CONV and UE (Fig. 1D). Also, all CD4⁺ T cells subsets were significantly reduced in numbers (Fig. 1E). With regard to CONV, we found elevated numbers of CD8⁺ T cells compared to UE, especially CCR7^-CD45RO⁺ effector memory (EM) and CCR7^-CD45RO^- TEMRA at the cost of naïve CD8⁺ T cells, which translated into a reduced CD4⁺/CD8⁺ ratio compared to UE and ICU (Fig. 1D). Thus, CD4⁺ and CD8⁺ T cells were both affected by the lymphopenia in COVID-19 ICU patients and particularly elevated CD8⁺ TEMRA T cell numbers in CONV suggest a pronounced T cell differentiation during SARS-CoV-2 infection.

In contrast to T cells, B cells were less affected by the lymphopenia in ICU (Fig. 1C, F). Detailed analyses of B cell subsets showed a broad range of IgD⁺CD27^- naïve, CD27⁺IgD^- memory (mem), CD27⁺IgD⁺ switch precursor (switch-pre) and IgD^-CD27^- effector memory (eff-mem) B cells in ICU accompanied by an unusual expansion of CD19⁺CD20^-CD27⁺CD38⁺ plasmablasts (Fig. 1F, S1E, S2C). Even though we did not observe significant differences in absolute numbers of these B cell subsets in ICU, frequencies of naïve B cells were increased and frequencies of eff-mem B cells were decreased in CONV compared to UE and ICU (Fig. 1F, S1E). Concluding, elevated absolute numbers and proportions of plasmablasts in ICU patients indicate a strong systemic B cell response during acute SARS-CoV-2 infection followed by a reverse transformation of the B cell population in CONV with elevated proportions of naïve B cells and reduced proportions of eff-mem B cells.

Plasmablasts contribute to antibody development and blocking activity is associated with disease severity

Based on the massive plasmablast expansion in ICU, we hypothesized that their presence in circulation might be accompanied by the appearance of SARS-CoV-2-specific antibodies. A Luminex-based multiplex assay was used for the detection of IgM, IgA and IgG specific for the S1, the receptor binding domain (RBD), the S2 regions of the spike protein (S) and the nucleocapsid (N) antigen of SARS-CoV-2 in plasma of ICU and CONV, while UE served as controls. In ICU, high IgG relative concentrations (MFI) were detected against all four SARS-CoV-2 antigens with the exception of two B cell-depleted patients, while IgM and IgA displayed broad variations of 45-100% positive samples (Fig. 2A-C). In contrast, CONV showed decreased proportions of seroconversion and lower concentrations of IgM and IgA (Fig. 2A-D). Since the ICU samples were obtained during acute infection, they were collected earlier in the course of disease than CONV samples, which was reflected by the different IgM, IgA and IgG kinetics and thus illustrated a time-dependent antibody class-switch from IgM to IgG and IgA (Fig. S3A, B). The decline of IgM and IgA levels 20 days after symptom onset (DASO) in ICU was accompanied by increased IgG levels that remained high up to 200 days after symptom onset in CONV (Fig. S3A, B). No major differences in IgM, IgA and IgG responses were detected between the S-and N-antigens arguing for humoral responses against multiple SARS-CoV-2 antigens (Fig. 2D).

• Download figure
• Open in new tab

Fig. 2: Plasmablasts contribute to antibody development and blocking activity is associated with disease severity

(A-E) Luminex-based multiplex assay was used to detect IgM, IgA and IgG antibodies against S1-, RBD, S2-, or N-antigen of SARS-CoV-2 in patient sera. (A-C) Percentage of seroconverted ICU (n=58) and CONV (n=28) patients for SARS-CoV-2-specific IgM (A), IgA (B) and IgG (C) antibodies. Threshold for positive samples was calculated based on the mean fluorescent intensity (MFI) of antibodies from 36 UEs + 2x standard deviation. (D) Antibody levels from UE (n=36), ICU (n=58) and CONV (n=28) are displayed as MFI. (E) Correlation analysis between S-and N-specific antibodies and plasmablasts proportions. (F) Percentage of ICU (n=58) and CONV (n=28) which developed blocking antibodies against SARS-CoV-2 RBD (left) and efficacy of blocking activity (right), which was assessed by competitive ELISA. Efficient blocking was expressed as the percent blocking at a 1:50 serum dilution relative to a UE serum control. (G) Correlation analysis between COVID-19 severity markers (CRP levels, SOFA-, WHO-score, PF ratio) and blocking efficiency.

Black triangles represent last samples from deceased patients.

Statistical analysis: Multigroup comparisons were performed using ANOVA test with Turkey multiple comparison test or Kruskal-Wallis with test with Dunn’s multiple comparison test. Two group comparison was performed using Mann-Whitney test; Spearman-Correlation. *p < 0.05,

**p<0.01, ***p<0.001, ****p<0.0001.

To unravel a potential relationship between plasmablast expansion and the occurrence of S-or N-specific antibodies, we correlated the proportions of plasmablasts with the corresponding relative antibody levels. While no correlation was found between plasmablasts and IgM antibodies, there was a weak positive correlation between plasmablasts and S-and N-specific IgA and IgG (Fig. 2E). Moreover, a competitive ELISA was performed to analyze whether antibodies in ICU were able to block in vitro binding of SARS-CoV-2 RBD to the ACE2 receptor. Interference was calculated as percent blocking in the presence of ICU plasma (1:50 dilution) relative to UE control plasma without SARS-CoV-2 antibodies. Interestingly, 85% of ICU and 75% CONV plasma samples contained interfering antibodies (Fig. 2F), but the blocking efficacy of ICU antibodies was significantly higher (median 40%) than CONV antibodies (median 24%) (Fig. 2F). Especially deceased patients exhibited high blocking activity (>60%) except for one B cell-depleted patient. Moreover, we correlated antibody blocking efficacy with disease severity makers: C-reactive protein (CRP) level as a marker for systemic inflammation, sepsis-related organ failure assessment (SOFA) score for the extent of systemic organ failure, World Health Organizations (WHO) eight-point scale for COVID-19 trial endpoints (https://www.who.int/publications/i/item/covid-19-therapeutic-trial-synopsis) and PaO2/FiO2 (PF) ratio representing pulmonary function. CRP levels, SOFA-and WHO-scores showed weak positive correlations with the antibody-blocking efficacy in ICU, whereas PF ratios were negatively correlated. Interestingly, deceased patients displayed remarkably high CRP-levels, SOFA-and WHO-scores and low PF ratios (Fig. 2G).

These results illustrate that the majority of ICU developed rapid and sustained IgG antibody responses specific for SARS-CoV-2 S-and N-proteins and their blocking efficacy seemed to increase with disease severity suggesting that antibodies were insufficient to protect from developing severe COVID-19.

Expansion of activated and terminally differentiated T and NK cells in COVID-19 ICU patients

Next, we aimed to identify the major immune cell populations significantly altered in ICU versus UE to determine ICU-associated cellular patterns. Unbiased volcano plot analysis depicted a strong decline in lymphocytes but high frequencies of HLA-DR⁺CD38⁺ activated CD4⁺ and CD8⁺ T cells together with HLA-DR⁺ NK cells and CD57⁺ terminally differentiated T and NK cells in ICU compared to UE (Fig. 3A, B). These high frequencies of activated HLA-DR⁺CD38⁺ T cells (Fig. 3C) together with HLA-DR⁺ and CD69⁺ activated CD56^dim NK cells (Fig. 3D) were associated with acute COVID-19 and decreased with recovery in CONV (Fig. 3E, F). Proportions of CD28^-CD27^- or CD57⁺CCR7^- terminally differentiated T and CD57⁺ NK cells were increased in ICU and remained elevated in CONV compared to UE (Fig. 3F). On the other hand, CD127⁺CD8⁺ T cells frequencies were decreased in ICU as well as in CONV compared to UE. Thus, severe SARS-CoV-2 infection is associated with strong activation and differentiation into late memory T and NK cells.

• Download figure
• Open in new tab

Fig. 3: Expansion of activated and terminally differentiated T and NK cells in COVID-19 ICU patients

(A) Volcano Plot visualizing two-group comparison of flow cytometry data including proportions of immune cells from ICU (n=58) and UE (n=36).

(B) Significantly (p ≥ 0.05) altered immune cell subsets between ICU (n=58) and UE (n=36) ordered based on their fold-change. (C, D) Representative flow cytometry plots (E, F) Representative flow cytometry data of immune cell frequencies from UE (n=36), ICU (n=58) and CONV (n=28).

Statistical analysis: Multigroup comparison was performed using ANOVA test with Turkey multiple comparison test or Kruskal-Wallis with test with Dunn’s multiple comparison test;

*p < 0.05, **p<0.01, ***p<0.001, ****p<0.0001.

In severe COVID-19 lymphocyte dynamics correlate with disease progression

In order to examine the implications of the significant changes within the lymphocyte compartment in ICU, they were correlated to disease severity (SOFA-, WHO-score, CRP level) and progression (days after symptom onset DASO). Several immune cell subsets significantly correlated with DASO, for example plasmablasts with a negative correlation, showing that plasmablasts developed early during severe disease but decreased with prolonged ICU hospitalization (Fig. 4A, B). In the myeloid compartment, proportions of CD14⁺CD16⁺ intermediate and CD14^-CD16⁺ non-classical monocytes showed a positive correlation to DASO. In parallel, NK cells increased with a simultaneous reduction of activated CD69⁺ or CD25⁺ NK cells over time (Fig. 4C, D). Moreover, different T cell subsets correlated with disease duration, indicating a time-dependent T cell differentiation caused by severe COVID-19 (Fig. 4E). While naïve CD8⁺, CD4⁺ T and CD25^hiCD127^lo Treg cells decreased after symptom onset, EM as well as HLA-DR⁺CD38⁺ activated CD8⁺ T cells accumulated over time (Fig. 4E, F). As expected, frequencies of naïve CD8⁺ T cells negatively and EM CD8⁺ T cells positively correlated with patient age in ICU (Fig. S4A, B). However, since no correlation was observed between COVID-19 disease duration and patient age we conclude that patient age had a minor impact on T cell differentiation during severe COVID-19 (Fig. S4C). Remarkably, CD4⁺ T cells and CD8⁺ TEMRA T cells did not follow the usual age correlation in ICU compared to UE (Fig. S4 A, B).

• Download figure
• Open in new tab

Fig. 4: In severe COVID-19, lymphocyte dynamics correlate with disease progression

(A, C, E) Waterfall-plot representing correlation coefficient (r) of Spearman-correlation analysis between immune cell proportions of B cells, monocytes, lymphocytes and granulocytes (A) or NK cells (C) or T cells (E) and disease duration as days after symptom onset (DASO) from ICUs (n=58). Red columns represent significant results. Dotted lines represent the minimal Spearman correlation coefficient (r) required for significant correlation (B, D, F) Spearman correlation analysis between proportions of plasmablasts (B), NK cells (D) or different T cell populations (F) and DASO from ICUs (n=58). (G) tSNE analysis of flow cytometry data from ICU (n=58). Variance-value cut off of 0.305 was used to identify significant differences within ICU cohort. Patients were classified into two groups. Blue: Late subgroup (n=32), yellow: early subgroup (n=26). (H) Heatmap of flow cytometry data including 98 immune cell populations from ICU (n=58). Variance-value cut off of 0.035 was used to identify significant differences within ICU cohort. Samples and immune cell subsets were ordered according to hierarchical clustering. Blue to yellow scale indicates the prevalence of each subset. Missing values are displayed in white. (I-J) Representative immune cell proportions from early (n=26) and late patients (n=32) within ICU cohort.

Statistical analysis: Unpaired t-test or Mann-Whitney test; Spearman-Correlation. *p < 0.05, **p<0.01, ***p<0.001, ****p<0.0001.

Unlike disease progression, severity markers (SOFA-, WHO-score, CRP level) correlated with only a few immune cell subsets, suggesting that disease severity had only a minor effect on the immune cell distribution (Fig. S5A). As expected, especially activated immune cells like CD69⁺CD8⁺ T cells and CD69⁺CD56^bright NK cells positively correlated with CRP levels.

Furthermore, t-distributed stochastic neighbor embedding (tSNE) (Fig. 4G) and unsupervised cluster analyses (Fig. 4H) revealed two distinct subgroups of ICU samples, which separated based on their DASO into a “late” and “early” subgroup (Fig. 4H, S5B, S5C). The “late” ICU subgroup consisted of samples obtained between day 7 to 68 (mean 34 days), whereas the “early” ICU subgroup consisted of samples from day 4 to 29 after symptom onset (mean 15 days). Neither of the two groups displayed significant differences regarding their SOFA-and WHO-scores nor PF ratios but the late subgroup exhibited decreased CRP levels compared to the early subgroup (Fig. S5D). However, detailed analyses of the immune cell compositions illustrated that these two subgroups could additionally be distinguished by several immune cell proportions (Fig. 4H-J). Particularly expression of the activation markers HLA-DR, CD38 and CD25 and frequencies of CD45RO⁺ or CD57⁺ memory T and NK cells increased with disease duration in ICU, whereas consequently, proportions of naïve T and NK cells declined over time (Fig. 4H-J).

These results show that the alterations within the lymphocyte compartment are strongly time-dependent during the individual ICU period but are less directly associated with disease severity.

Systemic inflammation and endothelial injury in COVID-19 patients

With regards to the systemic microenvironment, we quantified concentrations of 83 plasma proteins in ICU, CONV and UE including cytokines, chemokines, growth factors and endothelial factors. Samples from ICU clearly clustered separately from CONV and UE samples and displayed significantly increased levels of 65 plasma proteins (Fig. 5A). Deceased ICU patients clustered in separated patterns, indicating different plasma protein signatures associated with fatal outcome. To identify the most statistically significant plasma proteins discriminating ICU from UE, we analyzed their levels based on the fold change between ICU and UE (Fig. 5B). The volcano plot illustrates a COVID-19-associated plasma protein profile composed of pro-inflammatory cytokines and chemokines (e.g. IL-6, CXCL10) in addition to growth factors (e.g. granulocyte-macrophage colony-stimulating factor (GM-CSF)) and endothelial factors (e.g. osteopontin (OPN), angiopoietin-2 (Ang-2)). In particular, chemokines and endothelial activation markers correlated with each other in clusters (Fig. 5C), indicating a co-regulated secretion and suggesting that inflammation and endothelial injury may collectively contribute to severe COVID-19. Furthermore, we observed a high variability in the respective plasma proteins levels within the ICU cohort, as previously observed for immune cell populations (Fig. 5D, S6). In addition to classical inflammation-associated cytokines (e.g. IL-1β, IL-6, TNF-α, IFN-γ) ICU displayed increased levels of soluble interleukin-6 receptor α (sIL-6Rα), urokinase-type-plasminogen-activator (uPA) and its antagonist plasminogen activator inhibitor 1 (PAI-1) (Fig. 5D, S6). Whereas Ang-2, a central mediator of endothelial activation, was elevated, the concentrations of its soluble angiopoietin receptor 2 (sTIE-2) were stable in ICU, CONV and UE, suggesting a regulation at the ligand and not the receptor level. Endothelial activation in ICU was further detectable by significantly higher IGFBP-1, OPN and VCAM levels, accompanied by several others. Organ contribution could be seen by higher hepatocyte growth factor (HGF) levels, for instance (Fig. 5D). Of note, the plasma protein profile of CONV displayed no significant alterations compared to UE.

• Download figure
• Open in new tab

Fig. 5: Systemic inflammation and endothelial injury in COVID-19 patients

Cytokine concentrations in patient sera were measured by Luminex-based multiplex assay (A) Heatmap of cytokine data including 83 plasma proteins from UE (n=36), ICU (n=58) and CONV (n=28). P-value cut off of 0.01 was used to identify significant differences among the three cohorts. Samples and immune cell subsets were ordered using to hierarchical clustering. Blue to yellow scale indicates the prevalence of each subset. Missing values are displayed in white. (B) Volcano Plot visualizing a two-group comparison of plasma protein data from ICU (n=58) and UE (n=36). (C)Correlation matrix of 83 plasma proteins from ICUs. Spearman correlation was used to calculate correlation coefficients, which were displayed in circles. Strength of correlation was depicted by circle size. Cytokines were ordered using hierarchical clustering. Red to blue scale indicates the prevalence of each subset. (D) Representative cytokine concentrations from UE (n=36), ICU (n=58) and CONV (n=28).

Statistical analysis: Kruskal-Wallis with test with Dunn’s multiple comparison test were performed. *p < 0.05, **p<0.01, ***p<0.001, ****p<0.0001.

In summary, we observed a highly inflammatory profile of plasma proteins in ICU, indicating strong immune system activation due to SARS-CoV-2 infection. Furthermore, high levels of endothelial factors suggest that endothelial injury may be a crucial contributor to severe COVID-19 progression.

Core plasma protein signature of severe COVID-19 consists of eleven inflammatory mediators and endothelial factors

Based on the extensive secretion of immunomodulating proteins and endothelial factors in ICU, we aimed to examine their implications for disease severity (SOFA-, WHO-scores, CRP levels) and disease duration (DASO). Several plasma protein levels positively correlated with disease duration such as IL-13, IL-1β, GM-CSF and the vascular endothelial growth factor (VEGF) showing an extraordinary increase beyond day 20 after symptom onset (Fig. 6A, B). In contrast, CXCL10 and IL-12(p40) concentrations significantly decreased over time. Deceased ICU patients displayed particularly low IL-1β and IL-13 levels, while CXCL10 and IL-12(p40) plasma concentrations were increased (Fig. 6B). High SOFA scores correlated with high levels of chemokines like CXCL10 and endothelial factors like HGF, insulin-like growth factor-binding protein 1 (IGFBP-1) and uPA (Fig. 6C, D). Based on the highly inflammatory profiles of these ICU patients, it was not surprising that correlation analyses with CRP levels and WHO scores showed almost identical patterns (data not shown). Comparing plasma proteins positively correlating to either CRP levels, SOFA-or WHO-scores, we detected an intersection of seven shared plasma proteins, representing the core signature associated with severe COVID-19 (Fig. 6E). This signature was composed of the pro-inflammatory cytokines IL-12(p40) and IL-6, the growth factor stem cell growth factor beta (SCGF-β), the chemokines CXCL8-10 and the endothelial factor HGF. Furthermore IGFBP-1, IL-16, M-CSF and VCAM-1 positively correlated with SOFA-and WHO-scores, representing the contribution of the endothelium to severe COVID-19. Additional plasma proteins as IL-17F, IL-25 and CCL2 contributed to the unique WHO pattern, while IL-31, CCL27, CCL7, uPA and sCD25 were included to the unique SOFA score-associated pattern and CXCL1 completed the CRP pattern.

• Download figure
• Open in new tab

Fig. 6: Core plasma protein signature of severe COVID-19 consists of eleven inflammatory mediators and endothelial factors

Cytokine concentrations in patient sera were measured by Luminex-based multiplex assay (A, C) Waterfall-plot representing correlation coefficient (r) of Spearman-correlation analysis between cytokines and disease duration as DASO (A) or SOFA score (C) from ICUs (n=58). Red columns represent significant results. Dotted lines represent the minimal Spearman correlation coefficient (r) required for significant correlation. (B) Spearman correlation analysis between representative cytokines and DASO. (D) Spearman correlation analysis of representative cytokines with SOFA score. (E) Venn-diagram displaying significantly positive correlations between plasma proteins and CRP levels, SOFA-and WHO-score. Triangles represent last samples from deceased patients.

Statistical analysis: Spearman-Correlation. *p < 0.05, **p<0.01, ***p<0.001, ****p<0.0001.

Taken together, these signatures demonstrate that severe COVID-19 is characterized by extensive secretion of multifunctional interleukins, growth factors and endothelial factors into the circulation, suggesting that SARS-CoV-2 causes more than a systemic inflammation but multiple endothelial injuries.

Endothelial dysregulation as major contributor to severe COVID-19

Finally, since we observed distinct plasma protein patterns in deceased ICU patients (Fig. 5A), we wanted to assess whether severe COVID-19 progression could be associated with distinct profiles of ICU patients. Unsupervised t-SNE analysis of plasma proteins revealed three different subgroups within the ICU cohort (Fig. 7A) that clearly clustered separately, with the highest cytokine concentrations in subgroup 3 (Sub3) followed by subgroup 2 (Sub2) and subgroup 1 (Sub1) (Fig. 7B). Furthermore, Sub1-3 differed by their disease severity (SOFA-, WHO-scores, CRP levels, PF ratio) and disease progression (DASO) (Fig. 7C, D). Of note, Sub2 contained ICU patients (n=16) with the lowest CRP levels, SOFA-and WHO-scores but highest PF ratios and with the longest disease durations (DASO) plus one deceased patient with chronic renal failure after renal transplantation (Fig. 7C, D). On the contrary, Sub3 (n=28) displayed the highest SOFA-and WHO-scores among the three subgroups containing two deceased patients. Lowest PF ratios and highest CRP-levels were observed for Sub1 (n=14) with two deceased patients. The plasma protein profiles demonstrated clear differences among these three ICU subgroups (Fig. 7E). In general, lowest cytokine levels were observed for Sub1, which nevertheless still displayed significantly higher levels of key cytokines and chemokines like IL-6, IFN-γ and CXCL10 compared to UE (Fig. S7). Sub3 exhibited the strongest inflammatory signature with high levels of pro-inflammatory cytokines like IL-6, IL-18 and TNF-α in addition to sIL-6Rα and the chemokines CXCL8-10 (Fig. 7D, S7). Remarkably, Sub3 had the highest levels of uPA, Ang-2, HGF, IGFBP-1, OPN, ICAM-1 and VCAM-1 demonstrating endothelial activation, damage and multi-organ failure (Fig. 7D, S7). Most characteristic for Sub2 was an extensive secretion of IL-1β, IL-13, GM-CSF and VEGF (Fig. 7D, S7), which was already observed as “late” ICU signature (Fig. 6B). Compared to Sub3, Sub2 exhibited significantly lower IL-6 and IL-16 levels but concentrations of IL-18, TNF-α, Ang-2 and OPN were still elevated compared to Sub1 and UE (Fig. 7D, S7). However, concentrations of tissue factors like IGFBP-1, uPA and HGF were comparable to UE, suggesting less endothelial dysfunction in Sub2 (Fig. S7).

• Download figure
• Open in new tab

Fig. 7: Endothelial dysregulation as major contributor to severe COVID-19

Plasma protein concentrations in patient sera were measured by Luminex-based multiplex assay (A) tSNE analysis of 83 plasma protein data from ICU (n=58). Variance-value cut off of 0.146 was used to identify significant differences within ICU cohort. Patients were classified into three subgroups. Green: subgroup 1 (sub1), red: subgroup 2 (Sub2), blue: subgroup 3 (Sub3). (B) Heatmap of 83 plasma protein from ICU (n=58). Variance cut off of 0.146 was used to identify significant differences among the three cohorts. Samples and cytokines were ordered using hierarchical clustering. Blue to yellow scale indicates the prevalence of each subset. Missing values are displayed in white. (C) Distribution of severity markers SOFA, WHO-score, CRP levels, PF ratio and disease duration (DASO) among the three ICU subgroups. (D) Mean-values of SOFA-, WHO-score, CRP level, PF ratio and DASO for the three ICU subgroups. (E) Representative cytokine concentrations from ICU subgroups Sub1 (n=14), Sub2 (n=16) and Sub3 (n=28).

Statistical analysis: ANOVA test with Turkey multiple comparison test or Kruskal-Wallis with test with Dunn’s multiple comparison test were performed. *p < 0.05, **p<0.01, ***p<0.001, ****p<0.0001.

In summary, the plasma protein profile of Sub3 implies massive inflammation, endothelial injury and impairment of multiple organs, which is further supported by their high SOFA-and WHO-scores. Patients from Sub2 had longer disease durations and decreased SOFA-and WHO-scores compared to patients from Sub1 and Sub3. Together with decreased pro-inflammatory markers and endothelial factors, this indicates a gradual recovery from COVID-19 for Sub2 patients. Lastly, based on their plasma protein profile and clinical data, Sub1 patients are considered as patients at the beginning of a severe COVID-19 progression with massive inflammation but restricted endothelial damage and limited impairment of multiple organs.

Study design

In total 25 patients from intensive care unit (ICU) and 17 convalescent patients were recruited to this study between April and September 2020 at Hannover Medical School (MHH). Convalescent cohort included former ICU patients. ICU and CONV had a confirmed SARS-CoV-2 infection via viral PCR. For some patients, samples were taken at multiple time points leading to a sample size of 58 for ICU patients and 28 for convalescent patients. Five patients from ICU died and were visualized with the last sample obtained before decease. Clinical severity of COVID-19 was assessed on different parameters. World Health Organizations (WHO) eight-point scale for COVID-19 trail endpoints was used to classify severity of COVID-19 (https://www.who.int/publications/i/item/covid-19-therapeutic-trial-synopsis). C-reactive protein (CRP) level was used as inflammation marker, whereas sepsis-related organ failure assessment (SOFA) score determined the extent of systemic organ failure. PF ratio represented the individual pulmonary function. Clinical and demographical characteristics of study participants are summarized in Supplementary Table 1. Of the 25 ICU patients, 12 received Remdesivir and 6 received Tocilizumab treatment. 29 unexposed individuals (UE) were enrolled as control group, some also at several time points leading to a sample size of 36 for UE.

Quantification of cells from EDTA blood via TruCount™ analysis

Absolute cell numbers were calculated from whole blood using BD Trucount™ Tubes (BD Biosciences), following manufacturer’s instructions.

Flow cytometry

For flow cytometry analyses, whole blood EDTA samples (100 µl) were incubated with antibodies for surface staining in FACS Buffer (0,1% NaN3, 2,5% FCS in PBS) at 4°C for 30 min and followed by 15 min erythrocyte lysis using 1x BD Lysing Solution. Cells were washed with PBS prior to acquisition. All antibodies used for flow cytometry analyses are listed in Supplementary Table 2. Cells were acquired and analyzed on a LSRII flow cytometer (BD Biosciences, USA) using FACS Diva software (v8.0).

Multiplex assays

Luminex-based multiplex assays were used to quantify plasma proteins and SARS-CoV-2 S-and N-specific antibodies.

Plasma proteins were measured using the Bio-Plex Pro^TM Human Assays (Bio-Rad, Hercules, USA): cytokine screening panel plus ICAM-1 and VCAM-1 (12007283, 171B6009M, 171B6022M), Cancer Biomarker Panel 1 and Panel 2 (171AC500M, 171AC600M) and Th17 cytokine Panel 15-Plex (171AA001M) following manufacturer’s instructions with thawed plasma and which was diluted twofold with assay buffer.

SARS-CoV-2 S-and N-specific antibodies were detected using the SARS-CoV-2 Antigen Panel 1 IgG, IgM, IgA assay (Millipore HC19SERM1-85K-04, HC19SERA1-85K-04, HC19SERG1-85K-04) following manufacturer’s instructions. The assays were performed with thawed plasma, which was diluted 1:200 for ICU and CONV and 1:100 for UE with assay buffer.

The semi-quantitative readout is given as median fluorescence intensity (MFI) of > 50 beads for each antigen and sample, acquired by the Bio-Plex 200 machine and the Bio-Plex Manager^TM Version 6.0 software (Bio-Rad Hercules, USA).

Competitive ELISA

Competitive ELISA was performed to detect blocking antibodies against SARS-CoV-2 RBD. Plates were coated over night at 4°C with 50 µl/well of recombinant RBD (400 ng/ml). After 3 times washing, plates were blocked for 1-2 h at room temperature (RT) with PBS+2%FCS. 0,25 µg/ml biotinylated human ACE2 (Acro Biosystems) were added followed by 50 µl patient plasma (1:50 dilution). Plates were incubated at RT for 1-2 h and washed 3 times afterwards. HRP-bound streptavidin (Merck-Millipore) was added and incubated at RT for 1-2 h, then washed 3-4 times. TMB substrate (BD Biosciences) was added and reaction was stopped with 50 µl H2SO4 (0,5 M). Optical densities were measured at 450 nm. Efficient blocking was expressed as the percent neutralization at a 1:50 plasma dilution relative to a UE control serum.

Statistical analyses

Statistical analyses of the data were performed with GraphPad Prism v7.0/v9.0 software (GraphPad Software). D’Agostino-Pearson omnibus normality test was calculated to assess data distribution. Parametric tests were performed where data were normally distributed, otherwise non-parametric tests were used. The statistical test used in each analysis is indicated in the figure legends. Correlation analyses were performed using Spearman rank-order correlation. Results were considered significant if p<0.05.

Qlucore Omics Explorer (version 3.6, Qlucore) was used to generate principal component analysis (PCA) plots, heatmaps, tSNE plots and volcano plots. For each analysis, q-value or variance used as a cut-off are indicated in the figure legends.

For the cytokine correlation matrix Spearman rank-order correlation coefficient was visualized using the corrplot R package.

Study approval

The study was approved by the Hannover Medical School Ethics Committee. All patients or participants provided written informed consent before participation in the study (9001 BO K, 968-2011).