Intended for healthcare professionals

CCBY Open access
Research

Risk prediction of covid-19 related death and hospital admission in adults after covid-19 vaccination: national prospective cohort study

BMJ 2021; 374 doi: https://doi.org/10.1136/bmj.n2244 (Published 17 September 2021) Cite this as: BMJ 2021;374:n2244

Read our latest coverage of the coronavirus pandemic

  1. Julia Hippisley-Cox, professor of clinical epidemiology and general practice1,
  2. Carol AC Coupland, senior research fellow and professor of medical statistics in primary care1 2,
  3. Nisha Mehta, clinical director for digital primary care3,
  4. Ruth H Keogh, professor of biostatistics and epidemiology4,
  5. Karla Diaz-Ordaz, associate professor of biostatistics4,
  6. Kamlesh Khunti, professor of primary care diabetes and vascular medicine5,
  7. Ronan A Lyons, professor of public health6,
  8. Frank Kee, professor of public health medicine7,
  9. Aziz Sheikh, professor of primary care research and development8,
  10. Shamim Rahman, deputy director for mental health, disability and shielding analysis9,
  11. Jonathan Valabhji, national clinical director for diabetes and obesity10,
  12. Ewen M Harrison, professor of surgery and data science8,
  13. Peter Sellen, lead analyst for covid-19 clinically extremely vulnerable9,
  14. Nazmus Haq, senior analyst for covid-19 clinically extremely vulnerable9,
  15. Malcolm G Semple, professor of child health and outbreak medicine11,
  16. Peter W M Johnson, national clinical director for cancer10,
  17. Andrew Hayward, professor of infectious disease epidemiology and inclusion health12,
  18. Jonathan S Nguyen-Van-Tam, professor of health protection and deputy chief medical officer2 9
  1. 1Nuffield Department of Primary Health Care Sciences, University of Oxford, Oxford, UK
  2. 2Division of Primary Care, School of Medicine, University of Nottingham, Nottingham, UK
  3. 3NHS-X, London, UK
  4. 4Department of Medical Statistics and Centre for Statistical Methodology, London School of Hygiene and Tropical Medicine, London, UK
  5. 5Diabetes Research Centre, University of Leicester, Leicester, UK
  6. 6Population Data Science, Swansea University, Swansea, UK
  7. 7Queen’s University, Belfast, UK
  8. 8Usher Institute, University of Edinburgh, Edinburgh, UK
  9. 9Department of Health and Social Care, England, UK
  10. 10NHS England and Improvement, London, UK
  11. 11NIHR Health Protection Research Unit, Institute of Infection, Veterinary and Ecological Sciences, Faculty of Health and Life Sciences, University of Liverpool, Liverpool, UK
  12. 12UCL Institute of Epidemiology and Health Care, London, UK
  1. Correspondence to: J Hippisley-Cox julia.hippisley-cox{at}phc.ox.ac.uk (or @JuliaHCox on Twitter)
  • Accepted 13 September 2021

Abstract

Objectives To derive and validate risk prediction algorithms to estimate the risk of covid-19 related mortality and hospital admission in UK adults after one or two doses of covid-19 vaccination.

Design Prospective, population based cohort study using the QResearch database linked to data on covid-19 vaccination, SARS-CoV-2 results, hospital admissions, systemic anticancer treatment, radiotherapy, and the national death and cancer registries.

Settings Adults aged 19-100 years with one or two doses of covid-19 vaccination between 8 December 2020 and 15 June 2021.

Main outcome measures Primary outcome was covid-19 related death. Secondary outcome was covid-19 related hospital admission. Outcomes were assessed from 14 days after each vaccination dose. Models were fitted in the derivation cohort to derive risk equations using a range of predictor variables. Performance was evaluated in a separate validation cohort of general practices.

Results Of 6 952 440 vaccinated patients in the derivation cohort, 5 150 310 (74.1%) had two vaccine doses. Of 2031 covid-19 deaths and 1929 covid-19 hospital admissions, 81 deaths (4.0%) and 71 admissions (3.7%) occurred 14 days or more after the second vaccine dose. The risk algorithms included age, sex, ethnic origin, deprivation, body mass index, a range of comorbidities, and SARS-CoV-2 infection rate. Incidence of covid-19 mortality increased with age and deprivation, male sex, and Indian and Pakistani ethnic origin. Cause specific hazard ratios were highest for patients with Down’s syndrome (12.7-fold increase), kidney transplantation (8.1-fold), sickle cell disease (7.7-fold), care home residency (4.1-fold), chemotherapy (4.3-fold), HIV/AIDS (3.3-fold), liver cirrhosis (3.0-fold), neurological conditions (2.6-fold), recent bone marrow transplantation or a solid organ transplantation ever (2.5-fold), dementia (2.2-fold), and Parkinson’s disease (2.2-fold). Other conditions with increased risk (ranging from 1.2-fold to 2.0-fold increases) included chronic kidney disease, blood cancer, epilepsy, chronic obstructive pulmonary disease, coronary heart disease, stroke, atrial fibrillation, heart failure, thromboembolism, peripheral vascular disease, and type 2 diabetes. A similar pattern of associations was seen for covid-19 related hospital admissions. No evidence indicated that associations differed after the second dose, although absolute risks were reduced. The risk algorithm explained 74.1% (95% confidence interval 71.1% to 77.0%) of the variation in time to covid-19 death in the validation cohort. Discrimination was high, with a D statistic of 3.46 (95% confidence interval 3.19 to 3.73) and C statistic of 92.5. Performance was similar after each vaccine dose. In the top 5% of patients with the highest predicted covid-19 mortality risk, sensitivity for identifying covid-19 deaths within 70 days was 78.7%.

