Skip to main content

The dysregulated innate immune response in severe COVID-19 pneumonia that could drive poorer outcome

A Correction to this article was published on 08 March 2021

This article has been updated



Although immune modulation is a promising therapeutic avenue in coronavirus disease 2019 (COVID-19), the most relevant targets remain to be found. COVID-19 has peculiar characteristics and outcomes, suggesting a unique immunopathogenesis.


Thirty-six immunocompetent non-COVID-19 and 27 COVID-19 patients with severe pneumonia were prospectively enrolled in a single center, most requiring intensive care. Clinical and biological characteristics (including T cell phenotype and function and plasma concentrations of 30 cytokines) and outcomes were compared.


At similar baseline respiratory severity, COVID-19 patients required mechanical ventilation for significantly longer than non-COVID-19 patients (15 [7–22] vs. 4 (0–15) days; p = 0.0049). COVID-19 patients had lower levels of most classical inflammatory cytokines (G-CSF, CCL20, IL-1β, IL-2, IL-6, IL-8, IL-15, TNF-α, TGF-β), but higher plasma concentrations of CXCL10, GM-CSF and CCL5, compared to non-COVID-19 patients. COVID-19 patients displayed similar T-cell exhaustion to non-COVID-19 patients, but with a more unbalanced inflammatory/anti-inflammatory cytokine response (IL-6/IL-10 and TNF-α/IL-10 ratios). Principal component analysis identified two main patterns, with a clear distinction between non-COVID-19 and COVID-19 patients. Multivariate regression analysis confirmed that GM-CSF, CXCL10 and IL-10 levels were independently associated with the duration of mechanical ventilation.


We identified a unique cytokine response, with higher plasma GM-CSF and CXCL10 in COVID-19 patients that were independently associated with the longer duration of mechanical ventilation. These cytokines could represent the dysregulated immune response in severe COVID-19, as well as promising therapeutic targets. NCT03505281.


Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is responsible for pneumonia with peculiar characteristics [1]. A race against time has been launched to find effective therapies likely to improve outcome in coronavirus disease 2019 (COVID-19).

SARS-CoV-2 infects type I and type II alveolar epithelial cells as well as alveolar macrophages, through binding to angiotensin-converting enzyme 2 (ACE2), triggering a type I interferon (IFN) response, and the release of a myriad of pro-inflammatory cytokines (i.e. interferon (IFN)-γ, interleukin (IL)-1RA, IL-6, IL-8, IL-10, IL-19, monocyte chemoattractant protein (MCP)-1, MCP-2, MCP-3, C-X-C motif chemokine ligand (CXCL)9, CXCL10, CXCL5, tumor-necrosis factor (TNF)-α). Due to the massive T cell stimulation, lower levels of T lymphocytes are observed and all of these abnormalities being associated with disease severity [2,3,4,5,6]. Complement activation, and especially the C5a/C5aR1 axis was also implicated in COVID-19 lung pathology [7]. Accordingly, the description of the so-called cytokine storm has been advocated as the cause of organ dysfunction and death during COVID-19. For now, dexamethasone is the only treatment that has proven to be effective in reducing 28-day mortality in severe COVID-19 patients receiving mechanical ventilation in a randomized clinical trial [8]. These data support the existence of a unique dysregulated immune response that could be one of the most promising therapeutic targets to date.

However, several caveats must be underlined. First, immune pathogenesis is poorly understood, and comparisons of the immune response between COVID-19 patients and patients with pneumonia of other origins are scarce [4, 9]. Second, the relevance of the cytokine storm paradigm is being questioned [10, 11]. We recently reported that COVID-19 patients had lower concentrations of interleukin (IL)-6 compared to non-COVID-19 patients with severe pneumonia [12]. Others showed that mean IL-6 concentrations were nearly 100 times higher in patients with cytokine release syndrome, 27 times higher in patients with sepsis and 12 times higher in patients with ARDS unrelated to COVID-19 [13]. Third, Remy et al. showed that COVID-19 patients display a severe immunosuppressive phenotype [14]. Finally, a delayed type I IFN response is associated with an impeded viral clearance and could promote the accumulation of pathogenic inflammatory monocyte-macrophages resulting in cytokine/chemokine release within the lung [15, 16]. These results are of utmost importance, since we also recently reported that the alveolar viral load is tightly correlated with subsequent severity in COVID-19 acute respiratory distress syndrome (ARDS) [17]. Modulating immunity remains a challenge, and there is a compelling need to identify the dysregulated immune response driving the outcomes observed in COVID-19 patients, in order to find the most relevant therapeutic targets.

Thus, our study aimed to compare clinical and biological characteristics, immune response and outcomes between non-COVID-19 and COVID-19 patients with severe pneumonia.


Study design and participants

The present work is a prospective, exploratory substudy of the ongoing LYMPHONIE trial ( NCT03505281), initiated in November 2018 at the University Hospital of Dijon-Bourgogne (France). Patients were eligible if they had severe community-acquired pneumonia (CAP): 1) pneumonia (≥ 2 acute signs including cough, purulent sputum, dyspnea, chest pain, temperature < 35 °C or ≥ 38.5 °C, and new radiological pulmonary infiltrate); 2) 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 h 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. COVID-19 patients were eligible if they were 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).

Variables of interest, clinical outcomes, and data collection

