Type 2 Diabetes and COVID-19-Related Mortality in the Critical Care Setting: A National Cohort Study in England, March-July 2020.

OBJECTIVE: To describe the relationship between type 2 diabetes and all-cause mortality among adults with coronavirus disease 2019 (COVID-19) in the critical care setting. RESEARCH DESIGN AND METHODS: This was a nationwide retrospective cohort study in people admitted to hospital in England with COVID-19 requiring admission to a high dependency unit (HDU) or intensive care unit (ICU) between 1 March 2020 and 27 July 2020. Cox proportional hazards models were used to estimate 30-day in-hospital all-cause mortality associated with type 2 diabetes, with adjustment for age, sex, ethnicity, obesity, and other major comorbidities (chronic respiratory disease, asthma, chronic heart disease, hypertension, immunosuppression, chronic neurological disease, chronic renal disease, and chronic liver disease).
RESULTS: A total of 19,256 COVID-19-related HDU and ICU admissions were included in the primary analysis, including 13,809 HDU (mean age 70 years) and 5,447 ICU (mean age 58 years) admissions. Of those admitted, 3,524 (18.3%) had type 2 diabetes and 5,077 (26.4%) died during the study period. Patients with type 2 diabetes were at increased risk of death (adjusted hazard ratio [aHR] 1.23 [95% CI 1.14, 1.32]), and this result was consistent in HDU and ICU subsets. The relative mortality risk associated with type 2 diabetes decreased with higher age (age 18-49 years aHR 1.50 [95% CI 1.05, 2.15], age 50-64 years 1.29 [1.10, 1.51], and age ≥65 years 1.18 [1.09, 1.29]; P value for age-type 2 diabetes interaction = 0.002).
CONCLUSIONS: Type 2 diabetes may be an independent prognostic factor for survival in people with severe COVID-19 requiring critical care treatment, and in this setting the risk increase associated with type 2 diabetes is greatest in younger people.

Introduction

In early 2020, coronavirus disease 2019 (COVID-19), caused by the new highly infectious organism now known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was formally declared a public health emergency of international concern by the World Health Organization (1). With a global prevalence of 9.3% (2), it is unsurprising that diabetes is one of the most common comorbidities seen in people with COVID-19, second only to hypertension (3). People with type 2 diabetes are known to not only be more susceptible to infections in general (4) but also require hospitalization more often (5), resulting in an overall worse prognosis (6). As the COVID-19 pandemic has progressed, a similar pattern of a worse prognosis for people with diabetes who contract COVID-19 has been reported in population-based studies (7,8), with an increased risk of intensive care admission (9–11) and worse mortality outcome in some but not all hospital cohorts (12–14).

For people admitted to intensive care, studies published early in the course of the global pandemic suggested diabetes is an important prognostic factor even in severe COVID-19 (New York, n = 373, 40.2% with diabetes, unadjusted odds ratio [OR] for mortality in those with diabetes 2.31 [95% CI 1.34, 4.00], and Wuhan, China, n = 193, 24.9% with diabetes, adjusted hazard ratio [aHR] 1.53 [95% CI 1.02, 2.30]) (11,15). Two larger, more recent critical care studies have reported smaller effect sizes (Italy, n = 3,988, 12.9% with diabetes, adjusted OR 1.18 [95% CI 1.01, 1.39], and U.S., n = 2,215, 38.9% with diabetes, adjusted OR 1.14 [95% CI 0.91, 1.43]) (16,17). Heterogeneity across these studies likely reflects a combination of local variation in prevalence of diabetes, variation in the characteristics of people with diabetes, and the limited geo-temporal snapshots each study represents. Data on the additional mortality associated with diabetes in critical care across the entire first wave of COVID-19 are lacking in these studies, and potential heterogeneity in the risk associated with diabetes by age and other clinical features has yet to be fully explored.

We aimed to establish the clinical utility of type 2 diabetes as a prognostic factor for survival in people admitted to critical care with COVID-19, using national data over the course of the entire first wave of COVID-19 from hospitals in England.

Research Design and Methods

Data Source

