TY - JOUR AU1 - Elhadi, Muhammed AU2 - Alsoufi, Ahmed AU3 - Abusalama, Abdurraouf AU4 - Alkaseek, Akram AU5 - Abdeewi, Saedah AU6 - Yahya, Mohammed AU7 - Mohammed, Alsnosy AU8 - Abdelkabir, Mohammed AU9 - Huwaysh, Mohammed AU1 - Amkhatirah, Emad AU1 - Alshorbaji, Kamel AU1 - Khel, Samer AU1 - Gamra, Marwa AU1 - Alhadi, Abdulmueti AU1 - Abubaker, Taha AU1 - Anaiba, Mohamed AU1 - Elmugassabi, Mohammed AU1 - Binnawara, Muhannud AU1 - Khaled, Ala AU2 - Zaid, Ahmed AU2 - Msherghi, Ahmed AB - Background The coronavirus disease (COVID-19) pandemic has severely affected African countries, specifically the countries, such as Libya, that are in constant conflict. Clinical and laboratory information, including mortality and associated risk factors in relation to hospital settings and available resources, about critically ill patients with COVID-19 in Africa is not available. This study aimed to determine the mortality and morbidity of COVID-19 patients in intensive care units (ICU) following 60 days after ICU admission, and explore the factors that influence in‐ICU mortality rate. Methods This is a multicenter prospective observational study among COVID-19 critical care patients in 11 ICUs in Libya from May 29th to December 30th 2020. Basic demographic data, clinical characteristics, laboratory values, admission Sequential Organ Failure Assessment (SOFA) score, quick SOFA, and clinical management were analyzed. Result We included 465 consecutive COVID-19 critically ill patients. The majority (67.1%) of the patients were older than 60 years, with a median (IQR) age of 69 (56.5–75); 240 (51.6%) were male. At 60 days of follow-up, 184 (39.6%) were discharged alive, while 281 (60.4%) died in the intensive care unit. The median (IQR) ICU length of stay was 7 days (4–10) and non-survivors had significantly shorter stay, 6 (3–10) days. The body mass index was 27.9 (24.1–31.6) kg/m2. At admission to the intensive care unit, quick SOFA median (IQR) score was 1 (1–2), whereas total SOFA score was 6 (4–7). In univariate analysis, the following parameters were significantly associated with increased/decreased hazard of mortality: increased age, BMI, white cell count, neutrophils, procalcitonin, cardiac troponin, C-reactive protein, ferritin, fibrinogen, prothrombin, and d-dimer levels were associated with higher risk of mortality. Decreased lymphocytes, and platelet count were associated with higher risk of mortality. Quick SOFA and total SOFA scores increase, emergency intubation, inotrope use, stress myocardiopathy, acute kidney injury, arrythmia, and seizure were associated with higher mortality. Conclusion Our study reported the highest mortality rate (60.4%) among critically ill patients with COVID-19 60 days post-ICU admission. Several factors were found to be predictive of mortality, which may help to identify patients at risk of mortality during the ongoing COVID-19 pandemic. 1 Introduction Since the outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in December 2019, the coronavirus disease (COVID-19) has affected more than 137 million people and has resulted in the death of more than 2.96 million patients by April 14th, 2021 [1]. Symptoms range from mild to severe or life-threatening and may require intensive care unit (ICU) admission. The major cause of morbidity and mortality in COVID-19 is the evolution of acute respiratory distress syndrome (ARDS) and multi-organ dysfunction [2–4]. Several factors, such as increasing age, presence of other comorbidities, and abnormalities in specific laboratory results, have been identified as risk factors for COVID-19 severity [5]. Previous reports have suggested that up to 29% of patients develop ARDS and require hospitalization [2]. An early report suggested that 2–12% of patients will require mechanical ventilation and that mortality can range from 26% to 97% [6]. The rate of COVID-19 patients admitted to ICUs ranged between 7% and 26%, depending on the country. Of 138 hospitalized patients in China, 26.1% were transferred to the ICU [7]. In another study of 1099 patients from 522 hospitals in China, 5% were admitted to the ICU [4]. In Italy, 16% of hospitalized patients and 5–12% of all patients were transferred to the ICU [8, 9]. In the United States, a study of 5700 hospitalized patients in New York found that 20% required ICU admission and mechanical ventilation [10]. Another study in the United States found that 4.9% to 11.5% of 2449 hospitalized patients were admitted to the ICU [11]. Patients requiring ICU admission have been reported to have higher rates of mortality and morbidity during ICU admission. Although rates vary according to geographical and income categories, ICU mortality was 26% among 1591 admitted Italian patients [3]. Another study in China reported that among 344 ICU patients, the mortality was 38.7% within 28 days after ICU admission [12]. Another study in Spain found a mortality rate of 23.2% among 237 patients admitted to the ICU [13]. A European study conducted in France, Belgium, and Switzerland found a 28.1% mortality rate in the ICU in a study of 4224 patients [14]. In a meta-analysis of 24 studies outside of Africa, ICU mortality was 41.6% (34.0–49.7%) [15]. However, the only study from the African region was that of Bruce Biccard et al. [16], which found a 48.2% mortality rate among patients in several African countries. However, there have been no other reports on the survival and mortality rates of African ICU patients. Differences in local settings, as well as differences in predisposing factors in a war-torn country, such as Libya, play a major role in increasing mortality and morbidity. The COVID-19 pandemic has caused unprecedented challenges for ICUs worldwide. Africa was severely affected by the COVID-19 pandemic due to the lack of available resources and shortage of equipment, along with a lack of properly trained healthcare workers [17, 18]. These major challenges have resulted in the highest ICU mortality in Africa compared to that in the rest of the world, and Bruce et al. [16] reported a 54.7% mortality rate for African patients after 30 days of ICU admission. This high ratio reflects the severe shortage of healthcare infrastructure and the requirement for facilities in order to fight COVID-19. In Libya, the first laboratory-confirmed case of COVID-19 was detected on March 24, 2020 [19]. Subsequently, the pandemic increased steeply to approximately 170,000 cases and 2,830 deaths by April 14th, 2021. Libya was not prepared for the COVID-19 pandemic. It had a severe shortage of ICU beds and mechanical ventilators [20]. In addition, 56.7% of healthcare workers reported a shortage of personal protective equipment (PPE), and 70% reported purchasing PPE themselves [20]. In this context, assessing mortality and complications, along with ICU needs, is crucial to provide an overview of the current COVID-19 situation in some developing African countries where conducting research is very challenging. Our main objective was to determine the mortality and morbidity of COVID-19 patients 60 days after ICU admission, and to examine individual risk factors associated with ICU mortality. 2 Materials and methods 2.1 Study and selection of patients This is a prospective cohort study of COVID-19 patients consecutively admitted to ICUs in 11 hospitals and isolation centers that managed COVID-19 patients in Libya from the 29th of May,2020 to the 30th of December 2020. All hospitals are public hospitals funded by the government, as private hospitals do not admit COVID-19 patients. Patients recruited in the study were followed up for 60 days after ICU admission or until discharge. The hospitals that participated in the study were originally mixed surgical-medical specialized COVID-19 centers with daily multidisciplinary rounds and specific protocols for patient care during the pandemic. The study included only patients whose COVID-19 infection was confirmed by real-time PCR (RT-PCR) assay for SARS-CoV-2 testing of a nasopharyngeal sample. Inclusion criteria were as follows: adults aged > 18 years admitted to the ICU. Exclusion criteria were negative RT-PCR for COVID-19 or patients transferred to wards or isolation departments. The primary outcome was 60-day mortality after ICU admission. The COVID-19 patient data recorded included basic demographic data, essential medical history, pre-admission medication history, comorbidities, antibiotic and antifungal history, radiological results, applied treatment, organ support indications, intubation history and details, thromboembolic complications and treatment, temperature, basic laboratory results, admission sequential organ failure assessment (SOFA) score, ICU length of stay in days, and status at hospital discharge or 60 days after critical care admission. SOFA score is an objective score, which allows for measurement in six organ systems (respiratory, coagulatory, hepatic, cardiovascular, renal and neurological), with the score measuring individual or multiple organ dysfunction [21, 22]. Each organ has a score range from 0 to 4. SOFA scoring system is used to predict the clinical outcomes of patients admitted to intensive care units [21, 23]. The study adhered to the Strengthening the Reporting of Observational Studies in Epidemiology statement [24], and the recruitment strategy is summarized in a flow chart in Fig 1. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 1. STROBE flow chart. STROBE, Strengthening the Reporting of Observational Studies in Epidemiology. https://doi.org/10.1371/journal.pone.0251085.g001 In Libya, there is no provision for extracorporeal membrane oxygenation (ECMO); in addition, there are no guidelines or policies for withdrawal of life support or to make decisions to limit therapy. 2.2 Statistical analysis The Shapiro–Wilk normality test was used to determine the non-normal distribution of continuous variables. Median and interquartile ranges were used for continuous variables, while categorical variables were presented as frequencies and percentages. For continuous variables, the Mann–Whitney or Kruskal–Wallis H tests were used to determine statistically significant differences between two or more groups of an independent variable and a continuous variable. The chi-square test was used for categorical variables for association and comparison between survivor and non-survivor groups. Phi (φ) and Cramer’s V were used to determine the impact of association in the chi-square test for two or more variables, respectively. Cox regression was performed to calculate hazard ratios (HRs) and 95% confidence intervals for both the mortality and ICU admission outcome variables, where the binary outcome was (survival/death), and the length of stay in the ICU was used as a measurement of time. Univariate and multivariate Cox regression analyses were performed. Kaplan–Meier survival analysis (Kaplan & Meier, 1958) was conducted to compare the study variables and the impact on ICU mortality. A log-rank test was performed to determine the significance of differences in the survival distribution for the study variables and intervention suggested. All tests were two-tailed and considered significant when p < 0.05. Statistical analysis was performed using IBM SPSS Statistics for Windows (IBM SPSS Statistics for Windows, Version 25.0, IBM Corp., Armonk, NY, USA). 2.3 Ethical approval Ethical approval for this study was obtained from the Bioethics Committee at the Biotechnology Research Center in Libya (Reference number: BEC-BTRC-13-2020). Informed consent was not required due to the observational nature of the study and the pandemic situation. 2.1 Study and selection of patients This is a prospective cohort study of COVID-19 patients consecutively admitted to ICUs in 11 hospitals and isolation centers that managed COVID-19 patients in Libya from the 29th of May,2020 to the 30th of December 2020. All hospitals are public hospitals funded by the government, as private hospitals do not admit COVID-19 patients. Patients recruited in the study were followed up for 60 days after ICU admission or until discharge. The hospitals that participated in the study were originally mixed surgical-medical specialized COVID-19 centers with daily multidisciplinary rounds and specific protocols for patient care during the pandemic. The study included only patients whose COVID-19 infection was confirmed by real-time PCR (RT-PCR) assay for SARS-CoV-2 testing of a nasopharyngeal sample. Inclusion criteria were as follows: adults aged > 18 years admitted to the ICU. Exclusion criteria were negative RT-PCR for COVID-19 or patients transferred to wards or isolation departments. The primary outcome was 60-day mortality after ICU admission. The COVID-19 patient data recorded included basic demographic data, essential medical history, pre-admission medication history, comorbidities, antibiotic and antifungal history, radiological results, applied treatment, organ support indications, intubation history and details, thromboembolic complications and treatment, temperature, basic laboratory results, admission sequential organ failure assessment (SOFA) score, ICU length of stay in days, and status at hospital discharge or 60 days after critical care admission. SOFA score is an objective score, which allows for measurement in six organ systems (respiratory, coagulatory, hepatic, cardiovascular, renal and neurological), with the score measuring individual or multiple organ dysfunction [21, 22]. Each organ has a score range from 0 to 4. SOFA scoring system is used to predict the clinical outcomes of patients admitted to intensive care units [21, 23]. The study adhered to the Strengthening the Reporting of Observational Studies in Epidemiology statement [24], and the recruitment strategy is summarized in a flow chart in Fig 1. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 1. STROBE flow chart. STROBE, Strengthening the Reporting of Observational Studies in Epidemiology. https://doi.org/10.1371/journal.pone.0251085.g001 In Libya, there is no provision for extracorporeal membrane oxygenation (ECMO); in addition, there are no guidelines or policies for withdrawal of life support or to make decisions to limit therapy. 2.2 Statistical analysis The Shapiro–Wilk normality test was used to determine the non-normal distribution of continuous variables. Median and interquartile ranges were used for continuous variables, while categorical variables were presented as frequencies and percentages. For continuous variables, the Mann–Whitney or Kruskal–Wallis H tests were used to determine statistically significant differences between two or more groups of an independent variable and a continuous variable. The chi-square test was used for categorical variables for association and comparison between survivor and non-survivor groups. Phi (φ) and Cramer’s V were used to determine the impact of association in the chi-square test for two or more variables, respectively. Cox regression was performed to calculate hazard ratios (HRs) and 95% confidence intervals for both the mortality and ICU admission outcome variables, where the binary outcome was (survival/death), and the length of stay in the ICU was used as a measurement of time. Univariate and multivariate Cox regression analyses were performed. Kaplan–Meier survival analysis (Kaplan & Meier, 1958) was conducted to compare the study variables and the impact on ICU mortality. A log-rank test was performed to determine the significance of differences in the survival distribution for the study variables and intervention suggested. All tests were two-tailed and considered significant when p < 0.05. Statistical analysis was performed using IBM SPSS Statistics for Windows (IBM SPSS Statistics for Windows, Version 25.0, IBM Corp., Armonk, NY, USA). 2.3 Ethical approval Ethical approval for this study was obtained from the Bioethics Committee at the Biotechnology Research Center in Libya (Reference number: BEC-BTRC-13-2020). Informed consent was not required due to the observational nature of the study and the pandemic situation. 3 Results 3.1 Patient baseline characteristics The study included 465 critically ill COVID-19 patients in this study. The study excluded 69 patients, including patients who refused to participate in the study, two who were < 18 years of age, and 22 patients for whom > 30% of data were not available (study enrollment flowchart is shown in Fig 1). At 60 days of follow-up, 184 (39.6%) patients were discharged alive from the hospital, whereas 281 (60.4%) died in the ICU. The majority of the patients were older than 60 years, with a median (IQR) age of 69 (56.5–75) years; 240 (51.6%) were male and 225 (48.4%) were female. The body mass index of the participants was 27.9 (24.1–31.6). Of the patients included, only eight (1.7%) were healthcare workers and five (1.1%) were pregnant women. A total of 264 (56.8%) patients had diabetes, and 259 (55.7%) had arterial hypertension. Among the study participants, 225 (48.4%) were transferred from the emergency department to the ICU, 121 (26%) were transferred from the in-hospital ward, and 119 (25.6%) were transferred from another hospital or healthcare facility. None of the patients included in the study had a history of human immunodeficiency virus (HIV) infection. The characteristics of the study participants are listed in Table 1, with a comparison between surviving and non-surviving patients. Median laboratory values were significantly different between survivors and non-survivors; white blood cells, neutrophils, procalcitonin, cardiac troponin, ferritin, fibrinogen, prothrombin, and D-dimer levels were higher in non-survivors than in survivors, while lymphocyte counts and platelet counts were significantly lower in non-survivors than in survivors. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 1. Basic and laboratory characteristics of COVID-19 patients in the intensive care unit during admission. https://doi.org/10.1371/journal.pone.0251085.t001 3.2 Clinical status at admission At admission to the ICU, the quick SOFA median (IQR) score was 1 (1–2), and the total SOFA score was 6 (range, 4–7). The means (SD) of quick SOFA and SOFA were 1.48 (0.84) and 6.04 (2.53), with a range from 0 to 3 and 1 to 19, respectively. In the non-survivor group (281/465), the quick SOFA and SOFA scores were 1.62 (0.8) and 6.85 (2.63), whereas in the survivor group (184/465) were 1.27 (0.79) and 4.81 (1.75), respectively. In the non-survivor group, the quick SOFA and total SOFA scores ranged from 0 to 3 and 1 to 19, whereas in the survivor group, they ranged from 0 to 3 and 1 to 9, respectively. When analyzing the composition of the quick SOFA score, 412 (88.6%) patients presented with a respiratory rate ≥ 22, 157 (33.8%) presented with a Glasgow Coma Scale score < 15, and 120 (25.8%) presented with a systolic blood pressure ≤ 100. Regarding the total SOFA score categories, 427 (91.8%) patients presented with a SOFA score of ≤ 9, which predicts a risk of mortality of less than 33.3%, whereas 28 (6%) scored between 10 and 12, which predicts an approximately 50% risk of mortality, and 10 (2.2%) presented with a SOFA score > 12, which predicts a mortality rate of 95.2%. However, none of the survivors had a SOFA score of > 9. The study found a statistically significant association between quick and full SOFA scores and mortality in patients. Table 2 provides an overview of the SOFA score and quick SOFA score at admission. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 2. Sequential Organ Failure Assessment (SOFA) and quick SOFA (qSOFA) score at admission. https://doi.org/10.1371/journal.pone.0251085.t002 3.3 Intensive care management An overview of the detailed management of ICUs is presented in Table 3. The median (IQR) ICU length of stay was 7 days (range, 4–10 days), with a significantly shorter stay among non-survivors. At the end of the study, 460 (98.9%) patients were either discharged or died at 30 days. Only five (1.1%) patients had a hospital stay of more than 30 days; two of them died at 32 and 38 days, and three were discharged at 35, 39, and 45 days of ICU admission, as shown in Fig 2A, which illustrates the distribution of mortality according to age, sex, and length of ICU stay, and Fig 2B, which shows mortality. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 2. (A) Distribution of mortality of critically ill patients with COVID-19 stratified according to age and sex. (B) Distribution of length of stay in critical care, in days, for COVID-19 patients stratified according to age and mortality. https://doi.org/10.1371/journal.pone.0251085.g002 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 3. Admission management and characteristics of COVID-19 patients in the intensive care unit. https://doi.org/10.1371/journal.pone.0251085.t003 A total of 154 (33.1%) patients developed infection during their ICU stay. Inotrope use, prone position with mechanical ventilation, thromboembolic events during admission, stress myocardiopathy, acute kidney injury, arrhythmia, myocarditis, atelectasis, prolonged delirium, seizures, decubitus ulcers, and endotracheal tube obstruction were significantly higher in the non-survivor group than in the survivor group. Antibiotics were administered to 433 (93.1%) patients; antibiotics, especially quinolones, macrolides, and penicillin, were significantly associated with the survivor group. Steroids were not associated with any group, as both received steroids at similar rates. Other therapies, such as convalescent plasma, blood transfusion, Remdesivir, Tocilizumab, Famotidine, and high or low doses of vitamin C or vitamin D, did not show significant association with either group due to the low frequency of use of these drugs. Respiratory complications were the most frequent, accounting for 435 cases (93.5%) of intensive care complications. A total of 127 cases of acute kidney injury (27.3%), 93 cases of atelectasis (20%), and 88 cases of cardiac arrhythmia (18.9%) were also major events during ICU stay. Thromboembolic events were higher during admission than during pre-admission. An overview of the detailed complications and management of ICU patients is presented in Table 3. 3.4 Factors associated with survival Patients were categorized into survivors (184/465) and non-survivors (281/465) at the end of the study follow-up, or at discharge or death. Univariate and multivariate Cox proportional regression hazard analyses were performed on the study variables, as detailed in Table 4. Multivariate analysis showed that only the following variables were significantly associated with mortality: lower lymphocyte count (HR = 0.35 [0.18, 0.68]), higher procalcitonin (HR = 1.89 [1.03, 3.46]), cardiac troponin (HR = 1.004 [1.002, 1.006]), C-reactive protein (HR = 1.004 [1.001, 1.007]), D-dimer (HR = 1.02 [1.001, 1.03]), total SOFA score (HR = 1.16 [1.07, 1.25]), emergency intubation (HR = 1.97 [1.37, 2.85]), and stress cardiomyopathy (HR = 2.31 [1.27, 4.16]). In univariate analysis, the following were significantly associated with increase/decrease hazard of mortality: increased age (HR = 1.0 [1, 1.02]), BMI (HR = 1.08 [10.6, 1.11]), white cell count (HR = 1.12 [1.09, 1.14]), neutrophil (HR = 1.07 [1.05–1.1]), procalcitonin (HR = 7.83 [5.25, 11.65]), cardiac troponin (HR = 1.009 [1.007, 1.011]), C-reactive protein (HR = 1.011 [1.009–1.014]), ferritin (HR = 1.00 [1.00, 1.001]), fibrinogen (HR = 1.002 [1.001–1.003]), prothrombin (HR = 1.09 [1.04, 1.14]), and D-dimer (HR = 1.05 [1.04, 1.06]). Decreased lymphocyte (HR = 0.14 [0.08, 0.21]) and platelet count (HR = 0.993 [0.991, 0.995]) were both associated with increased mortality. However, increased antibiotic use reduced mortality (HR = 0.56 [0.38, 0.83]). Both increased quick SOFA score (HR = 1.05 [1.04, 1.06]), and increased total SOFA score (HR = 1.12 [1.08, 1.26]) increased the risk of mortality. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 4. Association of patient characteristics with 60-day post-ICU admission mortality using ICU length of stay as time variable (univariate/multivariate cox regression models). https://doi.org/10.1371/journal.pone.0251085.t004 Age > 70 years (HR = 1.52 [1.14, 2.03]), elective intubation (HR = 1.99 [1.33, 2.96]), emergency intubation (HR = 2.63 [2.03, 3.41]), inotrope use (HR = 1.72 [1.36, 2.19]), stress myocardiopathy (HR = 1.71 [1.203, 2.42]), acute kidney injury (HR = 1.37 [1.08, 1.75]), arrhythmia (HR = 1.38 [1.06, 1.81]), and seizure (HR = 1.52 [1.15, 1.99]) were associated with a higher risk of mortality. Kaplan–Meier survival analysis was performed to assess the association between the study variables and survival. The ICU length of stay was used as the mean time. The distribution of the study variables is summarized in Table 1. A log-rank test was conducted to determine if there were differences in the survival distributions of the study variables. The survival distributions for the following variables were significantly different: higher age groups were associated with lower survival rate (χ2(2) = 12.288, p = 0.002), BMI > 30 was associated with lower survival rate (χ2(2) = 20.372, p < 0.001) and antibiotic use was significantly associated with better survival (χ2(2) = 9.37, p = 0.002). The results found that the quick SOFA and SOFA scores were associated with mortality and survival. Lower quick SOFA scores were associated with better survival (χ2(2) = 24.467, p < 0.001). Lower total SOFA scores were significantly associated with better survival outcomes (χ2(2) = 50.852, p < 0.001). When dividing the SOFA score into three categories (≤ 9, 10–12, and > 12), the results showed a significant difference in survival based on the SOFA risk score category. High scores were associated with higher mortality (χ2(2) = 14.38, p = 0.001), and the inverse was also true. The type of intubation was another variable that could predict mortality; emergency intubation was associated with a high risk of mortality compared to elective or no intubation (χ2(2) = 22.45, p = 0.001). Patients who received deep vein thrombosis (DVT) prophylaxis during the first 12 hours of admission had better outcomes and lower mortality (χ2(2) = 24.689, p < 0.001). Patients who received inotropes or vasopressors had higher mortality rates than those who did not (χ2(2) = 22.45, p < 0.001). When comparing the effects of some newly used treatments, such as convalescent plasma, blood transfusion, remdesivir, tocilizumab, famotidine, and high or low doses of vitamin C or vitamin D; however, the results did not find a significant difference in survival distribution, which may be due to the limited number of patients who received these therapies. The majority of patients received steroids; however, there was no significant difference in the survival distribution (p = 0.43). Figs 3 and 4 provide Kaplan–Meier estimates of overall survival based on the previous significant variables. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 3. Kaplan–Meier survival estimates during 60 days following ICU admission. (A) Age grades (≤ 60, 60–70, > 70 years). (B) BMI (≤ 30, > 30). (C) Quick SOFA (0, 1, 2, 3). (D) Total SOFA Score. https://doi.org/10.1371/journal.pone.0251085.g003 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 4. Kaplan–Meier survival estimates during 60 days following ICU admission. (A) Antibiotic (yes, no). (B) SOFA score category (≤ 9, 10–12, >12). (C) Type of intubation (no intubation, elective, emergency). (D) Inotrope use (yes, no). (E) DVT prophylaxis during first 12 hours. https://doi.org/10.1371/journal.pone.0251085.g004 3.5 Hospital resources and infrastructure Regarding hospital facility infrastructure and availability of equipment among the participating hospitals, seven (63.5%) were teaching hospitals, whereas four (36.4%) were not teaching hospitals. The median number of ICU beds (IQR) was 15 (range, 6–16), with a range from 4 to 18 ICU beds. The number of mechanical ventilators was 9 (range, 7–12), with a range from 6 to 20 available devices. The patient-to-nurse ratio ranged from 2 to 12, with a median (IQR) of 4 (range, 2–4). The number of intensivist specialists ranged from 1 to 6 specialists, with a mean (IQR) of 2 (1–4), whereas the number of doctors in training ranged from 2 to 9, with a median (IQR) of 4 (2–6). All hospitals had non-invasive and invasive ventilation facilities. None of the hospitals had extracorporeal membrane oxygenation (ECMO). None of the hospitals had a 24 h visitor policy; only five (45.5%) had a visitor policy of limited hours, and two (18.2%) provided end-of-life visitation. Table 5 summarizes the hospital infrastructure and facilities. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 5. Hospital infrastructure and facilities (n = 11). https://doi.org/10.1371/journal.pone.0251085.t005 3.1 Patient baseline characteristics The study included 465 critically ill COVID-19 patients in this study. The study excluded 69 patients, including patients who refused to participate in the study, two who were < 18 years of age, and 22 patients for whom > 30% of data were not available (study enrollment flowchart is shown in Fig 1). At 60 days of follow-up, 184 (39.6%) patients were discharged alive from the hospital, whereas 281 (60.4%) died in the ICU. The majority of the patients were older than 60 years, with a median (IQR) age of 69 (56.5–75) years; 240 (51.6%) were male and 225 (48.4%) were female. The body mass index of the participants was 27.9 (24.1–31.6). Of the patients included, only eight (1.7%) were healthcare workers and five (1.1%) were pregnant women. A total of 264 (56.8%) patients had diabetes, and 259 (55.7%) had arterial hypertension. Among the study participants, 225 (48.4%) were transferred from the emergency department to the ICU, 121 (26%) were transferred from the in-hospital ward, and 119 (25.6%) were transferred from another hospital or healthcare facility. None of the patients included in the study had a history of human immunodeficiency virus (HIV) infection. The characteristics of the study participants are listed in Table 1, with a comparison between surviving and non-surviving patients. Median laboratory values were significantly different between survivors and non-survivors; white blood cells, neutrophils, procalcitonin, cardiac troponin, ferritin, fibrinogen, prothrombin, and D-dimer levels were higher in non-survivors than in survivors, while lymphocyte counts and platelet counts were significantly lower in non-survivors than in survivors. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 1. Basic and laboratory characteristics of COVID-19 patients in the intensive care unit during admission. https://doi.org/10.1371/journal.pone.0251085.t001 3.2 Clinical status at admission At admission to the ICU, the quick SOFA median (IQR) score was 1 (1–2), and the total SOFA score was 6 (range, 4–7). The means (SD) of quick SOFA and SOFA were 1.48 (0.84) and 6.04 (2.53), with a range from 0 to 3 and 1 to 19, respectively. In the non-survivor group (281/465), the quick SOFA and SOFA scores were 1.62 (0.8) and 6.85 (2.63), whereas in the survivor group (184/465) were 1.27 (0.79) and 4.81 (1.75), respectively. In the non-survivor group, the quick SOFA and total SOFA scores ranged from 0 to 3 and 1 to 19, whereas in the survivor group, they ranged from 0 to 3 and 1 to 9, respectively. When analyzing the composition of the quick SOFA score, 412 (88.6%) patients presented with a respiratory rate ≥ 22, 157 (33.8%) presented with a Glasgow Coma Scale score < 15, and 120 (25.8%) presented with a systolic blood pressure ≤ 100. Regarding the total SOFA score categories, 427 (91.8%) patients presented with a SOFA score of ≤ 9, which predicts a risk of mortality of less than 33.3%, whereas 28 (6%) scored between 10 and 12, which predicts an approximately 50% risk of mortality, and 10 (2.2%) presented with a SOFA score > 12, which predicts a mortality rate of 95.2%. However, none of the survivors had a SOFA score of > 9. The study found a statistically significant association between quick and full SOFA scores and mortality in patients. Table 2 provides an overview of the SOFA score and quick SOFA score at admission. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 2. Sequential Organ Failure Assessment (SOFA) and quick SOFA (qSOFA) score at admission. https://doi.org/10.1371/journal.pone.0251085.t002 3.3 Intensive care management An overview of the detailed management of ICUs is presented in Table 3. The median (IQR) ICU length of stay was 7 days (range, 4–10 days), with a significantly shorter stay among non-survivors. At the end of the study, 460 (98.9%) patients were either discharged or died at 30 days. Only five (1.1%) patients had a hospital stay of more than 30 days; two of them died at 32 and 38 days, and three were discharged at 35, 39, and 45 days of ICU admission, as shown in Fig 2A, which illustrates the distribution of mortality according to age, sex, and length of ICU stay, and Fig 2B, which shows mortality. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 2. (A) Distribution of mortality of critically ill patients with COVID-19 stratified according to age and sex. (B) Distribution of length of stay in critical care, in days, for COVID-19 patients stratified according to age and mortality. https://doi.org/10.1371/journal.pone.0251085.g002 Download: PPT PowerPoint slide PNG larger image TIFF original image Table 3. Admission management and characteristics of COVID-19 patients in the intensive care unit. https://doi.org/10.1371/journal.pone.0251085.t003 A total of 154 (33.1%) patients developed infection during their ICU stay. Inotrope use, prone position with mechanical ventilation, thromboembolic events during admission, stress myocardiopathy, acute kidney injury, arrhythmia, myocarditis, atelectasis, prolonged delirium, seizures, decubitus ulcers, and endotracheal tube obstruction were significantly higher in the non-survivor group than in the survivor group. Antibiotics were administered to 433 (93.1%) patients; antibiotics, especially quinolones, macrolides, and penicillin, were significantly associated with the survivor group. Steroids were not associated with any group, as both received steroids at similar rates. Other therapies, such as convalescent plasma, blood transfusion, Remdesivir, Tocilizumab, Famotidine, and high or low doses of vitamin C or vitamin D, did not show significant association with either group due to the low frequency of use of these drugs. Respiratory complications were the most frequent, accounting for 435 cases (93.5%) of intensive care complications. A total of 127 cases of acute kidney injury (27.3%), 93 cases of atelectasis (20%), and 88 cases of cardiac arrhythmia (18.9%) were also major events during ICU stay. Thromboembolic events were higher during admission than during pre-admission. An overview of the detailed complications and management of ICU patients is presented in Table 3. 3.4 Factors associated with survival Patients were categorized into survivors (184/465) and non-survivors (281/465) at the end of the study follow-up, or at discharge or death. Univariate and multivariate Cox proportional regression hazard analyses were performed on the study variables, as detailed in Table 4. Multivariate analysis showed that only the following variables were significantly associated with mortality: lower lymphocyte count (HR = 0.35 [0.18, 0.68]), higher procalcitonin (HR = 1.89 [1.03, 3.46]), cardiac troponin (HR = 1.004 [1.002, 1.006]), C-reactive protein (HR = 1.004 [1.001, 1.007]), D-dimer (HR = 1.02 [1.001, 1.03]), total SOFA score (HR = 1.16 [1.07, 1.25]), emergency intubation (HR = 1.97 [1.37, 2.85]), and stress cardiomyopathy (HR = 2.31 [1.27, 4.16]). In univariate analysis, the following were significantly associated with increase/decrease hazard of mortality: increased age (HR = 1.0 [1, 1.02]), BMI (HR = 1.08 [10.6, 1.11]), white cell count (HR = 1.12 [1.09, 1.14]), neutrophil (HR = 1.07 [1.05–1.1]), procalcitonin (HR = 7.83 [5.25, 11.65]), cardiac troponin (HR = 1.009 [1.007, 1.011]), C-reactive protein (HR = 1.011 [1.009–1.014]), ferritin (HR = 1.00 [1.00, 1.001]), fibrinogen (HR = 1.002 [1.001–1.003]), prothrombin (HR = 1.09 [1.04, 1.14]), and D-dimer (HR = 1.05 [1.04, 1.06]). Decreased lymphocyte (HR = 0.14 [0.08, 0.21]) and platelet count (HR = 0.993 [0.991, 0.995]) were both associated with increased mortality. However, increased antibiotic use reduced mortality (HR = 0.56 [0.38, 0.83]). Both increased quick SOFA score (HR = 1.05 [1.04, 1.06]), and increased total SOFA score (HR = 1.12 [1.08, 1.26]) increased the risk of mortality. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 4. Association of patient characteristics with 60-day post-ICU admission mortality using ICU length of stay as time variable (univariate/multivariate cox regression models). https://doi.org/10.1371/journal.pone.0251085.t004 Age > 70 years (HR = 1.52 [1.14, 2.03]), elective intubation (HR = 1.99 [1.33, 2.96]), emergency intubation (HR = 2.63 [2.03, 3.41]), inotrope use (HR = 1.72 [1.36, 2.19]), stress myocardiopathy (HR = 1.71 [1.203, 2.42]), acute kidney injury (HR = 1.37 [1.08, 1.75]), arrhythmia (HR = 1.38 [1.06, 1.81]), and seizure (HR = 1.52 [1.15, 1.99]) were associated with a higher risk of mortality. Kaplan–Meier survival analysis was performed to assess the association between the study variables and survival. The ICU length of stay was used as the mean time. The distribution of the study variables is summarized in Table 1. A log-rank test was conducted to determine if there were differences in the survival distributions of the study variables. The survival distributions for the following variables were significantly different: higher age groups were associated with lower survival rate (χ2(2) = 12.288, p = 0.002), BMI > 30 was associated with lower survival rate (χ2(2) = 20.372, p < 0.001) and antibiotic use was significantly associated with better survival (χ2(2) = 9.37, p = 0.002). The results found that the quick SOFA and SOFA scores were associated with mortality and survival. Lower quick SOFA scores were associated with better survival (χ2(2) = 24.467, p < 0.001). Lower total SOFA scores were significantly associated with better survival outcomes (χ2(2) = 50.852, p < 0.001). When dividing the SOFA score into three categories (≤ 9, 10–12, and > 12), the results showed a significant difference in survival based on the SOFA risk score category. High scores were associated with higher mortality (χ2(2) = 14.38, p = 0.001), and the inverse was also true. The type of intubation was another variable that could predict mortality; emergency intubation was associated with a high risk of mortality compared to elective or no intubation (χ2(2) = 22.45, p = 0.001). Patients who received deep vein thrombosis (DVT) prophylaxis during the first 12 hours of admission had better outcomes and lower mortality (χ2(2) = 24.689, p < 0.001). Patients who received inotropes or vasopressors had higher mortality rates than those who did not (χ2(2) = 22.45, p < 0.001). When comparing the effects of some newly used treatments, such as convalescent plasma, blood transfusion, remdesivir, tocilizumab, famotidine, and high or low doses of vitamin C or vitamin D; however, the results did not find a significant difference in survival distribution, which may be due to the limited number of patients who received these therapies. The majority of patients received steroids; however, there was no significant difference in the survival distribution (p = 0.43). Figs 3 and 4 provide Kaplan–Meier estimates of overall survival based on the previous significant variables. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 3. Kaplan–Meier survival estimates during 60 days following ICU admission. (A) Age grades (≤ 60, 60–70, > 70 years). (B) BMI (≤ 30, > 30). (C) Quick SOFA (0, 1, 2, 3). (D) Total SOFA Score. https://doi.org/10.1371/journal.pone.0251085.g003 Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 4. Kaplan–Meier survival estimates during 60 days following ICU admission. (A) Antibiotic (yes, no). (B) SOFA score category (≤ 9, 10–12, >12). (C) Type of intubation (no intubation, elective, emergency). (D) Inotrope use (yes, no). (E) DVT prophylaxis during first 12 hours. https://doi.org/10.1371/journal.pone.0251085.g004 3.5 Hospital resources and infrastructure Regarding hospital facility infrastructure and availability of equipment among the participating hospitals, seven (63.5%) were teaching hospitals, whereas four (36.4%) were not teaching hospitals. The median number of ICU beds (IQR) was 15 (range, 6–16), with a range from 4 to 18 ICU beds. The number of mechanical ventilators was 9 (range, 7–12), with a range from 6 to 20 available devices. The patient-to-nurse ratio ranged from 2 to 12, with a median (IQR) of 4 (range, 2–4). The number of intensivist specialists ranged from 1 to 6 specialists, with a mean (IQR) of 2 (1–4), whereas the number of doctors in training ranged from 2 to 9, with a median (IQR) of 4 (2–6). All hospitals had non-invasive and invasive ventilation facilities. None of the hospitals had extracorporeal membrane oxygenation (ECMO). None of the hospitals had a 24 h visitor policy; only five (45.5%) had a visitor policy of limited hours, and two (18.2%) provided end-of-life visitation. Table 5 summarizes the hospital infrastructure and facilities. Download: PPT PowerPoint slide PNG larger image TIFF original image Table 5. Hospital infrastructure and facilities (n = 11). https://doi.org/10.1371/journal.pone.0251085.t005 4 Discussion This study described the clinical characteristics and presentation of critically ill Libyan patients with COVID-19. In the analysis, the study assessed several basic, clinical, laboratory, and management characteristics and outcomes of severely ill patients who were admitted to ICUs, where followed them up to hospital discharge or mortality at 60 days. The main finding is the 60.4% mortality rate following ICU admission, which, to our knowledge, is the highest recorded compared to all countries’ previous studies. The study also found that age, BMI, laboratory findings, admission SOFA score, and quick SOFA score were strong predictive risk factors for mortality. Mortality among critically ill COVID-19 patients has been reported in several previous studies. However, there is a scarcity of information on COVID-19 status in Africa, which is at a higher risk due to resource limitations and several challenges that are endured by healthcare systems in African countries. One study conducted in several African countries found a mortality rate of 54.7% 30 days after ICU admission [16]. A previous study conducted in several countries revealed different mortality rates and associated risk factors. A study conducted in Wuhan, China, in February 2020, found that among 226 COVID-19 critical care patients, 87 (38.5%) died [25]. Another single-center study in China also found that the mortality rate among 344 critically ill patients was 38.7% (133) at or before 28 days [12]. The study also confirmed that lymphopenia is associated with mortality in critically ill patients with COVID-19. Another single-center study in Spain [13] of 237 critically ill COVID-19 patients found that 55 died and 116 remained in the ICU by April 19, 2020. In a large multicenter prospective study in three European countries including Belgium, France, and Switzerland that included 4224 critically ill patients with COVID-19, the mortality rate was 31% (1298) on day 90 after ICU admission [14]. Another multicenter study in Spain and Andorra [26] of 960 patients reported a mortality rate of 31% (203). Another multicenter prospective study in France [27] of 335 critically ill patients reported a mortality rate of 7.76% (26). In a single-center study in Sweden of 260 participants, mortality was found in 60 of those admitted to the ICU [28]. In Germany, a multicenter prospective study of 223 patients found a mortality rate of 35% (78) in ICU patients [29]. In a large multicenter prospective study in European countries, among 639 patients admitted to ICUs, ICU mortality was 24% [30]. A retrospective multicenter study in Georgia in the United States found a mortality rate of 28.6% (62 out of 217) among critically ill patients with COVID-19 [31]. Another retrospective single center in Ohio [32] of 495 critical care patients found that 91 (18.4%) died in the ICU. A study in Michigan of 141 patients found a mortality frequency of 40.4% (57) at 30 days of admission [33]. In Australia, a study of 18 patients admitted reported death of four (22%) patients in the ICU [34]. Our study demonstrated the utilization of both the SOFA score and quick SOFA score at admission. The results found a median (IQR) quick SOFA score of 1 (1–2) and a total SOFA score of 6 (4–7) on admission. When compared with other studies, a single-center retrospective study in Singapore in 22 patients with a 2/22 mortality rate found an initial admission median (IQR) SOFA score of 2.5 (1.25–7) [35]. A single-center study in Turkey between March and April 2020 found a SOFA score of 4 (3–6) at admission, with a mortality rate of 50.5% (52 out of 103) [36]. Another study of 55 patients in the ICU in the United Arab Emirates reported a SOFA score of 4 (3–7) with a mortality rate of 22% (12) [37]. A prospective multicenter study of 226 COVID-19 critical care patients in China reported a SOFA score of 4 (2–8) at admission [25]. In the COVID-ICU group, a study of 4224 European COVID-19 critical care patients found a SOFA score of 5 (3–8) [14]. In a Spanish study [26], of 960 critically ill patients had a SOFA score of 5 (3–7), and the overall mortality rate was 31% (203 out of 960). In a Swiss single-center prospective study of 129 patients, the SOFA score was 6 (4–7); however, the mortality rate was 19% (24 out of 129) [38]. Another multicenter study conducted in Germany found a SOFA score of 5 (3–9), with a mortality rate of 35% (78 out of 223) [29]. A retrospective study in Canada reported a SOFA score of 5 (3–8), and 18 patients died out of 106 (17%) admitted [39]. Therefore, both SOFA and quick SOFA scores are useful tools for predicting mortality among ICU patients with COVID-19. Our study showed that both were significantly associated with higher ICU mortality. When the quick SOFA score was divided into two categories, 0–1 vs. 2–3. The results showed a significant difference in mortality between the two categories. Similar results were obtained when the total SOFA score was ≤ 9 vs. 10–12 vs. > 12. Both tools indicated that mortality and non-survivorship were significantly associated with higher SOFA or quick SOFA scores. Thus, the quick SOFA score can be used in low-resource settings as an initial assessment of mortality risk to predict outcomes. Calculating the total SOFA score is challenging, especially in low-resource settings, as it requires extensive laboratory results, which can be difficult to obtain in countries with limited resources. Therefore, this recommend the use of quick SOFA as it only requires clinical variables such as Glasgow coma scale < 15, respiratory rate ≥ 22, and systolic blood pressure ≤ 100, in which each one is given one point per answer. Detecting patients with a high risk of mortality can help provide specific care by healthcare providers to reduce the risk of death in the ICU. It can also be a usable tool to prioritize the management of patients at increased risk for complications and mortality, given the treatment and management option limitations. Our study demonstrated high SOFA and quick SOFA scores at admission, which, when compared to previous reports, indicates the severity of study participants and explains the high mortality observed in this study. Furthermore, this study was conducted during the surge in cases in Libya, when many patients could not be admitted appropriately and in a timely manner to the ICU due to a lack of resources and shortage of ICU beds in comparison to the high number of patients. This resulted in some patients waiting for hours or days to find a suitable ICU bed at COVID-19 isolation centers. This delay in appropriate management carries a high risk of hypoxia and complications, as demonstrated by the high SOFA and quick SOFA scores at admission. The study also observed that some patients had to search for ICU beds in several cities, which has major consequences, especially in frail and elderly patients who cannot tolerate delays in the management of hypoxia and respiratory issues. Another possible explanation is the shortage of healthcare supplies and the lack of adequate training of healthcare workers. Libyan ICUs were not prepared in terms of equipment and facilities to fight COVID-19; there was a severe shortage of supplies and a low number of mechanical ventilators, as demonstrated in a previous report [20]. The lack of training and the low number of trained intensive care specialists is another obstacle, both in Libya and in the rest of Africa. In this study, it reported 2 (1–4) intensive care specialists per hospital and a patient to nurse ratio of 4 (2–4), which are not enough to cover the shifts and manage patients appropriately and accurately. Therefore, health authorities need to increase the resilience of hospitals, especially during the second wave. Hospitals need to be equipped with more ICU beds to avoid delays in admission and a higher number of mechanical ventilators to treat those patients who present with respiratory compromise. As in our report, the shortage of mechanical ventilators can significantly affect the ability of Libya and other African countries to treat critically ill patients with COVID-19. In addition, there is a need to recruit more ICU nurses and intensive care specialists to increase the ability of hospitals to manage COVID-19 patients. There is also a need to provide adequate training for healthcare providers and to provide management guidelines that can help utilize the resources and provide hospitals with newly developed drug therapies that have been proven to increase survival and reduce mortality, which are not adequately available in Libya due to their high prices and worldwide demands. However, the government should invest more in healthcare facilities and provide adequate supplies of these drugs in order to mitigate the catastrophic consequences of COVID-19 in critically ill patients, as demonstrated in our report. Increased age was found to be associated with an increased likelihood of mortality among critically ill patients. The study also found that those aged > 70 years had the highest mortality rate compared to the other age groups, which is in line with previous studies [40]. The results found that increased BMI was associated with a higher risk of mortality among COVID-19 critical care patients. This is in line with previous studies, as obesity and increased BMI carry a risk of higher hospital admission rates, morbidity, and mortality for COVID-19 patients [41, 42]. Increased blood laboratory values, such as white cell count and neutrophils, were significantly associated with higher mortality risk, which is in line with previous reports [43]. In contrast, a lower lymphocyte count was significantly associated with a higher risk of mortality, which has been reported previously [44]. In addition, a lower platelet count was significantly associated with mortality rate, which has also been reported in the literature as a risk factor [45]. Regarding other laboratory findings, procalcitonin levels were significantly associated with a sevenfold higher risk of mortality in critically ill patients with COVID-19, which has also been reported in several previous studies as a predictor of mortality [46]. Other biomarkers, such as cardiac troponin, C-reactive protein, ferritin, fibrinogen, prothrombin time, and D-dimer were shown to be associated with a higher risk of mortality among critically ill COVID-19 patients in our study. These biomarkers have been reported in several studies and meta-analyses regarding their role as risk factors and factors predicting mortality and worsening of COVID-19 severity [47–53]. Antibiotic use was significantly associated with lower mortality, as previously reported. These findings may support the use of antibiotics in COVID-19 patients in order to prevent secondary bacterial infection or infection acquired in the ICU, as patients with sepsis were more likely to die in critical care. However, antibiotic use among COVID-19 patients is still debatable, and some previous reports did not support their use [54]. DVT is considered a potential complication and life-threatening risk factor for COVID-19 patients, as reported previously [55]. This is supported by these finding; which found an increase in the risk of DVT and pulmonary embolism (PE) complications while comparing the frequency between admission and during admission, DVT (0.4% vs. 1.7%), and PE (4.9% vs. 6.9%), and found that 0.6% developed DVT and PE during admission. The number of patients who received low molecular weight heparin (LMWH), such as Fraxiparine and Enoxaparin, as deep venous thrombosis (DVT) prophylaxis during the first 12 hours, was higher in the survivor group than in the non-survivor group, in agreement with a previous report [56]. The incidence of COVID-19 major complications is another major concern. For cardiovascular complications, the study found stress cardiomyopathy in 44 (9.5%) patients, cardiac arrhythmia in 88 (18.9%), and myocarditis in 34 patients (7.3%), all of which were significantly higher in the non-survivor group. Stress cardiomyopathy has been reported to be associated with COVID-19, as its incidence has increased compared to that in the pre-pandemic era (7%–8% vs. 1%–2%) [57]. Cardiac arrhythmia has been reported in the literature as a complication of COVID-19, and a study in Wuhan in China found that 44.4% of critical care patients had arrhythmia [7]. Cardiac arrhythmia may have several underlying mechanisms, including electrophysiological or structural changes caused by viral infection [58]. Myocarditis has also been reported to be associated with COVID-19 infection, although a previous study reported that the majority of patients survived; however, higher proportion of patients with myocarditis in the non-survivor group [59]. Prolonged delirium as a neurocritical issue was observed in 75 (16.1%) of our patients and was significantly associated with the non-survivor group. Delirium as brain dysfunction has been reported to occur at a high rate in a previous international study; up to 54.9% of patients developed some type of delirium [60]. Seizure was observed in 42 (9%) cases and was found to be higher in the non-survivor group than in the survivor group. This was explained in the literature by several hypothesized mechanisms, such as the hypoxic effect of COVID-19 [61] and the ability of SARS-COV-2 to enter the central nervous system through the angiotensin-converting-enzyme-2 (ACE-2) receptor cells or through the olfactory nerve [62–64]. 5 Limitations This study has several limitations. First, the study was conducted in a country that has an ongoing civil war and severe financial crisis, which could explain the high mortality rate, which may need to be addressed in other African countries with similar settings. Second, detailed information about ethnicity, race, and other laboratory findings was not provided due to the limited availability of resources, which limited our ability to analyze these factors on mortality. Another limitation is that the study was conducted in a single country, and the results should not be generalized to other parts of the world. Finally, center-to-center changes are possible; however, our study included the main hospitals and healthcare facilities for COVID-19 patients in Libya. 6 Conclusion This is a large multicenter prospective study of COVID-19 patients admitted to the ICU in Libya with detailed laboratory, clinical, and SOFA score characteristics, complications, and mortality up to 60 days following ICU admission as the evaluated outcome. Our study also examined the differences in survival rates and factors contributing to mortality. The study showed one of the highest mortality rates in a country with limited resource settings and ongoing conflict, along with a financial crisis. Several challenges for both the healthcare system and healthcare workers in Libya underlie these results. However, this requires urgent intervention from policymakers and health authorities by financially supporting healthcare centers, providing sufficient mechanical ventilation devices, increasing ICU bed capacity, and increasing the nurse and specialist ratio to patients. In addition, there is an urgent need to support healthcare workers by providing financial compensation, psychological support, and adequate training. Libyan patients with COVID-19 in critical care are at a higher risk of mortality and complications than those in other regions of the world. Several factors were found to be predictive of mortality, which should be considered by clinicians and healthcare policy makers in future planning and support for hospitals and healthcare workers during the ongoing COVID-19 pandemic. Acknowledgments Declarations: All methods were carried out in accordance with relevant guidelines and regulations. TI - Epidemiology, outcomes, and utilization of intensive care unit resources for critically ill COVID-19 patients in Libya: A prospective multi-center cohort study JF - PLoS ONE DO - 10.1371/journal.pone.0251085 DA - 2021-04-30 UR - https://www.deepdyve.com/lp/public-library-of-science-plos-journal/epidemiology-outcomes-and-utilization-of-intensive-care-unit-resources-axJgOF9eum SP - e0251085 VL - 16 IS - 4 DP - DeepDyve ER -