Disparities in outcomes of COVID-19 hospitalizations in native American individuals

Objectives: This study aimed to investigate COVID-19-related disparities in clinical presentation and patient outcomes in hospitalized Native American individuals.

Methods: The study was performed within 30 hospitals of the Banner Health system in the Southwest United States and included 8,083 adult patients who tested positive for SARS-CoV-2 infection and were hospitalized between 1 March 2020 and 4 September 2020. Bivariate and multivariate analyses were used to assess racial and ethnic differences in clinical presentation and patient outcomes.

Results: COVID-19-related hospitalizations in Native American individuals were over-represented compared with non-Hispanic white individuals. Native American individuals had fewer symptoms at admission; greater prevalence of chronic lung disease in the older adult; two times greater risk for ICU admission despite being younger; and 20 times more rapid clinical deterioration warranting ICU admission. Compared with non-Hispanic white individuals, Native American individuals had a greater prevalence of sepsis, were more likely to require invasive mechanical ventilation, had a longer length of stay, and had higher in-hospital mortality.

Conclusion: Native American individuals manifested greater case-fatality rates following hospitalization than other races/ethnicities. Atypical symptom presentation of COVID-19 included a greater prevalence of chronic lung disease and a more rapid clinical deterioration, which may be responsible for the observed higher hospital mortality, thereby underscoring the role of pulmonologists in addressing such disparities.

Introduction

The coronavirus 2019 (COVID-19) pandemic has been a major threat to public health. As of August 2022, the number of confirmed SARS-CoV-2 infections in the United States has reached 93.7 million with a death toll surpassing 1 million individuals (1). COVID-19 has unevenly impacted different populations (2, 3). Vulnerable groups such as the older adults, disabled, and racial minorities have been more severely affected (3, 4). Studies have shown that African Americans were disproportionally affected by COVID-19 with greater rates of hospitalization and death than non-Hispanic white individuals (3, 5). Such health disparities in the African American community were associated with social determinants of health, including lower income, lack of health insurance, poor access to care, crowded housing, lack of independent transportation, and greater exposure to COVID-19 (6). Native Americans and Hispanic individuals suffer a greater burden of chronic medical illnesses, such as diabetes, obesity, and heart disease, which, in turn, may confer a greater risk for COVID-19 (7, 8). However, few studies have examined the COVID-19 hospitalization and mortality rates in Native American populations (9–11). A recent systematic review and meta-analysis stated that there was insufficient evidence in Native Americans to draw strong conclusions regarding disparities in COVID-19 and that more research was needed (12). Moreover, they concluded that the increased susceptibility to COVID-19 alone does not seem to explain the observed disparities in other minority populations (African American individuals and Hispanic individuals) (12). Perhaps barriers to healthcare access and delayed medical attention may be responsible for the observed disparities. The unique challenges faced by Native American populations are well documented and could explain the increased risk of COVID-19 morbidity and mortality (13–15). In addition to having a higher burden of risk factors for severe COVID-19 disease (15), poor access to tertiary care centers could lead to a delay in presentation, hence poorer outcomes (13). We explore the preventable differences in hospitalization including the burden of intensive care unit (ICU) admissions across such race/ethnicities and explore excess hospital mortality, which are defined as constituting health disparities by the National Institutes of Health for Minority Health and Disparities (NIMHD) and Centers for Disease Control (CDC) (16, 17).

Moreover, although African American and Hispanic populations experience disproportionately greater rates of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection and COVID-19-related mortality, case-fatality rates were similar in non-Hispanic white individuals (12). But inadequate adjustment for the acuity of illness or inter-hospital transfers—which are major determinants of in-hospital mortality—may have prevented the detection of differences in case-fatality rates as the studies in Native Americans were derived primarily from hospital reporting of deaths to the CDCl (9).

The objective of this study was to investigate COVID-19-related disparities in clinical presentation and patient outcomes in the Native American population. Additionally, the study aimed to assess the case-fatality rates after fully adjusting for the acuity of illness and inter-hospital transfers in minority populations with COVID-19 requiring hospitalization.

Methods

Study design, setting, and population

This was a retrospective observational cohort study conducted within the Banner Health system (30 acute care hospitals) in the Southwest United States. Banner Health is a large nationally recognized health system operating in six states in the United States: Arizona, California, Nebraska, Nevada, Wyoming, and Colorado. Most of the hospitals are located in Arizona, including three academic medical centers: Banner—University Medical Center Tucson, Banner—University Medical Center Phoenix, and Banner—University Medical Center South. Banner Health provides healthcare for a large segment of the Native American population in Arizona and neighboring states in the mountain west of the United States.

The study included adult (≥18) patients who tested positive for infection with SARS-CoV-2 using a qualitative polymerase-chain-reaction assay and were hospitalized between 1 March 2020 and 4 September 2020. We obtained data prior to the availability of the COVID-19 vaccine in order to avoid the effects of vaccine hesitancy, vaccination rates, monoclonal antibodies, and routine administration of steroids and remdesivir on observed disparities (18–20). The Institutional Review Board of Banner Health (IRB#483-20-0076) approved this study.

Data collection

Data were retrospectively collected from the electronic medical records (EMR) stored in a data warehouse (Banner Health Enterprise Data Warehouse). Data included demographics (age, sex, and race), health insurance, patient self-reported address, admission source, ICU utilization, timestamps of registration and ICU admission, timestamps of discharge from ICU and hospital, discharge disposition (alive or expired), and use of invasive mechanical ventilation. Self-reported ethnicity in the EMR system was used to identify the Hispanic population, and self-reported race information was used to identify other racial groups. Racial groups in this study include non-Hispanic white, Hispanic, Black or African American, Native American, Asian, Middle Eastern Indian, Native Hawaiian, or Other Pacific Islander. To ensure sufficient sample size in each comparator group, we collapsed cases into the following five race/ethnicity groups: non-Hispanic white, Hispanic, Black, Native American, and others. The median income for zip codes of patient residence was used as a surrogate for socioeconomic status. Insurance plans were categorized into five groups: Private, Medicaid, Medicare, Self-Pay, and others (Supplementary Table 1). Chronic conditions were collected based on previous ICD10 diagnosis codes related to any type of encounter in the last 5 years within the Banner delivery network (Supplementary Table 2). COVID-19 symptoms such as cough, fever, dyspnea, abdominal symptoms, and myalgia were extracted from an unstructured text field “chief complaint” from the patient's records at the beginning of visits.

Statistical analysis

We compared in-hospital outcomes of patients with COVID-19 by race/ethnicity groups. Outcome (dependent) variables included in-hospital mortality, invasive mechanical ventilation, ICU admission, and inter-hospital transfer vs. admission through ED, ICU LOS, and hospital LOS. Time (hours) from entering the hospital (“registration”) to ICU admission and ICU LOS were also deemed outcome variables in a subgroup of patients requiring ICU admission. Time from registration to ICU admission was examined to reflect how quickly a patient's condition deteriorated. Both hospital LOS and ICU LOS were investigated in live discharges. Covariates included various chronic medical conditions, insurance status, type of insurance, median income imputed from ZIP code of residence, age, sex, admission through ED or inter-hospital transfer, and symptom presentation (cough, fever, dyspnea, and abdominal symptoms). Bivariate analyses were used to characterize patients across the racial groups. Continuous variables were expressed using median and interquartile range (IQR). Categorical variables were described with frequencies and percentages. Univariate logistic regression was used to examine distributions of categorical variables among race/ethnicity groups, and ANOVA was used for examining continuous variables among race/ethnicity groups. The chi-square tests and Student's t-tests were used for post-hoc analyses. Two-sided exact P-values were reported, and P-values of ≤ 0.05 were considered statistically significant. P-values were adjusted by multiplicity with Bonferroni adjustment. Four strata by age group (18–39, 40–59, 60–79, and ≥80 years) were created to age-stratify distributions of chronic conditions, insurance plans, and outcome factors.

Multivariable logistic regression and linear regression were performed to examine factors associated with outcome variables in both hospitalized patients and ICU patients. Collinearity was determined among covariates using the Pearson correlation and variance inflation factor (VIF) methods (Supplementary Tables 3, 4 for inpatients, Supplementary Tables 5, 6 for ICU patients). The criteria (21) for selecting covariates of the model were as follows: (1) All covariates with P of < 0.25 in the bivariate analyses were excluded (22). (2) Pearson correlation coefficients <0.75 and VIF <5 for each covariate were included (23). The variable for the total number of chronic conditions was removed from the regression analysis based on the guidance for controlling confounders suggested by Lederer et al. (21). Regression assumptions such as linearity, homoscedasticity, and normality of error terms were checked before performing univariate and multivariate regressions. Median income values were missing in 12.3% of individuals due to a lack of residential address. An imputation with a sample median was used to replace the missing values. All analyses were performed using Python 3.6.2 (Python Software Foundation, Beaverton, OR).

Results

In total, there were 8,083 inpatients with COVID-19 (≥18 years), 7,515 (92.9%) were in Arizona. In such COVID-19 inpatients, when compared to the proportion of race/ethnicity populations in Arizona (2019 census) (24), the minority race/ethnicities—Native American (7.5%), Black (6.0%), and Hispanic individuals (40.6%)—were over-represented (all P < 0.0001), while non-Hispanic white individuals (41.7%) were under-represented (P < 0.0001, Supplementary Figure 1).

Characteristics of hospitalized patients by race using bivariate analysis

Age, sex, and socioeconomic status

Among 8,083 hospitalized patients with COVID-19, the median age of patients was 60 years (interquartile range [IQR] 46–73) and 50.4% were male patients. Native American, Black, and Hispanic patients were younger than the non-Hispanic white patients (Table 1). Of patients <60 years, 67.6% were Native American, 58.8% Hispanic, and 55.6% Black, compared with 34.8% of non-Hispanic white individuals (Supplementary Figure 2A). The median income of Native Americans was lower than non-Hispanic white individuals (Table 1; P < 0.0001). Native American patients were 3.5-fold more likely to be on Medicaid insurance than non-Hispanic white individuals (51.75 vs. 14.3%, Table 1; P < 0.0001). In age-stratified data, Native American patients aged 18–39 years had the highest proportion of Medicaid coverage among all age categories (Supplementary Figure 2B, P < 0.001). Hispanic (33.6%) and Black patients (28.5%) had higher rates of Medicaid coverage than non-Hispanic white individuals (14.2%) (Table 1; P < 0.0001).

TABLE 1

Table 1. Characteristics of hospitalized COVID-19 patients, overall and by race.

Comorbidities and clinical symptoms

Diabetes mellitus, obesity, and chronic liver disease were more prevalent in Native Americans than in non-Hispanic white individuals (Table 1; all P < 0.0001), and this was consistent across all age categories (all P < 0.05, Supplementary Figures 3A–C). The rate of chronic lung disease in Native Americans at age 60–79 years was 2-fold greater than in non-Hispanic white individuals (5.7 vs. 2.3%, P = 0.03, Supplementary Figure 3D). Interestingly, the percentage of individuals with multiple chronic medical conditions (>3 chronic medical conditions) was lower in the minority race/ethnicities such as Native American (36.8%), Hispanic (37.6%), and Black individuals (46.3%) than non-Hispanic white individuals (56.0%) (Table 1; P < 0.0001). The predominant presenting symptom was dyspnea (45%), and in descending order were cough (31%), fever (29%), abdominal complaints (11%), and myalgia (1.6%; Table 1). Excepting for fever, Native Americans had lower rates of other symptomatology than non-Hispanic white individuals (all P < 0.05; Table 1).

Characteristics of COVID-19 outcomes by race and age using bivariate analysis

Sepsis, ICU admission, mechanical ventilation, and hospital mortality

Overall, in-hospital mortality rate for COVID-19 patients was 16.5% (Table 2) and progressively increased across age (from 3% in 18–39 years of age to 36.4% in individuals ≥80 years; 18–29 vs. 40–59, chi-square = 88.9; 40–59 vs. 60–79, chi-square=88.9, 60–79 vs. ≥80, chi-square=88.9; all Bonferroni-adjusted P < 0.0001, Figure 1A). Native Americans had the highest rate of ICU admissions (43.3%) among all race/ethnicity groups (Native vs. Black, chi-square = 23.9; Native vs. Hispanic, chi-square = 23.9; Native vs. white, chi-square=23.9; Native vs. Others, chi-square=23.9; all Bonferroni-adjusted P < 0.0001, Table 2). This pattern was consistent across all age categories (all P < 0.0001, Figure 1B). Approximately 50% of Native American patients >60 years were admitted to the ICU. Greater than 15% of patients needed invasive mechanical ventilation, and both Native Americans (28.1%) and Hispanic individuals (15.9%) had higher rates of invasive mechanical ventilation than non-Hispanic white individuals (12.7%). Native Americans belonging to the age category (60–79 years) had the highest rate of mechanical ventilation compared with any other race/ethnicity age group categories (P < 0.05, Figure 1C; Supplemetnary Table 7). Native Americans had higher rates of sepsis than Black individuals (chi-square = 2.11, Bonferroni-adjusted P < 0.001), Hispanic individuals, or non-Hispanic white individuals (chi-square = 2.11, Bonferroni-adjusted P < 0.0001) (Table 2; Figure 1D). Approximately 50% of Native American patients >40 years developed sepsis (Figure 1D).

TABLE 2

Table 2. Characteristics of COVID-19 patients' outcomes by race.

FIGURE 1

Figure 1. Distributions of clinical outcomes of COVID-19 by age and race in hospitalized population. LOS was examined in live inpatients; A, Native American; B, Black and African American; H, Hispanic; O, Others; W, White. (A) Mortality by Age and Race. (B) ICU admission by Age and Race. (C) Ventilation by Age and Race. (D) Sepis by Age and Race. (E) Admission through ED by Age and Race. (F) LOS by Age and Race.

Admission through ED and hospital stay

Approximately 82% of hospitalized patients were admitted through the ED (85.2% Hispanic individuals, 87.7% Black individuals, 85.5% white individuals, and 43.5 Native American individuals) (Table 1). Native American individuals were less likely to be admitted through the ED or conversely more likely to undergo inter-hospital transfers (P-values <0.0001, Figure 1E). The median hospital LOS was 6.4 days for all inpatients (Table 1). LOS in Native American individuals was significantly longer than in non-Hispanic white individuals (Table 1). Native American individuals who were <80 years old had a longer LOS than non-Hispanic white individuals (all P < 0.0001; Table 3), and Native American patients who were between 40 and 69 years of age had the longest LOS of any age group (all P < 0.05, Figure 1F).

TABLE 3

Table 3. Multivariable logistic and linear regression results for factors associated with hospital death, mechanical intervention, ICU admission, acute complication as sepsis, LOS, admitted through ED in hospitalized patients.

Factors associated with COVID-19 outcomes using multivariable analyses

Sepsis, ICU admission, mechanical ventilation, and hospital mortality

Multivariable analyses indicated that after controlling for other factors, including a proportion of inter-hospital transfers, Native American patients had a greater chance of dying in the hospital (OR, 2.3; 95%CI, 1.6–3.2), receiving invasive mechanical ventilation (OR, 2.3; 95%CI, 1.8–3.0), being admitted to the ICU (OR, 2.4; 95%CI, 1.9–3.0), and developing sepsis (OR, 1.7; 95%CI, 1.4–2.8) (all P < 0.05; Table 3). Hispanic patients had a greater risk of hospital mortality (OR, 1.4; 95%CI, 1.1–1.7), mechanical ventilation (OR, 1.4; 95%CI, 1.2–1.6), and sepsis (OR, 1.2; 95%CI, 1.1–1.3) than non-Hispanic white individuals (all P < 0.05; Table 3). Male patients had an increased risk for hospital death (OR, 1.7; 95%CI, 1.4–2.0), invasive mechanical ventilation (OR, 2.0; 95%CI, 1.7–2.3), ICU admission (OR, 1.7; 95%CI, 1.5–1.9), and sepsis (OR, 1.4; 95%CI, 1.3–1.6) (all P < 0.0001; Table 3). Older age only significantly increased the odds of hospital death (OR, 1.02; 95%CI, 1.01–1.04; P < 0.0001; Table 3). Having Medicaid or Medicare compared to having private insurance was not significantly associated with any outcomes. However, residence in high-income areas was associated with lower odds of hospital death (OR, 0.4; 95%CI, 0.2–0.9), ventilation (OR, 0.2; 95%CI, 0.1–0.5), and ICU admission (OR, 0.3; 95%CI, 0.27–0.6) (all P < 0.05; Table 3). Chronic conditions such as lung diseases, diabetes, ventricular fibrillation, obesity, and symptoms such as dyspnea were associated with increased odds of hospital death, ventilation, ICU admission, and sepsis (all OR > 1 and P < 0.05). Patients admitted through the ED had lower odds of ventilation (OR, 0.63; 95%CI, 0.52–0.77; P < 0.0001) and ICU admission (OR, 0.69; 95%CI, 0.59–0.81; P < 0.0001).

Hospital stay and admission through ED

Native American race and Hispanic ethnicity were associated with longer LOS by 2.06 days (coefficient, 2.06; 95%CI, 1.16–2.96; P < 0.0001) and 0.48 days (coefficient, 0.48, 95%CI, 0.06–0.96; P < 0.0001) compared with non-Hispanic white individuals. Medicaid insurance was not associated with LOS when compared to commercial plans while increasing income was associated with shorter LOS (coefficient, −3.0; 95%CI, −5.2–0.8), P < 0.01; Table 3; P < 0.01). Male sex (coefficient,1.2; 95%CI, 0.8–1.6, P < 0.0001), older age (coefficient,0.1; 95%CI, 0.1–0.1; P < 0.0001), chronic lung conditions (coefficient, 11.9; 95%CI, 9.7–14.0), and kidney conditions (coefficient; 1.5; 95%CI, 0.9–2.1) increased LOS (all P < 0.0001; Table 3). However, admission through the ED was associated with shorter LOS (coefficient, −1.6; 95%CI, −2.2—1.0; P < 0.0001; Table 3).

