- Research
- Open access
- Published:
The relative expression of miR-31, miR-29, miR-126, and miR-17 and their mRNA targets in the serum of COVID-19 patients with different grades during hospitalization
European Journal of Medical Research volume 26, Article number: 75 (2021)
Abstract
Background
The aim of this study was to evaluate the expression of four up/down-regulated inflammatory miRNAs and their mRNA targets in the serum samples of COVID-19 patients with different grades. Also, we investigated the relative expression of these miRNAs and mRNAs during hospitalization.
Methods
In this cross-sectional study, 5 mL of blood sample were taken from COVID-19 patients with different grades and during hospitalization from several health centers of Yazd, Tehran, and Zahedan province of Iran from December 20, 2020 to March 2, 2021. The relative expression of miRNAs and mRNAs was evaluated by q-PCR.
Results
We found that the relative expression of hsa-miR-31-3p, hsa-miR-29a-3p, and hsa-miR-126-3p was significantly decreased and the relative expression of their mRNA targets (ZMYM5, COL5A3, and CAMSAP1) was significantly increased with the increase of disease grade. Conversely, the relative expression of hsa-miR-17-3p was significantly increased and its mRNA target (DICER1) was significantly decreased with the increase of disease grade. This pattern was exactly seen during hospitalization of COVID-19 patients who did not respond to treatment. In COVID-19 patients who responded to treatment, the expression of selected miRNAs and their mRNA targets returned to the normal level. A negative significant correlation was seen between (1) the expression of hsa-miR-31-3p and ZMYM5, (2) hsa-miR-29a-3p and COL5A3, (3) hsa-miR-126-3p and CAMSAP1, and (4) hsa-miR-17-3p and DICER1 in COVID-19 patients with any grade (P < 0.05) and during hospitalization.
Conclusions
In this study, we gained a more accurate understanding of the expression of up/down-regulated inflammatory miRNAs in the blood of COVID-19 patients. The obtained data may help us in the diagnosis and prognosis of COVID-19.
Trial registration: The ethics committee of Zahedan University of Medical Sciences, Zahedan, Iran. (Ethical Code: IR.ZAUMS.REC.1399.316) was registered for this project.
Introduction
Coronavirus 2019 (COVID-19) is a worldwide pandemic disease with near 2% mortality [1]. The disease is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and mainly affects the lungs and causes acute respiratory distress syndrome (ARDS) [2]. In COVID-19, there are significant changes in the function and morphology of lung endothelial cells, disruption of intercellular connections, and loss of the basal membrane [3]. From a molecular point of view, the S protein of SARS-CoV-2 binds to the angiotensin-converting enzyme 2 (ACE2), a vital protein in human cells, and then invigilated into target cells [4]. Of course, it must be said that not only the lungs are affected by the virus, other organs, such as the brain, can be targeted by it [5]. It has been shown that SARS-CoV-2 can enter the brain through the olfactory nerve and cause severe neuroinflammation [6]. With the entry of the virus and stimulation of the human immune system, a series of pro-inflammatory cytokines, such as interleukin (IL)-6, IL-1b, IL-8, IL-2R, and tumor necrosis factor (TNF-a), are produced and lead to cytokine storm [7]. In addition, hyperinflammatory reactions cause endothelial dysfunction [8] and pyroptosis [9].
By activation of different immune cells, several regulatory factors are also produced to control immunological events [10]. One of these regulators is MicroRNAs (miRNAs), which are small endogenous RNA molecules and control various physiological processes in cells [11]. Researchers had previously shown that the expression of some miRNAs is associated with respiratory syndromes and viral infections [12, 13]. In the case of COVID-19 disease, some researchers have compiled a list of up/down-regulated miRNAs by both bioinformatics [14,15,16,17] and experimental analysis [18,19,20]. The role of miRNAs in viral infections is very important and they can cause or modulate appropriate immune responses. Interestingly, some miRNAs can bind to viral mRNAs and block the expression of viral genes [21].
In this study, we evaluated the expression of four up/down-regulated inflammatory miRNAs and their mRNA targets in the blood of COVID-19 patients with different grades. We also wanted to see how the relative expression of these miRNAs and their mRNA targets differs during hospitalization. It is important to know the profile of miRNAs in COVID-19 because it can be useful in diagnosis and prognosis. In addition, new molecular drugs can be designed by miRNA profile [22].
Materials and methods
Bioinformatics
DIANA tools v.5.0 (http://diana.imis.athena-innovation.gr) was applied to identify miRNAs potentially related to COVID-19 and inflammatory pathways. Gene targets of the selected miRNAs were detected by miRDB (http://mirdb.org/) [23] and miRTarBase v.8.0 (http://mirtarbase.cuhk.edu.cn) [24], based on a target score. To find signaling pathways of the selected miRNAs, Kyoto Encyclopedia of Genes and Genomes (KEGG) [24] was used by P value (P < 0.05) [25].
COVID-19 patients
In this cross-sectional study, 5 mL of blood sample were taken from COVID-19 patients, obtained from several health centers of Yazd, Tehran, and Zahedan province of Iran. Sampling was from December 20, 2020 to March 2, 2021. We also took blood samples during the patients’ hospitalization to see how the relative expression of miRNAs and mRNAs changed during hospitalization. The age range of patients was 38–63 years old and 52 ± 2% of them was female and 48 ± 2% was male. In this study, COVID-19 patients with grade 1 (n = 21), grade 2 (n = 20), grade 3 (n = 20), grade 4 (n = 21), and grade 5 (n = 21) with informed consent were enrolled in the study. Positive real time-PCR, CT scan (depend on its grade), and clinical symptoms of COVID-19 were important inclusion criteria. Also, COVID-19 patients with severe comorbidities and inflammatory autoimmune diseases were excluded from this study. All COVID-19 patients were treated with remdesivir and favipiravir during hospitalization. Blood samples were also taken from 20 healthy individuals for control. All experiments were under the guidelines of the National Institute of Health, and the ethics committee of Zahedan University of Medical Sciences, Zahedan, Iran. (Ethical Code: IR.ZAUMS.REC.1399.316). After collecting peripheral blood samples, their serum was separated and stored at − 80 °C.
The relative expression of miRNAs and mRNAs
Briefly, total RNAs were extracted from serum samples by mirPremier microRNA isolation kit (Sigma-Aldrich, USA), according to manufacture protocol. After confirmation of RNA purity by NanoDrop ND-1000 UV–VIS spectrophotometer, cDNA was synthesized by Mir-X miRNA First-Strand Synthesis kit (Takara Bio Inc., USA), according to manufacture protocol. Then, real time-PCR was done in presence of Mir-X miRNA qPCR SYBR (Invitrogen, UK). Finally, the threshold cycle (CT) values were recorded for each miRNA and mRNA. The relative expression of miRNAs and mRNAs was calculated by delta-delta CT formula. The reference gene for miRNAs was RNU 48 and the reference gene for mRNAs was GAPDH.
Statistical analysis
All data were reported as the mean ± standard deviation. To find significant differences between groups, a one-way ANOVA was used. Also, Spearman’s correlation coefficient was used to correlate the expression of miRNAs and their mRNA targets. A P value of less than 0.05 was considered statistically significant.
Results
Computational analysis
Five target genes for each miRNA with the highest target score are shown in Table 1, obtained from miRDB and miRTarBase. For hsa-miR-31-3p, the best target gene was zinc finger MYM-type containing 5 (ZMYM5). For hsa-miR-29a-3p, the best target gene was collagen type V alpha 3 chain (COL5A3). For hsa-miR-126-3p, the best target gene was calmodulin-regulated spectrin-associated protein 1 (CAMSAP1). For hsa-miR-17-3p, the best target gene was dicer 1, ribonuclease III (DICER1). Figure 1 shows significant pathways for each miRNA, extracted from KEGG molecular pathway. As seen, important cellular pathways are Wnt signaling pathway and AMPK signaling pathway, PI3K–Akt signaling pathway, mRNA surveillance pathway, MAPK signaling pathway, and mTOR signaling pathway.
The expression of selected miRNAs and their mRNA targets
Figure 2 shows the relative expression of hsa-miR-31-3p (a), hsa-miR-29a-3p (b), hsa-miR-126-3p (c), and hsa-miR-17-3p (d) in COVID-19 patients with grades 1–5. As found, the relative expression of hsa-miR-31-3p, hsa-miR-29a-3p, and hsa-miR-126-3p was significantly decreased with the increase of disease grade. Conversely, the relative expression of hsa-miR-17-3p was significantly increased with the increase of disease grade. These results were in consistent with the relative expression of these miRNAs during hospitalization. That is, in COVID-19 patients who did not respond to treatment (e), the relative expression of hsa-miR-31-3p, hsa-miR-29a-3p, and hsa-miR-126-3p was significantly decreased and the relative expression of hsa-miR-17-3p was significantly increased during hospitalization. Conversely, in COVID-19 patients who respond to treatment (f), the relative expression of hsa-miR-31-3p, hsa-miR-29a-3p, and hsa-miR-126-3p was significantly increased and the relative expression of hsa-miR-17-3p was significantly decreased during hospitalization.
Figure 3 demonstrates the relative expression of ZMYM5 (a), COL5A3 (b), CAMSAP1 (c), DICER1 (d) in COVID-19 patients with grades 1–5. As seen, the relative expression of ZMYM5, COL5A3, and CAMSAP1 was significantly increased with the increase of disease grade. Conversely, the relative expression of DICER1 was significantly decreased with the increase of disease grade. In COVID-19 patients who did not respond to treatment (e), the relative expression of ZMYM5, COL5A3, and CAMSAP1 was significantly increased and the relative expression of DICER1 was significantly decreased during hospitalization. Conversely, in COVID-19 patients who respond to treatment (f), the relative expression of ZMYM5, COL5A3, and CAMSAP1 was significantly decreased and the relative expression of DICER1 was significantly increased during hospitalization. Based on Table 2, a negative significant correlation was seen between the expression of (1) hsa-miR-31-3p and ZMYM5, (2) hsa-miR-29a-3p and COL5A3, (3) hsa-miR-126-3p and CAMSAP1, and (4) hsa-miR-17-3p and DICER1 in COVID-19 patients with any grade (P < 0.05).
Discussion
This study aimed to evaluate the expression of four inflammatory miRNAs and their mRNA targets in the serum samples of COVID-19 patients with different grades. Also, we investigated the relative expression of these miRNAs and their mRNA targets during hospitalization. Since some COVID-19 patients respond to medication (e.g., remdesivir and favipiravir) and some of them do not respond well, the relative expression of miRNAs and their mRNA targets is different in the two groups. The first question is that which miRNAs must be selected and investigated from the human miRNA pool. Although it is a very difficult and complex question, there are two main solutions. The first solution is the evaluation of all human miRNAs with a microarray [26]. With this technique, although we can prepare a complete list of up/down-regulated miRNAs, it requires a complex and expensive device. Another solution is the selection of potential miRNAs using computational and bioinformatics methods. Most of the articles on miRNAs and COVID-19 have followed the same procedure [14,15,16,17]. Generally, miRNAs can be isolated from tissue [18] or blood samples [20]. Importantly, the miRNAs detected in the blood can help us to diagnose and prognoses COVID-19. Arguably, the main source of these miRNAs is from immune and respiratory cells. It is good to say that there are various destructive enzymes in the blood and the desired miRNAs may be rapidly broken down in it [27]. It is interesting to note other types of RNA, such as long non-coding RNAs (long ncRNAs, lncRNA), can also give us information about the presence of SARS-COV-2 [28].
In this study, we found that the relative expression of hsa-miR-31-3p, hsa-miR-29a-3p, and hsa-miR-126-3p was significantly decreased and the relative expression of their mRNA targets (ZMYM5, COL5A3, and CAMSAP1) was significantly increased with the increase of disease grade. Conversely, the relative expression of hsa-miR-17-3p was significantly increased and its mRNA target (DICER1) was significantly decreased with the increase of disease grade. This pattern was exactly seen during hospitalization of COVID-19 patients who did not respond to treatment. In COVID-19 patients who responded to treatment, the expression of selected miRNAs and their mRNA targets returned to the normal level. A negative significant correlation was seen between the expression of (1) hsa-miR-31-3p and ZMYM5, (2) hsa-miR-29a-3p and COL5A3, (3) hsa-miR-126-3p and CAMSAP1, and (4) hsa-miR-17-3p and DICER1 in COVID-19 patients with any grade (P < 0.05) and during hospitalization. Our findings are consistent with the results of other researchers. For example, Farr et al. found that the expression of miR-31 was altered in COVID-19 [29]. Centa et al. [18] showed that hsa-miR-29a-3p is associated with endothelial dysfunction in post-mortem lung biopsies of COVID-19 patients. Widiasta et al. [30] indicated the potential role of ACE2-related microRNAs (such as hsa-miR-126) in COVID-19-associated nephropathy. Li et al. [31] declared that the differential miRNA expression (such as hsa-miR-17-3p) found in COVID-19 patients could regulate the immune system. Based on review literature, all four selected miRNAs have an important role in inflammation. For example, Shi et al. [32] found that miR-31 could mediate inflammatory signaling. Eken et al. [33] indicated that miR-29b could modulate the chronic inflammatory response. Tang et al. [34] showed that miR-126 could suppress inflammation in endothelial cells. Tan et al. [35] found that inhibition of miR-17 reduces the inflammation. Among the limitations of this project was that since only four miRNAs have been analyzed, we missed other important miRNAs. Therefore, it is suggested that in another study, the expression of other inflammatory miRNAs be evaluated in COVID-19 patients.
Conclusions
Taken together, the relative expression of hsa-miR-31-3p, hsa-miR-29a-3p, and hsa-miR-126-3p was down-regulated and the relative expression of their mRNA targets was up-regulated with the increase of COVID-19 grade. This pattern was exactly seen during hospitalization for COVID-19 patients who did not respond to treatment. In COVID-19 patients who responded to treatment, the expression of selected miRNAs and their mRNA targets returned to the normal level. On the other hand, the relative expression of hsa-miR-17-3p was up-regulated and its mRNA target (DICER1) was down-regulated with the increase of COVID-19 grade.
Availability of data and materials
Not applicable.
References
Baud D, Qi X, Nielsen-Saines K, Musso D, Pomar L, Favre G. Real estimates of mortality following COVID-19 infection. Lancet Infect Dis. 2020;20(7):773.
Xu Z, Shi L, Wang Y, Zhang J, Huang L, Zhang C, et al. Pathological findings of COVID-19 associated with acute respiratory distress syndrome. Lancet Respir Med. 2020;8(4):420–2.
Nagashima S, Mendes MC, Camargo Martins AP, Borges NH, Godoy TM, dos Miggiolaro AFRS, et al. Endothelial dysfunction and thrombosis in patients with COVID-19—brief report. Arterioscler Thromb Vasc Biol. 2020;40(10):2404–7.
Hoffmann M, Kleine-Weber H, Schroeder S, Krüger N, Herrler T, Erichsen S, et al. SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor. Cell. 2020;181(2):271-80. e8.
Song E, Zhang C, Israelow B, Lu-Culligan A, Prado AV, Skriabine S, et al. Neuroinvasion of SARS-CoV-2 in human and mouse brain. J Exp Med. 2021;218(3):e20202135.
Yachou Y, El Idrissi A, Belapasov V, Benali SA. Neuroinvasion, neurotropic, and neuroinflammatory events of SARS-CoV-2: understanding the neurological manifestations in COVID-19 patients. Neurol Sci. 2020;10:1–13.
Song P, Li W, Xie J, Hou Y, You C. Cytokine storm induced by SARS-CoV-2. Clin Chim Acta. 2020;509:280.
Pober JS, Sessa WC. Evolving functions of endothelial cells in inflammation. Nat Rev Immunol. 2007;7(10):803–15.
Bergsbaken T, Fink SL, Cookson BT. Pyroptosis: host cell death and inflammation. Nat Rev Microbiol. 2009;7(2):99–109.
Okeke EB, Uzonna JE. The pivotal role of regulatory T cells in the regulation of innate immune cells. Front Immunol. 2019;10:680.
O’connell RM, Rao DS, Chaudhuri AA, Baltimore D. Physiological and pathological roles for microRNAs in the immune system. Nat Rev Immunol. 2010;10(2):111–22.
Bruscella P, Bottini S, Baudesson C, Pawlotsky J-M, Feray C, Trabucchi M. Viruses and miRNAs: more friends than foes. Front Microbiol. 2017;8:824.
Stolzenburg LR, Harris A. The role of microRNAs in chronic respiratory disease: recent insights. Biol Chem. 2018;399(3):219–34.
Bartoszewski R, Dabrowski M, Jakiela B, Matalon S, Harrod KS, Sanak M, et al. SARS-CoV-2 may regulate cellular responses through depletion of specific host miRNAs. Am J PhysiolLung Cell Mol Physiol. 2020;319(3):L444–55.
Khan M, Sany M, Us R, Islam M, Islam ABMM. Epigenetic regulator miRNA pattern differences among SARS-CoV, SARS-CoV-2, and SARS-CoV-2 world-wide isolates delineated the mystery behind the epic pathogenicity and distinct clinical characteristics of pandemic COVID-19. Front Genet. 2020;11:765.
Nersisyan S, Shkurnikov M, Turchinovich A, Knyazev E, Tonevitsky A. Integrative analysis of miRNA and mRNA sequencing data reveals potential regulatory mechanisms of ACE2 and TMPRSS2. PLoS ONE. 2020;15(7):e0235987.
Demirci MDS, Adan A. Computational analysis of microRNA-mediated interactions in SARS-CoV-2 infection. PeerJ. 2020;8:e9369.
Centa A, Fonseca AS, da Silva Ferreira SG, Azevedo MLV, De Paula CBV, Nagashima S, et al. Deregulated miRNA expression is associated with endothelial dysfunction in post-mortem lung biopsies of COVID-19 patients. Am J Physiol Lung Cell Mol Physiol. 2021;320(3):L405–12.
Donyavi T, Bokharaei-Salim F, Baghi HB, Khanaliha K, Janat-Makan MA, Karimi B, et al. Acute and post-acute phase of COVID-19: analyzing expression patterns of miRNA-29a-3p, 146a–3p, 155–5p, and let-7b-3p in PBMC. Int Immunopharmacol. 2021;97:107641.
Pinilla L, Benitez ID, González J, Torres G, Barbé F, de Gonzalo-Calvo D. Peripheral blood microRNAs and the COVID-19 patient: methodological considerations, technical challenges and practice points. RNA Biol. 2021;18(5):688–95.
Alam T, Lipovich L. miRCOVID-19: Potential targets of human miRNAs in SARS-CoV-2 for RNA-based drug discovery. Non Coding RNA. 2021;7(1):18.
Zhang S, Amahong K, Sun X, Lian X, Liu J, Sun H, et al. The miRNA: a small but powerful RNA for COVID-19. Brief Bioinform. 2021;22(2):1137–49.
Wong N, Wang X. miRDB: an online resource for microRNA target prediction and functional annotations. Nucleic Acids Res. 2015;43(D1):D146–52.
Hsu S-D, Lin F-M, Wu W-Y, Liang C, Huang W-C, Chan W-L, et al. miRTarBase: a database curates experimentally validated microRNA–target interactions. Nucleic Acids Res. 2011;39(suppl_1):D163–9.
Kanehisa M, Goto S, Kawashima S, Nakaya A. The KEGG databases at GenomeNet. Nucleic Acids Res. 2002;30(1):42–6.
Lim LP, Lau NC, Garrett-Engele P, Grimson A, Schelter JM, Castle J, et al. Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature. 2005;433(7027):769–73.
Zhang Z, Qin YW, Brewer G, Jing Q. MicroRNA degradation and turnover: regulating the regulators. Wiley Interdiscip Rev RNA. 2012;3(4):593–600.
Wu Y, Zhao T, Deng R, Xia X, Li B, Wang X. A study of differential circRNA and lncRNA expressions in COVID-19-infected peripheral blood. Sci Rep. 2021;11(1):1–14.
Farr R, Rootes C, Rowntree L, Nguyen T, Hensen L, Kedzierski L, et al. Altered microRNA expression in COVID-19 patients enables identification of SARS-CoV-2 infection. 2021. https://doi.org/10.21203/rs.3.rs-253459/v1.
Widiasta A, Sribudiani Y, Nugrahapraja H, Hilmanto D, Sekarwana N, Rachmadi D. Potential role of ACE2-related microRNAs in COVID-19-associated nephropathy. Non Coding RNA Res. 2020;5(4):153–66.
Li C, Hu X, Li L, Li JH. Differential microRNA expression in the peripheral blood from human patients with COVID-19. J Clin Lab Anal. 2020;34(10):e23590.
Shi J, Ma X, Su Y, Song Y, Tian Y, Yuan S, et al. MiR-31 mediates inflammatory signaling to promote re-epithelialization during skin wound healing. J Investig Dermatol. 2018;138(10):2253–63.
Eken SM, Christersdottir T, Winski G, Sangsuwan T, Jin H, Chernogubova E, et al. miR-29b mediates the chronic inflammatory response in radiotherapy-induced vascular disease. JACC Basic Transl Sci. 2019;4(1):72–82.
Tang S-T, Wang F, Shao M, Wang Y, Zhu H-Q. MicroRNA-126 suppresses inflammation in endothelial cells under hyperglycemic condition by targeting HMGB1. Vasc Pharmacol. 2017;88:48–55.
Tan L, Liu L, Jiang Z, Hao X. Inhibition of microRNA-17-5p reduces the inflammation and lipid accumulation, and up-regulates ATP-binding cassette transporterA1 in atherosclerosis. J Pharmacol Sci. 2019;139(4):280–8.
Acknowledgements
We thank the Reference Laboratory of Zahedan University of Medical Sciences. This article was financially supported by Zahedan University of Medical Sciences, Zahedan, Iran (Grant Number: 9936).
Funding
This article was financially supported by Zahedan University of Medical Sciences, Zahedan, Iran (Grant Number: 9936).
Author information
Authors and Affiliations
Contributions
Conception and design: RK and AJ. Administrative support: RK and AJ. Provision of study materials or patients: RK. Collection and assembly of data: AJ. Data analysis and interpretation: RK and AJ. Manuscript writing: all authors. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
All experiments were under the guidelines of the National Institute of Health, the provisions of the Declaration of Helsinki, and the ethics committee of Zahedan University of Medical Sciences, Zahedan, Iran. (Ethical Code: IR.ZAUMS.REC.1399.316). In this study, COVID-19 patients with grade 1 (n = 21), grade 2 (n = 20), grade 3 (n = 20), grade 4 (n = 21), grade 5 (n = 21), and healthy individuals (n = 20) with informed consent were enrolled in the study.
Consent for publication
Not applicable.
Competing interests
There are no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Keikha, R., Hashemi-Shahri, S.M. & Jebali, A. The relative expression of miR-31, miR-29, miR-126, and miR-17 and their mRNA targets in the serum of COVID-19 patients with different grades during hospitalization. Eur J Med Res 26, 75 (2021). https://doi.org/10.1186/s40001-021-00544-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40001-021-00544-4