key: cord-1017431-w5ui1ona authors: de Bruin, Sanne; Bos, Lieuwe D.; van Roon, Marian A.; Tuip-de Boer, Anita M.; Schuurman, Alex R.; Koel-Simmelinck, Marleen J.A.; Bogaard, Harm Jan; Tuinman, Pieter Roel; van Agtmael, Michiel A.; Hamann, Jörg; Teunissen, Charlotte E.; Wiersinga, W. Joost; (Koos) Zwinderman, A.H.; Brouwer, Matthijs C.; van de Beek, Diederik; Vlaar, Alexander P.J. title: Clinical features and prognostic factors in Covid-19: A prospective cohort study date: 2021-05-14 journal: EBioMedicine DOI: 10.1016/j.ebiom.2021.103378 sha: 6cae89a21725f788aaa5aeb266df7c4a35a749d9 doc_id: 1017431 cord_uid: w5ui1ona BACKGROUND: Mortality rates are high among hospitalized patients with COVID-19, especially in those intubated on the ICU. Insight in pathways associated with unfavourable outcome may lead to new treatment strategies. METHODS: We performed a prospective cohort study of patients with COVID-19 admitted to general ward or ICU who underwent serial blood sampling. To provide insight in the pathways involved in disease progression, associations were estimated between outcome risk and serial measurements of 64 biomarkers in potential important pathways of COVID-19 infection (inflammation, tissue damage, complement system, coagulation and fibrinolysis) using joint models combining Cox regression and linear mixed-effects models. For patients admitted to the general ward, the primary outcome was admission to the ICU or mortality (unfavourable outcome). For patients admitted to the ICU, the primary outcome was 12-week mortality. FINDINGS: A total of 219 patients were included: 136 (62%) on the ward and 119 patients (54%) on the ICU; 36 patients (26%) were included in both cohorts because they were transferred from general ward to ICU. On the general ward, 54 of 136 patients (40%) had an unfavourable outcome and 31 (23%) patients died. On the ICU, 54 out of 119 patients (45%) died. Unfavourable outcome on the general ward was associated with changes in concentrations of IL-6, IL-8, IL-10, soluble Receptor for Advanced Glycation End Products (sRAGE), vascular cell adhesion molecule 1 (VCAM-1) and Pentraxin-3. Death on the ICU was associated with changes in IL-6, IL-8, IL-10, sRAGE, VCAM-1, Pentraxin-3, urokinase-type plasminogen activator receptor, IL-1-receptor antagonist, CD14, procalcitonin, tumor necrosis factor alfa, tissue factor, complement component 5a, Growth arrest–specific 6, angiopoietin 2, and lactoferrin. Pathway analysis showed that unfavourable outcome on the ward was mainly driven by chemotaxis and interleukin production, whereas death on ICU was associated with a variety of pathways including chemotaxis, cell-cell adhesion, innate host response mechanisms, including the complement system, viral life cycle regulation, angiogenesis, wound healing and response to corticosteroids. INTERPRETATION: Clinical deterioration in patients with severe COVID-19 involves multiple pathways, including chemotaxis and interleukin production, but also endothelial dysfunction, the complement system, and immunothrombosis. Prognostic markers showed considerable overlap between general ward and ICU patients, but we identified distinct differences between groups that should be considered in the development and timing of interventional therapies in COVID-19. FUNDING: Amsterdam UMC, Amsterdam UMC Corona Fund, and Dr. C.J. Vaillant Fonds. Patients with severe COVID-19 are characterized by hyper-inflammation and coagulopathy. Studies have evaluated plasma biomarkers for inflammation and coagulopathy as prognostic factors for outcome in COVID-19. Many of these studies were retrospective, included a limited number of markers, or were limited to markers measured on admission only. We searched PubMed, EMBASE and Cochrane Reviews for research articles from January 1 2020 to February 1, 2021. More than 200 observational cohorts were published describing prognostic factors for severity of COVID-19. Only few of these studies reported dynamic changes of a very limited set of specific biomarkers; changes in interleukin (IL)-6 and D-dimer were related to clinical outcome. In this large prospective cohort of mild and severe hospitalized COVID-19 patients, serial prognostic biological factors were determined including blood molecular markers of endothelium activation, inflammation, neutrophil activation and neutrophil extracellular traps formation, complement activation, coagulopathy and epithelial barrier disruption. This study shows that endothelial dysfunction is a key characteristic in COVID-19 patients admitted to the ward. Patients with mild COVID-19 progressing to severe COVID-19 needing ICU admission show involvement of inflammation, coagulation, complement and epithelial barrier disruption on top of endothelial dysfunction. This study shows the key pathways involved in the progression of severe COVID-19 and may guide the way to interventional trials in severe COVID-19. Temporal differences between general ward and ICU patients should be taken account when designing interventional studies for severe COVID-19. On the general ward early interventions should focus on enhancing endothelial integrity and limiting chemotaxis. For intubated patients on the ICU, therapeutic interventions in multiple pathways may be needed. However, for single interventions those treatment with a day function within the crossroads of inflammation and coagulation might be most promising. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes COVID-19, a respiratory disease with high clinical variability [1] . Clinical risk factors for developing life-threatening disease are male gender, older age and comorbidities such as diabetes mellitus and hypertension [1] . Host genetic variants of interferon immunity and presence of auto-antibodies against type I interferon also have been identified as risk factors [2] [3] [4] . Patients hospitalized with COVID-19 predominantly present with hypoxemia caused by virus-induced lung inflammation characterized by lymphocyte infiltration and activation of the coagulation system [5] . As a result, 20-40% of the hospitalized COVID-19 patients require ICU admission, [6] and 35-50% of these patients have a fatal outcome [7] . The wide range between rate of ICU admission and death can be explained by various reasons, most importantly differences in the threshold for hospital admission. Biomarkers, mainly blood chemistry markers, [8] have been associated with outcome in COVID-19. However, many of these studies were retrospective, evaluated a limited number of markers, or only evaluated markers measured on admission [8] . So far, dynamic changes of a few biomarkers, including interleukin (IL)-6 and D-dimer, have been associated with outcome in large datasets [ 9 , 10 ] . While it has been postulated that secondary damage in COVID-19 is caused by a "cytokine storm", systemic levels of proinflammatory cytokines are generally lower than reported in patients with acute respiratory distress syndrome (ARDS) from aetiologies other than corona virus [ 11 , 12 ] . Autopsy reports in COVID-19 suggested that respiratory failure is also driven by other pathways than the "cytokine-storm-pathway", for example, endothelial activation, neutrophil activation and neutrophil extracellular trap (NET) formation, immunothrombosis, and epithelial barrier disruption [ 13 , 14 ] . So far, adjunctive dexamethasone therapy has been shown to prevent mortality in severe COVID-19, while monoclonal antibodies against IL-6 showed ambiguous results in clinical trials [ 15-19 , 20 ] . One study preprint reported that monoclonal antibodies against IL-6 improved survival and other clinical outcomes in hospitalised COVID-19 patients with hypoxia and systemic inflammation [21] . Among those not receiving invasive mechanical ventilation at baseline, patients allocated tocilizumab were less likely to reach the composite endpoint of invasive mechanical ventilation or death [21] . Optimal patient selection and timing of treatment must be considered crucial to optimize adjunctive treatments for COVID-19. We performed a prospective cohort study with serial sampling to evaluated the importance of multiple biomarkers in potential pathways over time in severe COVID-19. The Amsterdam UMC COVID-19 biobank is a prospective cohort study containing clinical data and archive material from adult patients admitted with COVID-19 in two academic hospitals of Amsterdam UMC. Based on previous studies, 64 blood molecular markers of endothelial activation, inflammation, neutrophil activation and NET formation, activation of the complement, coagulopathy and epithelial barrier disruption, were measured in serial blood archive samples from patients admitted from 23 rd of March until 26 th of May 2020. Patients or their legal representatives received written information about the study and were asked to give written informed consent for participation. If direct informed consent of patients was not feasible, patients could be included with a deferred consent procedure. In case of deferred consent, patients or their legal representatives were informed about the biobank as soon as possible. To ensure all patients wilfully participated in the biobank and provide a possibility to opt out from the biobank comprehensive information and an opt-out form was send to the patients three months after discharge. This study was approved by the biobank ethics committees of both Amsterdam UMC hospitals (2020_065). Patients admitted to one of the two academic hospitals within Amsterdam UMC were prospectively included if they were admitted on the COVID-19 ward (general ward cohort), COVID-19 ICU (ICU cohort) or on the regular ICU (control ICU cohort). A definite COVID-19 case was defined as a positive PCR from nasopharyngeal swab with clinical signs consistent with COVID-19. For the current analysis, patients were selected when a blood sample was available within 2 days of ward-and/or ICU-admission. So, patients were excluded if there was no blood sample available. Patients were treated according to the local protocol that included thromboprophylaxis, but did not routinely include the use of remdesivir, hydrochloroquine, azithromycin, convalescent plasma, corticosteroids or any other immunomodulatory therapy. High flow nasal cannulae were not used during the inclusion period due to initial concerns for the safety of healthcare workers. Invasively ventilated patients received low tidal volume ventilation with PEEP set to the lower PEEP-FiO2 table and mandatory prone positioning at a PaO2/FiO2 below 150 mmHg. Clinical data were collected using the case record form (CRF) of the World Health Organisation (WHO) [22] . The day of admission to the general ward or to the ICU was identified as timepoint 0. Subsequent timepoints were defined up to 28 days after admission (so day 0-2, followed by day 3-4, day 5-6, day 9-10, day 13-15, day 20-22, day [27] [28] [29] . Biomarkers were measured using a luminex platform or ELISA (online appendix e Table 1 ). Plate-to-plate variation was accounted for using negative and positive controls. Values below the detection limit were imputed with the lower limit of quantification given by the calibration curve for the univariate comparisons, but were not taken into account in joint model analysis (see statistical analysis paragraph). More than 50% of the ferritin, endothelin-1, interleukin (IL)-1 beta and IL-12-p70 measurements were judged to be unreliable because of stringent quality criteria (more than 25 beads counted and no extrapolation outside of the reference standard concentrations) and were therefore excluded for analysis. The analysis was stratified for patients admitted to the general ward and to the ICU. For patients admitted to the ward, the primary outcome was admission to the ICU or mortality (referred to as unfavourable outcome). For patients admitted to the ICU, the primary outcome was 12-week mortality. Kaplan-Meier curves for the risk of unfavourable outcomes are reported in the supplement material. Secondary outcomes were severe complications of COVID-19 infection: pulmonary embolism, vascular thrombotic events and acute kidney failure defined by the need for renal replacement treatment (only for ICU patients). Patients who clinically deteriorated routinely underwent pulmonary CT angiography to rule out pulmonary embolism. COVID-19 patients were screened for vascular thrombotic events with duplex ultrasound once weekly. All analyses were performed in R-statistics through the R-studio interface. Continuous data was summarized with mean and standard deviation or median with inter-quartile range (IQR) dependent on the distribution of values. Categorical data were summarized with numbers and percentages. Differences between groups at baseline were tested with a T-test, Mann-Whitney U or Fisher exact test where appropriate. The associations between outcome risk and biomarker-values were estimated using joint models that combine Cox regression and linear mixed-effects (LME) models [23] . With the LME part of the joint models we estimated the linear change pattern of Table 1 Baseline demographics. the biomarkers over follow-up time. We used the assumption that both the intercepts and the slopes of follow-up time varied between patients if indicated based on the Akaike Information Criterion (AIC). Absolute biomarker concentrations and fold changes in those concentrations over time were illustrated with boxplots, stratified per cohort and outcome. The log Hazard Ratios (HRs) of the time-dependent biomarker-values for outcome were estimated with the Cox regression part of the joint models and the Weibullspecification was used to describe the cumulative baseline hazard curve [23] . These analyses were done in the ward and ICU cohort separately for the 60 biomarkers with data of sufficient quality. For three biomarkers in the ward cohort (complement component 3a, C-reactive protein and Plasminogen activator inhibitor-1) the joint models did not converge and a Cox regression model with timedependent co-variates was used instead. To evaluate the baseline risk for outcomes, HRs of baseline biomarker values and of clinical risk factors at day of admission were estimated using Cox regression. All HRs were adjusted for age and gender of the patients. In order to evaluate the association between biomarker concentration at baseline and the occurrence of severe complications (pulmonary embolism, thrombotic vascular events and kidney failure) multinomial regression was performed and the odds ratios (ORs) were reported for occurrence of a complication and mortality. Missing baseline data was imputed using the linear mixed effect model for this analysis. All biomarkers were [2] log-transformed in order to normalize their distributions. As a consequence, the reported Hazard Ratios must be interpreted as effects on outcome-hazard upon a doubling of the biomarker. Statistical significance was defined after Bonferroni-correction for the number of analysed biomarkers (pvalue less than 0 • 05 divided by 60). For pathway analysis, the biomarkers were matched to their corresponding genes via the NCBI gene library. Biomarkers that were significantly associated with outcome in the joint model analyses were used for gene ontology enrichment analysis ( en-richGO algorithm from the clusterProfiler package), [24] using list of all measured biomarkers mapped to the human genome as reference. Significant pathways were defined by a P-value below 0 • 01 after Bonferroni correction with a Q-value below 0 • 05. This study was funded by the Amsterdam UMC, Amsterdam UMC Corona Fund, and Dr. C.J. Vaillant Fonds (to DB). Funders had no role in the in study design, data collection, data analyses, interpretation, or writing of this report. Between March 23 rd and May 26 th 2020, 416 patients with COVID-19 were included in the Amsterdam UMC COVID-19 Biobank (online appendix e Fig. 1 ). From 200 of the 417 patients (48%), no blood sample was present on the first day, resulting in 219 (52%) patients eligible for the current analysis. The 219 selected patients were were older and more frequently male and had a more severe clinical course defined by increased thromboembolic events and mortality (online appendix eTable 2). Of the 219 patients included, 136 (62%) were initially managed on the ward and 119 patients (54%) were admitted to the ICU. 36 patients (26%) were transferred from the ward to the ICU and were included in both cohorts. Median age of general ward COVID-19 patients was 64 years (SD 12 • 1; Table 1 ) of whom 86 out of 136 were men (63%). Median age of ICU COVID-19 patients was 62 years (SD 10 • 6) of whom 84 out of 119 were men (71%). 25 For patients on the general ward, we used 26 age-matched controls from the outpatient clinic, with a mean age of 64 years (SD 15 • 5; online appendix eTable 3) of whom 18 (69%) were male. For ICU COVID-19 patients, controls were patients admitted on the ICU who did not have COVID-19 disease or other viral respiratory illness, with a median age of 60 years (SD 17 • 6; online appendix eTable 4) of whom 17 were men (68%). Most common admission diagnoses were respiratory failure (6 of 25 [24%]) and post-surgical care (7 of 25 [28%]). Several plasma biomarkers were associated with outcome ( Figs. 1 and 2 ; online appendix eTable 3 and eFig. 4-9). On the general ward, changes in concentrations over time (comparison between the first blood with subsequent blood samples) of IL-6, IL-8, IL-10, soluble Receptor for Advanced Glycation End Products (sRAGE), vascular cell adhesion molecule 1 (VCAM-1) and Pentraxin-3 (PTX3) were associated with unfavourable outcome. On the ICU, changes in concentrations over time (first blood in ICU and then subsequent blood samples) of IL-6, IL-8, IL-10, sRAGE, VCAM-1, Pentraxin-3, urokinase-type plasminogen activator receptor (UPAR), IL-1-receptor antagonist (IL-1RA), CD14, procalcitonine, tumor necrosis factor alfa (TNFa), tissue factor (TF), complement component 5a (C5a), Growth arrest-specific 6 (GAS-6), angiopoietin 2 (ANG2), and lactoferrin were associated with death. Fig. 3 shows the individual biomarker trajectories of four individual patients and highlights the change in prognosis based on the dynamic biomarker profiles. The association between each biomarker and outcome at baseline alone is reported in the online appendix eFig. 10. Temporal changes in the plasma concentration of CRP and PAI-1 showed a stronger association with adverse outcome in patients admitted to the normal ward when corrected for BMI in a sensitivity analysis (online appendix eFig. 11). The same analysis for patients admitted to the ICU showed a stronger association for TFF3 but a weaker association for MPO (online appendix eFig. 12). Baseline biomarker concentrations were associated with pulmonary embolism, thrombotic vascular events and kidney failure, and with unfavourable outcome on the general ward and death on the ICU (online appendix eFig. 13, 14, and 15). The biomarkers showed a similar prognostic value for pulmonary embolism and thrombotic vascular events as for death. However, the association between kidney failure and baseline concentrations for Cystacin C (OR 3 • 89, 95%-CI: 1 • 80-8 • 39), Neutrophil gelatinaseassociated lipocalin (NGAL; OR 2 • 81, 95%-CI: 1 • 53-5 • 20), Trefoil factor 3 (TFF3; OR 12 • 2, 95%-CI: 3 • 93-38 • 0), TNF-receptor inhibitor (TNFRI; OR 6 • 51, 95%-CI 2 • 81-15 • 1) and Triggering receptors expressed on myeloid cells-1 (TREM-1; OR 5 • 52, 95%-CI: 2 • 46-12 • 4) was stronger than the associations between these markers and death. Pathway analyses revealed that on the general ward biomarkers associated with outcome were mainly mapped to chemotaxis and interleukin production (Fig. 4, online appendix eTable 4) . For the ICU patients, there was a much wider variety of pathways in- volved which included chemotaxis, cell-cell adhesion, innate host response mechanisms, including the complement system, viral life cycle regulation, angiogenesis, wound healing and response to corticosteroids (online appendix eTable 5) . In our prospective study with serial sampling, we evaluated the combination of multiple biomarkers in potential pathways involved in COVID-19. Our results show that deterioration in patients with hospitalized COVID-19 involves multiple pathways, including chemotaxis and interleukin production, but also endothelial dysfunction, the complement system, and immunothrombosis. Identi- Fig. 3 . Schematic overview of the pathophysiological processes involved in clinical deterioration in ward (left) and ICU patients with Covid-19 (right). Dynamic changes of important biomarkers in four patients that examplify the change in prognosis based on longitudinal biomarker information modeled by joint model analysis. legend: Four plots illustrating the change of four markers since day of ICU admission in four different patients. Panel A shows the seven available IL-6 blood levels measurements for the first patient (male, 63 yr). At day of admission IL-6 level was close to the average baseline IL-6 levels in our sample and therefore the predicted survival-curve at ICU admission for this patient (green line) was close to the observed Kaplan-Meier curve for the entire sample of patients. At days 3 and 5 of ICU admission IL-6 levels were increased compared to baseline and the predicted survival-curves at these days (dark and light blue curves) indicated substantial worsening of the prognosis. At days 9, 13, 20 and 27 IL-6 levels were decreased in this patient and the predicted survival-curves at these days (purple, yellow, grey and black curves) indicated clear improvement of the prognosis for this patient who was discharged alive from the ICU at day 38. Panels B and C illustrate similar change patterns of the C5a and RAGE markers in two other patients (male, 59 yr and male 74 yr), who were discharged alive from the ICU alive after 48 and 44 days, respectively. Panel D illustrates the change pattern of six measurements of the VCAM-1 blood levels in a 74 yr old male patient. In contrast to what was observed in other patients, the VCAM-1 measurements of this patient remained relatively stable at an increased level during his ICU stay, and a progressive worsening of the predicted survival-curves was seen (green, blue, purple, yellow, black curves). The patient died at day 47 after ICU admission. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.) fied prognostic markers showed considerable overlap between patients on the ward and ICU, but also distinct differences. On the general ward the strongest predictors were indicative of endothelial activation and chemotaxis. On the ICU, enhanced involvement of inflammation, activation of the complement system and coagulation, and further breaching of the epithelial barrier were found as markers of poor prognosis. Evaluated markers may not indicate a direct causal relationship, but our findings have implications for development and timing of interventional therapies in severe COVID-19. Our findings implicate that therapeutic interventions in patients with COVID-19 on the general ward, so early after hospital admission, should aim to improve endothelial integrity and limit chemotaxis towards the alveolar compartment. Concentrations of IL-6 were associated with unfavourable outcome for ward and ICU patients. Tocilizumab is a humanized IL-6 receptor-inhibiting monoclonal antibody for the rheumatoid arthritis giant cell arteritis, systemic juvenile idiopathic arthritis and for people with chimeric antigen receptor T cell-induced severe or life-threatening cytokine release syndrome. Ten randomized controlled studies on anti-IL6 treatment in hospitalized patients with COVID-19 have been reported [ 17-21 , 25-28 ] . These studies showed conflicting data regarding the use of anti-IL6 for COVID-19 patients, which may be explained by the additional use of dexamethasone, the timing of start of anti-IL6 in the disease course and patient selection. Our results may help to optimize the timing of anti-IL6 as well as the selection of patients which may benefit most. Imatinib is one of the drugs that reinforce the endothelial barrier and mitigate alveolar inflammatory responses through NFkB mediated chemotaxis in several models of acute lung injury and may therefore prevent clinical progression of COVID-19 [29] . On an individual level, stable VCAM-1 concentration over time were strongly associated with progressive worsening of the predicted survival-curves. At least two randomized controlled studies are ongoing that test this intervention in hospitalised patients with COVID-19 (EUDRACT2020-005447-23) [30] . Other drugs that could limit chemotaxis and endothelial dysfunction are blockers of RAGE and transient receptor potential vanilloid 4 (TRPV4) channel inhibitors [ 31 , 32 ] . These drugs are currently not registered as being tested in randomized controlled studies for COVID-19. Our study confirms the key role of hyperinflammation and immunothrombosis in severe COVID-19 on the ICU. Some patients are at risk to develop severe COVID-19, which can partly be explained by presence of the risk factors male gender, older age and comorbidities. Recent studies have provided insight why patients with severe COVID-19 develop hyperinflammation, among others due to improper IgG glycosylation, auto-antibodies against IFN and Annexin A2, and genetic polymorphisms [ 2 , 3 , 33 , 34 ] . The occurrence of hyperinflammation, illustrated by enhancing levels of IL-6, IL-8, and sRAGE, explains the biphasic pattern in COVID-19, with an early viral replication phase, followed by a hyper-inflammatory phase involving cytokine release [35] . Identified markers in our study go beyond cytokine release, with PTX3, GAS6, and C5a. PTX3 is an essential component of humoral innate immunity, involved in resistance to selected pathogens and in the regulation of inflammation [36] . RNA-sequencing analysis of peripheral blood mononuclear cells, single-cell bioinformatics analysis and immunohistochemistry of lung autopsy samples revealed that myelomonocytic cells and endothelial cells express high levels of PTX3 in patients with COVID-19 [36] . The cumulative response of the immune system to SARS-CoV-2, both through inflammation and immune cell expression of prothrombotic proteins, is likely to be a major contributor to hypercoagulability in COVID-19. Recognition of pathogen-associated molecular patterns through the toll like receptors and CD14 receptor of monocytes promotes the transcription and expression of TF [37] . Neutrophils of COVID-19 patients yielded high TF expression and released NETs carrying active TF [38] . Treatment of control neutrophils with COVID-19 platelet-rich plasma generated TF-bearing NETs that induced thrombotic activity of endothelial cells [38] . Thrombin or NETosis inhibition or C5a blockade attenuated platelet-mediated NET-driven thrombogenicity. The association of COVID-19 inflammation with activation of the C5a-C5aR1 axis has been reported [39] . The potent anaphylatoxin C5a attracts neutrophils and monocytes to the infection site, and strongly activates these cells, causing tissue damage by oxidative radical formation and enzyme release but also inducing release of tissue factor from endothelial cells and neutrophils thereby activating the coagulation system [40] . A phase 2 randomized controlled trial show that C5a inhibition was safe and associated with decreased risk for pulmonary embolisms in severe COVID-19 patients [40] . The GAS6-PROS1/TAM system has been suggested to play an important role in SARS-CoV-2 infection and progression complications [41] . Dysregulation of the urokinase receptor system (UPAR) have been associated with the development of immunothrombosis associated with respiratory failure in COVID-19 patients [45] . Our study has limitations. First, we studied hospitalized patients with COVID-19. Patients with mild complaints were not taken into consideration in this analysis. This limits our conclusions to patients with severe COVID-19 disease admitted to the hospital. Second, we evaluated selected patients from the Amsterdam UMC Covid-19 Biobank. Selection was made on the basis of plasma sample availability on the first admission day. Although selected patients were comparable to non-selected patients, they tended to have more comorbidities and a more severe clinical course. This could well be explained by inclusion of relatively higher proportion of severe patients, who had a relatively high frequency of blood withdrawals. Third, we included patients from two academic hospitals in the Netherlands. Nevertheless, we believe the cohort is representative of and for the wider population. Of course, in a "normal" academic setting a potential bias maybe introduced this way. For example, transfers often have missing data and they will be excluded this way. In a normal setting transfers are often representative for complex academic cases. In the setting of COVID-19 care no transfers were made for academic expertise. The only reason for inter-hospital transfers in the Netherlands (and the same holds true for other European countries) was capacity issues. So, we believe the cohort is representative of and for the wider population. Fourth, we did not evaluate ethnicity in the current study, so we could not include this in the prognostic model. Fifth, identified markers may not indicate a causal relationship. Higher levels of IL-10 were associated with unfavourable out-come. IL-10 is considered anti-inflammatory, but has been noted to be higher in sepsis patients, proposed the increase might be an attempt to moderate the immune response. Sixth, our patients were included in the first wave of COVID-19 in the Netherlands. Ever since treatment has involved, including the use of anticoagulants and adjunctive dexamethasone therapy. The impact on outcome of the evolving treatment of COVID19 has been limited, but nevertheless, the introduction of steroids might have impacted prognostic factors. Detailed data on treatment is being collected but currently not available for the analysis. Finally, some potential important biomarkers were not included. For example, we did not evaluate troponin that has been associated with complications in COVID-19 patients [ 42 , 43 , 44 ] . We aimed to provide insight in pathways associated with unfavourable outcome, and did evaluate the additive prognostic effect of biomarkers on known clinical prognostic factors in COVID-19. Our results may explain why single pathway intervention studies in severe COVID-19 so far remain negative or have limited impact on outcome, and only a general broad intervention on inflammation such as steroids showed benefit [ 15-20 , 24-27 ] . We did not evaluate which markers were the best independent biomarkers predicting unfavorable outcome. However, we identified multiple pathways to be important in the pathophysiology of severe COVID-19. On the general ward early interventions should focus on enhancing endothelial integrity and limiting chemotaxis. For intubated patients on the ICU, therapeutic interventions in multiple pathways may be needed. However, for single interventions those treatment with a day function within the crossroads of inflammation and coagulation might be most promising. Sanne de Bruin: study design, data collection, data analysis, data interpretation, writing Lieuwe D. Bos:s tudy design, data collection, data analysis, data interpretation, writing Marian A. van Roon: data collection, data analysis, data interpretation Anita M. All the authors had full access to the data, vouch for the completeness and accuracy of the data. The first draft of the manuscript was written by the first, second and last authors with input from authors. All the authors participated in reviewing and editing the manuscript and approved the submitted draft. Data will be shared on request for non-commercial purposes. To submit a request to the Amsterdam UMC COVID-19 Biobank please contact the corresponding author (d.vandebeek@amsterdamumc.nl). The authors have declared no potential conflicts of interest. Features of 20 133 UK patients in hospital with covid-19 using the ISARIC WHO clinical characterisation protocol: prospective observational cohort study Autoantibodies against type I IFNs in patients with life-threatening COVID-19 Inborn errors of type I IFN immunity in patients with life-threatening COVID-19 Genetic mechanisms of critical illness in Covid-19 Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study Epidemiology, clinical course, and outcomes of critically ill adults with COVID-19 in New York City: a prospective cohort study Ventilation management and clinical outcomes in invasively ventilated patients with COVID-19 (PRoVENT-COVID): a national, multicentre, observational cohort study Prognostic factors for severity and mortality in patients infected with COVID-19: a systematic review D-dimer levels on admission to predict in-hospital mortality in patients with Covid-19 Prognostic value of interleukin-6, C-reactive protein, and procalcitonin in patients with COVID-19 Prevalence of phenotypes of acute respiratory distress syndrome in critically ill patients with COVID-19: a prospective observational study Cytokine elevation in severe and critical COVID-19: a rapid systematic review, meta-analysis, and comparison with other inflammatory syndromes Viral presence and immunopathology in patients with lethal COVID-19: a prospective autopsy cohort study Insights into pathogenesis of fatal COVID-19 pneumonia from histopathology with immunohistochemical and viral RNA studies Effect of hydrocortisone on mortality and organ support in patients with severe COVID-19: the REMAP-CAP COVID-19 corticosteroid domain randomized clinical trial RECOVERY Collaborative Group Dexamethasone in hospitalized patients with Covid-19 -preliminary report Efficacy of tocilizumab in patients hospitalized with COVID-19 Tocilizumab in patients hospitalized with COVID-19 pneumonia Effect of tocilizumab on clinical outcomes at 15 days in patients with severe or critical coronavirus disease 2019: randomised controlled trial REMAP-CAP Investigators Interleukin-6 receptor antagonists in critically ill patients with Covid-19 Tocilizumab in patients admitted to hospital with COVID-19 (RECOVERY): preliminary results of a randomised, controlled, open-label, platform trial. medRxiv Epub World Health Organisation. Global COVID-19 clinical platform: rapid core case report form (CRF) Joint Models for Longitudinal and Time-to-Event Data clusterProfiler: an R package for comparing biological themes among gene clusters Effect of tocilizumab vs standard care on clinical worsening in patients hospitalized with COVID-19 pneumonia: a randomized clinical trial Effect of tocilizumab vs usual care in adults hospitalized with COVID-19 and moderate or severe pneumonia: a randomized clinical trial Tocilizumab in hospitalized patients with COVID-19 pneumonia Interleukin-6 blocking agents for treating COVID-19: a living systematic review Emerging pharmacological therapies for ARDS: COVID-19 and beyond Safety and efficacy of imatinib for hospitalized adults with COVID-19: a structured summary of a study protocol for a randomised controlled trial TRPV4 inhibition counteracts edema and inflammation and improves pulmonary function and oxygen saturation in chemically induced acute lung injury Urgent reconsideration of lung edema as a preventable outcome in COVID-19: inhibition of TRPV4 represents a promising and feasible approach Afucosylated IgG characterizes enveloped viral responses and correlates with COVID-19 severity Clinical effectiveness of drugs in hospitalized patients with COVID-19 infection: a systematic review and meta-analysis. medRxiv Life-threatening COVID-19: defective interferons unleash excessive inflammation Macrophage expression and prognostic significance of the long pentraxin PTX3 in COVID-19 Platelet activation and platelet-monocyte aggregate formation trigger tissue factor expression in patients with severe COVID-19 Complement and tissue factor-enriched neutrophil extracellular traps are key drivers in COVID-19 immunothrombosis Association of COVID-19 inflammation with activation of the C5a-C5aR1 axis Anti-C5a antibody IFX-1 (vilobelimab) treatment versus best supportive care for patients with severe COVID-19 (PANAMO): an exploratory, open-label, phase 2 randomised controlled trial Role of Vitamin K-dependent factors protein S and GAS6 and TAM receptors in SARS-CoV-2 infection and COVID-19-associated immunothrombosis Sanchis-Gomar F . Cardiac troponin I in patients with coronavirus disease 2019 (COVID-19): evidence from a meta-analysis. Prog Cardiovasc Characterization of myocardial injury in patients with COVID-19 Cardiac Injury in COVID-19-Echoing Prognostication Soluble urokinase plasminogen activator receptor (suPAR) as an early predictor of severe respiratory failure in patients with COVID-19 pneumonia This study was funded by the Amsterdam UMC, Amsterdam UMC Corona Reseach Fund, and Dr. C.J. Vaillant Fonds (to DvdB). We would like to thank all medical, paramedical, laboratory and nursing staff involved in the care of the COVD-19 patients for making it possible to build the Amsterdam UMC COVID-19 Biobank Cov in the middle of the COVID-19 outbreak in The Netherlands. Supplementary material associated with this article can be found in the online version at doi: 10.1016/j.ebiom.2021.103378 .