Introduction

Globally, healthcare workers encounter challenges in reducing coronavirus disease of 2019 (COVID‐19) severity and mortality1. Many subpopulations of patients with mild to non-severe COVID‐19 experience serious problems or even death, which is a growing concern1,2. According to reports, approximately 19% of COVID-19 patients experienced serious illness and 61.5% died within 28 days of admission, while 50% of hospitalized patients had no significant clinical and medical remission after 10 days3,4,5. Therefore, early diagnosis of patients with a possible serious COVID-19 infection and a high risk of mortality will relieve pressure on medical services, since treating a large number of patients places a significant burden on medical resources. The role of risk prediction is drastically shifting, and it helps to effectively determine if preventive protocols and treatment of positive cases are attempted. Therefore, early prognosis and care for this patient group are crucial to limiting disease progression and death1. In general, hematological predictors, risk factors, and potential complications of COVID-19 mortality and intensive care unit (ICU) referral need further investigation. Understanding the relevance of each risk factor to disease progression and mortality can assist in recognizing at-risk subpopulations and evaluating healthcare quality. Efforts should also be made to prepare for risk groups and estimate the risk of fatality in order to better understand the true patterns of mortality6. For example, age and gender have been identified as well-known risk factors for severe COVID19 outcomes: about 90% of COVID-19-related deaths in the United Kingdom (UK) have occurred in people over the age of 60, with 60% occurring in men7. Pre-existing conditions (such as tobacco use, coronary disease, hypertension, diabetes, respiratory and renal disease, and cancer) have also been related to an elevated risk of death6. For example, the Chinese Center for Disease Control and Prevention reported in a survey of 44,672 people (1023 deaths) that cardiovascular diseases, hypertension, diabetes, respiratory diseases, and cancers were associated with a high risk of death from COVID-198.

The aim of this meta-analysis was to assess the relevance of hematological/biochemical indices, pre-existing conditions, and complications to COVID-19 mortality and ICU referral and death, assuming an association between these factors. Furthermore, the association between hematological/biochemical indices, pre-existing conditions, and complications in COVID-19 non-survivors was investigated. The primary outcomes of this study were COVID-19 mortality and survival, and exposures included pre-existing conditions (such as age, gender, smoking, alcohol consumption, and so on) and comorbidities (i.e., hypertension, cerebrovascular diseases, diabetes, any comorbidities, cardiovascular diseases, renal diseases, chronic obstructive pulmonary disease (COPD), cancer, and liver diseases).

A network analysis was also used to estimate complex patterns of interaction9, to assess the specific structure of interrelationships between various factors, and to classify the degree of centrality and connectivity of patient characteristics/prognostic features to COVID-19 mortality.

Methods

Assessment of multicollinearity and correlation in influencing factors

The variance inflation factor (VIF) was used to determine whether any comorbidities or hematological and biochemical indices contribute to complications in COVID-19 patients. Indeed, VIF is used to calculate the amount of multicollinearity in variables, i.e., influencing factors. As a rough rule of thumb, VIF > 1 means that the variables are subject to collinearity, and VIF > 10 means that the correlation (collinearity) between the variables is very high. In addition, bivariate Pearson and Spearman correlations were used to assess linearity or nonlinearity between influencing factors (i.e., blood index and comorbidities) and complications. Correlation is significant at the 0.05 level (Bonferroni-corrected P value). Since the Anderson–Darling test confirmed that the data were normal, the Pearson correlation was used for network analysis.

Search strategy

A literature review was performed between December 2019 and April 2020 (no language restrictions were applied) in PubMed, Scopus, and various preprint servers (bioRxiv, medRxiv and SSRN) using MeSH keywords/terms, such as “COVID-19” AND “novel coronavirus” AND “new coronavirus” AND “coronavirus-2019” AND “COVID-2019” AND “SARS-COV-2” AND “2019-nCOV”.

Data extraction and analysis

The results of the search strategy were initially evaluated using abstracts and titles. Following that, the full text of the related articles was then checked using the inclusion and exclusion criteria. The final list of papers included was contrasted, and the disagreements were settled by a consensus discussion between authors. Three researchers (RK, AM, and AMR) independently extracted data such as the type and date of release, country, sample size, age, sex, blood indices and parameters, pre-existing health conditions, and complications. Using a structured spreadsheet, three authors (RK, AMR, and MS) tested the consistency of the data obtained. The Anderson–Darling test was used to determine if the resulting data was normal. For Meta-analysis and network analysis, META-MAR V2.7.0 and PAST were used, respectively. The Standardized Mean Difference (SMD) was used to identify the effect size of various hematological/biochemical indices, risk factors (pre-existing conditions), and complications in COVID-19 survivors and non-survivors, or in ICU-admitted patients (including both survivors and non-survivors) and non-ICU patients. Due to heterogeneity, random-effect models were used to calculate the weighted mean effect sizes and 95% confidence interval (CI).

Inclusion criteria

Figure 1 shows the flow diagram of studies assessed for inclusion. We collected data from research that recorded, at least, one of the following criteria, such as hematological/biochemical indices, pre-existing conditions, demographic factors, complications, and clinical outcomes (i.e., mortality or survival).

Figure 1
figure 1

PRISMA flow diagram of the studies identified, screened, reviewed, and included in the meta-analysis.

Pre-existing health conditions were described as conditions that existed prior to the patient’s arrival at the emergency department/hospital, and the decision was based on the medical record10. Complications (i.e., heart failure, septic shock, acidosis, respiratory failure, coagulopathy, acute renal injury, liver dysfunction, and secondary infection) were unanticipated events, diseases or symptoms triggered by the combination of a given disease11 and pre-existing conditions. Importantly, complications can occur even though appropriate precautions are taken. COVID-19 has been found to cause serious medical complications and death, especially in elderly patients and people with pre-existing health conditions.

Exclusion criteria

Data from over 200 reported clinical studies and preprints were screened. Following a comprehensive review of the data in figures and tables, 105 reports were excluded due to a lack of survivor or non-survivor sup-groups, examining infants/children/pediatrics, or having no DOI. This resulted in 97 eligible observational studies (19,014 patients)4,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107, with survivors (n = 14,359) and non-survivors (n = 4655) as strict subpopulations.

The reference list of papers and other systematic reviews were also scrutinized. The hematological and biochemical indices were checked to ensure they were the same, unless otherwise converted to the same unit. In terms of inclusion and exclusion criteria, paper selection resulted in different sample sizes for different parameters between survivors and non-survivors. These studies included patients from China, Italy, Scotland, the United States, UK, Japan, Singapore, South Korea, Iceland, Chile, the Netherlands, and Germany. Supplementary Table 1 describes the characteristics of the studies selected for meta-analysis.

Heterogeneity and risk of bias in individual studies

The I2 and Tau2 statistics were used to determine statistical heterogeneity108. According to the Cochrane's handbook for Systematic Reviews of Interventions, the I2 expresses the proportion of variance due to heterogeneity (i.e., 30% to 60%, 50% to 90% and 75% to 100% correspond to moderate, substantial and considerable degrees of heterogeneity, respectively)109. In addition, a meta-regression analysis was performed to assess the effect of variables on the effect size. The Z-test and its associated P-values determined whether the observed prevalence varied from zero percent.

Moreover, the Quality in Prognostic Studies (QUIPS) tool, a non-validated instrument with space for personal interpretation, was employed to assess the risk of bias in the studies110,111. The Cochrane Methods Prognosis group recommends the QUIPS tool for prognosis research because it tackles all common sources of bias110,111. Based on this, we concluded that the QUIPS method was sufficient for determining bias risk. Two team members (AM and RK) independently assessed the risk of bias in each study and classified it as low, moderate or high risk112. Any disagreements were resolved by consensus. The QUIPS tool includes domains on: study participation, study attrition, prognostic factor measurement, outcomes measurement, statistical analysis and reporting, and study confounding113.

Studies with a low risk of bias had thorough explanations of the population, design, and measures; as well as clear descriptions of how the measure was executed, equipments used, and how the data were interpreted. In studies with a medium risk of bias, there was some bias, but not enough to invalidate the data. These studies do not meet all of the requirements for a low risk of bias ranking, but no error is likely to result in significant bias. Studies with a high risk of bias had major flaws, indicating different forms of bias that might invalidate the findings. The high-risk study included one or more crucial or “fatal” flaws in its design, analysis, or reporting, as well as significant amounts of missing information112 The findings of the risk of bias assessment were identified in the narrative synthesis but were not included in the meta-analysis.

Network analysis

Data from variable, such as blood indices, pre-existing conditions (i.e., male, liver diseases, renal diseases, cerebrovascular diseases, diabetes, COPD, drinking, smoking, any comorbidities, cancer, cardiovascular diseases, hypertension), or complications (i.e., heart failure, respiratory failure, secondary infection, coagulopathy, acidosis, liver dysfunction, septic shock, acute cardiac injury, and acute kidney injury) were used for network analysis as the proportion of each parameter obtained from each study. In terms of parameters, such as age, time to hospital and BMI, the mean of each parameter obtained from each study was used for network analysis. Then data were classified in a binary manner for network analysis, with 0 and 1 representing mortality and survival, respectively. The PAST software (accessible at: http://folk.uio.no/ohammer/past) was used to implement the Circular and Fruchterman-Reingold algorithms as a force-directed layout algorithm. The Anderson–Darling test confirmed the normality of data, so the Pearson similarity index was used as a parametric index for network analysis. The Pearson correlation threshold of 50% (as a basic level) was used to establish the network of all variables. To be more specific, Pearson correlation thresholds of 72% and 93% were respectively selected to characterize the relationship between blood indices and outcomes (i.e., survivors or non-survivors, Fig. 2a–g). Pearson correlation thresholds of 79% and 97% were selected to determine the relationship between risk factors (pre-existing conditions) or complications and outcomes (i.e., survivors or non-survivors, Fig. 2f,g). These cutoff points were selected because there was no subsequent interaction among the variables after these thresholds. The node’s size represents its degree of connectivity, and the edges display the relationship between the two variables. The thicker edges indicate stronger correlations between variables. Nodes with more links are closer together. Small nodes and thin edges represent small values.

Figure 2
figure 2

Correlation-based network analysis. The Pearson correlation threshold of 50% was used to show the network of all variables (ac). More precisely, the Pearson correlation thresholds of 72% and 93% (d, e) were respectively selected to define the connection between survivors, non-survivors and blood parameters. The Pearson correlation thresholds of 79% and 97% were respectively chosen to assess the relationship between survivors, non-survivors and (f) risk factors or (g) complications. Circles of the network indicate the blood parameters (a, d, e), risk factors (b, f) and complications (c, g). The size of the node reflects the degree of connectivity of the node and the edges display the relationship between the two variables. The thicker edges reveal higher correlations between variables. Nodes with more links are close to each other. Network analysis and visualization was carried out using PAST and Fruchterman-Reingold algorithm or Circular algorithm as a force-directed layout algorithm. Abbreviations in panels (a), (d), and (e): Alb, albumin; HBG, hemoglobin; NEU, neutrophil; PLT, platelet; LYM, lymphocyte; WBC, white blood cells, PCT, procalcitonin; GGT, gamma-glutamyl transferase; CRP, C-reactive protein, CK, creatine kinase; Creat: creatinine, BUN, blood urea nitrogen; Bili, total bilirubin; AST, aspartate aminotransferase; ALT, alanine aminotransferase. Abbreviations in panels (b) and (f): BMI, body mass index; Time to H, time from symptoms appearance to hospitalization; Renal, renal disease; Cereb, cerebrovascular disease; Liver, liver disease; COPD, chronic obstructive pulmonary disease; Cardio, cardiovascular disease. Abbreviations in panels (c) and (g): Fail, failure; Res, respiratory; Liver, liver dysfunction; Sec, secondary; Kidney, acute kidney injury; Cardiac, acute kidney injury.

Results

Data characteristics and multicollinearity among the influencing factors

Table 1 presents the clinical outcomes, and Appendix 1 describes the data collection process for Table 1. The patient’s age (n = 9375) ranged from 25.3 to 80.0 years (49.8; CI95% [46.9, 52.7]). There were 5448 survivors (age: 46.6; CI95% [44.2, 48.9]) and 3927 non-survivors (age: 71.5; CI95% [66.4, 76.5]). Males outnumbered females in the non-survivor population (33.3 vs. 17.7%).

Table 1 COVID-19 patient characteristics, blood parameters, pre-existing conditions, and complications.

Both Pearson and Spearman analyzes were performed to examine linearity and nonlinearity between variables. Bivariate correlation analysis revealed that certain hematological and biochemical indicators were associated with complications in COVID-19 patients (see Supplementary Table 2).

Liver dysfunction had no linear relationship with blood index, but Spearman analysis showed a significant nonlinear relationship with gamma-glutamyltransferase (GGT) levels (r = 0.80), interleukin (IL)-6 (r = 0.80), c-reactive protein (CRP) levels (r = 0.70), creatinine (Cr) levels (r = 0.70), platelets (PLT) numbers (r = 0.70), alanine aminotransferase (ALT, r = 0.69), and higher white blood cells (WBC) counts (r = 0.68) (see Supplementary Table 2).

Using Pearson analysis, elevated gamma-glutamyltransferase (GGT) levels (r = 0.91) and decreased platelet count (PLT) (r = -0.75) were found to be significantly associated with respiratory failure. Interestingly, Spearman analysis revealed that respiratory failure was strongly associated with all blood indices except LYM, GGT, and albumin (see Supplementary Table 2).

Pearson analysis revealed that heart failure was significantly associated with higher GGT levels (r = 0.88), higher white blood cells (WBC) counts (r = 0.77), higher blood urea nitrogen (BUN) levels (r = 0.89), and lower PLT counts (r = − 0.83). Spearman found no significant non-linear relationship between heart failure and blood indices (see Supplementary Table 2).

Using Pearson analysis, we found that septic shock was significantly associated with higher Cr (r = 0.75) and lower PLT counts (r = − 0.76). However, Spearman analysis showed a positive association between septic shock and IL-6 (r = 0.90), BUN (r = 0.86), and creatine kinase (CK, r = 0.85) (see Supplementary Table 2).

Pearson correlation showed that acute cardiac injury was significantly correlated with higher GGT level (r = 0.75), higher Cr level (r = 0.70), and lower PLT counts (r = − 0.85). Spearman correlation identified an association between acute cardiac injury and BUN (r = 0.87) and CK (r = 0.87). Acute kidney injury was significantly correlated with higher GGT levels using Pearson (r = 0.77) and Spearman (r = 0.89) analysis (see Supplementary Table 2).

Pearson analysis showed that secondary infection was associated with higher neutrophil (NEU) counts (r = 0.76). Spearman analysis showed a trend (P = 0.06) to associate secondary infection with IL-6 (r = 0.76) and BUN (r = 0.76) (see Supplementary Table 2).

Additionally, we found associations with certain pre-existing health conditions and complications in COVID-19 patients (see Supplementary Table 3). According to the Pearson and Spearman analysis, there was no significant relationship between liver dysfunction and pre-existing health conditions in COVID-19 patients. Spearman analysis revealed that age (r = 0.77), comorbidities (r = 0.74), and diabetes (r = 0.69) were significantly associated with respiratory failure. Pearson analysis found no significant association between respiratory failure in COVID-19 patients and pre-existing health problems (see Supplementary Table 3).

Pearson analysis showed that age (r = 0.80), comorbidities (r = 0.78), cardiovascular disease (r = 0.79) and cerebrovascular disease (r = 0.75) were all significantly correlated with heart failure in COVID-19 patients. Spearman analysis showed that heart failure in COVID-19 patients was associated with cerebrovascular disease (r = 0.88), comorbidities (r = 0.82), diabetes (r = 0.81), cardiovascular disease (r = 0.71) and renal disease (r = 0.72) (see Supplementary Table 3).

Using Pearson analysis we found that septic shock was associated with cerebrovascular diseases (r = 0.73) in COVID-19 patients. Spearman analysis showed that septic shock was associated with diabetes (r = 0.77), comorbidities (r = 0.73), age (r = 0.70), hypertension (r = 0.70), and cardiovascular disease (r = 0.69) (see Supplementary Table 3).

Pearson analysis found no significant associations between acute cardiac injury, acute kidney injury, or secondary infections and pre-existing health conditions in COVID-19 patients. However, Spearman analysis detected a significant association between acute cardiac injury and COPD (r = 0.78), renal disease (r = 0.72), or cardiovascular disease (r = 0.70). Spearman analysis confirmed a significant association between renal disease (r = 0.70) or cardiovascular disease (r = 0.68) and COVID-19 complications, acute kidney injury. Using Spearman analysis, secondary infections were found to be associated with comorbidities (r = 0.83), cardiovascular disease (r = 0.88), cerebrovascular disease (r = 0.77), and diabetes (r = 0.70) (see Supplementary Table 3).

Furthermore, we found a very high VIF in blood indices as influencing variables for COVID-19 patient’s complications (see Supplementary Table 4). In addition, COPD and cerebrovascular disease were pre-existing health conditions that contributed to elevated collinearity (see Supplementary Table 5). This suggested that complications in COVID-19 patients could be caused by pre-existing health conditions or changes in blood parameters.

Heterogeneity and publication bias

Clinical heterogeneity can lead to statistical heterogeneity, and can be observed using techniques like the I2 index or meta-regression114. The I2 index was high (100%) across all meta-analysis (Figs. 3, 4 and 5). In an attempt to better explain heterogeneity, multivariate meta-regression analysis revealed substantial heterogeneity in the outcome of interest (i.e. mortality and survival), which may be attributed to heterogeneity in exposure sources, such as blood indices, pre-existing conditions, or complications (Table 2).

Figure 3
figure 3

Forest plot of blood parameters in survivors and non-survivors of COVID-19. The Standardized Mean Difference (SMD) and the 95% confidence intervals (CIs) were used to define the effect size of different blood indices in survivors and non-survivors. S, survivors; GGT, gamma-glutamyl transferase; NEU, neutrophil; WBC, white blood cell; CRP, C-reactive protein; AST, aspartate aminotransferase; CK, creatine kinase; IL-6, interleukin-6; BUN, blood urea nitrogen; ALT, alanine aminotransferase; PCT, procalcitonin; HBG, hemoglobin; PLT, platelet; LYM, lymphocyte; n, population size.

Figure 4
figure 4

Forest plot of pre-existing health conditions in survivors and non-survivors of COVID-19. The Standardized Mean Difference (SMD) and the 95% confidence intervals (CIs) were used to define the prevalence of various risk factors and complications for survivors and non-survivors of COVID-19. Time to hospital, time from symptoms appearance to hospitalization; Cerebrovascular, cerebrovascular disease; Cardiovascular, cardiovascular disease; Renal, renal disease; Liver, liver disease; BMI, body mass index; COPD, chronic obstructive pulmonary disease; S, survivors; n, population size.

Figure 5
figure 5

Forest plot of complications in survivors and non-survivors of COVID-19. The Standardized Mean Difference (SMD) and the 95% confidence intervals (CIs) were used to define the prevalence of various risk factors and complications for survivors and non-survivors of COVID-19. ARDS, acute respiratory distress syndrome; S, survivors; n, population size.

Table 2 The Multivariate meta-regression analysis.

Supplementary Table 6 shows the risk of bias assessment using QUIPS tool for all the observational studies used in the meta-analysis (n = 97). The majority of these 97 studies scored well (low risk of bias) in three of the QUIPS tool’s six domains: study attrition, prognostic factor measurement, and outcome measurement. The other three domains, namely study participation, study confounding and statistical analysis, were more problematic. According to the QUIPS tool, the majority of studies (74/97) used in this analysis had a moderate to high risk of bias.

Quantitative synthesis of data

As shown in Table 2 and Appendix 2, the multivariate meta-regression analysis revealed that the risk factors (exposures) (t, 4.77; CI95% [0.64, 1.68]; P < 0.001) were associated with the estimated intervention effect on COVID-19 mortality (outcome). Biochemical/hematological indices (t, 1.85; CI95% [− 0.11, 1.60]; P = 0.083) tended to be associated with the estimated intervention effect on COVID-19 mortality. Moreover, complications were linked to the estimated intervention effect on COVID-19 mortality (t, 3.80; CI95% [1.07, 4.36]; P = 0.005).

Meta-analysis of individual hematological indices

Figure 3 displays the individual Hedges’g for each blood parameter, along with the corresponding CI95%. Since there was substantial statistical heterogeneity (P < 0.001), a random-effect model was used to evaluate the effect sizes (see Appendix 2).

Non-survivor COVID-19 had higher number of NEU (2.81[2.70, 2.91]; Z = 53.97; P < 0.001) and WBC (2.38 [2.29, 2.47]; Z = 50.05; P < 0.001) counts, as well as higher GGT (4.10 [3.81, 4.39]; Z = 27.40; P < 0.001), Cr (2.40 [2.30, 2.49]; Z = 49.67; P < 0.001), c-reactive protein (CRP) (2.28 [2.19, 2.38]; Z = 46.23; P < 0.001), aspartate aminotransferase (AST) (1.44 [1.34, 1.54]; Z = 29.11; P < 0.001), CK (1.14 (1.03, 1.25); Z = 19.58; P < 0.001), interleukin (IL) − 6 (0.95 [0.82, 1.08]; Z = 14.04; P < 0.001), BUN (0.47 [0.38, 0.57]; Z = 9.62; P < 0.001), and bilirubin (0.20 [0.11, 0.29]; Z = 4.46; P < 0.001) levels.

Compared to survivors, non-survivors had a smaller number of lymphocyte (LYM) (− 1.74 [− 1.83, − 1.66]; Z = 41.36; P < 0.001) and PLT (− 1.55 [− 1.63, − 1.47]; Z = 36.89; P < 0.001), as well as lower hemoglobin (− 1.26 [− 1.35, − 1.17]; Z = 26.12; P < 0.001), albumin (− 0.80 [− 0.90, − 0.70]; Z = 15.50; P < 0.001), and procalcitonin (− 0.12 [− 0.20, − 0.03]; Z = 2.69; P = 0.007) levels.

COVID-19 mortality increases with age, hypertension, cerebrovascular disease, and diabetes

Figure 4 shows the meta-analysis forest plot in terms of pre-existing conditions based on the random effect model (see Appendix 2). Patient’s age (3.11 [3.05, 3.17]; Z = 100.70; P < 0.001); hypertension (2.30 (2.26, 2.35); Z = 100.00; P < 0.001), cerebrovascular disease (2.22 [2.13, 2.32]; Z = 45.95; P < 0.001), diabetes (2.11 [2.06, 2.15]; Z = 96.66; P < 0.001), any comorbidities (1.97 [1.99, 2.01]; Z = 84.99; P < 0.001), cardiovascular disease (1.55 [1.51, 1.59]; Z = 76.90; P < 0.001), COPD (1.16 [1.11, 1.20]; Z = 56.68; P < 0.001), renal disease (1.10 [1.06, 1.14]; Z = 52.59; P < 0.001), male sex (0.78 [0.75, 0.82]; Z = 44.59; P < 0.001), body mass index (BMI) (0.73 [0.46, 0.99]; Z = 5.38; P < 0.001), time from symptoms appearance to hospitalization (0.66 [0.61, 0.72]; Z = 23.17; P < 0.001), liver disease (0.52 [0.47, 0.56]; Z = 22.42; P < 0.001), cancer (0.45 [0.41, 0.48]; Z = 23.13; P < 0.001) and smoking history (0.13 [0.02, 0.24]; Z = 2.41; P = 0.016) was higher among non-survivors.

Meta-analysis identifies common complications among COVID-19 non-survivors

The prevalence of current complications was higher in COVID-19 non-survivors (2.71 [1.91, 3.51]; Z = 6.66; P < 0.001; I2 = 100.0%; Tau2 = 1.48) than survivors (see Fig. 5 and Appendix 2). Heart failure was the most common complication in COVID-19 non-survivors (7.40 [7.15, 7.64]; Z = 58.45; P < 0.001), followed by septic shock (4.49 [4.36, 4.63]; Z = 65.90; P < 0.001), acidosis (2.90 [2.64, 3.15]; Z = 22.24; P < 0.001), respiratory failure (2.80 [2.73, 2.87]; Z = 78.36; P < 0.001), acute cardiac injury (1.89 [1.83, 1.96]; Z = 54.84; P < 0.001), coagulopathy (1.79 [1.66, 1.93]; Z = 25.32; P < 0.001), acute kidney injury (1.64 [1.58, 1.69]; Z = 58.91; P < 0.001), secondary infection (1.31 [1.24, 1.37]; Z = 39.24; P < 0.001), and liver dysfunction (0.10 [0.01, 0.20]; Z = 2.08; P = 0.037).

Network analysis supports the results of the meta-analysis

At a cutoff point of 50%, the network correlation for blood indices (Fig. 2a), risk factors (Fig. 2b), and complication (Fig. 2c) was demonstrated. The number of PLT and LYM, as well as hemoglobin concentration, was associated with COVID-19 survivors at a maximum cutoff point of 72% where all parameters were disconnected after this cutoff point (Fig. 2d). At a maximum cutoff point of 93%, the number of NEU, GGT concentration, and the incidence of COVID-19 mortality were all related (Fig. 2e).

COVID-19 mortality was discovered to be linked to patient’s age, hypertension, cerebrovascular disease, diabetes, any comorbidities, cardiovascular disease (as pre-existing conditions, at a maximum cutoff point of 79%, Fig. 2f) and heart failure (as a complication, at a maximum cutoff point of 97%, Fig. 2g), according to the network analysis. These findings corroborated the results of the meta-analysis. The network analysis was able to map the association between all parameters, including mortality, blood indices, complications and pre-existing conditions, while the meta-analysis only ranked the potent factors involved in COVID-19 mortality.

Clinical outcomes and meta-analysis of COVID-19-infected ICU/non-ICU patients

The clinical outcomes of ICU-admitted COVID-19 patients and non-ICU COVID-19 patients are shown in Table 3. The average age of ICU survivors (n = 336) and non-survivors (n = 201) was 59.5 years (CI95%: 50.1, 68.9) and 70.7 years (CI95%: 65.5, 76.0), respectively. Non-ICU survivors (n = 2112) had an average age of 43.6 years (CI95%: 38.5, 48.7). Males outnumbered females in the ICU non-survivors (63.7%), ICU survivors (50.9%), and non-ICU survivors (52.3%). Diarrhea, nausea, coughing, and dyspnoea were more common in ICU-admitted patients than in non-ICU patients. In comparison to ICU non-survivors, a greater proportion of ICU survivors received antibiotics (83.5 vs. 62.5%), corticosteroids (83.0 vs. 54.5%), and intravenous immunoglobin (51.4 vs. 20.5%).

Table 3 Characteristics, blood parameters, symptoms, and pre-existing conditions among ICU and non-ICU patients infected with COVID-19.

Figure 6 shows the individual Hedges'g in terms of pre-existing conditions and blood parameters in ICU and non-ICU patients, as well as the corresponding CI95%. Because there was significant statistical heterogeneity (P < 0.001), the effect sizes were evaluated using a random-effect model. The prevalence of renal disease (0.72 [0.57, 0.86]; Z = 9.64; P < 0.001), COPD (0.29 [0.17, 0.41]; Z = 4.74; P < 0.001), and male percentage (0.23 [0.13, 0.32]; Z = 4.63; P < 0.001) were greater in COVID-19 patients admitted to ICU (Fig. 6a). Cancer was the most common in non-ICU patients (− 2.73 [− 2.88, − 2.57]; Z = − 34.55; P < 0.001), followed by age (− 1.75 [− 1.85, − 1.64]; Z = − 32.34; P < 0.001), and diabetes (− 0.97 [− 1.09, − 0.85]; Z = − 15.97; P < 0.001, Fig. 6a).

Figure 6
figure 6

Forest plot of (a) pre-existing health issues and (b) blood parameters in COVID-19-infected ICU and non-ICU patients. The Standardized Mean Difference (SMD) and the 95% confidence intervals (CIs) were used to define the effect size. n, population size; D, diseases; COPD, chronic obstructive pulmonary disease; S, survivors; NS, non-survivors; ALT, alanine aminotransferase; WBC, white blood cell; PCT, procalcitonin; BUN, blood urea nitrogen; NEU, neutrophil; AST, aspartate aminotransferase; CRP, C-reactive protein; PLT, platelet; LYM, lymphocyte; HBG, hemoglobin.

Platelet (0.63 [0.51, 0.76]; Z = 10.06; P < 0.001), WBC (0.50 [0.38, 0.62]; Z = 7.92 P < 0.001), and NEU counts (0.40 [0.27, 0.53]; Z = 6.03; P < 0.001) were observed to be linked with ICU referral (Fig. 6b). Creatinine levels (0.71 [0.58, 0.83]; Z = 10.97; P < 0.001), AST activity (0.51 [0.37, 0.64]; Z = 7.49; P < 0.001), albumin levels (0.47 [0.34, 0.61]; Z = 6.82; P < 0.001), and BUN levels (0.45 [0.32, 0.58]; Z = 6.75; P < 0.001) were also observed to be linked with ICU referral (Fig. 6b). It was found that CRP levels (− 0.69 [− 0.79, − 0.58]; Z = − 13.04; P < 0.001) and HBG levels (− 0.17 [− 0.33, − 0.01]; Z = − 2.10; P < 0.001) preferred non-ICU patients (Fig. 6b).

Age (− 3.83 [− 4.12, − 3.54]; Z = − 25.99; P < 0.001) favored survival among ICU-admitted patients, while cardiovascular disease (0.53 [0.34, 0.72]; Z = 5.55; P < 0.001) was frequent among ICU non-survivors (Fig. 7a). Compared to ICU survivors, ICU non-survivors had a greater number of WBC (− 1.36 [− 1.62, − 1.10]; Z = − 10.35; P < 0.001) and NEU (− 0.96 [− 1.22, − 0.70]; Z = − 7.26; P < 0.001). Also, ICU non-survivors had higher alanine aminotransferase (ALT, − 1.80 [− 2.09, − 1.51]; Z = − 12.25; P < 0.001), Cr (− 1.20 [− 1.45, − 0.94]; Z = − 9.29; P < 0.001), PCT (− 1.00 [− 1.25, − 0.76]; Z = − 7.97; P < 0.001), BUN (− 0.99 [− 1.25, − 0.73]; Z = − 7.46; P < 0.001), and AST (− 0.55 [− 0.79, − 0.30]; Z = − 4.33; P < 0.001) (Fig. 7b). However, lower levels of HBG (6.52 [5.74, 7.30]; Z = 16.30; P < 0.001), total bilirubin (0.76 [0.52, 1.00]; Z = 6.14; P < 0.001) and albumin (0.30 [0.05, 0.55]; Z = 2.38; P = 0.02), as well as LYM (0.93 [0.69, 1.18]; Z = 7.46; P < 0.001), favored mortality in ICU-admitted (ICU non-survivors; Fig. 7b).

Figure 7
figure 7

Forest plot of (a) pre-existing health conditions and (b) blood parameters in survivors and non-survivors of COVID-19 patients admitted to the intensive care unit (ICU). The Standardized Mean Difference (SMD) and the 95% confidence intervals (CIs) were used to define the effect size. n, population size; D, diseases; COPD, chronic obstructive pulmonary disease; S, survivors; NS, non-survivors; ALT, alanine aminotransferase; WBC, white blood cell; PCT, procalcitonin; BUN, blood urea nitrogen; NEU, neutrophil; AST, aspartate aminotransferase; CRP, C-reactive protein; PLT, platelet; LYM, lymphocyte; HBG, hemoglobin.

Discussion

According to the findings of this study, certain hematological/biochemical markers, pre-existing conditions, and complications were associated with mortality in COVID-19 patients, particularly in ICU-admitted patients, and should be taken into account for patient care. Furthermore, a multilevel approach based on network analysis and correlation analysis will aid in the deeper understanding of interactions within and between exposures and outcomes.

Our meta-analysis revealed that COVID-19 non-survivors, as well as ICU non-survivors, had higher NEU and greater LYM counts. According to the study of Qin et al115, lymphopenia (low LYM counts) and an increased NEU–LYM ratio were frequently observed in patients with severe COVID-19. This was also a more common characteristic in COVID-19-related death12. Inflammatory mediators, such as IL-2 and IL-6, can cause serious lymphopenia, resulting in LYM loss115. Qin et al.115 indicated that SARS-CoV-2 infection affects LYMs, resulting in secondary bacterial infections and an increased NEU count. Indeed, neutrophilia (NEU count > 7.5 × 109/L) has been linked to bacterial inflammation, cytokine storm, and hyper-inflammation116, all of which play significant pathogenetic roles in COVID-19 infection115,117. Consistent with previous studies32, we identified an increase in WBC and NEU counts and a decrease in LYM in COVID-19 non-survivors, particularly in ICU non-survivors. Therefore, changes in WBC, NEU, and LYM counts were associated with the risk of mortality in COVID-19 patients and patients admitted to the ICU.

In this meta-analysis, we found that GGT and AST concentrations were higher in COVID-19 non-survivors. Furthermore, ALT and AST levels were greater in ICU patients, as well as ICU non-survivors. Concentrations of ALT, AST, and GGT have been found to be markedly greater in dead patients than in recovered patients32,118. In a previous study119, GGT levels were shown to be elevated in COVID-19 patients. Higher GGT levels were associated with lower albumin and higher CRP, ALT, and ALP levels in the 82 COVID-19 patients who did not have chronic liver disease or an alcohol history. This elevation suggests that the liver is involved in COVID-19 patients120. Bernal-Monterde et al.119 reported that increased levels of GGT and ALP, as well as decreased albumin levels, were associated with increased risk of death in COVID-19. Indeed, viral infections that commonly affect the respiratory tract cause hypoxia119. In patients with pandemic H1N1 influenza infection, serum levels of ALT, AST, and GGT were found to be positively correlated with hypoxia121. These findings were consistent with our Pearson analysis, which found a strong relationship between GGT levels and respiratory failure (r = 0.91) or heart failure (r = 0.88), but not with liver dysfunction (r = 0.14). As mentioned above, heart failure was found to be the most common complications in non-survivors. Interestingly, Spearman correlation showed significant and non-linear associations between GGT and ALT and liver dysfunction. Increased ALT, AST, and GGT levels in COVID-19 patients, particularly in ICU non-survivors, appeared to be caused by heart failure-induced hypoxia, although further research is required to understand the details.

Since meta-analysis cannot establish relationships between variables, bivariate analysis, such as Spearman and Pearson correlation, and VIF were used to determine relationships and multicollinearity among influencing factors (i.e., blood indices and pre-existing conditions). Data showed that some complications found in COVID-19 patients, such as heart failure, were correlated with pre-existing conditions (i.e., age), as well as lower PLT and LYM counts or higher GGT levels. Moreover, network analysis was used to visualize the structure of relationships between factors affecting the COVID-19 outcome. Indeed, this approach will help in explaining the relationships between variables like blood indices, pre-existing conditions and complications, as well as the relationships between these variables and the outcome, i.e., mortality and survival. Network analysis identified a link between heart failure and increased mortality in COVID-19 patients. Therefore, combining a multi-level analysis with meta-analysis would help to achieve a better understanding of the relationships between patient characteristics and outcome.

This meta-analysis revealed a greater concentration of GGT, BUN, and Cr in COVID-19 non-survivors. Furthermore, COVID-19 patients admitted to the ICU had higher Cr levels and a greater prevalence of renal disease. Creatinine and BUN levels were also higher in died ICU patients. These findings suggest that SARS-CoV-2 has a clear impact on human kidneys122 and increases COVID-19 patient referral to the ICU, as well as ICU mortality. A study of 701 patients revealed that elevated serum Cr levels on admission were associated with severity due to severe coagulation pathway abnormalities123. Furthermore, increased urea levels had comparable, if not greater, impacts on hazard ratios. Another kidney failure marker is GGT124, which is a cell-surface enzyme that metabolizes extracellular glutathione, the primary antioxidant in mammalian cells125. A high level of GGT is often regarded as an early and marker of oxidative stress126 and it can be a source of reactive oxygen species in the presence of iron126,127. This, in turn, may result in renal vasoconstriction, salt retention, and subsequent kidney damage128. Abnormalities in the routine urine test performed on admission have been linked to disease progression and an increased risk of in-hospital death129. As a result, renal abnormalities on admission revealed a greater risk of deterioration, requiring proper triaging129; further research is needed.

Current evidence indicates that complications in COVID-19 patients may be caused by the virus’s direct effect, immune-mediated inflammation or drug-induced toxicity, assuming that the majority of patients were given high doses of antibiotics, antiviral drugs, and steroids118. Bivariate analyses showed that some complications in COVID-19 patients (e.g., heart failure, septic shock, respiratory failure, acute cardiac and kidney injury, and secondary infections) were mainly associated in a non-linear manner with both pre-existing conditions (such as age, comorbidities, diabetes, and cerebrovascular disease) and blood parameters (such as PLT, CK, and IL-6). However, it is unclear to what extent complications are exacerbated by COVID-19 infection. Zhou et al.12 found that sepsis was the most common complication, followed by respiratory failure, ARDS, and heart failure. In our meta-analysis, heart failure and septic shock were the most common complications diagnosed in dead patients. Sahu et al.130 found that COVID-19 patients who died from infection had a gradual increase in CRP levels. Li et al.131 suggested that direct viral disruption, hyper-inflammation, and hypoxemia may all contribute to cardiac injury. Serum CRP, as an inflammatory marker, has been linked to disease severity132, lung lesions133, acute kidney damage134, and cardiac injuries135 in COVID-19 patients. Furthermore, we found that CRP levels were greater in COVID-19 ICU patients than in non-ICU patients. Our findings imply that CRP may be a biomarker of ICU referral, as well as death in COVID-19 patients, emphasizing the need of regularly monitoring CRP changes.

According to the present meta-analysis, COVID-19 patients admitted to the ICU and COVID-19 non-survivors had a lower PLT counts, as well as lower hemoglobin and albumin concentrations. Our findings corroborated the previous study that showed a decrease in the number of PLTs in non-survivors but an increase in survivors136. Zhao et al.136 found that PLT count may dramatically reflect pathophysiological changes in COVID-19 patients, and an early decrease in PLT was associated with COVID-19 mortality. Viral infection appears to have damaged lung tissue, resulting in PLT activation, aggregation, and entrapment, which lead to thrombosis and increased PLT consumption136. PLTs have a short life cycle (8–10 days) and very few PLTs are preserved in bone marrow;137 it may be responsive to the severity of the patient’s conditions. Furthermore, viruses may cause a decrease in PLT production as a result of megakaryocyte infection, which may contribute to megakaryocyte apoptosis138. Therefore, PLT measurement may be beneficial in the care of COVID-19 patients and ICU referral, resulting in a much earlier and more effective prognosis.

Liu et al.139 reported that COVID-19 patients had the most consistent decreases in hemoglobin levels. The first case of COVID-19 in the United States revealed a minor decrease in hemoglobin on day 6 of illness30. Notably, patients with a composite outcome (i.e., ICU admission, invasive ventilation, and death) had lower hemoglobin levels18. The present meta-analysis identified a link between low hemoglobin levels in COVID-19 patients admitted to the ICU and ICU mortality (Figs. 6b and 7b). It is worth noting that the atypical type of ARDS in COVID-19 patients causes inadequate blood oxygenation and can be fatal140. Because hemoglobin concentration in the blood is one of the most significant factors of the oxygen-carrying capacity of the blood, hemoglobin levels in COVID-19 patients may play a key role in this respect140. Anemia in ICU patients has been linked to an increased risk of death. It has also been linked to a higher risk of complications such as acute kidney injury141. Non-survivors' hemoglobin levels upon admission and throughout the last ICU days are much lower than survivors'140. Similar to our findings, it was recently reported that patients with COVID-19 have considerably lower hemoglobin levels than patients that were not admitted to ICU142. Dinevari et al.143 reported that, the frequency of mortality (anemic: 23.9% vs. non-anemic: 13.8%) and ICU admission (anemic: 27.8% vs. non-anemic: 14.71%) were significantly greater in anemic COVID-19 patients than in non-anemic COVID-19 patients. It has been shown that inflammation caused by SARS-CoV-2 can disrupt erythropoiesis and decrease hemoglobin production. For example, IL-6 has been shown to be elevated in severe COVID-19 infection117 and disrupts hemoglobin production144. The present meta-analysis revealed that COVID-19 non-survivors had higher levels of IL-6. Our findings suggested that lower hemoglobin levels may be attributed to higher levels of IL-6, which requires further study in COVID-19 patients.

Age, hypertension, cerebrovascular disease, and diabetes have been found to be common risk factors among non-survivors in our meta-analysis. In accordance with the previous study12, we noticed that non-survivors were older than survivors (46.6 vs. 76.5 years), ICU survivors were older than non-ICU survivors (43.6 vs. 59.5 years), and ICU non-survivors were older than ICU survivors (70.7 vs. 59.5 years). Additionally, Spearman analysis showed that age had a linear relationship with heart failure but a nonlinear relationship with complications like respiratory failure, septic shock, and acute cardiac injury in COVID-19 patients. Furthermore, COVID-19 non-survivors had a greater proportion of hypertension and diabetes than survivors. ACE2 has been shown to be over-expressed in diabetic or hypertensive patient145,146. Diabetes and hypertension are treated with ACE inhibitors and angiotensin II type-I receptor blockers (ARBs), which causes an increase in ACE2 expression and infection with COVID-19146,147. According to previous research, the re-admission rate to the ICU for sepsis was twice as high in the ACE inhibitor (1.08%) and ARB (1.04%) groups than in the non-user group148. The findings of our Spearman analysis showed a significant correlation between diabetes and hypertension with septic shock as a complication of COVID-19 infection. Moreover, another study found a link between cerebrovascular disease and the risk of death in COVID-19 patients149, which was consistent with our findings. SARS-CoV-2 has been shown to have neuro-invasive properties and the ability to spread from the respiratory system to the central nervous system150. COVID-19 can also cause cerebrovascular complications as a result of inflammation, hypoxia, and diffuse intravascular coagulation151. Therefore, pre-existing conditions, such as age, cerebrovascular disease, diabetes and hypertension, may lead to a greater risk of death in COVID-19 patients. Furthermore, we found that cardiovascular disease, renal disease, age, and gender were connected to increased ICU referral and mortality. Thus, COVID-19 patients with these pre-existing conditions should be closely monitored.

We found that exposures, i.e., demographic factors (e.g., age, gender, smoking, and alcohol consumption), as well as pre-existing conditions or comorbidities, were the primary sources of heterogeneity in this study. This may be due to inconsistencies in study designs, large differences in sample size, and differences in study characteristics. In this study, we focused on a large particular subgroup (e.g., survivors or non-survivors, which included patients of various ages, genders, and pre-existing health conditions). This basically results in heterogeneity as confirmed by I2 index and multivariate meta-regression analysis. Moreover, the high heterogeneity in this meta-analysis may be explained by studies that reported either individual patient data or the mean for a cohort of patients. Other factors such as heterogeneity in survival group (which included mild to severe cases) and non-survival groups (who had various treatments) can also lead to publication bias152.

It should also be noted that the sample size of the studies used for this meta-analysis ranged from 1 to 3200 cases (mean: 132 [CI95%: 75.0, 189.0]). In terms of variability in blood sample collection, it was observed that most blood samples were collected around the time of admission The interval between the onset of symptoms and hospitalization ranged from 0 to 18 days, with an average of 7.0 (CI95%: 6.4, 7.5) days. As a result, obtaining blood samples and recording symptoms at various stages of the infection at the time of admission may obviously affect outcomes; this difference in data collection period may contribute to significant heterogeneity in samples. Furthermore, the period of follow-up in the studies used in this meta-analysis differed from 6 days to 4 weeks; such changes in follow-up times may impact disease development and, therefore, corresponding outcomes13.

According to the QUIPS assessment, the majority of the studies used had a moderate risk of bias. The majority of the studies included in this meta-analysis lacked data on the impact of blood parameters and pre-existing conditions on comorbidities, as well as the relationship between such comorbidities and mortality. QUIPS assessment suggested that future research should consider experiments with adequate statistical power and appropriate statistical methods to address the potential interrelationships between all prognostic factors, complications, and outcomes in COVID-19 patients.

In this study, we used a multi-level approach, including meta-analysis, bivariate analysis and network analysis, to establish potential associations between exposures (e.g., patient characteristics) and outcomes (e.g., mortality or survival). However, this meta-analysis has several limitations. We did not perform sensitivity and subgroup analyses, despite the inclusion of studies with patients at various stages of COVID-19. Moreover, the data were obtained from a variety of countries, including developed and developing nations, with varying levels of medical facilities, suggesting different management guidelines for related medical comorbidities.

In conclusion, various pre-existing conditions and biochemical/hematological indices were related with an increased risk of ICU referral and mortality in COVID-19 patients, particularly those admitted to the ICU. Also, the data showed that complications, such as heart failure and septic shock, were more common in COVID-19 non-survivors, which could be attributed to patient characteristics, emphasizing the importance of pre-screening at triage153.