Conclusion This population based risk algorithm performed well showing high levels of discrimination for identifying those patients at highest risk of covid-19 related death and hospital admission after vaccination.

Introduction

During the first waves of the covid-19 pandemic (March 2020 to August 2020), before the introduction of vaccines, it was essential to be able to identify people at highest risk of adverse outcomes if they were infected with SARS-CoV-2. The QCovid risk assessment tool for predicting risk of covid-19 related death or hospital admission based on individual characteristics was developed,1 independently externally validated,2 and found to have performed well at identifying those individuals at high risk of severe outcomes from covid-19. The tool was used in England to identify patients at high risk of severe covid-19 outcomes, adding an additional 1.5 million people to the national shielded patient list in February 2021 and, on a UK basis, prioritising them for vaccination (if they had not already been offered the vaccine on account of their age or occupation).3

Since then, clinical trials of covid-19 vaccinations have demonstrated safety and efficacy in healthy volunteers456 and have been rolled out to the adult UK population, beginning with the most elderly groups (aged ≥90 years) and those people most at risk. Although vaccines have been found to be highly effective in trials and observational studies, a residual risk of serious covid-19 outcomes (in particular, hospital admission or death) remains after vaccination, despite allowing adequate time for immunity to develop. The risk of a severe outcome in vaccinated groups includes the risk of exposure, the risk of a breakthrough infection if exposed, and the risk of a breakthrough infection becoming severe. However, the relevant risk factors are currently unknown because clinical trials have not included many people in whom vaccine response might be suboptimal (eg, elderly people, people with complex comorbidities (eg, in receipt of solid organ transplants or immunosuppressive treatment for autoimmune disorders), or patients with cancer receiving chemotherapy or radiotherapy7).

Therefore, vaccinated individuals at highest risk of consequent severe outcomes such as covid-19 related hospital admission or death need to be identified urgently. A risk stratification tool for the vaccinated population would enable identification of patients to prioritise for targeted, early interventions once these become available—including booster vaccination and preventive treatments such as passive antibody delivery (for either prophylactic or therapeutic use). Risk stratification tools also provide a robust pragmatic mechanism for avoiding unnecessary lifestyle precautions, investigations, and therapeutic interventions for those individuals whose risk is relatively low, but who might perceive it to be much higher.

We developed and validated two new QCovid risk algorithms, based on data from the second pandemic wave in England, to identify those groups at highest risk of severe covid-19 outcomes: QCovid2 (based on unvaccinated patients) and QCovid3 (based on vaccinated patients). Given the nature of the pandemic, the speed of the vaccination programme, the relaxation of lockdown measures, and the urgent need to develop national policy, this work had to be undertaken during the pandemic period and national vaccination programme, and therefore included people who had only one vaccination dose as well as those who were fully vaccinated.

Methods

Data sources

We used the QResearch database (version 46) of 12 million patients with personal, clinical, and drug data that have been used for clinical18 and drug safety research.910 QResearch is linked to multiple datasets at individual patient level. For this analysis, we used the National Immunisation Database of covid-19 vaccinations to identify data on vaccine date and doses for all people vaccinated in England. For hospital admissions, we used the linked Hospital Episode Statistics dataset supplemented by the more regularly updated Secondary Users Service data. We also obtained and linked the following datasets: national data for mortality; SARS-CoV-2 infection; systemic anticancer treatment; radiotherapy treatment datasets; and national cancer registry data.

Study design and period for vaccinated cohort

We undertook a prospective cohort study of vaccinated individuals from 8 December 2020 (the earliest vaccination date in England) to 15 June 2021 (the latest date for which data were available at the time of the analysis). We considered outcomes after the first and second vaccination doses. The cohort included people who received one or two doses of a covid-19 vaccine during the study period. Individuals were followed from 14 days after receiving each vaccine dose until they had the outcome of interest, died, or reached the end of the study period. Use of follow-up time after the first dose and after the second dose is described below.

Inclusion criteria for vaccinated cohort

We included all adults aged 19-100 years who had one or two doses of the ChAdOx1 nCoV-19 (Oxford-AstraZeneca) or BNT162b2 (Pfizer-BioNTech) vaccine during the study period. Both these vaccinations require two doses for full vaccination. People were excluded from the analysis of hospital outcomes if they had a covid-19 associated hospital admission before their start of follow-up (14 days after the first or second dose of vaccination).

Outcomes for vaccinated cohort

The primary outcome was time to covid-19 related death (either in or out of hospital) as recorded on the death certification, or death within 28 days of a SARS-CoV-2 infection confirmed by reverse transcription polymerase chain reaction (RT-PCR). The secondary outcome was time to hospital admission with covid-19, defined as either confirmed or suspected covid-19 on ICD-10 (international classification of diseases, 10th revision) codes U071 and U072, or new hospital admission associated with a confirmed SARS-CoV-2 infection in the preceding 14 days. These outcome definitions are also used for covid-19 death and hospital admission in the UK.11 Both outcomes were assessed from 14 days or more after the first and second doses of vaccination, from which time we considered it was reasonable to expect some immunity.

Predictor variables for vaccinated cohort