Clinical and biological parameters, severity scores (SOFA [18], Simplified Acute Physiology Score (SAPS II) [19] and Pneumonia Severity Index (PSI) [20]) were calculated at the time of inclusion. ARDS was defined according to the Berlin definition [21], 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: 30-day mortality, hospital- and ICU- length of stay, duration of MV and the occurrence of ventilator-acquired pneumonia (VAP). Dedicated clinical research assistants collected all data using a standardized electronic case report form. Automatic checks were generated for missing or incoherent data.

Sample collection

Ethylenediamine tetraacetic acid blood (plasma biomarker) and heparin anticoagulated blood (cell stimulation) were obtained after inclusion of the patient (within 48 h of hospital admission, with a diagnosis of severe community acquired pneumonia and according to the inclusion and non-inclusion criteria). Within 4 h following sampling, plasma was collected after centrifugation at 2000 x g for 10 min at 4 °C and stored at −80 °C until use, without freeze–thaw cycle. All samples were collected and stored in the biological resource center of Dijon University Hospital (CRB Ferdinand Cabanne;; NF S96-900 certification).

Lymphocyte phenotyping

Absolute counts for CD3 + , CD3 + CD4 + ,CD3 + CD8 + , CD3-CD19 + , CD3-CD56 + and/or CD16 + lymphocyte subsets were performed using an AQUIOS CL flow cytometer (Beckman Coulter, Hialeah, FL). The AQUIOS CL is a single platform, fully automated volumetric flow cytometry technology and uses a 488 nm solid state diode laser to measure light diffraction, fluorescence and electronic volume which estimates the relative size of the cells. Whole blood was incubated with the monoclonal antibody reagent followed by no-wash erythrocyte lysis. A ready-to-use mix of antibodies was used. AQUIOS Tetra-1 Panel CD45-FITC/CD4-RD1/CD8-ECD/CD3-PC5 reagent provides identification and enumeration of CD45 + , CD45 + Low SS, and CD3 +/CD4 + , CD3 +/CD8 + , and CD3 + lymphocyte percentages and absolute counts in peripheral whole blood. AQUIOS Tetra-2 + Panel CD45-FITC/(CD56 + CD16)-RD1/CD19-ECD/CD3-PC5 provides identification and enumeration of CD45 + , CD45 + Low SS, and CD3 + , CD3-/CD19 + and CD3-/CD56 + CD16 + lymphocyte percentages and absolute counts in peripheral whole blood. Additionally, both panels provide for CD45 + absolute count and CD45 + Low SS absolute count and percentage. The AQUIOS System Software includes the algorithms and test definitions that provide automated analysis and results for AQUIOS reagents. Normal range (5%–95% reference ranges) values of absolute counts for immune cells and lymphocyte subsets are indicated as Ref. [22, 23].

Measurement of cytokines

Thirty analytes were quantified in plasma using the Human XL Cytokine Magnetic Luminex® assay (R&D Systems, USA) according to the manufacturer’s instructions: C–C motif chemokine ligand (CCL)2, CCL3, CCL4, CCL5, CCL11, CCL19, CCL20, soluble CD40 ligand, fractalkine, CXCL1, CXCL2, CXCL10, FMS-like tyrosine kinase 3 ligand (FLT3L), granulocyte colony-stimulating factor (G-CSF), granulocyte–macrophage colony-stimulating factor (GM-CSF), granzyme B, interferon (IFN)-α, IL-1α, IL-1β, IL-1RA, IL-2, IL-6, IL-7, IL-8, IL-10, IL-15, IL-33, programmed death-ligand 1 (PDL1), transforming growth factor (TGF)-α, TNF-α, and TNF-related apoptosis inducing ligand (TRAIL). All samples were measured in the same experiment. Briefly, on the day of the assay, plasma was centrifuged again at 16,000xg for 4 min immediately prior to use. A twofold dilution with calibrator was used for all samples. The acquisition was performed using Bio-Plex 200 system and analyzed using Bio-Plex ManagerTM software (Bio-rad, Hercules, CA). Cytokine concentrations were automatically determined from standard curves and expressed as pg/ml. Samples with values above the ranges were tested again with a 40-fold dilution. All raw data were collected by a data-manager for further analyses in the SAS Software.

Whole blood leukocyte ex vivo stimulation (WBS)

The standardized functional immunoassay QuantiFERON Monitor® (QFM, Qiagen) was used according to the manufacturer’s instructions. Within 3 h after blood sampling, one milliliter of whole blood was incubated at 37 °C for 20 ± 1 h with a QFM LyoSphere containing anti-CD3 T-cell receptor ligand and R848 (TLR7/8 ligand), or without LyoSphere (non-stimulated blood). Plasma was harvested after centrifugation at 4000 rpm for 10 min and stored at −80 °C until use. Whole blood leukocyte production of IFN-γ (IU/ml) upon stimulation was measured using ELISA (Qiagen), and fifteen other analytes using the Human Th9/Th17/Th22 Discovery Luminex® assay (R&D Systems, USA): CD40 ligand, GM-CSF, IL-1β, IL-2, IL-4, IL-5, IL-6, IL-10, IL-12, IL-13, IL-15, IL-17A, IL-33, TNF-α, CCL20. All samples were measured the same day by the same person and using the same kit. Samples with values above the ranges were tested again with further dilution. The cytokine production after stimulation was expressed as the difference of concentrations between plasma from stimulated blood and those from non-stimulated blood.