Data were accessed from the COVID-19 Hospitalisation in England Surveillance System (CHESS), established by Public Health England (PHE) (18). Reporting to CHESS is daily via Web tool and is mandatory for hospitals in England. The data specification comprises epidemiological information for people with proven or a high likelihood of COVID-19 requiring hospitalization and high dependency unit (HDU) or intensive care unit (ICU) management, including demographics, comorbidities including diabetes and obesity, and outcome. Confirmation of infection was by RT-PCR of nasopharyngeal and/or oropharyngeal swabs, as it was the only available method during the study period. National testing policy changed throughout the study period; however, all people with a suspicion of COVID-19 and who were admitted to hospital were tested (19), and people were potentially retested multiple times if the initial results were negative but clinical suspicion remained high. At the time of data extraction (27 July 2020), data on 23,935 cases across 108 hospitals had been submitted to CHESS.

Study Population

Inclusion and Exclusion Criteria

The study population comprised people in CHESS aged 18–99 years inclusive who required HDU or ICU admission with confirmed or clinically diagnosed COVID-19 (91.4% with PCR swab positive result at some point in their admission) from 1 March 2020 to 27 July 2020 inclusive and were not pregnant (n = 22,082). Pregnancy is associated with altered metabolic physiology; due to the limited information about this group, an a priori decision was made to exclude pregnant patients from the analysis.

Cohorts

The primary cohort for analysis was all patients admitted to either HDU or ICU at any point during the study period (HDU/ICU cohort). We separately examined the subsets of all patients admitted to HDU but not ICU during the study period (HDU cohort) and all patients admitted to ICU during the study period (ICU cohort). People with type 1 diabetes were excluded from these primary cohorts and are analyzed separately.

Recorded Clinical Features

Available clinical features comprised demographic characteristics (age, sex, ethnicity), obesity status, and comorbidities, as defined in the CHESS data specification (18), recorded during critical care admission. We grouped ethnicity into major U.K. ethnic census groups: White, Black, Asian, mixed, other, or not recorded. We categorized obesity status as nonobese (BMI <30 kg/m2), obese (BMI ≥30 kg/m2), or not recorded, as derived from recorded BMI category (coded in CHESS as <18.5, 18.5–24.9, 25–29.9, 30–39.9, or >39.9 kg/m2 or no, missing, or unknown) where available or recorded clinical obesity (coded as no, borderline, yes, missing, or unknown) if BMI was not available. Coded comorbidities were diabetes (separated by type 1 and 2), chronic respiratory disease, asthma requiring medication, chronic/congenital heart disease, hypertension, immunosuppression (due to disease or treatment), chronic neurological disease, chronic renal disease, and chronic liver disease, all coded as no, yes, missing, or unknown.

Missing Data

We found variation in reporting of information on comorbidity at the hospital level, with a limited number of hospitals reporting very little information on comorbidity. For our primary analysis, we excluded all patients from 15 hospitals where comorbidity recording was missing for >50% of people (n = 2,151) (Supplementary Flowchart 1). After filtering the CHESS cohort by missingness, we defined binary variables for each comorbidity (yes vs. no/missing/unknown), assuming that, in this mostly complete subset, missing or unknown coded data represented the absence of that condition. We tested the validity of this assumption in the sensitivity analysis, described below.

Outcome

The outcome was in-hospital all-cause mortality. People were followed up from admission to hospital until the earliest of the following: death, hospital transfer, date of last follow-up in CHESS, or a length of stay of 30 days. To account for discharge from hospital as a competing risk, we assumed that people who were discharged were no longer at risk of death, and so these individuals were not censored at the day of discharge and instead had their follow-up time set to 30 days (12).

Statistical Analysis

Primary Analysis

