Abstract
American Indian and Alaska Native (AI/AN) individuals are more likely to die with COVID-19 than other groups, but there is limited empirical evidence to explain the cause of this inequity. The objective of this study was to determine whether medical comorbidities, area socioeconomic deprivation, or access to treatment can explain the greater COVID-19 related mortality among AI/AN individuals. The design was a retrospective cohort study of harmonized electronic health record data of all inpatients with COVID-19 from 21 United States health systems from February 2020 through January 2022. The mortality of AI/AN inpatients was compared to all Non-Hispanic White (NHW) inpatients and to a matched subsample of NHW inpatients. AI/AN inpatients were more likely to die during their hospitalization (13.2% versus 7.1%; odds ratio [OR] = 1.98, 95% confidence interval [CI] = 1.48, 2.65) than their matched NHW counterparts. After adjusting for comorbidities, area social deprivation, and access to treatment, the association between ethnicity and mortality was substantially reduced (OR 1.59, 95% CI 1.15, 2.22). The significant residual relation between AI/AN versus NHW status and mortality indicate that there are other important unmeasured factors that contribute to this inequity. This will be an important direction for future research.
Introduction
Health inequities between American Indians and Alaska Natives (AI/AN) and other groups in the United States (US) have persisted throughout the 500 years of colonization1. Throughout US history, AI/AN Tribes and tribal communities have been repeatedly and disproportionately impacted by epidemics of diseases such as smallpox and measles1,2, and more recently, the H1N1 virus3. AI/AN individuals are more likely than non-Hispanic White (NHW) individuals to die from diseases associated with lower respiratory infection such as influenza and pneumonia4. Similarly, AI/AN individuals appear to be more severely affected by the COVID-19 pandemic than the general US population. Age-adjusted COVID-19 relative mortality risk of AI/ANs compared to NHW individuals range from 1.8 to 3.85,6,7,8. Death from COVID-19 has led to a widening of the life expectancy gap between AI/ANs and NHWs from 7 years in 2019 to 11 years in 20219. Despite declines in COVID-19 mortality from 2021 to 2022, racial and ethnic disparities persist10. Provisional 2022 mortality data indicate that COVID-19 was the fourth leading cause of death in the US, and COVID-19 death rates remained high for AI/AN individuals compared to NHW individuals10.
Few empirical investigations have explored the potential causes of the disparities in COVID-19 outcomes between AI/AN and NHW individuals. One explanation that has been proposed is that AI/AN individuals experience more medical comorbidities, such as diabetes, liver, and kidney disease than NHW individuals11,12, and some of these have been linked to severe COVID-19 outcomes13,14. To our knowledge, no study has empirically evaluated whether comorbidities partly explain the COVID-19 mortality disparity for AI/AN individuals. Two studies demonstrated that, at every level of comorbidity, differences between race groups (including AI/AN, Black, and White) persisted, suggesting that factors other than comorbidity contributed to the disparity in COVID-19 mortality15,16.
For instance, social and economic conditions in which AI/AN populations live (i.e., area socioeconomic deprivation) may drive inequities in infection and mortality12. Historically, US government policies, including forced removal from lands, perpetuated the severity of infectious epidemics by causing starvation, severe crowding, and “historical trauma”17. Today, due to these inequities, AI/AN individuals are more likely than NHW individuals to live in neighborhoods characterized by poverty and by housing that is overcrowded and without plumbing18,19.
Area socioeconomic deprivation is associated with higher risk for COVID-19 infection20, hospitalization and mortality21. Wong et al.19 examined the extent to which area socioeconomic deprivation contributed to AI/AN disparities in COVID-19 infection in a large geographically diverse sample extracted from the United States Veteran’s Health Administration electronic health records (including 3045 AI/AN US Veterans). The investigators demonstrated that 17–35% of the COVID-19 infection disparity between AI/AN and NHW Veterans could be explained by area socioeconomic deprivation, and this did not differ as a function of whether the Veteran lived on or near a reservation. Studies have not tested the hypothesis that area socioeconomic deprivation explains AI/AN disparities in COVID-19 mortality.
Access to health care is an additional mechanism of disparities in COVID-19 mortality between AI/ANs and NHWs2,18,22. AI/ANs often need to travel a great distance to the nearest health care facility and many lack adequate health insurance coverage18. To our knowledge, the hypothesis that poorer access to health care may explain disparities between AI/ANs and NHWs in COVID-19 mortality has never been directly tested.
Although most AI/AN individuals (87%) do not reside on reservations or tribal land23, most research on socioeconomic deprivation has focused on the 13% minority of AI/AN individuals who do19. The experiences of urban versus tribal residing AI/AN individuals differ substantially. For example, urban AI/AN individuals experience less continuity in their health care compared to AI/AN individuals who reside on tribal land24,25. There is a critical need for more research on COVID-19 outcomes among urban-dwelling AI/ANs.
The dearth of COVID-19 data on AI/AN individuals in state-level public health surveillance systems has been of grave concern25,26,27. COVID-19 tracking systems have failed to measure or report the race of those affected or to include AI/AN as a distinct category26,28. Lack of quality data may be contributing to a widening of health inequities for AI/ANs.
The present study assembled a sample of AI/AN individuals affected with COVID-19 by mining electronic health records (EHRs) of all inpatients affected with COVID-19 from 21 health systems over 2 years. There were 145,944 geographically diverse, predominantly urban patients (87% with known race) hospitalized with COVID-19—546 (0.37%) were AI/AN. The AI/AN inpatients were from 29 different states, which makes this one of the most geographically diverse samples of AI/AN individuals assembled to examine mortality disparities. Importantly, most of the AI/AN patients were living in urban areas, which represents the understudied majority of the AI/AN population19,29.
The goals of this study were to (1) explore the potential inequity in COVID-19 mortality in a primarily urban, geographically diverse sample of AI/AN, and to (2) examine whether comorbidities, area social deprivation, and access to treatment might contribute to inequities in mortality. We hypothesized that we would observe the same mortality disparity as found in previous studies of AI/AN and that this disparity would be partially explained by comorbid medical conditions, area social deprivation, and access to treatment.
Methods
Study design
The COVID EHR Cohort at the University of Wisconsin (CEC-UW;30) is a retrospective cohort study supported by the National Cancer Institute (ClinicalTrials.gov NCT04506528) that included 21 health systems from across the US (see Figure S1 in Supplemental Materials). Data extractions were performed using customized extraction code altered to accommodate unique health system specific EHR features. Each data extraction captured data on new patients meeting inclusion criteria and follow-up data on existing cohort members. Participating health systems provided selected data elements from the EHR of all COVID-19 patients encountered during the study period (February 1, 2020 to January 31, 2022). Data were transferred to the CEC-UW Coordinating Center in Madison, Wisconsin, where they were harmonized and merged. Harmonization, merging, and data analysis occurred September 30, 2021 through July 3, 2023. This study follows the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) reporting guidelines31.
Ethics statement
The CEC-UW study was initially approved in May 2020 by the University of Wisconsin-Madison Health Sciences Minimal Risk Institutional Review Board (MR-IRB) with approval for the collection of de-identified EHR data from the 21 health systems. The MR-IRB also determined that the study met criteria for a human subjects research exemption and qualified for a waiver of informed consent under the Federal Common Rule. All participating health systems provided written notice of either their own institution’s IRB approval or determination of exemption status before sharing EHR data. In February 2021, the MR-IRB approved a change of protocol for a Limited Data Set, allowing the collection of additional information (e.g., death dates, five-digit zip codes) but excluding direct patient identifiers. Each patient in the data set from each health system was assigned an enduring cryptographically processed Patient ID based on the SHA256 algorithm, which yielded a 64-character unique and private hash-based message authentication code (HMAC).
Analysis samples
The full CEC-UW inpatient cohort included 145,944 adult patients hospitalized with COVID-19 who had prior contact with the health system and who completed their hospitalization from February 1, 2020 through January 31, 2022 at a participating health system30.
This study focused on three analysis samples. One sample comprised all CEC-UW inpatients who identified as AI/AN (N = 546) [78% were non-Hispanic]) or NHW (N = 78,128). The other two were matched samples comprised of all CEC-UW inpatients who identified as AI/AN (N = 546) along with five matched NHW inpatients (N = 2645) selected using the SAS GMATCH macro32. (More information is provided below under “Statistical Analysis”).
Primary outcome
The primary outcome was all-cause in-hospital mortality during the index COVID-19 hospitalization. Because cause of death was not extracted from the EHRs of the participating health systems, we could not definitively attribute it to COVID-19 as some patients could have died during their hospitalization from other causes. Also reported were rates of three severity indicators (admission to the intensive care unit, intubation for ventilator use, and days in the hospital [among patients who did not die prior to their discharge]). The mean age at death was reported among those who died.
Patient characteristics
The following patient characteristics were extracted from the electronic heath records: sex, age, cigarette smoking status, co-occurring medical conditions (obesity, Type 2 diabetes, kidney, liver, and heart disease, cancer, alcohol use disorder, and drug use disorder [see Supplemental Text S1 for the ICD-10 codes corresponding to the medical conditions]), insurance type, receipt of antiviral medication during the hospitalization (see Supplemental Text S2 for a list of the antiviral medications prescribed), and COVID-19 vaccination status (see Supplemental Text S3). The date of the index hospitalization was extracted to incorporate a variable indicating the year in which the hospitalization occurred. (Note that although it is included as a race category in the EHR, we do not refer to AI/AN as a race based on the recommendation of our Indigenous consultant. The rationale is that AI/ANs are represented by 574 federally recognized tribes that are sovereign political governments and are not a single unified group.)
Additional variables were linked to the EHR data based on the patient’s home ZIP code: region (Northeast, Midwest, South, and West), urbanicity33, area social deprivation34, and distance to hospital. Urbanicity was based on Rural–Urban Commuting Area (RUCA) codes33; codes 1–6 were classified as “urban”, and codes 7–10 were classified as “rural.” Area social deprivation was based on seven census-based indicators of deprivation (see Table S1 in the Supplemental Materials) that were combined via factor analysis34,35 (see Supplemental Table S2). The area social deprivation index was calculated at the ZIP code tabulation area level; results are presented for each quintile based on the full sample. The ZIP code of the facility in which treatment was received was used to compute the distance to treatment using a SAS function that calculated the geodetic distance in miles between the centroids of two ZIP code locations, in this case the home and treatment facility ZIP codes. Because distance to treatment was very skewed and kurtotic, it was dichotomized at 60 miles, a distance that represents a significant barrier to getting timely treatment (that is, about 60 min travel time) while still including adequate numbers in each group for analyses.
These variables were included as covariates in adjusted multivariable models (covariate categorizations are given in Table 2).
Statistical analysis
Matching was conducted to create samples that were aligned on key variables while allowing for variability in important predictors. Two different matched samples were created. One set of matched samples (Match 1) was matched on age, region, month/year, and sex. The other set of matched samples (Match 2) was matched on region, month/year, and sex (not age). Because age was confounded with AI/AN versus NHW status in the full sample (see Table 2), NHW patients were matched to the AI/AN patients based on age (± 5 years) in Match 1. AI/AN and NHW patients were not matched on age when the outcome was age at death because age at admission and age at death for those who died during their hospitalization were nearly perfectly correlated (Match 2). Matching on region of residence and time of hospitalization minimized potential differences in availability of SARS CoV-2 vaccinations or exposure to different COVID-19 strains for AI/AN and NHW (Match 1 and 2). Both samples were matched based on sex (Match 1 and 2) to maintain the nearly equal proportions of men and women.
Analyses were conducted in the matched samples to identify factors that may contribute to disparities in mortality. First, we examined whether there were differences between AI/AN and NHW inpatients for each of the patient characteristics to identify potential explanations for the inequities in mortality. Second, we examined the extent to which each of these patient characteristics predicted mortality individually and after controlling for all the other characteristics in adjusted models. Of particular interest was the extent to which the association between AI/AN versus NHW status and mortality were reduced after patient characteristics, such as smoking and medical comorbidities, access to treatment (as indicated by distance to treatment, insurance type, residence in a rural or urban region, and receipt of antiviral medication while hospitalized), and area social deprivation were considered. The impact of controlling for these characteristics was quantified by examining the percentage attenuation in the effect sizes obtained in unadjusted and adjusted models: 100 × (Bunadjusted model − Badjusted model)/(Bunadjusted model) cf.36,37. Third, in the fully adjusted models, we examined whether any of the patient characteristics differentially predicted mortality in AI/AN and NHW inpatients. The purpose was to determine whether a potential risk factor exerted a greater or lesser effect on mortality among AI/AN compared to NHW individuals. Analyses were conducted in SPSS38. Multilevel generalized linear models with a binomial distribution and a logit link were fit. Multilevel analysis was used to account for the clustering of patients within the 21 health systems.
Bootstrapped confidence intervals around percentages and means and t-tests were estimated in SPSS. The Benjamini–Hochberg procedure39 was applied to control the false discovery rate in the multivariable analyses and results of such correction are shown in relevant tables.
Results
Differences in mortality and other outcomes between AI/AN and NHW inpatients
AI/AN patients were more likely to die during their hospitalization than were NHW patients (Table 1 and Fig. 1). AI/AN patients died at a significantly younger age than NHW patients in the full and matched samples, which was not surprising given that AI/ANs were admitted at a younger age. AI/ANs were also more likely to require intubation, be admitted to the intensive care unit, and spend more days in the hospital than NHWs; Mann–Whitney tests indicated that AI/ANs and NHWs did not significantly differ in the number of days hospitalized in the full sample but did in the matched sample (Table 1).
Differences in patient characteristics between AI/AN and NHW inpatients
Differences in patient characteristics between AI/ANs and the full and matched NHW samples are presented in Table 2. As mentioned, AI/AN patients were significantly younger than the full NHW sample but, as expected, were not younger than the matched NHW sample. Some of the significant differences observed in the full sample (smoking status, obesity, heart disease, drug use disorder, insurance type, and vaccination status) were no longer significant in the matched samples, and one new difference (chronic renal failure) was revealed by the matching process.
In the matched samples, AI/AN patients who were hospitalized with COVID-19 were significantly more likely than NHW patients to be diagnosed with the comorbid conditions of Type 2 diabetes, chronic renal failure, liver disease, and alcohol use disorder. They also had significantly higher comorbidity burden as indicated by a count of the number of comorbid disorders (median 2.0) than NHW patients (median 1.0); Mann–Whitney test (U = 842,344.00, z = 6.47, p < 0.001).
Although the AI/AN patients were predominantly urban dwelling (86.4%), they were more likely than NHW patients to live in a rural area and to live 60 or more miles from the facility in which they received treatment. AI/AN patients were significantly more likely than NHW patients to live in socioeconomically deprived areas. For example, they were more likely to live in areas with a higher percentage of residents living below the federal poverty level, with less than a high school education, who did not own a car, who lived in a crowded housing unit and in households headed by a single parent (see Supplemental Table S1). In summary, these results suggest that comorbid medical conditions, access to treatment, and area social deprivation are potential contributors to the inequity in mortality between AI/AN and NHW individuals with COVID-19. AI/AN and NHW patients did not significantly differ in smoking status, some medical comorbidities (obesity, heart disease, cancer, and drug use disorder), type of insurance coverage, receipt of antiviral medication, year admitted to the hospital, and vaccination status. These results suggest that these are unlikely to be contributors to the inequity in mortality between AI/AN and NHW individuals with COVID-19.
Associations of AI/AN versus NHW status and patient characteristics with mortality
As presented in Table 3, several comorbidities were associated with in-hospital mortality (after covariate adjustment and control for false discovery rate): chronic renal failure, liver disease, and heart disease. Several indicators of health care access were associated with in-hospital mortality in unadjusted analyses: being a Medicare recipient, receipt of antiviral medication, living a greater distance from the treatment facility, and living in a rural area; being in the top quintile of area social deprivation was also associated with in-hospital mortality, but after covariate adjustment and control for false discovery rate, being a Medicare recipient and receipt of antiviral medication were the only ones that remained significant predictors of mortality.
The unadjusted odds ratio of the association between AI/AN versus NHW status and mortality was 1.98 (95% CI 1.48, 2.65). After adjusting for comorbidities, access to treatment, area social deprivation and vaccination status, the odds ratio was reduced to 1.59 (95% CI 1.15, 2.22), which represents a 32% reduction in the association between AI/AN versus NHW status and mortality.
Analyses were conducted to identify the source of this diminution in the association between AI/AN versus NHW status and mortality. In a model that adjusted for only comorbidities (excluding those that were inversely associated with mortality in the adjusted model), the odds ratio was reduced to 1.82 (95% CI 1.35, 2.46), which represents a 12% reduction in the association between AI/AN versus NHW status and mortality. This effect was further probed by individually examining each of the six comorbid disorders. The single disorder that accounted for most of the mortality disparity was any liver disease (odds ratio reduced to 1.87 [95% CI 1.40, 2.51], representing a 9% reduction; the others ranged from 1.92 to 1.99). When specific liver diseases were included in the adjusted model presented in Table 3, the single best predictor of mortality with COVID-19 was hepatic failure, not elsewhere classified (OR 8.05, 95% CI 4.22, 15.36). Given that liver disease is often secondary to hepatitis B and C and HIV infections40, we conducted post hoc analyses to determine that the relation between liver disease and mortality could not be explained by comorbid hepatitis or HIV (see Supplemental Text S4).
In a model that adjusted for only access to treatment (distance to treatment, insurance type, residence in a rural or urban region, and receipt of antiviral medication while hospitalized), the odds ratio was reduced to 1.82 (95% CI 1.34, 2.47), also representing a 12% reduction. In a model that adjusted for only area social deprivation, the odds ratio was reduced to 1.79 (95% CI 1.32, 2.43), representing a 15% reduction. These results suggest that comorbidities (especially comorbid liver disease), access to treatment and area social deprivation all contribute to the disparity in COVID-19 related mortality between AI/AN and NHW individuals. The significant residual relation between AI/AN versus NHW status and mortality after accounting for these candidate explanatory variables indicate that they do not fully account for this disparity.
Moderation of associations between AI/AN versus NHW status and mortality
We examined whether any of the patient characteristics differentially predicted mortality in AI/AN and NHW inpatients in fully adjusted models by including an interaction term between AIAN versus NHW status and each of the patient characteristics. After control for false discovery rate, there was no evidence that any of the potential risk factors exerted a greater or lesser effect on mortality among AI/AN compared to NHW individuals (see Supplemental Table S3). Before control for false discovery rate, however, there was a significant interaction between AIAN versus NHW status and alcohol use disorder (see Supplemental Text S5). The magnitude of the association between AI/AN versus NHW status and mortality did not change from 2020 to 2022.
Discussion
This study compared the outcomes of 546 AI/AN to 78,128 NHW individuals from 21 US health systems who were hospitalized with COVID-19 from February 1, 2020 to January 31, 2022. To better isolate differences, the AI/AN sample was also compared to a sample of 2645 NHW individuals matched on age, sex, region, and month/year of hospital admission. As expected, based on previous epidemics and recent studies4,5,6,7, the odds of in-hospital mortality were doubled in the AI/AN compared to the matched sample of NHW individuals with COVID-19. Novel to this study was the direct examination of the extent to which comorbidities, area social deprivation, and access to treatment could contribute to the AI/AN disparity in mortality.
Comorbidity and mortality differences between AI/AN and NHW
AI/AN inpatients were compared to NHW inpatients on eight diseases known to be more common among AI/AN12,25 and that have been linked to severe COVID-19 outcomes13,14. In this inpatient sample, AI/AN patients were more likely to have been diagnosed with three of these (diabetes, chronic renal failure, and liver disease); liver disease accounted for the largest portion of the comorbidity-related mortality disparity between AI/AN and NHW patients.
Area social deprivation and mortality differences between AI/AN and NHW
AI/AN inpatients were more than twice as likely to reside in the most socially deprived areas relative to NHW inpatients. Consistent with prior research19,41 higher area social deprivation increased the risk of COVID-19 inpatients’ dying in the hospital. Furthermore, adjusting for area social deprivation attenuated the association between AI/AN versus NHW status and hospital mortality. The present findings extend the limited number of studies that have linked area social deprivation with patient-level COVID-19 clinical outcomes.
Access to treatment and mortality differences between AI/AN and NHW
Higher rates of intubation, ICU admission, and days hospitalized suggest that AI/AN inpatients may have presented to inpatient care later in the COVID-19 disease course compared to NHW inpatients. In this study, 23% of AI/AN inpatients lived 60 or more miles from their treatment facility, compared to only 10% of the matched sample of NHW. It is likely that this greater distance to treatment may account for the greater severity of COVID-19 illness and ultimately higher rates of death in AI/AN than NHW patients. The association between AI/AN versus NHW status and mortality was reduced after adjusting for distance to treatment. More research is needed to identify other obstacles to timely treatment that AI/AN individuals are more likely to encounter than their NHW counterparts42. In particular, transportation barriers43 and lack of internet access17,18 may also impede timely treatment.
Differences between AI/AN and NHW in age at death
Previous studies have documented younger ages at death for AI/AN than NHW individuals with COVID-199,44. In the full sample of the present study, AI/AN patients were younger (d = 0.74) and had younger in-hospital ages at death (d = 0.74) than NHW patients. Of those who died during their hospitalization, 60% of the AI/ANs compared to only 34% of the NHWs were less than 70 years old. The contributors to mortality disparities previously discussed also likely contributed to disparities in age at death. AI/ANs were more likely than NHWs to be living in crowded housing, which may have promoted greater viral load. Living further from treatment may have led to greater delays in getting or seeking treatment and more advanced illness at hospitalization for AI/ANs compared to NHWs.
Disparities in age at death might also be explained by the concept of “weathering” that posits that the cumulative impact of repeated experience with social or economic adversity and political marginalization may lead to physiological deterioration45. This cumulative wear and tear has been termed “allostatic load”45. An empirical demonstration of weathering was conducted in a community-based sample of Black and White individuals showing that the allostatic load of a 40-year-old Black person was equal to that of a 50-year-old White person46. To our knowledge, studies of the causes of disparities in age at death have not been conducted among AI/AN individuals but should be a top priority for future research. The profound loss of lifespan among AI/AN has widespread effects given the value that Indigenous communities place on their elders28,47.
Other studies quantitatively explaining COVID-19 related mortality disparities
Several studies have systematically documented COVID-19 related mortality disparities among minority populations (e.g.,48,49). To our knowledge, only two studies, both conducted in the UK, have attempted to quantitatively explain the potential causes of these mortality disparities36,50. In one UK study36, comorbidities explained 10% of the association between Black/White status and mortality and 39% of association between South Asian/White status and mortality; social factors (educational attainment, occupational attainment, household size and area deprivation) explained 28% of the association between Black/White status and mortality and 4% of association between South Asian/White status and mortality. In the other UK study50, comorbidities and social factors (household size and area deprivation), but not lifestyle factors (smoking, body mass index), explained about 40% of the association between Black/White status and mortality and between South Asian/White status and mortality. (Access to treatment was not included in either of the UK studies.) In the present study, comorbidities explained 12% of the association between AI/AN status and mortality and area social deprivation explained 15% of the association between AI/AN status and mortality. Taken together, results from the UK studies and the present study suggest that comorbidities and area social deprivation account for similarly small fractions of the COVID-19 related mortality disparity among Black, South Asian, and AI/AN individuals. Much of the disparity was left unexplained for all three ethnic groups.
Limitations
First, the sample included only hospitalized AI/AN and NHW patients during their first hospitalization for COVID-19, so it does not reflect the course of COVID-19 and mortality differences in the broader AI/AN population. Second, outcomes occurring post-discharge or outcomes occurring at nonparticipating health systems were not captured. Third, results across time could not be linked with type of COVID-19 variant. The analyses were conducted over the first 2 years of the pandemic, suggesting that the data obtained were contemporaneous with high prevalence of alpha, delta, and early omicron variants51. Fourth, risk factors for COVID-19 related mortality were considered in isolation when they were more likely to act in concert or be stages in a causal chain. For example, area social deprivation and distance to treatment may be barriers to adequate prevention and intervention for comorbid liver disease52. Future research should model the process by which risk factors combine to influence COVID-19 related mortality.
Fifth, the CEC-UW inpatient cohort is not a representative sample. There was selection of the participants based on having a diagnosis of COVID-19 and for being hospitalized. It is well known that studies based on EHR data are plagued by collider bias53,54,55. In this case, COVID-19 may be a collider associated with both AI/AN versus NHW status and mortality and could have induced distorted or spurious findings56. It is reassuring that the association between AI/AN versus NHW status and mortality observed in this study was similar to results obtained from other sources, such as state-level surveillance systems5,7,8. However, the positive association of receipt of antiviral medication with mortality may have been due to collider bias (but also possibly due to sicker patients being more likely to be prescribed medication). Although the association with mortality may have been distorted, it was an important observation that AI/AN were not less likely than NHW to receive pharmacologic treatment for their COVID-19.
Conclusions
Notwithstanding limitations, this study is an important contribution to the literature because it represents the first attempt to explain COVID-19 related mortality disparities among AI/AN. Comorbidities, area social deprivation, and access to treatment were all important contributors to the mortality disparity between AI/AN and NHW inpatients with COVID-19. Nonetheless, the significant residual relation between AI/AN versus NHW status and mortality after accounting for the candidate explanatory variables of comorbidity burden, neighborhood socioeconomic deprivation and reduced access to health care indicate that there are other important unmeasured factors that contribute to this inequity. This likely includes living conditions, such as multigenerational and crowded housing57, being a frontline worker57, and having inadequate access to transportation43 and to the internet17,18. Accounting for the unexplained causes of disparities among AI/AN will be an important direction for future research.
Health disparities among AI/AN are not a new problem but reflect “legacies of failing to address historical and ongoing inequities” (58, p. 2739). Results of the present study likely extend beyond the current COVID-19 pandemic and may apply to many other past, current, and future health disparities experienced by AI/AN. The availability of quality data on the disparate impacts of health threats such as COVID-19 on AI/ANs is essential in the effort to reduce disparities and enhance health equity.
Data availability
The existing Data Transfer and Use Agreements negotiated with each of the participating health systems preclude the University of Wisconsin from sharing CEC-UW data with any entity. Information Management Services, Inc. (IMS), under contract with the National Cancer Institute, is responsible for housing the CEC-UW dataset. Investigators desiring access to CEC-UW data can apply to IMS (https://www.imsweb.com/).
References
-
Jones, D. S. The persistence of American Indian health disparities. Am J Public Health 96, 2122–2134 (2006).
-
Connolly, M., Jacobs, B. & Notzon, F. C. COVID-19 among American Indians and Alaska Natives in the United States: An early look. Stat. J. IAOS 37, 25–36 (2021).
-
Castrodale, L. et al. Deaths related to 2009 pandemic influenza A (H1N1) among American Indians—12 states, 2009. MMWR Morb. Mortal. Wkly. Rep. 58, 1341–1344 (2009).
-
Groom, A. V. et al. Pneumonia and influenza mortality among American Indian and Alaska Native people, 1990–2009. Am. J. Public Health 104, S460–S469 (2014).
-
Arrazola, J. et al. COVID-19 mortality among American Indian and Alaska native persons—14 states, January–June 2020. MMWR Morb. Mortal. Wkly. Rep. 69, 1853 (2020).
-
Feldman, J. M. & Bassett, M. T. Variation in COVID-19 mortality in the US by race and ethnicity and educational attainment. JAMA Netw. Open 4, e2135967 (2021).
-
Hatcher, S. M. et al. COVID-19 among American Indian and Alaska native persons—23 states, January 31–July 3, 2020. MMWR Morb. Mortal. Wkly. Rep. 69, 1166–1169 (2020).
-
Ward, L. A., Black, K. P., Britton, C. L., Tompkins, M. L. & Provost, E. M. COVID-19 cases, hospitalizations, and deaths among American Indian or Alaska Native persons—Alaska, 2020–2021. MMWR Morb. Mortal. Wkly. Rep. 71, 730 (2022).
-
Arias E, Tejada-Vera, B, Kochanek KD, Ahmad FB. Provisional life expectancy estimates for 2021. Vital Stat. Surveill. Rep. 23 (2022).
-
Ahmad, F. B., Cisewski, J. A., Xu, J. & Anderson, R. N. Provisional mortality data—United States, 2022. MMWR Morb. Mortal. Wkly. Rep. 72, 488–492 (2023).
-
Urban Indian Health Commission. Invisible Tribes: Urban Indians and Their Health in a Changing World. Seattle: Urban Indian Health Commission, 2007. https://www2.census.gov/cac/nac/meetings/2015-10-13/invisible-tribes.pdf. Accessed 17 Nov 2022.
-
Webb Hooper, M., Nápoles, A. M. & Pérez-Stable, E. J. COVID-19 and racial/ethnic disparities. JAMA 323, 2466–2467 (2020).
-
Bennett, T. D. et al. Clinical characterization and prediction of clinical severity of SARS-CoV-2 infection among US adults using data from the US National COVID Cohort Collaborative. JAMA Netw. Open 4, e2116901 (2021).
-
Harrison, S. L., Fazio-Eynullayeva, E., Lane, D. A., Underhill, P. & Lip, G. Y. Comorbidities associated with mortality in 31,461 adults with COVID-19 in the United States: A federated electronic medical record analysis. PLoS Med. 17, e1003321 (2020).
-
Musshafen, L. A. et al. In-hospital mortality disparities among American Indian and Alaska Native, Black, and White patients with COVID-19. JAMA Netw. Open 5, e224822–e224822 (2022).
-
Qeadan, F. et al. Racial disparities in COVID-19 outcomes exist despite comparable Elixhauser comorbidity indices between Blacks, Hispanics, Native Americans, and Whites. Sci. Rep. 11, 8738 (2021).
-
Ehrenpreis, J. E. & Ehrenpreis, E. D. A historical perspective of healthcare disparity and infectious disease in the Native American population. Am. J. Med. Sci. 363, 288–294 (2022).
-
Sequist, T. D. Improving the health of the American Indian and Alaska Native population. JAMA 325, 1035–1036 (2021).
-
Wong, M. S. et al. The role of community-level factors on disparities in COVID-19 infection among American Indian/Alaska Native veterans. J. Racial Ethn. Health Disparities 9, 1861–1872 (2022).
-
Madhav, K. C., Oral, E., Straif-Bourgeois, S., Rung, A. L. & Peters, E. S. The effect of area deprivation on COVID-19 risk in Louisiana. PLoS ONE 15, e0243028 (2020).
-
Zhang, Y. et al. Socioeconomic variation in characteristics, outcomes, and healthcare utilization of COVID-19 patients in New York City. PLoS ONE 6, e0255171 (2021).
-
Burki, T. COVID-19 among American Indians and Alaska Natives. Lancet 21, 325–326 (2021).
-
U.S. Department of Health and Human Services Office of Minority Health. Profile: American Indian/Alaska Native. https://www.minorityhealth.hhs.gov/omh/browse.aspx?lvl=3&lvlid=62. Accessed 10 Oct 2022.
-
Adakai, M. et al. Health disparities among American Indians/Alaska natives—Arizona, 2017. MMWR Morb. Mortal. Wkly. Rep. 67, 1314–1318 (2018).
-
Urban Indian Health Commission. Invisible Tribes: Urban Indians and Their Health in a Changing World. Seattle: Urban Indian Health Commission (2007).
-
Urban Indian Health Institute. Data Genocide of American Indians and Alaska Natives in COVID-19 Data. Seattle, WA: Urban Indian Health Institute (2021). https://www.uihi.org/projects/data-genocide-of-american-indians-and-alaska-natives-in-covid-19-data/. Accessed 11 Nov 2022.
-
Yellow Horse, A. J. & Huyser, K. R. Indigenous data sovereignty and COVID-19 data issues for American Indian and Alaska Native Tribes and populations. J. Popul. Res. 39, 527–531 (2021).
-
Mays, V. M., Echo-Hawk, A., Cochran, S. D. & Akee, R. Data equity in American Indian/Alaska Native populations: Respecting sovereign nations’ right to meaningful and usable COVID-19 data. Am. J. Public Health 112, 1416–1420 (2022).
-
Carroll, S. R., Rodriguez-Lonebear, D., Akee, R., Lucchesi, A. & Richards, J. R. Indigenous Data in the COVID-10 Pandemic: Straddling Erasure, Terrorism, and Sovereignty. Social Science Research Council, Items: Insights from the Social Sciences (2020). https://items.ssrc.org/covid-19-and-the-social-sciences/disaster-studies/indigenous-data-in-the-covid-19-pandemic-straddling-erasure-terrorism-and-sovereignty/. Accessed 11 Nov 2022.
-
Fiore, M. C. et al. The first 20 months of the COVID-19 pandemic: Mortality, intubation and ICU rates among 104,590 patients hospitalized at 21 United States health systems. PLoS ONE 17, e027451 (2022).
-
von Elm, E. et al. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: Guidelines for reporting observational studies. Lancet 370, 1453–1457 (2007).
-
Bergstralh, E. J., Kosanke, J. L. & Jacobsen, S. J. Software for optimal matching in observational studies. Epidemiology 7, 331–332 (1996).
-
U.S. Department of Agriculture Economic Research Service. (2020). Rural-Urban Commuting Area Codes (RUCA). https://www.ers.usda.gov/data-products/rural-urban-commuting-area-codes/documentation/. Accessed 9 Oct 2022.
-
Robert Graham Center. Social Deprivation Index (SDI). https://www.graham-center.org/rgc/maps-data-tools/sdi/social-deprivation-index.html. Accessed 9 Oct 2022.
-
Butler, D. C., Petterson, S., Phillips, R. L. & Bazemore, A. W. Measures of social deprivation that predict health care access and need within a rational area of primary care service delivery. Health Serv. Res. 48, 539–559 (2013).
-
Batty, G. D., Gaye, B., Gale, C. R., Hamer, M. & Lassale, C. Explaining ethnic disparities in COVID-19 mortality: Population-based, prospective cohort study. Am. J. Epidemiol. 191, 275–281 (2022).
-
Stringhini, S. et al. Association of socioeconomic position with health behaviors and mortality. JAMA 303, 1159–1166 (2010).
-
IBM Corp. IBM SPSS Statistics for Windows. Version 27.0. Armonk, NY: IBM Corp; 2020.
-
Benjamini, Y. & Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B 57, 289–300 (1995).
-
Joshi, D., O’Grady, J., Dieterich, D., Gazzard, B. & Agarwal, K. Increasing burden of liver disease in patients with HIV infection. Lancet 377, 1198–1209 (2011).
-
Hu, J. et al. Race, ethnicity, neighborhood characteristics, and in-hospital coronavirus disease-2019 mortality. Med. Care 59, 888–892 (2021).
-
Agency for Healthcare Research and Quality, National Healthcare Quality and Disparities Report. Rockville, MD: October 2022. AHRQ Pub. No. 22.
-
Syed, S. T., Gerber, B. S. & Sharp, L. K. Traveling towards disease: Transportation barriers to health care access. J. Community Health 38, 976–993 (2013).
-
Goldman, N. & Andrasfay, T. Life expectancy loss among Native Americans during the COVID-19 pandemic. Demogr. Res. 47, 233–246 (2022).
-
Geronimus, A. T., Hicken, M., Keene, D. & Bound, J. “Weathering” and age patterns of allostatic load scores among blacks and whites in the United States. Am. J. Public Health 96, 826–833 (2006).
-
Seeman, T. E., Singer, B. H., Rowe, J. W., Horwitz, R. I. & McEwen, B. S. Price of adaptation—Allostatic load and its health consequences: MacArthur studies of successful aging. Arch. Intern. Med. 157, 2259–2268 (1997).
-
Viscogliosi, C. et al. Importance of Indigenous elders’ contributions to individual and community wellness: Results from a scoping review on social participation and intergenerational solidarity. Can. J. Public Health 111, 667–681 (2020).
-
Irizar, P. et al. Ethnic inequalities in COVID-19 infection, hospitalisation, intensive care admission, and death: A global systematic review and meta-analysis of over 200 million study participants. EClinicalMedicine 57, 101877 (2023).
-
Raharja, A., Tamara, A. & Kok, L. T. Association between ethnicity and severe COVID-19 disease: A systematic review and meta-analysis. J. Racial Ethn. Health Disparities 8, 1563–1572 (2021).
-
Zaccardi, F. et al. Ethnic disparities in COVID-19 outcomes: A multinational cohort study of 20 million individuals from England and Canada. BMC Public Health 23, 399 (2023).
-
Iuliano, A. D. et al. Trends in disease severity and health care utilization during the early omicron variant period compared with previous SARS-CoV-2 high transmission periods—United States, December 2020–January 2022. MMWR Morb. Mortal. Wkly. Rep. 71, 146–152 (2022).
-
Kardashian, A., Wilder, J., Terrault, N. A. & Price, J. C. Addressing social determinants of liver disease during the COVID-19 pandemic and beyond: A call to action. Hepatology 73, 811–820 (2021).
-
Griffith, G. J. et al. Collider bias undermines our understanding of COVID-19 disease risk and severity. Nat. Commun. 11, 5749 (2020).
-
Weiskopf, N. G., Dorr, D. A., Jackson, C., Lehmann, H. P. & Thompson, C. A. Healthcare utilization is a collider: An introduction to collider bias in EHR data reuse. JAMIA 30, 971–977 (2023).
-
Bots, S. H., Groenwold, R. H. & Dekkers, O. M. Using electronic health record data for clinical research: A quick guide. Eur. J. Endocrinol. 186, E1–E6 (2022).
-
Holmberg, M. J. & Andersen, L. W. Collider bias. JAMA 327, 1282–1283 (2022).
-
Leggat-Barr, K., Uchikoshi, F. & Goldman, N. COVID-19 risk factors and mortality among Native Americans. Demogr. Res. 45, 1185–1218 (2021).
-
Power, T. et al. COVID-19 and indigenous peoples: An imperative for action. J. Clin. Nurs. 29, 2737–2741 (2020).
Acknowledgements
Sources of funding and support: The CEC-UW data collection was funded by a contract from the National Cancer Institute (CRDF Award #66590 to B. Rolland). Funding for this project was provided by the UW School of Medicine and Public Health from the Wisconsin Partnership Program (to M. Fiore).
Disclaimer
The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the Nation Cancer Institute or the Wisconsin Partnership Program. The University of Wisconsin–Madison acknowledges that it occupies ancestral Ho-Chunk land.
Author information
Authors and Affiliations
Contributions
All authors contributed to the work presented in this paper. M.C.F led the preservation of the original data on which the paper is based, with the assistance of K.L.C. S.S.S assembled and harmonized data used in the manuscript with assistance from W.S.S, K.L.C, and J.A.K. W.S.S. conducted the statistical analyses with assistance from J.A.K, S.S.S, and T.M.P. W.S.S., K.L.C, J.A.K., S.S.S., T.M.P., A.L.J., D.E.M., P.N.H. and M.C.F. wrote and provided feedback on the paper. All authors discussed the results and implications and commented on the manuscript at all stages.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Slutske, W.S., Conner, K.L., Kirsch, J.A. et al. Explaining COVID-19 related mortality disparities in American Indians and Alaska Natives.
Sci Rep 13, 20974 (2023). https://doi.org/10.1038/s41598-023-48260-9
-
Received:
-
Accepted:
-
Published:
-
DOI: https://doi.org/10.1038/s41598-023-48260-9
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.