Introductory paragraph
The mechanisms that underpin COVID-19 disease severity, and determine the outcome of infection, are only beginning to be unraveled. The host inflammatory response contributes to lung injury, but circulating mediators levels fall below those in classical ‘cytokine storms’. We analyzed serial plasma samples from 619 patients hospitalized with COVID-19 recruited through the prospective multicenter ISARIC clinical characterization protocol U.K. study and 39 milder community cases not requiring hospitalization. Elevated levels of numerous mediators including angiopoietin-2, CXCL10, and GM-CSF were seen at recruitment in patients who later died. Markers of endothelial injury (angiopoietin-2 and von-Willebrand factor A2) were detected early in some patients, while inflammatory cytokines and markers of lung injury persisted for several weeks in fatal COVID-19 despite decreasing antiviral cytokine levels. Overall, markers of myeloid or endothelial cell activation were associated with severe, progressive, and fatal disease indicating a central role for innate immune activation and vascular inflammation in COVID-19.
Main text
Fatal COVID-19 is associated with acute respiratory distress syndrome and raised systemic inflammatory markers including IL-6 and C-reactive protein, often accompanied by neutrophilia and lymphopenia 1. The beneficial effect of corticosteroid treatment in severe disease highlights the role of steroid-responsive inflammation in pathogenesis 2, 3, and post-mortem studies report pulmonary vessel vasculitis (most commonly myeloid cells) and microthrombosis in fatal COVID-19 4, 5, 6, 7. The virus-induced inflammatory state has laboratory features that resemble secondary haemophagocytic lymphohistiocytosis (sHLH) 8, 9, 10 but the exact pattern and severity of inflammatory responses has been only partially characterized. Levels of some inflammatory mediators, including IL-6, are elevated in COVID-19, but are typically ten times lower than those reported in acute respiratory distress syndrome (ARDS) and sepsis 11, 12, 13, suggesting that other factors may play a major role in COVID-19 severity. Host genetic factors may also influence disease severity, with polymorphisms in several regions, including the interferon pathway genes IFNAR2 and OAS1/2/3 recently associated with enhanced disease severity 14. Identification of such genetic of inflammatory factors may define a ‘treatable trait’ 15, allowing both stratification of patients likely to benefit from therapies such as dexamethasone and targeted biological anti-cytokine therapies, and design of novel therapeutics targeting causative pathways.
Early clinical studies of COVID-19 identified elevated neutrophil counts and lymphopenia in peripheral blood 1, 16, predominantly seen in late-stage disease and of limited prognostic value. Peripheral blood neutrophilia is also seen in other severe respiratory viral 17 and bacterial 18 infections, suggesting that this is not a unique feature of COVID-19. Elevated levels of D-dimer, a product of fibrin-degradation associated with thrombosis and inflammation, have also been observed in COVID-19 16, consistent with systemic inflammation and the high frequency of macrovascular thrombotic complications in severe cases 7, 19. Post-mortem studies show that thromboses and microthrombi within pulmonary vessels are common in fatal COVID-19 and are associated with endothelial responses distinct from those that occur during fatal influenza A virus infection 4, 5, 7, 20. However, the thrombotic aspects of life-threatening COVID-19, and the interaction of this process with cytokine release have hitherto been described in relatively small groups of cases, from single-center studies, or with a narrow range of disease severities.
Within the ISARIC4C study we obtained clinical data and 1,047 plasma samples from 619 hospitalized patients with COVID-19 21, 22. Given the large number of cases, patients from the ISARIC4C database could be stratified into five levels of severity according to their peak illness, in line with the World Health Organization COVID-19 ordinal scale 23 (Supplementary Table 1): (1) no oxygen requirement (Severity 3, n=169); (2) patients requiring oxygen by face mask (Severity 4, n=143); (3) patients requiring high-flow nasal cannulae, a continuous positive airway pressure mask or other non-invasive ventilation (Severity 5 n=99); (4) patients requiring invasive mechanical ventilation (Severity 6/7, n=113); and (5) fatal COVID-19 (Severity 8, n=95). The median duration of symptoms prior to study enrollment was similar in all groups: Severity 3, 7 days; Severity 4, 9 days; Severity 5, 11 days; Severity 6/7, 11 days; and Severity 8, 8 days. Some differences in routinely performed clinical hematology and biochemistry measures were evident between clinical outcome groups at the time of study enrolment: Lymphopenia was evident in groups 6/7 and 8, relative to 3, alongside neutrophilia in 6/7 and 8 relative to 3 and 4 (Supplementary Fig. 1a and 1b, respectively). No differences between groups were observed in ferritin levels, whilst LDH was elevated in groups 5, 6/7, and 8 relative to 3 (Supplementary Fig. 1c and 1d, respectively). Procalcitonin levels were elevated in group 8 relative to 3 and 4, and in group 6/7 relative to 4 (Supplementary Fig. 1e). Partial HScores 24 were calculated (fever, cytopenia, ferritin, triglycerides, and AST) but the only significant difference between groups was between 6/7 and 4, indicating that sHLH is unlikely to be the predominant pathophysiological mechanism in life-threatening COVID-19 (Supplementary Fig. 1f). The ISARIC4C mortality scores 25 for these patients demonstrated an elevated risk of mortality, calculated from admission data, in those that would progress to fatal disease (group 8) relative to all other groups (Supplementary Fig. 1g), though there was considerable overlap between all groups. Together, these data indicated limited clinical or biochemical differences between patient outcome groups at the time of hospital admission.
a) Peripheral blood lymphocyte count, b) neutrophil count, c) ferritin levels, d) lactate dehydrogenase (LDH) levels, e) procalcitonin levels, f) partial HScores, and g) ISARIC4C mortality scores at the time of enrolment in hospitalized patients with COVID-19 that would: not require oxygen support (‘3’, n=9-93); require an oxygen face mask (‘4’, n=22-71); require non-invasive ventilation or high-flow nasal cannulae (‘5’, n=15-63); require invasive mechanical ventilation (‘6/7’, n=19-91); or progress to fatal disease (‘8’, n=15-63). Data were analyzed for statistical significance using Kruskal-Wallis tests with Dunn’s tests for multiple comparisons between all groups. *P<0.05, **P<0.01, ***P<0.001, ****P<0.0001.
We hypothesized that differences in the levels of plasma inflammatory mediators would reflect the nature and scale of immunopathology in COVID-19 and would associate with different disease outcomes. We therefore quantified 33 mediators in all available plasma samples using panels designed to study a broad range of mediators that could be broadly categorized as having roles in antiviral immunity, inflammation, or coagulation 16, 19. Analysis of plasma mediator levels at the time of enrolment distinguished 3 clusters of patients that were associated with distinct patterns of mediator levels (Fig. 1). The first of these clusters was enriched in patients from groups 6/7 and 8 and was associated with higher levels of CXCL10, GM-CSF, D-dimer, and vWF-A2. The second cluster contained a more diverse mixture of severities and had a more pronounced pattern of coagulation factor XIV and angiopoietin-2 containing mediator clusters, but lower levels of the CXCL10 containing mediator cluster. The third patient cluster had lower levels of the CXCL10, D-dimer, and coagulation factor XIV containing mediator clusters, but had a more varied pattern of other mediators including IL-6Rα, VEGF-D, and IL-4. Interestingly, this analysis did not indicate any obvious patterns of age, symptom duration (onset), or sex in these plasma mediator levels. This analysis shows that, at the time of enrolment, different COVID-19 outcome groups were already identifiable and associated with distinct patterns of inflammatory mediators and that markers such as D-dimer, EN-RAGE, CXCL10, and GM-CSF were particularly associated with enhanced disease severity. However, entry to the study was determined by hospitalization which will be influenced by predisposing factors; these factors may therefore not be evident in the data that we accumulate.
Clustered heatmap of 33 immune mediators in plasma samples collected from patients hospitalized with COVID-19 at the time of study enrolment. Values for each mediator were scaled and rows and columns were split by K-means clustering. Each patients’ column is additionally annotated with data on disease outcome (“Severity”) as one of the following outcome groups: not requiring oxygen support (‘3’, n=128), requiring oxygen via a face mask (‘4’, n=103), requiring non-invasive ventilation or high-flow nasal canulae (‘5’, n=78), requiring invasive mechanical ventilation (‘6/7’, n=87) or fatal disease (‘8’, n=69). Columns are additionally annotated with patient age, sex and duration of illness at the time of sample collection (“Onset”).
To further explore the relationship between the mediator levels and severity we analyzed plasma from 15 healthy controls (7 males, median age 55, range 45-71) and 39 individuals recruited 7 days after a SARS-CoV-2 positive PCR test who did not require hospitalization (15 males, median age 43, range 27-62, termed group ‘1/2’ as per the WHO scale 23) and related these to hospitalized patients. At the time of enrollment, numerous differences were evident between hospitalized COVID-19 patients and the control groups, along with many differences across the clinical outcome groups in hospitalized patients (Fig. 2 and Supplementary Fig. 2). In contrast to other reports 26 we found no evident deficiency in IFN-α levels in those with severe disease (Fig. 2a). IFN-γ was elevated in hospitalized COVID-19 patients relative to both healthy controls (HC) and group 1/2 (Fig. 2b) and was elevated in the most severe outcome groups, relative to lower severity grades. The interferon-induced chemokine CXCL10 was also substantially elevated in all hospitalized COVID-19 cases relative to the control groups, with the most pronounced increases in groups 6/7 and 8 (Fig. 2c). These results are in contrast to the decreased ISG gene expression in peripheral blood samples from patients with severe COVID-19 26, showing that the gene expression pattern from blood does not necessarily reflect the directly measured levels of gene product. We speculate that the abundance of IFN-γ and CXCL10 results from release from the site of disease rather than from circulating cells, though anti-IFN autoantibodies 27 and polymorphisms in IFN signaling 14 may influence this pathway in some patients.
Mediator levels were quantified from plasma collected at the point of study enrolment from hospitalized patients with COVID-19 that would: not require oxygen support (‘3’, n=128), require an oxygen face mask (‘4’, n=103), require non-invasive ventilation or high-flow nasal cannulae (‘5’, n=78), require invasive mechanical ventilation (‘6/7’, n=87) or progress to fatal disease (‘8’, n=69). Data were analyzed for statistical significance using Kruskal-Wallis tests with Dunn’s tests for multiple comparisons between all groups. *P<0.05, **P<0.01, ***P<0.001, ****P<0.0001.
Plasma samples from the time of study enrolment were analyzed for levels of the antiviral cytokines a) IFN-α, b) IFN-γ, and c) the interferon-induced chemokine CXCL10 in healthy control (HC, n=15), patients with COVID-19 not requiring hospitalization (‘1/2’, n=39), and hospitalized patients with COVID-19 that would: not require oxygen support (‘3’, n=32-128), require an oxygen face mask (‘4’, n=23-103), require non-invasive ventilation or high-flow nasal cannulae (‘5’, n=19-78), require invasive mechanical ventilation (‘6/7’, n=19-87) or progress to fatal disease (‘8’, n=14-69). Mediators associated with coagulation and endothelial injury were also quantified in these plasma samples; d) D-dimer, e) Angiopoietin-2, and f) von-Willebrand factor A2 (vWF-A2). Similarly, mediators associated with inflammation were quantified: g) IL-6; h) GM-CSF; and i) EN-RAGE/S100A12. Data were analyzed for statistical significance using Kruskal-Wallis tests with Dunn’s tests for multiple comparisons between all groups. *P<0.05, **P<0.01, ***P<0.001, ****P<0.0001.
The fibrin degradation product D-dimer has been reported to be elevated in severe COVID-19 16, implicating thrombosis in disease severity 4, 5, 20. In agreement with these reports, D-dimer was elevated in all hospitalized groups, but little difference was observed between the severity groups at the time of enrolment (Fig. 2d). Given reports of the association between COVID-19 mortality and pulmonary vasculitis 4, we hypothesized that endothelial injury may be a feature of COVID-19, potentially triggering coagulation and the thrombotic complications common in severe disease 19, 28. Indeed, levels of angiopoietin-2, a marker of endothelial injury, were elevated in all hospitalized patients relative to both control groups (Fig. 2e), with levels 5.6-fold higher in the mildest hospitalized patients (group 3, median=1983pg/ml) than HCs (median=352pg/ml). Angiopoetin-2 levels were also significantly elevated in groups 6/7 and 8 relative to all other hospitalized COVID-19 outcome groups (Fig. 2e). As both angiopoietin-2 and vWF-A2 can enter the blood plasma through exocytosis of endothelial cell Weibel-Palade bodies 29, we also quantified vWF-A2, which was similarly elevated in hospitalized COVID-19 patients (Fig. 2f). In line with these markers of endothelial injury and thrombosis, thrombomodulin, vWF-A2, and endothelin-1 were also elevated in COVID-19, predominantly in those most severe patient outcome groups (Supplementary Fig. 2). Elevations in these prothrombotic mediators were not counteracted by levels of the inhibitors angiopoietin-1 or soluble Tie2, which were not significantly different between the tested groups (Fig. S2). These results suggested that endothelial injury and coagulation are common features of patients hospitalized with COVID-19 and that these are most pronounced in severe and fatal COVID-19.
In line with other reports 1, 12, we found that IL-6 was also significantly elevated in most hospitalized groups relative to the controls (Fig. 2g), with a stepwise increase in levels with escalating severity. IL-6 levels in groups 6/7 and 8 were significantly elevated above all other groups (all P<0.0001, Fig. 2g). In agreement with the association of a strong inflammatory response with COVID-19 severity, GM-CSF was similarly elevated in all hospitalized groups, relative to controls and was most pronounced in the groups 6/7 and 8 (Fig. 2h). Numerous other inflammatory cytokines and chemokines showed similar results including TNF-α, IL-2, GDF-15, G-CSF, and VEGF-D (Supplementary Fig. 2). EN-RAGE/S100A12 has previously been characterized as a marker of respiratory damage in ARDS 30 and indeed was elevated in groups 6/7 and 8 relative to most others (Fig. 2i). The neutrophil chemokine IL-8 (CXCL8) was similarly elevated in severe disease, as was the neutrophil gelatinase associated lipocalin (LCN-2/NGAL) (Supplementary Fig. 2), in line with the reported association between blood neutrophilia and severity 16 also seen in this cohort (Supplementary Fig. 1b).
Other immunological mediators (IL-6Rα, IL-13, IL-17) were not significantly different between groups, indicating that only limited aspects of the immune repertoire were active in COVID-19. Interestingly, IL-4 levels were lower in the non-severe disease outcome groups (3, 4, and 5) relative to both the control groups and the severe disease groups 6/7 and 8 (Supplementary Fig. 2), indicating that suppression of the normal levels of type-2 cytokines may be associated with milder COVID-19 disease, and that this mechanism is lost in severe disease. Similarly, IL-12p70, commonly released by antigen presenting cells (APCs) 31, was decreased in all hospitalized cases relative to the HCs and group 1/2 (Figure S2), possibly owing to the trafficking of APCs to the site(s) of viral infection.
To determine the strength of the relationships between these individual plasma mediators we performed a hierarchical correlation matrix analysis of mediators from plasma samples collected at the time of study enrolment. This identified a strongly correlated cluster of inflammatory mediators including GM-CSF, CXCL10, vWF-A2, and IL-6 (Fig. 3a); increases in which were commonly associated with the most severe COVID-19 outcome groups. Given the strong association between age and COVID-19 severity 22, and reports of increased inflammatory responses in males relative to females with COVID-19 32 we investigated the influence of these demographic factors on plasma mediators levels in hospitalized patients. As the major effect in our cluster analysis was severity (Fig. 1), we further stratified each of these severity groups by age (≥ or < 70 years of age) and sex, to better account for the influence of disease severity on plasma mediator levels. Following adjustment for multiple testing, no mediator was found to be statistically different between males and females within each outcome group (Supplementary Fig. 3). By contrast, several differences were evident between those aged ≥70 and <70 years, including elevated levels of D-dimer, CXCL10, and GM-CSF in those aged ≥70 years; IFN-γ levels were, by contrast, greater in younger patients within severity group 4 (Fig. 3b and Supplementary Fig. 3).
Heatmap of false-discovery rate adjusted P-values for each plasma mediator between males and females (“Sex”) and those aged ≥70 years and <70 years (“Age”) within each disease outcome group (‘8’=Red, ‘6/7’=Orange, ‘5’=Purple, ‘4’=Dark blue, ‘3’=Cyan). Data were analyzed using Mann-Whitney U tests with P-value adjustment for false discovery rate.
a) Correlogram of the association between plasma mediator levels at the time of enrolment in all patients hospitalized with COVID-19 (n=465). b) Inflammatory mediator levels within an outcome group, stratified as those ≥ or < than 70 years of age. Data in panel a were analyzed using Spearman’s rank correlations with correction for multiple testing; significant correlations are denoted by a circle, the color of which denotes the Spearman’s R value. Data in panel b were analyzed using Mann-Whitney U tests with P-value adjustment for false discovery rate.
We next sought to determine the changes in levels of some key plasma mediators from the time of enrolment over the course of disease, by relating data to the patient reported duration of symptoms at the time of each sample collection, including consecutive samples collected from individual patients. This analysis indicated that many mediators were stable over the time-course of hospitalization, supporting the validity of using samples from the time of enrolment to study the immunologic basis of COVID-19. However, some mediators did change over time; for example, there was a gradual decrease in IFN-γ and CXCL10 over time in most groups (Supplementary Fig. 4), including group 8 (Fig. 4a and 4b, respectively). By contrast some other mediators remained elevated or appeared to increase over the duration of symptoms in group 8, including angiopoietin-2 and D-dimer (Fig. 4c and 4d, respectively). Similarly, the inflammatory mediators GM-CSF and EN-RAGE remained elevated or increased in group 8 in the latter stages of disease (Fig. 4e and 4f, respectively). Together, these results indicated that the most severe outcomes of COVID-19 disease were associated with persistent coagulation and inflammation, even as IFN levels declined.
All data within each severity group was related to the duration of symptoms at the time of sample collection (“Onset to sample”, measured in days) for each plasma mediator. Generalized additive modelling was used to fit a restricted cubic spline which is plotted together with the standard error (grey).
Plasma levels of a) IFN-γ, b) CXCL10, c) Angiopoietin-2, d) D-dimer, e) GM-CSF, and f) EN-RAGE/S100A12 over the course of disease in patients with fatal COVID-19. Plasma mediator levels of g) IL-2, h) IL-6, i) GM-CSF, j) D-dimer, and k) von-Willebrand factor A2 (vWF-A2) within the first 4 days of symptom onset in patients in severity groups 6/7 or 8 (“Severe”, n=22) and groups 3, 4, or 5 (“Non-Severe”, n=54). Linear regressions with 95% confidence intervals are shown in panels a-f. Data in panels g-k were analyzed for statistical significance using Mann-Whitney U tests, where thick horizontal dashed lines denote the median values and thin horizontal dashed lines denote the interquartile ranges.
Finally, we hypothesized that differences in plasma mediator levels between patients with Severe (groups 6/7 and 8) and Non-severe (groups 3, 4, and 5) COVID-19 would be apparent within the first few days of symptoms. Indeed, within the first 4 days of symptoms several mediators were significantly elevated in the Severe group, relative to Non-severe, including IL-2, IL-6, and GM-CSF (all P<0.0001, Fig. 4g-i, respectively), indicating a pronounced inflammatory response early in Severe disease. Similarly, many markers of coagulation and endothelial injury were elevated in Severe, relative to Non-severe, including D-dimer and vWF-A2 (P<0.0001, Fig. 4j and 4k, respectively), in addition to angiopoietin-2 and IL-1α (which can be activated by thrombin 33) (Supplementary Fig. 5). By comparison the lung damage-associated marker EN-RAGE 30 was not significantly different between the Severe and Non-severe groups in the first 4 days of symptoms (P=0.098, Supplementary Fig. 5). Together, these data indicated that severe COVID-19 is associated with elevated levels of plasma mediators indicative of coagulation, endothelial activation and a broad inflammatory response including CXCL10, GM-CSF, and IL-6. These differences were apparent within the first days of symptoms, while markers of lung damage may only become elevated later in disease, potentially indicating a pathological role for these processes and a window of opportunity for early immunomodulation to prevent significant lung damage.
Levels of immune mediators collected within the first 4 days of symptom onset in patients in the groups 6/7 or 8 (“Severe”, n=22) and groups 3, 4, or 5 (“Non-Severe”, n=54). Data were analyzed for statistical significance using Mann-Whitney U tests, where thick horizontal dashed lines denote the median values and thin horizontal dashed lines denote the interquartile ranges.
While markers of fibrinolysis have previously been associated with disease severity 16 and thrombosis is common in severe and fatal COVID-19 4, 5, 20 the causes of this manifestation of severe disease are not known. We demonstrate that increasing disease severity is associated with broad elevations in inflammatory mediator levels, alongside a signature of endothelial injury. This signal was most pronounced in fatal COVID-19 and was apparent even in the early stages of disease.
The elevation of angiopoietin-2, thrombomodulin, and vWF-A2 in fatal COVID-19 cases provides evidence for the involvement of endothelial injury in COVID-19 severity. Endothelial injury following inflammatory damage, including the increasingly recognized pulmonary artery vasculitis 4, 20 in COVID-19, may result in the initiation of a pro-coagulant role for these cells 34. Alternatively, this response could be triggered by direct viral infection of vascular cells (though this has yet to be conclusively determined 34 viral replication in non-respiratory tissues is commonly observed at post-mortem 4, 7); or thrombin mediated activation of IL-1α 33. This pro-coagulant role could lead to the deposition of microthrombi, evident in COVID-19 4, the development of features of disseminated intravascular coagulopathy (DIC) and ultimately elevated levels D-dimer through the degradation of fibrin rich thrombi 28. Neutrophilic inflammation could have an etiological role in endothelial injury though neutrophilia is predominantly a feature of the later phases of COVID-19 1, while endothelial injury was evident in the first days of symptoms. However, the continued thrombosis in late stage fatal COVID-19 may result from neutrophil mediated coagulation, observed in other settings 35, 36, 37 and recently demonstrated in COVID-19 38. Combined, these results indicate a multiplicity of possible pro-coagulant triggers that may contribute to pathology at different stages of disease.
We found that the antiviral immune mediator CXCL10 and the myeloid cell growth factor GM-CSF, were strikingly elevated in fatal cases of COVID-19. This is confirmed by a recent report describing the potential utility of CXCL10 as an early prognostic marker of COVID-19 severity 39. An influx of monocytes/macrophages has been described in the lung parenchyma in fatal COVID-19, combined with a mononuclear cell pulmonary artery vasculitis 6, and presence of pro-inflammatory monocyte-derived macrophages in bronchoalveolar lavage fluid from patients with severe COVID-19 4, 40. The elevation of CXCL10 and GM-CSF in severe disease reported here could contribute to monocyte recruitment and activation leading to this vasculitis, alongside the role of GM-CSF in the recruitment of neutrophils to the pulmonary vasculature 41.
Large scale randomized clinical trials for IL-6 signaling antagonists are on-going, though early results of the COVACTA trial of Tocilizumab found no improvement in clinical status or mortality 42. Small scale studies of anti-GM-CSF have shown promising results 43, 44 but require formal testing in a clinical trial. Given the role of GM-CSF in granulopoiesis and enhancement of neutrophil survival, alongside the neutrophil activation observed in late stage fatal COVID-19, these trials may inform our understanding of the importance of this pathway in COVID-19 immunopathogenesis 45. While early studies demonstrated elevated GM-CSF levels in both ICU and non-ICU treated COVID-19 patients 1, we now demonstrate a positive association with disease severity and outcome, in agreement with reports of elevated frequencies of GM-CSF+ Th1 cells in patients with COVID-19 requiring ICU treatment 46.
While many cytokines and other inflammatory mediators were most significantly elevated in fatal and critical COVID-19, these data do not necessarily support the concept of a “cytokine storm” in COVID-19 12, 13. While some elements, such as elevated IL-6 and ferritin levels (reported in other studies, but not seen here) 8, 9, 10, are reminiscent of sHLH, the relatively gradual clinical progression and persistent elevation of some cytokines, even during the early stages of symptomatic disease, are uncommon amongst conditions associated with cytokine storms such as toxic-shock syndrome and bacterial sepsis.
To our knowledge, this is to date the largest study of inflammatory responses in COVID-19. The multicenter nature of ISARIC4C adds to the ability to interpret and apply these results to other settings. However, further studies are needed to determine the prognostic value of these key plasma biomarkers, including multivariable analyses of biological data alongside clinical and demographic data. This detailed level of analysis may also enable the phenotyping of patients most likely to respond to individual therapies. Future analyses should focus on the biological features of patients that respond to therapeutic interventions, such as dexamethasone 2, 3, to enable mechanistic insight and targeting of treatment. The clear distinction between patients that would progress to severe COVID-19 and those that would not, even in the earliest stages of disease, indicates that early therapeutic intervention may be crucial to limit mortality. Overall, these data indicate an early inflammatory response in COVID-19, most prominent in those who will later suffer severe or fatal disease. These responses may enable the development of prognostic biomarkers, inform our understanding of immunopathogenesis in COVID-19 and enable novel approaches for therapeutic intervention.
Supplementary Methods
Study design and setting
The ISARIC WHO Clinical Characterization Protocol for Severe Emerging Infections in the UK (CCP-UK) is an ongoing prospective cohort study of hospitalized patients with COVID-19, which is recruiting in 258 hospitals in England, Scotland, and Wales (National Institute for Health Research Clinical Research Network Central Portfolio Management System ID: 14152) 47. The protocol, revision history, case report form, patient information leaflets, consent forms and details of the Independent Data and Material Access Committee are available online 48 and published previously 22.
Participants
Hospitalized patients with PCR-proven or high likelihood of SARS-CoV-2 infection were recruited, including both patients with community- and hospital-acquired COVID-19. This study analyzed plasma from blood samples obtained on the day of enrolment to the study (day 1, Tier 1) and additional serial samples obtained following a sampling schedule (Tier 2) harmonized with international investigators to allow meaningful comparison of results between studies 21. Healthy controls were recruited prior to December 2019 under approval from the London – Fulham Research Ethics Committee (REC) (reference 14/LO/1023) or from healthy donors following informed consent from a sub-collection of the Imperial College Healthcare NHS Trust National Institute for Health Research Imperial Biomedical Research Centre Tissue Bank. Use of the sub-collection was approved by the Tissue Bank Ethics Committee (Approval R12023). Samples from community managed COVID-19 cases were collected through a subproject of Imperial College London Communicable Disease Research Tissue Bank, under approval from the south central Oxford REC (reference 15/SC/0089).
HScores
To calculate partial HScores 24, ferritin, triglyceride and AST measurements from this study were combined with recorded results from case report forms for temperature and routine hemoglobin, white cell counts, and platelet counts.
Immunoassays
IFN-γ, TNF-α, IL-1β, IL-2, IL-4, IL-6, CXCL8/IL-8, IL-10, IL-12p70 and IL-13 were quantified using MSD (Mesoscale Diagnostics, Rockville, Maryland, USA) V-Plex proinflammatory plates on a SQ120 Quickplex instrument. IL-1α, IL-1ra, IL-6Rα, angiopoetin-1, angiopoetin-2, endothelin-1, VEGF-D, D-dimer, thrombomodulin, Tie2, von-Willebrand Factor-A2 (vWF-A2), G-CSF, GM-CSF, IL-17A, LCN2/NGAL, CXCL10/IP-10, CCL2, CCL3, CCL4 and CCL5 were quantified using a Bio Plex 200 instrument (Bio-Rad, Hercules, California, USA) with custom Luminex panel kits from Biotechne (Minneapolis, Minnesota, USA) and MilliporeSigma (Burlington, Massachusetts, USA). IFN-α was quantified using Quanterix (Billerica, Massachusetts, USA) IFN-α assay kits on the SIMOA platform. All values at or below the lower limit of detection (LLOD) were replaced with the geometric mean of the lower limits of detection across plates for each assay.
Statistical analyses
Statistical analyses used GraphPad Prism v8.3.0 (GraphPad, La Jolla, California, USA) R version 3.6.1 and Python 3.7.3 with Pandas 1.0.3 and Seaborn 0.10.0. Non-parametric mediator data (as determined by D’Agostino and Pearson normality test) were analyzed by ANOVA using Kruskal-Wallis tests with Dunn’s test for multiple comparisons of patient groups within in time group. Non-parametric two-way analyses were performed using Mann-Whitney U tests. Correlation matrix analysis was performed using the R packages ggplot2 and ggcorrplot and Spearman’s test for correlation of non-parametric data, after P-value adjustment for multiple testing. The false discovery rate, or expected proportion of discoveries which are falsely rejected, was controlled using the methods of Benjamini and Hochberg. Heatmaps of scaled plasma mediator data were generated using the ComplexHeatmap package in R with rows and columns split by K-means clustering and dendrograms based on Ward’s minimum variance method (ward.D2) and Spearmans rank correlations. For heatmap analyses missing values were imputed by predictive mean matching using the Multivariate Imputation by Chained Equations (MICE) package 49.
The ISARIC4C Investigators
Consortium Lead Investigator
J Kenneth Baillie, Chief Investigator: Malcolm G Semple, Co-Lead Investigator: Peter JM Openshaw. ISARIC Clinical Coordinator: Gail Carson.
Co-Investigators
Beatrice Alex, Benjamin Bach, Wendy S Barclay, Debby Bogaert, Meera Chand, Graham S Cooke, Annemarie B Docherty, Jake Dunning, Ana da Silva Filipe, Tom Fletcher, Christopher A Green, Ewen M Harrison, Julian A Hiscox, Antonia Ying Wai Ho, Peter W Horby, Samreen Ijaz, Saye Khoo, Paul Klenerman, Andrew Law, Wei Shen Lim, Alexander J Mentzer, Laura Merson, Alison M Meynert, Mahdad Noursadeghi, Shona C Moore, Massimo Palmarini, William A Paxton, Georgios Pollakis, Nicholas Price, Andrew Rambaut, David L Robertson, Clark D Russell, Vanessa Sancho-Shimizu, Janet T Scott, Thushan de Silva, Louise Sigfrid, Tom Solomon, Shiranee Sriskandan, David Stuart, Charlotte Summers, Richard S Tedder, Emma C Thomson, AA Roger Thompson, Ryan S Thwaites, Lance CW Turtle, Maria Zambon.
Project Managers
Hayley Hardwick, Chloe Donohue, Ruth Lyons, Fiona Griffiths, Wilna Oosthuyzen.
Data Analysts
Lisa Norman, Riinu Pius, Tom M Drake, Cameron J Fairfield, Stephen Knight, Kenneth A Mclean, Derek Murphy, Catherine A Shaw.
Data and Information System Manager
Jo Dalton, Michelle Girvan, Egle Saviciute, Stephanie Roberts, Janet Harrison, Laura Marsh, Marie Connor, Sophie Halpin, Clare Jackson, Carrol Gamble.
Data integration and presentation
Gary Leeming, Andrew Law, Murray Wham, Sara Clohisey, Ross Hendry, James Scott-Brown.
Material Management
William Greenhalf, Victoria Shaw, Sarah McDonald.
Patient engagement
Seán Keating
Outbreak Laboratory Staff and Volunteers
Katie A. Ahmed, Jane A Armstrong, Milton Ashworth, Innocent G Asiimwe, Siddharth Bakshi, Samantha L Barlow, Laura Booth, Benjamin Brennan, Katie Bullock, Benjamin WA Catterall, Jordan J Clark, Emily A Clarke, Sarah Cole, Louise Cooper, Helen Cox, Christopher Davis, Oslem Dincarslan, Chris Dunn, Philip Dyer, Angela Elliott, Anthony Evans, Lorna Finch, Lewis WS Fisher, Terry Foster, Isabel Garcia-Dorival, Willliam Greenhalf, Philip Gunning, Catherine Hartley, Antonia Ho, Rebecca L Jensen, Christopher B Jones, Trevor R Jones, Shadia Khandaker, Katharine King, Robyn T. Kiy, Chrysa Koukorava, Annette Lake, Suzannah Lant, Diane Latawiec, L Lavelle-Langham, Daniella Lefteri, Lauren Lett, Lucia A Livoti, Maria Mancini, Sarah McDonald, Laurence McEvoy, John McLauchlan, Soeren Metelmann, Nahida S Miah, Joanna Middleton, Joyce Mitchell, Shona C Moore, Ellen G Murphy, Rebekah Penrice-Randal, Jack Pilgrim, Tessa Prince, Will Reynolds, P. Matthew Ridley, Debby Sales, Victoria E Shaw, Rebecca K Shears, Benjamin Small, Krishanthi S Subramaniam, Agnieska Szemiel, Aislynn Taggart, Jolanta Tanianis-Hughes, Jordan Thomas, Erwan Trochu, Libby van Tonder, Eve Wilcock, J. Eunice Zhang.
Local Principal Investigators
Kayode Adeniji, Daniel Agranoff, Ken Agwuh, Dhiraj Ail, Ana Alegria, Brian Angus, Abdul Ashish, Dougal Atkinson, Shahedal Bari, Gavin Barlow, Stella Barnass, Nicholas Barrett, Christopher Bassford, David Baxter, Michael Beadsworth, Jolanta Bernatoniene, John Berridge, Nicola Best, Pieter Bothma, David Brealey, Robin Brittain-Long, Naomi Bulteel, Tom Burden, Andrew Burtenshaw, Vikki Caruth, David Chadwick, Duncan Chambler, Nigel Chee, Jenny Child, Srikanth Chukkambotla, Tom Clark, Paul Collini, Catherine Cosgrove, Jason Cupitt, Maria-Teresa Cutino-Moguel, Paul Dark, Chris Dawson, Samir Dervisevic, Phil Donnison, Sam Douthwaite, Ingrid DuRand, Ahilanadan Dushianthan, Tristan Dyer, Cariad Evans, Chi Eziefula, Chrisopher Fegan, Adam Finn, Duncan Fullerton, Sanjeev Garg, Sanjeev Garg, Atul Garg, Effrossyni Gkrania-Klotsas, Jo Godden, Arthur Goldsmith, Clive Graham, Elaine Hardy, Stuart Hartshorn, Daniel Harvey, Peter Havalda, Daniel B Hawcutt, Maria Hobrok, Luke Hodgson, Anita Holme, Anil Hormis, Michael Jacobs, Susan Jain, Paul Jennings, Agilan Kaliappan, Vidya Kasipandian, Stephen Kegg, Michael Kelsey, Jason Kendall, Caroline Kerrison, Ian Kerslake, Oliver Koch, Gouri Koduri, George Koshy, Shondipon Laha, Steven Laird, Susan Larkin, Tamas Leiner, Patrick Lillie, James Limb, Vanessa Linnett, Jeff Little, Michael MacMahon, Emily MacNaughton, Ravish Mankregod, Huw Masson, Elijah Matovu, Katherine McCullough, Ruth McEwen, Manjula Meda, Gary Mills, Jane Minton, Mariyam Mirfenderesky, Kavya Mohandas, Quen Mok, James Moon, Elinoor Moore, Patrick Morgan, Craig Morris, Katherine Mortimore, Samuel Moses, Mbiye Mpenge, Rohinton Mulla, Michael Murphy, Megan Nagel, Thapas Nagarajan, Mark Nelson, Igor Otahal, Mark Pais, Selva Panchatsharam, Hassan Paraiso, Brij Patel, Justin Pepperell, Mark Peters, Mandeep Phull, Stefania Pintus, Jagtur Singh Pooni, Frank Post, David Price, Rachel Prout, Nikolas Rae, Henrik Reschreiter, Tim Reynolds, Neil Richardson, Mark Roberts, Devender Roberts, Alistair Rose, Guy Rousseau, Brendan Ryan, Taranprit Saluja, Aarti Shah, Prad Shanmuga, Anil Sharma, Anna Shawcross, Jeremy Sizer, Manu Shankar-Hari, Richard Smith, Catherine Snelson, Nick Spittle, Nikki Staines, Tom Stambach, Richard Stewart, Pradeep Subudhi, Tamas Szakmany, Kate Tatham, Jo Thomas, Chris Thompson, Robert Thompson, Ascanio Tridente, Darell Tupper-Carey, Mary Twagira, Andrew Ustianowski, Nick Vallotton, Lisa Vincent-Smith, Shico Visuvanathan, Alan Vuylsteke, Sam Waddy, Rachel Wake, Andrew Walden, Ingeborg Welters, Tony Whitehouse, Paul Whittaker, Ashley Whittington, Meme Wijesinghe, Martin Williams, Lawrence Wilson, Sarah Wilson, Stephen Winchester, Martin Wiselka, Adam Wolverson, Daniel G Wooton, Andrew Workman, Bryan Yates, and Peter Young.
Funding
This work is supported by grants from: the National Institute for Health Research [award CO-CIN-01], the Medical Research Council [grant MC_PC_19059] and by; a National Institute for Health Research (NIHR) Health Protection Research Senior Investigator award to PO; the Imperial Biomedical Research Centre (NIHR Imperial BRC, grant P45058) and the Health Protection Research Unit (HPRU) in Respiratory Infections at Imperial College London and NIHR HPRU in Emerging and Zoonotic Infections at University of Liverpool, both in partnership with Public Health England, [NIHR award 200907], Wellcome Trust and Department for International Development [215091/Z/18/Z], and the Bill and Melinda Gates Foundation [OPP1209135], Liverpool Experimental Cancer Medicine Centre for providing infrastructure support (Grant Reference: C18616/A25153); a Wellcome Trust-University of Edinburgh Institutional Strategic Support Fund award (to YJC and DH) and the UK Coronavirus Immunology Consortium (UK-CIC). LT is supported by the Wellcome Trust (grant number 205228/Z/16/Z). The views expressed are those of the authors and not necessarily those of the NIHR, BMGF, MRC, Wellcome Trust or PHE.
Acknowledgements
This work uses data provided by patients and collected by the NHS as part of their care and support #DataSavesLives. We are extremely grateful to the 2,648 frontline NHS clinical and research staff and volunteer medical students, who collected this data in challenging circumstances; and the generosity of the participants and their families for their individual contributions in these difficult times.
We also acknowledge the support of Jeremy J Farrar, Nahoko Shindo, Devika Dixit, Nipunie Rajapakse, Lyndsey Castle, Martha Buckley, Debbie Malden, Katherine Newell, Kwame O’Neill, Emmanuelle Denis, Claire Petersen, Scott Mullaney, Sue MacFarlane, Nicole Maziere, Julien Martinez, Oslem Dincarslan, Annette Lake, Lucy Cook, Robert Hammond, Rachael Quinlan, Benjamin Harris, Mohamed Zuhair, Michael Fertleman, Claudia Selck, Danai Koftori, Alex Cocker, Zainab Saeed, Isaac Day-Weber, Aris Aristodemou, Patricia Watber, and Yunchuan Ding.