Recorded clinical features were described, overall and by recorded type 2 diabetes status. We estimated the risk of 30-day mortality for each cohort by type 2 diabetes status using survival analysis with days at risk as the timescale. Thirty-day mortality covered 98% of the deaths observed over the maximum cohort follow-up. Kaplan-Meier survival functions and Cox proportional hazards models were used to evaluate the association between type 2 diabetes and mortality, with proportional hazards assumptions tested through examination of the Schoenfeld residuals. All models included type 2 diabetes status and were sequentially increased in complexity with adjustment firstly for age (with exploration of age as linear, category, and a three-knot nonlinear restricted cubic spline), before addition of sex, ethnicity, obesity, and comorbidities. The primary, fully adjusted, model comprised all these features. For the ICU cohort, we started follow-up on the date of ICU admission and adjusted for days from hospitalization to ICU admission.

Subgroup Analysis

Analysis of the primary outcome for the HDU/ICU cohort was repeated in subgroups defined by major clinical characteristics; age (18–49, 50–64, ≥65 years), sex, and ethnicity; the presence or absence of obesity, and the presence or absence of at least one other comorbidity. For exploration of temporal trends, subgroups were defined by calendar month of hospital admission (March, April, and May–July 2020).

Sensitivity Analysis

To explore the sensitivity to the decision to filter patients by degree of missing comorbidity data at hospital level, we described recorded characteristics of patients in included and excluded hospitals and assessed whether the mortality associated with type 2 diabetes differed in included and excluded hospitals. To explore sensitivity of our findings to potential temporal changes in mortality over the study period, calendar week of hospital admission was included as an additional linear continuous covariate in multivariable analysis.

We conducted a propensity score–matched analysis to evaluate whether the presence of type 2 diabetes increased the risk of death compared with risk in patients with the same clinical characteristics but without type 2 diabetes. The full set of recorded clinical features previously described was considered as potential confounders and used for matching. Missing information for ethnicity and obesity was coded as missing categories, providing a full patient set (20). Pairs were matched without replacement on the logit of the propensity score, and a caliper size of 0.05 was applied for all matched pairs. Separate propensity score models were fitted for each cohort. Cox proportional hazards models were stratified by matched set to account for the matched nature of the propensity score–defined cohorts.

We repeated the primary analysis with the following changes to formally assess the impact of the assumptions made: 1) excluding transferred patients and patients who remained in hospital, 2) including only patients with PCR swab positive results during admission, 3) including only hospitals with at least 75% complete comorbidity data, 4) including all patients in all hospitals and assuming that missing/unknown represented the absence of a comorbidity, 5) extending the Cox proportional hazards model to incorporate hospital-level random effects to account for within-hospital homogeneity in outcomes (frailty model with a random effect per hospital trust) (21), and 6) using a Fine and Gray model to check the validity of our approach by accounting for discharge as a competing risk.

Evaluation of Missing Data

To evaluate the appropriateness of equating diabetes status coded as missing in CHESS to the absence of diabetes, we used multiple imputation by chain equations (five imputations) to fill missing values for diabetes status and all other missing values for ethnicity, obesity, and other comorbidities. This approach assumes that missing values are missing at random (22). We then repeated the full multivariable analysis for each cohort using the imputation model, to determine whether this model gave different results to our primary analysis.

COVID-19 Mortality Risk in People With Type 1 Diabetes

Mortality in people with recorded type 1 diabetes was compared with that in people with type 2 diabetes with adjusted Cox regression, as described above. A propensity-matched analysis was also used, as previously described, to compare the mortality risk for people with type 1 diabetes with that for similar patients without type 1 diabetes. Analyses were limited to the combined HDU/ICU cohort due to sample size.

Computational Resources

Analyses were conducted with R (version 3.6.2) (23), including use of the packages survival, rms, coxme, paf, mlr3, mlr3tuning, mlr3proba, xgboost, and mice.

Ethics and Governance

The study was reviewed and approved by the Warwick Biomedical & Scientific Research Ethics Committee (BSREC 119/19-20), and sponsorship for this study was provided by University of Warwick (SOC.28/19-20).

Data and Resource Availability

Data cannot be shared publicly, as it was collected by PHE as part of their statutory responsibilities, which allows them to process patient confidential data without explicit patient consent. Data used in this study were made available through an agreement between the University of Warwick and PHE. Individual requests for access to CHESS data are considered directly by PHE (contact via covid19surv@phe.gov.uk).