As a reference, we used samples from 7 control patients included in the Pneumochondrie study ( NCT03955887) and who underwent QuantiFERON Monitor® assay in the same conditions. The control population consisted of outpatients without fever during the previous 15 days and who underwent bronchoalveolar lavage for a non-infectious condition [24]. Samples were conditioned and measured in the same way as for the LYMPHONIE study.

Statistical analysis

Data were described according to COVID-19 status (i.e. non-COVID-19 vs COVID-19). Continuous variables were expressed as mean ± standard deviation (SD) or median and inter-quartile range (IQR), according to their distribution, and categorical variables as frequencies and percentage. Univariate comparisons were performed using Student’s test for means, Wilcoxon Mann–Whitney test for medians and IQRs and Chi square test (or Fisher’s exact test when appropriate) for percentages. Cytokines with p < 0.05 were presented by boxplots to visualize potential associations with COVID-19 status.

Principal component analysis (PCA) was used to identify potentially significant patterns of 64 variables: clinical characteristics and outcomes (n = 6), biological findings (n = 13), plasma cytokines (n = 30), cytokine production upon ex vivo stimulation (n = 15). PCA identifies factors, called principal components, that induce the most variation in the overall data [25]. 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 are interpreted as correlation coefficients. Moreover, each factor describes a percent of variation in the OVs. The number of factor to retain was determined using the scree plot [26] and the clinical interpretability of factors [27]. Finally, patients OVs data can be projected on the plans defined by the retained factors, which allows observing patient’s patterns in a two-dimensional plot.

Spearman correlations were computed between cytokines and the most pertinent clinical outcomes associated with Covid-19 status in univariate analyses and PCA patterns. To account for potential confounders, we constructed multivariable linear regressions, with the MV duration as an outcome, for each selected cytokine. adjusted for age, COVID-19 status and either SOFA score (model 1) or PaO2:FiO2 (model 2). The interaction between COVID-19 status and cytokines was systematically tested. Absence of serial correlation and heteroscedasticity were assessed using the DW statistic [28] and the White test [29] respectively. The proportion of variance explained by the models was quantified using the R2 coefficient. Measures of association are expressed as mean differences ± standard error (SE). A p-value < 0.05 was considered statistically significant. Analyses were performed using SAS version 9.4 (SAS Institute Inc., Cary, NC, USA).


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 (Additional file 1: Table S1). 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 (Table 1). Although PaO2:FiO2 ratio and SOFA score were not different between groups (p = 0.35 and p = 0.52, respectively), fewer COVID-19 patients had septic shock (0 vs. 11 (31%); p = 0.0015), and arterial lactates (p = 0.01), serum creatinine (p = 0.02), and NT-proBNP (p = 0.0003) were all lower in COVID-19 patients (Table 1). White blood cells counts did not differ significantly between groups (Table 1).

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

COVID-19 patients had a significantly longer duration of mechanical ventilation

