Skip to main content
Advertisement
  • Loading metrics

Changes in the associations of race and rurality with SARS-CoV-2 infection, mortality, and case fatality in the United States from February 2020 to March 2021: A population-based cohort study

  • George N. Ioannou ,

    Roles Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Writing – original draft, Writing – review & editing

    georgei@medicine.washington.edu

    Affiliations Divisions of Gastroenterology, Veterans Affairs Puget Sound Healthcare System and University of Washington, Seattle, Washington, United States of America, Research and Development, Veterans Affairs Puget Sound Health Care System, Seattle, Washington, United States of America

  • Jacqueline M. Ferguson,

    Roles Conceptualization, Writing – review & editing

    Affiliation Center for Innovation to Implementation, VA Palo Alto Healthcare System, US Department of Veterans Affairs, Palo Alto, California, United States of America

  • Ann M. O’Hare,

    Roles Conceptualization, Investigation, Methodology, Writing – review & editing

    Affiliation Nephrology, Veterans Affairs Puget Sound Healthcare System and University of Washington, Seattle, Washington, United States of America

  • Amy S. B. Bohnert,

    Roles Conceptualization, Writing – review & editing

    Affiliation Department of Psychiatry, University of Michigan Medical School, Ann Arbor, Michigan, United States of America

  • Lisa I. Backus,

    Roles Conceptualization, Writing – review & editing

    Affiliation Department of Veterans Affairs, Population Health, Palo Alto Healthcare System, Palo Alto, California, United States of America

  • Edward J. Boyko,

    Roles Conceptualization, Writing – review & editing

    Affiliation General Internal Medicine, Veterans Affairs Puget Sound Healthcare System and University of Washington, Seattle, Washington, United States of America

  • Thomas F. Osborne,

    Roles Conceptualization, Writing – review & editing

    Affiliation Veterans Affairs Palo Alto Healthcare System, Palo Alto, and Department of Radiology, Stanford University School of Medicine, Stanford, California, United States of America

  • Matthew L. Maciejewski,

    Roles Writing – review & editing

    Affiliations Center of Innovation to Accelerate Discovery and Practice Transformation (ADAPT), Durham Veterans Affairs Health Care System, Durham, North Carolina, United States of America, Department of Population Health Sciences, Duke University School of Medicine, Durham, North Carolina, United States of America, Duke-Margolis Center for Health Policy, Duke University School of Medicine, Durham, North Carolina, United States of America, Division of General Internal Medicine, Duke University School of Medicine, Durham, North Carolina, United States of America

  • C. Barrett Bowling,

    Roles Conceptualization, Writing – review & editing

    Affiliation Durham Veterans Affairs Geriatric Research Education and Clinical Center, Durham Veterans Affairs Medical Center (VAMC), Durham, NC and Department of Medicine, Duke University, Durham, North Carolina, United States of America

  • Denise M. Hynes,

    Roles Conceptualization, Writing – review & editing

    Affiliations Center of Innovation to Improve Veteran Involvement in Care, VA Portland Healthcare System, Portland, Oregon, United States of America, Health Management and Policy, School of Social and Behavioral Health Sciences, College of Public Health and Human Sciences, Health Data and Informatics Program, Center for Genome Research and Biocomputing, Oregon State University, Corvallis, Oregon, United States of America

  • Theodore J. Iwashyna,

    Roles Conceptualization, Writing – review & editing

    Affiliations Center for Clinical Management Research, VA Ann Arbor Health System, Ann Arbor, Michigan, United States of America, Division of Pulmonary and Critical Care Medicine, Department of Internal Medicine, University of Michigan, Ann Arbor, Michigan, United States of America

  • Melody Saysana,

    Roles Conceptualization, Project administration, Writing – review & editing

    Affiliation Research and Development, Veterans Affairs Puget Sound Health Care System, Seattle, Washington, United States of America

  • Pamela Green,

    Roles Conceptualization, Data curation, Formal analysis, Software

    Affiliation Research and Development, Veterans Affairs Puget Sound Health Care System, Seattle, Washington, United States of America

  • Kristin Berry

    Roles Conceptualization, Formal analysis, Methodology, Writing – original draft, Writing – review & editing

    Affiliation Research and Development, Veterans Affairs Puget Sound Health Care System, Seattle, Washington, United States of America

Abstract

Background

We examined whether key sociodemographic and clinical risk factors for Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infection and mortality changed over time in a population-based cohort study.

Methods and findings

In a cohort of 9,127,673 persons enrolled in the United States Veterans Affairs (VA) healthcare system, we evaluated the independent associations of sociodemographic and clinical characteristics with SARS-CoV-2 infection (n = 216,046), SARS-CoV-2–related mortality (n = 10,230), and case fatality at monthly intervals between February 1, 2020 and March 31, 2021. VA enrollees had a mean age of 61 years (SD 17.7) and were predominantly male (90.9%) and White (64.5%), with 14.6% of Black race and 6.3% of Hispanic ethnicity. Black (versus White) race was strongly associated with SARS-CoV-2 infection (adjusted odds ratio [AOR] 5.10, [95% CI 4.65 to 5.59], p-value <0.001), mortality (AOR 3.85 [95% CI 3.30 to 4.50], p-value < 0.001), and case fatality (AOR 2.56, 95% CI 2.23 to 2.93, p-value < 0.001) in February to March 2020, but these associations were attenuated and not statistically significant by November 2020 for infection (AOR 1.03 [95% CI 1.00 to 1.07] p-value = 0.05) and mortality (AOR 1.08 [95% CI 0.96 to 1.20], p-value = 0.21) and were reversed for case fatality (AOR 0.86, 95% CI 0.78 to 0.95, p-value = 0.005). American Indian/Alaska Native (AI/AN versus White) race was associated with higher risk of SARS-CoV-2 infection in April and May 2020; this association declined over time and reversed by March 2021 (AOR 0.66 [95% CI 0.51 to 0.85] p-value = 0.004). Hispanic (versus non-Hispanic) ethnicity was associated with higher risk of SARS-CoV-2 infection and mortality during almost every time period, with no evidence of attenuation over time. Urban (versus rural) residence was associated with higher risk of infection (AOR 2.02, [95% CI 1.83 to 2.22], p-value < 0.001), mortality (AOR 2.48 [95% CI 2.08 to 2.96], p-value < 0.001), and case fatality (AOR 2.24, 95% CI 1.93 to 2.60, p-value < 0.001) in February to April 2020, but these associations attenuated over time and reversed by September 2020 (AOR 0.85, 95% CI 0.81 to 0.89, p-value < 0.001 for infection, AOR 0.72, 95% CI 0.62 to 0.83, p-value < 0.001 for mortality and AOR 0.81, 95% CI 0.71 to 0.93, p-value = 0.006 for case fatality). Throughout the observation period, high comorbidity burden, younger age, and obesity were consistently associated with infection, while high comorbidity burden, older age, and male sex were consistently associated with mortality. Limitations of the study include that changes over time in the associations of some risk factors may be affected by changes in the likelihood of testing for SARS-CoV-2 according to those risk factors; also, study results apply directly to VA enrollees who are predominantly male and have comprehensive healthcare and need to be confirmed in other populations.

