https://orcid.org/0000-0002-9549-567X
ABSTRACT
Background
MicroRNAs (miRNAs) have a crucial role in regulating immune response against infectious diseases, showing changes early in disease onset and before the detection of the pathogen. Thus, we aimed to analyze the plasma miRNA profile at COVID-19 onset to identify miRNAs as early prognostic biomarkers of severity and survival.
Methods and results
Plasma miRNome of 96 COVID-19 patients that developed asymptomatic/mild, moderate and severe disease was sequenced together with a group of healthy controls. Plasma immune-related biomarkers were also assessed. COVID-19 patients showed 200 significant differentially expressed (SDE) miRNAs concerning healthy controls, with upregulated putative targets of SARS-CoV-2, and inflammatory miRNAs. Among COVID-19 patients, 75 SDE miRNAs were observed in asymptomatic/mild compared to symptomatic patients, which were involved in platelet aggregation and cytokine pathways, among others. Moreover, 137 SDE miRNAs were identified between severe and moderate patients, where miRNAs targeting the SARS CoV-2 genome were the most strongly disrupted. Finally, we constructed a mortality predictive risk score (miRNA-MRS) with ten miRNAs. Patients with higher values had a higher risk of 90-days mortality (hazard ratio = 4.60; p-value < 0.001). Besides, the discriminant power of miRNA-MRS was significantly higher than the observed for age and gender (AUROC = 0.970 vs. 0.881; p = 0.042).
Conclusions
SARS-CoV-2 infection deeply disturbs the plasma miRNome from an early stage of COVID-19, making miRNAs highly valuable as early predictors of severity and mortality.
Introduction
The pathogenesis of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and the virus-host interactions during COVID-19 are still not entirely understood. Therefore, it is essential to continue unraveling the biological and molecular mechanisms behind the disparity in disease severity. Omics approaches are indispensable in this scenario, generating crucial information for decision-making in public health policies [Citation1] and allowing the identification of new biomarkers for predicting severity in COVID-19 patients.
MicroRNAs (miRNAs) are small non-coding RNAs that play a fundamental role in regulating gene expression by binding the 3' untranslated region of target RNAs. This post-transcriptional regulation is essential in both innate and adaptive immunity [Citation2], and the fine-tune of the immune response. Several miRNAs are key regulators of inflammation-related mediators, being essential in some inflammatory diseases [Citation2]. In addition, evidences point to an essential function of miRNAs in the pathogenesis and therapeutics of different viral diseases, including those related to the respiratory tract [Citation3], such as the respiratory syncytial virus (RSV) [Citation4], middle east respiratory syndrome coronavirus (MERS-CoV) [Citation5], SARS-CoV [Citation6], and SARS-CoV-2 [Citation7].
In this regard, several host miRNAs have been identified to target the SARS-CoV-2 genome, repressing its expression and mitigating the pathogenesis of COVID-19. Therefore, the interaction between SARS-CoV-2 and the cellular miRNome could modify the outcome of the infection, as some mutations in the 3´UTR region of the viral genome leads to a viral escape from the host immune system [Citation8]. SARS-CoV-2 infection induces a robust host miRNA response, as it was observed in patients admitted to intensive care units (ICU) [Citation6]. However, it is unknown the miRNA profile of other severity outcomes, specially at disease onset, which could give new insights to improve COVID-19 detection and management.
Therefore, miRNAs could have a potential role in counteracting the SARS-CoV-2 induction of inflammatory response [Citation9]. In this setting, miRNAs are deeply involved in regulating the expression of cytokines, chemokines, and growth factors [Citation10–12], which is directly related with the cytokine release syndrome (CRS) of COVID-19. To date, the plasma miRNome of COVID-19 patients of different severity grades has not been massively assessed. The integrative analysis of the miRNome data with plasma cytokine levels could unravel the miRNA-cytokine interplay in COVID-19, identifying potential miRNAs as new therapeutic approaches against CRS and putative survival biomarkers at disease onset.
This study aimed to characterize the plasma miRNome at the onset of COVID-19 of different severity statuses and to identify a miRNA signature of mortality.
Methods
Design and study population
A multicenter observational study was performed in 96 patients with SARS-CoV-2 infection enrolled from March to August 2020 at three public hospitals from Madrid: University Hospital Infanta Leonor, Del Tajo University Hospital, and Príncipe de Asturias University Hospital. The study protocol was approved by the Ethics Committee of the Institute of Health Carlos III (PI 33_2020-v3) and the Ethics Committee of each hospital. Written, oral, or delegated informed consent was obtained from patients or legal representatives whenever possible. The ethics committee, in some cases, authorized an informed consent waiver.
Patients were classified according to their highest severity grade during the evolution of COVID-19 as follows: 1) Severe: i) death during hospitalization, ii) ICU admission, iii) invasive mechanical ventilation, or iv) presence of bilateral pulmonary infiltrates, mechanical ventilation and oxygen saturation (Sat02) ≤ 93%; 2) Moderate: the remaining patients admitted to hospital who did not fulfil severe COVID-19 criteria; 3) Asymptomatic/Mild (AM): individuals with minor or no COVID-19 symptoms; 4) Control: healthy donors recruited before COVID-19 outbreak. The STROBE-ID checklist was used to strengthen the design and conduct the study.
Clinical data and sample collection
Epidemiological and clinical variables were exhaustively collected from medical records using an electronic case report form (eCRF) built using REDCap electronic data capture tools [Citation13]. Plasma samples were collected at hospital entry or within the first days after hospitalization and before treatment with specific therapies for COVID-19 such as immunotherapy (tocilizumab, interferon-beta (IFNβ), corticoids, or ribavirin, among others).
Analysis of miRNome by high throughput sequencing
Total RNA, including small RNAs, was isolated from 400μl of plasma with the miRNeasy Serum Plasma Advanced kit (Qiagen, Hilden, Germany) following manufacturer's instructions. RNA quality and quantity were evaluated by the Bioanalyzer 2100 with Agilent RNA 6000 Nano Kit. Small RNA libraries were constructed using the NEBNext Multiplex Small RNA Library Prep for Illumina (New England Biolabs, Ipswich, MA, USA) at Parque Científico of Madrid (Spain). Next, small RNAs were sequenced in NovaSeq 6000, Illumina, at the National Centre for Microbiology (Majadahonda, Spain), estimating over 50 million reads per sample (Supplementary Data 1).
Raw data were analyzed using a specific bioinformatic pipeline detailed in Supplementary Data 2. Briefly, reads were quality checked with FastQC (v.0.11.3), and adapter sequences were trimmed with cutadapt (v.1.13). Reads were then processed with miRDeep2 (v.0.0.7) to identify and quantify known miRNAs from the reference human genome (GRCh38) and miRBase (v2.0) [Citation14].
Cytokines and chemokines quantification
Twenty-six plasma markers previously related to COVID-19 (Supplementary Data 3) were evaluated with a custom Procartaplex multiplex immunoassay (Invitrogen) using the Bio-plex 200TM system (BioRad Laboratories, Hercules, California, USA) following the manufacturer’s specifications. We used the raw fluorescence intensity (FI) values as a relative quantification of the analyte abundances, as previously described [Citation15].
Outcome variables
We analyzed several outcomes: i) infection: all SARS-CoV-2 infected patients were compared to healthy controls; ii) Symptomatology: symptomatic patients (moderate plus severe) versus asymptomatic; iii) Severity: severe versus moderate patients; iv) Mortality: 90-days mortality of SARS-CoV-2 infected patients, where the baseline was the date of diagnosis of SARS-CoV-2 infection.
Statistical analysis
The statistical analysis was performed by R statistical package (v4.0.3). Detailed statistical analysis information is available in Supplementary Data 4.
The Kruskal–Wallis test and Pearson’s chi-squared test were used for descriptive data for continuous and discrete variables. miRNA counts were normalized with the DESeq2 method (v1.30.0), and expression differences between severity groups were analyzed using a generalized linear model (GLM) with negative binomial distribution adjusted by age and gender. miRNAs with fold change (FC) ≥ 1.5 and q-value ≤ 0.05 (p-value corrected for the false discovery rate (FDR) by Benjamini-Hochberg correction) were considered significant.
Mortality analysis for 90 days after hospitalization was performed by a Cox Proportional-Hazard regression model adjusted by age and gender. Variance stabilizing transformation (VST) normalized counts of all miRNAs were included in the least absolute shrinkage and selection operator (LASSO) multivariate Cox regression model to select the miRNAs with the best predictive value and construct the miRNA mortality risk score (miRNA-MRS) within 90 days of hospital admission (see Supplementary Data 4 for extended information). The predictive performance for miRNA-MRS was assessed using the area under the receiver-operating characteristic curve (AUROC) and by calculating sensitivity, specificity and positive and negative predictive value.
Plasma cytokine and chemokine FI values were preprocessed with a weighted Box–Cox, followed by quantile normalization [Citation16]. Normalized data were analyzed with a GLM with gamma distribution. Correlation analysis between SDE miRNAs of each contrast and cytokines/chemokines was performed using the Spearman correlation test (q-value < 0.1).
miRNA-based target prediction and pathway enrichment analysis
Significant differentially expressed (SDE) miRNAs were analyzed for experimentally validated miRNA-target interactions with MIENTURNET (MicroRNA ENrichment TURned NETwork) [Citation17]. Gene targets with at least three validated interactions with SDE miRNAs were selected to perform functional enrichment analysis with the Reactome Pathway database [Citation18]. Enrichment p-values (Fischer’s exact test with hypergeometric distribution) were corrected for false discovery ratio (FDR) (q-value ≤ 0.1).
Results
Patient characteristics
Characteristics of the four groups of patients are shown in . There were no differences in the time point at which the plasma samples were obtained among groups, neither between symptoms onset and sampling. Ethnicity, smoker status, arterial hypertension, chronic pulmonary disease, COVID-19-related symptoms, and basal therapy with nonsteroidal anti-inflammatory drugs showed significant differences between groups. COVID-19 hospitalization requirements were found significantly different between groups, as they were used for severity definition (see methods): i.e. oxygenotherapy, mechanical ventilation, presence of bilateral pulmonary infiltrates, ICU and exitus.
Table 1. Clinical, epidemiological, and virological characteristics of SARS-CoV-2 infected patients.
miRNome characterization and correlation with cytokine release by severity groups
We obtained 50 million reads per sample, roughly three times over the minimum depth required for miRNA expression analysis. A total of 2656 miRNAs (91.97% of miRBase miRNAs) were identified, and 767 miRNAs were successfully filtered with sufficiently large counts for the differential expression analysis.
The partial least squares discriminant analysis (PLS-DA) showed differences between severity groups according to the miRNome expression (A-D), where each group showed a specific expression pattern (E).
Figure 1. Exploratory analysis. A to D. Multivariate analysis was performed by supervised partial least squares discriminant analysis (PLS-DA) from normalized log transformed and scaled miRNA expression data: A) all individuals (healthy controls, asymptomatic, moderate and severe); B) COVID+ patients vs. healthy controls; C) AM vs symptomatic; D) Severe vs moderate. Each symbol represents the miRNA profile of each participant of the study. E) Heatmap and hierarchical clustering of the SDE miRNAs. Top 20 SDE miRNAs for each comparison are shown: Healthy vs COVID+, AM vs Symptomatic and Moderate vs Severe. Study subjects are represented in columns and SDE miRNAs in rows, with clustering dendograms on the left for miRNAs and at the top for samples. The colour scale shows the relative expression level of SDE miRNAs. Red colour indicates a higher expression level and blue a lower expression level. Patients are grouped by severity. AM, asymptomatic; SDE, significant differentially expressed.
We also evaluated the change of cytokine and chemokine levels for all comparisons (Supplementary Data 5). Patients infected with SARS-CoV-2 showed higher levels of proinflammatory cytokines, which were increasing with COVID-19 severity. Later, we analyzed the correlation between SDE miRNAs and cytokines for each severity group contrast (). We observed that several cytokines, chemokines, and growth factors significantly correlated with the dysregulated miRNAs resulting from each comparison (A, COVID-19 patients vs. healthy controls; B, asymptomatic vs. symptomatic; C, severe vs. moderate).
Figure 2. Pearson correlation plot showing the correlation between SDE miRNAs and plasma cytokines/chemokines considering all patients for each comparison: A) COVID+ patients and healthy controls; B) Asymptomatic and symptomatic; C) Severe and moderate. The size of the circles is proportional to the strength of the correlation, and the colour represents the direction (colour legends are shown on the right), where large dark blue represents a strong negative correlation, and a large dark red circle represents a strong positive correlation. SDE miRNAs are on the horizontal axis and cytokines/chemokines on the vertical axis.
Characterization of COVID-19 infection
We first studied the miRNAs expression differences between SARS-CoV-2 infected patients and healthy controls, which showed the presence of 200 SDE miRNAs (A; Supplementary Data 6). In COVID-19 patients, 142 miRNAs were upregulated and 58 downregulated. SARS-CoV-2 infected patients had a strong upregulation of miRNAs such as the hsa-miR-4665-5p, hsa-miR-3190-3p, hsa-miR-331-3p, hsa-miR-4525, hsa-miR-431-5p, hsa-miR-6721-5p, hsa-miR-4661-5p, hsa-miR-548a-3p, hsa-miR-4745-5p, and hsa-miR-3150b-3p, which showed a FC > 100.
Figure 3. Differential expression analysis of miRNAs. Volcano plots showing the SDE miRNAs for each comparison: A. COVID+ vs healthy controls, B. AM vs symptomatic and C. Moderate vs severe. Red dots show miRNAs with a FDR corrected p-value ≤ 0.05 and a FC ≥ 1.5; Green dots show miRNAs with a FDR corrected p-value ≤ 0.05 and a FC ≤ −1.5; Blue dots represent miRNAs with a |FC| ≥ 1.5 that do not present statistical significance; Gray dots show miRNAs without statistical significance and a |FC| < 1.5. AM, asymptomatic; FDR, false discovery rate; FC, fold change.
The functional analysis reported 29 Reactome biological pathways. The most relevant pathways were influenza infection, influenza viral RNA transcription and replication, viral mRNA translation, and estrogen-dependent gene expression, among others (Supplementary Data 7; Supplementary Data 8).
SDE miRNAs in COVID-19 patients were significantly correlated with proinflammatory cytokines such as IL-6, IL-12, IP10, and TNFα (A), while healthy and infected patients showed different correlation patterns (Supplementary Data 9A and B, respectively). The proinflammatory hsa-miR-320 family was highly expressed in SARS-CoV-2 infected patients, positively correlated with proinflammatory cytokines such as IL-12, IL-6, IL-8, IP10, MIG, TIM3, and TNFα, especially the hsa-miR-320b. These correlations were not observed in healthy controls, where most of the significant correlations with proinflammatory cytokines were negative (Supplementary Data 9A).
Asymptomatic COVID-19 patients
We found 75 SDE miRNAs between asymptomatic (asymptomatic and mild) and symptomatic (moderate and severe) COVID-19 patients (B; Supplementary Data 10). Nineteen of them showed a higher expression in asymptomatic patients, and the hsa-miR-1291 was the most upregulated (FC = 14.6). In contrast, most of the SDE miRNAs were downregulated, showing nine of them over a 10-fold decrease (hsa-miR-8485, hsa-miR-1229-3p, hsa-miR-412-5p, hsa-miR-937-3p, hsa-miR-4750-5p, hsa-miR-296-5p, hsa-miR-26a-2-3p, hsa-miR-939-5p and hsa-miR-424-5p).
The functional pathway analysis showed that some of the most deregulated pathways were platelet aggregation (plug formation), IL-3, IL-5, and GM-CSF signalling, VEGFA-VEGFR2 pathway, and signalling by VEGF, among others (Supplementary Data 11; Supplementary Data 12). In line with these findings, asymptomatic patients showed lower levels of proangiogenic markers such as FGF2 and HGF than symptomatic patients, which were positively correlated with hsa-miR-4854 and hsa-miR-548t-5p, respectively (). This latter miRNA showed a high positive correlation with IFNγ, IL-1b, IL-12, IL-1RA, and IL-4 exclusively in asymptomatic participants (Supplementary Data 13A). Symptomatic patients showed a different correlation pattern, showing a high positive correlation of the hsa-miR-320 with ferritin, HGF, IL-12, IL-4, IL-6, IL-8, IP10, MIG, TIM3, and TNFα (Figure 13B).
Severe COVID-19 patients
We found 137 SDE miRNAs between severe and moderate COVID-19 patients (Supplementary Data 14). Seventy-seven miRNAs were upregulated in severe patients (C), being the hsa-miR-3976, hsa-miR-4488, hsa-miR-3150b-3p, hsa-miR-7704 and hsa-miR-3168 the most upregulated (FC > 10). Moreover, 60 miRNAs showed a downregulation in severe patients (C) and the most downregulated miRNAs with a 10-fold decrease (FC < 0.1) where the hsa-miR-7848-3p, hsa-miR-7110-3p, hsa-miR-1197, hsa-miR-4686, hsa-miR-374b-3p, hsa-miR-4473, hsa-miR-758-3p, hsa-miR-5004-3p, hsa-miR-146a-3p, hsa-miR-329-5p, hsa-miR-744-3p, hsa-miR-1285-5p and hsa-miR-6516-3p.
The pathway analysis indicated specific deregulation of several epigenetic pathways (PRC2 methylate histones and DNA, epigenetic regulation of gene expression, and SUMOylation of DNA methylation proteins) (Supplementary Data 15; Supplementary Data 16).
Severe patients showed an increase of proinflammatory cytokines, where the most elevated were IL-8, IL-6, IP10, MIG, MIP1, and TIM3, among others (Supplementary Data 5). The proinflammatory hsa-miR-320 family was highly expressed in severe patients (Supplementary Data 14), showing a high positive correlation with proinflammatory cytokines.
Conversely, the hsa-miR-26a-5p (a negative regulator of the inflammatory response) and the miR-98-5p (a member of the let-7 miRNA family) showed lower levels in severe patients and were negatively correlated with all proinflammatory cytokines (C), (Supplementary Data 17B).
Mortality predictive model
Characteristics of patients by mortality status are shown in Supplementary Data 18. We identified 12 miRNAs as predictors of 90-days mortality (hsa-let-7f-1-3p, hsa-let-7g-5p, hsa-miR-1255a, hsa-miR-140-3p, hsa-miR-20a-5p, hsa-miR-22-3p, hsa-miR-3180-3p, hsa-miR-3180, hsa-miR-363-5p, hsa-miR-4510, hsa-miR-548h-3p and hsa-miR-6130) (Supplementary Data 19, 20A and 20B), but two of them (hsa-miR-548h-3p and hsa-miR-3180) did not meet the filtering criteria (Supplementary Data 19). Among the ten selected miRNAs, only two were associated with a lower probability of survival at 90 days (hsa-miR-22-3p and hsa-miR-3180-3p), whereas eight were associated with a better prognosis (hsa-let-7f-1-3p, hsa-let-7g-5p, hsa-miR-1255a, hsa-miR-140-3p, hsa-miR-20a-5p, hsa-miR-363-5p, hsa-miR-4510, and hsa-miR-6130).
This ten miRNA-expression signature (Supplementary Data 21) was used to generate the miRNA-MRS, whose median value was −3.992, which was used as a threshold to stratify the study subjects into low-risk (risk score < −3.992, n = 42) and high-risk (risk score ≥ −3.992, n = 42) (A). Kaplan-Meier analysis confirmed the significant differences in the overall mortality (OS) (p-value <0.001) (B). Moreover, ROC curve analysis confirmed the predictive power of the miRNA-MRS, which was improved by adjusting with age and gender, showing a significantly higher AUROC value than the basic model considering only age and gender (0.968 vs. 0.881, p-value = 0.042), (C). Additionally, we analyzed the diagnostic accuracy of adding the miRNA-MRS to the model with age and gender by calculating sensitivity, specificity, and positive and negative predictive value. Considering the cut-off for maximum sensitivity and specificity, the model showed sensitivity and specificity values above 90% (Supplemental Data 22). Additionally, results were confirmed by the multivariate Cox-regression model, highlighting the independent prognostic value of the risk-score (hazard ratio (HR) = 4.599; 95% confidence interval (CI) (1.977− 10.700), p-value <0.001) (D).
Figure 4. Survival analysis of COVID-19 patients and 10-miRNA signature risk score. A. Distribution of risk score values for low-risk and high-risk subjects. B. KM survival curve between the high-risk and low-risk groups stratified by median risk score. C. ROC curve of survival showing the prognostic ability (AUC) of age + gender (red) and miRNA risk score + age + gender (blue). D. Forest plot of the multivariate Cox regression analysis for age, gender and 10-miRNA risk score. The HR of each variable is indicated with 95% CI and statistical significance (p-value). KM, Kaplan-Meier; ROC, receiver operating characteristic; AUC, area under the curve; HR, hazard ratio; CI, confidence interval.
Discussion
Our results showed that SARS-CoV-2 infection produces a strong disturbance of miRNA profile from the onset of the disease, exhibiting at this time a specific miRNA signature for each COVID-19 severity. The identified miRNAs could be used as early severity markers at hospital admission. To our knowledge, this is the first large-scale study that characterizes the expression profile of plasma miRNAs in COVID-19 patients of different clinical statuses, identifying specific miRNA signature for each severity group and a specific miRNA mortality risk score (miRNA-MRS).
miRNA altered profile by SARS-CoV-2 infection
We identified 200 SDE miRNAs between COVID-19 patients and healthy controls, most of them upregulated in SARS-CoV-2 infected patients. Although scarce studies on human miRNAs during COVID-19 have been published, some of the most upregulated miRNAs in COVID-19 patients have been identified to target the SARS-CoV-2 genome, which supports their functionality. This is the case of hsa-miR-4665-5p (predicted binding site in SARS-CoV-2 genome) [Citation19], hsa-miR-3190-3p (predicted targeting of ORF1ab and ORF8 genes) [Citation20], hsa-miR-4745-5p (expressed in lung tissue with affinity to the 3′UTR of SARS-COV-2 genome) [Citation21] and hsa-miR-3150b-3p (predicted binding site in the leader sequence of SARS-CoV-2) [Citation22]. Remarkably, other strongly dysregulated miRNAs, such as hsa-miR-431-5p, hsa-miR-548a-3p, and hsa-miR-331-3p have been previously linked to several viral infections, including ZIKV and HBV [Citation23,Citation24]. It has also been suggested that hsa-miR-431-5p plays a role in restricting HCV replication after its induction by IFNβ [Citation25]. Due to the known antiviral activity of IFNβ, constituting the first line of defence against viral infections, IFNβ-regulated miRNAs could be used as early biomarkers of COVID-19 disease, being detected before the traditional markers of immune activation or inflammation.
Disrupted miRNAs in COVID-19 patients were strongly involved in biological processes related to viral infections, including influenza viral RNA transcription and replication or viral mRNA translation, suggesting a direct binding/effect of the SDE miRNAs on the viral genome. Remarkably, biological pathways involving eukaryotic mRNA translation also appeared as dysregulated, including cap-dependent translation initiation, translation initiation complex formation, elongation, and termination, among others. These functional alterations agree with in vitro studies that point to the inhibition of cellular mRNA expression in cells infected with different coronaviruses, including SARS-CoV-2 [Citation26]. Coronaviruses seem to hijack the host translational machinery, thus favouring viral protein production. In fact, Kim et al. showed that host gene expression is almost suppressed in the first hours of infection [Citation27], which remarks the early disturbance of miRNA machinery, as miRNAs are the master regulators of gene expression.
Additionally, certain trace elements pathways such as metabolism and transport of copper (Cu) and selenium (Se) appeared also disrupted. These elements are essential for a normal immune response, and their serum levels have been proposed as good predictors of COVID-19 survival [Citation28].
However, the most dysregulated was the estrogen signalling pathway. Estrogen plays a crucial role in innate and adaptive immune responses, showing an immunoenhancing effect. Estrogen receptors are involved in tissue repair processes during respiratory virus infection [Citation29]. A weaker estrogen receptor signalling is linked to a poor evolution of influenza, while estrogen treatment help to control the inflammatory reactions and improve survival [Citation30]. An anti-inflammatory effect of estrogens in COVID-19 has been suggested, reducing SARS-CoV-2 infectivity through the regulation of proinflammatory signalling pathways [Citation31]. This evidence, along with our results and ongoing clinical trials [Citation32], point to the necessity of continuing to investigate the effects of these molecules in COVID-19, to discover the potential use of estrogens as therapy or preventive strategy against COVID-19.
miRNA profile of asymptomatic infection
We identified a signature of 75 SDE miRNAs in asymptomatic with respect to symptomatic patients that could be potentially used as early predictors of a favourable clinical course. Some of the most upregulated miRNAs in symptomatic patients have previously been related to viral infections, including influenza A viruses (IAVs), HIV1, and HIV2 viruses [Citation33,Citation34]. The hsa-miR-424-5p targets different RNA viruses, including SARS-CoV-2, and the hsa-miR-939-5p targets the SARS-CoV-2 genome 3′UTR [Citation21]. Additionally, the hsa-miR-939-5p regulates several anti-inflammatory genes, including IL-6, VEGFA, TNFα, NFκB2, and NOS2A, showing a key role as a mediator in inflammatory networks [Citation35].
We also observed up to 19 upregulated miRNAs in asymptomatic patients. Among the most upregulated miRNAs, hsa-miR-708-3p and hsa-miR-1291 were recently found by Chen et al. as differentially expressed between patients with mild and severe symptoms of COVID-19 [Citation36]. The hsa-miR-1291 was more than 14 times expressed in our asymptomatic patients. This miRNA acts upstream of the chemokine receptor CCR2, which is known to induce monocyte and macrophage recruitment to sites of inflammation and regulates the immune response by controlling the proportion of effector and regulatory T cells [Citation36]. An increased CCR2 expression has been reported in patients with severe COVID-19 [Citation37], and similarly, hsa-miR-1291 has been pinpointed as a potential early biomarker of COVID-19 severity [Citation36]. Therefore, tight control of inflammation by hsa-miR-1291 could be behind the COVID-19 control in asymptomatic patients.
Remarkably, asymptomatic patients showed a higher expression of hsa-miR-150-5p, a master regulator of the immune system with crucial roles in the development of immune response, including respiratory viral infections [Citation38]. Hsa-miR-150-5p is involved in the modulation of B cell differentiation, development of natural killer (NK), and invariant NK T (iNKT) cells by targeting the transcription factor c-Myb [Citation39], among others. Critically ill patients showed decreased levels of hsa-miR-150-5p, and its expression is inversely correlated with the number of days of ICU stay [Citation37]. Thus, a higher expression of hsa-miR-150-5p could lead to a fine-tune immune response to successfully clear the SARS-CoV-2 virus without an exaggerated immune response.
The functional analysis further confirmed these results, which showed inflammatory and immune processes as the most dysregulated biological events by the SDE miRNAs. In fact, we observed disruption of platelet aggregation, interleukin-3, −5 and GM-CSF signalling, VEGFA-VEGFR2, and signalling by VEGF. Although we could not evaluate VEGF, higher plasma levels of other angiogenic growth factors such as HGF and FGF were observed in asymptomatic patients. Moreover, a strong enrichment of phagocytic processes, including FcγRIIIa-mediated phagocytosis, regulation of actin dynamics for phagocytic cup formation, and Fcgamma receptor (FCγR) dependent phagocytosis was observed. Phagocytosis is one of the main processes in innate immune response, contributing to the fight against infectious agents, as it has been detected in COVID-19 severe patients [Citation40]. The FCγR is also involved in cell activation and increased proinflammatory cytokine production [Citation41], which could be promoting the outbreak of COVID-19 cytokine storm.
Early miRNA signature predicting the severity of COVID-19
Certain risk factors linked to a more severe course of COVID-19 have already been identified, such as age or male gender. However, these factors cannot predict the evolution of the disease by themselves, and it remains essential to find additional biological factors decisive to trigger COVID-19 severity. In this sense, we identified the molecular players in the early stages that could be behind the progression to a severe COVID-19.
We identified 137 SDE miRNAs between severe patients and those with moderate symptoms. Seventy-seven of these miRNAs were upregulated in severe patients, with hsa-miR-3976, hsa-miR-4488, hsa-miR-3150b-3p, hsa-miR-7704, and hsa-miR-3168 showing the strongest disruption. Most of them showed some evidence of playing a role in different viral infections [Citation42,Citation43]. Moreover, hsa-miR-3150b-3p and hsa-miR-7704 are putative targets to the leader sequence and ORF3a gene of the SARS-CoV-2 genome, respectively [Citation20,Citation22]. Meanwhile, 60 miRNAs showed downregulation in severe patients. Again, some of the strongly disrupted miRNAs potentially target the SARS-CoV-2 genome, such as hsa-miR-1197, hsa-mir-5004-3p, and hsa-miR-1285-5p [Citation19,Citation44], which could be involved in inhibiting virus replication. In this regard, our functional enrichment results suggest a strong involvement of SDE miRNAs in biological processes such as viral RNA transcription, replication and translation, further supporting this hypothesis. Several examples of host miRNA implication on RNA viruses replication and pathogenesis have been described, as these viruses may modify miRNAs profile within the cell to establish a proviral environment [Citation45]. Our data suggest that the altered miRNA profile observed in COVID-19 patients could lead to a differential viral load according to the severity of the disease. Similarly, RNAemia and viral load on admission have been related to mortality [Citation46]. However, future studies would be necessary to clarify the interplay between miRNAs, RNAemia and viral load.
Additionally, we observed a strong disturbance of miRNAs that modulate the immune response, suggesting that the differential severity status could be a result of both the altered regulation of viral translation/replication processes but also the differential modulation of immune / inflammatory responses following infection.
In fact, other strongly downregulated miRNAs in severe patients directly regulate antiviral response, such as the hsa-miR-744-3p, against IAV and RSV through targeting MAPK-activated protein kinase 2 [Citation47]. On the contrary, hsa-miR-758-3p downregulates the expression of the viral RNA receptors Toll-like receptors 3 (TLR3) and 7 (TLR7) following HCV infection [Citation48]. Additionally, we observed a deep downregulation of hsa-miR-146a-3p, involved in several autoimmune and inflammatory diseases such as cystic fibrosis, where it controls the IL-6 production, contributing to the restriction of the inflammatory response [Citation49]. Thus, the reduced expression of this miRNA in severe patients is related to an increase in IL-6 production, which has been reported in severe and critical COVID-19 patients, where IL-6 acts as a vital amplifier [Citation50]. Besides, we also observed an increase in additional proinflammatory cytokines such as IL-8, IP10, MIG, MIP-1, and TIM3, among others.
The dysregulated miRNAs in severe patients were mainly involved in epigenetic pathways, which have been reported as altered in COVID-19 [Citation51].
Ten miRNA-based risk score as a predictor of COVID-19 mortality
The overall mortality rate of COVID-19 is around 1-3% in most countries, and it is hugely affected by certain epidemiological factors, such as age, gender, and comorbidities [Citation52]. However, these estimates are highly dependent on the countries healthcare systems and vaccination progress. While no differences in comorbidities or treatments by mortality status were observed in this study, we identified early biomarkers that allowed predicting an increased risk of death. We screened 767 miRNAs using a LASSO Cox regression model and identified a signature of ten SDE miRNAs that can be used as good predictors of COVID-19 mortality.
While little is known about the potential role of hsa-miR-1255a and hsa-miR-6130 in COVID-19 disease, some of the selected miRNAs have already been predicted to target the SARS-CoV-2 genome, including hsa-let-7g-5p, hsa-miR-363-5p, and hsa-miR-4510 [Citation20]. Moreover, a decreased hsa-miR-20a-5p expression in peripheral blood was observed in patients with COVID-19 [Citation53]. This miRNA has been computationally predicted to target several proteins that may play a role in SARS-CoV-2 induced cell death, including SMAD3, SMAD4, and TGFBR1 [Citation54]. In fact, the binding of SARS-CoV N protein to SMAD3 has been shown to interfere with the complex formation between SMAD3 and SMAD4, thus inhibiting the apoptosis of SARS-CoV-infected cells mediated by the TGF-signalling pathway [Citation55], in favour of virus packaging and replication.
A link between COVID-19 disease and the miRNAs hsa-let-7f-1-3p and hsa-miR-140-3p has also been reported, showing a downregulation in cells infected with SARS-CoV-2 [Citation19,Citation56]. Kim et al. proposed that hsa-miR-140-3p targets the serine protease TMPRSS2, which processes the viral S protein that binds to the ACE2 host receptor [Citation56]. Moreover, several studies have described its potential to modulate Bcl-2 expression, which has an essential role in regulating several cellular processes that could have an impact on viral infection progression [Citation57].
Although both age and gender have been widely reported as risk factors for the greater severity and mortality of COVID-19, we found that the sensitivity and specificity in the prediction of COVID-19 mortality significantly increased when including the ten-miRNA signature-based risk score in our model. Moreover, we confirmed the independent predictive value of the risk-score by a multivariate Cox-regression model including both epidemiological variables.
Although further research is still needed to decipher the biological mechanisms that could be behind the potential effect of these human miRNAs in COVID-19 mortality, our results provide a new promising diagnostic tool that could contribute to the improved management of patients with a higher risk of death at disease onset.
This miRNA-MRS could be implemented in routine diagnostics through an amplification-based detection method such as multiplex real time quantitative PCR, or even through emerging technologies such as hybridization-based detection methods [Citation58]. In this latter case, low amount of starting material and no reverse transcription or cDNA amplification is needed, which make this medium-throughput technology highly suitable for a scalable use [Citation59]. These systems can be easily combined with an analysis software to give a simple readout indicating the disease probability score based on the concentration of each miRNA [Citation60].
Nevertheless, this study has several limitations. First, our cohort included a limited number of deaths among participants. Thus, we decided to include all the study samples as a training set in order to obtain more reliable and representative results. Therefore, our findings would need to be validated in an external cohort of COVID-19 patients. Nevertheless, at the time of manuscript submission, there were no high-throughput miRNA sequencing data available from studies of COVID-19 patients reporting mortality analysis results to perform the validation of our miRNA-based model. Moreover, most of the participants in our cohort were Caucasian and Hispanic, while other ethnic groups were underrepresented. Due to the differences in the genetic background and miRNA expression, the predictive value of the miRNA-based risk score should also be confirmed in studies with ethnically more diverse populations. Additionally, it would also be interesting to study whether miRNA profiles correlate with SARS-CoV-2 viral load or RNAemia. However, we were unable to analyze it because no remaining sample was available. Finally, further studies would be needed to explore how the duration of illness influences the microRNA profiles.
Conclusion
SARS-CoV-2 infection deeply disturbs the plasma miRNA expression profile from an early stage of COVID-19, making miRNAs highly valuable as early predictors of severity and mortality.
Author contribution
Funding body: AFR and MAJS.
Study concept and design: AFR and MAJS
Patients’ selection and clinical data acquisition: PR, RB, MMV, AVB, OBK, FCC, FPG, OMG, AM, NBL, JT, DA, LMP, IRMA, NBL.
Sample preparation, and plasma biomarker analysis: MMV, AVB, EJVA, SBS, OBK.
Statistical analysis and interpretation of data: AFP, AVB, OBK, DVM, AFR, and MAJS.
Writing of the manuscript: AFP, AFR, and MAJS.
Critical revision of the manuscript for relevant intellectual content: SR, PR, RB, FPG, OMG.
Supervision and visualization: AFR and MAJS.
All authors read and approved the final manuscript.
Supplemental Material
Download MS Word (3.2 MB)Acknowledgements
This study would not have been possible without the collaboration of all the patients, their families, medical and nursery staff, and data managers who have taken part in the project. We would like also to acknowledge the Genomics Unit, and the Bioinformatics Unit at Instituto de Salud Carlos III, whose assistance was essential.
Disclosure statement
No potential conflict of interest was reported by the author(s).
Additional information
Funding
References
- Oude Munnink BB, Nieuwenhuijse DF, Stein M, et al. Rapid SARS-CoV-2 whole-genome sequencing and analysis for informed public health decision-making in the Netherlands. Nat Med. 2020 09;26(9):1405–1410.
- Chandan K, Gupta M, Sarwat M. Role of Host and Pathogen-Derived MicroRNAs in Immune Regulation During Infectious and Inflammatory Diseases. Front Immunol. 2019;10:3081.
- Hum C, Loiselle J, Ahmed N, et al. MicroRNA Mimics or Inhibitors as Antiviral Therapeutic Approaches Against COVID-19. Drugs. 2021 Apr;81(5):517–531.
- Feng S, Zeng D, Zheng J, et al. MicroRNAs: Mediators and Therapeutic Targets to Airway Hyper Reactivity After Respiratory Syncytial Virus Infection. Front Microbiol. 2018;9:2177.
- Hasan MM, Akter R, Ullah MS, et al. A Computational Approach for Predicting Role of Human MicroRNAs in MERS-CoV Genome. Adv Bioinformatics. 2014;2014:967946.
- de Gonzalo-Calvo D, Benítez ID, Pinilla L, et al. Circulating microRNA profiles predict the severity of COVID-19 in hospitalized patients. Transl Res. 2021 Oct;236:147–159..
- Ying H, Ebrahimi M, Keivan M, et al. miRNAs; a novel strategy for the treatment of COVID-19. Cell Biol Int. 2021 Oct;45(10):2045–2053.
- Chan AP, Choi Y, Schork NJ. Conserved Genomic Terminals of SARS-CoV-2 as Coevolving Functional Elements and Potential Therapeutic Targets. mSphere. 2020 11 25;5(6):e00754–20.
- Narożna M, Rubiś B. Anti-SARS-CoV-2 Strategies and the Potential Role of miRNA in the Assessment of COVID-19 Morbidity, Recurrence, and Therapy. Int J Mol Sci. 2021 Aug 12;22(16):8663.
- Chen L, Xie W, Wang L, et al. miRNA-133a aggravates inflammatory responses in sepsis by targeting SIRT1. Int Immunopharmacol. 2020 Nov;88:106848.
- Goodwin AJ, Li P, Halushka PV, et al. Circulating miRNA 887 is differentially expressed in ARDS and modulates endothelial function. Am J Physiol Lung Cell Mol Physiol. 2020 Jun 1;318(6):L1261–L1269.
- Song X, Wang CT, Geng XH. MicroRNA-29a promotes apoptosis of monocytes by targeting STAT3 during sepsis. Genet Mol Res. 2015 Oct 29;14(4):13746–13753.
- Harris PA, Taylor R, Minor BL, et al. The REDCap consortium: Building an international community of software platform partners. J Biomed Inform. 2019 07;95:103208.
- Kozomara A, Birgaoanu M, Griffiths-Jones S. miRBase: from microRNA sequences to function. Nucleic Acids Res. 2019 01 08;47(D1):D155–D162.
- Breen EJ, Tan W, Khan A. The Statistical Value of Raw Fluorescence Signal in Luminex xMAP Based Multiplex Immunoassays. Sci Rep. 2016 05;6:26996.
- Rausch TK, Schillert A, Ziegler A, et al. Comparison of pre-processing methods for multiplex bead-based immunoassays. BMC Genomics. 2016 Aug 11;17(1):601.
- Licursi V, Conte F, Fiscon G, et al. MIENTURNET: an interactive web tool for microRNA-target enrichment and network-based analysis. BMC Bioinformatics. 2019 Nov 04;20(1):545.
- Jassal B, Matthews L, Viteri G, et al. The reactome pathway knowledgebase. Nucleic Acids Res. 2020 01 08;48(D1):D498–D503.
- Chow JT, Salmena L. Prediction and Analysis of SARS-CoV-2-Targeting MicroRNA in Human Lung Epithelium. Genes (Basel). 2020 08 26;11(9):1002.
- Saçar Demirci MD, Adan A. Computational analysis of microRNA-mediated interactions in SARS-CoV-2 infection. PeerJ. 2020;8:e9369.
- Balmeh N, Mahmoudi S, Mohammadi N, et al. Predicted therapeutic targets for COVID-19 disease by inhibiting SARS-CoV-2 and its related receptors. Inform Med Unlocked. 2020;20:100407.
- Natarelli L, Parca L, Mazza T, et al. MicroRNAs and Long Non-Coding RNAs as Potential Candidates to Target Specific Motifs of SARS-CoV-2. Noncoding RNA. 2021 Feb 18;7(1):14.
- Seong RK, Lee JK, Cho GJ, et al. mRNA and miRNA profiling of Zika virus-infected human umbilical cord mesenchymal stem cells identifies miR-142-5p as an antiviral factor. Emerg Microbes Infect. 2020 Dec;9(1):2061–2075.
- Cao Y, Chen J, Wang D, et al. Upregulated in Hepatitis B virus-associated hepatocellular carcinoma cells, miR-331-3p promotes proliferation of hepatocellular carcinoma cells by targeting ING5. Oncotarget. 2015 Nov 10;6(35):38093–38106.
- Pedersen IM, Cheng G, Wieland S, et al. Interferon modulation of cellular microRNAs as an antiviral mechanism. Nature. 2007 Oct 18;449(7164):919–922.
- Banerjee AK, Blanco MR, Bruce EA, et al. SARS-CoV-2 Disrupts Splicing, Translation, and Protein Trafficking to Suppress Host Defenses. Cell. 2020 11 25;183(5):1325–1339.e21.
- Kim D, Lee JY, Yang JS, et al. The Architecture of SARS-CoV-2 Transcriptome. Cell. 2020 05 14;181(4):914–921.e10.
- Hackler J, Heller RA, Sun Q, et al. Relation of Serum Copper Status to Survival in COVID-19. Nutrients. 2021 May 31;13(6):1898.
- Kadel S, Kovats S. Sex Hormones Regulate Innate Immune Cells and Promote Sex Differences in Respiratory Virus Infection. Front Immunol. 2018;9:1653.
- Robinson DP, Lorenzo ME, Jian W, et al. Elevated 17β-estradiol protects females from influenza A virus pathogenesis by suppressing inflammatory responses. PLoS Pathog. 2011 Jul;7(7):e1002149.
- Al-Kuraishy HM, Al-Gareeb AI, Faidah H, et al. The Looming Effects of Estrogen in Covid-19: A Rocky Rollout. Front Nutr. 2021;8:649128.
- Adel Ganaw MPF, Madhani M, Montgomery H, et al. Oestrogen Treatment for COVID-19 Symptoms. In: Corporation HM, editor. Clinical Trial. Qatar: Hamad Medical Corporation; 2021.
- Gao J, Gao L, Li R, et al. Integrated analysis of microRNA-mRNA expression in A549 cells infected with influenza A viruses (IAVs) from different host species. Virus Res. 2019 04 02;263:34–46.
- Biswas S, Chen E, Haleyurgirisetty M, et al. Comparison of miRNA Expression Profiles between HIV-1 and HIV-2 Infected Monocyte-Derived Macrophages (MDMs) and Peripheral Blood Mononuclear Cells (PBMCs). Int J Mol Sci. 2020 Sep 22;21(18):6970.
- McDonald MK, Ramanathan S, Touati A, et al. Regulation of proinflammatory genes by the circulating microRNA hsa-miR-939. Sci Rep. 2016 08 08;6:30976.
- Chen Z, Wang X, Li L, et al. Construction of an autophagy interaction network based on competitive endogenous RNA reveals the key pathways and central genes of SARS-CoV-2 infection in vivo. Microb Pathog. 2021 Jun 19;158:105051.
- Kim CW, Oh JE, Lee HK. Single Cell Transcriptomic Re-analysis of Immune Cells in Bronchoalveolar Lavage Fluids Reveals the Correlation of B Cell Characteristics and Disease Severity of Patients with SARS-CoV-2 Infection. Immune Netw. 2021 Feb;21(1):e10.
- Morán J, Ramírez-Martínez G, Jiménez-Alvarez L, et al. Circulating levels of miR-150 are associated with poorer outcomes of A/H1N1 infection. Exp Mol Pathol. 2015 Oct;99(2):253–261.
- Bezman NA, Chakraborty T, Bender T, et al. miR-150 regulates the development of NK and iNKT cells. J Exp Med. 2011 Dec 19;208(13):2717–2731.
- Chakraborty S, Gonzalez J, Edwards K, et al. Proinflammatory IgG Fc structures in patients with severe COVID-19. Nat Immunol. 2021 Jan;22(1):67–73.
- Vogelpoel LT, Baeten DL, de Jong EC, et al. Control of cytokine production by human fc gamma receptors: implications for pathogen defense and autoimmunity. Front Immunol. 2015;6:79.
- Shabani M, Nasr Esfahani B, Sadegh Ehdaei B, et al. Inhibition of herpes simplex virus type 1 replication by novel hsa-miR-7704. Res Pharm Sci. 2019 Apr;14(2):167–174.
- Chahar HS, Corsello T, Kudlicki AS, et al. Respiratory Syncytial Virus Infection Changes Cargo Composition of Exosome Released from Airway Epithelial Cells. Sci Rep. 2018 01 10;8(1):387.
- Mohammadi-Dehcheshmeh M, Moghbeli SM, Rahimirad S, et al. A Transcription Regulatory Sequence in the 5’ Untranslated Region of SARS-CoV-2 Is Vital for Virus Replication with an Altered Evolutionary Pattern against Human Inhibitory MicroRNAs. Cells. 2021 02 04;10(2):319.
- Trobaugh DW, Klimstra WB. MicroRNA Regulation of RNA Virus Replication and Pathogenesis. Trends Mol Med. 2017 01;23(1):80–93.
- Bermejo-Martin JF, González-Rivera M, Almansa R, et al. Viral RNA load in plasma is associated with critical illness and a dysregulated host response in COVID-19. Crit Care. 2020 12 14;24(1):691.
- McCaskill JL, Ressel S, Alber A, et al. Broad-Spectrum Inhibition of Respiratory Virus Infection by MicroRNA Mimics Targeting p38 MAPK Signaling. Mol Ther Nucleic Acids. 2017 Jun 16;7:256–266.
- Yang Q, Fu S, Wang J. Hepatitis C virus infection decreases the expression of Toll-like receptors 3 and 7 via upregulation of miR-758. Arch Virol. 2014 Nov;159(11):2997–3003.
- Luly FR, Lévêque M, Licursi V, et al. MiR-146a is over-expressed and controls IL-6 production in cystic fibrosis macrophages. Sci Rep. 2019 11 07;9(1):16259.
- Sabbatinelli J, Giuliani A, Matacchione G, et al. Decreased serum levels of the inflammaging marker miR-146a are associated with clinical non-response to tocilizumab in COVID-19 patients. Mech Ageing Dev. 2021 01;193:111413.
- Herbein G. An epigenetic signature to fight COVID-19. EBioMedicine. 2021 May;67:103385.
- Wiersinga WJ, Rhodes A, Cheng AC, et al. Pathophysiology, Transmission, Diagnosis, and Treatment of Coronavirus Disease 2019 (COVID-19): A Review. JAMA. 2020 Aug 25;324(8):782–793.
- Li C, Hu X, Li L, et al. Differential microRNA expression in the peripheral blood from human patients with COVID-19. J Clin Lab Anal. 2020 Oct;34(10):e23590.
- Yousefi H, Poursheikhani A, Bahmanpour Z, et al. SARS-CoV infection crosstalk with human host cell noncoding-RNA machinery: An in-silico approach. Biomed Pharmacother. 2020 Oct;130:110548.
- Zhao X, Nicholls JM, Chen YG. Severe acute respiratory syndrome-associated coronavirus nucleocapsid protein interacts with Smad3 and modulates transforming growth factor-beta signaling. J Biol Chem. 2008 Feb 08;283(6):3272–3280.
- Kim WR, Park EG, Kang KW, et al. Expression Analyses of MicroRNAs in Hamster Lung Tissues Infected by SARS-CoV-2. Mol Cells. 2020 Nov 30;43(11):953–963.
- Zhu ZR, He Q, Wu WB, et al. MiR-140-3p is Involved in In-Stent Restenosis by Targeting C-Myb and BCL-2 in Peripheral Artery Disease. J Atheroscler Thromb. 2018 Nov 01;25(11):1168–1181.
- Gustafson D, Tyryshkin K, Renwick N. microRNA-guided diagnostics in clinical samples. Best Pract Res Clin Endocrinol Metab. 2016 10;30(5):563–575.
- Kappel A, Keller A. miRNA assays in the clinical laboratory: workflow, detection technologies and automation aspects. Clin Chem Lab Med. 2017 May 01;55(5):636–647.
- Gliddon HD, Herberg JA, Levin M, et al. Genome-wide host RNA signatures of infectious diseases: discovery and clinical translation. Immunology. 2018 02;153(2):171–178.