Results

Clinical Characteristics

A total of 19,256 patients were included in the primary cohort admitted to HDU (n = 13,809) or ICU (n = 5,447) (Supplementary Flowchart 1 and Supplementary Table 1). Mean age was 67 years, 60.1% of patients were male, 27.0% of patients were of non-White ethnicity, and 18.3% had type 2 diabetes. Supplementary Table 1 shows the overall characteristics of patients in the separate HDU and ICU cohorts. ICU patients were younger and more commonly male and of non-White ethnicity (mean age 58 years, 70.6% male, 37.1% non-White) than HDU-only admitted patients (mean age 70 years, 55.9% male, 23.1% non-White).

In the primary HDU/ICU cohort, patients with type 2 diabetes (n = 3,524 [18.3%]) were of similar age to those without type 2 diabetes (mean 67 vs. 66 years, respectively) and were more commonly recorded as being of non-White ethnicity (41.5% vs. 23.5%) and obese (45.1% vs. 31.4%) (Table 1). All comorbidities were more common in patients with type 2 diabetes.

Table 1

Recorded characteristics of HDU/ICU cohort, by recorded type 2 diabetes status

	No type 2 diabetes (n = 15,732)	Type 2 diabetes (n = 3,524)
Age (years)	66 (17.4)	67 (14.1)
Age-group (years)
18–24	164 (1.0)	8 (0.2)
25–34	612 (3.9)	40 (1.1)
35–44	1,010 (6.4)	142 (4.0)
45–54	2,039 (13.0)	429 (12.2)
55–64	2,814 (17.9)	848 (24.1)
65–74	2,978 (18.9)	872 (24.7)
75–84	3,334 (21.2)	718 (20.4)
≥85	2,781 (17.7)	467 (13.3)
Sex
Female	6,440 (40.9)	1,243 (35.3)
Male	9,292 (59.1)	2,281 (64.7)
Ethnicity*
White	9,624 (76.5)	1,815 (58.5)
Asian	1,356 (10.8)	680 (21.9)
Black	737 (5.9)	360 (11.6)
Mixed	163 (1.3)	48 (1.5)
Other	694 (5.5)	201 (6.5)
Obesity**
Nonobese	6,034 (68.6)	1,509 (54.9)
Obese	2,764 (31.4)	1,240 (45.1)
Comorbidity
Any comorbidity	6,228 (39.6)	2,945 (83.6)
Chronic respiratory disease	1,250 (7.9)	423 (12.0)
Asthma	1,200 (7.6)	428 (12.1)
Hypertension	3,439 (21.9)	2,218 (62.9)
Chronic heart disease	1,828 (11.6)	834 (23.7)
Chronic renal disease	996 (6.3)	831 (23.6)
Chronic liver disease	252 (1.6)	138 (3.9)
Chronic neurological disease	1,098 (7.0)	310 (8.8)
Immunosuppressive disease	374 (2.4)	94 (2.7)

Data are N (%), except for Age (years), which is mean (SD). *Ethnicity not recorded for 3,578 patients.

Obesity not recorded for 7,709 patients.

Type 2 Diabetes Is an Independent Prognostic Factor for Mortality in Severe COVID-19

A total of 5,077 (26.4%) patients in the primary HDU/ICU cohort died within the 30-day study period. Mean (SD) follow-up was 21.2 (10.9) days in people with type 2 diabetes and 21.2 (19.4) days in those without.

Risk of mortality was consistently higher in people with type 2 diabetes in all analyses performed. The overall unadjusted cumulative incidence of 30-day mortality was higher in people with diabetes (34.7% [95% CI 33.1, 36.3]) than in people without diabetes (25.5% [24.8, 26.2]) (Fig. 1). Mortality did not differ between people with diabetes recorded as missing/unknown diabetes status and those recorded as no diabetes (Supplementary Fig. 1), supporting the assumption that missing coding represents the absence of diabetes. Unadjusted hazard ratios and aHRs showed higher mortality risk in people with diabetes, with a 23% increase in mortality risk observed in the fully adjusted model (Table 2). Full multivariable model outputs are provided in Supplementary Table 2. A 23% increase in mortality risk for diabetes was also observed in the propensity score–matched analysis (Table 2).