Candidate predictor variables likely to be associated with increased risk of covid-19 death or hospital admission were identified from the original QCovid protocol12 and from previous studies.1813 The variables were vaccine dose (first or second), age, sex, ethnic origin, Townsend deprivation score (an area level score based on postcode where higher scores indicate higher levels of deprivation14), body mass index,13 domicile (care home, homeless, neither), chronic kidney disease, chemotherapy in previous 12 months, type 1 or type 2 diabetes (with glycated haemoglobin (HbA1c) levels <59 or ≥59 mmol/mol), blood cancer, bone marrow transplantation in past six months, respiratory cancer, radiotherapy in past six months, solid organ transplantation, chronic obstructive pulmonary disease, asthma, rare lung diseases (cystic fibrosis, bronchiectasis, or alveolitis), pulmonary hypertension or pulmonary fibrosis, coronary heart disease, stroke, atrial fibrillation, heart failure, venous thromboembolism, peripheral vascular disease, congenital heart disease, dementia, Parkinson’s disease, epilepsy, rare neurological conditions (motor neurone disease, multiple sclerosis, myasthenia gravis, or Huntington’s chorea), cerebral palsy, osteoporotic fracture, rheumatoid arthritis or systemic lupus erythematosus, liver cirrhosis, bipolar disorder or schizophrenia, inflammatory bowel disease, sickle cell disease, HIV/AIDS, severe combined immunodeficiency, and record of a SARS-CoV-2 positive test result before cohort entry.

To account for changing infection rates during the study period (since the vaccination programme was started during the second pandemic wave in England), we calculated a seven-day moving average of the background rates of positive SARS-CoV-2 tests per 100 000 people, using published English national data.15 We linked the rate to the date of cohort entry for each individual (that is, 14 days after each vaccine dose).

We defined predictors using information recorded in primary care electronic health records at the start of follow-up at 14 days after the first dose, except for chemotherapy, radiotherapy, and transplantations, which were based on linked data related to systemic anticancer and radiotherapy treatment, Hospital Episode Statistics data, and Secondary Users Service data. For all predictor variables, we used the most recently available value at the cohort entry date.

Model development

To maximise the number of events after second dose of vaccine, we used all 1336 practices with linked data available up to 15 June 2021 at the time of model development. We subsequently validated it in the remaining 182 practices once the updated linked data to 15 June 2021 became available a few weeks later.

People entered the cohort at 14 days after their first vaccination dose. We used a landmarking approach16 to handle the time dependent dose variable, because some people contributed follow-up time after their second dose as well as after their first dose. For people with only one dose, we followed them up until they had the event of interest, died, or reached the study end. For those with two vaccination doses, we split follow-up time into two periods. Period 1 included the time from 14 days after their first vaccination dose until 14 days after their second dose (therefore, outcomes during the first 14 days after the second dose were attributed to the first dose). Period 2 included time from 14 days after their second dose until they had the event of interest, died, or reached the study end. We fitted all models using combined data from follow-up after the first and second doses, with dose number entered into the model as a predictor.

We developed the risk models using cause specific Cox proportional hazard models to calculate hazard ratios and develop risk scores accounting for the competing risk of death due to other causes. A hazard ratio is a measure of the rate at which a particular outcome happens in one group relative to the rate at which it happens in another group over time. We fitted two cause specific Cox models to derive a risk algorithm for our primary outcome—one for covid-19 deaths and one for deaths due to other causes, censoring patients with the respective competing event. For our secondary outcome, we fitted one model for covid-19 admission and another model for all cause mortality (excluding deaths occurring after a covid-19 admission).

We used second degree fractional polynomials to model non-linear associations for continuous variables including SARS-Cov-2 infection rates, age, body mass index, and Townsend deprivation score.14 We fitted the models to the complete cases (that is, with no missing values for predictor variables) to derive the fractional polynomial terms. We used multiple imputation with chained equations to impute missing values for ethnic origin, Townsend score, body mass index, and HbA1c. We carried out five imputations and fitted the prediction models in each imputed dataset, and used Rubin’s rules to combine the model parameter estimates across the imputed datasets.17

We retained variables in the final models that were significant at the 5% level (taking account of the clustered nature of the data) or when adjusted cause specific hazard ratios for categorical variables were more than 1.1. Clinically similar variables with low numbers of events, such as bone marrow and solid organ transplantation, were combined. We examined interactions between predictor variables and age, as well as interactions between vaccine dose and age, body mass index, ethnic origin, deprivation, and each comorbidity. Furthermore, we derived estimates of the cumulative incidence function for covid-19 mortality accounting for the competing risk of death from other causes by combining estimates obtained from the two cause specific Cox models using an appropriate formula.18 The same method was used to derive the cumulative incidence function for covid-19 hospital admission, accounting for competing risk of death. These final algorithms for predicting absolute risk in vaccinated individuals are referred to as QCOVID3.18

We developed an additional model restricted to vaccinated patients with a positive SARS-CoV-2 test result after vaccination. This model separately quantified the risk of severe outcomes (mortality and admission) in individuals with a record of infection.

Model evaluation

We evaluated model performance in the separate validation cohort. We used multiple imputation to replace missing values for ethnic origin, body mass index, and Townsend score; the imputation model used was the same as that used in the derivation cohort. We applied the final risk equations to calculate the risk scores for each outcome accounting for competing risks, and calculated a C index accounting for competing risks using R.19 We also calculated R2 values and D statistics20 although these statistics were only available for the cause specific outcomes.

We assessed model calibration in the validation cohort accounting for competing risks by comparing mean predicted risks with the observed cumulative incidence function by twentieths of predicted risk.21 A model is well calibrated if predicted risks closely approximate the observed risks. We calculated each metric in the whole validation cohort, separately for individuals who had received one and two vaccination doses, and in subgroups for age and sex (ethnic groups had too few patients to undertake analyses).

