Hostname: page-component-848d4c4894-x5gtn Total loading time: 0 Render date: 2024-06-12T04:20:44.439Z Has data issue: false hasContentIssue false

Changes in COVID-19-related outcomes, potential risk factors and disparities over time

Published online by Cambridge University Press:  10 August 2021

Youfei Yu
Affiliation:
Department of Biostatistics, University of Michigan School of Public Health, Ann Arbor, MI 48109, USA
Tian Gu
Affiliation:
Department of Biostatistics, University of Michigan School of Public Health, Ann Arbor, MI 48109, USA
Thomas S. Valley
Affiliation:
Division of Pulmonary and Critical Care Medicine, Department of Internal Medicine, University of Michigan Medical School, Ann Arbor, MI 48109, USA Institute for Healthcare Policy and Innovation, University of Michigan, Ann Arbor, MI 48109, USA
Bhramar Mukherjee*
Affiliation:
Department of Biostatistics, University of Michigan School of Public Health, Ann Arbor, MI 48109, USA Rogel Cancer Center, University of Michigan Medicine, Ann Arbor, MI 48109, USA Center for Statistical Genetics, University of Michigan School of Public Health, Ann Arbor, MI 48109, USA Department of Epidemiology, University of Michigan School of Public Health, Ann Arbor, MI 48109, USA
Lars G. Fritsche
Affiliation:
Department of Biostatistics, University of Michigan School of Public Health, Ann Arbor, MI 48109, USA Rogel Cancer Center, University of Michigan Medicine, Ann Arbor, MI 48109, USA Center for Statistical Genetics, University of Michigan School of Public Health, Ann Arbor, MI 48109, USA
*
Author for correspondence: Bhramar Mukherjee, E-mail: bhramar@umich.edu
Rights & Permissions [Opens in a new window]

Abstract

To investigate temporal trends in coronavirus disease 2019 (COVID-19)-related outcomes and to evaluate whether the impacts of potential risk factors and disparities changed over time, we conducted a retrospective cohort study with 249 075 patients tested or treated for COVID-19 at Michigan Medicine (MM), from 10 March 2020 to 3 May 2021. Among these patients, 26 289 were diagnosed with COVID-19. According to the calendar time in which they first tested positive, the COVID-19-positive cohort were stratified into three-time segments (T1: March–June, 2020; T2: July–December, 2020; T3: January–May, 2021). Potential risk factors that we examined included demographics, residential-level socioeconomic characteristics and preexisting comorbidities. The main outcomes included COVID-19-related hospitalisation and intensive care unit (ICU) admission. The hospitalisation rate for COVID-positive patients decreased from 36.2% in T1 to 14.2% in T3, and the ICU admission rate decreased from 16.9% to 2.9% from T1 to T3. These findings confirm that COVID-19-related hospitalisation and ICU admission rates were decreasing throughout the pandemic from March 2020 to May 2021. Black patients had significantly higher (compared to White patients) hospitalisation rates (19.6% vs. 11.0%) and ICU admission rates (6.3% vs. 2.8%) in the full COVID-19-positive cohort. A time-stratified analysis showed that racial disparities in hospitalisation rates persisted over time and the estimates of the odds ratios (ORs) stayed above unity in both unadjusted [full cohort: OR = 1.98, 95% confidence interval (CI) (1.79, 2.19); T1: OR = 1.70, 95% CI (1.36, 2.12); T2: OR = 1.40, 95% CI (1.17, 1.68); T3: OR = 1.55, 95% CI (1.29, 1.86)] and adjusted analysis, accounting for differences in demographics, socioeconomic status, and preexisting comorbid conditions (full cohort: OR = 1.45, 95% CI (1.25, 1.68); T1: OR = 1.26, 95% CI (0.90, 1.76); T2: OR = 1.29, 95% CI (1.01, 1.64); T3: OR = 1.29, 95% CI (1.00, 1.67)).

Type
Original Paper
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © The Author(s), 2021. Published by Cambridge University Press

Introduction

Since the World Health Organization declared coronavirus disease 2019 (COVID-19) a pandemic on 11 March 2020 [Reference Cucinotta and Vanelli1], there have been over 173.1 million confirmed cases worldwide, leading to 3.73 million deaths in the subsequent 15 months [2]. The United States has surpassed 33 million cases, with the state of Michigan reporting 994 935 confirmed cases and over 20 000 deaths as of 5 June 2021 [3]. Studies have found that racial and ethnic minority groups have been disproportionately affected by COVID-19 [Reference Webb Hooper, Nápoles and Pérez-Stable4Reference Xu15]. For example, higher incidence of confirmed cases [Reference Millett5, Reference Muñoz-Price6, Reference Tirupathi8, Reference Moore12], increased risk of hospitalisation [Reference Muñoz-Price6, Reference Azar7, Reference Gu9, Reference Gianfrancesco13], greater fatality burden [Reference Millett5, Reference Tirupathi8, Reference Ojinnaka10, Reference Alcendor11], and higher burden of years of potential life lost [Reference Bassett, Chen and Krieger14, Reference Xu15] were identified for Black individuals when compared to White individuals. Our previous study found certain risk factors, such as specific comorbidities, to affect outcomes in White and Black patients differently [Reference Gu9]. For example, preexisting type 2 diabetes was associated with higher risk of COVID-19-related hospitalisation in White patients but not in Black patients.

