Skip to main content

A blood microRNA classifier for the prediction of ICU mortality in COVID-19 patients: a multicenter validation study



The identification of critically ill COVID-19 patients at risk of fatal outcomes remains a challenge. Here, we first validated candidate microRNAs (miRNAs) as biomarkers for clinical decision-making in critically ill patients. Second, we constructed a blood miRNA classifier for the early prediction of adverse outcomes in the ICU.


This was a multicenter, observational and retrospective/prospective study including 503 critically ill patients admitted to the ICU from 19 hospitals. qPCR assays were performed in plasma samples collected within the first 48 h upon admission. A 16-miRNA panel was designed based on recently published data from our group.


Nine miRNAs were validated as biomarkers of all-cause in-ICU mortality in the independent cohort of critically ill patients (FDR < 0.05). Cox regression analysis revealed that low expression levels of eight miRNAs were associated with a higher risk of death (HR from 1.56 to 2.61). LASSO regression for variable selection was used to construct a miRNA classifier. A 4-blood miRNA signature composed of miR-16-5p, miR-192-5p, miR-323a-3p and miR-451a predicts the risk of all-cause in-ICU mortality (HR 2.5). Kaplan‒Meier analysis confirmed these findings. The miRNA signature provides a significant increase in the prognostic capacity of conventional scores, APACHE-II (C-index 0.71, DeLong test p-value 0.055) and SOFA (C-index 0.67, DeLong test p-value 0.001), and a risk model based on clinical predictors (C-index 0.74, DeLong test-p-value 0.035). For 28-day and 90-day mortality, the classifier also improved the prognostic value of APACHE-II, SOFA and the clinical model. The association between the classifier and mortality persisted even after multivariable adjustment. The functional analysis reported biological pathways involved in SARS-CoV infection and inflammatory, fibrotic and transcriptional pathways.


A blood miRNA classifier improves the early prediction of fatal outcomes in critically ill COVID-19 patients.


The COVID-19 pandemic has exerted dramatic pressure on the health care system globally, especially in intensive care units (ICUs). COVID-19 is characterized by substantial heterogeneity in progression rates [1], posing a considerable challenge to triage critically ill patients, inform early intervention and guide both ICU capacity and resource allocation. In this scenario, reliable biomarkers would be invaluable to improve risk stratification and allow for more effective clinical decision-making. These biomarkers may also be legitimate targets for therapeutic intervention to prevent ICU complications and adverse outcomes.

MicroRNAs (miRNAs) are small noncoding RNAs involved in posttranscriptional gene regulation that have emerged as innovative biomarkers for a number of conditions [2,3,4], including viral respiratory infections [5] and the management of critical patients [6]. Circulating miRNAs are cost-effective biomarkers that can be quantified through a minimally invasive blood draw using techniques already employed in clinical laboratories [7]. Furthermore, recent clinical trials have demonstrated that miRNA-based therapies seem to be well tolerated and show promising effects [8].

Here, we used samples from a large multicenter cohort of critically ill patients with COVID-19, i.e., the CIBERESUCICOVID study (NCT04457505), to construct a blood miRNA classifier that could be used to predict all-cause in-ICU mortality. We focused our attention on a 16-miRNA panel previously associated by our group with the severity of the disease and its adverse clinical outcomes in a miRNA biomarker discovery study [9]. To the best of our knowledge, our study is the largest on miRNAs as biomarkers for the clinical management of critically ill COVID-19 patients.

Patients and methods

Study design and data collection

This is a substudy of the CIBERESUCICOVID study registered at with the identification NCT04457505. CIBERESUCICOVID is a multicenter, observational, prospective/retrospective cohort study that enrolled critically ill COVID-19 patients admitted to the ICUs of 55 Spanish hospitals [10]. CIBERESUCICOVID started in May 2020 by collecting the retrospective data of patients admitted to participating ICUs (from February 2020) and continued prospectively until February 2021. After enrollment, comprehensive demographic, clinical, pharmacological and laboratory data were exhaustively collected at hospital and ICU admission, as previously described in Torres et al. [10]. The pharmacologic treatments administered and interventions performed during hospital admission until either discharge from hospital or death were also collected. Definitions have also been previously published [11]. Deidentified patient data were abstracted manually from the electronic medical records and stored in a Research Electronic Data Capture (REDCap) database hosted in the Centro de Investigación Biomédica en Red (CIBER, Madrid, Spain). Data from patients’ medical records were incorporated into the database by trained local researchers. The study coordinators ensured the integrity and timely completion of data collection. Prior to statistical analyses, incoherent or missing data were checked by independent experienced data collectors trained in critical care.