Conclusions

In this study, we found that strongly positive associations of Black and AI/AN (versus White) race and urban (versus rural) residence with SARS-CoV-2 infection, mortality, and case fatality observed early in the pandemic were ameliorated or reversed by March 2021.

Author summary

Why was this study done?

  • As the Coronavirus Disease 2019 (COVID-19) pandemic continues to evolve, some risk factors for infection with COVID-19 and death due to COVID-19 that were described early in the pandemic may be changing.
  • Recognizing such changes is important in informing population-based approaches to prevent infection and reduce mortality.

What did the researchers do and find?

  • We investigated how the associations of key sociodemographic and clinical factors with COVID-19 infection, mortality, or case fatality changed between February 2020 and March 2021 among a cohort of approximately 9.1 million persons enrolled in the national US Veterans Affairs (VA) healthcare system, including 216,046 who tested positive and 10,230 who died of COVID-19 during the study period.
  • Black (versus White) race was strongly associated with a 5-fold higher risk of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infection, a 4-fold higher risk of mortality, and a 2.5-fold higher risk of case fatality in February to March 2020, but these associations attenuated over time and were no longer statistically significant by November 2020 for infection and mortality and were reversed for case fatality.
  • American Indian/Alaska Native (AI/AN versus White) race was associated with SARS-CoV-2 infection early in the pandemic, but this association declined over time and reversed by March 2021.
  • Urban (versus rural) residence was associated with 2-fold higher risk of infection, a 2.5-fold higher risk of mortality, and 2.2-fold higher risk of case fatality in February to April 2020, but these associations attenuated over time and reversed by September 2020.
  • Throughout the observation period, high comorbidity burden, younger age, Hispanic ethnicity, and obesity were consistently associated with infection, while high comorbidity burden, older age, Hispanic ethnicity, and male sex were consistently associated with mortality.

What do these findings mean?

  • Early in the pandemic, there were strongly positive associations of Black and AI/AN (versus White) race and urban (versus rural) residence with SARS-CoV-2 infection, mortality, and case fatality, but these were ameliorated or even reversed by March 2021.
  • Our results apply directly to VA enrollees who are predominantly male and have access to universal healthcare; they need to be confirmed in other populations.

Introduction

Sociodemographic factors and comorbidity burden have emerged as major risk factors for Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infection and mortality [113]. Black, American Indian/Alaska Native (AI/AN), and Hispanic persons have been reported to have higher risk of SARS-CoV-2 infection and mortality than White and non-Hispanic persons [711,14]. Obese persons [12,15] and those with higher comorbidity burden [9] have also been reported to have higher risk of SARS-CoV-2 infection and mortality, while older age is one of the strongest risk factors for SARS-CoV-2–related mortality [2,9]. Over the course of the pandemic, residing in a geographical region with high incidence of SARS-CoV-2 infection at a given time has proven to be a strong risk factor for SARS-CoV-2 infection and mortality [2,5,16].

Three “waves” of the SARS-CoV-2 pandemic have been described in the US, with peaks in cases in April 2020 (first wave), July 2020 (second wave), and December 2020 to January 2021 (third wave). Risk factors for SARS-CoV-2–related infection and mortality may be changing over time since the pandemic began, especially in relation to these waves. Although it is clear that geographical regions with high infection, hospitalization, and mortality rates changed over time as the pandemic surged in different parts of the country, prior studies have not examined whether risk factors such as age, comorbidity burden, race, and ethnicity also varied. Over the course of the pandemic, there have been marked changes in use of prophylactic measures and access to care as well as viral characteristics (e.g., emergence of new variants) and availability of treatments (e.g., use of different pharmacotherapies). These changes may have affected the associations of sociodemographic and other risk factors with SARS-CoV-2 infection or mortality. Understanding changing patterns of risk factors could be important in informing population-based approaches to prevent infection and reduce mortality by targeting those at highest risk at any given time during the time course of an evolving pandemic. Also, changing patterns over time in sociodemographic risk factors for SARS-CoV-2 may provide insights that are more broadly applicable to disparities research in general.

We identified a cohort of approximately 9.1 million persons who were enrolled in the national US Department of Veterans Affairs (VA) healthcare system at the beginning of the pandemic in February 2020 and followed this cohort over the ensuing 14-month period to determine whether the magnitude and direction of the associations of established risk factors with SARS-CoV-2 infection, mortality, and case fatality changed over time.

Methods

Study population and data source

The VA supports the largest integrated national healthcare system in the US, providing care at 168 medical centers and 1,112 outpatient clinics throughout the country. We identified a cohort of all persons aged ≥18 years who were alive and enrolled in the VA healthcare system on February 1, 2020 (n = 9,127,673). We followed this cohort for the development of SARS-CoV-2 infection and SARS-CoV-2–related mortality until March 31, 2021.

The VA employs a comprehensive, nationwide electronic health record (EHR) system. EHR data from all VA facilities are transferred to the VA’s centralized relational database, the Corporate Data Warehouse (CDW), on a nightly basis to support research and clinical operations [17]. The CDW includes the “COVID-19 Shared Data Resource,” a set of analytic variables and datasets related to Coronavirus Disease 2019 (COVID-19) developed and maintained by the VA Informatics and Computing Infrastructure (VINCI) specifically to facilitate COVID-19 research and operations, which we used in combination with other CDW data.

The study was approved by the VA Puget Sound Institutional Review Board (IRB# 01885), which waived the requirement to obtain informed consent because this was a retrospective study of data from an existing database.

Definition of SARS-CoV-2 infection and SARS-CoV-2–related death

Cohort members who tested positive for SARS-CoV-2 within the VA system based on polymerase chain reaction (PCR) tests were defined as having infection. The earliest date of a documented positive test was taken as each patient’s date of infection. Patients who died of any cause within 30 days of infection were defined as having a SARS-CoV-2–related death consistent with prior studies [2,3]. Deaths occurring both within and outside the VA are comprehensively captured in CDW from a variety of VA and non-VA sources including VA inpatient files, VA Beneficiary Identification and Records Locator System (BIRLS), Social Security Administration (SSA) death files, and the Department of Defense [18]. Deaths occurring from February 1, 2020 to February 28, 2021 were ascertained from files that were updated on June 25, 2021 to allow time for deaths to be electronically recorded in CDW.

Study outcomes: Monthly SARS-CoV-2 infection rates, mortality rates, and case fatality rates

SARS-CoV-2 infection and mortality rates were calculated monthly as a proportion of all cohort members who were still alive and at risk on the first day of each month. For the monthly analysis of SARS-CoV-2 infection rates, we excluded persons who died of any cause or were infected with SARS-CoV-2 before the beginning of each monthly period. For the analysis of SARS-CoV-2 mortality rates, we excluded persons who died of any cause before the beginning of each monthly period or were infected with SARS-CoV-2 more than 30 days before the beginning of the month (since SARS-CoV-2–related mortality was defined as death within 30 days of infection). The aim of removing from our cohort persons who were no longer at risk at the beginning of each observation period was to avoid a spurious attenuation of risk factors over time due to “depletion” of at-risk persons.

Monthly case fatality rates were calculated as the proportion of the patients who tested positive each month who died of any cause within 30 days of the earliest positive test.

Sociodemographic factors and comorbidity burden

We ascertained the following 8 characteristics for all cohort members on the first day of each monthly observation period: age (categorized as shown in Table 1 or modeled as restricted cubic splines with 5 knots at ages 30, 49, 64, 73, and 88 years corresponding to 5, 27.5, 50, 72.5, and 95 percentiles as recommended [19]), sex, self-reported race (White, Black, Asian, AI/AN, Pacific Islander/Native Hawaiian [PI/NH], and Other) and ethnicity (Hispanic and non-Hispanic), urban versus rural residence (based on zip codes, using data from the VA Office of Rural Health [20], which uses the Secondary Rural–Urban Commuting Area [RUCA] for defining rurality), geographic location (divided into 10 standard US Federal Regions [21]), body mass index (BMI, categorized using the World Health Organization groups [22] as shown in Table 1 or modeled as restricted cubic splines with 5 knots at BMIs of 21.3, 26.0, 29.0, 32.5, and 39.9 kg/m2 corresponding to 5, 27.5, 50, 72.5, and 95 percentiles), and Charlson comorbidity index (CCI) [23] calculated using the Deyo modification [24], which takes into account 19 comorbid conditions reported in the 2 years prior to each observation period. We focused on these sociodemographic factors, obesity, and the CCI because they are some of the most important risk factors for SARS-CoV-2 infection or mortality reported in the literature [113] and because trends in the associations over time could be plausibly hypothesized.

thumbnail
Table 1. Cohort characteristics and incidence of SARS-CoV-2 infection presented by month in a cohort of 9.1 million VA enrollees followed from February 2020 to March 2021.

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

Missing values for BMI (20.9%) were multiply imputed using values of the other covariates included in multivariable analyses. “Missing” values for race (18.7%) or ethnicity (15.6%) included persons who refused to declare their race/ethnicity or reported unknown or mixed race/ethnicity and did not self-identify as belonging to one of the prespecified racial or ethnic group. For these persons, we did not perform imputation of race/ethnicity, but rather included them in a “missing/unknown/refused” category.

Statistical analysis

We calculated the number of new SARS-CoV-2 infections and deaths each month and the monthly incidence as a proportion of all persons in our cohort who were still alive and at risk on the first day of the month. Monthly case fatality was calculated as the proportion of patients who tested positive each month who died within 30 days. For each monthly observation period (for infections) or bimonthly observation period (for mortality and case fatality), we used a separate multivariable logistic regression model to simultaneously adjust for the 8 characteristics listed above to determine the associations between each risk factor and SARS-CoV-2 infection, mortality, or case fatality reported as an adjusted odds ratio (AOR). We used bimonthly periods for trends in mortality and case fatality because there were too few deaths in some subgroups to generate reasonably precise monthly analyses of time trends.

We used a Wald test with cluster-robust standard errors to formally evaluate whether the associations between risk factors and each outcome (infection, mortality, or case fatality) changed over time by creating another model that combined all time periods and included an interaction term separately for each risk factor (risk factor * time period) where time period was an ordinal variable (see S1 and S2 Tables).

Almost all analyses were proposed a priori as shown in the S1 Analytic Plan. Notable exceptions were the analyses of case fatality as an outcome and additionally modeling age and BMI as restricted cubic splines that were performed in response to reviewers’ comments.

We followed the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) reporting guidelines as summarized in S1 STROBE Checklist.

Results

Baseline characteristics of VA enrollees

The baseline characteristics of the cohort of VA enrollees (n = 9,127,673) in February 2020 are shown in Table 1. Mean and median age were 61.0 and 64 years, respectively, with a substantial proportion in categories 65 to 74 (26.7%), 75 to 84 (13.8%), and ≥85 (8.4%) years of age. Most cohort members were male (90.9%) and White (64.5%), with 14.6% of Black race and 6.3% of Hispanic ethnicity. Cohort members were distributed almost evenly between urban versus rural locations, and the cohort included Veterans residing in all US Federal Regions and states, with the greatest contribution from region 4 (Southeast US, states of AL, FL, GA, KY, MS, NC, SC, and TN). Overweight (36.3%) and obesity (42.8%) were common among cohort members and 33.6% had a CCI ≥2.

Trends over time in the associations of risk factors with SARS-CoV-2 infection

During the entire follow-up period (February 1, 2020 to March 31, 2021), 216,046 out of 9,127,673 VA enrollees (2.4%) in our cohort tested positive for SARS-CoV-2 (Table 1). Infection rates peaked in accordance with the 3 well-described national waves of the pandemic in April 2020 (first wave), July 2020 (second wave), and December 2020 (third wave) as shown in Fig 1 (for infection rates), Fig 2 (for mortality rates), and Fig 3 (for case fatality rates). Characteristics independently associated with testing positive over the entire period included Black (versus White) race (AOR 1.39, 95% CI 1.37 to 1.41, p-value = <0.001), Hispanic ethnicity (AOR 1.64, 95% CI 1.62 to 1.67, p-value = <0.001), higher BMI (e.g., BMI 40 versus BMI 25 kg/m2, AOR 1.51, 95% CI 1.49 to 1.53, p-value < 0.001), and higher CCI (e.g., CCI ≥ 6 versus CCI = 0 to 1, AOR 3.16, 95% CI 3.12 to 3.21, p-value < 0.001) (Table 2). The highest risk of infection was observed at age 45 with progressively lower risk at both older ages (e.g., 85 versus 45 years old, AOR 0.50, 95% CI 0.43 to 0.58, p-value < 0.001) and younger ages (e.g., 25 versus 45 years old, AOR 0.77, 95% CI 0.76 to 0.79, p-value < 0.001). Men had significantly lower risk of infection than women (AOR 0.88, 95% CI 0.87 to 0.90, p-value = <0.001). Compared to Federal Region 4 (Southeast US), which served as the reference category, some regions were associated with significantly higher risk of positive SARS-CoV-2 test (e.g., Federal Regions 5 to 9), while others were associated with significantly lower risk (Federal Regions 1 to 3 and 10).