Emerging evidence shows an appreciable decline in mortality rates among patients with severe COVID-19 outcomes after the initial wave in March and April 2020 [Reference Dennis16Reference Griffin18]. A study conducted in a single health system in New York City found that the hospital mortality rate dropped from 25.6% at the start of the pandemic in March to 7.6% by mid-August 2020, adjusting for demographic and clinical factors [Reference Horwitz17]. The same trend of improved survival over time was seen among COVID-19 patients requiring hospitalisation [Reference Griffin18] or critical care management [Reference Dennis16] in England. Racial and ethnic disparities still remain, but recent data showed reduced differences in age-adjusted case fatality rate across ethnic groups [Reference Fabic, Choi and Bishai19]. Introduction of vaccines in the beginning of the year 2021 has also impacted COVID outcomes [Reference Tenforde20, Reference Cook and Roberts21]. To understand the dynamics of potential risk factors on COVID-19 outcomes over time, time-stratified analyses are required. Understanding how risk factors changed over the course of the pandemic can enhance the coordination of healthcare resources and help protect the most vulnerable. It can also help us understand whether targeted outreach/prevention efforts in the early months of the pandemic to specific vulnerable subgroups had an effect on improving their outcomes or if disparities persisted over time.

As a follow-up to the smaller cross-sectional study of Gu et al. [Reference Gu9], this time-stratified retrospective cohort study aims to examine the sociodemographic and clinical characteristics that are associated with various COVID-19 outcomes by race/ethnicity in a much larger COVID-19 cohort over a 15-month time period, as well as to characterise the risk factor trajectories over time, using electronic health records (EHRs) from Michigan Medicine (MM), a large academic health care system in the State of Michigan.

Methods

Study cohort and COVID-19 testing

The study sample consisted of 249 075 patients tested or treated for COVID-19 at MM, the University of Michigan Health System, from 10 March 2020 to 3 May 2021, presenting 26 289 testing positive. The descriptive statistics of the tested cohort are summarised in Table 1. Patients in the tested cohort received one of the five types of diagnostic tests: an in-house polymerase chain reaction (PCR) test (127 615 patients (52.1%)), a point of care PCR test (5781 patients (2.4%)), a commercial PCR test (Viracor; 418 patients (0.2%)), COVID-19 nasopharynx or oropharynx PCR tests deployed by the Michigan Department of Health and Human Services (52 patients (0.02%)), and a small fraction of ribonucleic acid (RNA) tests (three patients (0.001%)); 111 098 tested patients (45.4%) were transferred, tested elsewhere, or had no information on the type of testing they received.

Table 1. Descriptive characteristicsa of the COVID-19 tested or diagnosed cohort (n = 249 075)

Abbreviations: BMI, body mass index (calculated as weight in kilograms divided by height in metres squared); COVID-19, coronavirus disease 2019; ICU, intensive care unit; IQR, interquartile range; NDI, 2010 Neighbourhood Socioeconomic Disadvantage Index; MM, Michigan Medicine.

a Unless specified otherwise, numbers shown are numbers (%) of individuals in each category. Percentages are reported as fraction of column totals excluding missing entries.

b At least one visit with a primary care provider at MM within the last 2 years.

c Includes White Hispanic or unknown; Black Hispanic or unknown; Asian Hispanic, non-Hispanic, or unknown; Native American Hispanic, non-Hispanic, or unknown; Pacific Islander Hispanic, non-Hispanic, or unknown; and other Hispanic, non-Hispanic, or unknown.

d Includes missing race and/or ethnicity.

COVID-19 outcomes

In this study, we primarily focused on three COVID-19-related outcomes among individuals who tested positive: hospitalisation, intensive unit care (ICU) admission, and all-cause mortality. All patients who died post COVID-19 were included in the mortality counts. We performed sensitivity analysis by only considering outcomes that happened within 6 months of a COVID-19 diagnosis. We also investigated the factors associated with being tested (a question related to both access to testing and infection rates) and testing positive for COVID-19 (using a randomly selected untested comparison group described in the Supplementary Material). The corresponding results for these two secondary outcomes are largely reported in the Supplementary Material. The definitions of this ensemble of COVID-19 outcomes are listed in Supplementary Table S1, and the sample sizes for each outcome category are presented in Supplementary Figure S1.

Definition of demographics, socioeconomic status, comorbidities, and other covariates