The study protocol was approved by the respective ethics committee of each participating hospital. The study was designed and conducted in compliance with the Declaration of Helsinki and national and international law on data protection. Participants, or their legal representatives, provided informed consent, when possible, for the use of the samples and data. In the remaining cases, an informed consent waiver was authorized by the ethics committee.

Primary and secondary outcomes

The primary outcome was all-cause in-ICU mortality. Secondary outcomes included all-cause 28-day mortality and all-cause 90-day mortality from ICU admission and length of hospital stay, ICU stay and invasive mechanical ventilation (IMV).

Study sample

Patients admitted to the ICU at the participating hospitals were enrolled in the current substudy if they fulfilled the following inclusion criteria: age over 18, laboratory-confirmed SARS-CoV-2 infection according to a standardized test, admission to the ICU and blood sample collected during the first 48 h available. Patients were excluded if they had unconfirmed SARS-CoV-2 infection, lacked data at baseline or hospital discharge or were admitted to an ICU for other causes. A list of participating hospitals is provided in Supplemental Figure S1 (Additional File 2).

Sample size calculation

The sample size calculation was based on comparisons of miRNA levels between study groups (survivors vs. nonsurvivors) using a two-sample t test. Relevant biological differences were considered as a fold change of 1.2 (or 0.83 in downregulated miRNAs). According to our previous data [9], we assumed a coefficient of variation of 0.5. The significance level was fixed at 0.05. A minimum sample size of 143 per group was necessary to achieve 90% statistical power. Expecting a mortality rate of 33%, a total sample size of 429 samples was necessary for the validation study. The sample size used was higher to ensure the necessary number of patients. Ultimately, 503 patients from the CIBERESUCICOVID consortium with blood samples available were included.

Experimental methods

Details regarding miRNA quantification and prediction of miRNA target regulation are provided in the Supplemental Methods (Additional File 1).

Statistical analysis

The characteristics of the study population were summarized by descriptive statistics. Data are presented as the medians [P25; P75] for continuous variables and as frequencies (percentage) for categorical variables. Continuous variables were compared using the Mann‒Whitney U test. Categorical variables were compared using the Fisher’s exact test. Linear models with Empirical Bayes statistic were used to evaluate differences in miRNA levels between survivors and nonsurvivors [12]. miRNAs with a significant difference [false discovery rate (FDR) < 0.05] between study groups and a fold change (FC) higher than 1.2 (or lower than 0.83 for downregulated levels) were considered differentially detected. Adjusted models were evaluated, including potential confounding factors that could affect the association of the miRNA and the outcome. The Pearson correlation coefficient was used to assess the correlation between continuous variables. Correlations between validated miRNAs and laboratory parameters were performed in the whole population. Correlations between the miRNAs included in the 4-blood miRNA classifier and length of hospital stay, ICU stay or IMV were performed in survivors of ICU stay. A blood miRNA classifier for all-cause in-ICU mortality was constructed using a relaxed least absolute shrinkage and selection operator (LASSO) model. This approach was used to reduce the collinearity of the multivariable model due the high correlation between the plasma levels of miRNAs. Fivefold cross-validation was carried out to determine the lambda parameter of the LASSO model. Lambda and gamma parameters were selected as the values associated with one standard error greater than the minimum mean square error (MSE). The miRNA levels were standardized prior to fitting the LASSO model. For levels of selected miRNAs, a cutoff point was established for fitted mortality risk using a maximally selected log-rank statistic [13]. For the blood miRNA classifier, in order to improve the clinical interpretation, patients were classified in high- or low-risk based on their probability of mortality, estimated by the final model, using the same method. The hazard ratio (HR) was estimated using Cox regression models including the dichotomized levels of individual miRNAs and the blood miRNA classifier (high- and low-risk groups) [13]. Kaplan‒Meier curves were used to illustrate differences among groups in the time-to-event outcome, and the log-rank test was performed to assess statistical significance. The clinical model for fatal outcomes included demographic characteristics, baseline comorbidities, laboratory tests, organ support and oxygenation variables: age (years), sex (female vs. male), hypertension (yes vs. no), chronic respiratory disease (yes vs. no), chronic kidney disease (yes vs. no), diabetes (yes vs. no) and at ICU admission: PaO2/FiO2 ratio, pH, lymphocyte count (x109/L), platelet count (x109/L), D-dimer (ng/mL), urea (mg/dL) and serum creatinine (mg/dL). The variables were included based on clinical relevance and bibliography [11]. Logistic regression models were used to analyze the association between the miRNA classifier and 28-day and 90-day mortality. The C-index and area under the ROC curve (AUC) were calculated for survival and logistic models, respectively. The incremental gain in model discrimination was evaluated using DeLong’s test. The p-value threshold defining statistical significance was set at < 0.05. All statistical analyses were performed using R software, version 4.0.2.