Factors associated with ED admission were also investigated. Native American patients had lower odds of admission through the ED than white individuals (OR, 0.17; 95%CI, 0.13–0.22; P < 0.0001), while Black individuals had higher odds (OR,1.69; 95%CI, 1.21–2.36, P < 0.01) than non-Hispanic white individuals (Table 3).

ICU stay and time to ICU

Within the ICU subgroup, chronic lung diseases (coefficient, 4.4; 95%CI, 0.0–8.8) and presenting symptoms such as fever (coefficient,1.7; 95%CI, 0.1–3.3) were associated with a longer ICU LOS (all P < 0.05; Table 4). Having Medicare (coefficient, −2.3; 95%CI, −4.2–0.4; P < 0.05) and being admitted through the ED (coefficient, −2.6; 95%CI, −3.9–1.47; P < 0.0001) were associated with a shorter ICU LOS (Table 4). Native American race was highly associated with a shorter time from registration to ICU admission after controlling for all covariates including being admitted to the ICU through the ED (coefficient, −20.4; 95%CI, −35.4–5.3; P < 0.01; Table 4).

TABLE 4

Table 4. Multivariable linear regression results for factors associated with ICU length of stay and time from registration to ICU admission.

Discussion

This retrospective cohort study characterized the COVID-19 health disparities and patient outcomes in Native American individuals within a large healthcare system in the Southwest US. Although many small-scale studies have investigated similar health disparities, to our knowledge, this is the first comprehensive study with sufficient representation of the Native American population that was performed to measure case-fatality rates. Our results showed that Native American, Black individuals, and Hispanic patients were disproportionately burdened with a greater proportion of hospitalization when compared to race/ethnicity distribution based upon census data (24). These findings are consistent with findings in other studies and may be attributable to high infection rates in these minorities secondary to social determinants of health and other reasons (25, 26). However, in such studies mortality was based upon data from the National Death Index of the CDC, and case-fatality rates and hospital course of these individuals were not assessed (12). Conceivably, a greater number of infections in the minority populations could have explained the greater mortality in Native American individuals. Our findings align with a recent report that used Agency for Healthcare Research and Quality (AHRQ) 2020 NIS dataset for hospitalizations from 1 January 2020 to 31 December 2020 to show that Native American men had the highest inpatient mortality due to COVID-19 (11). To our knowledge, our study is the first to demonstrate a greater case-fatality rate in a hospitalized setting among Native Americans irrespective of gender after adjusting for multiple covariates.

