ABSTRACT
Background There is concern about medium to long-term adverse outcomes following acute COVID-19, but little relevant evidence exists. We aimed to investigate whether risks of hospital admission and death, overall and by specific cause, are raised following discharge from a COVID-19 hospitalisation.
Methods and Findings Working on behalf of NHS-England, we used linked primary care and hospital data in OpenSAFELY to compare risks of hospital admission and death, overall and by specific cause, between people discharged from COVID-19 hospitalisation (February-December 2020), and (i) demographically-matched controls from the 2019 general population; (ii) people discharged from influenza hospitalisation in 2017-19. We used Cox regression adjusted for personal and clinical characteristics.
24,673 post-discharge COVID-19 patients, 123,362 general population controls, and 16,058 influenza controls were followed for ≤315 days. Overall risk of hospitalisation or death (30968 events) was higher in the COVID-19 group than general population controls (adjusted-HR 2.23, 2.14-2.31) but similar to the influenza group (adjusted-HR 0.94, 0.91-0.98). All-cause mortality (7439 events) was highest in the COVID-19 group (adjusted-HR 4.97, 4.58-5.40 vs general population controls and 1.73, 1.60-1.87 vs influenza controls). Risks for cause-specific outcomes were higher in COVID-19 survivors than general population controls, and largely comparable between COVID-19 and influenza patients. However, COVID-19 patients were more likely than influenza patients to be readmitted/die due to their initial infection/other lower respiratory tract infection (adjusted-HR 1.37, 1.22-1.54), and to experience mental health or cognitive-related admission/death (adjusted-HR 1.36, 1.01-2.83); in particular, COVID-19 survivors with pre-existing dementia had higher risk of dementia death. One limitation of our study is that reasons for hospitalisation/death may have been misclassified in some cases due to inconsistent use of codes.
Conclusions People discharged from a COVID-19 hospital admission had markedly higher risks for rehospitalisation and death than the general population, suggesting a substantial extra burden on healthcare. Most risks were similar to those observed after influenza hospitalisations; but COVID-19 patients had higher risks of all-cause mortality, readmissions/death due to the initial infection, and dementia death, highlighting the importance of post-discharge monitoring.
Introduction
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged in early 2020 and rapidly spread around the world, infecting >140 million people globally.1 Acute infection can be asymptomatic or mild,2 but a substantial minority of infected people experience severe COVID-19 requiring hospitalisation,3 with age being a major risk factor, along with male sex, non-White ethnicity and certain comorbidities.4-6 Early in the pandemic, the proportion surviving hospitalisation was around 50-70%,7 though improved treatment guidelines and the identification of effective therapies such as dexamethasone helped to improve survival rates.8,9 There is now a large and growing population of people who have survived a COVID-19 hospitalisation, but little is known about their longer-term health outcomes.
Emerging evidence suggests that a subset of people infected with SARS-CoV-2 can experience health problems for at least several months after the acute phase of their infection, with fatigue, pain, respiratory and cardiovascular symptoms, and mental health and cognitive disturbances being among the problems frequently described under the term “post-acute COVID-19 syndrome”;10 however, epidemiological characterisation of such sequelae remains immature. Small descriptive studies of COVID-19 survivors have been suggestive of high incidence of a range of outcomes including respiratory, cardiovascular, and mental-health related, but firm conclusions are difficult due to lack of comparison groups.11,12 There remains limited evidence about post-COVID sequelae across the full range of health outcomes. One recent study of US Department of Veterans Affairs (VA) data examined a wide range of diagnoses, prescriptions and laboratory abnormalities among 30-day survivors of COVID-19, showing excess risks of several health outcomes in the 6 months following infection, compared with the general VA population.13 Whether these findings generalise to the entire US population or other settings remains unclear. Another US study limited to people aged <65 years also found excess risks of a range of clinical outcomes ascertained from health insurance data among people with a record of SARS-CoV-2 infection.14 A UK study of routinely-collected primary care and hospitalisation data described raised rates of all-cause hospital admission and death among patients discharged following a COVID-19 hospitalisation; the authors also noted raised risks of adverse respiratory and cardiovascular sequalae among the selected outcomes investigated.15 Only a general population comparator was used, making it difficult to disentangle risks specific to COVID-19 from those associated with hospitalisation more generally; furthermore, a hospitalised cohort is likely to have been more prone to health problems at the outset than the general population comparator group.
To strengthen the evidence base in this important emerging area, we therefore aimed to investigate the incidence of subsequent hospital admission and death, both overall and from a wide range of specific causes, following a COVID-19 hospitalisation in England. We aimed to compare post-COVID risks with two separate comparison groups: (i) the general population, and (ii) people hospitalised for influenza prior to the current pandemic.
Methods
Study design and study population
A cohort study was carried out within the OpenSAFELY platform, which has been described previously.6 We used routinely-collected electronic data from primary care practices using TPP SystmOne software, covering approximately 40% of the population in England, linked at the individual patient level to NHS Secondary Uses Service (SUS) data on hospitalisations, and Office of National Statistics death registration data (from 2019 onwards). We selected all individuals discharged between 1st February and 30th December 2020 from a hospitalisation that lasted >1 day and where COVID-19 was coded as the primary diagnosis (based on the International Classification of Diseases (ICD)-10 codes U07.1 “COVID-19 - virus identified” and U07.2 “COVID-19 - virus not identified”) and who were alive and under follow-up one week after discharge (to avoid a focus on hospital transfers and immediate readmissions/deaths). We excluded a small number of people with missing age, sex, or index of multiple deprivation, which are likely to indicate poor data quality. Two comparison groups were also selected: (i) people under follow-up in the general population in 2019, individually matched 5:1 to the COVID-19 group on age (within 3 years), sex, Sustainability and Transformation Plans (STP, a geographical area used as in NHS administration, of which there were 32 in our data), and calendar month (e.g. a patient discharged from a COVID-19 hospitalisation in April 2020 was matched to 5 individuals of the same age, sex and STP who were under follow-up in general practice on 1st April 2019); (ii) all individuals discharged from hospital in 2017-2019 where influenza was coded as the primary reason for hospitalisation and who were alive and under follow-up one week after discharge.
Outcomes and covariates
The outcomes were (i) time to first hospitalisation or death (composite outcome); (ii) all-cause mortality; and (iii) time to first cause-specific hospitalisation or death. Hospitalisations were identified from linked SUS data. All-cause mortality was identified using date of death in the primary care record so that deaths before 2019 (in the influenza group) could be included (as linked ONS data were not available prior to 2019); concordance of death dates between primary care and linked ONS data has been shown to be high.16 Cause-specific outcomes were categorised based on ICD-10 codes into infections (ICD-10 codes beginning with “A”), cancers except non-melanoma skin cancer (C, except C44), endocrine/nutritional/metabolic (E), mental health and cognitive (F, G30 and X60-84), nervous system (G, except G30), circulatory (I), COVID-19/influenza/pneumonia/other lower respiratory tract infections (J09-22, U07.1/2), other respiratory (J23-99), digestive (K), musculoskeletal (M), genitourinary (N), and external causes (S-Y, except X60-84). For each of these, the outcome was time to the earliest of hospitalisation with the relevant outcome listed as primary diagnosis, or death with the relevant outcome listed as the underlying cause on the death certificate.17 The influenza control group was restricted to those discharged in 2019 for analyses of these cause-specific outcomes, because we did not have linked death registration data (and thus cause of death) for earlier years.
Other covariates considered in the analysis were factors that might be associated with both risk of severe COVID-19 and subsequent outcomes, namely age, sex, ethnicity, obesity, smoking status, index of multiple deprivation quintile (derived from the patient’s postcode at lower super output area level), and comorbidities considered potential risk factors for severe COVID-19 outcomes (see Table 1 and footnotes for full specification of covariate categories and comorbidities).
Information on all covariates was obtained by searching TPP SystmOne records for specific coded data, based on a subset of SNOMED-CT mapped to Read version 3 codes. Covariates were identified using data prior to the patient’s hospital admission date (for the COVID-19 and influenza groups) or the index date (for the matched control group, i.e. 1st day of the matched calendar month in 2019). For the COVID-19 and influenza hospitalised groups, primary care data on ethnicity was supplemented with information from the hospitalisation record, to improve completeness. All codelists, along with detailed information on their compilation are available at https://codelists.opensafely.org for inspection and re-use by the wider research community.
Statistical analysis
Follow-up began on the 8th day after hospital discharge for the COVID-19 and influenza groups, and on the 1st of the same calendar month in 2019, for the general population control group. Follow-up ended at the first occurrence of the analysis-specific outcome, or the earliest relevant censoring date for data availability/coverage for the outcome being analysed; the control groups were additionally censored after the maximum follow-up time of the COVID-19 group (315 days). For outcomes involving hospital admissions the administrative censoring date (for SUS data) was 30th December 2020; for outcomes involving cause of death, the administrative censoring date (for ONS mortality data) was 11th March 2021; for the all-cause death outcome which was ascertained in primary care data, patients were censored at date of deregistration if they had left the TPP general practice network. Cumulative incidence of the composite hospitalisation/death outcome and all-cause mortality were calculated using Kaplan-Meier methods. Hazard ratios comparing COVID-19 and controls were estimated using Cox regression models. Separate models were fitted for the comparisons with matched 2019 general population controls (models stratified by matched set), and with influenza controls (models adjusted for age, sex, STP and calendar month). The additional covariates noted above were then added to the models. People with missing data on ethnicity, body mass index or smoking were excluded from models that used these variables (“complete case analysis”, which is valid under the assumption that missingness is conditionally independent of the outcome);18 imputation was not used because these variables were thought to be missing not at random in primary care (e.g. smokers more likely to have smoking status recorded). Cumulative incidence of cause-specific hospitalisation/death outcomes were calculated with deaths from other causes treated as a competing risk. Hazard ratios for these outcomes were then estimated from a Cox model targeting the cause-specific hazard, with deaths from competing risks censored.
The study was approved by the Health Research Authority (REC reference 20/LO/0651) and by the LSHTM Ethics Board (ref 21863).
Results
24,673 individuals discharged after a COVID-19 hospitalisation were included, alongside 123,362 matched controls from the 2019 general population, and 16,058 individuals discharged after influenza hospitalisation in 2017-19 (Appendix Figure A1-A2).
At entry, the COVID-19 group had similar age and sex distribution to the general population groups due to matching, but had younger median age and were more likely to be male than the influenza group (Table 1). Body mass index, smoking and ethnicity data were 93-99% complete in all groups, except that ethnicity was 25% missing in the matched control group (no hospital-based ethnicity records were available for this group). The COVID-19 group were more likely to be obese, non-White and less likely to be current smokers than both comparison groups. Pre-existing comorbidities were more common in both COVID-19 and influenza-discharged patients than general population controls. COVID-19 patients had longer median duration of hospital stay and were more likely to have received critical care during their admission than influenza patients.
Numbers of outcome events are shown in Appendix Table A1. Cumulative incidence of subsequent hospital admission or death after study entry in the COVID-19 group was higher than in general population controls, but similar to that in the influenza group (cumulative incidence at 6 months = 34.8, 15.2 and 37.8% in the three groups respectively; fully adjusted HR = 2.23, 2.14-2.31 for COVID-19 vs general population; 0.94, 0.91-0.98 for COVID-19 vs influenza, Figures 1a, 2). Cumulative all-cause mortality was higher in the COVID-19 group than both the general population and influenza groups (7.5, 1.4 and 4.9% at 6 months in the three groups respectively; fully adjusted HR = 4.97, 4.58-5.40 for COVID-19 vs general population; 1.73, 1.60-1.87 for COVID vs influenza, Figures 1b, 2). To further explore this, causes of death were examined (appendix Table A2). A substantial proportion of deaths in the COVID-19 group had COVID-19 listed as the underlying cause (500/2022, 24.7%), while in the influenza group ≤5 deaths were coded with influenza as the underlying cause.
Cumulative incidences of cause-specific hospital admissions or deaths are shown in Figure 3. After adjustment for matching factors and other covariates, risks of all cause-specific outcomes were substantially higher in COVID-19 groups than general population controls (Figure 4). Compared with influenza patients, people in the COVID-19 group had similar or lower risk of admission or death from most causes, but higher risks of admission/death from COVID-19/influenza/lower respiratory tract infection (LRTI, adjusted HR 1.37, 1.22-1.54); in the post-COVID-19 group these outcomes were dominated by codes for COVID-19 itself (515/1122 [46%] of hospitalisations and 342/368 [93%] of deaths) and pneumonia (461/1122 [41%] of hospitalisations). The COVID-19 group also had higher risks than the influenza group for mental health or cognitive outcomes (adjusted HR 1.36, 1.01-1.83). This was further explored in a post-hoc analysis of specific outcomes within the mental health and cognitive category (Table 2). Raised risks in the COVID-19 group appeared to be driven by dementia hospitalisations/deaths (age/sex-adjusted HR 2.32, 1.48-3.64), particularly among those with pre-existing dementia at baseline (HR 2.47, 1.37-4.44) and/or resident in care homes (HR 2.53, 0.99-6.41). Of note, 129/161 dementia outcome events (80.1%) were deaths (rather than hospitalisations). Higher rates of hospitalisations/deaths due to mood disorders, and neurotic/stress-related/somatoform disorders were also observed in COVID-19 patients, but confidence intervals were too wide to be conclusive.
Discussion
Key findings
Patients discharged from a COVID-19 hospitalisation had more than double the risk of hospitalisation or death, and a 5-fold higher risk of all-cause mortality than controls from the general population, after adjusting for baseline personal and clinical characteristics. Risks were higher for all categories of disease-specific hospital admissions/deaths after a COVID-19 hospitalisation than in general population controls. Most outcomes were similar between people discharged from a COVID-19 hospitalisation, and people discharged from an influenza hospitalisation in 2017-19, but the COVID-19 group had higher subsequent all-cause mortality; higher rates of respiratory infection admissions and deaths (predominantly COVID-19), and more adverse mental health and cognitive outcomes (particularly deaths attributed to dementia among people with pre-existing dementia) compared with the influenza group.
Findings in context of literature and explanations
A recent study of VA data on US veterans examined a wide range of diagnostic and other outcomes in 30-day COVID-19 survivors, compared with the general VA population.13 Among veterans where COVID-19 had led to a hospitalisation, hazard ratios of every category of outcome were raised. This concurs with findings from our study, despite different characteristics of the VA population. In the UK, an earlier study found an 8-fold higher risk of death in post-acute COVID-19 patients compared with general population controls, and raised risks of respiratory disease, diabetes and cardiovascular disease.15 Interestingly, recent data from Denmark suggest limited post-acute complications following non-hospitalised COVID-19;19 this is in contrast to a recent study using US health insurance data which found raised risks of a range of outcomes among a relatively young cohort with mostly (92%) non-hospitalised COVID-19 disease, compared with both the general population and people with a record of other viral lower respiratory tract infections.14
Our data showed that COVID-19 hospitalised patients were more likely to have baseline comorbidities than general population controls, reflecting known associations between comorbidities and risks of severe COVID-19 outcomes.6 Differences in outcomes between hospitalised patients and general population controls might therefore reflect baseline differences not fully captured in our adjustment models, and might also reflect a generic adverse effect of hospitalisation.20 This is supported by the more similar risks we observed when COVID-19 survivors were compared with people who had experienced influenza hospitalisation. However, all-cause mortality was substantially higher after COVID-19 compared with influenza. A quarter of deaths after a COVID-19 hospitalisation had COVID-19 listed as the underlying cause, but it is not clear from our data whether patients experienced specific complications after hospital discharge that were then attributed to COVID-19.
Our analysis of cause-specific outcomes also suggested a disproportionate rate of dementia deaths post-COVID-19, particularly among those with pre-existing dementia. Cognitive decline after hospitalisation and critical illness have been previously described;21,22 acute COVID-19 and associated hospital admission, social isolation, and medications may have accelerated progression of patients’ dementia; it is unclear whether post-discharge care was adequate for this vulnerable group. However, it is possible that deaths where the underlying cause was recorded as dementia may have been due to progression of underlying health problems following an acute illness as well as difficulty in managing these due to dementia. Due to small numbers, we could not confirm whether higher rates of mood disorders and neurotic/stress-related/somatoform disorders after COVID-19 compared with influenza were due to chance, but a number of previous studies outside the pandemic context have found that critical illness is associated with raised risks of depression, anxiety and post-traumatic stress.23-25 It will be important to continue to monitor these outcomes as more follow-up accumulates.
Strengths and Limitations
We identified COVID-19 hospitalisations from a base population covering around 40% of the population of England, giving us high statistical power. We examined a broad range of hospitalisation and mortality outcomes, and were able to describe and adjust for a wide range of personal and clinical characteristics using rich primary care data.
However, our study has some limitations. We relied on ICD-10 codes entered as the primary reason for hospitalisation or underlying cause of death to define our cause-specific outcomes, but these fields may not have been used consistently.26 In particular, there might have been a tendency for clinicians aware of a recent COVID-19 hospitalisation to code COVID-19 for a range of clinical complications, masking more specific sequelae. Outcomes were classified in broad categories to obtain an overview of post-COVID-19 disease patterns; more granular disease categories would be of future interest but will require more follow-up to maintain statistical power. Our main comparisons may have been affected by time-related factors. We compared post-COVID patients in 2020 with controls from 2019 and earlier; consultations for non-COVID-19 conditions in 2020 are known to have been subdued in the general population,27 perhaps due to lockdown, or public reluctance to seek care, potentially affecting comparison with earlier years. On the other hand, patients with recent a COVID-19 hospitalisation may assume immunity from reinfection and be less reticent in seeking care than the general population. The comparison with influenza may also have been affected by seasonality, since the first wave of COVID-19 in England happened outside the typical influenza season. Lack of overlap in the data meant that we could not incorporate seasonal adjustment into our statistical models for this comparison; any confounding by seasonality is likely to have led to underestimation of hazard ratios comparing COVID-19 and influenza patients, since cases of the former were under-represented in the winter months (which typically confer higher health risks). We did not have detailed data on disease severity, though descriptive data showed that COVID-19 patients tended to have longer hospital stays and more critical care than those hospitalised for influenza. Finally, data were unavailable on new/emerging COVID-19 variants during the study period.
Conclusions
Patients surviving a COVID-19 hospitalisation were at substantially higher risk than the general population for a range of subsequent adverse outcomes over a period of up to 10 months’ follow-up included in this study. Risks for most outcomes were broadly comparable to those experienced by influenza hospitalisation survivors prior to the pandemic, but in the period following hospital discharge COVID-19 patients had higher risks of all-cause mortality, readmission or death attributed to their initial infection, and adverse mental health and cognitive outcomes; in particular, among people with pre-existing dementia, we observed an excess of deaths where dementia was recorded as the underlying cause. These findings suggest a need for services to support and closely monitor people following discharge from hospital with COVID-19.
Data Availability
All data were linked, stored and analysed securely within the OpenSAFELY platform https://opensafely.org/. Data include pseudonymized data such as coded diagnoses, medications and physiological parameters. No free text data are included. All code is shared openly for review and re-use under MIT open license (https://github.com/opensafely/post-admission-admissions-research). Detailed pseudonymised patient data is potentially re-identifiable and therefore not shared. We rapidly delivered the OpenSAFELY data analysis platform without prior funding to deliver timely analyses on urgent research questions in the context of the global Covid-19 health emergency: now that the platform is established we are developing a formal process for external users to request access in collaboration with NHS England; updates on this process are available on OpenSAFELY.org.
Contributors
KB, CTR, WJH led the study design. GH led on platform support. KB did the analysis and wrote the first draft. BG conceived the OpenSAFELY platform and the approach. BG and LS lead the project overall and are guarantors. SB led on software development. AM led on information governance. Other contributions were – data curation - CB, JP, JC, SH, SB, DE, PI and CEM; disease category conceptualisation and codelists - KB CTR BM CB JC CEM AJW HIM IJD HJC JP; statistical analysis code: KB, EW; ethical approvals - HJC EW LS BG; software - SB, DE, PI, AJW, WJH, CEM, CB, FH, JC; writing (reviewing and editing) - KB, CTR, CWG, LT, BM, AS, AM, RME, CEM, IJD, SJWE, LS, BG. All authors were involved in design and conceptual development and reviewed and approved the final manuscript.
Information Governance and Ethics
NHS England is the data controller; TPP is the data processor; and the key researchers on OpenSAFELY are acting on behalf of NHS England. OpenSAFELY is hosted within the TPP environment which is accredited to the ISO 27001 information security standard and is NHS IG Toolkit compliant;28,29 patient data are pseudonymised for analysis and linkage using industry standard cryptographic hashing techniques; all pseudonymised datasets transmitted for linkage onto OpenSAFELY are encrypted; access to the platform is via a virtual private network (VPN) connection, restricted to a small group of researchers who hold contracts with NHS England and only access the platform to initiate database queries and statistical models. All database activity is logged. No patient-level data leave the platform; only aggregate statistical outputs leave the platform environment following best practice for anonymisation of results such as statistical disclosure control for low cell counts.30 The OpenSAFELY platform adheres to the data protection principles of the UK Data Protection Act 2018 and the EU General Data Protection Regulation (GDPR) 2016. In March 2020, the Secretary of State for Health and Social Care used powers under the UK Health Service (Control of Patient Information) Regulations 2002 (COPI) to require organisations to process confidential patient information for the purposes of protecting public health, providing healthcare services to the public and monitoring and managing the COVID-19 outbreak and incidents of exposure.31 Taken together, these provide the legal bases to link patient datasets on the OpenSAFELY platform. The study was approved by the Health Research Authority (REC reference 20/LO/0651) and by the LSHTM Ethics Board (ref 21863).
Declaration of interests
This work was jointly funded by UKRI, NIHR and Asthma UK-BLF [COV0076; MR/V015737/] and the Longitudinal Health and Wellbeing strand of the National Core Studies programme. The OpenSAFELY data science platform is funded by the Wellcome Trust. TPP provided technical expertise and infrastructure within their data centre pro bono in the context of a national emergency. BG has received research funding from the Laura and John Arnold Foundation, the NHS National Institute for Health Research (NIHR), the NIHR School of Primary Care Research, the NIHR Oxford Biomedical Research Centre, the Mohn-Westlake Foundation, NIHR Applied Research Collaboration Oxford and Thames Valley, the Wellcome Trust, the Good Thinking Foundation, Health Data Research UK, the Health Foundation, the World Health Organisation, UKRI, Asthma UK, the British Lung Foundation, and the Longitudinal Health and Wellbeing strand of the National Core Studies programme; he also receives personal income from speaking and writing for lay audiences on the misuse of science. LS reports grants from Wellcome, MRC, NIHR, UKRI, British Council, GSK, British Heart Foundation, and Diabetes UK outside this work. KB holds a Senior Research Fellowship from Wellcome (220283/Z/20/Z). CWG is supported by a Wellcome Intermediate Clinical Fellowship (201440/Z/16/Z). HIM is funded by the National Institute for Health Research (NIHR) Health Protection Research Unit in Immunisation, a partnership between Public Health England and LSHTM. EW holds grants from MRC. ID golds grants from NIHR and GSK. HF holds a UKRI fellowship. RME is funded by HDR UK (grant: MR/S003975/1) and MRC (grant: MC_PC 19065). The views expressed are those of the authors and not necessarily those of the NIHR, NHS England, Public Health England or the Department of Health and Social Care. Funders had no role in the study design, collection, analysis, and interpretation of data; in the writing of the report; and in the decision to submit the article for publication.
Open access statement
This research was funded in whole or in part by the Wellcome Trust. For the purpose of Open Access, the author has applied a CC BY public copyright licence to any Author Accepted Manuscript (AAM) version arising from this submission.
Data availability statement
All data were linked, stored and analysed securely within the OpenSAFELY platform https://opensafely.org/. Data include pseudonymized data such as coded diagnoses, medications and physiological parameters. No free text data are included. All code is shared openly for review and re-use under MIT open license (https://github.com/opensafely/post-admission-admissions-research). Detailed pseudonymised patient data is potentially re-identifiable and therefore not shared. We rapidly delivered the OpenSAFELY data analysis platform without prior funding to deliver timely analyses on urgent research questions in the context of the global Covid-19 health emergency: now that the platform is established we are developing a formal process for external users to request access in collaboration with NHS England; updates on this process are available on OpenSAFELY.org.
Patient and Public Involvement
Patients were not formally involved in developing this specific study design that was developed rapidly in the context of a global health emergency. We have developed a publicly available website https://opensafely.org/ through which we invite any patient or member of the public to contact us regarding this study or the broader OpenSAFELY project.
Appendix
Footnotes
↵+ joint senior authors