thumbnail
Fig 1. Monthly results in the VA healthcare system from February 2020 to March 2021 of SARS-CoV-2 infection rates among a cohort of VA enrollees.

SARS-CoV-2, Severe Acute Respiratory Syndrome Coronavirus 2; VA, Veterans Affairs.

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

thumbnail
Fig 2. Monthly results in the VA healthcare system from February 2020 to March 2021 of SARS-CoV-2 mortality rates among a cohort of VA enrollees.

SARS-CoV-2, Severe Acute Respiratory Syndrome Coronavirus 2; VA, Veterans Affairs.

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

thumbnail
Fig 3. Monthly results in the VA healthcare system from February 2020 to March 2021 of SARS-CoV-2 case fatality rates among VA enrollees testing positive for SARS-CoV-2.

SARS-CoV-2, Severe Acute Respiratory Syndrome Coronavirus 2; VA, Veterans Affairs.

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

thumbnail
Table 2. Trends over time in the associations (AORs*) of sociodemographic characteristics and comorbidity burden with risk of SARS-CoV-2 infection in a cohort of 9.1 million VA enrollees from February 2020 to March 2021.

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

The magnitude of the association between Black (versus White) race and SARS-CoV-2 infection declined steadily from February/March 2020 (AOR 5.10, 95% CI 4.65 to 5.59, p-value < 0.001) to November 2020 (AOR 1.03, 95% CI 1.00 to 1.07, p-value = 0.05) when it was no longer significant (p-value < 0.001 for interaction term of [black race * time period] testing trends over time, see S1 and S2 Tables). However, during the last 4 months of the observation period from December 2020 to March 2021 (corresponding to the third wave of the pandemic), Black race was again significantly associated with infection with AORs ranging from 1.15 to 1.30, although the magnitude of this association was still much lower than in the early months of the pandemic (Table 2, Fig 4). When we categorized our cohort by age (<65 and ≥65 years), the associations and their trends over time between Black (versus White) race and SARS-CoV-2 infection were very similar for persons aged <65 and ≥65 years (S3 Table).

thumbnail
Fig 4.

Trends over time in the associations of the following factors with risk of SARS-CoV-2 infection and mortality: (A and B) Black versus White race. (C and D) Urban versus rural location. (E and F) CCI categories. (G and H) Age. CCI, Charlson comorbidity index; SARS-CoV-2, Severe Acute Respiratory Syndrome Coronavirus 2.

https://doi.org/10.1371/journal.pmed.1003807.g004

The magnitude of the association between AI/AN (versus White) race and SARS-CoV-2 infection declined steadily over time (p-value for trends over time <0.001; see S1 and S2 Tables) and shifted from a positive association during the period from February to June (AORs ranging from 1.13 to 1.54) to a negative association in March 2021 (AOR 0.66, 95% CI 0.51 to 0.85, p-value = 0.004).

The magnitude of the association between urban versus rural location also declined steadily over time (p-value for trends over time <0.001; see S1 and S2 Tables) and shifted from a positive association in February/March 2020 (AOR 2.02, 95% CI 1.83 to 2.22, p-value = <0.001) to a negative association in September to October 2020 and a nonsignificant association in March 2021 (AOR 0.98, 95% CI 0.94 to 1.02, p-value = 0.32) (Table 2, Fig 4).

The magnitude of the associations between CCI and SARS-CoV-2 infection attenuated early in the pandemic until July 2020 (p-value < 0.001) and then appeared to plateau between July 2020 and January 2021, before declining again from January to March 2021, after the introduction of vaccination. However, CCI was still strongly associated with infection even in March 2021. Geographical regions at higher risk of infection fluctuated over time reflecting surges in different parts of the country. For example, region 2 (NY, NJ, and PR), which represented the earliest epicenter of the pandemic in the US, had the highest risk of infection in February to March 2020, but one of the lowest risks of infection in July to August 2020 and approximately average risk by March 2021.

The magnitude of the association of the following characteristics with SARS-CoV-2 infection did not vary appreciably over the observation period: age, sex, PI/NH race, Hispanic ethnicity, and BMI.

Trends over time in the associations of risk factors with SARS-CoV-2–related mortality

From February 2020 to March 2021, 10,230 SARS-CoV-2–related deaths were identified among our cohort of 9,127,673 VA enrollees (Table 3). Monthly mortality rates are shown in Table 3 and Fig 2. Significant, independent risk factors for mortality over the entire observation period included male sex (AOR 1.59, 95% CI 1.38 to 1.82, p-value < 0.001), older age (e.g., AOR 20.37, 95% CI 17.48 to 23.74, p-value < 0.001, comparing ages 85 year versus 45 year (the reference age)), Black (AOR 1.55, 95% CI 1.47 to 1.64, p-value = <0.001), and AI/AN (AOR 1.66, 95% CI 1.39 to 1.98, p-value = <0.001) race relative to White race, Hispanic ethnicity (AOR 1.57, 95% CI 1.45 to 1.70, p-value = <0.001), higher BMI (e.g., AOR 1.48, 95% CI 1.40 to 1.57, p-value < 0.001, comparing BMI 40 kg/m2 to BMI 25 kg/m2), and higher CCI (AOR 9.46, 95% CI 8.89 to 10.07, p-value = <0.001 for CCI ≥6 relative to CCI 0 to 1) (Table 4).

thumbnail
Table 3. SARS-CoV-2–related mortality presented by month in a cohort of 9.1 million VA enrollees followed from February 2020 to February 2021.

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

thumbnail
Table 4. Trends over time in the associations (AORs*) of sociodemographic characteristics and comorbidity burden with risk of SARS-CoV-2–related mortality in a cohort of 9.1 million VA enrollees from February 2020 to February 2021.

https://doi.org/10.1371/journal.pmed.1003807.t004

The magnitude of the association between Black (versus White) race and SARS-CoV-2–related mortality declined steadily over time from February to April (AOR 3.85, 95% CI 3.30 to 4.50, p-value < 0.001) and became nonsignificant by September to October (AOR 1.16, 95% CI 0.95 to 1.43, p-value = 0.15) (Table 4, Fig 2); however, Black race was again positively associated with SARS-CoV-2–related mortality in January to February 2021 (AOR 1.28, 95% CI 1.16 to 1.42, p-value = <0.001), albeit at a much lower magnitude of association.

The magnitude of the association between urban versus rural location declined steadily over time (p-value for trends over time <0.001, see S1 and S2 Tables) and shifted from a positive association in February/April (AOR 2.48, 95% CI 2.08 to 2.96, p-value = <0.001) to a negative association in all the time periods after September 2020. Geographical regions with the highest risk of SARS-CoV-2–related mortality in February to April (Federal Regions 1 and 2) had the lowest risk of SARS-CoV-2–related mortality by September to October with some increase in risk thereafter.

Associations between sex, CCI, or age and SARS-CoV-2–related mortality appeared to be stable over time.

Trends over time in the associations of risk factors with SARS-CoV-2 case fatality

Among a total of 206,789 persons who tested positive for SARS-CoV-2 infection from February 1, 2020 to February 28, 2021, 10,429 (5.0%) SARS-CoV-2–related deaths occurred, i.e., deaths within 30 days of infection (Table 5). Case fatality declined progressively over time from 12.6% in persons who tested positive in February to March 2020 to 3.1% in persons who tested positive in February 2021 (Table 5, Fig 3). Significant, independent risk factors for case fatality over the entire observation period included male sex (AOR 1.73, 95% CI 1.50 to 2.00, p-value < 0.001), older age (e.g., AOR 50.09, 95% CI 42.87 to 58.52, p-value < 0.001, for age 85 versus age 45 years), Black (AOR 1.13, 95% CI 1.06 to 1.19, p-value < 0.001), and AI/AN (AOR 1.76, 95% CI 1.44 to 2.14, p-value < 0.001) race relative to White race, Hispanic versus non-Hispanic ethnicity (AOR 1.21, 95% CI 1.11 to 1.31, p-value < 0.001), obesity/high BMI (e.g., AOR 1.35, 95% CI 1.27 to 1.43, p-value < 0.001 comparing BMI 40 versus 25 kg/m2) or underweight/low BMI (e.g., AOR 1.44, 95% CI 1.35 to 1.53, p-value < 0.001 comparing BMI 20 versus 25 kg/m2), and progressively higher CCI (e.g., AOR 2.51, 95% CI 2.35 to 2.68, p-value = <0.001 for CCI ≥6 relative to CCI 0 to 1) (Table 6).

thumbnail
Table 5. SARS-CoV-2 case fatality presented for VA enrollees who tested positive for SARS-CoV-2 each monthly period from February 2020 to February 2021.

https://doi.org/10.1371/journal.pmed.1003807.t005

thumbnail
Table 6. Trends over time in the associations (AORs*) of sociodemographic characteristics and comorbidity burden with risk of SARS-CoV-2 case fatality among 206,789 VA enrollees who tested positive for SARS-CoV-2 from February 1, 2020 to February 28, 2021, presented overall and divided into time periods.

https://doi.org/10.1371/journal.pmed.1003807.t006

The magnitude of the association between Black (versus White) race and case fatality declined over time and shifted from a positive association in February to April (AOR 2.56, 95% CI 2.23 to 2.93, p-value < 0.001) to a negative association in all time periods after September 2020 (Table 6). Also, the magnitude of the association between urban versus rural location and case fatality declined steadily over time and shifted from a positive association in February/April (AOR 2.24, 95% CI 1.93 to 2.60, p-value = <0.001) to a negative association in all the time periods after September 2020. Geographical regions with some of the highest case fatality rates in February to April (e.g., Federal Regions 1, 2, and 5) had some of the lowest case fatality rates by January/February 2021. Associations between sex, CCI, or age and case fatality appeared to be stable over time.

Discussion

Our study of a national cohort of 9.1 million VA enrollees followed from February 2020 to March 2021 demonstrates that the strongly positive associations between Black (versus White) race and SARS-CoV-2 infection, mortality, and case fatality that were observed in the early months of the pandemic attenuated over time and were no longer significant by November 2020. Positive associations between AI/AN (versus White) race and risk of infection noted early in the pandemic also attenuated over time and reversed by March 2021. Similarly, strongly positive associations between urban (versus rural) location and SARS-CoV-2 infection, mortality, and case fatality that were present early in the pandemic attenuated over time and were no longer significant by March 2021. Throughout the observation period, high comorbidity burden, younger age, Hispanic ethnicity, and obesity were consistently associated with infection, while high comorbidity burden, older age, Hispanic ethnicity, obesity, and male sex were consistently associated with mortality.

Multiple studies reported higher risk of SARS-CoV-2 infection and mortality in Black versus White persons [8,9,14]. In a nationally representative US study, Black persons accounted for 18.7% of SARS-CoV-2–related deaths from May to August 2020 despite representing just 12.5% of the US population. During this time period, the percentage of decedents who were Black decreased from 20.3% in May to 17.4% in August but was still higher than the percentage of Black persons in the US population [8]. Our study extends these prior findings in significant ways. First, we demonstrate that the association between Black race and SARS-CoV-2–related mortality continued to decline in the VA healthcare system even after August and became nonsignificant by November 2020. Second, we simultaneously adjusted for comorbidity burden, age, sex, ethnicity, BMI, geographic region, and rural/urban location to identify the associations of Black race and mortality that were not explained by differences in these factors. Third, we show that the associations of Black race with infection also declined over time and became nonsignificant by November 2020, while the positive associations with case fatality attenuated over time and actually reversed after September 2020. However, we also show that during the third wave of the pandemic, from December 2020 to March 2021, Black race was again associated with higher risk of SARS-CoV-2 infection and mortality, albeit at a much lower level than in the early period of the pandemic.

It is unclear why such a dramatic reduction in risk of infection in Black persons relative to White persons occurred in the space of 9 months between February and November 2020 [25,26]. Although the pandemic shifted from early urban epicenters with high percentages of Black persons to a broad, nationwide distribution in both urban and rural communities with lower percentages of Black persons, our results persisted even after adjustment for both urban/rural location and geographic region. Some factors that mediate racial inequities and contributed to increased risk in Black persons are unlikely to be reversed quickly, such as occupation (e.g., disproportionate representation in essential work settings with high exposure risk such as healthcare facilities, farms, factories, grocery stores, and public transportation), inability to work from home, housing (e.g., crowded housing, multigenerational households, and residence in densely populated neighborhoods with high infection rates), and reliance on public transportation. Factors that can change quickly include those related to prophylactic measures (masking, handwashing, physical distancing, and adhering to stay-at-home orders), disease awareness and behavior (e.g., limiting unnecessary exposure to crowds), and increased rates of testing (leading to early identification and reduction in transmission rates [27]).

Our comparisons of Black versus White VA enrollees need to be interpreted with caution when trying to extrapolate to differences by race in the US population as a whole. Black Veterans report higher median incomes (US$44,000 versus US$26,000), lower unemployment (3% versus 5%), and lower poverty levels (10% versus 21%) than non-Veteran Black adults[28]. However, Black VA enrollees have lower socioeconomic status than White VA enrollees including much higher unemployment (22.5% versus 6.8%) [29]. Levels of perceived discrimination in healthcare have been shown to be equally high among Black Veteran and non-Veteran groups in the US, despite access to comprehensive healthcare through the VA [30]. Furthermore, racial discrimination experienced during military service has been shown to be common among Black Veterans and associated with long-term impacts on self-reported physical health [31,32].

AI/AN persons had significant higher risk of infection than White persons early in the pandemic, but this association attenuated over time and actually reversed during the third wave of the pandemic between December 2020 and March 2021. However, SARS-CoV-2 mortality remained higher in AI/AN persons throughout the observation period.

Hispanic ethnicity was associated with increased risk of SARS-CoV-2 infection and mortality during almost every time period that we investigated. In contrast to Black race, the associations of Hispanic ethnicity with infection or mortality did not attenuate over time. This suggests that the factors listed above that mediate inequities in SARS-CoV-2 infection and mortality were not improved over time in Hispanic communities.

Big cities and metropolitan areas constituted the initial epicenters of SARS-CoV-2 infection in the US. However, the virus quickly spread throughout the US to both rural and urban areas, likely explaining the observed attenuation in the association between urban residence and SARS-CoV-2 infection and mortality over time. Unique challenges exist in both urban (dense population, housing distress, overcrowding, and public transport) and rural (geographic inaccessibility to SARS-CoV-2–related screening and care, higher disability, lack of social capital, and high-risk occupations such as meat and poultry processing [33]) that tend to raise risk of SARS-CoV-2 infection and mortality.

Comorbidity burden, estimated by the CCI, was one of the strongest risk factors for both SARS-CoV-2 infection and mortality. Furthermore, after May 2020, we did not observe any attenuation over time in the associations of high CCI with SARS-CoV-2–related mortality. It is possible that increased likelihood of testing persons with high comorbidity burden even if minimally symptomatic or asymptomatic may have contributed to the observed association with risk of infection. However, this does not explain the strong associations between high CCI and SARS-CoV-2–related mortality and case fatality that we observed. Most of the component comorbidities that make up the CCI have been individually associated with adverse COVID-19 outcomes, especially diabetes, congestive heart failure, chronic pulmonary disease, cerebrovascular disease, liver disease, kidney disease, and malignancy [1,3439]. The cumulative burden of these comorbidities appears to have a dramatic effect on SARS-CoV-2–related mortality. The associations between high CCI and risk of infection appear to decline steadily from January to March 2021. This may be related to initiation of vaccination that was restricted to high-risk groups in that time period.

Compared to age 45 years (the reference age in our analysis), both older and younger persons had significantly lower risk of testing positive in all time periods after adjusting for sex, race, ethnicity, geographical region, urban/rural location, BMI, and CCI. It is interesting that risk of testing positive in older persons declined even further in March 2021, which may be related to initiation of vaccination in older age groups during that time period. Similar associations between age ≥65 and lower risk of infection have been described with influenza virus [40] and are likely related to lower risk of exposure to SARS-CoV-2 in older persons. Despite the lower incidence of infection, older age was the strongest risk factor for SARS-CoV-2–related mortality (driven by much higher case fatality in older persons), an association that did not change in magnitude over time.

Our results should be interpreted in the context of some important limitations. It is impossible to identify all cases of SARS-CoV-2 infection in a population, since many are asymptomatic and even some symptomatic cases do not get confirmed by testing. We identified only the cases who were tested and identified as positive within the VA system (positive tests performed outside the VA were only captured if documented in the EHR). Therefore, changes in the associations of certain risk factors over time that we observed may also be caused by changes in likelihood of testing for SARS-CoV-2 and identifying positive cases. “Attenuation” over time in risk associated with a high-risk group may also occur if a sufficiently high proportion of that group has already acquired infection and has developed immunity. To avoid this problem, we removed from our cohort persons with known infection prior to the beginning of each monthly/bimonthly observation period. However, persons infected but not tested would still be retained in our cohort and might contribute to attenuation in high-risk groups. Although we adjusted for federal region and urban/rural location, it is possible that some residual confounding by geographic location persists [41]. We did not have data on educational attainment, poverty rates, or income. It is possible that race or ethnicity may be serving as a proxy for low socioeconomic status; the extent to which the associations of race or ethnicity are confounded by socioeconomic status should be examined in future research. Finally, our results apply to US VA enrollees who are predominantly male and have access to comprehensive healthcare. Therefore, they need to be confirmed in other populations. In particular, access to the VA’s comprehensive healthcare system has been shown to attenuate some racial and geographic disparities in adverse outcomes in other contexts [42].

In conclusion, strongly positive associations of Black race, AI/AN race, and urban residence with SARS-CoV-2 infection, mortality, and case fatality that were observed early in the pandemic attenuated over time in the VA system. On the other hand, other risk factors, such age, comorbidity burden, Hispanic ethnicity, and obesity, were strongly associated with infection and mortality throughout the observation period. Recognizing the potentially dynamic nature of associations between known risk factors for infection and mortality can help to inform ongoing population-based approaches to prevention and treatment of SARS-CoV-2 as well as providing insights for disparities research in other fields.

Disclaimers

The contents do not represent the views of the US Department of Veterans Affairs or the US Government.

Supporting information

S1 STROBE Checklist. STROBE Statement—Checklist of items that should be included in reports of cohort studies.

STROBE, Strengthening the Reporting of Observational Studies in Epidemiology.

https://doi.org/10.1371/journal.pmed.1003807.s001

(DOCX)

S1 Table. AOR* for interaction term of risk factor and monthly time period of SARS-CoV-2 infection treated as an ordinal variable (risk factor * time period) among 9.1 million VA enrollees from February 2020 to March 2021, performed as a test of trends over time.

AOR, adjusted odds ratio; SARS-CoV-2, Severe Acute Respiratory Syndrome Coronavirus 2; VA, Veterans Affairs.

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

(DOCX)