We followed the same process described in Gu et al. [Reference Gu9] to extract the relevant variables from the EHRs. The sociodemographic variables considered included age, self-reported sex, race/ethnicity, smoking status, alcohol consumption, body mass index (BMI), Neighbourhood Socioeconomic Disadvantage Index (NDI) [22], and population density (in persons per square mile) [22]. NDI is a scale ranging from 0 to 1, with larger values indicating more disadvantaged communities. To characterise the underlying medical conditions of the patients, we considered seven comorbid conditions: respiratory conditions, circulatory conditions, any cancer, type 2 diabetes, kidney diseases, liver diseases and autoimmune diseases, all of which were coded binary. The general health status of patients was represented by a comorbidity score ranging from 0 to 7 that sums up the aforementioned seven conditions. In addition, all patients receiving medical care at least once at any MM primary care location since 1 January 2018 were considered as having MM as their primary care provider.

Statistical analysis

For each COVID-19-related outcome Y COVID, we fit a Firth's bias-corrected logistic regression to avoid the potential separation issues in the traditional maximum likelihood estimates. The analysis model is given by

(1)$$\matrix{ {{\rm logit}\;P( { {Y_{{\rm COVID}} = 1} \vert X, \;{\rm \;}Covariate} ) = \beta _0 + \beta _XX + \beta _{Cov}Covariate, \;\;\;} \cr } $$

where X and Covariate denote the predictors of interest (including Race) and the vector of adjustment covariates, respectively. To assess the sensitivity of our analysis to the choice of confounders/adjustment variables, we explored four nested sets of variables for Covariate, respectively labelled adjustment 0–3 (defined in Supplementary Table S1 and results reported in Supplementary Table S2). Our final model (adjustment 3) adjusted for age, sex, race/ethnicity, NDI and the composite comorbidity score. We excluded the cumulative comorbidity score when the factor of interest X was one of the individual comorbid conditions. The models for being tested and testing positive further adjusted for population density as that is likely to be related to disease transmission.

We further conducted a time-stratified analysis to assess the changes in racial/ethnic disparities and the effects of other risk factors for severe COVID-19-related outcomes over time since the start of the pandemic. According to the time period in which they first tested positive, the COVID-19-positive cohort were stratified into three groups: Time Period 1 (T1), 10 March 2020–30 June 2020; Time Period 2 (T2), 1 July 2020–31 December 2020; Time Period 3 (T3), 1 January 2021–3 May 2021. Similarly, the COVID-19-negative cohort were stratified based on the time period in which they first received the tests. The logistic regression model (1) was fitted using data for each time period.

The Race variable was self-reported and consisted of four categories: non-Hispanic White (White), non-Hispanic Black (Black), other known race/ethnicity, and unknown race/ethnicity. While the disproportionate burden of disease for Hispanics/Latinos is an important public health issue [Reference Alcendor11, Reference Weng23], our cohort did not have a sufficiently large sample size to afford the analysis. Therefore, we focus on comparing White and Black patients to evaluate racial/ethnic differences.

To evaluate the possibility of effect modification by race/ethnicity on the effect of potential risk factors, we fit a set of models with interactions by race:

$$\eqalign{& {\rm logit}\;P( { {Y_{{\rm COVID}} = 1} \vert X, \;\;Race, \;{\rm \;}Covariate} ) \cr & = \beta _0 + \beta _XX + \beta _{Race}Race + \beta _{int}X \times Race + \beta _{Cov}Covariate.} $$

For each model, we report the Firth's bias-corrected estimates of odds ratios along with their associated 95% Wald confidence intervals (CI) and P values. For the interaction models, we also obtained the P values of the differences in subgroup effects for White and Black patients by testing the null hypothesis H 0: β int = 0. All analyses were performed using complete cases in R statistical software version 4.0.3 (R Project for Statistical Computing).

Results

Descriptive statistics and unadjusted analysis

Table 1 reports the summary statistics for the tested cohort and each of the outcome subgroups. The tested cohort included 249 075 patients, with mean (s.d.) age of 44.2 (23.1), mean (s.d.) BMI of 29.0 (7.4), and 108 656 (43.6%) males. During the study period, more White patients (n = 182 475 (73.3%)) were tested than the other races/ethnicities, with 26 829 (73.7%), 112 912 (73.4%), and 39 799 (72.8%) White patients tested in T1, T2, and T3, respectively (Supplementary Table S3). The untested comparison group, which was used only for the secondary outcomes, consisted of 26 363 individuals, with mean (s.d.) age of 43.1 (24.6), mean (s.d.) BMI of 28.4 (7.4), and 12 285 (46.6%) males.

Among the 26 289 patients who tested positive, 3071 (11.7%) were hospitalised, 844 (3.2%) were admitted to the ICU, and 485 (1.8%) died. Older patients, male sex, ever smokers, and increased comorbidity burden tended to be associated with worsened disease severity among COVID-19-positive patients (Table 1). From T1 to T3, the mean and median ages decreased over time for patients who were hospitalised or admitted to the ICU. For example, the mean ages for hospitalised patients were 58.2, 54.2, and 49.0 years in T1, T2, and T3, respectively (Supplementary Table S3). This could potentially be attributed to the increased availability of vaccines for the elderly, the younger population in working-age groups returning to more in-person activities, and a harvesting effect.

Stratifying the descriptive statistics by White and Black patients, we noted differences in covariate distributions between the two groups (Supplementary Table S4). For example, in the tested cohort, Black patients tended to be younger (mean (s.d.) age, 41.4 (21.7) vs. 45.7 (23.2)), live in higher populated areas (mean (s.d.) population density in persons per square mile, 3439.4 (2270.7) vs. 1985.5 (2201.6)), and be of lower socioeconomic status (mean (s.d.) NDI, 0.21 (0.13) vs. 0.08 (0.06)) than White patients. Similar patterns existed in the subgroups (e.g. hospitalised patients) of the tested cohort. We assumed that the data were missing completely at random and used observed complete cases for the analysis. The fractions of missingness for each covariate are presented in Supplementary Table S5.

Figure 1 presents the changes in the overall and race-stratified COVID-19 outcomes over time. The overall hospitalisation rate among patients testing positive decreased from 36.2% in T1 to 14.2% in T3, and the overall ICU admission rate decreased from 16.9% to 2.9%. In the sensitivity analysis that only considers outcomes that happened within 6 months of a COVID-19 diagnosis, similar patterns were noted for the hospitalisation and ICU admission rates (Supplementary Figure S2). In the full cohort, Black patients had a significantly higher (unadjusted) hospitalisation rate (603/3078 (19.6%) vs. 2012/18 358 (11.0%)), ICU admission rate (193/3078 (6.3%) vs. 512/18 358 (2.8%)), and mortality rate (78/3101 (2.5%) vs. 317/18 499 (1.7%)) compared to White patients. The unadjusted odds ratios (ORs) for hospitalisation, ICU admission, and death were 1.98 (95% CI (1.79, 2.19)), 2.33 (95% CI (1.96, 2.77)), and 1.48 (95% CI (1.14, 1.91)), respectively. When we examined each time period separately, the point estimates of the hospitalisation rate remained higher for Black patients compared to White patients (Fig. 1) across the study period, but the (unadjusted) ORs were numerically lower in T2 (OR = 1.40, 95% CI (1.17, 1.68)) and T3 (OR = 1.55, 95% CI (1.29, 1.86)) than in T1 (OR = 1.70, 95% CI (1.36, 2.12)). The unadjusted OR for the ICU admission decreased from 1.54 (95% CI (1.16, 2.05)) in T1 to 0.99 (95% CI (0.64, 1.47)) in T2, followed by a rebound in T3 (OR = 1.80, 95% CI (1.25, 2.54)). The mortality rates for Black patients were lower than those of White patients in T1 (40/577 (6.9%) vs. 68/900 (7.6%)) and T2 (18/1136 (1.6%) vs. 161/9135 (1.8%)), and higher in T3 (11/971 (1.1%) vs. 41/5540 (0.7%)), though the differences were not statistically significant due to the small number of patients. At the end of the analysis, 110 patients were still hospitalised (n = 62) or still in ICU (n = 48).

Fig. 1. COVID-19 Outcomes Stratified by Race/Ethnicity in Each Time Period. Abbreviations: COVID-19, coronavirus disease 2019; ICU, intensive care unit; OR, odds ratio; T1, 10 March 2020, to 30 June 2020; T2, 1 July 2020, to 31 December 2020; T3, 1 January 2021, to 3 May 2021. + Logistic regression with Firth's correction. a Multivariable logistic regression with adjustment 1 (age + sex + race/ethnicity). b Multivariable logistic regression with adjustment 2 (adjustment 1 + Neighbourhood Socioeconomics Disadvantage Index). c Multivariable logistic regression with adjustment 3 (adjustment 2 + comorbidity score).

Covariate-adjusted analysis using multivariable logistic regression

Factors associated with COVID-19 testing

Overall, younger patients, female sex, higher BMI, ever smoker, alcohol consumption, Black race/ethnicity, lower NDI, areas with higher population densities, and accumulation of comorbidities were associated with an increased chance of getting tested (Supplementary Table S6). These associations were fairly consistent over time. For example, Black patients had a higher chance of getting tested across the study period (T1: OR = 1.35, 95% CI (1.26, 1.45); T2: OR = 1.06, 95% CI (0.99, 1.13); T3: OR = 1.27, 95% CI (1.18, 1.36)) than White patients.

Factors associated with testing positive

In the full cohort, Black patients had significantly higher risk of testing positive than White patients (OR = 1.48, 95% CI (1.37, 1.61)) (Supplementary Figure S2 and Table S6). Furthermore, age, male sex, ever smoker, and lower NDI were inversely associated with the odds of testing positive. All the comorbidity conditions considered were associated with an increased risk of testing positive. The time-stratified analysis revealed a sharp decrease in the odds of testing positive for Black patients relative to White patients from T1 to T3 (T1: OR = 4.17, 95% CI (3.52, 4.92); T2: OR = 1.19, 95% CI (1.07, 1.31); T3: OR = 1.52, 95% CI (1.36, 1.70)). A complete set of results for testing positive for COVID-19 are reported in the Supplementary Materials.

Factors associated with COVID-19-related outcomes among COVID-19-positive patients

Table 2 shows that among the patients who tested positive for COVID-19 through 3 May 2021, Black patients had significantly higher (covariate adjusted) odds of being hospitalised (OR = 1.45, 95% CI (1.25, 1.68)) and admitted to the ICU (OR = 1.37, 95% CI (1.05, 1.79)) compared to White patients. The racial/ethnic differences in mortality (OR = 1.31, 95% CI (0.93, 1.85)) were not statistically significant after covariate adjustment (Supplementary Table S6). The racial/ethnic disparities in hospitalisation rates persisted and the estimates of the ORs stayed above unity over time (T1: OR = 1.26, 95% CI (0.90, 1.76); T2: OR = 1.29, 95% CI (1.01, 1.64); T3: OR = 1.29, 95% CI (1.00, 1.67)). The racial/ethnic disparities in the ICU admission were not significant in each time period (T1: OR = 1.00, 95% CI (0.64, 1.58); T2: OR = 0.84, 95% CI (0.51, 1.39); T3: OR = 1.34, 95% CI (0.81, 2.22)), possibly due to the reduced sample size after stratification.

Table 2. Odds ratios (95% confidence intervals)a of COVID-19-related outcomes from logistic regression

Abbreviations: BMI, body mass index; ICU, intensive unit care; NDI, Neighbourhood Socioeconomic Disadvantage Index; NA, not applicable.

Odds ratios that are significantly greater than 1 at the level of 0.05 are shaded orange. Time period 1, 10 March 2020, to 30 June 2020; Time period 2, 1 July 2020, to 31 December 2020; Time period 3, 1 January 2021, to 3 May 2021.

a Results were obtained from Firth's multivariable logistic regression model logit P(Y COVID = 1|X, Covariate) = β 0 + β XX + β CovCovariate, where Y COVID is the COVID-19 outcomes (i.e. positive test results, hospitalisation, admission to ICU, or mortality). Covariate = age + sex + race + NDI + comorbidity score.

b Includes White Hispanic or unknown; Black Hispanic or unknown; Asian Hispanic, non-Hispanic, or unknown; Native American Hispanic, non-Hispanic, or unknown; Pacific Islander Hispanic, non-Hispanic, or unknown; and other Hispanic, non-Hispanic, or unknown.

c Not adjusted for composite comorbidity score.

Additionally, age, male sex, smoking, living in densely populated areas or disadvantaged neighbourhoods, and heavier overall comorbidity burden were positively associated with hospitalisation and ICU admission in the full COVID-positive cohort (Table 2). Most of the associations in each time period stayed in the same direction as in the full cohort analysis (Table 2). Specifically, older patients (every 10-year increase in age) were consistently associated with increased risk of being hospitalised (T1: OR = 1.20, 95% CI (1.11, 1.29); T2: OR = 1.12, 95% CI (1.08, 1.16); T3: OR = 1.13, 95% CI (1.08, 1.18)). Among the underlying conditions, kidney diseases in general had the largest ORs of hospitalisation (T1: OR = 3.53, 95% CI (2.55, 4.88); T2: OR = 4.46, 95% CI (3.80, 5.25); T3: OR = 4.55, 95% CI (3.74, 5.53)) and ICU admission (T1: OR = 2.83, 95% CI (1.91, 4.19); T2: OR = 4.67, 95% CI (3.48, 6.27); T3: OR = 5.04, 95% CI (3.47, 7.31)) across the study periods.

Interaction analysis by race/ethnicity within the COVID-19-positive cohort

Figure 2 displays the ORs of potential risk factors stratified by race/ethnicity for the hospitalisation outcome (Fig. 2A) and ICU admission (Fig. 2B). The main effects of the seven comorbid conditions considered were significant in both Black and White patients for hospitalisation, while none of the comorbidities exhibited significant interaction effects. For example, the OR of kidney diseases for Black and White patients were 3.88 (95% CI (3.07, 4.91), P < 0.001) and 3.86 (95% CI (3.38, 4.40), P < 0.001), respectively, and the interaction P value was 0.96. Higher NDI (i.e. lower socioeconomic status) had significant positive associations with hospitalisation in White patients (OR = 9.29, 95% CI (4.15, 20.8), P < 0.001), but not in Black patients (OR = 1.64, 95% CI (0.69, 3.87), P = 0.263), with the interaction being significant (Pint = 0.004). For the ICU admission (Fig. 2B), significant interaction effects were observed in patients with obesity (BMI ≥30) (White: OR = 0.94, 95% CI (0.69, 1.26), P = 0.662; Black: OR = 2.20, 95% CI (1.08, 4.46), P = 0.029; Pint = 0.030), alcohol drinkers (White: OR = 0.74, 95% CI (0.58, 0.96), P = 0.020; Black: OR = 1.34, 95% CI (0.86, 2.07), P = 0.193; Pint = 0.023), and patients with underlying liver diseases (White: OR = 1.83, 95% CI (1.35, 2.48), P < 0.001; Black: OR = 0.82, 95% CI (0.40, 1.70), P = .592; Pint = 0.046) and autoimmune diseases (White: OR = 1.09, 95% CI (0.82, 1.46), P = 0.550; Black: OR = 1.89, 95% CI (1.25, 2.86), P = 0.002; Pint = 0.031).

Fig. 2. Hospitalisation (A) and ICU Admission (B) for Black and White Patients in the Full Cohort. Abbreviations: BMI, body mass index; NDI, Neighbourhood Socioeconomic Disadvantage Index. The results were from the model logit P(Y COVID = 1|X, Covariate) = β 0 + β XX + β RaceRace + β intX × Race + β CovCovariate, where Y COVID denotes hospitalisation (A) or ICU admission (B), and Covariate = age + sex + NDI ( + comorbidity score in the demographic and socioeconomic status models). Results that are statistically significant at the level of 0.05 are in bold.

Discussion

In this study, we evaluated the potential risk factors and racial/ethnic disparities for COVID-19 prognosis using a larger cohort from the same healthcare system that covers a longer time frame than Gu et al. [Reference Gu9]. We observed that Black patients and patients with higher comorbidity burden were more susceptible to COVID-19 and more likely to be hospitalised, a result consistent with previous findings in the literature. On the other hand, the association between age and positive test results reversed compared to what was reported in Gu et al. [Reference Gu9], such that younger people, especially those aged 18–35 years, were more likely to be tested positive. The finding aligns with recent data from the Centers for Disease Control and Prevention, which shows an increase in the COVID-19 incidence in persons aged <30 years and a decrease in the median age for confirmed cases in May−August 2020 [Reference Boehmer24].

A fraction of patients who were hospitalised (n = 249 of 3071 (8.1%)) in our sample were transferred from external hospitals, mainly from the Detroit Metro area. These transferred patients tended to have more severe symptoms and different distributions in sociodemographic characteristics than non-transferred patients. It is important to note the temporal changes in the context of care at MM (Supplementary Table S7). There was an enrichment in transferred patients in T1 (n = 123 (19.3%)) compared to T2 (n = 76 (5.6%)) and T3 (n = 50 (4.7%)) for those hospitalised. However, a sensitivity analysis that excluded patients who did not receive their primary care at MM was consistent with the overall analysis, namely improved outcomes and reduced but persisting health disparities over time (Supplementary Table S8).

Our definition of COVID-19 related hospitalisation, ICU admission, and mortality may be liberal. We performed a sensitivity analysis restricting the time period to six months since the first positive diagnosis. For example, of the 485 deaths, 458 (94.4%) happened within six months. The results of the sensitivity analysis (Supplementary Figure S3) showed consistent patterns with the analysis based on unfiltered outcomes (Fig. 1) with respect to the changes in the rates and ORs of the three outcomes. Specifically, the overall hospitalisation rate, ICU admission rate, and mortality rate decreased from 34.2%, 16.5%, and 7.1% to 9.6%, 1.9%, and 0.8%, respectively.

Our findings suggest that the COVID-19-related hospitalisation, ICU admission, and mortality rates declined from T1 to T2, which may be partly explained by the rapid increase in testing and adoption of preventive measures. The rates of these severe outcomes remained relatively steady from T2 to T3. However, racial/ethnic disparities persisted over time. This study also confirms a number of potential risk factors for COVID-19-related outcomes: older age, male sex, Black race/ethnicity, preexisting conditions, and lower socioeconomic status. Among the examined comorbidities, underlying kidney diseases in general had the largest odds ratios for the COVID-19-related outcomes. The interaction analysis identified several risk factors that affected the hospitalisation and ICU admission differently in Black and White patients. Specifically, underlying respiratory diseases and liver diseases increased the risk of hospitalisation in the subgroup of White patients, while the effects on Black patients were not significant. On the other hand, obesity and autoimmune diseases were positively associated with ICU admission in Black patients but not in White patients.

There are several limitations of our study. First, the sample cannot be considered a population-based random sample, as MM prioritised testing to those presenting COVID-19 symptoms or at greatest risk of exposure, particularly during the early stages of the pandemic when the availability of tests was limited. Second, since the hospitalisation records were available only for those treated at MM, not all hospitalised patients were captured by our data. Patients hospitalised elsewhere may be classified as non-hospitalised, which may lead to biased estimates. Finally, 486 out of the 52 325 people who were first tested in T3 received at least one vaccination dose by May 3, but we did not consider vaccination as a covariate, which may explain part of the reduction in hospitalisation and mortality. Despite the limitations, our study assessed the temporal trend in COVID-19-related outcomes and evaluated whether the impacts of potential risk factors and disparities changed over time. Elderly adults, Black patients, and patients with underlying health conditions were disproportionately affected by COVID-19 over the course of the pandemic, which calls for vaccine prioritisation and targeted outreach in these population groups.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0950268821001898.

Author contributions

Conceptualisation, B.M., L.G.F., T.G., Y.Y., and T.S.V; methodology, B.M., L.G.F., T.G., and Y.Y.; software, L.G.F. and T.G.; validation, L.G.F. and T.G.; formal analysis, L.F.G. and T.G; investigation, all authors; resources, B.M. and T.S.V.; data curation, L.G.F; writing—original draft preparation, Y.Y., L.G.F., B.M., T.S.V.; visualisation, L.G.F. and T.G.; supervision, B.M.; project administration, B.M.; funding acquisition, B.M. All authors have read and agreed to the published version of the manuscript.

Financial support

This study was funded by the University of Michigan Precision Health Initiative, University of Michigan Rogel Cancer Center, and Michigan Institute for Data Science. Dr Mukherjee's research was funded by grant No. NSF DMS 1712933 from the National Science Foundation, and Dr Fritsche's research was supported by grant No. CA 046592 from the National Cancer Institute, National Institutes of Health.

Institutional Review Board Statement

The University of Michigan Medical School institutional review board reviewed the study and determined that it is exempt. Ethical review and approval were waived for this study per institutional policy.

Informed Consent Statement

The institutional review board waived the need to obtain consent for the collection, analysis, and publication of the anonymised COVID-19 data, per institutional policy.

Conflict of interest

The authors declare no conflict of interest.

Data

The individual level data cannot be shared publicly due to patient confidentiality. The data underlying the results presented in the study are available from University of Michigan Data Office for Clinical & Translational Research for researchers who meet the criteria for access to confidential data.

References

Cucinotta, D and Vanelli, M (2020) WHO Declares COVID-19 a pandemic. Acta Bio Medica Atenei Parmensis 91, 157160.Google ScholarPubMed
Johns Hopkins Coronavirus Resource Center. COVID-19 Map. Available at: (https://coronavirus.jhu.edu/map.html). Accessed 6 June 2021.Google Scholar
Coronavirus - Michigan Data. Available at: (https://www.michigan.gov/coronavirus/0,9753,7-406-98163_98173---,00.html). Accessed 6 June 2021.Google Scholar
Webb Hooper, M, Nápoles, AM and Pérez-Stable, EJ (2020) COVID-19 and racial/ethnic disparities. JAMA 323, 2466.CrossRefGoogle ScholarPubMed
Millett, GA et al. (2020) Assessing differential impacts of COVID-19 on Black communities. Annals of Epidemiology 47, 3744.CrossRefGoogle ScholarPubMed
Muñoz-Price, LS et al. (2020) Racial disparities in incidence and outcomes among patients with COVID-19. JAMA Network Open 3, e2021892.CrossRefGoogle ScholarPubMed
Azar, KMJ et al. (2020) Disparities in outcomes among COVID-19 patients in a large health care system in California. Health Affairs Health Affairs 39, 12531262.CrossRefGoogle Scholar
Tirupathi, R et al. (2020) COVID-19 disparity among racial and ethnic minorities in the US: a cross sectional analysis. Travel Medicine and Infectious Disease 38, 101904.CrossRefGoogle ScholarPubMed
Gu, T et al. (2020) Characteristics associated with racial/ethnic disparities in COVID-19 outcomes in an academic health care system. JAMA Network Open 3, e2025197.CrossRefGoogle Scholar
Ojinnaka, CO et al. (2020) Factors associated with COVID-related mortality: the case of Texas. Journal of Racial and Ethnic Health Disparities Nov 9, 19. doi: 10.1007/s40615-020-00913-5.Google Scholar
Alcendor, DJ (2020) Racial disparities-associated COVID-19 mortality among minority populations in the US. Journal of Clinical Medicine Multidisciplinary Digital Publishing Institute 9, 2442.Google ScholarPubMed
Moore, JT (2020) Disparities in incidence of COVID-19 among underrepresented racial/ethnic groups in counties identified as hotspots during June 5–18, 2020 — 22 states, February–June 2020. MMWR. Morbidity and Mortality Weekly Report 69(33), 11221126. doi: 10.15585/mmwr.mm6933e1.CrossRefGoogle ScholarPubMed
Gianfrancesco, MA et al. (2020) Race/ethnicity association with COVID-19 outcomes in rheumatic disease: data from the COVID-19 global rheumatology alliance physician registry. Arthritis & Rheumatology (Hoboken, N.J.) 73(3), 374380. doi: 10.1002/art.41567.CrossRefGoogle Scholar
Bassett, MT, Chen, JT and Krieger, N (2020) Variation in racial/ethnic disparities in COVID-19 mortality by age in the United States: a cross-sectional study. PLOS Medicine Public Library of Science 17, e1003402.Google ScholarPubMed
Xu, JJ et al. (2021) Racial and ethnic disparities in years of potential life lost attributable to COVID-19 in the United States: an analysis of 45 states and the District of Columbia. International Journal of Environmental Research and Public Health 18(6), 2921. doi: 10.3390/ijerph18062921.CrossRefGoogle ScholarPubMed
Dennis, JM et al. (2020) Improving survival of critical care patients with coronavirus disease 2019 in England: a national cohort study, March to June 2020. Critical Care Medicine 49(2), 209214. doi: 10.1097/CCM.0000000000004747.CrossRefGoogle Scholar
Horwitz, LI et al. (2020) Trends in COVID-19 risk-adjusted mortality rates. Journal of Hospital Medicine 16(2), 9092. doi: 10.12788/jhm.3552.CrossRefGoogle Scholar
Griffin, S (2020) Covid-19: continued outbreaks in care homes risk extending pandemic, say experts. BMJ British Medical Journal Publishing Group 369(m2530), 1. doi: 10.1136/bmj.m2530.Google ScholarPubMed
Fabic, MS, Choi, Y and Bishai, D (2020) Deaths among COVID Cases in the United States: Racial and Ethnic Disparities Persist. MedRxiv. doi: 10.1101/2020.11.15.20232066.Google Scholar
Tenforde, MW (2021) Effectiveness of Pfizer-BioNTech and Moderna vaccines against COVID-19 among hospitalized adults aged ⩾65 years — United States, January–March 2021. MMWR. Morbidity and Mortality Weekly Report 70(18), 674679. doi: 10.15585/mmwr.mm7018e1.CrossRefGoogle ScholarPubMed
Cook, TM and Roberts, JV (2021) Impact of vaccination by priority group on UK deaths, hospital admissions and intensive care admissions from COVID-19. Anaesthesia 76, 608616. doi: 10.1111/anae.15442.CrossRefGoogle Scholar
Weng, C-H et al. (2021) Persistently high SARS-CoV-2 positivity rate and incidence for Hispanic/Latinos during state reopening in an urban setting: a retrospective cohort study. Epidemiology and Infection 149, e25.CrossRefGoogle Scholar
Boehmer, TK (2020) Changing age distribution of the COVID-19 pandemic — United States, May–August 2020. MMWR. Morbidity and Mortality Weekly Report 69(39), 14041409. doi: 10.15585/mmwr.mm6939e1.CrossRefGoogle ScholarPubMed
Figure 0

Table 1. Descriptive characteristicsa of the COVID-19 tested or diagnosed cohort (n = 249 075)

Figure 1

Fig. 1. COVID-19 Outcomes Stratified by Race/Ethnicity in Each Time Period. Abbreviations: COVID-19, coronavirus disease 2019; ICU, intensive care unit; OR, odds ratio; T1, 10 March 2020, to 30 June 2020; T2, 1 July 2020, to 31 December 2020; T3, 1 January 2021, to 3 May 2021. + Logistic regression with Firth's correction. a Multivariable logistic regression with adjustment 1 (age + sex + race/ethnicity). b Multivariable logistic regression with adjustment 2 (adjustment 1 + Neighbourhood Socioeconomics Disadvantage Index). c Multivariable logistic regression with adjustment 3 (adjustment 2 + comorbidity score).

Figure 2

Table 2. Odds ratios (95% confidence intervals)a of COVID-19-related outcomes from logistic regression

Figure 3

Fig. 2. Hospitalisation (A) and ICU Admission (B) for Black and White Patients in the Full Cohort. Abbreviations: BMI, body mass index; NDI, Neighbourhood Socioeconomic Disadvantage Index. The results were from the model logit P(YCOVID = 1|X, Covariate) = β0 + βXX + βRaceRace + βintX × Race + βCovCovariate, where YCOVID denotes hospitalisation (A) or ICU admission (B), and Covariate = age + sex + NDI ( + comorbidity score in the demographic and socioeconomic status models). Results that are statistically significant at the level of 0.05 are in bold.

Supplementary material: File

Yu et al. supplementary material

Yu et al. supplementary material

Download Yu et al. supplementary material(File)
File 314 KB