Figure 1

Kaplan-Meier plots for in-hospital COVID-19 death in 19,256 patients admitted to critical care (HDU or ICU) in CHESS by time since hospital admission. Plots show the proportion of individuals at risk who were still alive at regular intervals up to 30 days from admission, stratified by the presence of type 2 diabetes. People discharged from hospital prior to 30 days were assumed to survive to 30 days and are included in the number at risk until 30 days, in keeping with the standard practice for time-to-event modeling in analysis of critical care patients.

Table 2

Hazard ratios for diabetes for the primary outcome of 30 day in-hospital mortality for people with type 2 diabetes admitted to critical care with COVID-19

Adjustment	Cohort 1: all patients (n = 19,256)	Cohort 2: HDU-only patients (n = 13,809)	Cohort 3: ICU patients (n = 5,447)
Unadjusted	1.44 (1.35, 1.54), P < 0.001	1.38 (1.26, 1.51), P < 0.001	1.34 (1.22, 1.48), P < 0.001
Age	1.42 (1.33, 1.51), P < 0.001	1.31 (1.20, 1.44), P < 0.001	1.28 (1.16, 1.41), P < 0.001
Age and sex	1.40 (1.31, 1.49), P < 0.001	1.30 (1.19, 1.42), P < 0.001	1.28 (1.16, 1.41), P < 0.001
Age, sex, ethnicity	1.36 (1.27, 1.45), P < 0.001	1.29 (1.18, 1.41), P < 0.001	1.23 (1.11, 1.36), P < 0.001
Age, sex, ethnicity, obesity	1.32 (1.23, 1.41), P < 0.001	1.27 (1.16, 1.39), P < 0.001	1.29 (1.17, 1.43), P < 0.001
Full covariate set adjustment model (age, sex, ethnicity, obesity, comorbidity*)	1.23 (1.14, 1.32), P < 0.001	1.19 (1.08, 1.31), P < 0.001	1.24 (1.11, 1.38), P < 0.001
Propensity score–matched model, adjustment for full covariate set	1.25 (1.14, 1.36), P < 0.001 (n = 3,289 with diabetes matched)	1.17 (1.04, 1.32), P = 0.01 (n = 2,007 with diabetes matched)	1.19 (1.05, 1.35), P = 0.009 (n = 1,244 with diabetes matched)
Competing risk model, adjustment for full covariate set (Fine and Gray)**	1.27 (1.18, 1.36), P < 0.001	1.24 (1.13, 1.36), P < 0.001	1.24 (1.12, 1.38), P < 0.001

Age fitted as a restricted cubic spline with three knots in all models.

Chronic respiratory disease, asthma, chronic heart disease, hypertension, immunosuppression, chronic neurological disease, chronic renal disease, chronic liver disease.

Fine and Gray model subdistribution hazard ratios taking into account the competing event of being discharged alive.

The mortality increase in people with type 2 diabetes was consistent in the HDU cohort (type 2 diabetes 30-day mortality 29.2% [95% CI 27.2, 31.1] vs. no diabetes 22.4% [21.6, 23.2]) and ICU cohort (30-day mortality 45.2% [42.4, 47.9] vs. 36.3% [34.7, 37.8]) (Table 2 and Supplementary Fig. 2).

The mortality risk for type 1 diabetes (n = 203) was similar to that observed for type 2 diabetes (aHR 0.98 [95% CI 0.75, 1.28], P = 0.89), and in propensity score–matched analysis the hazard ratio for type 1 diabetes was similar to that observed for matched control subjects without diabetes (aHR 1.25 [95% CI 0.86, 1.84], P = 0.25).

Subgroup Analysis Suggests Lower Excess Mortality Risk Associated with Type 2 Diabetes at Older Ages and in Males Compared with Females