S2 Table. AOR* for interaction term of risk factor and monthly time period of SARS-CoV-2–related mortality treated as an ordinal variable (risk factor * time period) among 9.1 million VA enrollees from February 2020 to February 2021, performed as a test of trends over time.

AOR, adjusted odds ratio; SARS-CoV-2, Severe Acute Respiratory Syndrome Coronavirus 2; VA, Veterans Affairs.

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

(DOCX)

S3 Table. Trends over time in the associations (AORs*) of race with risk of SARS-CoV-2 infection presented separately for persons aged <65 versus ≥65 year among 9.1 million VA enrollees from February 2020 to March 2021.

AOR, adjusted odds ratio;SARS-CoV-2, Severe Acute Respiratory Syndrome Coronavirus 2; VA, Veterans Affairs.

https://doi.org/10.1371/journal.pmed.1003807.s004

(DOCX)

S1 Analytic Plan. Changes in risk factors for SARS-CoV-2 infection and mortality over time in a national healthcare system.

SARS-CoV-2, Severe Acute Respiratory Syndrome Coronavirus 2.

https://doi.org/10.1371/journal.pmed.1003807.s005

(DOCX)

References

  1. 1. 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(7821):430–6. Epub 2020/07/09. pmid:32640463.
  2. 2. Ioannou GN, Locke E, Green P, Berry K, O’Hare AM, Shah JA, et al. Risk Factors for Hospitalization, Mechanical Ventilation, or Death Among 10131 US Veterans With SARS-CoV-2 Infection. JAMA Netw Open. 2020;3(9):e2022310. Epub 2020/09/24. pmid:32965502; PubMed Central PMCID: PMC7512055.
  3. 3. King JT Jr, Yoon JS, Rentsch CT, Tate JP, Park LS, Kidwai-Khan F, et al. Development and validation of a 30-day mortality index based on pre-existing medical administrative data from 13,323 COVID-19 patients: The Veterans Health Administration COVID-19 (VACO) Index. PLoS ONE. 2020;15(11):e0241825. Epub 2020/11/12. pmid:33175863; PubMed Central PMCID: PMC7657526.
  4. 4. Rentsch CT, Kidwai-Khan F, Tate JP, Park LS, King JT Jr, Skanderson M, et al. Patterns of COVID-19 testing and mortality by race and ethnicity among United States veterans: A nationwide cohort study. PLoS Med. 2020;17(9):e1003379. Epub 2020/09/23. pmid:32960880; PubMed Central PMCID: PMC7508372.
  5. 5. Fan VS, Dominitz JA, Eastment MC, Locke E, Green P, Berry K, et al. Risk Factors for testing positive for SARS-CoV-2 in a national US healthcare system. Clin Infect Dis. 2020. Epub 2020/10/27. pmid:33105485; PubMed Central PMCID: PMC7665412.
  6. 6. Dai CL, Kornilov SA, Roper RT, Cohen-Cline H, Jade K, Smith B, et al. Characteristics and Factors Associated with COVID-19 Infection, Hospitalization, and Mortality Across Race and Ethnicity. Clin Infect Dis. 2021. Epub 2021/02/21. pmid:33608710; PubMed Central PMCID: PMC7929051.
  7. 7. Bassett MT, Chen JT, Krieger N. Variation in racial/ethnic disparities in COVID-19 mortality by age in the United States: A cross-sectional study. PLoS Med. 2020;17(10):e1003402. Epub 2020/10/21. pmid:33079941; PubMed Central PMCID: PMC7575091.
  8. 8. Gold JAW, Rossen LM, Ahmad FB, Sutton P, Li Z, Salvatore PP, et al. Race, Ethnicity, and Age Trends in Persons Who Died from COVID-19—United States, May-August 2020. MMWR Morb Mortal Wkly Rep. 2020;69(42):1517–21. Epub 2020/10/23. pmid:33090984; PubMed Central PMCID: PMC7583501.
  9. 9. Escobar GJ, Adams AS, Liu VX, Soltesz L, Chen YI, Parodi SM, et al. Racial Disparities in COVID-19 Testing and Outcomes: Retrospective Cohort Study in an Integrated Health System. Ann Intern Med. 2021. Epub 2021/02/09. pmid:33556278; PubMed Central PMCID: PMC7893537.
  10. 10. Mackey K, Ayers CK, Kondo KK, Saha S, Advani SM, Young S, et al. Racial and Ethnic Disparities in COVID-19-Related Infections, Hospitalizations, and Deaths: A Systematic Review. Ann Intern Med. 2021;174(3):362–73. Epub 2020/12/01. pmid:33253040; PubMed Central PMCID: PMC7772883.
  11. 11. 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(26):2534–43. Epub 2020/05/28. pmid:32459916; PubMed Central PMCID: PMC7269015.
  12. 12. Eastment MC, Berry K, Locke E, Green P, O’Hare A, Crothers K, et al. Body mass index (BMI) and outcomes of SARS-CoV-2 among US veterans. Obesity (Silver Spring). 2020. Epub 2020/12/19. pmid:33336934.
  13. 13. Ioannou GN, Green P, Fan VS, Dominitz JA, O’Hare AM, Backus LI, et al. Development of COVIDVax Model to Estimate the Risk of SARS-CoV-2-Related Death Among 7.6 Million US Veterans for Use in Vaccination Prioritization. JAMA Netw Open. 2021;4(4):e214347. Epub 2021/04/07. pmid:33822066.
  14. 14. Ferguson JM, Abdel Magid HS, Purnell AL, Kiang MV, Osborne TF. Differences in COVID-19 Testing and Test Positivity Among Veterans, United States, 2020. Public Health Rep. 2021:333549211009498. Epub 2021/04/08. pmid:33826875.
  15. 15. Tartof SY, Qian L, Hong V, Wei R, Nadjafi RF, Fischer H, et al. Obesity and Mortality Among Patients Diagnosed With COVID-19: Results From an Integrated Health Care Organization. Ann Intern Med. 2020;173(10):773–81. Epub 2020/08/14. pmid:32783686; PubMed Central PMCID: PMC7429998.
  16. 16. Asch DA, Sheils NE, Islam MN, Chen Y, Werner RM, Buresh J, et al. Variation in US Hospital Mortality Rates for Patients Admitted With COVID-19 During the First 6 Months of the Pandemic. JAMA Intern Med. 2020. Epub 2020/12/23. pmid:33351068; PubMed Central PMCID: PMC7756246.
  17. 17. Veterans Affairs Corporate Data Warehouse. [cited 2021 Mar 22]. Available from: http://www.hsrd.research.va.gov/for_researchers/vinci/cdw.cfm
  18. 18. Sohn MW, Arnold N, Maynard C, Hynes DM. Accuracy and completeness of mortality data in the Department of Veterans Affairs. Popul Health Metr. 2006;4:2. pmid:16606453; PubMed Central PMCID: PMC1458356.
  19. 19. Harrell FE Jr. Regression Modeling Strategies: With Applications to Linear Models, Logistic Regression, and Survival Analysis. New York: Springer; 2001.
  20. 20. Veterans Health Administration, Office of Rural Health. [cited 2021 Mar 28]. Available from: https://vaww.vashare.vha.va.gov/sites/ruralhealth/GSOD/default.aspx
  21. 21. Federal Emergency Management (FEMA) US Regions. [cited 2021 Feb 12]. Available from: https://www.fema.gov/about/organization/regions
  22. 22. World Health Organization. Body Mass Index. [cited 2021 Mar 21]. Available from: https://www.euro.who.int/en/health-topics/disease-prevention/nutrition/a-healthy-lifestyle/body-mass-index-bmi
  23. 23. Charlson ME, Pompei P, Ales KL, MacKenzie CR. A new method of classifying prognostic comorbidity in longitudinal studies: development and validation. J Chronic Dis. 1987;40(5):373–83. pmid:3558716
  24. 24. Deyo RA, Cherkin DC, Ciol MA. Adapting a clinical comorbidity index for use with ICD-9-CM administrative databases. J Clin Epidemiol. 1992;45(6):613–9. Epub 1992/06/01. [pii]. pmid:1607900.
  25. 25. Dorn AV, Cooney RE, Sabin ML. COVID-19 exacerbating inequalities in the US. Lancet. 2020;395(10232):1243–4. Epub 2020/04/20. pmid:32305087; PubMed Central PMCID: PMC7162639.
  26. 26. Baptiste DL, Commodore-Mensah Y, Alexander KA, Jacques K, Wilson PR, Akomah J, et al. COVID-19: Shedding light on racial and health inequities in the USA. J Clin Nurs. 2020;29(15–16):2734–6. Epub 2020/05/31. pmid:32472591; PubMed Central PMCID: PMC7300762.
  27. 27. Zhang M. Estimation of differential occupational risk of COVID-19 by comparing risk factors with case data by occupational group. Am J Ind Med. 2021;64(1):39–47. Epub 2020/11/20. pmid:33210336; PubMed Central PMCID: PMC7753309.
  28. 28. Profile of Veterans:2016, Data from the American Community Survey, Prepared by the National Center for Veterans Analysis and Statistics, US Department of Veterans Affairs, February 2018. [cited 2021 Aug 1]. Available from: https://www.va.gov/vetdata/veteran_population.asp
  29. 29. 2019 Survey of Veteran Enrollees’ Health and Use of Health Care, March 2020. [cited 2021 Aug 1]. Available from: https://www.va.gov/HEALTHPOLICYPLANNING/SOE2019/2019_Enrollee_Data_Findings_Report-March_2020_508_Compliant.pdf
  30. 30. Hausmann LR, Jeong K, Bost JE, Kressin NR, Ibrahim SA. Perceived racial discrimination in health care: a comparison of Veterans Affairs and other patients. Am J Public Health. 2009;99 Suppl 3:S718–24. Epub 2009/05/16. pmid:19443818; PubMed Central PMCID: PMC2774170.
  31. 31. Sohn L, Harada ND. Effects of racial/ethnic discrimination on the health status of minority veterans. Mil Med. 2008;173(4):331–8. Epub 2008/05/14. pmid:18472621.
  32. 32. Foynes MM, Smith BN, Shipherd JC. Associations between race-based and sex-based discrimination, health, and functioning: a longitudinal study of Marines. Med Care. 2015;53(4 Suppl 1):S128–35. Epub 2015/03/15. pmid:25767966.
  33. 33. Waltenburg MA, Victoroff T, Rose CE, Butterfield M, Jervis RH, Fedak KM, et al. Update: COVID-19 Among Workers in Meat and Poultry Processing Facilities—United States, April-May 2020. MMWR Morb Mortal Wkly Rep. 2020;69(27):887–92. Epub 2020/07/10. pmid:32644986; PubMed Central PMCID: PMC7732361.
  34. 34. Ioannou GN, Liang PS, Locke E, Green P, Berry K, O’Hare AM, et al. Cirrhosis and SARS-CoV-2 infection in US Veterans: risk of infection, hospitalization, ventilation and mortality. Hepatology. 2020. Epub 2020/11/22. pmid:33219546.
  35. 35. Izzi-Engbeaya C, Distaso W, Amin A, Yang W, Idowu O, Kenkre JS, et al. Adverse outcomes in COVID-19 and diabetes: a retrospective cohort study from three London teaching hospitals. BMJ Open Diabetes Res Care. 2021;9(1). Epub 2021/01/08. pmid:33408084; PubMed Central PMCID: PMC7789097.
  36. 36. Perez-Belmonte LM, Torres-Pena JD, Lopez-Carmona MD, Ayala-Gutierrez MM, Fuentes-Jimenez F, Huerta LJ, et al. Mortality and other adverse outcomes in patients with type 2 diabetes mellitus admitted for COVID-19 in association with glucose-lowering drugs: a nationwide cohort study. BMC Med. 2020;18(1):359. Epub 2020/11/17. pmid:33190637; PubMed Central PMCID: PMC7666969.
  37. 37. 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. Epub 2020/02/25. pmid:32091533.
  38. 38. Freaney PM, Shah SJ, Khan SS. COVID-19 and Heart Failure With Preserved Ejection Fraction. JAMA. 2020. Epub 2020/10/02. pmid:33001179.
  39. 39. Carlson N, Nelveg-Kristensen KE, Freese Ballegaard E, Feldt-Rasmussen B, Hornum M, Kamper AL, et al. Increased vulnerability to COVID-19 in chronic kidney disease. J Intern Med. 2021. Epub 2021/01/17. pmid:33452733; PubMed Central PMCID: PMC8014284.
  40. 40. Tokars JI, Olsen SJ, Reed C. Seasonal Incidence of Symptomatic Influenza in the United States. Clin Infect Dis. 2018;66(10):1511–8. Epub 2017/12/06. pmid:29206909; PubMed Central PMCID: PMC5934309.
  41. 41. Bilal U, Tabb LP, Barber S, Diez Roux AV. Spatial Inequities in COVID-19 Testing, Positivity, Confirmed Cases, and Mortality in 3 U.S. Cities: An Ecological Study. Ann Intern Med. 2021. Epub 2021/03/30. pmid:33780289.
  42. 42. 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–e11. Epub 2018/02/08. pmid:29412713; PubMed Central PMCID: PMC5803811.