Risk stratification

We applied the algorithms to the validation cohort to define the centile thresholds based on absolute risk using the prevailing SARS-CoV-2 rate 14 days after the date of each vaccination dose. Sensitivity was calculated as the total cumulative number of patients with a risk score above the risk threshold with a covid-19 death by 70 days divided by the total cumulative number of patients with a covid-19 death by 70 days.

QCovid2 model in the unvaccinated cohort

We also developed and evaluated two additional models (QCovid2) based on a cohort of unvaccinated people aged 19-100 years and observed between 1 September 2020 and 31 May 2021 but censoring people who were subsequently vaccinated on the date of their first vaccination. Additional variables not included in the original QCovid model were used, such as inflammatory bowel disease and levels of diabetes control according to HbA1c measurements. We also examined separate variables for sickle cell disease, HIV/AIDs, immunodeficiency conditions, and a refined definition of severe mental illness (to determine the contribution of moderate and severe depression). The first model included all unvaccinated patients (restricting to the time before vaccination for those who were subsequently vaccinated). The second model was restricted to unvaccinated patients with a positive SARS-CoV-2 test result, to separately quantify the risk of SARS-CoV-2 infection from the risk of severe outcomes (covid-19 mortality and admission) in those people with a positive test result. These final algorithms for predicting absolute risk in unvaccinated individuals are referred to as QCovid2.

Reporting

Stata (version 17) and R were used for analyses. The study adhered to the TRIPOD (transparent reporting of a multivariable prediction model for individual prognosis or diagnosis) statement for reporting.22

Patient and public involvement

Patients were involved in framing the research question, identifying predictors, and in developing plans for design and implementation of the study.

Results

Characteristics of the vaccinated cohort

Table 1 shows the characteristics of the 6 952 440 vaccinated patients in the derivation cohort, of whom 4 026 592 (57.9%) had the Oxford-AstraZeneca vaccine and 2 925 848 (42.1%) had the Pfizer-BioNTech vaccine. Overall, the mean age was 52 years (standard deviation 17.7), 3 321 247 (47.8%) were men, and 5 150 310 (74.1%) had two vaccine doses. The median follow-up time was 72 days (interquartile range 59-77) after the first dose and 35 (18-53) days after the second dose. Of 2031 covid-19 related deaths and 1929 covid-19 related hospital admissions, 81 deaths and 71 admissions occurred 14 days or more after the second vaccine dose. Of the 1929 patients in hospital, 446 (23.1%) subsequently died. Supplementary table 1 shows corresponding results for the 626 656 vaccinated patients in the validation cohort, of whom 174 had a covid-19 death and 179 had a covid-19 hospital admission. Of these, 10 deaths and seven admissions occurred 14 days or more after the second vaccine dose.

Table 1

Personal and medical characteristics for the derivation cohort and covid-19 related death or hospital admission 14 days or more after vaccination. Data are number (%) unless stated otherwise

View this table:

QCovid3: associations of outcomes with predictor variables

The final risk algorithms for covid-19 mortality included age, sex, ethnic origin, Townsend deprivation, body mass index, a range of comorbidities, and SARS-CoV-2 infection rate. We did not find any evidence of interactions between the dose variable and other predictors (although we did find small numbers for some pre-existing health conditions). Figure 1 and figure 2 show the adjusted hazard ratios and 95% confidence intervals for predictor variables included in the final cause specific models for covid-19 related deaths and hospital admissions. Supplementary figures 1-3 show the adjusted hazard ratios for the fractional polynomial terms for age, body mass index, and background SARS-CoV-2 infection rate, respectively. Hazard ratios increased with increasing age, Townsend deprivation, and background rates of SARS-CoV-2 infection. We saw a J shaped association between body mass index and rates of both hospital admission and mortality outcomes. Supplementary figures 4 and 5 show the adjusted hazard ratios for the competing events of non-covid-19 death and all cause deaths.

Fig 1
Fig 1

Use of QCovid3 model showing adjusted cause specific hazard ratios for covid-19 death after vaccination, mutually adjusted and adjusted for fractional polynomial terms for age, body mass index, vaccination dose, and background infection rate at time of vaccination. CKD=chronic kidney disease; HbA1c=glycated haemogoblin

Fig 2
Fig 2

Use of QCovid3 model showing adjusted cause specific hazard ratios for covid-19 hospital admission after vaccination, mutually adjusted and adjusted for fractional polynomial terms for age, body mass index, vaccination dose, and background infection rate at time of vaccination. CKD=chronic kidney disease; HbA1c=glycated haemogoblin

Covid-19 mortality incidence increased with age and deprivation, male sex, and Indian and Pakistani ethnic origin. Hazard ratios were highest for those with Down’s syndrome (12.7-fold increase), kidney transplantation (8.1-fold), sickle cell disease (7.7-fold), care home residency (4.1-fold), group B (3.6-fold) and group C chemotherapy (4.3-fold), recent bone marrow transplantation or a solid organ transplantation ever (2.5-fold), HIV/AIDS (3.3-fold), dementia (2.2-fold), Parkinson’s disease (2.2-fold), neurological conditions (2.6-fold), and liver cirrhosis (3.0-fold). Other conditions associated with increased covid-19 mortality included chronic kidney disease, blood cancer, epilepsy, chronic obstructive pulmonary disease, coronary heart disease, stroke, atrial fibrillation, heart failure, thromboembolism, peripheral vascular disease, and type 2 diabetes (with highest risks among those with HbA1c ≥59 mmol/mol (>7.5%)). The adjusted hazard ratio for covid-19 related death was 0.17 (95% confidence interval 0.13 to 0.22) after the second vaccine dose (plus 14 days) compared with after the first vaccine dose (plus 14 days).