Results of subgroup analyses for the primary outcome are shown in Fig. 2. The association between type 2 diabetes and mortality was attenuated by age at hospital admission (age 18–49 years aHR 1.50 [95% CI 1.05, 2.15]), age 50–64 years 1.29 [1.10, 1.51], and age ≥65 years 1.18 [1.09; 1.29], P = 0.002 for age category–type 2 diabetes interaction), and there was a numerically lower excess mortality risk in males (1.15 [1.05, 1.26]) than females (1.36 [1.20, 1.54], P = 0.06 for sex-diabetes interaction). There was no evidence of a difference in the mortality associated with type 2 diabetes in subgroups defined by ethnicity (P = 0.08), obesity (P = 0.28), or presence or absence of other recorded comorbidity (P = 0.41). There was no evidence of a temporal trend in the mortality associated with type 2 diabetes by calendar month (hazard ratio range 1.13–1.34 over March–July) (Supplementary Table 3) or with adjustment for calendar week of hospital admission as an additional covariate (aHR for diabetes in the primary HDU/ICU cohort 1.21 [95% CI 1.13, 1.31]) (Supplementary Table 3).

Figure 2

Major subgroup analysis of the primary outcome of 30-day in-hospital mortality. Hazard ratios represent the mortality risk associated with the presence of type 2 diabetes in each subgroup. Bars represent 95% CIs. *Hazard ratios for other subgroups: mixed ethnicity 1.91 (95% CI 1.03, 3.56), other ethnicity 0.96 (0.67, 1.39), missing ethnicity 1.34 (1.11, 1.62), missing obesity 1.28 (1.08, 1.46).

Sensitivity and Multiple Imputation Analyses Were Consistent With the Primary Analysis

Results were consistent in all sensitivity analyses (Supplementary Table 3), including inclusion of only patients with a positive PCR swab result during admission; repetition of the primary analysis with hospital added as a random effect in the multiple imputation model; and, despite differences in clinical characteristics (Supplementary Table 4), inclusion of patients from all hospitals currently reporting to CHESS irrespective of completeness of comorbidity coding.

Conclusions

Our analysis of >19,000 patients admitted to critical care over the entire first wave of COVID-19 in England demonstrates that type 2 diabetes is associated with ∼20% increase in mortality risk in people with severe COVID-19, independent of age, sex, ethnicity, obesity, and other major comorbidity. Risk is similar in both the HDU and ICU cohorts, despite the fact that in England these groups differed markedly, with those admitted to ICU being younger, more commonly male, and more commonly of non-White ethnicity. Additional mortality risk appears similar in type 1 diabetes compared with type 2 diabetes in this setting.

Our national-level study provides evidence for a relatively small, but significant, effect of diabetes on mortality risk in one of the largest samples of COVID-19–associated critical care admissions to date. The association between diabetes and COVID-19 that mortality we observed is substantially smaller than those reported in studies conducted early in the pandemic (11,15). Notably, we extend previous contributions by evaluating both HDU and ICU admissions, which we demonstrate to be very different populations in terms of clinical characteristics. In England, these differences between HDU and ICU populations likely result from individualized care strategies based on the patient’s best interests (which one UK-based center implemented for 61% of 429 consecutive admissions) (24) rather than sustained saturation of ICUs, for which there was limited evidence (25).

Given the large sample available, we were able to evaluate heterogeneity in mortality risk across subgroups defined by clinical features and demonstrate that, in relative terms, the additional mortality risk associated with type 2 diabetes attenuates markedly in older people. In contrast to age, we found no evidence of heterogeneity in the mortality risk associated with diabetes in subgroups defined by the presence or absence of obesity, comorbidity, and ethnicity. There was weak evidence of heterogeneity by sex, with a numerically higher excess risk associated with diabetes in females than males (P = 0.06), which may warrant future investigation. Importantly, we demonstrate a consistent strength of association between type 2 diabetes and mortality accounting for time trends in in-hospital mortality, and accounting for geographical clustering of hospitals, which is notable, as recent studies have demonstrated substantial variation in COVID-19–specific critical care outcomes between institutions (26), as well as a marked reduction in critical care mortality over the course of the pandemic (27).