In our study, the greater prevalence of chronic lung disease (OR 6.14; 95%CI; 3.91–9.63; Table 3) may be responsible for the observed greater hospital mortality in Native American individuals when compared to non-Hispanic white individuals or even other minorities. Previously, investigators have found that patients with chronic lung disease (e.g., chronic obstructive pulmonary disease) had the highest odds of COVID-19 mortality compared with patients with other pre-existing cardiometabolic conditions, such as obesity, diabetes, and hypertension (27). Moreover, the Native American populations have a higher adjusted prevalence of chronic respiratory disease compared with the non-Hispanic white population (28). Such differences in the presence of chronic lung disease as comorbidity are important information that healthcare providers and systems need to consider when caring for Native American individuals or even other patients.

Information on racial disparities in COVID-19 outcomes has largely come from studies on the Black population, with minimal research on the impact on Native American individuals (12). Some studies reveal that Native American individuals experienced disproportionate infection and mortality rates during the COVID-19 pandemic (10). Our study furthers our understanding of such findings by demonstrating the greater likelihood of hospital death, invasive ventilation, ICU admission, and sepsis in Native American individuals, even after adjusting for differences in sociodemographic factors (insurance status and median household income) and other clinical characteristics. Interestingly, while older age and male sex are commonly reported risk factors for adverse COVID-19 outcomes, the median age of Native American patients in our study was 17 years younger than in non-Hispanic white individuals with a near 1:1 ratio between men and women (29). Thus, age and sex did not contribute to the observed adverse COVID-19 outcomes in Native American individuals.

Observed associations between Native American race and COVID-19 outcomes may be multifactorial. One factor may be racial differences in immune response (30). Studies have found that Native American individuals have genetic susceptibility and subdued innate response to novel viral infection (31). Another factor may be related to the lower rates of common COVID-19 symptoms in Native American individuals than non-Hispanic white individuals. The racial difference in symptom presentation may reflect a difference in the underlying immune response or end-organ response to COVID-19 and may lead to delays in seeking medical attention. Such atypical clinical presentation and symptomatology of COVID-19 in Native American patients may increase the difficulty for patients and healthcare providers to identify COVID-19 and monitor clinical changes. Such factors may also underlie greater case-fatality rates in Native American individuals, highlighting the necessity of increased precautions when caring for Native American patients with COVID-19. As COVID-19 is still a novel viral respiratory infection, its epidemiology and clinical responses are yet unknown among different racial groups, and more studies are needed to assess racial differences in immune response and clinical presentation.

Factors such as social determinants of health may also contribute to the observed racial differences in COVID-19 outcomes. One in four Native American individuals is in poverty with a median annual income below the national average and is more likely to have chronic diseases (7, 32). The median income of residence areas for Native American patients was significantly lower than white individuals, and a greater share of Native American individuals had Medicaid than any other racial group, with 79% of Native American individuals between ages 18 and 39 on Medicaid. Consistent with national surveys, Native American inpatients had higher rates of diabetes, obesity, and chronic lung disease than non-Hispanic white individuals (7). Studies in other races have reported associations between social inequalities in housing, education level, and jobs, and the susceptibility to COVID-19 (33). Many long-standing socioeconomic issues such as living conditions, transportation, employment, access to health services, and running water affect Native American individuals as well (34). These health and socioeconomic determinants may have made Native American individuals more vulnerable to COVID-19. However, no strong association was found between the Native American race and the number of comorbidities in our study, although individual conditions such as chronic lung and liver disease were associated with greater risk for mortality. Such findings further underscore the importance of understanding organ-specific mechanisms—such as lung and liver—that may be responsible for the observed difference. Interestingly, certain single nucleotide polymorphisms (SNPs; rs2070788 and rs383510) that were associated with severe pulmonary damage caused by influenza A(H7N9) in 2014 and influenza A(H1N1) in 2009 are surprisingly the highest in Native American individuals in a worldwide study (35).

In our study, only 43% of Native American individuals were admitted through the emergency department (ED), whereas more than 80% of Black, Hispanic and non-Hispanic white individuals were admitted through ED. The low rate of hospitalization through the ED was carried across all age cohorts of Native American individuals and remained independent of sociodemographic factors and clinical presentation. Previous studies have found that older individuals, men, Black individuals, and low-income populations had higher rates of visiting the ED (36). In keeping with such findings, our study showed that the older adults, men, Black individuals, individuals from lower socioeconomic status, and individuals with chronic liver disease, hypertension, dementia, depression, and symptoms were more likely to be admitted through the ED. But the finding that a majority of Native American individuals were inter-hospital transfers may be reflective of barriers to accessing higher level care and adjusted for in our analyses of other outcomes such as mortality and time to ICU admission. However, it is well known that inter-hospital transfer is independently associated with greater in-hospital mortality, especially from rural areas (37). While such findings suggest issues with access to care, delayed presentation, and destabilization during limited treatment availabilities during transit, we observed that the greater in-hospital mortality in Native American individuals when compared to non-Hispanic white individuals remained even after adjusting for inter-hospital transfers (Table 3). Such findings suggest that other factors may be responsible for the observed greater mortality in Native American individuals.

We investigated the relationship between time from hospital admission to ICU transfer as a surrogate measure of the rate of COVID-19 progression and clinical deterioration. Our multivariable analysis revealed that the time to ICU transfer in Native American individuals was shorter than in other racial groups across all age groups. These findings suggest that Native American individuals with COVID-19 progressed more rapidly than others. We recognize that we do not have time spent by patients at outside hospitals when they underwent inter-hospital transfers, and that is a limitation. However, despite adjusting for inter-hospital transfer status, we found that the time to deterioration was greater in Native American individuals. Moreover, an increased odd of admission from the ED to the ICU was noted in Native American individuals (Table 4), which was not affected by inter-hospital transfer status or time spent on the floor in the outside hospital. To our knowledge, our study is the first to investigate the clinical deterioration rates in COVID-19 inpatients and higher case-fatality rates for various race and ethnicity groups. In combination with the lower rate of commonly reported symptoms of COVID-19 in Native American individuals, these results may emphasize the importance of closer monitoring of Native American patients upon entry or admission to the hospital.

Our study has several limitations, the first being the retrospective observational nature of the data. Therefore, adjustment for all confounders was not possible. Our study sample has a better representation of Native American populations because Banner Health hospitals serve as the tertiary care referral centers for many critical access hospitals in the southwestern United States. Therefore, it is possible that we missed a large proportion of the Native American population, hence limiting the generalizability of our findings. However, our large sample size and detailed phenotyping allowed adjustment for multiple covariates.

Public health implications

The greater case-fatality rates following hospitalization in Native American individuals emphasize the systematic disadvantage experienced by Native American individuals and point out the need to address health inequalities by developing programs that improve overall health outcomes. Atypical symptom presentation of COVID-19, greater prevalence of chronic lung disease, and more rapid clinical deterioration than other races/ethnicities underscore the role of pulmonologists in addressing such disparities. The association between greater COVID-19-related mortality and chronic lung disease is being addressed by improving COVID-19 vaccination rates in adults in the pulmonary ambulatory setting (38–40).

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.

Ethics statement

The studies involving human participants were reviewed and approved by Banner Health IRB 483-20-0076. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author contributions

CB, YW, GC, DS, SK, PT, and SP participated in the conception of the work, planning, data collection, analysis and manuscript writing and review of final draft. VK participated in data analysis, manuscript writing and review of final draft. All authors contributed to the article and approved the submitted version.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2023.1220582/full#supplementary-material

References

1. COVID-19. (2021). Available online at: https://www.cdc.gov/coronavirus/2019-ncov/index.html (accessed July 8, 2021).

Google Scholar

