Skip to main content

Racial/ethnic disparities on inflammation and response to methylprednisolone in severe COVID-19 pneumonia

Abstract

Background

Racial/ethnic minorities are at higher risk for severe COVID-19. This may be related to social determinants that lead to chronic inflammatory states. The aims of the study were to determine if there are racial/ethnic disparities with inflammatory markers and association of methylprednisolone to in hospital survival.

Methods

This was a secondary analysis of a retrospective cohort study of patients ≥ 18 years of age and admitted for severe COVID-19 pneumonia between March and June 2020 in 13 Hospitals in New Jersey, United States. Patients who received other formulation of corticosteroids were not included. Area under the receiver operating characteristics curves were performed to test for discriminatory ability of each inflammatory makers. Univariate and multivariate Cox regression assessed the association of variables to in hospital survival.

Results

Propensity matched sample (n = 759) between no methylprednisolone (n = 380) and methylprednisolone (n = 379) had 338 Whites, 102 Blacks, 61 Asian/Indians, and 251 non-Black non-White Hispanics. Compared to CRP, area under receiving operating characteristic curve for d-dimer in Hispanics (0.742) was statistically different (DeLong Test P = 0.0041). Multivariate cox regression showed that different variables in Blacks [age ≥ 60 years (HR = 3.71, P = 0.0281), mechanical ventilation (HR = 5.07, P = 0.0281) and creatinine ≥ 1.5 mg/dL (HR = 3.61, P = 0.0007)], Whites [cancer (HR = 1.68, P = 0.0213), qSOFA score of 1 (HR = 1.81, P = 0.0213), qSOFA score of 2 (HR = 5.16, P < 0.0001), qSOFA score of 3 (HR = 11.81, P < 0.0001) and creatinine ≥ 1.5 mg/dL (HR = 2.16, P = 0.0006)], Hispanics [hypertension (HR = 2.52, P = 0.0007), cancer (HR = 2.99, P = 0.0244 and D-dimer ≥ 2 mcg/mL (HR = 2.22, P = 0.0077)], and Asian/Indians [ chronic kidney disease (HR = 6.36, P = 0.0031) and CRP > 20 mg/L (HR = 5.02, P = 0.0032)] were statistically significant for mortality. Low dose and high dose methylprednisolone were significantly associated with prolonged survival in Whites [low dose (HR = 0.37, P < 0.0001) and high dose (HR = 0.48, P < 0.0183)] and Asian/Indians [low dose (HR = 0.13, P = 0.0101) and high dose (HR = 0.15, P = 0.01)]. However, high dose was not associated with improved survival compared to low dose. Methylprednisolone was not associated with prolonged survival in Blacks and Hispanics.

Conclusion

Racial/Ethnic disparities with inflammatory markers preclude the use of one marker as a predictor of survival. Methylprednisolone is associated with prolonged survival in Asian/Indians and Whites.

Peer Review reports

Introduction

As of July 2, 2021, there has been 182,319,261 coronavirus disease 2019 (COVID-19) cases and 3,954,324 deaths in worldwide [1], and the main cause of mortality is hyperinflammatory acute respiratory distress syndrome (ARDS). Advanced age, diabetes, cardiovascular disease, chronic lung disease, chronic kidney disease (CKD), and racial/ethnic minorities are among the factors that appear to increase the risk for severe COVID-19 [2,3,4]. Racial/ethnic minorities include non-black non-white Hispanics, Blacks, Native Americans, Native Hawaiians and other Pacific Islanders [5]. After age-related adjustments, mortality in Blacks, non-white non-black Hispanics, and Asians are higher compared to Whites in the United States, United Kingdom, and Brazil [1, 6,7,8].

One explanation is social determinants (multi-generational homes, essential workers, low socio-economic status, lack of access to quality health care) predispose these minorities to higher COVID-19 exposure, and behaviors (depression, anxiety, smoking, alcoholism, high sugar, salt and fatty diet) that lead to chronic inflammatory states. These might influence the severe clinical presentation [5]. Inflammatory markers, such as IL-6 (> 25 pg/mL), D-dimer (≥ 2.0 mcg/mL), CRP (≥ 20 mg/L) and/or ferritin ≥ 1400 μg/L, are believed directly correlate with mortality in COVID-19 [9,10,11,12,13,14,15,16]. There are racial/ethnic variability with inflammatory markers [17] but there has been no study comparing the variability between survivors and non-survivors.

Histologic studies show the sequela of this inflammation, with severe endothelial damage, diffuse alveolar damage, thrombosis in situ, intussusceptive angiogenesis, and steroid responsive histologic patterns of organizing pneumonia (OP) and/or acute fibrosing organizing pneumonia (AFOP) [18, 19]. Dexamethasone has been shown to improve mortality in patients requiring oxygen support including invasive mechanical ventilation [20]. Our initial study suggested that low dose methylprednisolone (< 1.36 mg/kg/day) given > 7 days from onset of symptoms for 7 days were associated with improved mortality and no additional benefit with duration > 14 days or high dose (≥ 1.36 mg/kg/day) [21]. Racial/ethnic minorities have a higher incidence of chronic inflammation.

It is unknown if these groups require higher doses of methylprednisolone.

Methodology

Eligibility criteria

Real world data was collected from Hackensack Meridian Health (HMH), a NJ health network comprising of 13 hospitals on patients ≥ 18 years of age, and hospitalized for at least 2 days between March 1, 2020 and June 15, 2020 with severe COVID-19 Pneumonia [21]. These patients had a positive SARS-CoV-2 PCR and had SpO2 < 94% on room air at sea level, a respiratory rate > 30 breaths/min, PaO2/FiO2 < 300 mm Hg, or lung infiltrates > 50%. We excluded patients who had different corticosteroid regimen other than methylprednisolone. Approval was obtained by the Hackensack Meridian Health Institutional Review Board (study #Pro2020-0485) and the study was also registered on ClinicalTrials.Gov as a prospective observational database (NCT04347993).

Data collection process and data items

Demographic data such as age, gender, race, ethnicity, comorbidities, and sex were self-reported. Weight and height were measured. SARS-CoV-2 was detected in nasal swabs by RT-PCR. Routine blood tests included complete blood count (CBC), coagulation profile, complete metabolic profile (CMP), inflammatory markers [interleukin-6 (IL-6), C reactive protein (CRP), d-dimer, and ferritin], and arterial blood gas(ABG). Data was entered into Redcap and abstracted from June to December 2020.

Outcomes

The primary outcomes are the levels of the inflammatory markers of survivors and non-survivors and the association of methylprednisolone dose to in hospital survival for each racial/ethnic group.

Statistical analysis

The data in this study is primarily the propensity score matched cohort of comparison of in-hospital survival of COVID-19 admitted patients that were treated with methylprednisolone and patients not treated with methylprednisolone. [21] In the one-to-one propensity score matched design of the parent study, patients from the methylprednisolone administration groups were matched based on variables associated with mortality such as age (age ≥ 60 years vs. age < 60 years), obesity (BMI ≥ 30.0 kg/m2 vs. BMI < 30.0 kg/m2), sex (M/F), diabetes (Yes/No), hypertension (Yes/No), cancer (Yes/No), respiratory Rate (respiratory rate > 22 bpm vs. ≤ 22 bpm), chronic kidney disease or chronic renal failure (Yes/No), low oxygen (oxygenation > 94% vs oxygenation ≥ 94%), CRP (CRP > 20 mg/dL vs. CRP ≤ 20 mg/dL), and qSOFA (score: 0,1,2,3) [21]. A nearest-neighbor method (greedy match) was employed using a caliper of 0.20 to obtain the propensity matched sample. A post-match assessment of the distribution of propensity scores (or logit of propensity scores) and balance in the adjusted variables between the no methylprednisolone (NMP) and methylprednisolone (MP) using standardized difference and variance ratio were performed. This analysis and graphical displays were obtained by the ASSESS statement of PROC MATCH in SAS 9.4 [21]. The aims of this subgroup analysis were pursued between and within the racial/ethnic groups that exhausted the span of racial/ethnic groups of the propensity score-matched cohort from primary study [21].

In this study descriptive statistics was used to summarize the data. Categorical variables were presented as frequencies and percentages. Continuous variables were presented as the median and interquartile range (IQR) for non-normally distributed data and mean and standard deviation for normally distributed data. The normality assumption of continuous variables was assessed using Shapiro–Wilk test. Comparison of matched continuous variables was performed using two-sided paired t-test or Wilcoxon signed-rank test, as appropriate. Comparison of paired categorical variables was performed using McNemar’s test or Bowker test of symmetry, as appropriate.

Comparison of levels of inflammatory markers, creatinine and age between racial/ethnic groups, which are independent, was performed using Kruskal–Wallis test followed by pairwise comparisons utilizing a two-sided Wilcoxon rank sum test.

Comparison of these continuous variables between survivors and non-survivors within each racial/ethnic group was performed using two-sided Wilcoxon rank sum test. Boxplots were used to illustrate the different distribution between survivor and non-survivors in the inflammatory markers. To assess the discriminative ability of inflammatory markers for in-hospital mortality, area under the receiver operating characteristics (ROC) curve of the inflammatory markers in each racial/ethnic group using ROC statements in PROC LOGISTIC SAS 9.4. Comparison of area under ROCs within each racial/ethnic group were evaluated by DeLong test, with ROC closest to 0.5 as the reference. Area under ROCs were reported as area (95% confidence interval). AUC closer to 1 with a 95% CI that excludes 1 was considered statistically significant.

Within each racial/ethnic group, this study also examined the association between in-hospital survival and inflammatory markers, methylprednisolone dose (none, low dose, high dose), older age obesity, sex, diabetes, hypertension cancer, respiratory rate (bpm), chronic kidney disease (or renal failure), supplemental oxygen support (None, non-mechanical ventilation, mechanical ventilation), CRP ≥ 20 mg/dL, qSOFA and creatinine ≥ 1.5 mg/dL. Methylprednisolone (MP) dose cutoff value of 1.36 mg/kg/day had been previously determined as optimal based on Youden index method such that low dose (LD) MP was defined as MP dose < 1.36 mg/kg/day and high dose (HD) MP was defined as MP dose ≥ 1.36 mg/kg/day. Dichotomized inflammatory marker we employed in all the time to in-hospital mortality event analysis as follows: IL-6 > 25 pg/mL, D-dimer ≥ 2 mcg/mL, CRP ≥ 20 mg/L and ferritin ≥ 1400 μg/L as high levels of each respective marker. Time to in-hospital mortality was estimated using Kaplan–Meier product limit method. Comparison of in-hospital survival between independent groups was estimated by a two-sided log-rank test and matched cohort were compared using stratified log rank. To examine the association of risk/benefit of in-hospital survival with the factor of interest such as methylprednisolone treatment, in a univariate model, the Cox proportional hazard regression analysis with robust covariance [22] (sandwich estimator) to account for paired observations was conducted. The proportional hazard (PH) assumption, critical in Cox regression modelling, was evaluated using a Kolmogorov-type supremum test [23] in the ASSESS statement of PROC PHREG. If the PH assumption was violated for any covariate, then a continuous variable which also violated the PH assumption and its interaction with time were included in the model to adjust for the significant interaction with time to the risk of mortality [24]. After conducting univariate PH Cox regression or NPH Cox regression, any covariate that reported P < 0.15 were entered into the initial full model of the multivariable Cox model to initiate backward elimination procedure. If any of covariates that satisfied the criterion of P < 0.15 above had been found to have violated PH assumption then a multivariate NPH Cox regression model was analysis, otherwise if all the variables satisfied the PH assumption then a standard multivariable PH Cox model was fit. Results of all univariate and multivariable Cox regression (PH/NPH) within each race/ethnic group were presented as hazard ratio (HR), 95% Cis and P-value. All data analysis was performed using SAS 9.4 (SAS Institute Inc., Cary, North Carolina, USA). Unless otherwise specified, any P < 0.05 was considered statistically significant.

Results

Between March 4 and June 15, 2020, 2041 patients were flagged in the electronic health record with a diagnosis of COVID-19 and pneumonia. A total of 539 patients were excluded based on eligibility criteria (< 18 years of age, pregnant, received other formulations of corticosteroids, or hospitalized for less than 2 days) Thus, 1072 patients had their data abstracted. A propensity score matched sample was constructed out of 759 patients (380 in NMP and 379 in MP). After an examination of the proportional hazard assumption, MP and Fractional inspired oxygen (FiO2) significantly violated it (both with P < 0.0001). [21] Data on P/F ratio was lacking; and FiO2 was used since 95% of patients had this data. The supremum test also indicated that non-proportionality was observed in other variables such as nursing home, lack of taste or smell, WBC < 11,000 cells/mL, creatinine > 1.5 ng/mL, respiratory rate > 22 bpm, hydroxychloroquine (HCQ), MP, HD or LD MP, calcium, and initial diastolic blood pressure.

All variables with non-proportional hazard were adjusted using FiO2, as indicated above. The Youden Index method yielded a MP dose cut-off 1.36 mg/kg/day. Low dose methylprednisolone (LDMP) was defined as < 1.36 mg/kg/day and high dose methylprednisolone (HDMP) was defined as ≥ 1.36 mg/kg/day. 215 received LDMP and 164 received HDMP [21]. In the original study, there was another category for racial and ethnic group in REDCAP. We manually reviewed the charts further for racial and ethnic group and other variables, and found there were 102 Blacks, 338 Whites, 61 Asian/Indians, 251 non-white/non-black Hispanics and 7 unknown racial/ethnic group (Table 1, Additional file 1: Tables S1, S2, Figs. S1, S2).

Table 1 Characteristics of hospitalized COVID-19 patients treated with or without methylprednisolone (n = 759)

Box plots were created to determine the distribution of each inflammatory marker between survivors, or those who are alive, and non-survivors, or those who expired (Fig. 1). For ferritin, Asian/Indians (median = 1048 μg/L) had elevated levels compared to Blacks (median = 743.0 μg/L), Whites (median = 777.5 μg/L) and Hispanics (median = 831.0 μg/L), however, this was not statistically significant (Kruskal–Wallis P = 0.1689). Ferritin values in Blacks (n = 72) were not significantly different between survivors (n = 21) and non-survivors (n = 51) (median: 745 vs 575 μg/L, P = 0.9280). Ferritin values in Whites (n = 338) were significantly different between survivors (n = 243) and non-survivors (n = 95) (median: 641 vs 987 μg/L, P = 0.0139). Ferritin values in Asian/Indians (n = 56) were significantly different between survivors (n = 38) and non-survivors (n = 18) (median: 1265 vs 418 μg/L, P = 0.0211). The ferritin values in Hispanics (n = 226), were not significantly different between survivors (n = 158) and non-survivors (n = 68) (median: 826.5 vs 897.5 μg/L, P = 0.0854).

Fig. 1
figure 1

Box plots of inflammatory markers between those who expired and those who are alive. Panel shows red corresponds to those who are alive or are survivors and blue corresponds to those who have expired or are non-survivors

D-dimer values in Blacks (median = 1.64 mcg/mL), Whites (median = 1.33 mcg/mL), Asian/Indians (median = 1.21 mcg/mL) and Hispanics (median = 1.11 mcg/mL) were not significantly different (Kruskal–Wallis P = 0.1240). D-dimer values in Blacks (n = 62) were not significantly different between survivors (n = 44) and non-survivors (n = 18) (median: 1.91 mcg/mL vs 1.53 mcg/mL, P = 0.9119). D-dimer values in Whites (n = 193), were significantly different between survivors (n = 133) and non-survivors (n = 60) was significant (median: 1.19 vs 1.76 mcg/mL, P = 0.0034). D-dimer values in Asian/Indians (n = 41), were not significantly different between survivors (n = 26) and non-survivors (n = 15) (median: 1.10 vs. 1.51 mcg/mL, P = 0.0784). D-dimer values in Hispanics (n = 172), were significantly different between survivors (n = 116) and non-survivors (n = 56) (median: 0.89 vs 2.24 mcg/mL, P < 0.0001).

CRP values in Blacks (median = 11.40 mg/mL), Asian/Indians (median = 12.85 mg/mL) and Hispanics (median = 14.32 mg/mL), were statistically significantly different compared to Whites (median = 9.47 mg/mL), (Kruskal–Wallis P < 0.0001). This result was driven by the significant difference between Whites and Hispanics adjusting for multiple testing using Hochberg method. CRP values in Blacks (n = 85) were not significantly different between survivors (n = 63) and non-survivors (n = 22) (median: 10.2 mg/mL vs 12.07 mg/mL, P = 0.3586). CRP values in Whites (n = 276) were significantly different between survivors (n = 183) and non-survivors (n = 93) (median: 7.71 mg/mL vs 12.30 mg/mL, P = 0.0004). CRP values in Asian/Indians (n = 60) were significantly different between survivors (n = 38) and non-survivors (n = 22) (median: 9.05 vs. 20.75 mg/mL, P < 0.0001). CRP values in Hispanics (n = 228) were not significantly different between survivors (n = 158) and non-survivors (n = 70) (median: 13.70 mg/mL vs 16.70 mg/mL, P = 0.0514).

IL-6 values were not significantly different between the groups (Kruskal–Wallis P = 0.9843). Except for Asian/Indians with median IL-6 of 11 pg/mL, all the other groups had a median IL-6 of 14.0 pg/mL. IL-6 values in Blacks (n = 51) were not significantly different between survivors (n = 34) and non-survivors (n = 17) (median: 13 pg/mL vs 11 pg/mL, P = 0.7868). IL-6 values in Whites (n = 157) were significantly different between survivors (n = 94) and non-survivors (n = 63) (median: 11.0 pg/mL vs 27.0 pg/mL, P = 0.0005). IL-6 values in Asian/Indians (n = 61) were not statistically different between survivors (n = 38) and non-survivors (n = 23) (median 30.0 pg/mL vs 8.0 pg/mL, P = 0.0924). IL-6 values in Hispanics (n = 110) were statistically different between survivors (n = 73) and non-survivors (n = 37) (median = 19.0 pg/mL vs 11.0 pg/mL, P = 0.0168).

ROC analysis was used to identify the candidate inflammatory markers to include in the subsequent univariate and multivariable analysis of in-hospital survival in each racial/ethnic group (Fig. 2). Since IL-6 had the most missing except for the Asian/Indian cohort, we excluded from comparative ROC analysis as the Proc LOGISTIC used to perform the AUC analysis always seeks the most common subjects between measurements in a DeLong test. Thus this analysis was limited to D-dimer, CRP and Ferritin. In the DeLong testing, we used the AUC closest to 0.5 and/or the AUC whose 95% CI included 0.5 as a reference inflammatory marker. The inflammatory marker which was significantly different from the reference by DeLong test was automatically selected as a candidate for Univariate/Multivariate Cox models. A second marker that was the second largest in the AUC value was also included regardless of lack of statistical difference from the reference marker, In cases where none of the AUCs of the markers were statistically different from the reference, we forced the markers with two largest AUCs to initiate the model fitting steps.

Fig. 2
figure 2

Area under receiver operator characteristic curves of inflammatory markers in each racial/ethnic group. Panel shows AUC for each inflammatory marker

In Blacks, the AUCs for D-dimer (Area = 0.543, 95% CI 0.351–0.735), CRP (Area = 0.514, 95% CI 0.312–0.716), and Ferritin (Area = 0.585, 95% CI 0.372–0.797) were not significantly different, taking CRP as the reference (Delong Test P = 0.8378).

In Whites, the AUCs of D-dimer (Area = 0.630, 95% CI 0.543–0.718), CRP (Area = 0.633, 95% CI 0.547–0.719), and Ferritin (Area = 0.599, 95% CI 0.507–0.691), taking ferritin as the reference, were not statistically different (DeLong Test P = 0.8265).

In Hispanics, the AUC of D-dimer (Area = 0.742, 95% CI 0.651–0.833), CRP (Area = 0.567, 95% CI 0.472–0.662), and Ferritin (Area = 0.575, 95% CI 0.447–0.673), taking CRP as the reference, were significant (Delong Test P = 0.0038). Compared to CRP, D-dimer was statistically different (DeLong Test P = 0.0041).

In Asian/Indian, the difference between AUC of D-dimer (Area = 0.659, 95% CI 0.487–0.830), CRP (Area = 0.746, 95% CI 0.571–0.920), and Ferritin (Area = 0.709, 95% CI 0.524–0.894) is not significant (P = 0.7076). Compared to D-dimer, none of the ROC were statistically different (DeLong Test P = 0.8014).

In Blacks, univariate Cox regression revealed that increased risk of in-hospital mortality (poor in-hospital survival) was significantly associated with older age ≥ 60 years vs < 60 years (HR = 4.11, 95% CI 1.31–12.90; P = 0.0152), mechanical ventilation compared vs non-mechanical ventilation (HR = 4.36, 95% CI 1.29–14.76, P = 0.0180), qSOFA score of 2 vs 0 (HR = 2.61, 95% CI 1.16–5.89; P = 0.029), qSOFA score of 3 vs 0 (HR = 21.33, 95% CI 8.84–51.44; P < 0.0001) and creatinine ≥ 1.5 mg/dL (HR = 3.59, 95% CI 1.61–8.02, P = 0.0018) (Table 2, Additional file 1: Table S3). Univariate analysis indicated that MP dose did not violate the proportional hazard assumption NPH supremum test P = 0.7120, hence a multivariable proportional hazard Cox model was employed. As candidates in the full multivariate model, these variables along with P-values < 0.15 in the univariate analysis were entered into the subsequent multivariate’s backward elimination procedure. The extra values that met criterion for be entered in this multivariable stage included: MP dose, obesity status (BMI ≥ 30 kg/m2), respiratory rate > 22 bpm, and qSOFA. After the variable selection step of multivariable concluded, it was indicated that risk of poor in-hospital survival was significantly associated with older age ≥ 60 years (HR = 3.71, 95% CI 1.15–11.95, P = 0.0281), mechanical ventilation compared (HR = 5.07, 95% CI 1.42–18.07, P = 0.0123) and creatinine ≥ 1.5 mg/dL (HR = 3.61, 95% CI 1.71–7.59, P = 0.0007).

Table 2 Univariable and multivariate analysis of in-hospital survival in admitted COVID-19 Black patients (n = 102)

In Whites, univariate Cox regression indicated that significant benefit in in-hospital survival was associated with LDMP vs NMP (HR = 0.35, 95% CI 0.22–0.56, P < 0.0001) and HDMP vs NMP (HR = 0.55, 95% CI 0.34–0.88, P = 0.0132). (Table 3, Additional file 1: Table S4) The analysis also revealed that increased risk of in-hospital mortality (poor in-hospital survival) was significantly associated with hypertension (HR = 2.03, 95% CI 1.36–3.05, P = 0.0006), respiratory rate > 22 bpm (HR = 1.66, 95% CI 1.07–2.57, P = 0.0232), chronic kidney disease (HR = 1.66, 95% CI 1.01–2.73, P = 0.0468), older age ≥ 60 years (HR = 2.22, 95% CI 1.31–3.75; P = 0.0030), cancer (HR = 1.90, 95% CI 1.24–2.92, P = 0.0032), receiving mechanical ventilation compared to no supplemental oxygen (HR = 1.72, 95% CI 1.09–2.71, P = 0.0194), qSOFA score of 1 vs 0 (HR = 2.28, 95% CI 1.53–3.38; P < 0.0001), qSOFA score of 2 vs 0 (HR = 6.83, 95% CI 4.13–11.28; P < 0.0001), qSOFA score of 3 vs 0 (HR = 8.21, 95% CI 5.64–11.95; P < 0.0001) and creatinine ≥ 1.5 mg/dL (HR = 2.02, 95% CI 1.39–2.94, P = 0.0002). In the White cohort, univariate analysis indicated that MP dose did not satisfy the proportional hazard assumption NPH supremum test P < 0.0001, hence a multivariable non-proportional hazard Cox model was employed. As candidates in the full multivariate model, these variables along with P-values < 0.15 in the univariate analysis were entered into the subsequent multivariate’s backward elimination procedure. The extra values that met criterion for be entered in this multivariable stage included sex. After the variable selection step of multivariable concluded, it was reported that significant benefit in in-hospital survival was associated with LD MP vs NMP (HR = 0.37, 95% CI 0.23–0.59, P<0.0001) and HD MP vs NMP (HR = 0.56, 95% CI 0.35–0.91, P = 0.0183). In addition, the risk of poor in-hospital survival was significantly associated with cancer (HR = 1.68, 95% CI 1.03–2.69, P = 0.0376), qSOFA score of 1 vs 0 (HR = 1.81, 95% CI 1.09–2.98, P = 0.0213), qSOFA score of 2 vs 0 (HR = 5.16, 95% CI 3.07–8.66; P < 0.0001), qSOFA score of 3 vs 0 (HR = 11.81, 95% CI 6.46–21.57; P < 0.0001) and creatinine ≥ 1.5 mg/dL (HR = 2.16, 95% CI 1.39–3.34, P = 0.0006).

Table 3 Univariable and multivariable analysis of in-hospital survival in admitted COVID-19 white patients (n = 338)

In Hispanics, univariate Cox regression indicated that significant benefit in in-hospital survival was associated with LD MP vs NMP (HR = 0.35, 95% CI 0.18–0.67, P = 0.0015) (Table 4, Additional file 1: Table S5). In addition, the analysis indicated that increased risk of in-hospital mortality (poor in-hospital survival) was significantly associated with HD MP vs LD MP (HR = 2.08, 95% CI 1.20–3.54, P = 0.0073), hypertension (HR = 2.65, 95% CI 1.66–4.24, P < 0.0001), chronic kidney disease (HR = 2.65, 95% CI 1.38–5.09, P = 0.0034), older age ≥ 60 years (HR = 2.56, 95% CI 1.55–4.23; P = 0.0002), cancer (HR = 2.51, 95% CI 1.19–5.27, P = 0.0154), D-dimer ≥ 2 mcg/mL (HR = 1.79, 95% CI 1.05–3.04, P = 0.0320), qSOFA score of 2 vs 0 (HR = 1.99, 95% CI 1.11–3.55; P = 0.0202), and creatinine ≥ 1.5 mg/dL (HR = 2.16, 95% CI 1.20–3.89, P = 0.0100). The univariate analysis indicated that MP dose did not violate the proportional hazard assumption NPH supremum test P = 0.1120. Since creatinine > 1.5 mg/dL, which did not satisfy the proportional hazards assumption (supremum test P = 0.0190), was selected as one of the variables to initiate the multivariate model fitting process, a multivariable Non-proportional hazard Cox model was employed. As candidates in the full multivariate model, these variables were entered into the subsequent multivariate’s backward elimination procedure. No extra variables reported P < 0.15 in the univariate stage. At the conclusion of variable selection step of multivariable analysis, it was indicated that risk of poor in-hospital survival was significantly associated with hypertension (HR = 2.52, 95% CI 1.14–4.47, P = 0.0012), cancer (HR = 3.04, 95% CI 1.16–8.02, P = 0.0244) and D-dimer ≥ 2 mcg/mL (HR = 2.21, 95% CI 1.23–3.96, P = 0.0077).

Table 4 Univariable and multivariate analysis of in-hospital survival in admitted COVID-19 Hispanic patients (n = 251)

In Asian/Indians, univariate Cox regression indicated that significant benefit in in-hospital survival was associated with LDMP vs NMP (HR = 0.16, 95% CI 0.04–0.62, P = 0.0080) and HDMP vs NMP (HR = 0.24, 95% CI 0.06–0.92, P = 0.0370). (Table 5, Additional file 1: Table S6) In addition, the analysis indicated that increased risk of in-hospital mortality (poor in-hospital survival) was significantly associated with chronic kidney disease (HR = 4.77, 95% CI 1.96–11.62, P = 0.0006), older age ≥ 60 years (HR = 9.45, 95% CI 1.22–72.98; P = 0.0313), and CRP > 20 mg/L (HR = 4.46, 95% CI 1.64–12.10, P = 0.0033). In the Asian/Indian cohort, univariate analysis indicated that MP dose did not satisfy the proportional hazard assumption NPH test P = 0.0320, hence a multivariable non-proportional hazard Cox model was employed. As candidates in the full multivariate model, these variables along with variables reporting P-values < 0.15 in the univariate analysis were entered into the subsequent multivariate’s backward elimination procedure. The extra values that met criterion for being entered in this multivariable stage included: hypertension, Ferritin ≥ 1400 μg/L and creatinine ≥ 1.5 mg/dL. At the conclusion of variable selection step of multivariable analysis, it was indicated that significant benefit in in-hospital survival was associated with LD MP vs NMP (HR = 0.13, 95% CI 0.03–0.61, P = 0.0101) and HD MP vs NMP (HR = 0.15, 95% CI 0.04–0.64, P = 0.0106). In addition, risk of poor in-hospital survival was significantly associated with chronic kidney disease (HR = 6.36, 95% CI 1.87–21.67, P = 0.0031) and CRP > 20 mg/L (HR = 5.02, 95% CI 1.72–14.66, P = 0.0032).

Table 5 Univariable analysis of in-hospital survival in admitted COVID-19 Asian/Indian patients (n = 61)

Overall, in hospital survival was significantly difference between the racial/ethnic groups (Wilcoxon P = 0.0320). This result was driven by the significant difference between Whites and Blacks (P = 0.0249) and Whites and Asian/Indians (P = 0.0463), after adjusting for multiple testing. The 30-day in hospital survival rates amongst Asia/Indians, Blacks, Hispanics, and Whites were 54.4% (95% CI 35.6–72.5%), 37.9% (95% CI 19.8–57.8%), 33.4% (95% CI 22.9–44.7%), and 41.0% (95% CI 32.6–49.6%), respectively. These differences become less significant after 30 days. There were 23 (7.42%) nosocomial infections in the no methylprednisolone group and 38 (13.72%) nosocomial infections in the methylprednisolone group (P 0.0145). The main cause of death was ARDS. (Additional file 1: Fig. S3).

Due to the varying frequencies, we were not able to compare the association of NMP, LDMP, and HDMP between racial/ethnic groups. We did compare NMP, LDMP, and HDMP and their association with survival between Whites and Hispanics and they were not significant. (P > 0.05).

(Additional file 1: Figures S4, S5, S6) Discussion

Social determinants, through weathering and/or allostatic load theory, can lead to low-grade chronic inflammation and persistently elevated inflammatory markers at baseline [25, 26]. With lack of access to quality health care, co-morbidities such as diabetes, CKD, hypertension, and coronary artery disease can remain undiagnosed in these minority groups. This may explain why in our study some inflammatory makers were above the cut-off values for mortality even in survivors [median CRP in Black (10.2 mg/mL) and Hispanic (13.70 mg/mL)].

However, there may also be genetic influences to elevated inflammatory markers. F3 and sickle cell variant (HBB rs334) are associated with higher d-dimer levels in Blacks and thalassemia, iron overload, or high iron (HFE) mutations are associated with elevated ferritin levels in Asian/Indians [27, 28]. This can explain the elevated median d-dimer in Black survivors (1.91 pg/mL) and median Ferritin in Asian/Indian survivors (> 1265 μg/L).

In the multivariate Cox regression, D-dimer, hypertension and cancer were independently associated with higher mortality in Hispanics. D-dimer, a fragment of fibrin and marker of coagulation, has also been implicated in angiogenesis, tumor cell invasion, and metastatic spread in cancer [29, 30]. Three studies with predominantly Hispanic population on COVID-19 showed an association with d-dimer and mortality [29, 31, 32]. One study on Hispanics showed that d-dimer and hypertension were associated with worse survival [32]. A study on COVID-19 showed that d-dimer was higher in patients with hypertension compared to no hypertension [29]. Non-COVID-19 studies showed that D-dimer was associated with hypertension in Blacks [33, 34]. This was not seen in our study and maybe due to smaller sample.

In Blacks and Whites, initial creatinine ≥ 1.5 mg/dL and higher qSOFA were independently associated with higher mortality. Prehospitalization creatinine was not available in all patients to confirm a diagnosis of acute kidney injury (AKI) [35]. Mortality can reach up to 50% in hospitalized patients with COVID-19 who developed AKI and COVID-19 patients were more likely to develop AKI than those without COVID-19 [36, 37]. Causes of AKI are glomerular injury from cytopathogenic effect of SARS-COV-2 and tubular injury from right heart failure, hypovolemia, endothelial dysfunction, hypercoagulability, and mechanical ventilation settings [38,39,40,41,42].

CKD was independently associated with mortality in Asian/Indians. Other studies have suggested that CKD has the highest risk for mortality out of all comorbidities, and has been associated with hypercoagulopathy and elevated inflammatory markers, particularly CRP [36, 43,44,45,46].

In the univariate analysis, age ≥ 60 years was associated with worse survival in all racial/ethnic groups. However, after multivariate Cox regression, age was noted to be independently associated with mortality in Blacks only. Mechanical ventilation has been shown to have higher mortality in COVID-19 [21] and was also associated with increased mortality in Blacks only.  In this study, comorbidities may have more prognostic significance rather than age in the other racial/ethnic groups.

Chronic inflammatory states are hypothesized to cause steroid desensitization, requiring higher doses of methylprednisolone to mount an effective response. LDMP and HDMP was associated with prolonged in hospital survival in Whites and Asian/Indians. There was no additional benefit of HDMP over LDMP. Lack of benefit of HDMP maybe due to a dose dependent increased severity of critical illness polyneuropathy, increased incidence of secondary infections or practice of higher doses for patients who are very sick [47]. Methylprednisolone, regardless of dose, was not associated with improved survival in Hispanics and Blacks. This may be due to vitamin D deficiency that is common in Blacks and Hispanics [48,49,50,51,52,53]. The National Health and Nutrition Examination Survey suggested that prevalence of Vitamin D deficiency is 65.4% in Blacks, 28.9% in Hispanics, and 14% for Whites in the United States [52]. Skin pigmentation is an evolutionary adaptation to intensive solar ultraviolet radiation. Ultraviolet B (UVB) is needed to convert 7-dehydrocholesterol to Vitamin D3. The abundant melanin which contributes to brown-black pigmentation, absorbs UVR and protects from skin damage and ensures adequate Vitamin D production in low latitudes [54]. The protective benefits of abundant melanin in higher latitudes, such as in the United States, are offset by decreased production of Vitamin D [50,51,52,53,54]. Therefore, patients with darker skin pigmentation living in higher latitudes like United States have lower vitamin D levels compared to patients with darker skin pigmentation living in lower latitudes like Africa [54]. They would require higher amounts of ultraviolet radiation (UVR) to make the same amount of Vitamin D as those with lighter pigmentation in higher latitudes [55,56,57,58,59,60]. Vitamin D sufficiency is often associated with efficacy of steroid response [6]. This may be due to Vitamin D’s upregulation of steroid receptor GR-α [61] or downregulation of IL-23R driven glucorticoid resistant MDR1+ proinflammatory Th17 cells [62]. Non-steroidal anti-inflammatories such as tocilizumab maybe indicated sooner for these racial/ethnic minorities [63]. Baricitinib is promising although the studies had proportionally fewer racial/ethnic minorities [64, 65].

Strengths of this study include the median days in both methylprednisolone and no methylprednisolone were 5 days, which is at the beginning stages of inflammatory phase of the disease and there were minimal number of patients on remdesivir, which has become standard of care for COVID-19. Therefore, focus there has been on the inflammatory phase and on anti-inflammatory medications.

Limitations

Our study has several limitations. First, since it is an observational study and there may be known and unknown confounders. However, propensity matching was employed to limit the known confounders. Second, misclassification of data is possible due to manual extraction of structured and unstructured data from medical health records. Third, there was a higher prevalence of Whites and non-white/non-black Hispanics, which might have skewed the analysis. Fourth, methylprednisolone was used as a rescue, given to patients who were at a higher risk of death. During the initial pandemic surge, there were reservations on the use of methylprednisolone due to extrapolated data on prolonging viral shedding in SARS and MERS and worse mortality in Influenza. Therefore, it was used as a rescue and reserved for patients who are already on high oxygen supplementation requirements or on mechanical ventilation.

Despite the possibility of corticosteroid resistance in certain racial/ethnic groups, the effect of corticosteroids in patients on lesser amounts of oxygen supplementation such as nasal cannula was not available.

Conclusions

Racial and ethnic disparities in inflammatory markers preclude the use of one marker as a solitary measure of mortality. Low dose and high dose methylprednisolone were associated with prolonged survival in Whites and Asian/Indians. However, high dose was not superior to low dose to prolonging survival. Methylprednisolone, regardless of dose, was not associated with prolonged survival in Blacks and Hispanics. Large, prospective studies are needed to confirm these conclusions.

Availability of data and materials

The datasets are available from the corresponding author upon reasonable request.

Abbreviations

ABG:

Arterial blood gas

AFOP:

Acute fibrinous organizing pneumonia

AKI:

Acute kidney injury

ARDS:

Acute respiratory distress syndrome

AUROC:

Area under receiver operating characteristic

BMI:

Body mass index

BNP:

Beta natriuretic peptide

CBC:

Complete blood count

CI:

Confidence interval

CMP:

Complete metabolic profile

Covid-19:

Coronavirus disease 2019

CRP:

C-reactive protein

EHR:

Electronic health record

FiO2 :

Fraction of inspired oxygen

GR-α :

Glucocorticoid receptor alpha

HCQ:

Hydroxychloroquine

HMH:

Hackensack Meridian Health

HR:

Hazard ratio

his:

In hospital survival

IL-6:

Interleukin 6

IL-23R:

Interleukin 23R

IQR:

Interquartile range

HDMP:

High dose methylprednisolone

HR:

Hazard ratio

LDMP:

Low dose methylprednisolone

MDR1:

Multidrug resistant mutation 1

MV:

Mechanical ventilation

NMV:

Non-mechanical ventilation

NMP:

No methylprednisolone

OP:

Organizing pneumonia

PaO2 :

Partial pressure of arterial oxygen

PH:

Proportional hazard

qSOFA:

Quick sequential failure assessment

RR:

Relative risk

SARS-COV-2:

Severe acute respiratory syndrome coronavirus 2

SpO2 :

Oxygen saturation

Th17:

T helper 17 cell

UVR:

Ultraviolet radiation

UVB:

Ultraviolet B

References

  1. WHO Coronavirus (COVID-19) Dashboard. World Health Organization. https://covid19.who.int/. Updated July 2, 2021. Accessed 4 July 2021.

  2. Wu Z, McGoogan JM. Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) outbreak in China: summary of a report of 72 314 cases from the Chinese Center for Disease Control and Prevention. JAMA. 2020;323(13):1239–42.

    CAS  PubMed  Google Scholar 

  3. Huang C, Wang Y, Li X, et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet. 2020;395:497–506.

    CAS  PubMed  PubMed Central  Google Scholar 

  4. Price-Haywood EG, Burton J, Fort D, Seoane L. Hospitalization and mortality among black patients and white patients with COVID-19. NEJM. 2020;382:2534–43.

    CAS  PubMed  Google Scholar 

  5. COVID-19 Racial and Ethnic Health Disparities. Centers for Disease Control and Prevention. Updated December 10, 2020. https://www.cdc.gov/coronavirus/2019-ncov/community/health-equity/racial-ethnic-disparities/index.html. Accessed 14 Apr 2021.

  6. Karan A, Ali K, Rambaran K, Del Corral P, Sakhamuri S, Teelucksingh S. COVID-19 and ethnicity: does reduced responsiveness to glucocorticoids explain the more aggressive nature of disease among minorities? Med Hypotheses. 2020. https://doi.org/10.1016/j/mehy.2020.110220.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Harrison E. Ethnicity and outcomes from COVID-19 in UK hospital patients using the ISARIC WHO Clinical Characterisation Protocol: prospective observational cohort study. (2020). https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/895695/S0452_CO-CIN_Ethnicity_and_Outcomes_from_COVID-19_in_UK_Hospital_Patients.pdf.

  8. Ribeiro H, Lima VM, Waldman EA. In the COVID-19 pandemic in Brazil, do brown lives matter? Lancet Glob Health. 2020;8(8):e976–7.

    PubMed  PubMed Central  Google Scholar 

  9. Leisman DE, Ronner L, Pinotti R, et al. Cytokine elevation in severe and critical COVID-19: a rapid systematic review, meta-analysis, and comparison with other inflammatory syndromes. Lancet Respir Med. 2020;8:1233–44.

    CAS  PubMed  PubMed Central  Google Scholar 

  10. Manson JJ, Crooks C, Naja M, et al. COVID-19 associated hyperinflammation and escalation of patient care: a retrospective longitudinal cohort study. Lancet Rheumatol. 2020;10:E594–602.

    Google Scholar 

  11. Zhang L, Yan X, Fan Q, Liu H, Liu X, Liu Z, Zhang Z. D-dimer levels on admission to predict in-hospital mortality in patients with Covid-19. J Thromb Haemost. 2020;18(6):1324–9.

    CAS  PubMed  Google Scholar 

  12. Gómez-Pastora J, Weigand M, Kim J, et al. Hyperferritinemia in critically ill COVID-19 patients—is ferritin the product of inflammation or a pathogenic mediator? Clin Chim Acta. 2020;509:249–51.

    PubMed  PubMed Central  Google Scholar 

  13. Sabaka P, Koščálová A, Straka I, et al. Role of interleukin 6 as a predictive factor for a severe course of Covid-19: retrospective data analysis of patients from a long-term care facility during Covid-19 outbreak. BMC Infect Dis. 2021;21:308. https://doi.org/10.1186/s12879-021-05945-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Rose-John S, Winthrop K, Calabrese L. The role of IL-6 in host defense against infections: immunobiology and clinical implications. Nat Rev Rheumatol. 2017;13(7):399–409.

    CAS  PubMed  Google Scholar 

  15. Santa Cruz A, Mendes-Frias A, Oliveira AI, et al. Interleukin-6 is a biomarker for the development of fatal severe acute respiratory syndrome coronavirus 2 pneumonia. Front Immunol. 2021. https://doi.org/10.3389/fimmu.2021.613422.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Ali N. Elevated level of C-reactive protein may be an early marker to predict risk for severity of COVID-19. J Med Virol. 2020;92(11):2409–11. https://doi.org/10.1002/jmv.26097.

    Article  CAS  PubMed  Google Scholar 

  17. Taylor DR, Buchanan D, Al-Hasani W, et al. Biochemical abnormalities in COVID-19: a comparison of white versus ethnic minority populations in the UK. J Clin Pathol. 2021. https://doi.org/10.1136/jclinpath-2021-207446.

    Article  PubMed  Google Scholar 

  18. Ackermann M, Verleden SE, Kuehnel M, et al. Pulmonary vascular endothelialitis, thrombosis, and angiogenesis in COVID-19. NEJM. 2020;383:120–8.

    CAS  PubMed  Google Scholar 

  19. Hariri LP, North CM, Shih AR, et al. Lung histopathology in coronavirus disease 2018 as compared with severe respiratory syndrome and H1N1 influenza. Chest. 2021;159(1):73–85.

    CAS  PubMed  Google Scholar 

  20. Recovery Group. Dexamethasone in hospitalized patients with COVID-19. NEJM. 2021;384:693–704.

  21. Go RC, Shah R, Nyirenda T, et al. Methylprednisolone and 60 days in hospital survival in coronavirus disease 2019 pneumonia. Crit Care Explor. 2019. https://doi.org/10.1097/CCE.0000000000000493.

    Article  Google Scholar 

  22. Austin PC. The use of propensity score methods with survival or time-to-event outcomes: reporting measures of effect similar to those used in randomized experiments. Stat Med. 2014;33(7):1242–58.

    PubMed  Google Scholar 

  23. Lin DY, Wei LJ, Ying Z. Checking the Cox Model with cumulative sums of martingale-based residuals. Biometrika. 1993;80(3):557–72.

    Google Scholar 

  24. Allison PD. Survival analysis using SAS®: a practical guide. 2nd ed. Cary: SAS Institute Inc.; 1995.

  25. Schulz AJ, Mentz G, Lachance L, et al. Associations between socioeconomic status and allostatic load: effects of neighborhood poverty and tests of mediating pathways. Am J Public Health. 2012;102(9):1706–14.

    PubMed  PubMed Central  Google Scholar 

  26. Schmeer KK, Tarrence J. Racial–ethnic disparities in inflammation: evidence of weathering in childhood? J Health Soc Behav. 2018;59(3):411–28.

    PubMed  PubMed Central  Google Scholar 

  27. Rastogi D, Jung M, Strizich G, et al. Association of systemic inflammation, adiposity, and metabolic dysregulation with asthma burden among Hispanic adults. Respir Med. 2017;125:72–81.

    PubMed  PubMed Central  Google Scholar 

  28. Raffield LM, Zakai NA, Duan Q, et al. D-dimer in African Americans: whole genome sequence analysis and relationship to cardiovascular disease risk in the Jackson Heart Study. Arterioscler Thromb Vasc Biol. 2017;37:2220–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. Andina T, Tampubolon D. 54. Comparation of D-dimer levels in hypertensive and non-hypertensive adult COVID-19 patients. J Hypertens. 2021;39:e14. https://doi.org/10.1097/01.hjh.0000752564.68879.af.

    Article  Google Scholar 

  30. Simes J, Robledo KP, White HD, et al. D-dimer predicts long-term cause-specific mortality, cardiovascular events, and cancer in patients with stable coronary heart disease: LIPID Study. Circulation. 2018;138(7):712–23.

    CAS  PubMed  Google Scholar 

  31. Song X, Ji J, Reva B, et al. Post-anticoagulant D-dimer is a highly prognostic biomarker of COVID-19 mortality. ERJ Open Res. 2021;7:00018–2021.

    PubMed  PubMed Central  Google Scholar 

  32. Chilimuri S, Sun H, Alemam A, Mantri N, Shehi E, Tejada J, Yugay A, Nayudu SK. Predictors of mortality in adults admitted with COVID-19: retrospective cohort study from New York City. West J Emerg Med. 2020;21(4):779–84.

    PubMed  PubMed Central  Google Scholar 

  33. Khaleghi M, Saleem U, McBane RD, Mosley TH Jr, Kullo IJ. African-American ethnicity is associated with higher plasma levels of D-dimer in adults with hypertension. J Thromb Haemost. 2009;7(1):34–40.

    CAS  PubMed  Google Scholar 

  34. Mukaz DK. D-dimer and the risk of incident hypertension in black and white Americans. J Stroke Cardiovasc Dis. 2020;29:5.

    Google Scholar 

  35. Summary of recommendation statements. Kidney Int Suppl (2011). 2012;2(1):8–12.

  36. Chan L, Chaudhary K, Saha A, et al. AKI in hospitalized patients with COVID-19. J Am Soc Nephrol. 2021;32(1):151–60.

    CAS  PubMed  Google Scholar 

  37. Fisher M, Neugarten J, Bellin E, Yunes M, Stahl L, Johns TS, Abramowitz MK, Levy R, Kumar N, Mokrzycki MH, Coco M, Dominguez M, Prudhvi K, Golestaneh L. AKI in hospitalized patients with and without COVID-19: a comparison study. J Am Soc Nephrol. 2020;31(9):2145–57.

    CAS  PubMed  PubMed Central  Google Scholar 

  38. Pan X-W, Da X, Zhang H, Zhou W, Wang L-H, Cui X-G. Identification of a potential mechanism of acute kidney injury during the COVID-19 out-break: a study based on single-cell transcriptome analysis. Intensive Care Med. 2020. https://doi.org/10.1007/s00134-020-06026-1.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Ronco C, Reis T, Husain-Syed F. Management of acute kidney injury in patients with COVID-19. Lancet Respir Med. 2020;8:738–42.

    CAS  PubMed  PubMed Central  Google Scholar 

  40. Ackermann M, Verleden SE, Kuehnel M, Haverich A, Welte T, Laenger F, et al. Pulmonary vascular endothelialitis, thrombosis, and angiogenesis in Covid-19. N Engl J Med. 2020;383:120–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  41. Su H, Yang M, Wan C, Yi L-X, Tang F, Zhu H-Y, et al. Renal histopathological analysis of 26 postmortem findings of patients with COVID-19 in China. Kidney Int. 2020. https://doi.org/10.1016/j.kint.2020.04.003.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Geri G, Darmon M, Zafrani L, et al. Acute kidney injury in SARS-CoV2-related pneumonia ICU patients: a retrospective multicenter study. Ann Intensive Care. 2021;11:86. https://doi.org/10.1186/s13613-021-00875-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Williamson EJ, Walker AJ, Bhaskaran K, et al. Factors associated with COVID-19-related death using OpenSAFELY. Nature. 2020;584:430–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  44. Chen TK, Knicely DH, Grams ME. Chronic kidney disease diagnosis and management: a review. JAMA. 2019;322(13):1294–304.

    CAS  PubMed  PubMed Central  Google Scholar 

  45. Muslimovic A, Rasic S, Tulumovic D, Hasanspahic S, Rebic D. Inflammatory markers and procoagulants in chronic renal disease stages 1–4. Med Arch. 2015;69(5):307–10.

    PubMed  PubMed Central  Google Scholar 

  46. Annuk M, Soveri I, Zilmer M, Lind L, Hulthe J, Fellström B. Endothelial function, CRP and oxidative stress in chronic kidney disease. J Nephrol. 2005;18(6):721–6.

    CAS  PubMed  Google Scholar 

  47. The National Heart, Lung, and Blood Institute Acute Respiratory Distress Syndrome (ARDS) Clinical Trials Network. Efficacy and safety of corticosteroids for persistent acute respiratory distress syndrome. NEJM. 2006;354:1671–84.

  48. Meltzer DO, Best TJ, Zhang H, Vokes T, Arora VM, Solway J. Association of vitamin D levels, race/ethnicity, and clinical characteristics with COVID-19 test results. JAMA NetwOpen. 2021;4(3): e214117. https://doi.org/10.1001/jamanetworkopen.2021.4117.

    Article  Google Scholar 

  49. Libon F, Cavalier E, Nikkels A. F: skin color is relevant to vitamin D synthesis. Dermatology. 2013;227:250–4.

    CAS  PubMed  Google Scholar 

  50. Olerod G, Hulten LM, Hammersten O, Klingber G. The variation in free 25-hydroxy vitamin D and vitamin D binding protein with season and vitamin D status. Endocr Connect. 2017;6(2):111–20.

    PubMed  PubMed Central  Google Scholar 

  51. Güven M, Gültekin H. The effect of high-dose parenteral vitamin D3 on COVID-19-related inhospital mortality in critical COVID-19 patients during intensive care unit admission: an observational cohort study. Eur J Clin Nutr. 2021;75:1383–8.

    PubMed  Google Scholar 

  52. Taksler GB, Cutler DM, Giovannucci E, Keating NL. Vitamin D deficiency in minority populations. Public Health Nutr. 2015;18(3):379–91.

    PubMed  Google Scholar 

  53. Meltzer DO, Best TJ, Zhang H, Vokes T, Arora V, Solway J. Association of vitamin D status and other clinical characteristics with COVID-19 test results. JAMA Netw Open. 2020;3(9):e2019722.

    PubMed  PubMed Central  Google Scholar 

  54. Ames BN, Grant WB, Willett WC. Does the high prevalence of vitamin D deficiency in African Americans contribute to health disparities? Nutrients. 2021;13:499. https://doi.org/10.3390/nu13020.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Luxwolda MF, Kuipers RS, Kema IP, Dijck-Brouwer DA, Muskiet FA. Traditionally living populations in East Africa have a mean serum 25-hydroxyvitamin D concentration of 115 nmol/l. Br J Nutr. 2012;108:1557–61.

    CAS  PubMed  Google Scholar 

  56. Luxwolda MF, Kuipers RS, Kema IP, van der Veer E, Dijck-Brouwer DA, Muskiet FA. Vitamin D status indicators in indigenous populations in East Africa. Eur J Nutr. 2013;52:1115–25.

    CAS  PubMed  Google Scholar 

  57. Chiang D, Kramer H, Luke A, Cooper R, Aloia J, Bovet P, Plange-Rhule J, Forrester T, Lambert V, Camacho P, et al. 25-Hydroxyvitamin D and blood pressure: a plateau effect in adults with African ancestry living at different latitudes. J Hypertens. 2017;35:968–74.

    CAS  PubMed  Google Scholar 

  58. Ginde AA, Liu MC, Camargo CA Jr. Demographic differences and trends of vitamin D insufficiency in the US population, 1988–2004. Arch Intern Med. 2009;169:626–32.

    CAS  PubMed  PubMed Central  Google Scholar 

  59. Enechukwu N, Cockburn M, Ogun G, Ezejiofor OI, George A, Ogunbiyi A. Higher vitamin D levels in Nigerian albinos compared with pigmented controls. Int J Dermatol. 2019;58:1148–52.

    CAS  PubMed  Google Scholar 

  60. Webb AR, Kazantzidis A, Kift RC, Farrar MD, Wilkinson J, Rhodes LE. Colour counts: sunlight and skin type as drivers of vitamin D deficiency at UK latitudes. Nutrients. 2018;10(4):457.

    PubMed Central  Google Scholar 

  61. Mahboub B, Al Heialy S, Hachim MY, et al. Vitamin D regulates the expression of glucocorticoid receptors in blood of severe asthmatic patients. J Immunol Res. 2021;2021:9947370. https://doi.org/10.1155/2021/9947370.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Mann EH, Ho TR, Pfeffer PE, et al. Vitamin D counteracts an IL-23-dependent IL-17A+IFN-γ+ response driven by urban particulate matter. Am J Respir Cell Mol Biol. 2017;57(3):355–66.

    CAS  PubMed  PubMed Central  Google Scholar 

  63. Salama C, Han J, Yau L, et al. Tocilizumab in patients hospitalized with COVID-19 pneumonia. NEJM. 2021;384:20–30.

    CAS  PubMed  Google Scholar 

  64. Marconi VC, Ramanan AV, de Bono S, et al. Efficacy and safety of baricitinib for the treatment of hospitalized adults with COVID-19 (COV-BARRIER): a randomized double-blind, parallel-group, placebo-controlled phase 3 trial. Lancet Respir Med. 2021. https://doi.org/10.1016/S2213-2600(21)00331-3.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Kalil AC, Patterson TF, Mehta AK, et al. Baricitinib plus remdesivir for hospitalized adults with COVID-19. NEJM. 2021;384:795–807.

    CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Dr. Andrew Ip and Dr. Stuart L. Goldberg.

Funding

There was no funding. The corresponding author had final responsibility for decision to submit for publication.

Author information

Authors and Affiliations

Authors

Contributions

RCG, TN, EGP, ACG, ZH, and SJS were responsible for study design and planning. RCG, MB, DKH, KK, and MR contributed to data collection. RCG, TN, EGP, AG, ACG, ZH, and SJS contributed to data analysis and interpretation. All authors contributed to writing the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Ronaldo C. Go.

Ethics declarations

Ethics approval and consent to participate

Ethical approval and clearance were obtained from IRB board from participating hospitals. Approval was obtained by the Hackensack Meridian Health Institutional Review Board (#Pro2020-0485) and the study was also registered on ClinicalTrials.Gov as a prospective observational database (NCT04347993). Waiver of consent [45cfr46.116(d)] was approved by Hackensack Meridian Health Institutional Review Board and was applied due to minimal risk, waiver will not adversely affect the rights and welfare of subjects, the research could not be carried out without the waiver, and whenever possible the subjects were provided with additional pertinent information after participation. At that time, the risk is reasonable due to importance of knowledge gained. The study was conducted according to the International Conference on Harmonization (ICH), Good Clinical Practice (GCP), the Declaration of Helsinki, Institutional Review Boards (IRB) and in accordance with the U.S. Code of Federal Regulations on Protection of Human Rights (21 CFR 50).

Consent for publication

Not applicable.

Competing interests

Dr. Ronaldo C. Go participated in a long-term follow-up committee for long term follow-up in COVID-19 patients sponsored by Hoffmann-LaRoche. All the other authors have 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

Additional file 1: Table S1.

Baseline demographics disease characteristics of unmatched population and propensity matched population. Table S2. Standardized mean differences (methylprednisolone − no methylprednisolone). Table S3. Non-proportional hazard assumption of covariates for In-hospital mortality in COVID-19 black patients. Table S4. Non-proportional hazard assumption of covariates for In-hospital mortality in COVID-19 White patients. Table S5. Non-proportional hazard assumption of covariates for In-hospital mortality in COVID-19 Hispanic patients. Table S6. Non-proportional hazard assumption of covariates for In-hospital mortality in COVID-19 Asian/Indian patients. Figure S1. Plot of differences methylprednisolone − no methylprednisolone in hospitalized COVID-19 patients. Figure S2. LPS cloud plots showing distributions of logit of propensity scores for MP and NMP treated COVID-19 patients. Figure S3. Kaplan Meier plot for overall in-hospital survival (IhS) for COVID-19. Figure S4. Kaplan Meier plot of In hospital survival of Whites versus Hispanics who received no methylprednisolone. Figure S5. Kaplan Meier plot of In hospital survival of Whites versus Hispanics who received low dose methylprednisolone. Figure S6. Kaplan Meier plot of In Hospital survival of Whites versus Hispanics who received high dose methylprednisolone.

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/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Go, R.C., Nyirenda, T., Bojarian, M. et al. Racial/ethnic disparities on inflammation and response to methylprednisolone in severe COVID-19 pneumonia. BMC Infect Dis 22, 254 (2022). https://doi.org/10.1186/s12879-022-07237-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12879-022-07237-1

Keywords