We found similar patterns of associations between predictors and the cause specific hazard for covid-19 admission (fig 2) except for conditions with too few events for analysis (that is, sickle cell disease, severe combined immunodeficiency, and HIV/AIDS). Similarly, the adjusted hazard ratio of covid-19 related hospital admission was 0.21 (95% confidence interval 0.16 to 0.27) after the second dose compared with after the first dose.

Supplementary figure 6 shows the corresponding results for risk of covid-19 death among the subgroup of patients with a SARS-CoV-2 positive test result. The associations for each factor in the restricted model were similar to those of the main QCovid3 model apart from ethnic origin (for which no significant associations were seen) and conditions with too few events for analysis (sickle cell disease, severe combined immunodeficiency, HIV/AIDS). All associations with pre-existing health conditions reported are conditional on the other predictors in the model and do not necessarily have a causal interpretation.

QCovid3 model evaluation of performance

Table 2 shows the performance of the risk equations in the validation cohort. The QCovid3 algorithm for covid-19 related death explained 74.1% (95% confidence interval 71.1% to 77.0%) of the variation in time to covid-19 death, the Royston’s D statistic was 3.46 (3.19 to 3.73) and the Harrell’s C statistic was 92.5. Results were similar in men and women. Corresponding results restricted to the first vaccine dose were 71.3% (67.9% to 74.7%), 3.23 (2.96 to 3.50), and 93.6. The results restricting to the validation cohort after the second vaccine dose were similar but with wider confidence intervals owing to smaller numbers. The values for the R2, D, and C statistics were similar for the hospital admission equation. Supplementary table 2 shows the corresponding results for covid-19 death and hospital admission by age band where performance tended to be lower in the higher age bands.

Table 2

Performance of QCovid3 risk model in the validation cohort for covid-19 related death and hospital admission

View this table:

Figure 3 shows the calibration plot for covid-19 related deaths and figure 4 shows the corresponding results for covid-19 related hospital admission, both accounting for competing risks. Model calibration was assessed by comparing mean predicted risks with the observed cumulative incidence function by twentieths of predicted risk21; a model is well calibrated if predicted risks closely approximate the observed risks. These plots showed reasonable correspondence between observed and predicted cumulative incidences at 70 days of follow-up. However, numbers of events were small in several groups, and in the higher twentieths we saw slight under-prediction for covid-19 death (fig 3) and in twentieths 17-19 for the admissions outcome (fig 4). For example, in the top twentieth of predicted risks for covid-19 death, the observed cumulative incidence was 0.28% over 70 days (95% confidence interval 0.24% to 0.33%) and the mean predicted risk was 0.25%.

Fig 3
Fig 3

Calibration of the QCovid3 risk model to predict covid-19 related death after vaccination. Data source: QResearch England, 8 December 2020 to 15 June 2021, https://www.qresearch.org/

Fig 4
Fig 4

Calibration of the QCovid3 risk model to predict covid-19 related hospital admission after vaccination. Data source: QResearch England, 8 December 2020 to 15 June 2021, https://www.qresearch.org/

QCovid3 risk stratification

Table 3 shows the percentage of covid-19 related deaths identified by the QCovid3 mortality equation at different thresholds based on centiles of predicted absolute risk in the validation cohort, using the background SARS-CoV-2 infection rate associated with 14 days after vaccination. For example, it shows that 78.7% of deaths occurred in individuals in the top 5% for predicted absolute risk of covid-19 death (predicted absolute risks at 70 days above 0.06%). Individuals in the top 20% for predicted absolute risk of death accounted for 98.9% of deaths. Table 4 summarises the characteristics of individuals at the highest predicted absolute risk of covid-19 death (top 5%). Box 1 lists clinical examples of patients and their predicted covid-19 risks (https://bmjSept2021.qcovid.org).

Table 3

Sensitivity for covid-19 related death at 70 days’ follow-up in the validation cohort (consisting of 626 656 vaccinated patients with 174 covid-19 related deaths at different QCovid3 thresholds of absolute risk)

View this table:
Table 4

Characteristics of patients at highest risk of covid-19 related death (top 5%) from 14 days after covid-19 vaccination in the validation cohort using QCovid3 model. Data are numbers (percentages) unless stated otherwise

View this table:
Box 1

Clinical examples of patients and their predicted covid-19 risks over a 70 day period, based on QCovid3 risk algorithms (https://bmjSept2021.qcovid.org)

Example 1

72 year old white man with a first vaccine dose, atrial fibrillation, and body mass index of 30 (background daily infection rate of 22 positive reverse transcription polymerase chain reaction (RT-PCR) test results per 100 000 people) would have:

  • 0.04% risk of covid-19 related hospital admission over a 70 day period

  • 0.02% risk of covid-19 death over a 70 day period

  • 5.15% risk of covid-19 related death after a SARS-CoV-2 positive test result.

Example 2

62 year old Pakistani woman with two vaccine doses, chronic kidney disease stage 5 with transplantation, and body mass index of 24 (background daily infection rate of 20 positive RT-PCR test results per 100 000 people) would have:

  • 0.04% risk of covid-19 related hospital admission over a 70 day period

  • 0.003% risk of covid-19 death over a 70 day period

  • 0.10% risk of covid-19 related death after a SARS-CoV-2 positive test result.

Example 3

60 year old white man with a first vaccine dose, stroke, epilepsy, well controlled type 2 diabetes, Down’s syndrome, and body mass index of 41 (background daily infection rate of 60 positive RT-PCR test results per 100 000 people) would have:

  • 0.56% risk of covid-19 related hospital admission over a 70 day period

  • 0.46% risk of covid-19 death over a 70 day period

  • 24.3% risk of covid-19 related death after a SARS-CoV-2 positive test result.

Example 4

67 year-old Caribbean woman with a first vaccine dose, liver cirrhosis, and body mass index of 41 (background daily infection rate of 40 positive RT-PCR test results per 100 000 people) would have:

  • 0.08% risk of covid-19 related hospital admission over a 70 day period

  • 0.04% risk of covid-19 death over a 70 day period

  • 7.29% risk of covid-19 related death after a SARS-CoV-2 positive test result.

RETURN TO TEXT

QCovid2 model in unvaccinated patients

Supplementary figure 7 shows the adjusted hazard ratios for risk of covid-19 death for men and women in the comparison cohort of unvaccinated individuals using the QCovid2 model. Supplementary figure 8 shows the corresponding results for covid-19 hospital admission. These models are updated versions of the original QCovid model from our earlier publication.1 The new models include additional variables for inflammatory bowel disease, levels of diabetes control (according to measured HbA1c values), separate variables for sickle cell disease, HIV/AIDS, and severe combined immunodeficiency, and a refined definition of severe mental illness (which now includes only schizophrenia or bipolar disease but does not include moderate and severe depression).

The hazard ratios for QCovid2 were generally similar in magnitude and direction for the subset of variables included in our main model (QCovid3). However, some additional variables included in the QCovid2 model did not reach statistical significance or resulted in hazard ratios lower than 1.1, and hence were not included in the main QCovid3 model. These variables were type 1 diabetes, asthma, rare lung conditions, pulmonary fibrosis or pulmonary hypertension, congenital heart disease, cerebral palsy, inflammatory bowel disease, and severe mental illness (schizophrenia or bipolar disorder).

Supplementary table 3 shows similar performance statistics for discrimination and explained variation for the main QCovid2 model when compared with the original QCovid model evaluated in the validation cohort. Supplementary figures 9 and 10 show adjusted hazard ratios from a similar analysis for risk of covid-19 death and hospital admission in unvaccinated patients but restricted to individuals with a positive SARS-CoV-2 test result.

Discussion

Principal findings

We have identified a range of important clinical risk factors for severe covid-19 outcomes in people in the UK, 14 days or more after covid-19 vaccination (first or second dose) when some immunity is expected to have developed. We have used national linked datasets from general practice, national immunisation and SARS-CoV-2 testing, death registry, and hospital episode data for a population representative sample of more than 6.9 million adults. Risk ratios were highest for people with Down’s syndrome, kidney transplantation, sickle cell disease, care home residency, chemotherapy, recent bone marrow transplantation or solid organ transplantation ever, HIV/AIDS, dementia, Parkinson’s disease, neurological conditions, and liver cirrhosis. We also developed and evaluated novel clinical risk prediction models to estimate the absolute risks of covid-19 related hospital admission and mortality in the general population of vaccinated people as well as in the subset of people with positive SARS-CoV-2 test results. The risk models showed high levels of discrimination (C statistics ≥0.88 for the primary outcome, covid-19 death) and good calibration.

For many of the predictors included in the original QCovid model (contributing to risk prediction in an unvaccinated population during wave 1), the magnitude of the relative risks is broadly comparable in both QCovid2 (risk prediction in an unvaccinated population, wave 2) and QCovid3 (risk prediction in a vaccinated population, wave 2). Although these associations cannot be given a causal interpretation, individual characteristics such as age,23 obesity, pre-existing medical conditions, and socioeconomic disadvantage24 are known to affect immune competence25 and, at least for certain diseases, affect the response to some vaccines2627282930 or to immunosuppressive drugs.731 The associations with Down’s syndrome in all the models are likely to reflect increased susceptibility to infection and genetic predisposition.8 Compared with the white ethnic group, the Pakistani and Indian groups had up to twofold increased hazards of covid-19 death and hospital admission after vaccination in the full QCovid3 model. These ethnic disparities in covid-19 outcomes could represent residual differential exposure (eg, linked to behaviour, lifestyle, household size, and occupation) more than differential susceptibility mechanisms,32 although we also acknowledge that being vaccinated could change behaviour (and exposure) in some groups more than in others.

These risk models can be deployed in several health and care settings, either during the current phase of the pandemic or in subsequent waves of infection (with recalibration as required); however, absolute risk for individuals will always depend on disease prevalence and personal exposures. Uses of QCovid3 could include supporting targeted recruitment for clinical trials, prioritisation of vaccine boosters, future preventive treatments such as prophylactic passive monoclonal antibody protection, shielding, and discussions between individuals and clinicians on workplace or health risk mitigation (eg, through improved glycaemic control, weight reduction,33or general risk avoidance behaviours). Our QCovid3 model provides absolute risks conditional on patient characteristics, including whether they have received one or two doses of a covid-19 vaccine, and on the underlying prevailing infection levels. It also enables individuals to be ranked in terms of their risk. The deployment of drug and non-drug interventions to protect individuals with residual vulnerability after vaccination needs to be considered in the context of absolute risks of severe outcomes at the time of making predictions. Absolute risks are related to both the prevalence of SARS-CoV-2 infection in the population and the likelihood of SARS-CoV-2 exposure in a vaccinated adult population. Although these algorithms have been designed to inform UK health policy and interventions to manage covid-19 related risks, they also have international potential, subject to local validation. Previous similar risk prediction models have been validated internationally and shown to have good performance outside of the UK.3435

Strengths and limitations of this study

Our study has some major strengths but also some important limitations, which include specific issues related to covid-19 along with factors similar to those for several other widely used clinical risk prediction algorithms developed using the QResearch database.363738 Key strengths included the use of large, validated, representative, population based contemporaneous data sources that have been used to develop other widely used risk prediction tools3637; the wealth of candidate risk predictors; the prospective recording of outcomes and their ascertainment using linkage of multiple national databases; lack of selection, recall, and respondent biases; and robust statistical analysis. We have used non-linear terms to model body mass index, age, and background SARS-CoV-2 infection rates. The inclusion of infection rates is a substantial improvement compared with the original QCovid model, because it enables risks to be updated according to the background infection, which is important given the nature of pandemic waves. Our analysis has been able to separately quantify the risk of severe outcomes among those people with a positive SARS-CoV-2 test result, which was not possible in the original QCovid model1 owing to a lack of testing data. Therefore, the analysis could be used at the point of testing to identify those who might benefit from additional interventions such as monoclonal antibodies once these treatments become available.

Limitations included a relatively short duration of follow-up, a partially vaccinated population, and small numbers of events in some subgroups—which are inevitable consequences of undertaking an analysis during rapid deployment of a national vaccination programme. Our analysis incorporated information on whether an individual had received one or two vaccination doses in our prediction models. We saw relatively few deaths in individuals who had received the second dose of the vaccine (4% of all covid-19 related deaths); therefore, most information about associations between predictors and mortality came from individuals who had received only one dose. Results from individuals who have received the second vaccination dose are likely to be most relevant for UK adults as full vaccine coverage increases. Although we examined for interactions, our study might have lacked power to detect whether certain associations differed according to whether one or two doses had been received. Our models also incorporated information on prevailing positive SARS-CoV-2 infection rates, as a proxy for a person’s risk of covid-19 infection at the start of the follow-up period.

We did not include information on different variants that emerged during the study period owing to incomplete data, particularly in those patients admitted to hospital.39 While we accounted for many risk factors for covid-19 mortality, some risks could remain, such as those conferred by rare medical conditions or other factors associated with exposure (eg, occupation) that are poorly recorded in general practice or hospital records and that might be being proxied to some extent by the covariates included. We did not distinguish vaccination type because this study was not designed to compare vaccine effectiveness. Younger patients without underlying health conditions had limited data, because the vaccination programme in England prioritised elderly patients and those people at highest risk. Furthermore, those patients who had two vaccines early in the pandemic were judged to be at highest risk of infection or severe outcomes.

Although we have reported a validation using practices from QResearch, these practices were separate to those used to develop the model. Previously we have used this approach to develop and validate other widely used prediction models. When these models have been validated on different clinical computer systems, the results have been similar.404142 Work is already underway to evaluate the new models in external datasets (such as the English national dataset hosted by the Office for National Statistics) including data from other general practice computer systems that have not been used to derive the algorithm. These data offer a fully independent dataset including data from general practice computer systems not included in the derivation of the dataset. They also offer a larger sample size for validation because clinical and demographic subgroups will have more events. Work is also underway to consider integration of this new algorithm within NHS clinical software systems.

Policy implications and conclusions

This study presents a robust risk prediction model (QCovid3) that can be used to stratify risk populations to identify those who are at highest risk of severe covid-19 outcomes despite covid-19 vaccination, and who might therefore benefit from further interventions to reduce risk or boost immunity once these become available. The model can be used in conjunction with QCovid2, which updates and replaces the original algorithm (QCovid1) and is designed for use in unvaccinated patients. We anticipate that these algorithms will be updated as the vaccination programme progresses and is extended to younger age groups, as understanding of covid-19 increases, as more post-vaccination follow-up data become available, as new variants of concern emerge, and in response to new policy interventions.

What is already known on this topic

  • The original QCovid tool for predicting risk of covid-19 related death or hospital admission based on individual characteristics was used in England to identify patients at high risk of severe covid-19 outcomes

  • Identification of these high risk patients added an additional 1.5 million people to the national shielded patient list in February 2021

  • On a UK basis, these patients would be prioritised for vaccination, if they had not already been offered the vaccine on account of their age or occupation

What this study adds

  • Commissioned by the Chief Medical Officer for England on behalf of the UK government, two new risk prediction algorithms have been derived and validated to estimate the risk of covid-19 related mortality and hospital admission in UK adults, 14 days or more after vaccination when some immunity is expected to have developed

  • Several clinical risk factors for severe covid-19 outcomes despite vaccination have been identified: Down’s syndrome, kidney transplantation, sickle cell disease, care home residency, chemotherapy, recent bone marrow transplantation or a solid organ transplantation ever, HIV/AIDS, dementia, Parkinson’s disease, neurological conditions, and liver cirrhosis

  • The QCovid3 risk algorithms (https://bmjSept2021.qcovid.org) showed high levels of discrimination for identifying adults at highest risk of covid-19 related death and hospital admission after vaccination; these risk stratification tools can help support public health policy and prioritise patients for targeted, early interventions

Ethics statements

Ethical approval

The QResearch ethics approval was provided on 8 June 2020 by the East Midlands-Derby research ethics committee (reference 18/EM/0400).

Data availability statement

To guarantee the confidentiality of personal and health information, only the authors have had access to the data during the study in accordance with the relevant licence agreements. Access to the QResearch data is according to the information on the QResearch website (https://www.qresearch.org/). The full model, model coefficients, functional form and cumulative incidence function are published on the qcovid.org website.

Acknowledgments

We acknowledge the contribution of Egton Medical Information Systems practices who contribute to QResearch and EMIS Health and the Universities of Nottingham and Oxford for expertise in establishing, developing, or supporting the QResearch database. This project involves data derived from patient level information collected by the NHS, as part of the care and support of patients with cancer. The data are collated, maintained, and quality assured by the National Cancer Registration and Analysis Service, which is part of Public Health England. The Hospital Episode Statistics data, SARS-Cov-2 results, and civil registration mortality data are used by permission from NHS Digital who retain the copyright in those data. NHS Digital and Public Health England bear no responsibility for the analysis or interpretation of the data.

Footnotes

  • Contributors: Study conceptualisation was led by JH-C, CACC, and JSN-V-T. JH-C specified the data, organised data approvals, and data linkage. JH-C, CACC, RHK, and KD-O designed the statistical analysis plan. JH-C and CACC undertook the analyses. JH-C wrote the first draft of the paper and developed the software for the web calculator. The Office of the Chief Medical Officer contributed to the development of the study question and facilitated access to relevant national datasets, contributed to interpretation of data, and contributed to the drafting of the report. All authors contributed to the interpretation of the results and revision of the manuscript and approved the final version of the manuscript. JH-C had full access to all data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. JH-C is the guarantor for the study. The corresponding author attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted.

  • Funding: This study was funded by the National Institute for Health Research (NIHR) following a commission by the Chief Medical Officer for England. The researchers are independent from the NIHR. QResearch was supported by funds from the John Fell Oxford University Press Research Fund, grants from Cancer Research UK (grant C5255/A18085), through the Cancer Research UK Oxford Centre, grants from the Oxford Wellcome Institutional Strategic Support Fund (204826/Z/16/Z), during the conduct of the study. The funders had no role in considering the study design or in the collection, analysis, interpretation of data, writing of the report, or decision to submit the article for publication.

  • Competing interests: All authors have completed the ICMJE uniform disclosure form at www.icmje.org/disclosure-of-interest/ and declare: support from the NIHR for the submitted work; JH-C reports grants from NIHR Biomedical Research Centre, Oxford, grants from John Fell Oxford University Press Research Fund, grants from Cancer Research UK, through the Cancer Research UK Oxford Centre, grants from the Oxford Wellcome Institutional Strategic Support Fund and other research councils, during the conduct of the study; JH-C is an unpaid director of QResearch, a not-for-profit organisation that is a partnership between the University of Oxford and EMIS Health, which supplied the QResearch database used for this work; JH-C is a founder and shareholder of ClinRisk and was its medical director until 31 May 2019 (ClinRisk produces open and closed source software to implement clinical risk algorithms (outside this work) into clinical computer systems); JH-C is chair of the NERVTAG risk stratification subgroup and a member of SAGE covid-19 groups and the NHS group advising on prioritisation of use of monoclonal antibodies in covid-19 infection; CACC reports receiving personal fees from ClinRisk outside this work, and is a member of the NERVTAG risk stratification subgroup; KK is supported by the NIHR Applied Research Collaboration-East Midlands and the Leicester BRC and is a member of SAGE; RHK was supported by a UKRI Future Leaders Fellowship (MR/S017968/1); KD-O was supported by a grant from the Alan Turing Institute Health Programme (EP/T001569/1); AS is a member of the Scottish Government Chief Medical Officer’s covid-19 advisory group and a member of AstraZeneca’s thrombotic thrombocytopenic advisory group (both roles are unremunerated); RAL is a member of the Welsh government covid-19 technical advisory group (unrenumerated); MGS reports grants from Department of Health and Social Care NIHR UK, grants from the UK Medical Research Council, grants from Health Protection Research Unit in Emerging and Zoonotic Infections, University of Liverpool, during the conduct of the study, and other funds from Integrum Scientific, Greensboro, NC, USA, outside the submitted work; MGS is a member of NERVTAG and attends SAGE covid-19; AH is a member of NERVTAG and the NHS group advising on prioritisation of use of monoclonal antibodies in covid-19 infection; JV is national clinical director for diabetes and obesity at NHS England and Improvement; FK is a member of the Northern Ireland Chief Medical Officer’s pandemic modelling group and strategic intelligence group; JSN-V-T is seconded to the Department of Health and Social Care, England. The views expressed in this manuscript are those of the authors and not necessarily those of Department of Health and Social Care or the UK government.

  • The lead author affirms that the manuscript is an honest, accurate, and transparent account of the study being reported; that no important aspects of the study have been omitted; and that any discrepancies from the study as planned (and, if relevant, registered) have been explained.

  • Dissemination to participants and related patient and public communities: Patients will be invited to advise on disseminating the results.

  • Provenance and peer review: Not commissioned; externally peer reviewed.

http://creativecommons.org/licenses/by/4.0/

This is an Open Access article distributed in accordance with the terms of the Creative Commons Attribution (CC BY 4.0) license, which permits others to distribute, remix, adapt and build upon this work, for commercial use, provided the original work is properly cited. See: http://creativecommons.org/licenses/by/4.0/.

References