Study sample characteristics

The main baseline characteristics of the study cohort are summarized in Table 1. The substudy included patients admitted from March 2020 to February 2021. The median age was 65.0 [56.0;73.0] years, and 31.6% were females. Of the total subcohort, 29.3% died during the ICU stay. The 28-day and 90-day mortality rates were 21.4% and 28.9%, respectively. Nonsurvivors were older and had a higher prevalence of comorbidities, including hypertension, diabetes mellitus and chronic kidney disease. At ICU admission, this group had significantly higher APACHE-II and SOFA scores. At the same time, nonsurvivors had a lower PaO2/FiO2 and a higher PaCO2 than survivors. Laboratory parameters highlighted disparities between survivors and nonsurvivors in neutrophil counts, lymphocyte counts, creatinine, urea and D-dimer. Patients who died during the ICU stay required a higher use of IMV and prone positioning.

Table 1 Characteristics of study sample

Differential microRNA profiles among ICU survivors and nonsurvivors: validation study

First, we interrogated our candidate miRNAs in an independent cohort of critically ill COVID-19 patients (Supplemental Table S1; Additional File 2). Linear models with the empirical Bayes statistic were used to identify differentially expressed miRNAs between survivors and nonsurvivors. From the sixteen candidates, six miRNAs (miR-27a-3p, miR-27b-3p, miR-148a-3p, miR-199a-5p, miR-214-3p and miR-491-5p) showed contradictory results when compared with previous findings (Supplemental Figure S2; Additional File 2). Furthermore, the previous association between miR-150-5p and COVID-19 was not observed in the current study (FDR 0.102) (Fig. 1A). Therefore, these seven miRNAs were discarded from subsequent analysis. Ultimately, nine miRNAs reached a statistically significant signal (FC < 0.83, FDR < 0.05) and showed similar results to those in the previous study in terms of the size effect and direction of the association: miR-16-5p, miR-92a-3p, miR-93-5p, miR-98-5p, miR-132-3p miR-192-5p, miR-323a-3p, miR-451a and miR-486-5p (Fig. 1A). No great impact of confounding factors was observed in the association between miRNA levels and all-cause in-ICU mortality (Supplemental Table S2; Additional File 2). Poor correlations were observed between validated miRNAs and laboratory parameters (rho < 0.3) (Supplemental Figure S3; Additional File 2). The nine validated miRNAs entered the next phase.

Fig. 1
figure 1

Construction of the 4-blood microRNA classifier. (A) Violin plot comparing validated microRNA levels between survivors and nonsurvivors to ICU stay. Between-group differences were analyzed using linear models for arrays. P-values describe the significance level for each comparison; (B) hazard ratio (HR) and 95% confidence interval (90% CI) for each validated microRNA. Cox regression models include the dichotomized levels of individual miRNAs; (C) Blood microRNA classifier constructed using a selection process based on relaxed least absolute shrinkage and selection operator (LASSO) model; (D) Kaplan‒Meier estimations for the 4-blood microRNA classifier. (E-G) Combination of the 4-blood microRNA classifier with established clinical predictors or contemporaneous prognostic scores. (E) Clinical model (n = 373); (F) APACHE-II (n = 352); (G) SOFA (n = 379). The graph displays the hazard ratio (HR) and 95% confidence interval (90% CI) for each variable. Cox regression models included the dichotomized levels of the 4-blood microRNA classifier. The hazard ratio (95% CI) is displayed as a 1-SD change for continuous predictors