At baseline respiratory severity, COVID-19 patients had a significantly longer duration of MV (15 [7–22] vs. 4 [0–14.5]; p = 0.0049), and ICU stay (p = 0.0274) and higher rate of ventilator-acquired pneumonia (p = 0.001) than non-COVID-19 patients (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 cytokine response pattern

COVID-19 patients had higher plasma CXCL-10 and CCL5 and marginally higher GM-CSF (Fig. 1). However, we observed lower plasma levels of FLT3L, G-CSF, CXCL1, IL-1β, IL1-RA, IL-2, IL-6, IL-8, IL-15, CCL2, CCL4, CCL19, CCL20, TGF-α and TNF-α, and a non-significant difference for plasma levels of sCD40-Ligand, CX3CL1, Granzyme B, CXCL2, INF-α, IL-1α, IL-7, IL-10, IL-33, PD-L1 and TRAIL (Fig. 1, Additional file 1: Table S2).

Fig. 1
figure 1

Box plot showing plasma concentrations of cytokines in non-COVID-19 and COVID-19 patients. Plasma concentration of cytokines was measured within 48 h of hospitalization in 36 non-COVID-19 and 27 COVID-19 patients with severe pneumonia. For each cytokine, p-values from both Student t and Wilcoxon-Mann–Whitney U tests are indicated and the difference was considered significant if at least one was < 0.05. CCL C–C motif chemokine ligand, COVID-19 coronavirus disease 2019, IL interleukin, GM-CSF Granulocyte–macrophage colony-stimulating factor, TNF tumor necrosis factor, TGF transforming growth factor (LYMPHONIE study, 2018–2020)

COVID-19 patients displayed similar T-cell exhaustion to non-COVID-19 patients

On ex vivo stimulation with CD3 and TLR7/8, T-cells from patients with severe CAP (non-COVID-19 and COVID-19), as compared to non-infectious controls, displayed a significantly decreased production of INF-γ, TNF-α, IL-1β, IL-2, IL-4, IL-5, IL-6, IL-10, IL-12, IL-13, IL-15, IL-17A and IL-33, with no difference between COVID-19 and non-COVID-19 groups (Fig. 2a–d, Additional file 1: Table S3). In addition, the IL-6:IL-10 and TNF-α:IL-10 ratios were significantly lower in COVID-19 patients as compared to non-COVID-19 patients (p = 0.0042 and p = 0.0001, respectively), driven by lower IL-6 and TNF-α levels in COVID-19 patients with similar levels of IL-10 between groups (Fig. 2e–f).

Fig. 2
figure 2

Immune-suppression phenotype and inflammatory/anti-inflammatory balance in COVID-19 and non-COVID-19 patients. Boxplot showing cytokine production (IFN-γ (a), IL-1β (b), IL-6 (c), TNF-α (d)) of blood leukocytes on ex vivo stimulation (CD3/TLR7-8 agonists), using a standardized test (QuantiFERON Monitor®) within 48 h of hospitalization in non-COVID-19 (n = 36) and COVID-19 (n = 27) patients. As a reference, the test was performed in 7 non-infected control patients included in the Pneumochondrie study (NCT03955887) [24]. Boxplot depicting IL-6:IL10 (d) and TNF-α:IL-10 ratios in non-COVID-19 and COVID-19 patients. COVID-19 coronavirus disease 2019, IL interleukin, TNF tumor necrosis factor, WBS whole blood stimulation (LYMPHONIE study, 2018–2020)

PCA identifies two patterns linking immune response and outcomes, with a clear distinction between non-COVID-19 and COVID-19 patients

In PCA, four factors were retained for interpretation, which together accounted for 53.81% of all the information obtained by using the whole 64 available variables (clinical characteristics and outcomes (n = 6), biological findings (n = 13), plasma cytokines (n = 30), cytokine production upon ex vivo stimulation (n = 15)). Two of the four factors in the final pattern were associated with outcomes (Table 2, Additional file 1: Fig. S1). Factor 1 associated [1] baseline severity and extra-respiratory organ dysfunction (SOFA and PSI scores, lactate, creatinine and NT-ProBNP levels which were higher in non-COVID-19, except for SOFA score), [2] “classical” inflammatory mediators (i.e. CXCL1, CXC2, IL-1β, IL-2, IL-6, IL-8, TNF-α) and (3) “T-cell exhaustion (inverse correlation with lymphocyte count, and production of 15 cytokines upon ex vivo stimulation of whole blood with anti-CD3 and TLR7/8 ligands). Factor 4 associated: [1] outcomes (duration of MV, ICU and hospital length of stay, which were higher in COVID-19 patients at 30 days) and (2) only GM-CSF, CXCL10 and INF-α levels (Table 2). The projection of patient data onto the directions defined by Factors 1 and 4 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 4 and moderate Factor 1 values (Fig. 3). Planes defined by factors 2 and 3 did not enable any such discrimination between patients according to COVID-19 status.

Table 2 Principal component analysis (PCA) factor pattern correlation (LYMPHONIE study, 2018-2020)
Fig. 3
figure 3

Two-dimensional score plot of principal component analysis according to pneumonia etiology. Principal component analysis (PCA) was used to identify potentially significant patterns of 65 variables (clinical (n = 8), biological (n = 12), plasma cytokines (n = 30), cytokine production on ex vivo stimulation (n = 15)) from 63 patients with severe pneumonia (non-COVID-19 (n = 36), COVID-19 (n = 27)). Factors 1 and 4 were used to build a two-dimensional score plot of PCA and COVID-19 patients (red circles) and non-COVID-19 patients (bacterial (blue triangles), mixed (blue diamonds), viral (grey crosses) and other non-documented severe community acquired pneumonia (grey Xs)) were represented. (LYMPHONIE study, 2018-2020)

Plasma GM-CSF, CXCL10 and IL-10 were 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 immune response could explain this poor outcome. We observed a significant correlation between the duration of MV and GM-CSF (p < 0.0001), IL-10 (p < 0.0001), CXCL10 (p < 0.0001), CCL-2 (p = 0.001), CX3CL1 (p = 0.0233), and Granzyme B (p = 0.0143) (Additional file 1: Table S4). Based on all these results, we performed a multivariate linear regression to identify factors associated with the duration of MV in the first 30 days, using two models (i.e. SOFA score (model 1) or PaO2:FiO2 ratio (model 2) as the variable of adjustment to account for severity) (Table 3). Only GM-CSF was independently associated with a longer duration of MV in both models. We estimated an excess of 22.11 ± 8.36 min of MV per increase of 1 pg/mL of GM-CSF (p = 0.0105, model 1; and p < 0.0001, model 2). Interleukin-10 and CXCL10 were independently associated with a longer duration of MV only when adjusted for respiratory severity (model 2; p = 0.0359 and p = 0.049 respectively) (Table 3). No significant interaction was found between these cytokines and COVID-19 status. Autocorrelation and heteroscedasticity were non-significant in all models.

Table 3 Multivariate 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 (LYMPHONIE study, 2018-2020)


In this study, we identify a dysregulated cytokine production of GM-CSF and CXCL-10 in COVID-19 patients that was independently associated with the duration of MV which represents the distinctive poorer outcome observed in severe COVID-19 patients.

COVID-19 pneumonia is unique in comparison with pneumonia of other origins, with, in particular, sudden deterioration 7–9 days after onset of symptoms, severity of hypoxemia that contrasts with the relatively well-preserved lung mechanics, and the protracted nature of ARDS [30, 31], as observed in our study. The beneficial effect of dexamethasone in the most severe forms of COVID-19 argues for a dysregulated immune response that mediates lung injury and outcome [8]. To date, the characteristics of the immune response in COVID-19 have not been completely elucidated. The terms “cytokine storm” and “macrophage activation syndrome” have been widely adopted to explain the immunopathogenesis, since the release of myriad inflammatory mediators is correlated with disease severity [4, 32]. In our study, we first showed that despite similar respiratory severity, plasma concentrations of numerous cytokines characterizing the “cytokine storm” (i.e. IL-1β, IL-6, IL-8, IL-15, TNF-α, CCL2, CCL4, CCL19, CCL20, TGF-α) were lower in COVID-19 compared to non-COVID-19 patients. These findings contrast with the results of McElvaney et al., showing that plasma IL-6 levels were higher in COVID-19, compared to non-COVID-19 patients (4), but are in line with Sinha’s retrospective observations (10). Conversely, based on a standardized functional immune-assay, we found that patients with severe pneumonia (whether COVID-19-related or not) displayed severe alterations of T-cell functionality on ex vivo CD3 and TLR7/8 stimulation. However, we observed a more unbalanced inflammatory/anti-inflammatory cytokine response in COVID-19 patients, as reflected by the IL-6:IL-10 and TNF-α:IL-10 ratios. Our results and those of Remy et al. [14] clearly challenge the classical paradigm of “cytokine storm”-mediated inflammation and show a markedly immunosuppressive phenotype in COVID-19 patients rather than hyperinflammation.

However, we identified a dysregulated immune response that was independently associated with the peculiar longer duration of MV observed in severe COVID-19 patients. PCA analysis, including a comprehensive study of inflammatory and anti-inflammatory immune responses and outcomes, identified two interesting patterns that clearly distinguish non-COVID-19 from COVID-19 patients. COVID-19 patients presented a unique phenotype associating higher levels of GM-CSF and CXCL10 and a longer duration of MV. In addition, multivariate regression analysis confirmed that GM-CSF, CXCL10, and also IL-10 were all independently associated with the duration of MV after adjustment for potential confounders. Based on this comprehensive analysis, we hypothesize that these cytokines could represent part of the dysregulated immune response driving the prolonged need for MV in COVID-19 patients.

GM-CSF is secreted by epithelial cells from injured tissue or leukocytes, to induce survival, proliferation and/or differentiation of myeloid cells [33], playing a critical role in regulating microbial defense [34]. However, as a consequence of aberrant Th-1 cell activation and inflammatory monocytes [35], aberrant production of GM-CSF may result in excessive inflammation and tissue damage, mainly by macrophage M1 polarization and overactivation [36]. CXCL10 is a pro-inflammatory Th1-chemokine driving migration to the site of infection of Th-1 T-cells, monocytes and neutrophils that express its receptor CXCR3 [37]. Production of CXCL10 has already been shown to be increased in SARS-CoV-1 [38]. Plasma concentrations of CXCL10 were recently reported to predict disease progression in COVID-19 [39, 40]. Finally, Ichikawa et al. showed that blocking the CXCL10-CXCR3 signaling pathway in viral and non-viral ARDS preclinical models improved survival [41]. High levels of interleukin-10 were also observed in our study, whether COVID-19-related or not, and were independently associated with the duration of MV. However, no statistically significant association between the COVID-19 status and IL-10 was observed, which may be explained either by the fact that IL-10 is not the only or main driver of the length of MV in COVID-19 patients, either by a lack of statistical power, or both. As recently described, concurrent immune suppression and hyperinflammation are a hallmark of the pathogenesis in non-COVID-19 CAP, and argue against two distinct phases of host response [42]. Nevertheless, in COVID-19 patients, we observed higher IL-6:IL-10 and TNF-α:IL-10 ratios that could reflect an unbalanced overproduction of IL-10. IL-10 can be produced by most cells of innate and adaptative immune system, including lymphoid and myeloid cells, resulting in pleiotropic immunosuppressive functions (e.g. inhibition of the release of pro-inflammatory mediators, inhibitory effects on T cells and monocytes/macrophages…) [43]. As for CXCL10, it can be advocated that an initial more robust TH1 response with monocytes/macrophage over activation in COVID-19 patients, as compared to non-COVID-19 patients, could lead to a subsequent IL-10 overproduction in order to limit T cell responses. The role of IL-10, either beneficial or deleterious remains a difficult issue. However, IL-10-mediated immune suppression could drive the onset of secondary infectious complications and morbi-mortality, especially since 59% of COVID-19 patients presented at least 1 VAP event.

Our results are in line with those of Hue et al. that recently identified a plasma chemokine signature in COVID-19 ARDS patients (CXCL10, GMCSF and IL-10) which was associated with mortality (9). In addition, we also recently reported that both plasma and alveolar CXCL10 concentrations were independently associated with the duration of mechanical ventilation in COVID-19 ARDS patients [24].

Corticosteroids were shown to improve survival in severe COVID-19 patients in a recent therapeutic trial [8], and this may be linked to a decrease in CXCL10 levels [40], via the inhibition of the Th-1 pathway. Since corticosteroids have many side effects, targeted therapies likely to dampen the dysregulated immune response in COVID-19 are urgently needed. Specifically, GM-CSF blockade (e.g. with lenzilumab) is increasingly being considered as a promising therapy in COVID-19 [33, 36] and is under investigation in a phase III clinical trial (NCT04351152). In addition, CXCL10 blockade (e.g. Eldelumab/MDX-1100) may also represent an attractive therapy likely to dampen the dysregulated immune response that could be driving the duration of MV.

Dampening inflammation in a context of high immune suppression is not always a hazardous route. During chimeric antigen receptor T (CAR-T) cell therapy, GM-CSF inhibition reduces cytokine release syndrome and neuro-inflammation, but enhances antitumoral CAR-T cell function [44]. In addition, an IL-6 blocker could partially rescue immune dysregulation caused by SARS-CoV-2 (2). These considerations are of utmost importance, since we reported that COVID-19 ARDS patients had a persistent alveolar SARS-CoV-2 viral load that correlates with severity [17]. Combined therapies associating immunomodulatory and antiviral agents are the most promising strategy likely to improve outcome in COVID-19 patients.

This study has several limitations. The statistical analysis suffers from a lack of power given the large number of variables studied and the small sample size. Then, we only used one approach (cytokine production after stimulation with anti-CD3 and TLR7/8 ligand) to measure T-cells exhaustion phenotype. However, it would have been also important to measure several phenotypic markers of exhaustion, namely PD-1, LAG-3, TIM-3 and CTLA-4. However, it was an exploratory study in the context of a pandemic emergency. Additionally, we used several statistical methods to assess the association between COVID-19 status and immune targets. Comparisons of immune response between non-COVID-19 and COVID-19 severe pneumonia are still scarce in the literature, even though they are mandatory to understand the distinctive pathogenesis of severe forms of COVID-19.


Higher plasma GM-CSF and CXCL10 were reported in COVID-19 patients and could represent the dysregulated immune response in COVID-19 patients driving the longer duration of mechanical ventilation in severe pneumonia. These cytokines could be considered as promising therapeutic targets.

Availability of data and materials

all data are available on demand.

Change history



Angiotensin-converting enzyme 2


Acute respiratory distress syndrome


Aspartate aminotransferase


C-C motif chemokine ligand


Cluster of differentiation


Coronavirus disease 2019


C-X-C motif chemokine ligand


CXC chemokine receptor 3


FMS-like tyrosine kinase 3 ligand


Granulocyte colony-stimulating factor


Granulocyte-macrophage colony-stimulating factor


Intensive care medicine


Interquartile range






INF-γ-induced protein 10


Mechanical ventilation


N-Terminal Fragment of the Prohormone Brain-Type Natriuretic Peptide

PaO2 :

FiO2: arterial pressure of oxygen/oxygen inspiratory fraction


Programmed death-ligand 1


Reverse transcriptase-polymerase chain reaction


Simplified acute physiology score II


Severe acute respiratory syndrome coronavirus 2


Sequential Organ Failure Assessment


transforming growth factor


Tumor necrosis factor


TNF-related apoptosis inducing ligand


Ventilator-acquired pneumonia


Whole blood stimulation


  1. 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(16):1574–81.

    Article  CAS  Google Scholar 

  2. Giamarellos-Bourboulis EJ, Netea MG, Rovina N, Akinosoglou K, Antoniadou A, Antonakos N, et al. Complex immune dysregulation in COVID-19 patients with severe respiratory failure. Cell Host Microbe. 2020;27(6):992-1000.e3.

    Article  Google Scholar 

  3. Minejima E, Bensman J, She RC, Mack WJ, Tuan Tran M, Ny P, et al. A dysregulated balance of proinflammatory and anti-inflammatory host cytokine response early during therapy predicts persistence and mortality in Staphylococcus aureus bacteremia. Crit Care Med. 2016;44(4):671–9.

    Article  CAS  Google Scholar 

  4. McElvaney OJ, McEvoy N, McElvaney OF, Carroll TP, Murphy MP, Dunlea DM, et al. Characterization of the inflammatory response to severe COVID-19 illness. Am J Respir Crit Care Med. 2020;202(6):812–21.

    Article  CAS  Google Scholar 

  5. Sims JT, Krishnan V, Chang C-Y, Engle SM, Casalini G, Rodgers GH, et al. Characterization of the cytokine storm reflects hyperinflammatory endothelial dysfunction in COVID-19. J Allergy Clin Immunol. 2020;S0091–6749(20):31242-2.

    Google Scholar 

  6. Del Valle DM, Kim-Schulze S, Huang H-H, Beckmann ND, Nirenberg S, Wang B, et al. An inflammatory cytokine signature predicts COVID-19 severity and survival. Nat Med. 2020;26(10):1636–43.

    Article  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  8. RECOVERY Collaborative Group, Horby P, Lim WS, Emberson JR, Mafham M, Bell JL, et al. Dexamethasone in Hospitalized Patients with Covid-19 - Preliminary Report. N Engl J Med. 2020.

    Article  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  10. Sinha P, Matthay MA, Calfee CS. Is a « Cytokine Storm » Relevant to COVID-19? JAMA Intern Med. 2020;180(9):1152–4.

    Article  CAS  Google Scholar 

  11. Chen LYC, Hoiland RL, Stukas S, Wellington CL, Sekhon MS. Confronting the controversy: interleukin-6 and the COVID-19 cytokine storm syndrome. Eur Respir J. 2020;56(4):2003006.

    Article  CAS  Google Scholar 

  12. Blot M, Bourredjem A, Binquet C, Piroth L. Is Interleukin 6 the right target in COVID-19 severe pneumonia? Am J Respir Crit Care Med. 2020.

    Article  Google Scholar 

  13. Leisman DE, Ronner L, Pinotti R, Taylor MD, Sinha P, Calfee CS, et al. Cytokine elevation in severe and critical COVID-19: a rapid systematic review, meta-analysis, and comparison with other inflammatory syndromes. Lancet Respir Med. 2020;S2213–2600(20):30404–5.

    Google Scholar 

  14. Remy KE, Mazer M, Striker DA, Ellebedy AH, Walton AH, Unsinger J, et al. Severe immunosuppression and not a cytokine storm characterize COVID-19 infections. JCI Insight. 2020;5(17):e140329.

    Article  Google Scholar 

  15. Hadjadj J, Yatim N, Barnabei L, Corneau A, Boussier J, Smith N, et al. Impaired type I interferon activity and inflammatory responses in severe COVID-19 patients. Science. 2020;369(6504):718–24.

    Article  CAS  Google Scholar 

  16. Channappanavar R, Fehr AR, Vijay R, Mack M, Zhao J, Meyerholz DK, et al. Dysregulated type I interferon and inflammatory monocyte-macrophage responses cause lethal pneumonia in SARS-CoV-infected mice. Cell Host Microbe. 2016;19(2):181–93.

    Article  CAS  Google Scholar 

  17. Blot M, Jacquier M, Manoha C, Piroth L, Charles P-E, Pneumochondrie study group. Alveolar SARS-CoV-2 viral load is tightly correlated with severity in COVID-19 ARDS. Clin Infect Dis. 2020.

    Article  PubMed  Google Scholar 

  18. 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(7):707‑10.

  19. 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(24):2957–63.

    Article  Google Scholar 

  20. 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(4):243–50.

    Article  CAS  Google Scholar 

  21. 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(23):2526‑33.

  22. Santagostino A, Garbaccio G, Pistorio A, Bolis V, Camisasca G, Pagliaro P, et al. An Italian national multicenter study for the definition of reference ranges for normal values of peripheral blood lymphocyte subsets in healthy adults. Haematologica. 1999;84(6):499–504.

    CAS  PubMed  Google Scholar 

  23. Troussard X, Vol S, Cornet E, Bardet V, Couaillac J-P, Fossat C, et al. Full blood count normal reference values for adults in France. J Clin Pathol. 2014;67(4):341–4.

    Article  Google Scholar 

  24. 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(1):632.

    Article  Google Scholar 

  25. Sainani KL. Introduction to principal components analysis. PM R. 2014;6(3):275–8.

    Article  Google Scholar 

  26. Cattell RB. The Scree Test For The Number Of Factors. Multivariate Behav Res. 1966;1(2):245–76.

    Article  CAS  Google Scholar 

  27. Cattell RB, Vogelmann S. A Comprehensive Trial Of The Scree And Kg Criteria For Determining The Number Of Factors. Multivariate Behav Res. 1977;12(3):289–325.

    Article  CAS  Google Scholar 

  28. Durbin J, Watson GS. Testing for serial correlation in least squares regression. II. Biometrika. 1951;38(1–2):159–78.

    Article  CAS  Google Scholar 

  29. White H. A heteroskedasticity-consistent covariance matrix estimator and a direct test for heteroskedasticity. Econometrica. 1980;48(4):817–38.

    Article  Google Scholar 

  30. Bhatraju PK, Ghassemieh BJ, Nichols M, Kim R, Jerome KR, Nalla AK, et al. Covid-19 in Critically Ill Patients in the Seattle Region-Case Series. N Engl J Med. 2020;382(21):2012–22.

    Article  CAS  Google Scholar 

  31. Gattinoni L, Coppola S, Cressoni M, Busana M, Rossi S, Chiumello D. COVID-19 Does Not Lead to a « Typical » Acute Respiratory Distress Syndrome. Am J Respir Crit Care Med. 2020;201(10):1299–300.

    Article  CAS  Google Scholar 

  32. Maggi E, Canonica GW, Moretta L. COVID-19: unanswered questions on immune response and pathogenesis. J Allergy Clin Immunol. 2020;146(1):18–22.

    Article  CAS  Google Scholar 

  33. Lang FM, Lee KM-C, Teijaro JR, Becher B, Hamilton JA. GM-CSF-based treatments in COVID-19: reconciling opposing therapeutic approaches. Nat Rev Immunol. 2020;20(8):507-514.

  34. Stanley E, Lieschke GJ, Grail D, Metcalf D, Hodgson G, Gall JA, et al. Granulocyte/macrophage colony-stimulating factor-deficient mice show no major perturbation of hematopoiesis but develop a characteristic pulmonary pathology. Proc Natl Acad Sci USA. 1994;91(12):5592–6.

    Article  CAS  Google Scholar 

  35. Zhou Y, Fu B, Zheng X, Wang D, Zhao C, Qi Y, et al. Aberrant pathogenic GM-CSF + T cells and inflammatory CD14 + CD16 + monocytes in severe pulmonary syndrome patients of a new coronavirus. bioRxiv. 2020;2020.02.12.945576.

  36. Mehta P, Porter JC, Manson JJ, Isaacs JD, Openshaw PJM, McInnes IB, et al. Therapeutic blockade of granulocyte macrophage colony-stimulating factor in COVID-19-associated hyperinflammation: challenges and opportunities. Lancet Respir Med. 2020;8(8):822–30.

    Article  CAS  Google Scholar 

  37. Oliviero A, de Castro F, Coperchini F, Chiovato L, Rotondi M. COVID-19 Pulmonary and Olfactory Dysfunctions: Is the Chemokine CXCL10 the Common Denominator? Neuroscientist. 2020;1073858420939033.

  38. Wong CK, Lam CWK, Wu AKL, Ip WK, Lee NLS, Chan IHS, et al. Plasma inflammatory cytokines and chemokines in severe acute respiratory syndrome. Clin Exp Immunol. 2004;136(1):95–103.

    Article  CAS  Google Scholar 

  39. Yang Y, Shen C, Li J, Yuan J, Wei J, Huang F, et al. Plasma IP-10 and MCP-3 levels are highly associated with disease severity and predict the progression of COVID-19. J Allergy Clin Immunol. 2020;146(1):119-127.e4.

    Article  Google Scholar 

  40. Lev S, Gottesman T, Levin GS, Lederfein D, Berkov E, Diker D, et al. Real-time IP-10 measurements as a new tool for inflammation regulation within a clinical decision support protocol for managing severe COVID-19 patients. medRxiv. 22 juill 2020;2020.07.21.20158782.

  41. Ichikawa A, Kuba K, Morita M, Chida S, Tezuka H, Hara H, et al. CXCL10-CXCR3 enhances the development of neutrophil-mediated fulminant lung injury of viral and nonviral origin. Am J Respir Crit Care Med. 2013;187(1):65–77.

    Article  CAS  Google Scholar 

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

    Article  CAS  Google Scholar 

  43. Neumann C, Scheffold A, Rutz S. Functions and regulation of T cell-derived interleukin-10. Semin Immunol. 2019;44:101344.

    Article  CAS  Google Scholar 

  44. Sterner RM, Sakemura R, Cox MJ, Yang N, Khadka RH, Forsman CL, et al. GM-CSF inhibition reduces cytokine release syndrome and neuroinflammation but enhances CAR-T cell function in xenografts. Blood. 2019;133(7):697–709.

    Article  CAS  Google Scholar 

Download references


The authors thank all the investigators of the Lymphonie study group (Pascal Andreu, François Aptel, Marie Labruyère, Sébastien Prin (Department of Intensive Care, Dijon Bourgogne University Hospital, Dijon, France); Guillaume Beltramo, Philippe Bonniaud (Department of Pneumology, Dijon Bourgogne University Hospital, Dijon, France); Philip Bielefeld, Hervé Devilliers (Department of Internal Medicine and Systemic Diseases, Dijon Bourgogne University Hospital, Dijon, France); Bernard Bonnotte (Department of Internal Medicine and Clinical Immunology, Dijon Bourgogne University Hospital, Dijon, France); Marielle Buisson (Infectious Diseases Department, Dijon Bourgogne University Hospital, Dijon, France); Alain Putot (Geriatrics Internal Medicine Department, Dijon Bourgogne University Hospital, Dijon, France). We thank 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;, 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.


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), and by crowdfunding (

Author information

Authors and Affiliations




Concept and design: MB, AB, CB, LP; Recruitment of patients: MB, JPQ, MN, MG, AL, AD,AG, SMM, JB, BB, PEC, PC; Acquisition, analysis, or interpretation of data: MB, JBB, JG, SM, AB, CB, LP; Drafting of the manuscript: MB, AB; Critical revision: MB, JBB, JPQ, AB, MN, JG, SM, MG, AL, AD, AG, SMM, JB, BB, PEC, PC, CB, LP; Supervision: CB, LP. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Mathieu Blot.

Ethics declarations

Ethics approval and consent to participate

Approval was obtained from the ethics committee (Comité de Protection des Personnes SUD MEDITERRANEE V; 2017-A03404-49). Oral consent was obtained from the patient or their legal representatives.

Consent for publication

all authors consent for publication.

Competing interests

we declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

The original version of this article was revised: the members of the LYMPHONIE study group were added as collaborators in the metadata of the article.

Supplementary information

Additional file 1: Table S1.

Microbiological etiologies of severe pneumonia from 63 patients. Table S2. Plasma cytokine concentrations. Table S3. Blood leukocyte cytokine production on ex vivo stimulation. Table S4. Spearman correlation between plasma cytokine concentrations and severity (PaO2:FiO2 ratio, SOFA score) or outcome (duration of mechanical ventilation). Table S5. Spearman correlation between blood leukocyte cytokine production on ex vivo stimulation (Whole Blood Stimulation (WBS)) and severity (PaO2:FiO2 ratio, SOFA score) or outcome (duration of mechanical ventilation). Figure S1. Principal component analysis screen plot (retained dimension = 4) Cattell scree plot criterion retained the first four independent factors; which were clinically interpretable and together preserved 53.47% of all the information contained in the 65 correlated original variables.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Blot, M., Bour, JB., Quenot, J.P. et al. The dysregulated innate immune response in severe COVID-19 pneumonia that could drive poorer outcome. J Transl Med 18, 457 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: