Abstract
The use of erythropoietin-stimulating agents (ESAs) as adjunctive therapy in critically ill patients with COVID-19 may have a potential benefit. This study aims to evaluate the effect of ESAs on the clinical outcomes of critically ill COVID-19 patients. A multicenter, retrospective cohort study was conducted from 01-03-2020 to 31-07-2021. We included adult patients who were ≥ 18 years old with a confirmed diagnosis of COVID-19 infection and admitted to intensive care units (ICUs). Patients were categorized depending on ESAs administration during their ICU stay. The primary endpoint was the length of stay; other endpoints were considered secondary. After propensity score matching (1:3), the overall included patients were 120. Among those, 30 patients received ESAs. A longer duration of ICU and hospital stay was observed in the ESA group (beta coefficient: 0.64; 95% CI: 0.31-0.97; P = < .01, beta coefficient: 0.41; 95% CI: 0.12-0.69; P = < .01, respectively). In addition, the ESA group's ventilator-free days (VFDs) were significantly shorter than the control group. Moreover, patients who received ESAs have higher odds of liver injury and infections during ICU stay than the control group. The use of ESAs in COVID-19 critically ill patients was associated with longer hospital and ICU stays, with no survival benefits but linked with lower VFDs.
Keywords
Introduction
The COVID-19 pandemic has resulted in approximately 626 million cases and 6.5 million fatalities worldwide. 1 Although the majority of patients with COVID-19 had minor symptoms,2,3 around 14% of the reported cases of COVID-19 cases progressed to serious illness and developed complications. 3 The severity of the illness may be attributed to the inflammatory response known as the “Cytokine Storm” caused by COVID-19. 4 One of the major respiratory complications of COVID-19 is acute respiratory distress syndrome (ARDS) which exacerbates other acute inflammatory disorders like sepsis, trauma, and pneumonia. 4 Patients with ARDS frequently develop respiratory failure, requiring invasive mechanical ventilation (MV) managed at the intensive care unit (ICU). 5
Thus, immunomodulator drugs that suppress cytokine storm and hyperinflammatory response could be a potential therapeutic choice for critically ill COVID-19. 6 Many immunomodulators have been studied for the treatment of hospitalized patients with COVID-19.7–10 These options include interleukin-6 inhibitors (IL-6), corticosteroids, and Janus kinase (JAK) inhibitors.7–10 In addition, the use of erythropoietin-stimulating agents (ESAs) for the treatment of severe anemia and reduced blood transfusion has been reported to have additional anti-inflammatory, immunomodulatory properties, and anti-apoptotic effects. 4 Therefore, ESA has been used in the treatment of acute lung injury induced by renal ischemia, ischemic stroke, and head trauma.11–17
A prospective observational study reported a lower amount of endogenous erythropoietin in mechanically ventilated critically ill patients with PaO2/FiO2 <300 which might be a risk factor for difficult extubation and poor response to anemia. 18 ESAs resistance, low level of circulating erythropoietin, and high inflammatory mediator, all can inhibit erythropoietic cellular proliferation and impairs tissue oxygenation, in critically ill patients who present with cytokines storm and hypoxia paradox.11,14,15,17 ESA can raise hemoglobin count and, thus, strengthens the oxygen reservoir, countering hypoxia. 4 Thus, ESAs have demonstrated a beneficial role in the management of ARDS.19,20 This effect was attained by ESAs pleiotropic effect, which protects lung integrity, improves respiration, restores homeostasis, antioxidant effect, and reduced inflammation.17,19–22
A case report by Hadadi et al. evaluated the use of ESAs in critically ill COVID-19 patients with severe anemia and reported improved anemia symptoms in addition to viral suppression. 23 Moreover, a case series by Begemann et al. reported that four critically ill patients with COVID-19 admitted to the ICU with respiratory failure who received ESA without blood transfusion showed improvement during hospital stay or rehabilitation. 11 On the other hand, the lack of funding has prevented the publication of the Ehrenreich et al. 2020 randomized controlled trial that examined the use of erythropoietin as an adjuvant therapy for the treatment of severe COVID-19. 4 The evidence proposed that in addition to anemia management, ESA can be used as supportive therapy in critically ill patients including patients with COVID-19. 4 The theory behind the use of erythropoietin in critically ill COVID-19 patients is the ability of ESAs to work in 3 different domains that control respiration through the brainstem, lung and phrenic nerve.11,12
The administration of ESAs in critically ill COVID-19 patients remains questionable due to insufficient evidence and the potential adverse effect.14,17,24–26 Recombinant human erythropoietin used as adjunctive therapy may have potential survival benefits in critically ill patients with COVID-19. 4 Therefore, the objective of this study is to evaluate the effect of ESAs on the clinical outcomes of critically ill COVID-19 patients.
Methods
Study Design
A multicenter retrospective cohort study was conducted from 01-03-2020 till 31-07-2021 and included all adult patients who were ≥ 18 years old with a confirmed diagnosis of COVID-19 infection and admitted to ICUs. COVID-19 diagnosis was confirmed by reverse transcriptase-polymerase chain reaction (RT-PCR) nasopharyngeal or throat swab testing. Critically ill patients were categorized into two sub-cohorts based on receiving ESA during ICU stay. There were no pre-defined criteria at the participating centers for ESA use; ESA is mainly used as either a continuation of ongoing therapy or a new initiation which is based on physician clinical judgment. ESA therapy was mainly administered via the subcutaneous route of administration. All patients were followed until they were discharged or died during the in-hospital stay. The King Abdullah International Medical Research Center (KAIMRC) - Institutional Review Board approved the study in January 2021 (Ref.# NRC21R-004-01) and approved by all participating centers’ internal IRB committees. All participating centers followed the same standards of care for case triaging, admission tiers, and transferring criteria. All methods were performed in accordance with relevant guidelines and regulations. Given the retrospective observational nature of the study, informed consent from the study patients was waived by the IRB committees.
Settings
This study was conducted at five hospitals in Saudi Arabia. The primary site for this study was King Abdulaziz Medical City (Riyadh), a tertiary care center, and the initial IRB approval originated from this center. Additional centers were added based on the availability of data, the center's willingness to participate, and geographic distribution, including King Abdulaziz Medical City (Jeddah), King Abdulaziz University Hospital (Jeddah), King Abdullah bin Abdulaziz University Hospital (KAAUH) (Riyadh), and King Salman Specialist Hospital (Hail).
Participants
The eligibility criteria included adult patients with confirmed COVID-19 infection through RT-PCR testing and exclusion of patients who received ESA therapy either prior to ICU admission without ICU continuation or after ICU discharge. In addition, patients were excluded if deceased within 24 h of ICU admission, had an ICU length of stay less than one day, were designated as do-not-resuscitate, or had unknown medical history. The flowchart of the study population is presented in Figure 1.

Flowchart for eligibility criteria.
Outcomes
This study aimed to assess the association between the administration of erythropoiesis-stimulating agents and clinical outcomes in critically ill COVID-19 patients.. The primary endpoint was the LOS, while secondary endpoints were ventilator-free days (VFDs), 30-day and in-hospital mortality, and complication(s) that occurred during ICU stay (ie, new onset Afib., liver injury, all thrombosis cases, and hospital-acquired infection(s)).
Data Collection
The data collection tool was standardized among all centers. The data was collected using Research Electronic Data Capture (REDCap®) form hosted by the primary institution. The form included demographic factors, comorbidities, and clinical data, including vitals and laboratory tests (eg, lactic acid, serum creatinine, liver function tests, and inflammatory markers). The severity and clinical scores were captured within 24 h of ICU admission, such as Acute Physiology and Chronic Health Evaluation II (APACHE II), Sequential Organ Failure Assessment (SOFA), and Multiple Organ Dysfunction Score (MODS). Moreover, Glasgow Coma Score (GCS), acute kidney injury, proning position status, inotropes use, and the need for MV were documented within 24 h of admission. Additionally, the study recorded the early use of tocilizumab, methylprednisolone, and dexamethasone in these patients.
Statistical Analysis
Descriptive analysis was used to present the categorical variables as a number with percentage and continuous data as mean with standard deviation (SD) or median with lower (Q1) and upper (Q3) quartiles based on the distribution of data. The normality assumptions for all numerical variables were evaluated using a statistical test (the Shapiro–Wilk test) and graphical representation (ie, histograms and Q-Q plots).
The baseline characteristics of the two study groups were compared using appropriate statistical tests. Categorical variables were compared using either the Chi-square test or the Fisher exact test. On the other hand, continuous variables were compared using the Student's t-test for normally distributed data [1]. In cases where the continuous variables were not normally distributed, the Mann-Whitney U-test was utilized.
Regression analysis was done by considering the propensity score (PS) as one of the covariates in the model. Multivariable Cox proportional hazards regression analyses were performed for the 30-day and in-hospital mortality and reported using the hazard ratios (HR). The proportionality assumption was assessed before fitting the Cox model. On the other hand, multivariable logistic and negative binomial regression analysis were used for the outcomes considered in this study and reported using the odds ratios (ORs) or estimates with 95% confidence intervals (CIs) as appropriate. Model fit was assessed using the Hosmer–Lemeshow goodness-of-fit test.
In the PS-matched analysis, we selected five matching covariates, including age, baseline serum creatinine, chronic kidney disease as a comorbid condition, APACHE II score, and hematocrit at admission. Those factors were selected for their possible association with the study outcomes. A greedy nearest-neighbor matching method was used in which one patient who received ESA was matched with three patients in the control group (1:3 ratio). This eventually produces the smallest within-pair difference among all available pairs with treated patients. These patients were matched only if the difference in the logits of the PSs for pairs of patients from the two groups was less than or equal to 0.1 times the pooled estimate of the standard deviation. We considered a P-value of < .05 statistically significant and used SAS version 9.4 for all statistical analyses.
Results
A total of 1592 patients with confirmed COVID-19 were admitted to the ICU during the study and evaluated for inclusion. Of these, 135 were excluded from the analysis due to various reasons (Figure 1). The total included number was 1457 patients, distributed as 1424 in the control group and 33 patients in the ESA group. After PS score matching (1:3) based on the predefined criteria, the number of patients included overall was 120. Of the total number, 30 patients were administered ESA (with 47% receiving Darbepoetin alfa and 53% receiving Erythropoietin). The mean weekly dose of Darbepoetin alfa was 40 mcg with a standard deviation of 20.4, while the mean weekly dose for Erythropoietin was 10,500 units with a standard deviation of 2582. Among these patients, 23% had previously received ESA therapy before their admission to the ICU and continued to receive it during their stay, whereas 76.7% received ESA therapy as a new initiation during their ICU stay.
Demographic and Clinical Characteristics
Table 1 summarizes all the baseline characteristics before and after PS matching. In the whole cohort and before the PS matching, the mediam age of the study population was 63.0 (53.0, 73.0) years, and the majority were males (62.3%). The most observed comorbidities were diabetes mellitus (61.6%), followed by hypertension (59.1%), dyslipidemia (22%), and CKD (10.6%). Moreover, patients who received ESA therapy were older, had a higher baseline MODS, serum creatinine, acute kidney injury, INR, aPTT, and a higher prevalence of CKD and IHD compared with the control group. On the other hand, the control group had a higher baseline hematocrit, albumin, and early use of tocilizumab within 24 h of ICU admission than patients who received ESA therapy before PS matching. The two groups were comparable after conducting the PS matching, except for baseline albumin and blood glucose levels within 24 h of ICU admission.
Summary of Demography and Baseline Characteristics.
*T Test / ^ Wilcoxon rank sum test is used to calculate the P-value.
^^ Chi square/ ** Fisher's exact teat is used to calculate P-value.
*$ Nephrotoxic medications/material included IV Vancomycin, Gentamicin, Amikacin, Contrast, Colistin, Furosemide, and/or Sulfamethoxazole/trimethoprim.
Length of Stay and Ventilator Free Days
In crude analysis, patients who received ESA therapy had a statistically significantly longer median of ICU (18 vs. 10) and hospital LOS (20 vs. 16) compared with the control group. Moreover, generalized linear regression analysis revealed a longer duration of ICU stay as well as hospital stay in patients who received ESA compared to the control group (beta coefficient: 0.64; 95% CI: 0.31-0.97; P = < .01, and beta coefficient: 0.41; 95% CI: 0.12-0.69; P = < .01, respectively). On the other hand, the ESA group's median VFDs was statistically significantly shorter than the control group (Table 2).
The Outcomes of Critically ill Patients With COVID-19 After Propensity Score Matching.
Denominator of the percentage is the total number of patients.
^ Wilcoxon rank sum test is used to calculate the P-value.
Chi-square test is used to calculate the P-value.
Cox proportional hazards regression analysis used to calculate HR and P-value.
Generalized linear model is used to calculate estimates and P-value.
30-day and in-Hospital Mortality
The 30-day and in-hospital mortality were higher in patients who received ESA therapy in comparison with the control at crude analysis (71.4% vs 50.6%; P = .06 and 78.6% vs 55.6%; P = .03, respectively). However, Cox proportional hazards regression analysis did not show statistically significant differences between the two groups in terms of 30-day and in-hospital mortality (HR: 1.22; 95% CI: 0.7-2.12, P = .48 and HR: 1.07; 95% CI: 0.62-1.83, P = .82, respectively) (Table 2).
Complications During ICU Stay
Among the ESA group, statistically significantly higher odds of liver injury were observed than the control group (OR: 3.49; 95% CI: 1.03-11.81, P = .04). Moreover, patients who received ESA therapy also had higher odds of hospital-acquired infections (OR: 2.93; 95% CI: 1.03-8.33, P = .04). Other ICU complications were not statistically significant between the two groups (Table 3).
The ICU Complications During Stay.
Denominator of the percentage is patients who don’t have atrial fibrillation as comorbidity.
Denominator of the percentage is the total number of patients.
^^ Chi square/ ** Fisher's exact teat is used to calculate P-value.
$* Logistic regression is used to calculate the OR and P-value.
Discussion
In this multicenter retrospective cohort study, the impact of erythropoiesis-stimulating agents (ESA) on the clinical outcomes of critically ill patients with COVID-19 was examined. After propensity score matching for age, APACHE II score, chronic kidney disease, baseline serum creatinine, and hematocrit at admission, we found that using ESA was associated with longer lengths of stay in the ICU and hospital. This could be attributed to the development of liver injury (20.0% vs 6.7%) and the occurrence of acquired infections during the patient's stay.27,28
In contradiction to our findings, Benjamin et al reported four cases of anemic COVID-19 patients, two were critically ill, and all four patients received EPO analogs, changes from baseline hemoglobin and hematocrit were variable. However, it is good to note that the four patients were smokers in addition to one patient with CKD, which subsequently suggests possible ESA administration pre-COVID-19 infection. Most importantly, the length of ICU and hospital stay were variable (72 h, 8 weeks for non-ICU patients, and 22 days, 3 weeks for ICU patients). Hence, it is not possible to generalize the aforementioned results due to the small sample size, lack of comparable groups, and patients’ baseline characteristics. 11 Moreover, in a recent interventional study conducted by Samimagham et al, the mean hospital LOS was significantly less than that of control group (P = .002). 29 These results were inconsistent with our findings. It is worth mentioning that our cohort was comparable in terms of anemia variables and CKD at baseline.
Liu et al investigated 245 patients with COVID-19 and found a non-significant difference between all-cause mortality during hospitalization and baseline hemoglobin levels, and the OR of death with increasing serum hemoglobin level was 0.98 (95% CI: 0.96, 1.00, P = .05). Our findings came in line with this result; both 30-day mortality and in-hospital mortality were not statistically significant after PS matching (P = .48 and .82, respectively). The number of VFDs during the ICU stay was significantly shorter in the crude and regression analysis in the ESA group (beta coefficient (95% CI): −1.26 (−2.45-0.07), P < .04). Whereas, in a double-blind multicenter randomized control trial for 1302 critically ill patients, ESA therapy was found to reduce blood transfusion without any favorable effect on VFD, hospital stay, or mortality. 13
There was a significant difference between the two groups’ baseline albumin. Hypoalbuminemia has been prominent in the ESA group (mean26.3 (±5.84) in the ESA group vs 31.3 (±4.27) in the control group, P = .0002). A meta-analysis has reported that hypoalbuminemia in COVID-19 patients was found to be associated with a greater risk of hospital-acquired infections, higher viral load, and pronounced organ dysfunction in this group of patients.30,31 he development of hypoalbuminemia in severe infectious diseases can be attributed to several factors, including increased capillary permeability, decreased protein synthesis, and elevated levels of pro-inflammatory cytokines in the bloodstream. These mechanisms contribute to the leakage of albumin from the blood vessels into the interstitial space, resulting in reduced levels of albumin in the blood. 32 In line with this, our study findings also showed significantly higher odds of hospital-acquired infections in the ESA group (P = .04). Moreover, our data shows that there were higher baseline blood glucose levels in both groups and a significant increase in the treatment group, which may be related to an increase in the risk for other hospital-acquired infections and increased hospital stay as well as mortality.
Based on the baseline data from the current cohort, it was observed that both groups had higher blood glucose levels at the baseline, with a significant increase in the ESA group, which may be related to an increase in the risk for other hospital-acquired infections and increased hospital stay as well as mortality. Hyperglycemia is one of the normal physiological responses for critically ill patients. Maintaining the blood glucose within the target ranges was found to significantly decrease mortality as well as morbidity in critically ill patients. 33 At the same time, hyperglycemia was found to decrease immune response and increase susceptibility to infection. 34
Lastly, we found that the ESA group had a significantly higher incidence of liver injury. However, preclinical, and observational studies didn’t report any instances of elevation in liver enzymes or acute liver injury due to ESA therapy, also the mechanism of ESA in causing liver injury is not well known. 35 However, we can attribute this to the low serum albumin in COVID-19 critically ill patients as discussed before which may contribute to end organ damage including hepatic dysfunction.
Our study had some limitations: First, some variables were not possibly to be controlled after PS matching such as albumin level, and blood glucose which might affect the external validity of the results. Second, having a small sample size and a short follow-up period may interpret study outcomes. Lastly, not being willing to assess inflammatory markers made it challenging to assess disease and outcome progression. In summary, our study is unique in the lack of well-controlled studies investigating the association between ESA administration in critically ill patients with COVID-19 in terms of clinical outcomes. Hence, controlled trials with larger sample sizes are warranted to illustrate a more pronounced correlation.
Conclusion
The use of ESAs in critically ill patients with COVID-19 was found to be associated with prolonged hospital and ICU stays. However, there were no observed survival benefits, and the use of ESAs was linked to lower ventilator-free days (VFDs). To validate these findings, it is recommended to conduct randomized controlled trials with larger sample sizes.
Footnotes
Acknowledgments
We would like to thank the investigators who participated in this project from the Saudi critical care pharmacy research (SCAPE) platform.
Author Contributions
All authors contributed to data collection, analysis, drafted, revised, and approved the final version of the manuscript.
Availability of Data and Material
The datasets used and/or analyzed during the current study are available from corresponding author on reasonable request.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Ethics Approval and Consent to Participate
The study was approved in January 2021 by King Abdullah International Medical Research Center (KAIMRC)–Institutional Review Board, Riyadh, Saudi Arabia (Ref.# NRC21R-004-01). Informed consent from the study patients was waived by KAIMRC due to the retrospective observational nature of the study. Participants’ confidentiality was strictly observed throughout the study by using anonymous unique serial number for each subject and restricting data only to the investigators.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