A blood microRNA classifier for predicting ICU mortality risk in critically ill COVID-19 patients

Cox regression analysis revealed that low expression of eight miRNAs predicted the highest risk of all-cause in-ICU mortality (HR from 1.54 to 2.61) (Fig. 1B & Supplemental Table S3; Additional File 2). No association was observed for miR-92a-3p [HR 1.37, p-value 0.06]. We next constructed a blood miRNA classifier of all-cause in-ICU mortality using the LASSO algorithm for feature selection. A 4-blood miRNA signature consisting of miR-16-5p, miR-192-5p, miR-323a-3p and miR-451a was selected (Fig. 1C). Using the classifier, the group of patients at higher risk of in-ICU mortality showed an HR (95% CI) of 2.53 (1.81–3.54) compared to the low-risk group. The Kaplan–Meier curve revealed clear separation of survival between the high- and low-risk subgroups (log-rank χ2 30.38, p-value < 0.001) (Fig. 1D). The 4-blood miRNA classifier was also significantly associated with 28-day mortality [OR (95% CI) 3.50 (2.34–5.29), p-value < 0.001] and 90-day mortality [OR (95% CI) 3.65 (2.44–5.51), p-value < 0.001]. No significant correlations were found between the classifier and the length of hospital stay, ICU stay or IMV in critically ill survivors (Supplemental Figure S4; Additional File 2).

Then, we explored whether the 4-blood miRNA classifier cooperatively discriminates all-cause in-ICU mortality in combination with established clinical predictors or contemporaneous prognostic scores. To this end, we combined a model based on clinical predictors and the prognostic scores APACHE-II and SOFA with the 4-blood miRNA classifier (Fig. 1E & 1F & 1G). The addition of the miRNA signature to available clinical parameters improved risk discrimination: miRNA classifier-clinical model pair: C-index 0.74, DeLong test p-value 0.035; miRNA classifier-APACHE-II pair: C-index 0.71, DeLong test p-value 0.055; and miRNA classifier-SOFA pair: C-index 0.67, DeLong test p-value 0.001. Concerning 28-day and 90-day mortality, the addition of the miRNA classifier also improved the prognostic capacity of the clinical model, APACHE-II and SOFA (Supplemental Figure S5; Additional File 2).

Of note, the ability of the 4-blood miRNA classifier to predict mortality risk remained significant after extensive adjustment for clinical predictors, APACHE-II or SOFA (Fig. 1E & 1F & 1G). Indeed, we analyzed the relationship between the 4-blood miRNA classifier and established clinical features of the disease, including clinically available blood biomarkers. Again, poor correlations (rho < 0.3) were observed between the miRNA score with blood cell counts and biochemical parameters (Supplemental Figure S6; Additional File 2).

Functional analysis of the 4-blood microRNA classifier

The components of the classifier were subjected to Reactome, GO and KEGG analyses. The functional analysis reported 32 Reactome biological pathways, 36 GO terms and 9 KEGG molecular pathways. We identified downstream targets and mechanisms associated with viral infections, such as pathways implicated in transcriptional regulation and SARS-CoV infections, and processes implicated in COVID-19 physiopathology, including the NF-κB, VEGF and TGF signaling pathways. (Fig. 2).

Fig. 2
figure 2

Functional enrichment analysis of the 4-blood microRNA signature. The microRNAs that composed the classifier were included in the analyses. Graph representing the p-value versus the number of target genes for each microRNA using the Reactome (A), GO (B) and KEGG (C) databases. miRWalk2.0 (accessed date July 26th, 2022) was used to predict the interaction of microRNAs with their targets (TargetScan filter)


The actual prognosis of ICU mortality in COVID-19 patients constitutes a challenge because of the lack of risk assessment metrics [1]. Determining early predictors of mortality is mandatory to guide ICU capacity and resource allocation. To identify novel biomarkers, the plasma of critically ill COVID-19 patients obtained within the first 48 h of ICU admission was subjected to miRNA profiling.

In the current study, we first validated the potential use of nine miRNA candidates selected from a previous discovery study as advanced mortality predictors in the ICU. Our differentially expressed miRNAs quantified in an independent multicenter cohort show a high concordance with our own previous data [9]. Then, we explored the predictive potential of the validated candidates. Accordingly, we showed that eight miRNAs were associated with in-ICU survival. Furthermore, we constructed a 4-blood miRNA classifier that provides an innovative estimate of early in-ICU mortality prediction. The association between the classifier and mortality persisted even after multivariable adjustment that included clinical history and established biomarkers. In this sense, neither individual components of the 4-blood miRNA classifier nor the signature were correlated with prominent features of SARS-CoV-2 infection and disease severity, such as leukocyte counts, including neutrophil and lymphocyte counts, D-dimer or creatinine, among others, suggesting that the miRNA classifier may provide novel and independent information for prognostication. No correlation was observed between the classifier and length of hospital stay, ICU stay or IMV in critically ill survivors. The explanation of this result is unclear. It may be associated with differences in the study sample or the pathobiological mechanisms implicated in the different clinical outcomes. Finally, we demonstrated that the use of the 4-blood miRNA classifier in concert with conventional scores of adverse outcomes, i.e., APACHE-II and SOFA, and a model based on clinical predictors is more accurate than either alone for the prognostication of fatal outcomes. Indeed, the combination of the host miRNA classifier with the clinical predictors computes prognosis for all-cause in-ICU, 28-day and 90-day mortality more precisely than any information currently available at ICU admission.

In addition, to support the previous association reported by our group between miR-16-5p, miR-98-5p, miR-132-3p, miR-192-5p and miR-323a-3p and in-ICU mortality [9], these results are in line with findings from the literature. For instance, Fernández-Pato et al. recently demonstrated that plasma miR-98-5p is reduced in severe COVID-19 patients [14]. In a study from Wilson et al., both miR-323-3p and miR-451a were shown to be downregulated in plasma samples from individuals with severe disease [15]. Overall, our miRNA signature introduces a new horizon for a combinatorial clinical data–transcriptomic biomarker systems for mortality prognostication in critically ill COVID-19 populations. This advancement is especially relevant since only a few potential blood-based biomarkers, mainly inflammatory and thrombotic mediators [16,17,18], have been proposed for the prognostication of critically ill COVID-19 patients. Furthermore, recent findings suggest the superior specificity of miRNA for COVID-19 mortality compared to protein biomarkers [19]. Risk stratification based on progression rate can optimize triage, inform early intervention, improve the allocation of hospital resources and allow patient selection for clinical trials. Additional work will be necessary to determine whether the integration of electronic health record data with the host miRNA classifier, which can be quantified with relatively low cost through techniques already available in clinical laboratories such as qPCR, is suitable for use as a prognostic test in the ICU.

The development of novel therapies to improve outcomes in critically ill COVID-19 patients remains a challenge [20]. Rational therapeutic approaches could be inferred from the current results. Disease- and pathway-specific biomarkers that can predict patients’ evolution in the ICU may also constitute therapeutic targets. In this context, a second notable finding from our study is that the functional analysis of the miRNA classifier components has identified relative enrichment of pathways not only implicated in inflammatory mechanisms but also in VEGF, TGF and transcriptional regulation signaling, which are hallmarks of viral infection and disease severity [14, 21, 22]. Notably, a biological pathway related to SARS-CoV infection was enriched in the targets of the miRNA signature. The biological function reported in independent miRNA-based studies was also captured in our analysis. miR-16-5p and miR-98-5p have been predicted to target the SARS-CoV-2 genome and host factors that mediate viral infectivity [23, 24]. Interestingly, miR-16-5p is downregulated in macrophages exposed to SARS-CoV‐2 virion spike 1 glycoprotein (S1) [25]. The levels of the miRNA are upregulated in response to pro‐resolving mediators, which suggests a role of miR-16-5p in the resolution of inflammation and the return to homeostasis after viral stimuli. The downregulation of miR-451a levels has been described as a possible mechanism implicated in cytokine storms in COVID-19 patients [26]. Therefore, the association of these miRNAs as biomarkers for fatal outcomes may not be coincidental. Potential miRNA-based antiviral therapeutics warrant consideration [27].

Strengths and limitations

The strengths of the study are the use of a multicenter design including patients from 19 hospitals from different geographical locations, the use of a priori sample size calculation, the high number of patients analyzed, the rigorous control of data quality and the evaluation of the miRNA classifier in conjunction with electronic health history and contemporaneous tests. The study population is a subcohort of critically ill COVID-19 patients collected in the same time frame as the nationwide cohort CIBERESUCICOVID [10]; and therefore, it has similar characteristics being females (29.6%), with a median age of 63 [54–71] and hypertension (49.7%), obesity (36.3%) and diabetes mellitus (24.6%) as the most common comorbidities (Supplemental Table S4; Additional File 2). In addition, the subcohort is similar to populations described in other studies, such as Dongelmans et al. [28] and Carbonell et al. [29]. The comparison of the biomarker performance with clinical predictors is also fundamental to obtaining robust evidence for the potential clinical applicability of the miRNA classifier. The low number of miRNAs constitutes an advantage for its potential translation to patient management.

Some limitations must be acknowledged. First, although the study population seems to be representative of critically ill COVID-19 patients and the study design resembles clinical routine, which suggests that our findings may be generalizable, further validation in larger and international cohorts is encouraged. Second, possible confounders cannot be ruled out despite adjustment. Third, since we used a real-world clinical practice setting, the sociogeographical context should be considered. For instance, the impact of therapeutic effort limitation on the outcome, SARS-CoV-2 variants and modifications in the treatments during the pandemic cannot be ruled out. Fourth, the role of the miRNAs in the causal pathway of the disease cannot be inferred from our study design, although this knowledge is not a prerequisite for a biomarker. In addition, the role of circulating miRNAs as endocrine genetic signals remains to be fully explored [30]. Fifth, the statistical significance of functional analysis should be treated cautiously since the analysis was based on selected miRNAs. Further studies should explore whether miRNAs participate mechanistically in the pathophysiology of adverse outcomes. Sixth, the current results suggest the need for studies that are more comprehensive and less biased than the panel profiling used here.

To conclude, we validated the use of host miRNA profiles as a source of predictors of fatal outcomes in critically ill COVID-19 patients. In addition, we constructed a 4-blood miRNA classifier that stratifies patients according to their risk of fatal outcome at early stages of ICU admission. This classifier can be integrated into clinical information currently available to improve prognostication in the ICU and therefore to prospectively inform health care management and clinical decision-making.

Data Availability

The dataset that supports the findings of this study is available from the corresponding author on reasonable request.



area under the ROC curve


Centro de Investigación Biomédica en Red


fold change


false discovery rate


hazard ratio


Intensive Care Unit


Invasive Mechanical Ventilation


least absolute shrinkage and selection operator




mean square error