Although direct evaluation of patient characteristics other than diabetes was not the focus of this study, the validity of our results is bolstered by demographic congruence with COVID-19 data from the U.K.’s Intensive Care National Audit & Research Centre (ICNARC) (28), where diabetes status is not available. Notably, both studies found similar effects for males (i.e., no difference) and people of Asian ethnicity (i.e., increased risk) in critical care settings (Supplementary Table 2). The exception is that we did not find an association between obesity and mortality, in contrast to ICNARC, which reported an association of small magnitude for BMI but without adjustment for diabetes or hyperglycemia.

Previous studies of other viral respiratory infections suggest that findings of population-level studies should be extrapolated to the critical care setting with caution. For influenza, national surveys suggested increased influenza-specific mortality in people with diabetes (29); however, large multistate European studies of ICU admissions for influenza showed no difference in mortality by diabetes status (30). Our results, contrasted with the numerically much greater associations for diabetes reported in U.K. population-based studies of COVID-19 mortality, mirror this pattern of substantially reduced effect size (7,8). Such setting-specific differences likely reflect different frames of reference, wherein each study is capturing a different portion of a complicated pathway including the risk of being infected and then accessing services, subsequently being admitted, hospital treatments, and finally, dying in hospital.

To our knowledge, CHESS is currently the largest critical care database of people with COVID-19 and coded diabetes status, providing near-complete capture of severe COVID-19 cases across England during the first wave of the pandemic (18). Our use of a methodologically robust modeling framework, with adjustment for potential confounders and multiple sensitivity analyses to evaluate the potential influence of missing data, lends additional credibility to our mortality estimates. Despite this, we cannot rule out unmeasured confounding as an explanation of our findings. Any interpretation of our results should be bounded by the knowledge that in the primary analysis we adjusted for variables along the causal pathway (e.g., cardiovascular comorbidities), which precludes any claim that associations relate to the causal role of type 2 diabetes in COVID-19 mortality (31). However, the consistent size of association for type 2 diabetes observed across models with sequential adjustment of clinical features supports the interpretation that type 2 diabetes is an important prognostic factor for critical care patients.

A limitation of the CHESS critical care data set is the lack of validated case definitions for recorded comorbidities, and it is possible that misclassification of diabetes status may have attenuated the mortality risk estimates observed in this study. A further important limitation is that the lack of standardized measurement of height and weight meant we had to adjust risk models for obesity status (BMI categorized as < or >30 kg/m2) rather than with more granular categorization of BMI. A recent study of 6,916 U.S. patients diagnosed with COVID-19 suggests an independent J-shaped association between BMI category and mortality, with patients of BMI <18.5 kg/m2 and ≥40 kg/m2 at higher risk of mortality compared with those of normal BMI (18.5–24 kg/m2) (32), and an analysis of 1,687 hospitalized adults in New York observed a similar J-shaped pattern between BMI and mortality, although with overweight patients at lowest risk (33). These analyses raise the possibility of residual confounding, due to the limited recording of BMI in CHESS, being a potential explanation of our findings. We lacked data to assess the impact of potentially modifiable diabetes-specific risk factors, in particular, hyperglycemia prior to, at, or during, hospital admission (7,14,34–37), as well as possible modulatory effects of specific antihyperglycemic medications (38–40). Clinical information on duration of diabetes and the presence of diabetes complications was similarly not available. Robustly evaluating the added value of such potential diabetes-specific prognostic factors is an important area for ongoing (35,37), and future, research to inform understanding of the mechanisms by which diabetes modifies outcome in severe COVID-19 (41), to establish which risk factors are most useful to identify which people with diabetes are most vulnerable to COVID-19 (42), and to inform individualized care strategies (38).

Conclusion

Type 2 diabetes may be an independent prognostic factor for in-hospital survival in people with severe COVID-19 admitted to a critical care setting. The additional mortality risk associated with type 2 diabetes is attenuated by age. These results can help inform in-hospital decision-making on appropriate care escalation and treatment provision for people with type 2 diabetes and severe COVID-19.