Skip to main content
Advertisement
  • Loading metrics

Patterns of COVID-19 testing and mortality by race and ethnicity among United States veterans: A nationwide cohort study

  • Christopher T. Rentsch ,

    Roles Conceptualization, Data curation, Formal analysis, Methodology, Project administration, Visualization, Writing – original draft, Writing – review & editing

    christopher.rentsch@lshtm.ac.uk

    Affiliations VA Connecticut Healthcare System, US Department of Veterans Affairs, West Haven, Connecticut, United States of America, Faculty of Epidemiology and Population Health, London School of Hygiene & Tropical Medicine, London, United Kingdom

  • Farah Kidwai-Khan ,

    Contributed equally to this work with: Farah Kidwai-Khan, Janet P. Tate, Lesley S. Park, Joseph T. King Jr, Melissa Skanderson

    Roles Data curation, Writing – review & editing

    Affiliations VA Connecticut Healthcare System, US Department of Veterans Affairs, West Haven, Connecticut, United States of America, Department of Internal Medicine, Yale School of Medicine, New Haven, Connecticut, United States of America

  • Janet P. Tate ,

    Contributed equally to this work with: Farah Kidwai-Khan, Janet P. Tate, Lesley S. Park, Joseph T. King Jr, Melissa Skanderson

    Roles Data curation, Methodology, Visualization, Writing – review & editing

    Affiliations VA Connecticut Healthcare System, US Department of Veterans Affairs, West Haven, Connecticut, United States of America, Department of Internal Medicine, Yale School of Medicine, New Haven, Connecticut, United States of America

  • Lesley S. Park ,

    Contributed equally to this work with: Farah Kidwai-Khan, Janet P. Tate, Lesley S. Park, Joseph T. King Jr, Melissa Skanderson

    Roles Data curation, Methodology, Visualization, Writing – review & editing

    Affiliation Stanford Center for Population Health Sciences, Stanford University School of Medicine, Stanford, California, United States of America

  • Joseph T. King Jr ,

    Contributed equally to this work with: Farah Kidwai-Khan, Janet P. Tate, Lesley S. Park, Joseph T. King Jr, Melissa Skanderson

    Roles Data curation, Writing – review & editing

    Affiliations VA Connecticut Healthcare System, US Department of Veterans Affairs, West Haven, Connecticut, United States of America, Department of Neurosurgery, Yale School of Medicine, New Haven, Connecticut, United States of America

  • Melissa Skanderson ,

    Contributed equally to this work with: Farah Kidwai-Khan, Janet P. Tate, Lesley S. Park, Joseph T. King Jr, Melissa Skanderson

    Roles Data curation, Writing – review & editing

    Affiliation VA Connecticut Healthcare System, US Department of Veterans Affairs, West Haven, Connecticut, United States of America

  • Ronald G. Hauser,

    Roles Data curation, Methodology, Writing – review & editing

    Affiliations VA Connecticut Healthcare System, US Department of Veterans Affairs, West Haven, Connecticut, United States of America, Department of Laboratory Medicine, Yale School of Medicine, New Haven, Connecticut, United States of America

  • Anna Schultze,

    Roles Data curation, Writing – review & editing

    Affiliation Faculty of Epidemiology and Population Health, London School of Hygiene & Tropical Medicine, London, United Kingdom

  • Christopher I. Jarvis,

    Roles Visualization, Writing – review & editing

    Affiliation Faculty of Epidemiology and Population Health, London School of Hygiene & Tropical Medicine, London, United Kingdom

  • Mark Holodniy,

    Roles Data curation, Writing – review & editing

    Affiliations VA Palo Alto Health Care System, US Department of Veterans Affairs, Palo Alto, California, United States of America, Division of Infectious Diseases and Geographic Medicine, Stanford University School of Medicine, Stanford, California, United States of America

  • Vincent Lo Re III,

    Roles Methodology, Writing – review & editing

    Affiliations Division of Infectious Diseases, Department of Medicine, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America, Center for Clinical Epidemiology and Biostatistics, Department of Biostatistics, Epidemiology, and Informatics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America

  • Kathleen M. Akgün,

    Roles Methodology, Writing – review & editing

    Affiliations VA Connecticut Healthcare System, US Department of Veterans Affairs, West Haven, Connecticut, United States of America, Department of Internal Medicine, Yale School of Medicine, New Haven, Connecticut, United States of America

  • Kristina Crothers,

    Roles Methodology, Writing – review & editing

    Affiliations VA Puget Sound Health Care System, US Department of Veterans Affairs, Seattle, Washington, United States of America, Department of Medicine, University of Washington School of Medicine, Seattle, Washington, United States of America

  • Tamar H. Taddei,

    Roles Writing – review & editing

    Affiliations VA Connecticut Healthcare System, US Department of Veterans Affairs, West Haven, Connecticut, United States of America, Department of Internal Medicine, Yale School of Medicine, New Haven, Connecticut, United States of America

  • Matthew S. Freiberg,

    Roles Methodology, Writing – review & editing

    Affiliations Geriatric Research Education and Clinical Center, Tennessee Valley Healthcare System, US Department of Veterans Affairs, Nashville, Tennessee, United States of America, Department of Medicine, Vanderbilt University Medical Center, Nashville, Tennessee, United States of America

  •  [ ... ],
  • Amy C. Justice

    Roles Conceptualization, Funding acquisition, Methodology, Project administration, Writing – original draft, Writing – review & editing

    Affiliations VA Connecticut Healthcare System, US Department of Veterans Affairs, West Haven, Connecticut, United States of America, Department of Internal Medicine, Yale School of Medicine, New Haven, Connecticut, United States of America, Center for Interdisciplinary Research on AIDS, Yale School of Public Health, New Haven, Connecticut, United States of America

  • [ view all ]
  • [ view less ]

Abstract

Background

There is growing concern that racial and ethnic minority communities around the world are experiencing a disproportionate burden of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and coronavirus disease 2019 (COVID-19). We investigated racial and ethnic disparities in patterns of COVID-19 testing (i.e., who received testing and who tested positive) and subsequent mortality in the largest integrated healthcare system in the United States.

Methods and findings

This retrospective cohort study included 5,834,543 individuals receiving care in the US Department of Veterans Affairs; most (91%) were men, 74% were non-Hispanic White (White), 19% were non-Hispanic Black (Black), and 7% were Hispanic. We evaluated associations between race/ethnicity and receipt of COVID-19 testing, a positive test result, and 30-day mortality, with multivariable adjustment for a wide range of demographic and clinical characteristics including comorbid conditions, health behaviors, medication history, site of care, and urban versus rural residence. Between February 8 and July 22, 2020, 254,595 individuals were tested for COVID-19, of whom 16,317 tested positive and 1,057 died. Black individuals were more likely to be tested (rate per 1,000 individuals: 60.0, 95% CI 59.6–60.5) than Hispanic (52.7, 95% CI 52.1–53.4) and White individuals (38.6, 95% CI 38.4–38.7). While individuals from minority backgrounds were more likely to test positive (Black versus White: odds ratio [OR] 1.93, 95% CI 1.85–2.01, p < 0.001; Hispanic versus White: OR 1.84, 95% CI 1.74–1.94, p < 0.001), 30-day mortality did not differ by race/ethnicity (Black versus White: OR 0.97, 95% CI 0.80–1.17, p = 0.74; Hispanic versus White: OR 0.99, 95% CI 0.73–1.34, p = 0.94). The disparity between Black and White individuals in testing positive for COVID-19 was stronger in the Midwest (OR 2.66, 95% CI 2.41–2.95, p < 0.001) than the West (OR 1.24, 95% CI 1.11–1.39, p < 0.001). The disparity in testing positive for COVID-19 between Hispanic and White individuals was consistent across region, calendar time, and outbreak pattern. Study limitations include underrepresentation of women and a lack of detailed information on social determinants of health.

Conclusions

In this nationwide study, we found that Black and Hispanic individuals are experiencing an excess burden of SARS-CoV-2 infection not entirely explained by underlying medical conditions or where they live or receive care. There is an urgent need to proactively tailor strategies to contain and prevent further outbreaks in racial and ethnic minority communities.

Author summary

Why was this study done?

  • There is growing concern that racial and ethnic minority communities around the world are experiencing a disproportionate burden of morbidity and mortality from symptomatic severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection or coronavirus disease 2019 (COVID-19).
  • Most studies investigating racial and ethnic disparities to date have focused on those who tested positive for SARS-CoV-2 or hospitalized patients.
  • No single study to our knowledge has yet investigated racial and ethnic disparities in testing patterns (i.e., who received testing and who tested positive) as well as COVID-19 outcomes in a nationwide cohort with adequate adjustment for potential confounders.

What did the researchers do and find?

  • We used electronic health records from the largest integrated healthcare system in the US to investigate racial and ethnic disparities in testing and subsequent COVID-19 mortality.
  • Non-Hispanic Black and Hispanic individuals were twice as likely as non-Hispanic White individuals to test positive for COVID-19, even after accounting for other demographics, geographic location, and underlying health conditions.
  • The racial disparity between Black and White individuals in testing positive for COVID-19 slightly decreased over the study period, and was highest in the Midwest compared to all other regions. The ethnic disparity between Hispanic and White individuals in testing positive for COVID-19 was consistent across time, geographic region, and outbreak pattern; the disparity was consistently observed across all strata.
  • Among those who tested positive for COVID-19, there was no observed difference in 30-day mortality by race/ethnicity group.

What do these findings mean?

  • Our findings highlight the urgent need for improved strategies to contain and prevent further outbreaks in racial and ethnic minority communities.

Introduction

The United States has the highest number of reported symptomatic severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections and related deaths in the world, accounting for one-fourth of global totals as of July 22, 2020 [1]. There is growing concern that racial and ethnic minority communities are experiencing a disproportionate burden of morbidity and mortality from symptomatic SARS-CoV-2 infection or coronavirus disease 2019 (COVID-19) [28]. One statewide study investigating racial disparities followed 3,481 COVID-19 cases in Louisiana and found that non-Hispanic Black individuals represented 77% of hospitalizations and 71% of deaths despite only making up 31% of the total source population [9]. Thus, the potential for racial and ethnic disparities in COVID-19 have been deemed an urgent public health research priority [10]. However, most studies investigating racial and ethnic disparities have focused on hospitalized patients or have not characterized who received testing or tested positive for COVID-19 [9,1115]. Given that COVID-19 testing was not performed at random, particularly in the early phases of the pandemic, evaluating underlying testing patterns and changes over time may provide important context for interpreting findings from models of COVID-19 outcomes. In addition, it is not yet known whether disparities in COVID-19 infection or severe outcomes are explained, at least in part, by differences in underlying health conditions, smoking and alcohol use, geographic location, or urban versus rural residence—essential information if we are to design effective interventions.

The electronic health record database of the Department of Veterans Affairs (VA) offers the single largest nationwide data resource available with the necessary information on system-wide testing and detailed medical histories to examine racial and ethnic disparities in the US. We evaluated associations between race/ethnicity and receipt of COVID-19 testing, a positive test result, and 30-day mortality, conditioning each analysis on the previous outcome and accounting for a wide range of demographic and clinical characteristics through July 22, 2020.

Methods

Data source

The VA is the largest integrated healthcare system in the US and comprises over 1,200 points of care (i.e., sites) nationwide including hospitals, medical centers, and community outpatient clinics. All care is recorded in an electronic health record with daily uploads into the VA Corporate Data Warehouse. Available data include demographics, outpatient and inpatient encounters, diagnoses, smoking and alcohol health behaviors, and pharmacy dispensing records.

This study was approved by the institutional review boards of VA Connecticut Healthcare System and Yale University. It has been granted a waiver of informed consent and is Health Insurance Portability and Accountability Act compliant. The analyses herein were not pre-specified in a formal protocol, rather were informed by hypotheses drawn from prior work [16]. This study is reported as per the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guidelines (S1 STROBE Checklist).

Sample, follow-up, and outcomes

All individuals in clinical care (defined as having at least 1 clinical encounter between January 1, 2018, and December 31, 2019, and alive as of January 1, 2020) were included in this analysis. We identified individuals tested for COVID-19 from the date of the VA’s first recorded test, on February 8, 2020, through July 22, 2020, by using text searching of laboratory results containing terms consistent with SARS-CoV-2 or COVID-19. Nearly all tests utilized nasopharyngeal swabs; 1% were from other sources. Testing was performed in VA, state public health, and commercial reference laboratories using FDA Emergency Use Authorization–approved SARS-CoV-2 assays. We did not include antibody tests in this analysis.

If an individual had more than 1 test and all were negative, we selected the date of the first negative test; otherwise we used the date of the first positive test. Baseline for individuals tested for COVID-19 was defined as the date of specimen collection unless testing occurred during hospitalization, in which case baseline was defined as the date of admission. If the admission began more than 14 days prior to testing, which may indicate hospital-acquired infection, we set baseline to 14 days prior to testing to better capture health status prior to SARS-CoV-2 infection. We examined 3 outcomes: (1) receipt of COVID-19 testing among all in care, (2) receipt of a positive test result among individuals tested for COVID-19, and (3) 30-day mortality among COVID-19 cases. Deaths were ascertained using inpatient records and VA death registry data to capture deaths outside of hospitalization. The choice of 30-day mortality as the outcome was guided by the distribution of mortality events by time since testing positive for COVID-19 (50th, 75th, 90th percentile time to death: 12, 20, 30 days) and to allow for sufficient follow-up within the study period. While there were some deaths beyond 30 days after testing positive for COVID-19, we were less certain that these deaths could be attributed to COVID-19. Given the low number of deaths after 30 days, 30-day mortality may be a reasonable proxy for case fatality rate. However, until longer follow-up has accrued, it remains to be seen whether those who develop symptomatic COVID-19 experience longer term excess mortality.

Variables

The primary exposure variable was self-reported race/ethnicity (non-Hispanic White [White], non-Hispanic Black [Black], and Hispanic). Analyses of other racial and ethnic backgrounds were underpowered at the time of this analysis, and therefore individuals who self-reported race/ethnicity other than White, Black, or Hispanic were excluded from the study population.

We selected demographic and clinical characteristics that have been evaluated in prior COVID-19 reports and could potentially mediate or explain racial/ethnic disparities in COVID-19 positivity and mortality. Demographics included age at baseline, sex, and rural/urban residence. Rural/urban residence was defined using geographic information system coding based upon established criteria [17]. Clinical characteristics were based on diagnostic codes for asthma, any cancer, chronic obstructive pulmonary disease (COPD), chronic kidney disease, diabetes mellitus, hypertension, liver disease, vascular disease, and alcohol use disorder (definitions provided in S1 Table). Presence of conditions was determined by 1 inpatient or 2 outpatient diagnoses in the 2 years prior to baseline, except for cancer, which was considered present if diagnosed ever prior to baseline. Diagnoses made in the 7 days prior to baseline were not included. We used a validated algorithm to capture smoking status [18] and alcohol consumption [19]. We collected pharmacy fills for angiotensin converting enzyme (ACE) inhibitors and angiotensin II receptor blockers (ARBs) and identified individuals with active prescriptions in the 30 days prior to baseline. Missing data for smoking and alcohol consumption affected only 5% of individuals included in multivariable models; thus, complete case analysis was performed.

We also created variables to assess potential variation in racial/ethnic disparities by calendar time, region, and outbreak pattern. We split the population into 3 groups based on date of COVID-19 test: February 8 to April 21, April 22 to June 21, and June 22 to July 22. States were grouped into 4 US Census regions (i.e., West, South, Midwest, and Northeast) [20]. Outbreak patterns were based on site-level percentage of positive tests per month among sites with at least 100 positive COVID-19 tests: early (≥10% in March or April), late (≥10% in June or July), resurgent (≥10% in March or April and June or July), steady (<10% in all months), and other (sites with <100 positive tests).

Statistical analysis

We calculated COVID-19 testing rate per 1,000 individuals in care and Clopper–Pearson 95% confidence intervals (CIs) by race/ethnicity category. Among those tested for COVID-19, we calculated percent testing positive and 95% CIs by race/ethnicity category. Logistic regression models were used to estimate associations between race/ethnicity and COVID-19 positivity and mortality, adjusting for sets of potential mediators of such disparities, moving from more distal to more proximate determinants of health. Age-adjusted models included race/ethnicity and age. Demographic-adjusted models additionally included sex and rural/urban residence, and were conditioned on site of care. Fully adjusted models additionally included all clinical covariates, substance use, and medication history. We report the estimates of each individual adjustment as well as those arising from a fully adjusted model. We repeated this modeling strategy to estimate odds ratios (ORs) and 95% CIs between race/ethnicity and 30-day mortality among those who tested positive for COVID-19 on or prior to June 21, 2020, to allow all individuals 30 days of follow-up. We evaluated variation in racial/ethnic disparities in testing positive for COVID-19 by stratifying the fully adjusted model by calendar time, geographic region, and site-level outbreak pattern. In sensitivity analyses, we restricted ascertainment of 30-day mortality to only include inpatient deaths to test the robustness of the associations found in the primary models. Analyses were performed using SAS version 9.4 (SAS Institute, Cary, NC, US). R version 3.6.3 was used to map COVID-19 cases nationwide.

Results

There were 5,834,543 individuals in care prior to the COVID-19 pandemic. Most (91%) were men, 74% were White, 19% were Black, and 7% were Hispanic (Table 1). Age distributions were similar by race/ethnicity, and age ranged from 20 to 105 years, with 24% less than 50 years, 35% 50–69 years, 28% 70–79 years, and 13% ≥70 years of age. Of these, 254,595 (43.6 per 1,000 individuals in care) were tested for COVID-19, of whom 65% were White, 26% were Black, and 9% were Hispanic. There were 16,317 (6.4%) individuals who tested positive for COVID-19 between February 8 and July 22, 2020, of whom 44% were White, 40% were Black, and 16% were Hispanic. While 66% of all individuals in care resided in urban areas, 76% of those tested and 87% of those testing positive for COVID-19 resided in urban areas. The geographic distribution of COVID-19 cases in the VA was similar to the pattern of known hotspots in the general population, including in the Northeast, South, and some Midwestern states (Fig 1). Several VA sites with the highest proportion of positive COVID-19 tests also performed a higher volume of tests and had the highest proportion of Black individuals in care, including New York City, New Orleans, and Chicago (Fig 1).

thumbnail
Fig 1. Distribution of 16,317 laboratory-confirmed COVID-19 cases in the US Department of Veterans Affairs as of July 22, 2020.

(a) Distribution of all COVID-19 laboratory-confirmed cases in the US Department of Veterans Affairs between February 8 and July 22, 2020, included in the current study. (b) Proportion of positive COVID-19 test results by the proportion of Black individuals in care by site. Map created using R library USMAP (v0.5.0) and RStudio (v3.6.3). COVID-19, coronavirus disease 2019.

https://doi.org/10.1371/journal.pmed.1003379.g001

thumbnail
Table 1. Characteristics of all individuals who were in care, were tested, tested positive, and died as of July 22, 2020.

https://doi.org/10.1371/journal.pmed.1003379.t001

Rate of testing and testing positive

Of the 254,595 patients tested for COVID-19, 73% received 1 test, 16% received 2 tests, 6% received 3 tests, and the remaining 5% received 4 or more tests. After reducing to 1 test per patient as described in Methods above, testing rates for COVID-19 were higher among Black (rate per 1,000 individuals: 60.0, 95% CI 59.6–60.5) and Hispanic individuals (52.7, 95% CI 52.1–53.4) compared to White individuals (38.6, 95% CI 38.4–38.7). Testing rates also varied by age, sex, rural/urban residence, region, and outbreak pattern (Table 2).

thumbnail
Table 2. COVID-19 testing by race/ethnicity among all individuals in care as of July 22, 2020.

https://doi.org/10.1371/journal.pmed.1003379.t002

Among individuals tested for COVID-19, the proportion with a positive test varied by race/ethnicity (Table 2); 4.4% (95% CI 4.3%–4.5%) of White, 10.2% (95% CI 10.0%–10.4%) of Black, and 11.4% (95% CI 11.0%–11.9%) of Hispanic individuals tested positive for COVID-19. For White and Black individuals, the proportion of positive COVID-19 tests was highest at ages under 60 years and at or over 80 years. For Hispanic individuals, the proportion of positive test results was highest among lower age groups (15.6%, 95% CI 14.6%–16.6%, for 20–39 years) and continuously decreased with increasing age.

Regression modeling of testing positive

Unadjusted associations between race/ethnicity and testing positive for COVID-19 yielded OR 2.49 (95% CI 2.40–2.58, p < 0.001) for Black and OR 2.80 (95% CI 2.67–2.94, p < 0.001) for Hispanic individuals compared to White individuals. After adjusting for age, the odds of testing positive did not change among Black individuals (OR 2.48, 95% CI 2.40–2.57, p < 0.001) and somewhat attenuated among Hispanic individuals (OR 2.56, 95% CI 2.44–2.69, p < 0.001) (Table 3). These associations further attenuated after additionally accounting for sex, rural/urban residence, and site of care among Black (OR 1.92, 95% CI 1.85–2.00, p < 0.001) and Hispanic individuals (OR 1.96, 95% CI 1.86–2.07, p < 0.001). These estimates were robust to any individual (S2 Table) or combined adjustment for comorbidities, substance use, and medication history. In fully adjusted models, Black (OR 1.93, 95% CI 1.85–2.01, p < 0.001) and Hispanic (OR 1.84, 95% CI 1.74–1.94, p < 0.001) individuals remained at increased odds of testing positive for COVID-19 (Fig 2).

thumbnail
Fig 2. Adjusted associations of demographic characteristics with testing positive for COVID-19 and subsequent 30-day mortality as of July 22, 2020.

(a) Positive test result among tested; (b) 30-day mortality among cases. Both models were conditioned on site of care and adjusted for baseline comorbidity (asthma, cancer, chronic kidney disease, chronic obstructive pulmonary disease, diabetes mellitus, hypertension, liver disease, vascular disease), substance use (alcohol consumption, alcohol use disorder, smoking status), and medication history (angiotensin converting enzyme inhibitor, angiotensin II receptor blocker). *Low number of mortality events in age groups 20–39 and 40–49 thus grouped with 50–59. CI, confidence interval; COVID-19, coronavirus disease 2019; OR, odds ratio.

https://doi.org/10.1371/journal.pmed.1003379.g002

thumbnail
Table 3. Associations with testing positive and subsequent 30-day mortality, February 8 to July 22, 2020.

https://doi.org/10.1371/journal.pmed.1003379.t003

The disparity in testing positive for COVID-19 between Black and White individuals decreased between the calendar periods February 8–April 21 (OR 2.16, 95% CI 1.98–2.36, p < 0.001) and June 22–July 22 (OR 1.74, 95% CI 1.64–1.85, p < 0.001) (Fig 3). By region, the disparity between Black and White individuals was highest in the Midwest (OR 2.66, 95% CI 2.41–2.95, p < 0.001) than any other region, and lowest in the West (OR 1.24, 95% CI 1.11–1.39, p < 0.001). By outbreak pattern, the disparity between Black and White individuals was highest at VA sites that experienced an early (OR 2.11, 95% CI 1.95–2.28, p < 0.001) or resurgent (OR 2.06, 95% CI 1.81–2.35, p < 0.001) outbreak and lowest at VA sites that experienced a late outbreak (OR 1.66, 95% CI 1.54–1.80, p < 0.001). In contrast, there was no variation observed in the disparity between Hispanic and White individuals by calendar time, region, or outbreak pattern.

thumbnail
Fig 3. Racial and ethnic disparities in testing positive for COVID-19, by calendar time, region, and outbreak pattern.

(a) Black versus White individuals; (b) Hispanic versus White individuals. All p < 0.001. Region based on US Census groupings. Outbreak pattern based on site-level percentage of positive tests per month among sites with at least 100 positive COVID-19 tests: early (≥10% in March or April), late (≥10% in June or July), resurgent (≥10% in March or April and June or July), steady (<10% in all months), other (sites with <100 positive tests). Models were conditioned on site of care and adjusted for baseline comorbidity (asthma, cancer, chronic kidney disease, chronic obstructive pulmonary disease, diabetes mellitus, hypertension, liver disease, vascular disease), substance use (alcohol consumption, alcohol use disorder, smoking status), and medication history (angiotensin converting enzyme inhibitor, angiotensin II receptor blocker). CI, confidence interval; COVID-19, coronavirus disease 2019; OR, odds ratio.

https://doi.org/10.1371/journal.pmed.1003379.g003

Regression modeling of 30-day mortality

There were 8,625 individuals who tested positive for COVID-19 on or before June 21, 2020, of whom 931 (457 [49%] White; 392 [42%] Black; and 82 [9%] Hispanic) died within 30 days. Unadjusted associations between race/ethnicity and mortality within 30 days of a positive test yielded OR 0.76 (95% CI 0.66–0.88, p < 0.001) for Black and OR 0.60 (95% CI 0.47–0.77, p < 0.001) for Hispanic individuals compared to White individuals. This association was not observed after adjusting for age among Black (OR 1.08, 95% CI 0.93–1.26, p = 0.32) and Hispanic individuals (OR 1.11, 95% CI 0.85–1.44, p = 0.46) (Table 3). The null association was robust to any further adjustment (Fig 2). Most deaths (n = 603, 65%) occurred in hospital. In sensitivity analyses, results from a model of 30-day mortality restricted to inpatient deaths did not alter the conclusions from the primary analyses (OR 1.12, 95% CI 0.90–1.40, p = 0.32, among Black individuals; OR 1.08, 95% CI 0.76–1.54, p = 0.65, among Hispanic individuals).

Discussion

This study examined racial and ethnic disparities in testing and subsequent COVID-19 mortality among approximately 6 million individuals receiving care in the US. We found that Black and Hispanic individuals were more likely to be tested and to test positive for COVID-19 than White individuals, even after comprehensive adjustment for underlying health conditions, other demographics, and geographic location. Among the variables assessed in this study, age, rural/urban residence, and site of care explained more of the racial/ethnic disparity in testing positive for COVID-19 than comorbidities, substance use, or medication history. While the disparity between Black and White individuals decreased over time, the disparity was strongest in the Midwest and at VA sites that experienced an early or resurgent outbreak. There was no variation observed in the disparity between Hispanic and White individuals by calendar time, region, or outbreak pattern. While individuals from minority backgrounds appeared to experience excess burden of COVID-19, among those infected, there was no observed difference in 30-day mortality by race/ethnicity group. The apparent racial/ethnic disparity in mortality in unadjusted data was principally explained by differing age structures between the populations.

Key strengths and limitations

This study elucidated racial and ethnic disparities in testing patterns of COVID-19 independent of underlying health status and other key factors in a nationwide sample. Strengths of this study included that it was based on well-annotated electronic health record data from a team with decades of experience using VA data, enabling a rapid and reliable analysis of COVID-19 outcomes by race and ethnicity. This analysis utilized patients’ records from an entire healthcare system, which made it less prone to collider bias (i.e., non-random selection of individuals into a study) than other COVID-19 studies limited to individuals testing positive or admitted to hospital [21]. Unlike other nationwide healthcare systems, linkage to COVID-19 testing data or outcomes was not required as the integrated nature of VA healthcare provided at over 1,200 sites allows all information to be stored in its Corporate Data Warehouse. We used validated algorithms to accurately extract information on and adjust models for a wide range of clinical, behavioral, and geographic factors, with very little missingness in the data. The scale of VA data also allowed us to assess the impact of COVID-19 separately across multiple racial and ethnic minority groups; combining or limiting analyses to a single minority group would have masked important differences between Black and Hispanic individuals. We continue to monitor COVID-19 outcomes for individuals of other minority backgrounds and plan to follow up these analyses when there are sufficient numbers for analysis.

While this analysis adds information, its limitations must be kept in mind. First, this study was conducted on veterans currently receiving care in the VA, who are older and have a higher prevalence of chronic health conditions and risk behaviors than the general US population [2224]. However, prior research has established that after adjusting for age, sex, race/ethnicity, region, and rural/urban residence, all of which were included in this study, there is no difference in total disease burden between veterans and non-veterans [24]. Our key finding of no observed racial disparity of COVID-19 mortality has also been shown in a smaller non-veteran population [9]; thus, associations reported in this study are likely generalizable to the wider US population. Second, while individuals in VA care represent a diversity of backgrounds, women represented a small proportion of individuals in the sample. Thus, our analysis was not powered to assess interactions between sex and race/ethnicity. Third, beyond adjusting for rural/urban location and site of testing, we were not able to explore likely social determinants of the pronounced differential burden of COVID-19 among minority individuals. More detailed information on nursing home residence and socioeconomic status (e.g., type of employment, income, number of individuals in household) were unavailable or not consistently recorded in VA data, as is the case in most other electronic health record data sources. Fourth, as is true outside the VA, only a small proportion of individuals have been tested (~5%), and rates of testing vary by site and within important subgroups. However, while initial testing was limited, by mid-April the VA began testing all individuals admitted to hospital and before any inpatient or outpatient procedures, even in those not suspected to have COVID-19. Our models for testing positive should be cautiously interpreted as a proxy of odds of infection since those with mild symptoms were unlikely to have received testing, particularly in the early stages of the outbreak.

Findings in context

Our findings of racial and ethnic disparities in COVID-19 provide important distinctions from previous reports in the US and other countries with ethnically diverse populations. To our knowledge, one of the largest studies to date on racial disparities in COVID-19 outcomes in the US followed 3,481 COVID-19 cases in the state of Louisiana and found that non-Hispanic Black individuals represented 77% of hospitalizations and 71% of deaths despite only making up 31% of the total source population [9]. However, this study was based on patients who tested positive for COVID-19 in a statewide healthcare system and was underpowered to investigate ethnic minorities. We were able to expand the scope of this finding nationally and to include Hispanic individuals. In the UK, which was the first country with a broadly ethnically diverse population to experience a COVID-19 outbreak [25], a study of 17 million individuals showed that those from minority backgrounds had a substantially higher risk of mortality from COVID-19, which was not fully explained by underlying health conditions or social deprivation [26]. While our study also found racial and ethnic disparities, we found that these disparities occurred primarily at a stage prior to hospitalization (i.e., testing positive for COVID-19). We found no evidence of racial or ethnic disparities in 30-day mortality once models were restricted to those who tested positive for COVID-19. Our findings may be an underestimate of the US population risk as health disparities in the VA tend to be smaller than in the private sector [27]. Nevertheless, at a population level the substantial excess burden of SARS-CoV-2 infection among Black and Hispanic individuals inevitably translates to excess COVID-19 mortality in these communities.

We demonstrated that Black and Hispanic individuals were more likely to test positive than their White counterparts even after accounting for underlying health conditions, other demographics, rural/urban residence, and site of care. Based on experience with the 1918 Spanish flu and the 2009 H1N1 epidemic, public health experts have warned that racial and ethnic minority populations may be at higher risk during infectious disease outbreaks due to underlying health conditions, lower access to care, and socioeconomic conditions [28,29]. Notably, our analysis found that underlying health conditions did not explain any of the disparity between racial/ethnic groups in the odds of testing positive for COVID-19 or subsequent mortality in models already accounting for demographics, principally age, rural/urban residence, and VA site of care—essential information to help guide effective interventions. Prior reports have also highlighted that members of racial and ethnic minorities are more likely to live in densely populated areas or multigenerational households, and minority groups are overrepresented in jails, prisons, and detention centers, all of which lead to reduced capacity to implement physical distancing [3034]. Similarly, Black and Hispanic workers are more likely than their White counterparts to be workers in essential industries, who continue to work outside the home despite outbreaks in their communities, making them more prone to exposure and therefore infection [3436].

We found substantial variation in the disparity between Black and White individuals in testing positive for COVID-19 by geographic region, with stronger disparity observed in the Midwest than all other regions, and disparity most attenuated in the West. Further breakdown of groups within the Black community (e.g., African American, Afro-Caribbean, African), which could potentially reveal additional variation, is not captured in VA data. The observed disparities may be due to differential social determinants of health between Black and White individuals across regions. A US Census Bureau report showed that while racial residential segregation has diminished over time nationally, communities in the Midwest remained less integrated than in the West [37]. If community-level exposure is driving risk of SARS-CoV-2 infections, then the disparity in testing positive for COVID-19 may be lower in regions with greater integration between White and Black residents, as is the case in the West. We also found that the disparity between Black and White individuals in testing positive slightly decreased over the study period and was highest at VA sites that experienced an early outbreak of COVID-19. This finding may be partially explained by the increased attention on racial disparities in COVID-19 in the media [35] that may have impacted behaviors like wearing face coverings in public to reduce the spread of infection [38].

Interestingly, the ethnic disparity between Hispanic and White individuals in testing positive for COVID-19 was consistent across time, geographic region, and outbreak pattern; the disparity was consistently observed across all strata. The lack of variation over time may be explained, in part, by less nationwide media coverage and epidemiological investigations of outbreaks of COVID-19 in Hispanic communities. Importantly, the Hispanic population in the US comprises a wide array of ethnic communities (e.g., Mexican, Puerto Rican, Cuban). However, these distinctions are not captured in VA data. The umbrella grouping may mask any potential variation within such a heterogeneous population. Further research on the impact of COVID-19 in Hispanic and Latinx communities is urgently needed.

Testing rates for COVID-19 in the VA were higher among Black and Hispanic individuals compared to White individuals. Local reporting from metropolitan areas with large minority populations, including New York City [39] and Chicago [40], has highlighted the disproportionate impact of COVID-19 in minority communities. We showed that VA facilities in these cities and others around the country that conducted the highest number of COVID-19 tests also had the highest proportion of Black individuals in care. There were also differences in the rate of COVID-19 testing and the proportion testing positive by age, sex, and type of residence across race/ethnicity groups. These findings demonstrate the need for epidemiological investigations to characterize testing patterns in the underlying population as they provide important context for interpretations of models of COVID-19 outcomes. To our knowledge, the largest medical record study to date analyzed COVID-19-related mortality in a population of 17 million residents in the UK, the vast majority of whom never tested for COVID-19 [26]. While the authors identified ethnic disparities in COVID-19-related mortality, the estimates reported can be interpreted as the overall burden of mortality by ethnicity without accounting for underlying testing patterns. We found a similar disparity in the overall burden of COVID-19-related death in the full source population of approximately 6 million individuals in care at the VA. However, when the model was restricted to individuals testing positive—which inherently accounts for factors related to access to testing, non-random testing, and odds of infection—racial and ethnic disparities in mortality were no longer observed.

Policy implications

These findings underscore the urgent need to proactively tailor strategies to contain and prevent further outbreaks in the US, principally focused on testing and getting individuals into care. Black and Hispanic communities are at increased risk of infection, justifying increased intensity of intervention. Our findings of variation in disparities over time and across geographic regions highlight the important need for community-based interventions at a state and local level to contain further exposure and outbreaks of COVID-19, particularly tailored to minority communities. Other interventions may include clinical decision support tools to prompt educational and testing interventions based upon an individualized risk assessment of testing positive for COVID-19. Outreach products about COVID-19 testing and disparities should also be distributed to patient advocates and groups at all points of care.

Future research

We appeal to other researchers investigating racial and ethnic disparities to perform analyses on the entire population at risk for COVID-19 where data are available, and to compare findings associated at each stage in the clinical course of COVID-19, from testing to outcomes. In this paper, we focused only on 30-day mortality among COVID-19 cases. We plan to explore other outcomes, including hospitalization, intensive care, and intubation, in subsequent analyses to examine whether racial and ethnic disparities exist in the clinical course of COVID-19 after testing positive and before death. Among other factors, future research should consider the role of other social determinants of health, including employment type, number of individuals in household, nursing home residence, and incarceration. Other racial and ethnic minorities in the US deserve attention, and while we did not have enough statistical power to include other groups in this analysis, we will continue to monitor these numbers for future research.

Supporting information

S1 Table. International Classification of Diseases–10th Revision, Clinical Modification (ICD-10-CM) diagnosis codes.

https://doi.org/10.1371/journal.pmed.1003379.s002

(XLSX)

S2 Table. Individual adjustments for the association between race/ethnicity and COVID-19 positivity and mortality.

https://doi.org/10.1371/journal.pmed.1003379.s003

(DOCX)

Acknowledgments

The authors wish to recognize Dr. Kendall Bryant as the NIAAA Scientific Collaborator for this study.

The views and opinions expressed in this paper are those of the authors and do not necessarily represent those of the Department of Veterans Affairs or the United States Government.


Transparency declaration

CTR and ACJ affirm that this paper 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 have been explained.

References

  1. 1. Center for Systems Science and Engineering. COVID-19 dashboard by the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University (JHU). Baltimore: Johns Hopkins University; 2020 [cited 2020 Aug 31]. Available from: https://coronavirus.jhu.edu/map.html.
  2. 2. Vahidy FS, Nicolas JC, Meeks JR, Khan O, Pan A, Jones SL, et al. Racial and ethnic disparities in SARS-CoV-2 pandemic: analysis of a COVID-19 observational registry for a diverse US metropolitan population. BMJ Open. 2020;10:e039849. pmid:32784264
  3. 3. Reyes C, Husain N, Gutowski C, St. Clair S, Pratt G. Chicago’s coronavirus disparity: black Chicagoans are dying at nearly six times the rate of white residents, data show. Chicago Tribune. 2020 Apr 7.
  4. 4. Thebault R, Tran A, Williams V. The coronavirus is infecting and killing black Americans at an alarmingly high rate. Washington Post. 2020 Apr 7.
  5. 5. Deslatte M. Louisiana data: virus hits blacks, people with hypertension. AP News. 2020 Apr 7.
  6. 6. Garg S, Kim L, Whitaker M, O’Halloran A, Cummings C, Holstein R, et al. Hospitalization rates and characteristics of patients hospitalized with laboratory-confirmed coronavirus disease 2019—COVID-NET, 14 states, March 1–30, 2020. MMWR Morb Mortal Wkly Rep. 2020;69(15):458–64. pmid:32298251
  7. 7. Rast J, Martinez YC, Williams LH. Milwaukee’s coronavirus racial divide: a report on the early stages of COVID-19 spread in Milwaukee County 2020. Milwaukee: Center for Economic Development; 2020 [cited 2020 Aug 31]. Available from: https://dc.uwm.edu/cgi/viewcontent.cgi?article=1053&context=ced_pubs.
  8. 8. Rodriguez-Diaz CE, Guilamo-Ramos V, Mena L, Hall E, Honermann B, Crowley JS, et al. Risk for COVID-19 infection and death among Latinos in the United States: examining heterogeneity in transmission dynamics. Ann Epidemiol. 2020 Jul 22. pmid:32711053
  9. 9. Price-Haywood EG, Burton J, Fort D, Seoane L. Hospitalization and mortality among black patients and white patients with Covid-19. N Engl J Med. 2020;382:2534–43. pmid:32459916
  10. 10. Pareek M, Bangash MN, Pareek N, Pan D, Sze S, Minhas JS, et al. Ethnicity and COVID-19: an urgent public health research priority. Lancet. 2020;395:1421–2. pmid:32330427
  11. 11. Zhou F, Yu T, Du R, Fan G, Liu Y, Liu Z, et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. Lancet. 2020;395(10229):1054–62. pmid:32171076
  12. 12. Grasselli G, Zangrillo A, Zanella A, Antonelli M, Cabrini L, Castelli A, et al. Baseline characteristics and outcomes of 1591 patients infected with SARS-CoV-2 admitted to ICUs of the Lombardy Region, Italy. JAMA. 2020;323(16):1574–81. pmid:32250385
  13. 13. CDC COVID-19 Response Team. Preliminary estimates of the prevalence of selected underlying health conditions among patients with Coronavirus Disease 2019—United States, February 12–March 28, 2020. MMWR Morb Mortal Wkly Rep. 69(13):382–6. pmid:32240123
  14. 14. Bhatraju PK, Ghassemieh BJ, Nichols M, Kim R, Jerome KR, Nalla AK, et al. Covid-19 in critically ill patients in the Seattle region—case series. N Engl J Med. 2020;382:2012–22. pmid:32227758
  15. 15. Wu Z, McGoogan JM. Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) outbreak in China: summary of a report of 72314 cases from the Chinese Center for Disease Control and Prevention. JAMA. 2020;323(13):1239–42. pmid:32091533
  16. 16. Rentsch CT, Kidwai-Khan F, Tate JP, Park LS, King JT, Skanderson M, et al. Covid-19 testing, hospital admission, and intensive care among 2,026,227 United States veterans aged 54–75 years. medRxiv. 2020 Apr 14. pmid:32511595
  17. 17. Abrams TE, Vaughan-Sarrazin M, Kaboli PJ. Mortality and revascularization following admission for acute myocardial infarction: implication for rural veterans. J Rural Health. 2010;26(4):310–7. pmid:21029165
  18. 18. McGinnis KA, Brandt CA, Skanderson M, Justice AC, Shahrir S, Butt AA, et al. Validating smoking data from the Veteran’s Affairs Health Factors dataset, an electronic data source. Nicotine Tob Res. 2011;13(12):1233–9. pmid:21911825
  19. 19. Fiellin DA, Reid MC, O’Connor PG. Screening for alcohol problems in primary care: a systematic review. Arch Intern Med. 2000;160(13):1977–89. pmid:10888972
  20. 20. United States Census Bureau. 2010 Census regions and divisions of the United States. Suitland (MD): United States Census Bureau; 2018 [cited 2020 Jul 30]. Available from: https://www.census.gov/geographies/reference-maps/2010/geo/2010-census-regions-and-divisions-of-the-united-states.html.
  21. 21. Griffith G, Morris TM, Tudball M, Herbert A, Mancano G, Pike L, et al. Collider bias undermines our understanding of COVID-19 disease risk and severity. medRxiv. 2020 May 20.
  22. 22. Hoerster KD, Lehavot K, Simpson T, McFall M, Reiber G, Nelson KM. Health and health behavior differences: U.S. Military, veteran, and civilian men. Am J Prev Med. 2012;43(5):483–9. pmid:23079170
  23. 23. Lehavot K, Hoerster KD, Nelson KM, Jakupcak M, Simpson TL. Health indicators for military, veteran, and civilian women. Am J Prev Med. 2012;42(5):473–80. pmid:22516487
  24. 24. RAND Health. Current and projected characteristics and unique health care needs of the patient population served by the Department of Veterans Affairs. Santa Monica (CA): RAND Corporation; 2015.
  25. 25. Khunti K, Singh AK, Pareek M, Hanif W. Is ethnicity linked to incidence or outcomes of covid-19? BMJ. 2020;369:m1548. pmid:32312785
  26. 26. Williamson EJ, Walker AJ, Bhaskaran K, Bacon S, Bates C, Morton CE, et al. Factors associated with COVID-19-related death using OpenSAFELY. Nature. 2020;584:430–6. pmid:32640463
  27. 27. Peterson K, Anderson J, Boundy E, Ferguson L, McCleery E, Waldrip K. Mortality disparities in racial/ethnic minority groups in the Veterans Health Administration: an evidence review and map. Am J Public Health. 2018;108(3):e1–11. pmid:29412713
  28. 28. Centers for Disease Control and Prevention. Flu-related hospitalizations and deaths in the United States from April 2009–January 30, 2010. Atlanta: Centers for Disease Control and Prevention; 2010 [cited 2020 Apr 2]. Available from: https://www.cdc.gov/H1N1flu/hosp_deaths_ahdra.htm.
  29. 29. Hutchins SS, Fiscella K, Levine RS, Ompad DC, McDonald M. Protection of racial/ethnic minority populations during an influenza pandemic. Am J Public Health. 2009;99(Suppl 2):S261–70. pmid:19797739
  30. 30. Fothergill A, Maestas EG, Darlington JD. Race, ethnicity and disasters in the United States: a review of the literature. Disasters. 1999;23(2):156–73. pmid:10379098
  31. 31. Thomas D, Phillips B, Lovekamp W, Fothergill A. Social vulnerability to disasters. Boca Raton (FL): Taylor & Francis Group; 2013.
  32. 32. Bolin B, Kurtz LC. Race, class, ethnicity, and disaster vulnerability. In: Rodríguez H, Donner W, Trainor J, editors. Handbook of disaster research. Berlin: Springer; 2018. pp. 181–203.
  33. 33. Resnick A, Galea S, Sivashanker K. Covid-19: the painful price of ignoring health inequities. BMJ Opinion. 2020 Mar 18 [cited 2020 Jul 27]. Available from: https://blogs.bmj.com/bmj/2020/03/18/covid-19-the-painful-price-of-ignoring-health-inequities/.
  34. 34. Hooper MW, Nápoles AM, Pérez-Stable EJ. COVID-19 and racial/ethnic disparities. JAMA. 2020;323(24):2466–7. pmid:32391864
  35. 35. Centers for Disease Control and Prevention. Health equity considerations & racial & ethnic minority groups. Atlanta: Centers for Disease Control and Prevention; 2020 [cited 2020 Sep 2]. Available from: https://www.cdc.gov/coronavirus/2019-ncov/need-extra-precautions/racial-ethnic-minorities.html.
  36. 36. Berchick E, Barnett J, Upton R. Health insurance coverage in the United States: 2018. Current Population Reports, P60-267(RV). Suitland (MD): United States Census Bureau; 2019.
  37. 37. United States Census Bureau. Racial and ethnic residential segregation in the United States: 1980–2000. Suitland (MD): United States Census Bureau; 2002 [cited 2020 Jul 31]. Available from: https://www.census.gov/prod/2002pubs/censr-3.pdf.
  38. 38. Fisher KA, Barile JP, Guerin RJ, Vanden Esschert KL, Jeffers A, Tian LH, et al. Factors associated with cloth face covering use among adults during the COVID-19 pandemic—United States, April and May 2020. MMWR Morb Mortal Wkly Rep. 2020;69(28):933–7. pmid:32673303
  39. 39. NYC Health. COVID-19: data. New York: NYC Health; 2020 [cited 2020 May 12]. Available from: https://www1.nyc.gov/site/doh/covid/covid-19-data.page.
  40. 40. Chicago Department of Public Health. Latest data. Chicago: City of Chicago; 2020 May 11 [cited 2020 May 12]. Available from: https://www.chicago.gov/city/en/sites/covid-19/home/latest-data.html.