Research Electronic Data Capture


  1. Wynants L, van Calster B, Collins GS, Riley RD, Heinze G, Schuit E, et al. Prediction models for diagnosis and prognosis of covid-19: systematic review and critical appraisal. BMJ. 2020;369:m1328.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Magen I, Yacovzada NS, Yanowski E, Coenen-Stass A, Grosskreutz J, Lu CH, et al. Circulating miR-181 is a prognostic biomarker for amyotrophic lateral sclerosis. Nat Neurosci. 2021;24:1534–41.

    Article  CAS  PubMed  Google Scholar 

  3. Blanco-Domínguez R, Sánchez-Díaz R, de la Fuente H, Jiménez-Borreguero LJ, Matesanz-Marín A, Relaño M, et al. A Novel circulating MicroRNA for the detection of Acute Myocarditis. N Engl J Med. 2021;384:2014–27.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Raut JR, Schöttker B, Holleczek B, Guo F, Bhardwaj M, Miah K, et al. A microRNA panel compared to environmental and polygenic scores for colorectal cancer risk prediction. Nat Commun. 2021;12:4811.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Poore GD, Ko ER, Valente A, Henao R, Sumner K, Hong C, et al. A miRNA host response signature accurately discriminates Acute respiratory infection etiologies. Front Microbiol. 2018;9:2957.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Molinero M, Benítez ID, González J, Gort-Paniello C, Moncusí-Moix A, Rodríguez-Jara F, et al. Bronchial aspirate-based profiling identifies MicroRNA Signatures Associated with COVID-19 and Fatal Disease in critically ill patients. Front Med (Lausanne). 2022;8:756517.

    Article  PubMed  Google Scholar 

  7. Pinilla L, Barbé F, de Gonzalo-Calvo D. MicroRNAs to guide medical decision-making in obstructive sleep apnea: a review. Sleep Med Rev. 2021;59:101458.

    Article  CAS  PubMed  Google Scholar 

  8. Täubel J, Hauke W, Rump S, Viereck J, Batkai S, Poetzsch J, et al. Novel antisense therapy targeting microRNA-132 in patients with heart failure: results of a first-in-human phase 1b randomized, double-blind, placebo-controlled study. Eur Heart J. 2021;42:178–88.

    Article  PubMed  Google Scholar 

  9. de Gonzalo-Calvo D, Benítez ID, Pinilla L, Carratalá A, Moncusí-Moix A, Gort-Paniello C, et al. Circulating microRNA profiles predict the severity of COVID-19 in hospitalized patients. Transl Res. 2021;236:147–59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Torres A, Motos A, Ceccato A, Bermejo-Martin J, de Gonzalo-Calvo D, Pérez R, et al. Methodology of a large Multicenter Observational Study of patients with COVID-19 in spanish intensive care units. Arch Bronconeumol. 2022;58(Suppl 1):22–31.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Torres A, Motos A, Riera J, Fernández-Barat L, Ceccato A, Pérez-Arnal R, et al. The evolution of the ventilatory ratio is a prognostic factor in mechanically ventilated COVID-19 ARDS patients. Crit Care. 2021;25:331.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43:e47–7.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Hothorn T, Lausen B. On the exact distribution of maximally selected rank statistics. Comput Stat Data Anal. 2003;43:121–37.

    Article  Google Scholar 

  14. Fernández-Pato A, Virseda-Berdices A, Resino S, Ryan P, Martínez-González O, Pérez-García F, et al. Plasma miRNA profile at COVID-19 onset predicts severity status and mortality. Emerg Microbes Infect. 2022;11:676–88.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Wilson JC, Kealy D, James SR, Plowman T, Newling K, Jagger C, et al. Integrated miRNA/cytokine/chemokine profiling reveals severity-associated step changes and principal correlates of fatality in COVID-19. iScience. 2022;25:103672.

    Article  CAS  PubMed  Google Scholar 

  16. de Bruin S, Bos LD, van Roon MA, Tuip-de Boer AM, Schuurman AR, Koel-Simmelinck MJA, et al. Clinical features and prognostic factors in Covid-19: a prospective cohort study. EBioMedicine. 2021;67:103378.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Gorog DA, Storey RF, Gurbel PA, Tantry US, Berger JS, Chan MY, et al. Current and novel biomarkers of thrombotic risk in COVID-19: a Consensus Statement from the International COVID-19 thrombosis biomarkers Colloquium. Nat Rev Cardiol. 2022;19:475–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Battaglini D, Lopes-Pacheco M, Castro-Faria-Neto HC, Pelosi P, Rocco PRM. Laboratory biomarkers for diagnosis and prognosis in COVID-19. Front Immunol. 2022;13:857573.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Gustafson D, Ngai M, Wu R, Hou H, Schoffel AC, Erice C, et al. Cardiovascular signatures of COVID-19 predict mortality and identify barrier stabilizing therapies. EBioMedicine. 2022;78:103982.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Hermine O, Mariette X, Porcher R, Resche-Rigon M, Tharaux P-L, Ravaud P. Effect of interleukin-6 receptor antagonists in critically ill adult patients with COVID-19 pneumonia: two randomised controlled trials of the CORIMUNO-19 collaborative group. Eur Respir J. 2022;60:2102523.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Evangelou K, Veroutis D, Paschalaki K, Foukas PG, Lagopati N, Dimitriou M, et al. Pulmonary infection by SARS-CoV-2 induces senescence accompanied by an inflammatory phenotype in severe COVID-19: possible implications for viral mutagenesis. Eur Respir J. 2022;60:2102951.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Giuliani A, Matacchione G, Ramini D, di Rosa M, Bonfigli AR, Sabbatinelli J, et al. Circulating miR-320b and mir-483-5p levels are associated with COVID-19 in-hospital mortality. Mech Ageing Dev. 2022;202:111636.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Kim WR, Park EG, Kang KW, Lee SM, Kim B, Kim HS. Expression analyses of MicroRNAs in Hamster lung tissues infected by SARS-CoV-2. Mol Cells. 2020;43:953–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Matarese A, Gambardella J, Sardu C, Santulli G. miR-98 regulates TMPRSS2 expression in human endothelial cells: key implications for COVID-19. Biomedicines. 2020;8:1–10.

    Article  Google Scholar 

  25. Recchiuti A, Patruno S, Mattoscio D, Isopi E, Pomilio A, Lamolinara A, et al. Resolvin D1 and D2 reduce SARS-CoV-2-induced inflammatory responses in cystic fibrosis macrophages. FASEB J. 2021;35:e21441.

    Article  CAS  PubMed  Google Scholar 

  26. Yang P, Zhao Y, Li J, Liu C, Zhu L, Zhang J, et al. Downregulated miR-451a as a feature of the plasma cfRNA landscape reveals regulatory networks of IL-6/IL-6R-associated cytokine storms in COVID-19 patients. Cell Mol Immunol. 2021;18:1064–6.

    Article  PubMed  Google Scholar 

  27. McDonald JT, Enguita FJ, Taylor D, Griffin RJ, Priebe W, Emmett MR, et al. Role of miR-2392 in driving SARS-CoV-2 infection. Cell Rep. 2021;37:109839.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Dongelmans DA, Termorshuizen F, Brinkman S, Bakhshi-Raiez F, Arbous MS, de Lange DW, et al. Characteristics and outcome of COVID-19 patients admitted to the ICU: a nationwide cohort study on the comparison between the first and the consecutive upsurges of the second wave of the COVID-19 pandemic in the Netherlands. Ann Intensive Care. 2022;12:5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Carbonell R, Urgelés S, Rodríguez A, Bodí M, Martín-Loeches I, Solé-Violán J, et al. Mortality comparison between the first and second/third waves among 3,795 critical COVID-19 patients with pneumonia admitted to the ICU: a multicentre retrospective cohort study. Lancet Reg Health-Eur. 2021;11:100243.

    PubMed  PubMed Central  Google Scholar 

  30. Devaux Y, Pinet F, de Gonzalo-Calvo D, Editorial. The non-coding transcriptome as a New Player in Intercellular Communication. Front Mol Biosci. 2022;9:121.

    Article  Google Scholar 