2. Cardemil CV, Dahl R, Prill MM, Cates J, Brown S, Perea A, et al. COVID-19-related hospitalization rates and severe outcomes among veterans from 5 veterans affairs medical centers: hospital-based surveillance study. JMIR Public Health Surveill. (2021) 7:e24502. doi: 10.2196/24502

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Peek ME, Simons RA, Parker WF, Ansell DA, Rogers SO, Edmonds BT. COVID-19 among african Americans: an action plan for mitigating disparities. Am J Public Health. (2021) 111:286–92. doi: 10.2105/AJPH.2020.305990

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Bae S, Kim SR, Kim MN, Shim WJ, Park SM. Impact of cardiovascular disease and risk factors on fatal outcomes in patients with COVID-19 according to age: a systematic review and meta-analysis. Heart. (2021) 107:373–80. doi: 10.1136/heartjnl-2020-317901

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Yancy CW. COVID-19 and African Americans. JAMA. (2020) 323:1891–2. doi: 10.1001/jama.2020.6548

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Khazanchi R, Beiter ER, Gondi S, Beckman AL, Bilinski A, Ganguli I. County-level association of social vulnerability with COVID-19 cases and deaths in the USA. J Gen Intern Med. (2020) 35:2784–7. doi: 10.1007/s11606-020-05882-3

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Adamsen C, Schroeder S, LeMire S, Carter P. Education, income, and employment and prevalence of chronic disease among American Indian/Alaska native elders. Prev Chronic Dis. (2018) 15:E37. doi: 10.5888/pcd15.170387

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Forrest KYZ, Leeds MJ, Ufelle AC. Epidemiology of obesity in the Hispanic adult population in the United States. Fam Community Health. (2017) 40:291–7. doi: 10.1097/FCH.0000000000000160

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Burki T. COVID-19 among American Indians and Alaska Natives. Lancet Infect Dis. (2021) 21:325–6. doi: 10.1016/S1473-3099(21)00083-9

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Arrazola J, Masiello MM, Joshi S, Dominguez AE, Poel A, Wilkie CM, et al. COVID-19 mortality among American Indian and Alaska Native Persons - 14 States, January-June 2020. MMWR Morb Mortal Wkly Rep. (2020) 69:1853–6. doi: 10.15585/mmwr.mm6949a3

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Pal S, Gangu K, Garg I, Shuja H, Bobba A, Chourasia P, et al. Gender and Race-Based Health Disparities in COVID-19 outcomes among hospitalized patients in the united states: a retrospective analysis of a national sample. Vaccines. (2022) 10:2036. doi: 10.3390/vaccines10122036

PubMed Abstract | CrossRef Full Text | Google Scholar

12. 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:362–73. doi: 10.7326/M20-6306

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Mangla A, Agarwal N. Clinical Practice Issues in American Indians and Alaska Natives. StatPearls. Treasure Island (FL) ineligible companies. Disclosure: Nikki Agarwal declares no relevant financial relationships with ineligible companies. StatPearls Publishing Copyright © 2023. Treasure Island: StatPearls Publishing LLC. (2023).

PubMed Abstract | Google Scholar

14. Skewes MC, Gonzalez VM, Gameon JA, FireMoon P, Salois E, Rasmus SM, et al. Health disparities research with American Indian communities: the importance of trust and transparency. Am J Community Psychol. (2020) 66:302–13. doi: 10.1002/ajcp.12445

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Adakai M, Sandoval-Rosario M, Xu F, Aseret-Manygoats T, Allison M, Greenlund KJ, et al. Health disparities among American Indians/Alaska Natives - Arizona, 2017. MMWR Morbid Mortality Weekly Report. (2018) 67:1314–8. doi: 10.15585/mmwr.mm6747a4

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Disparities NIoMHa,. Minority Health Health Disparities: Definitions Parameters. (2022). Available online at: https://www.nimhd.nih.gov/about/strategic-plan/nih-strategic-plan-definitions-and-parameters.html (accessed August 28, 2022).

Google Scholar

17. Control CfD,. Health Disparities. (2022). Available online at: https://www.cdc.gov/healthyyouth/disparities/index.htm (accessed August 28, 2022).

Google Scholar

18. Ignacio M, Oesterle S, Mercado M, Carver A, Lopez G, Wolfersteig W, et al. Narratives from African American/Black, American Indian/Alaska Native, and Hispanic/Latinx community members in Arizona to enhance COVID-19 vaccine and vaccination uptake. J Behav Med. (2022) 46:140–152. doi: 10.1007/s10865-022-00300-x

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Group A-3L-CS, Lundgren JD, Grund B, Barkauskas CE, Holland TL, Gottlieb RL, et al. A Neutralizing Monoclonal Antibody for Hospitalized Patients with Covid-19. N Engl J Med. (2021) 384:905–14. doi: 10.1056/NEJMoa2033130

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Beigel JH, Tomashek KM, Dodd LE, Mehta AK, Zingman BS, Kalil AC, et al. Remdesivir for the Treatment of Covid-19 - Final Report. N Engl J Med. (2020) 383:1813–26. doi: 10.1056/NEJMoa2007764

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Lederer DJ, Bell SC, Branson RD, Chalmers JD, Marshall R, Maslove DM, et al. Control of confounding and reporting of results in causal inference studies. Guidance for Authors from Editors of Respiratory, Sleep, and Critical Care Journals. Ann Am Thorac Soc. (2019) 16:22–8. doi: 10.1513/AnnalsATS.162erratum

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Hosmer DW, Lemeshow S, Sturdivant RX. Applied Logistic Regression, 3rd Edition. NY: Wiley, Inc. (2013). doi: 10.1002/9781118548387

CrossRef Full Text | Google Scholar

23. Neter J, Nachtscheim CJ, Li. W. Applied Linear Statistical Models, 5th Ed. New York: Mc-Graw Hill. (2004).

Google Scholar

24. United States Census Bureau Quick Facts (2022). Available online at: https://www.census.gov/quickfacts/US (accessed July 7, 2022).

Google Scholar

25. Muñoz-Price LS, Nattinger AB, Rivera F, Hanson R, Gmehlin CG, Perez A, et al. Racial Disparities in Incidence and Outcomes Among Patients With COVID-19. JAMA Netw Open. (2020) 3:e2021892. doi: 10.1001/jamanetworkopen.2020.21892

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Price-Haywood EG, Burton J, Fort D, Seoane L. Hospitalization and mortality among black patients and white patients with COVID-19. N Engl J Med. (2020) 382:2534–43. doi: 10.1056/NEJMsa2011686

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Conway FM, Bloom CI, Shah PL. Susceptibility of patients with airways disease to SARS-CoV-2 infection. Am J Respir Crit Care Med. (2022) 206:696–703. doi: 10.1164/rccm.202111-2547PP

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Laffey KG, Nelson AD, Laffey MJ, Nguyen Q, Sheets LR, Schrum AG. Chronic respiratory disease disparity between American Indian/Alaska Native and white populations, 2011-2018. BMC Public Health. (2021) 21:1466. doi: 10.1186/s12889-021-11528-8

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Rosenthal N, Cao Z, Gundrum J, Sianis J, Safo S. Risk factors associated with in-hospital mortality in a us national sample of patients with COVID-19. JAMA Netw Open. (2020) 3:e2029058. doi: 10.1001/jamanetworkopen.2020.29058

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Nédélec Y, Sanz J, Baharian G, Szpiech ZA, Pacis A, Dumaine A, et al. Genetic ancestry and natural selection drive population differences in immune responses to pathogens. Cell. (2016) 167:657–669.e621. doi: 10.1016/j.cell.2016.09.025

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Kakol M, Upson D, Sood A. Susceptibility of Southwestern American Indian tribes to coronavirus disease 2019 (COVID-19). J Rural Health. (2021) 37:197–9. doi: 10.1111/jrh.12451

PubMed Abstract | CrossRef Full Text | Google Scholar

32. New New Detailed Statistics on Race Hispanic Origin Ancestry Tribal Groups. (2017). Available online at: https://www.census.gov/newsroom/press-releases/2017/acs-selected-population-tables-aian.html (accessed July 7, 2022).

Google Scholar

33. Adegunsoye A, Ventura IB, Liarski VM. Association of black race with outcomes in COVID-19 disease: a retrospective cohort study. Ann Am Thorac Soc. (2020) 17:1336–9. doi: 10.1513/AnnalsATS.202006-583RL

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Rodriguez-Lonebear D, Barcelo NE, Akee R, Carroll SR. American Indian reservations and COVID-19: correlates of early infection rates in the pandemic. J Public Health Manag Pract. (2020) 26:371–7. doi: 10.1097/PHH.0000000000001206

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Kehdy FSG, Pita-Oliveira M, Scudeler MM, Torres-Loureiro S, Zolini C, Moreira R, et al. Human-SARS-CoV-2 interactome and human genetic diversity: TMPRSS2-rs2070788, associated with severe influenza, and its population genetics caveats in Native Americans. Genet Mol Biol. (2021) 44:e20200484. doi: 10.1590/1678-4685-gmb-2020-0484

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Dufour I, Chouinard MC, Dubuc N, Beaudin J, Lafontaine S, Hudon C. Factors associated with frequent use of emergency-department services in a geriatric population: a systematic review. BMC Geriatr. (2019) 19:185. doi: 10.1186/s12877-019-1197-9

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Mohr NM, Harland KK, Shane DM, Ahmed A, Fuller BM, Torner JC. Inter-hospital transfer is associated with increased mortality and costs in severe sepsis and septic shock: An instrumental variables approach. J Crit Care. (2016) 36:187–94. doi: 10.1016/j.jcrc.2016.07.016

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Ortiz JR, Hartert TV. The Role American Thoracic Society Healthcare Providers Have in Immunization. Am J Respir Crit Care Med. (2021) 203:944–5. doi: 10.1164/rccm.202010-3990ED

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Hastert FD, Hein S, von Rhein C, Benz NI, Husria Y, Oberle D, et al. The SARS-CoV-2 variant omicron is able to escape vaccine-induced humoral immune responses, but is counteracted by booster vaccination. Vaccines (Basel). (2022) 10:794. doi: 10.3390/vaccines10050794

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Hastert FD, Henss L, von Rhein C, Gerbeth J, Wieters I, Borgans F, et al. Longitudinal analysis of coronavirus-neutralizing activity in COVID-19 patients. Viruses. (2022) 14:882. doi: 10.3390/v14050882

PubMed Abstract | CrossRef Full Text | Google Scholar