Background MicroRNAs (miRNAs) are a class of small non-coding RNA that can downregulate their targets by selectively binding to the 3′ untranslated region (3′UTR) of most messenger RNAs (mRNAs) in the human genome. MiRNAs can interact with other molecules such as viruses and act as a mediator for viral infection. In this study, we examined whether, and to what extent, the SARS-CoV-2 virus can serve as a “sponge” for human miRNAs. Results We identified multiple potential miRNA/target pairs that may be disrupted during SARS-CoV-2 infection. Using miRNA expression profiles and RNA-seq from published studies, we further identified a highly confident list of 5 miRNA/target pairs that could be disrupted by the virus’s miRNA sponge effect, namely hsa-miR-374a-5p/APOL6, hsa-let-7f-1-3p/EIF4A2, hsa-miR-374a-3p/PARP11, hsa-miR-548d-3p/PSMA2 and hsa-miR-23b-3p/ZNFX1 pairs. Using single-cell RNA-sequencing based data, we identified two important miRNAs, hsa-miR-302c-5p and hsa-miR-16-5p, to be potential virus targeting miRNAs across multiple cell types from bronchoalveolar lavage fluid samples. We further validated some of our findings using miRNA and gene enrichment analyses and the results confirmed with findings from previous studies that some of these identified miRNA/target pairs are involved in ACE2 receptor network, regulating pro-inflammatory cytokines and in immune cell maturation and differentiation. Conclusion Using publicly available databases and patient-related expression data, we found that acting as a “miRNA sponge” could be one explanation for SARS-CoV-2-mediated pathophysiological changes. This study provides a novel way of utilizing SARS-CoV-2 related data, with bioinformatics approaches, to help us better understand the etiology of the disease and its differential manifestation across individuals.

The SARS-CoV-2 (COVID-19) pandemic has caused millions of deaths worldwide. Early risk assessment of COVID-19 cases can help direct early treatment measures that have been shown to improve the prognosis of severe cases. Currently, circulating miRNAs have not been evaluated as canonical COVID-19 biomarkers, and identifying biomarkers that have a causal relationship with COVID-19 is imperative. To bridge these gaps, we aim to examine the causal effects of miRNAs on COVID-19 severity in this study using two-sample Mendelian randomization approaches. Multiple studies with available GWAS summary statistics data were retrieved. Using circulating miRNA expression data as exposure, and severe COVID-19 cases as outcomes, we identified ten unique miRNAs that showed causality across three phenotype groups of COVID-19. Using expression data from an independent study, we validated and identified two high-confidence miRNAs, namely, hsa-miR-30a-3p and hsa-miR-139-5p, which have putative causal effects on developing cases of severe COVID-19. Using existing literature and publicly available databases, the potential causative roles of these miRNAs were investigated. This study provides a novel way of utilizing miRNA eQTL data to help us identify potential miRNA biomarkers to make better and early diagnoses and risk assessments of severe COVID-19 cases.