Download references


On behalf of the CIBERESUCICOVID Project (COV20/00110, ISCIII). The authors are indebted to María Arguimbau, Raquel Campo, Natalia Jarillo, Javier Muñoz and Manuel Sánchez (†) for their extensive support with project management and article preparation. This work was supported by IRBLleida Biobank (B.0000682) and “Plataforma Biobancos PT20/00021”. The human sample manipulation was performed in the Cell Culture Technical Scientific Service of the Universitat de Lleida (Lleida, Catalonia, Spain).


DdGC has received financial support from Instituto de Salud Carlos III (Miguel Servet 2020: CP20/00041), co-funded by the European Union. MM is the recipient of a predoctoral fellowship (PFIS 2021: FI21/00187) from Instituto de Salud Carlos III and co-funded by the European Union. MCGH is the recipient of a predoctoral fellowship from the University of Lleida. AC is supported by Instituto de Salud Carlos III (Sara Borrell 2021: CD21/00087), co-funded by the European Union. Financial support was provided by the Instituto de Salud Carlos III de Madrid (COV20/00110), co-funded by the European Union and Centro de Investigación Biomedica En Red – Enfermedades Respiratorias (CIBERES). CIBERES is an initiative of the Instituto de Salud Carlos III. Further supported by: Programa de donaciones “estar preparados”; UNESPA (Madrid, Spain); La Fundació La Marató de TV3, projecte amb codi 202108-30/-31.

Author information

Authors and Affiliations



Conception and design: IDB, AT, FB, DdGC. Data acquisition: All authors. Statistical analysis: IDB. Data analysis and interpretation: MM, MP, IDB, JFBM, AT, FB, DdGC. Manuscript drafting: MM, MP, DdGC. Review & editing: All authors. Final approval of the submitted version: All authors.

Corresponding author

Correspondence to Ferran Barbé.

Ethics declarations

Ethics approval and consent to participate

The study protocol was approved by the respective ethics committee of each participating hospital. The study was designed and conducted in compliance with the Declaration of Helsinki and national and international law on data protection.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Supplementary Material 2

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

de Gonzalo-Calvo, D., Molinero, M., Benítez, I.D. et al. A blood microRNA classifier for the prediction of ICU mortality in COVID-19 patients: a multicenter validation study. Respir Res 24, 159 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: