Introduction
COVID-19, the illness caused by infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has caused a global pandemic and has placed an unprecedented strain on healthcare systems worldwide.
Early studies from China found that older age, male sex and a number of comorbidities (notably hypertension, cardiovascular disease and respiratory disease) [Reference Cao1–Reference Chen6] were associated with poorer outcomes. Furthermore, a number of blood test abnormalities such as raised inflammatory markers [Reference Zhang3, Reference Giacomelli7–Reference Qiu11] and lymphopoenia [Reference Li10, Reference Huang and Pranata12] have been associated with worse outcomes.
In the UK, two large cohort studies have corroborated the associations between older age, male sex and pre-existing comorbidities with poorer outcomes [Reference Williamson13, Reference Docherty14]. However, these studies included demographic and morbidity characteristics only, with no data to indicate disease severity at the time of presentation to hospital. The inclusion of investigations and observations that indicate clinical condition at presentation allows for increased acuity of predictive models. The aim of this large UK inpatient cohort study is to explore the association between results of admission laboratory tests and clinical observations, alongside demographic and morbidity characteristics, with clinical outcomes of patients hospitalised with COVID-19. All patients admitted with COVID-19 in the study period were enrolled thus avoiding the selection bias observed in other cohorts.
Methods
Study population
Northwick Park Hospital is a 658 bedded district general hospital serving an ethnically and socio-economically diverse population in North West London [15]. This observational cohort study retrospectively identified all adult patients (age ≥18 years) admitted to this hospital between 12 March and 15 April 2020 who were confirmed to be positive for SARS-CoV-2 by reverse transcription-polymerase chain reaction of pooled nasal and pharyngeal swabs. Patients were tested for COVID-19 in accordance with contemporaneous Public Health England guidelines [16]. Advice on which patients should be tested remained constant throughout the study period. Patients were admitted to hospital at the discretion of the admitting physician and managed using best supportive treatment, including high flow oxygen, antibiotics at the discretion of the treating physician, continuous positive airway pressure (CPAP), invasive ventilation and organ support as needed. This study had National Research Ethics approval (REC 20/NW/0218, IRAS 282630) and was carried out in accordance with the Declaration of Helsinki and the principles of Good Clinical Practice.
Variables
Data were collected from electronic patient records and entered into a securely stored database. This was pseudonymised for data analysis. Age, sex and ethnicity were recorded. Ethnicity was categorised as White, Black, Asian or Other (including mixed ethnicity). Comorbidities and medications which previous research had suggested might be associated with prognosis were extracted, including: hypertension, heart failure (HF), ischaemic heart disease (IHD), diabetes mellitus (DM), active malignancy and respiratory disease (including asthma, chronic obstructive pulmonary disease (COPD), pulmonary fibrosis and obstructive sleep apnoea). Use of angiotensin converting enzyme inhibitors (ACEi), angiotensin-receptor blockers (ARBs), metformin and immunosuppressive medications (including for autoimmune disease and anti-rejection drugs for solid organ transplant) were recorded. Information on comorbidities from admission documentation was cross-referenced with recent inpatient and specialist communications.
The following observations were collected on arrival to hospital: heart rate, blood pressure, temperature, respiratory rate, Glasgow Coma Scale (GCS) and hypoxia. Hypoxia was defined as oxygen saturations of <94% on air or the need for supplemental oxygen to maintain saturations ≥94%. Exceptions were made for two patients with established COPD who were known to retain carbon dioxide and had pre-COVID-19 saturations <94%.
The following blood test results were collected: haemoglobin, platelet count, lymphocyte count, neutrophil count, potassium, sodium, creatinine, urea, alkaline phosphatase (ALP), alanine aminotransferase (ALT), bilirubin, albumin, C-reactive protein (CRP) and lactate. As lactate levels may change rapidly in response to clinical interventions lactate was only recorded if taken within 4 h of admission; all other blood test results represent the first recorded on admission.
The presence or absence of formally reported pulmonary infiltrates on plain chest radiographs taken on admission was recorded.
Medications potentially active against SARS-CoV-2 as part of research trials including remdesivir, dexamethasone, hydroxychloroquine, lopinavir-ritonavir, azithromycin and tocilizumab were recorded.
Outcomes
Outcomes were ascertained for all patients on the census date of 13 May 2020, 28 days following the latest admission date. For all patients who were discharged from hospital alive, the end of follow-up date was considered to be the census date.
The primary outcome was death during the study period. Two secondary outcomes were considered. First, the composite of death or invasive mechanical ventilation (IMV). Second, the composite of death, IMV or CPAP.
Statistical analysis
Descriptive statistical analyses were performed on all of the population characteristics. Kaplan–Meier methods were used to calculate and graph cumulative survival.
Analyses examining factors associated with survival times were performed using a Cox proportional hazards model. Variable selection was executed by a multi-stage procedure. First, independent variables were examined in a series of univariable analyses. Where necessary, higher-order terms were used to better describe the relationship between continuous variables. Continuous variables found to have a heavily positively skewed distribution were analysed on the log10 scale.
In the second stage, only factors found to show association with survival times were included. A backwards selection procedure was used to retain the significant variables in the final model.
To visually display the effect of the combinations of risk factors on outcomes, patients were allocated to three equal-sized groups based on their predicted risk from the multivariable regression model. These outcomes were displayed as tertiles on a Kaplan–Meier plot.
Forty-one patients did not have a positive swab result until >5 days after admission. We were unable to accurately establish the reason for this, but potential explanations include initial false-negative swabs and presentations that did not initially fit the contemporaneous case definition of COVID-19. To explore this further, sensitivity analyses were performed with these patients excluded.
Statistical analysis was performed using STATA v15.1 (StataCorp LLC, Texas, USA).
Results
Demographics and characteristics
Nine-hundred-and-eighty-one patients were included in the study. The median age was 69 (interquartile range 56–80) years, 64.3% were male and 776 (79.1%) had at least one documented comorbidity. The cohort's full demographic and comorbidity data are shown in Table 1.
Demographics, comorbidities and medication usage
ACEi/ARB, angiotensin converting enzyme inhibitor/angiotensin II receptor blockers; COPD, chronic obstructive pulmonary disease; IQR, interquartile range; OSA, obstructive sleep apnoea.
Values are presented as number (percentage) unless otherwise stated.
Thirty patients received trial drugs for COVID-19: one received azithromycin alone; one received azithromycin with tocilizumab; 10 received dexamethasone; 10 received hydroxychloroquine; four received lopinavir-ritonavir and four received remdesivir.
Laboratory findings
A raised CRP was observed in 96.2%, and 60% had raised platelets, consistent with inflammation. A significant number of patients had hypernatraemia (58.0%), raised urea levels (40.7%) and lymphopoenia (76%).
Almost all patients (97.9%) had a chest radiograph taken on admission and of these 75.0% were reported as having pulmonary infiltrates. (The blood test results are available in the Supplementary material.)
Cohort outcomes
Median follow-up was 37 days (range 10–46). By the censoring date, 354 (36.0%) had died. A Kaplan–Meier survival curve is shown in Figure 1. Of the 627 patients alive at the end of the study period, 566 (90.3%) had been discharged from hospital, 13 (2.1%) had been discharged to rehabilitation facilities and 48 (7.7%) remained on acute wards.
Kaplan–Meier plot of time to death (reduced y-axis scale). Blue shading corresponds to 95% confidence intervals.
For the 579 discharged alive the median duration of stay was 8 days (range 0–46).
Out of 981 patients admitted in the study period, 114 (11.6%) required CPAP and 160 (16.3%) IMV. Notably, 28 patients (24.6%) who required CPAP survived without requiring IMV. The cohort outcomes are shown in Figure 2.
Cohort outcomes. CPAP, continuous positive airway pressure; IMV, invasive mechanical ventilation.
Primary outcome
The association between baseline demographic, clinical, laboratory and radiological features with death, the primary outcome was examined. In multivariable analysis, 12 variables retained statistical significance. The final model was based on data from 741 of the 981 patients, due to missing values for some variables. These results are shown in Table 2.
Primary outcome – time to death
ACEi/ARB, angiotensin converting enzyme inhibitor/angiotensin II receptor blockers.
*Hazard ratios given for a 5-unit increase in variable.
**Hazard ratios given for a 10-unit increase in variable.
***Hazard ratios given for a 50-unit increase in variable.
****Hazard ratios given for a 100-unit increase in variable.
+Variable analysed on the log scale (base 10).
Variables found to be significant in univariable analysis but did not retain significance in multivariable analysis were excluded from the final model.
Age was found to be associated with the risk of death, with a 53% higher chance for each 10-year increase in age (adjusted hazard ratio (aHR) 1.53 [1.37–1.71]). No association was found with either sex or ethnicity. The absence of association between ethnicity and death persisted when Black, Asian and minority ethnic (BAME) groups were combined (see Supplementary material).
Respiratory disease (aHR 1.37 [1.03–1.81]) was associated with a higher risk of death as was immunosuppression (aHR 2.23 [1.23–4.05]). Conversely, no association was found with DM, hypertension, IHD, HF or use of ACEi, ARBs, or metformin.
Of the initial bedside observations, respiratory rate, hypoxia and low GCS retained a statistically significant association with death in multivariable analysis. Each increase in respiratory rate by 10 breaths/min was associated with a 28% increase risk of death (aHR 1.28 [1.08–1.52]) and GCS <15 had a 92% increased risk of death (aHR 1.92 [1.41–2.62]). Hypoxia conferred a 36% increased risk of death (aHR 1.36 [1.03–1.82]).
Higher urea levels were associated with an increased risk of death (aHR 2.67 [1.64–4.32]), as were higher ALP (aHR 2.53 [1.37–4.67]), CRP (aHR 1.15 [1.08–1.22]) and lactate (aHR 2.67 [1.65–4.33]). Higher platelet counts were associated with a reduced risk of death (aHR 0.77 [0.67–0.89]).
The presence of pulmonary infiltrates on chest radiograph was associated with an increased risk of death (aHR 1.89 [1.38–2.60]).
In order to illustrate the difference in outcome between patients with different characteristics, patients were categorised into three tertiles based on the multivariable model. Survival at 60 days was 90% for the patients in the tertile with the ‘best’ combination of risk factors compared to 33% for the patients in the tertile with the ‘worst’ combination (Fig. 3).
Kaplan–Meier plot of time to death split into tertiles by predicted risk from the multivariable model.
A sensitivity analysis was performed excluding patients with a first positive swab >5 days after admission. The same variables were included as in the original model and yielded comparable results (see Supplementary material).
Secondary outcomes
Time to death or IMV was examined as a secondary composite endpoint using the same variables as for the primary outcome. Of the 981 patients in the cohort, 426 died or required IMV during the study period. The variables associated with this secondary outcome were similar to those for the primary endpoint, although ALP, respiratory disease and immunosuppressive medication were not found to be significantly associated. In addition, lower sodium and albumin levels were associated with higher risk of death or requiring IMV.
Four-hundred-and-sixty-one patients incurred the secondary endpoint of death, IMV or CPAP. The same variables were significantly associated with this endpoint as with the time to death or IMV, although the size of effects varied for some factors.
The full description of these secondary outcomes is shown in Table 3. Kaplan–Meier plots and sensitivity analyses are shown in the Supplementary material.
Secondary outcomes – multivariable analyses
CI, confidence interval.
*Hazard ratios given for a 5-unit increase in variable.
**Hazard ratios given for a 10-unit increase in variable.
***Hazard ratios given for a 50-unit increase in variable.
****Hazard ratios given for a 100-unit increase in variable.
Variables found to be significant in univariable analysis but did not retain significance in multivariable analysis were excluded from the final model. Only variables found to be associated with outcomes are shown in this table. The variables considered are the same as shown in Table 2 and the full results are available in the Supplementary material.
Discussion
This large cohort study of patients hospitalised with COVID-19 describes the association between a number of clinical factors and the risk of severe disease. The study was conducted in a major London hospital that serves an ethnically diverse population, with 69.6% of study patients from BAME groups. Overall, the mortality rates were comparable with those in other parts of the UK [Reference Docherty14]. Although the results are from a single hospital, our catchment area covers a part of London with a very wide range of income, social status and ethnicity. We have no reason to believe that it is not representative of the city generally.
This study demonstrates a dramatic difference in outcomes between those who are in the tertile with the best combination of risk factors (90% survival at 60 days) vs. those in the tertile with the worse combination of risk factors (33% survival). This provides clear evidence that risk-stratification can be made early in hospital presentation.
The results of this study generally corroborate a number of previously published findings. However, some differences also exist. Increasing age was strongly associated with poorer outcomes, corroborating the findings of others [Reference Zhang3, Reference Chen6–Reference Petrilli8, Reference Wang17, Reference Luo18]. In line with previous research [Reference Zhou19, Reference Sanchez-Ramirez and Mackey20] a strong association between respiratory diseases and death was found. An association was also found between the use of immunosuppression and death which has been reported separately [Reference Vaid21].
This population had a high incidence of hypertension and although this was a risk factor for severe disease in univariable analysis, it was not associated in multivariable analysis. This implies that although hypertension itself is not a risk factor for severe disease, it is associated with conditions which are. This is in line with other studies [Reference Giacomelli7–Reference Bhargava9].
Unlike the recent UK cohort reported by the OpenSAFELY Collaborative [Reference Williamson13], this study did not find an association between ethnicity and worse outcomes. Although BAME groups have been noted to have poor outcomes from COVID-19 at a population level, the explanation of which is likely to be a complex interaction between socio-economic and medical factors, the outcomes between different ethnic groups in hospitalised cohorts has often been found to be similar [Reference Bhargava9, Reference Gold22–24] and this was the case in our study.
In this cohort, 64.3% were male, which is notably higher than the hospital's catchment area (50%) [15]. In line with a number of cohort studies of hospitalised patients, we failed to demonstrate an association with sex when adjusting for other factors [Reference Chen6, Reference Bhargava9, Reference Cummings25]. One large UK cohort did find sex to be a significant factor when examined at the population level [Reference Docherty14]. It may be that BAME and male patients are prone to more severe disease and therefore more of them appear in hospital-based cohorts but once in hospital they do not have higher mortality. Further research is needed to understand this.
Seventy-five percent of this cohort had pulmonary infiltrates reported on their admission chest radiograph which were associated with poorer outcomes as corroborated by other studies [Reference Liang26–Reference Borghesi28].
Several of the laboratory investigations taken on admission were associated with worse outcomes including CRP, urea and ALP. Although several studies have found an association with raised hepatic enzymes [Reference Garrido, Liberal and Macedo29–Reference Chen31], an independent association between death and ALP has not previously been reported. ALT was not associated with outcomes, suggesting that hepatocellular injury is not the principal explanation for this. Notably, the relationship between outcomes and ALP has been demonstrated in other infectious diseases and this finding merits further research [Reference Hwang32].
Unlike some previous studies, reduced lymphocyte counts were not found to be associated with death [Reference Huang and Pranata12]. An association was found between raised lactate and adverse outcomes. Although a few other studies have included this variable, it has been previously reported [Reference Price-Haywood23, Reference Li33]. Hyperlactataemia may represent hypoperfusion in patients with a sepsis-like presentation or a high adrenergic state secondary to immune activation and it should be explored in future research.
The results of the analyses on secondary endpoints were broadly comparable with the results of the primary endpoints, although neither the use of immunosuppression nor comorbidity with respiratory disease was significantly associated with the secondary outcomes. The dilution of this association may be explained by selection bias from clinicians on the suitability for either mechanical ventilation or CPAP for patients with these comorbidities.
This study has a number of strengths. Crucially, it included all sequentially admitted patients during a predetermined time period thus reducing the selection bias which has been observed in many studies. This study includes admission observations and investigations, accounting for disease severity in the model. It also benefits from a high degree of data completeness, with patient characteristics cross-validated with discharge summaries and other information sources to ensure accuracy.
A limitation of this study is that for the purpose of the time to event analyses, it had to be assumed that those who had clinically improved and were discharged prior to censoring remained alive until the censoring date. Although it is possible that some of these patients died post-discharge, this risk was minimised by ensuring that patients were screened for readmission and by ensuring that none of the cohort were discharged for community end-of-life care.
A second limitation is that a significant minority of patients (157/981) did not have an ethnicity documented. The majority of these patients had opted not to have their ethnicity recorded.
In conclusion, this large cohort study included all patients who were hospitalised with COVID-19 during the study period. It was conducted in an area of London with high levels of ethnic diversity [15]. When considering the high mortality among those hospitalised with COVID-19, early identification of those most susceptible to a severe disease is of paramount importance. This study found that predominantly respiratory features such as respiratory rate, oxygen saturations, infiltrates on chest radiograph and a history of respiratory disease were associated with severe disease. Comparatively few non-respiratory features were also identified as being associated with severe disease (age, reduced GCS and immunosuppression) along with a number of biochemical markers. For this new disease, understanding at the point of admission which patients are more likely to have a severe disease course is key to making appropriate and timely treatment escalation decisions, and ultimately reducing avoidable mortality.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S0950268820002472
Acknowledgements
We thank Dr Poppy Utting, Dr Megan Miscampbell, Dr Dorina Condurache, Dr Nathan Eager and Dr Krisha Shah for their help in collecting patient data. We are grateful to all our colleagues at London North West Hospitals University Healthcare NHS Trust for their outstanding commitment to maintaining good clinical care during the COVID-19 pandemic and for their support of this project.
Author contributions
NV, DLC and TANR designed the study. MA, JWG and TANR collected the data with assistance from those acknowledged below. PB and JWG conducted the data analysis. JWG wrote the first draft of the article and conducted the literature search. TANR, MA, PB, AMW, DLC and NV all reviewed and approved the final report.
Conflict of interest
All authors declare no conflict of interest.
Data availability statement
Data on all COVID-19 admissions to the Hospital have already been shared with Public Health England. Application for further patient level data should be addressed to the Hospital's Research Committee. STATA codes can be obtained from the authors.
Open access
