Introduction

Since emerging in Wuhan, China, in late 2019, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes coronavirus disease 2019 (COVID-19), unfortunately has spread rapidly worldwide [1,2,3]. On March 11, 2020, the World Health Organization (WHO) declared COVID-19 a global pandemic (www.who.int) [4, 5]. So far, in all countries of the Western world, we are still fighting against this terrible infective disease even whether vaccines are now available.

The primary mechanism of transmission of SARS-CoV-2 is through infected respiratory droplets occurring by direct or indirect contact with nasal, conjunctival, or oral mucosa. Typical clinical manifestations of COVID-19 include fever, cough, myalgia, headache, and taste and smell dysfunctions [6]. Severe pneumonia as well as acute respiratory distress syndrome (ARDS) are pathological complications of COVID-19 [6]. Regarding the clinical adverse outcomes, including mortality, these are influenced by host factors, such as age, male sex, and comorbidities, which among the most significant include hypertension, diabetes, coronary heart disease, cerebrovascular disease, chronic obstructive pulmonary disease (COPD), and kidney diseases [6]. The immune response seems to have a dual effect: whilst in the initial phases it is pivotal and beneficial [7], in the later advanced stages it may be counterproductive [8], causing lymphopenia [9]. However, also factors directly related to virus, such as viral load kinetics, host immune-system defense mechanisms, and cross-reactive immune memory from previous exposure to other coronaviruses, may influence history of the disease that ranges from absence of symptoms to a mild infection, to a severe disease accompanied by high mortality [10]. While the underlying molecular event is still the focus of intense investigation [3, 11, 12], the monitoring [13, 14] and therapy [15,16,17,18] are still waiting innovative effective approach. Nevertheless, COVID-19 worse outcomes are all exacerbated by the SARS-CoV-2 ability to induce vessels thrombosis [19] leading to disseminated intravascular coagulation (DIC) and the cytokines storm [20]. In COVID-19 patients with high rate of mortality, a state of hyperfibrinolysis, characterized by increased fibrin degradation products, such as D-dimer and reduced platelets, is significantly present [21].

Fibrinolysis is a complex process whereby a fibrin-rich thrombus is degraded and remodeled by the protease plasmin, which belongs from the conversion of plasminogen (zymogen) underlying the action of plasminogen activators (tissue-type (tPA) and urokinase-type (uPA)) and plasminogen activator inhibitor-1 (PAI-1) [22]. A role in enhancing virulence and pathogenicity of plasmin(ogen) has been suggested for viruses containing a furin site in their envelope proteins, such as SARS-CoV-2 [23], since plasmin may be able to cleave the S protein of the virus, increasing its ability to bind with angiotensin converting enzyme 2 receptors of host cells, and facilitating virus entry and fusion [23]. Recently, a role for tPA and PAI-1 as biomarkers for COVID-19 has been also proposed since both found were elevated among COVID-19 hospitalized patients [24].

Based on this double role of plasmin(ogen) in facilitating both virus infection and in regulating fibrinolysis, which the deregulation is the most fatal complication of COVID-19, interest it has been focused on this mechanism, even if the association between levels of plasminogen and COVID-19 outcomes was merely speculative and inconclusive [20]. Only a study conducted in 20 patients with COVID-19 and 20 non-COVID-19 sick controls enrolled in the emergency department (ED) reported significant lower plasminogen plasma concentration in hospitalized COVID-19 patients, and in patients with COVID-19 who required intensive care unit (ICU), compared with those were discharged from the ED [25].

Therefore, in the present study, we aimed to analyze the association between plasminogen levels and the main COVID-19-related outcomes, including biological parameters, in a population of infected Caucasian patients.

Results

Characteristic of the study population

All characteristics of the study population are reported in Table 1. Results are reported as mean ± standard deviation and in percentage for categorical variables. The mean age of the whole study group was 69.8 ± 14.3 years, and they were mostly men (32, 58%). The overall mortality was 24% and the day of hospitalization 22.8 ± 17.8. As regards comorbidities, 24% had diabetes mellitus, 52% hypertension, 30% cardiovascular disease, 11% stroke, 15% COPD, and 11% were obese. Forty-nine patients (92%) were treated with hydroxychloroquine, 13% with tocilizumab/sarilumab antibodies, and 72% with antivirals.

Table 1 Demographic and clinical characteristics of the study population by level of plasminogen.

In the univariate analysis reported in Table 1, significant association was present between low levels of plasminogen and high rate of mortality (p < 0.001), higher levels of pancreatic amylase (p < 0.049), higher C reactive protein (CRP) levels (p < 0.002), high procalcitonin (PCT) (p < 0.019), and high fasting blood glucose (FBG) (p < 0.047). Moreover, significant associations were present for lower level of plasminogen and lower level of calcium (p < 0.014), higher levels of interleukin-6 (IL-6) (p < 0.017), and high levels of D-dimer (p < 0.011).

Low plasminogen levels are associated with higher inflammatory parameters in COVID-19 patients

Based on the rationale previously reported, and on the results shown in Table 1, we investigated further the association between plasminogen levels and the levels of the most relevant parameters associated to COVID-19 pathogenesis. As reported in the scatterplots in Fig. 1, lower levels of plasminogen resulted in a significant increase of CRP (r2 = −0.17, p = 0.003, Fig. 1A), FBG (r2 = −0.10, p = 0.001, Fig. 1B), and creatinine (r2 = −0.03, p = 0.005, Fig. 1C), a decrease in the glomerular filtration rate (GFR) (r2 = −0.03, p = 0.04, Fig. 1D), and a variation in coagulation parameters, including increase in INR (r2 = −0.06, p = 0.03, Fig. 1E) and in D-dimer (r2 = −0.06, p = 0.003, Fig. 1F), and significant decrease in activated partial thromboplastin time (APTT) (r2 = −0.001, p = 0.03, Fig. 1G). Even the association between plasminogen and levels of PCT, calcium, and IL-6 was inversely related, i.e., r2 = −0.405, p = 0.003; r2 = −0.333, p = 0.014; and r2 = −0.366, p = 0.006, respectively.

Fig. 1: Analysis of the correlation between level of plasminogen and markers related to COVID-9 pathogenesis.
figure 1

Levels of plasminogen were correlated with several markers linked with COVID-19 pathogenesis such as CRP (A), FBG (B), creatinine (C), GFR (D), PT/INR (E), D-D (F), and APTT (G). CRP C reactive protein, FBG fasting blood glucose, GFR glomerular filtration rate, PT/INR prothrombin time and international normalized ratio, D-D D-dimer, APTT activated partial thromboplastin time.

To evaluate the relationships between different inflammatory parameters and to take into account the novel findings that COVID-19-related inflammatory circuits formed by cytokines storm are mediated by an interaction between IL-6 and plasminogen mediators [26], we designed a descriptive multidimensional model including the correlation of the expression of coagulation and inflammatory parameters (Figs. 2 and 3). It was clear as plasminogen, IL-6, and APTT tended to cluster together compared to D-dimer and CRP, which tended to be dissociated, and therefore to possibly act independently on different pathways of proinflammation and DIC.

Fig. 2: Multidimensional scaling analysis reporting the correlation of the expression of coagulation and inflammatory parameters with plasminogen.
figure 2

IL-6 and APTT tended to cluster together with plasminogen compared to D-dimer and CRP, which tended to be dissociated, and therefore to possibly act independently on different pathways of proinflammation. D-D D-dimer, CRP C reactive protein, IL-6 interlukin-6, APTT activated partial thromboplastin time.

Fig. 3: Schematic representation of COVID-19 impact on subjects with lower levels of plasminogen.
figure 3

Subjects with plasminogen ≤ 80 showed poor prognosis following SARS-CoV-2 infection due to alteration in several inflammatory and coagulation-related parameters. CRP C reactive protein, FBG fasting blood glucose, IL-6 interlukin-6, GFR glomerular filtration rate. Created by BioRender.com.

Taken together, these data suggest a potential relationship between levels of inflammation, mediated by IL-6, levels of circulating plasminogen, as marker of fibrinolysis, and the shorter time in forming clots, as per APTT levels, typical of COVID-19 patients.

Low level of plasminogen is one of the most significant predictors of mortality in patients with COVID-19

Since the mortality in COVID-19 patients increases in logarithmic way and comorbidities and older age have a great impact on it [27], we aimed to evaluate whether factors, such as alteration in plasminogen levels, that generate the “blood clots storm” may be predictive markers.

Then, by using univariate logistic regression analysis, we tested the most significant parameters found associated with levels of plasminogen in univariate analysis (Table 2). COVID-19 patients with low levels of plasminogen showed a 12-fold significant increase in mortality compared to COVID-19 patients with normal or high levels of plasminogen (OR 12.57, 95% CI 2.46–64.0, β = 2.53, p = 0.002). By using the median values as cutoff, we also found that in patients with COVID-19, high level of CRP, as marker of inflammation, was associated with up to sevenfold risk of mortality (OR 7.353, 95% CI 1.29–41.7, β = 1.99, p = 0.02); high FBG, as marker of metabolic dysfunction, was greatly associated until up to 17-fold risk of mortality (OR 17.153, 95% CI 1.72–170.4, β = 2.84, p = 0.01); high creatinine, as marker of kidney failure, tended to be associated with up to fourfold risk of mortality (OR 4.707, 95% CI 0.95–23.1, β = 1.54, p = 0.05); and D-dimer, as marker of procoagulation, was associated up to sevenfold risk of mortality (OR 7.655, 95% CI 1.32–44.2, β = 2.03, p = 0.02) (Table 2).

Table 2 Risk of death in COVID-19 patients by level of selected inflammatory parameters (univariate logistic regression analysis).

The relative contribution to the risk of death for all parameters investigated was evaluated with a multiple logistic regression, in a model adjusted also for presence of comorbidities, age, and sex. As reported in Table 3, males and COVID-19 patients older than 70 years of age showed a borderline increase of the mortality risk compared to female and younger patients. In the univariate analysis, similar increases were found in patients with high CRP and D-dimer levels, while after multiple adjustments, the only parameter that remained significantly associated to mortality, with essentially the same risk, was the plasminogen (OR 12.68, 95% CI 1.69–95.0, β = 2.02, p = 0.01) (Table 3). Taken together these results suggest that in the complex pathway of COVID-19, levels of plasminogen may be indicated as independent predictor of mortality in these patients.

Table 3 Risk of death in COVID-19 patients (multiple logistic regression analysis).

Discussion

In the present study, in an aged population of patients with diagnosis of COVID-19 and with high rate of comorbidity, we demonstrated that low levels of plasminogen were significantly associated with several prognostic parameters of complications, such as inflammatory markers (CRP, PCT, and IL-6), markers of coagulation (D-dimer, INR, and APTT), and markers of organ dysfunctions (FBG and GFR). In the present population, by using a multidimensional model, we showed the tendency of lower levels of plasminogen to cluster with IL-6 and APTT, compared to other parameters of inflammation and coagulation, suggesting a possible pathway that links these three factors in triggering the disseminate intravasal coagulation, as complication of COVID-19. In addition, low levels of plasminogen were associated with 12-fold higher risk for mortality even after multiple adjustments, suggesting its power as independent predictor of worse outcomes in subjects with COVID-19.

Since the beginning it was clear that COVID-19 was a particular viral infection characterized by various complications, often fatal. One of the most important complication, leading to death, as we previously highlighted, is the ability of SARS-CoV-2 to induce disseminated procoagulations [20]. The mechanisms beyond this pathological process are not understood yet, and even in the era of vaccination, scientific committees have an open discussion worldwide on how to prevent and manage this virus-related side effect. Therefore, biomarkers that might predict this complication would be invaluable, especially in triaging patients on hospital admission and in frail and elderly subjects receiving vaccination.

Recently, among blood biomarkers measured from 187 patients admitted at hospital with diagnosis of COVID-19, also including IL-6, Krebs von den Lungen 6, troponin, ferritin, lactate dehydrogenase, B-type natriuretic peptide, PCT, and CRP, the strongest association with negative outcomes, defined as noninvasive ventilation or ICU admission, was the high levels of soluble urokinase plasminogen activator receptor (suPAR) and IL-6 [28]. These results align with our findings since high levels of suPAR lead to lower levels of circulating plasminogen for its degradation [29]. Moreover, in the present study, we found how plasminogen and IL-6 cluster together, along with APTT, in a multidimensional analysis, suggesting that cytokines storm, mainly mediated by IL-6, activates the process of fibrinolysis and the intrinsic and common pathways of the coagulation cascade. This link was already suggested in 91 patients with ARDS where among the proinflammatory cytokines, IL-6 was the most clinically suitable biomarker associated with high levels of PAI-1, and this pathway was beyond the cytokines storm activation in endothelial cells [26].

The role of levels of plasminogen activator and inhibitor in the association with severity of COVID-19 progression has been already reported [30]. However, here, we demonstrated the direct association between the consumption of plasminogen and increased markers of inflammation, coagulation, and organ dysfunctions. To the best of our knowledge, we are the first demonstrating the direct impact of plasminogen levels with COVID-19 outcomes and prognostic parameters. A previous study, in line with our data, demonstrated that lower levels of plasminogen resulted in highest rate of admission to ICU in COVID-19 patients but without reporting the association with other biomarkers of the disease [25].

This association may be explained since the shift in the procoagulant pathways has been observed in the later stages of COVID-19 and the switch of plasmin(ogen) can promote further infection, other than disseminate intravasal coagulation, by cleaving proteins that allow cell infection [31]. These would further justify the association we found in the present study between low levels of plasminogen and high rate of mortality. In univariate logistic regression analysis, we found low levels of plasminogen increased risk of death up to 12-fold, along with other prognostic markers, such as CRP, high FBG, high creatinine, and high D-dimer. These findings, excluding for plasminogen, further confirmed evidence from other studies [32]. However, after multiple logistic regression, which included in the model all comorbidities, the only parameter that maintained significance with the same level of strength was plasminogen. In agreement with our results, a previous study demonstrated as high PAI-1 tracked most closely with impaired oxygenation efficiency in COVID-19 patients, and tPA was the best predictor of death, as activator of plasminogen, then leading to its consume [24]. These data are of particular importance since understanding the mechanisms beyond the COVID-19 complications, especially those linked with the disseminate vascular coagulation that kills these patients, would be helpful both as biomarker of disease to prevent mortality and also to promote or adjust therapies modulating fibrinolysis, like the trial already running with alteplase (https://clinicaltrials.gov/ct2/show/NCT04640194?term=tpa&cond=Covid19&draw=2&rank=6).

We need to acknowledge some limitations for the present study over the current search for the underlying molecular mechanism of the disease [33,34,35]. Considering the emerging role of predictive medicine [36,37,38] and personalized therapy [39,40,41,42], this will soon become essential. Currently, we are far from this goal, as, for example, we do not know the role of Bcl-2 family [43, 44], degradation pathways [45, 46], autophagy [47,48,49,50,51], and, as described in this paper, plasminogen. First, the number of subjects enrolled in the study was small due to the research restrictions during the pandemic. Then, serial plasminogen measurements at different time-points will be also necessary in future studies to delineate this potential mechanism.

We would like to highlight that atomization inhalation of 10 mg plasminogen dissolved in 2-ml sterile water given twice daily in five clinically moderate patients with COVID-19 has been shown to improve the conditions of lung lesions, heart rate, and oxygen saturation, about 1 h after the first inhalation, suggesting that plasminogen may be effective and efficient in treating several complications during COVID-19 infections. The levels of plasminogen were not measured in the patients [52]. In this study, we did report a strong association between low levels of plasminogen and high rate of mortality in patients with COVID-19. Further studies are imperative to better understand this mechanism and to provide a timely identification of patients at higher risk of complications.

Material and methods

Study population

Fifty-five SARS-CoV-2-positive adults were enrolled in an open study by the Internal Medicine Unit, Departments of Systems Medicine and Experimental Medicine, University of Rome Tor Vergata from May to September 2020. The diagnosis was made according to the relevant diagnostic or classification criteria, and COVID-19 severity was classified according to WHO guidance (https://apps.who.int/iris/handle/10665/330893). Only patients with SARS-CoV-2 confirmed by polymerase chain reaction from nasopharyngeal swabs were included in the study. All demographic and clinical data were collected from medical charts and are reported in Table 1. The ethical committee of Tor Vergata University/Hospital approved this study (protocol no. 48.20, version 2020). Informed written consent was obtained from each patient. The study was done in accordance with the ethical principles of the Declaration of Helsinki and the Guidelines for Good Clinical Practice.

Blood parameters and plasminogen measurement

Blood samples were collected via routine blood draws for clinical indications at Emergency Department and Internal Medicine Unit of University of Rome Tor Vergata at the time of admission to Internal Medicine Unit, according to standard procedure for cellular [53,54,55] and clinical biochemistry [56] or redox analysis [56,57,58,59,60], as previously described [61, 62]. Plasminogen blood levels were measured at the same time. Within 3 h of collection, samples were centrifuged at 2000 × g for 15 min at 4 °C and subsequently frozen at −80 °C until analysis. Circulating (plasma) plasminogen activity was measured by performing an automated chromogenic assay (HemosIL, Werfen, Italy) and analyzed by using a Behring Nephelometer II System (Siemens Medical Solutions, Inc, Malvern, PA, USA) according to the manufacturer’s protocol and literature [63]. Based on the plasminogen kit protocols and reference range [63], we classified low level as the patients with values of plasminogen ≤ 80 (n = 19), and normal/high as those with values > 80 (n = 36).

Statistical analysis

Main characteristics of the study group were explored using descriptive statistics. Patients with levels of plasminogen below 80% were compared with those over this threshold, using Student’s t test and χ2 test for continuous and categorical variables, respectively. Bidimensional descriptive graphics were based on the Pearson’s r correlation coefficient, while the multidimensional plot was generated using the procedure Multidimensional Scaling Model by SPSS (IBM SPSS Statistics for Windows, Version 26.0), a procedure which helps understanding data clustering assigning observations to specific locations in a conceptual space. To provide a quantitative estimate of the risk of death associated to those variables, which from descriptive analyses resulted more associated with inflammatory parameters, a univariate logistic regression model was run. Confounding effect was taken subsequently into account including in the logistic model, an appropriate set of actual and potential confounders related to individual characteristics, lifestyle, and dietary patterns. For each model, the occurrence of overdispersion was checked by comparing the residual deviance with its degrees of freedom. The likelihood ratio test was applied to assess the significance of each variable in the logistic model. STATA (Stata Statistical Software: Release 11, StataCorp LP) software was used for modeling.