Cardiology Research and Practice
Volume 2021 (2021), Article ID 6682183, 10 pages
https://doi.org/10.1155/2021/6682183

Research Article

Identifying a Serum Exosomal-Associated lncRNA/circRNA-miRNA-mRNA Network in Coronary Heart Disease

Jia Mao,¹ Yufei Zhou,² Licheng Lu,³ Ping Zhang,¹ Runhua Ren,⁴ Yaqiong Wang,⁴ and Jing Wang⁴

¹Emergency Department, The Affiliated Wuxi No. 2 People’s Hospital of Nanjing Medical University, Wuxi 214000, Jiangsu, China
²Department of Cardiology, The First Affiliated Hospital of Nanjing Medical University, Nanjing 210029, Jiangsu, China
³Department of Cardiology, Kunshan Hospital of Traditional Chinese Medicine, Kunshan 215300, Jiangsu, China
⁴Department of Geriatric Medicine, The Affiliated Jiangning Hospital with Nanjing Medical University, Nanjing 211100, Jiangsu, China

Correspondence should be addressed to Jing Wang; roubaobao09@163.com

Received 1 January 2021; Revised 23 April 2021; Accepted 16 June 2021; Published 24 June 2021

Academic Editor: Domenico Della Rocca

Copyright © 2021 Jia Mao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Background. Accumulating evidence supports the importance of noncoding RNAs and exosomes in coronary heart disease (CHD). However, exosomal-associated competing endogenous RNA- (ceRNA-) mediated regulatory mechanisms in CHD are largely unexplored. The present study aimed to explore exosomal-associated ceRNA networks in CHD. Methods. Data from 6 CHD patients and 32 normal controls were downloaded from the ExoRBase database. CHD and normal controls were compared by screening differentially expressed mRNAs (DEMs), lncRNAs (DELs), and circRNAs (DECs) in serum exosomes. MicroRNAs (miRNAs) targeting DEMs were predicted using the Targetscan and miRanda databases, and miRNAs targeted by DELs and DECs were predicted using the miRcode and starBase databases, respectively. The biological functions and related signaling pathways of DEMs were analyzed using the David and KOBAS databases. Subsequently, a protein-protein interaction (PPI) network was established to screen out on which hub genes enrichment analyses should be performed, and a ceRNA network (lncRNA/circRNA-miRNA-mRNA) was constructed to elucidate ceRNA axes in CHD. Results. A total of 312 DEMs, 43 DELs, and 85 DECs were identified between CHD patients and normal controls. Functional enrichment analysis showed that DEMs were significantly enriched in “chromatin silencing at rDNA,” “telomere organization,” and “negative regulation of gene expression, epigenetic.” PPI network analysis showed that 25 hub DEMs were closely related to CHD, of which ubiquitin C (UBC) was the most important. Hub genes were mainly enriched in “cellular protein metabolic process” functions. The exosomal-associated ceRNA regulatory network incorporated 48 DEMs, 73 predicted miRNAs, 10 DELs, and 15 DECs. The LncRNA/circRNA-miRNA-mRNA interaction axes (RPL7AP11/hsa-miR-17-5p/UBC and RPL7AP11/hsa-miR-20b-5p/UBC) were obtained from the network. Conclusions. Our findings provide a novel perspective on the potential role of exosomal-associated ceRNA network regulation of the pathogenesis of CHD.

1. Introduction

Coronary heart disease (CHD) is a complex biological process accompanied by a wide range of transcriptional changes, but the mechanism of CHD is still complex and unclear [1].

Noncoding RNAs mainly consist of microRNAs (miRNAs/miRs), long noncoding RNAs (lncRNAs), and circular RNAs (circRNAs). MiRNAs are a class of small noncoding RNAs that can block protein translation or induce degradation by targeting specific regions of messenger RNA (mRNA) [2]. LncRNAs have more diverse functions as epigenetic regulators, molecular scaffolds, and decoys [3]. CircRNAs can function as templates for viroid and viral replication, as intermediates in RNA processing, as regulators of transcription, as small nucleolar RNAs, and as miRNA sponges [4]. With the development of sequencing technology and bioinformatics, it has been found that noncoding RNAs are involved in the pathophysiology of cardiovascular diseases [5]. Ahmadi et al. [6] demonstrated that miR-342-5p could be a biomarker for diagnosis of CHD associated with inflammatory cytokines. Wang et al. [7] revealed that BRAF-activated lncRNA is associated with CHD. Moreover, circRNAs or lncRNAs have been found to interact with miRNAs as competing endogenous RNAs (ceRNAs) to regulate target mRNA activity in CHD. For example, YOD1 deubiquitinase might be a novel target for diagnosing CHD from the lncRNA/circRNA-miRNA-mRNA ceRNA network [8].

However, serum RNAs can often be degraded by RNA enzymes and, thus, may not accurately reflect pathological differences; exosomes can protect serum RNAs from degradation [9]. Exosomes are small vesicles with a diameter of approximately 30–150 nm containing proteins, nucleic acids, and lipids [10] and are related to diverse regulatory processes in cardiovascular disorders, including myocardial injury, repair, and regeneration [11].

To better understand the underlying molecular regulatory mechanisms of CHD, we aimed to identify differentially expressed exosomal-related lncRNAs, circRNAs, and miRNAs and used these to construct a ceRNA network to discover accurate and reliable diagnostic biomarkers and therapeutic targets for CHD.

2. Materials and Methods

2.1. Data Collection

A flowchart of the study is shown in Figure 1. The exoRBase database (http://www.exorbase.org/) is a repository of circRNAs, lncRNAs, and mRNAs derived from RNA-seq data, including analysis of human blood exosomes. These samples come from different biological conditions, including normal people (NP), CHD, colorectal cancer, hepatocellular carcinoma, pancreatic adenocarcinoma, and breast cancer [12]. In this study, data from NP and CHD blood samples were downloaded, including 6 patients with CHD and 32 normal controls.

Figure 1: Study flowchart. Abbreviations: CHD, coronary artery disease; mRNAs, messenger RNAs; miRNAs, microRNAs; lncRNAs, long noncoding RNAs; circRNAs, circular RNAs, DEMs, differentially expressed mRNAs; DELs, differentially expressed lncRNAs; DECs, differentially expressed circRNAs; and ceRNA, competing endogenous RNA.

2.2. Identification of Differentially Expressed mRNAs, lncRNAs, and circRNAs

The lists of differentially expressed circRNAs (DECs), lncRNAs (DELs), and mRNAs (DEMs) between controls and patients with CHD were generated using the LIMMA package in R software. The values of |log₂ (fold change (FC))| > 0 and value < 0.05 were selected as cutoff criteria.

2.3. Integration of a PPI Network and Module Analysis

A protein-protein interaction (PPI) network of DEMs was constructed by STRING (https://string‐db.org) and visualized with Cytoscape software [13]. Furthermore, the Molecular Complex Detection (MCODE) application in Cytoscape was used to select the PPI network modules, with cutoff = 2, node score cutoff = 0.2, k-core = 2, and maximum depth = 100 as selection criteria. In addition, nodes with degree ≥5 were identified as hub nodes in the PPI network.

2.4. Functional Enrichment Analyses

Gene ontology (GO) analysis was used to annotate the DEMs and hub genes based on biological processes (BP), cellular components (CC), and molecular functions (MF) [14]. To investigate the biological function of DEMs and hub genes, the database for annotation, visualization, and integrated discovery (DAVID) online tool (version 6.8; http://david.abcc.ncifcrf.gov) was utilized to perform GO analysis [15]. In addition, the KOBAS 3.0 online analysis database was used to perform pathway enrichment analysis [16]. The significant enrichment threshold for GO and KEGG (Kyoto Encyclopedia of Genes and Genomes) analyses was a value <0.05 and count ≥ 2.

2.5. Prediction of microRNAs Targeting mRNAs

The TargetScan (http://www.targetscan.org/vert_72/) [17] and miRanda (http://www.mirdb.org/) [18] databases were used to predict miRNAs that target the DEMs. To increase the accuracy of the predictions, the targeting miRNAs predicted by both databases were used. The miRcode (http://www.mircode.org/) [19] and starBase (http://starbase.sysu.edu.cn/) [20] databases were used to predict miRNAs targeted by the DELs and DECs, respectively. Pairs of miRNAs-mRNAs, lncRNA-miRNA, and circRNA-miRNA were subsequently constructed.

2.6. Construction of the lncRNA/circRNA-miRNA-mRNA ceRNA Network

CeRNA regulation has been reported to serve important roles in human disease; the circRNA or lncRNA-miRNA-mRNA interaction network was constructed to explore the associations among circRNAs, lncRNAs, miRNAs, and mRNAs [21]. Finally, Cytoscape was used to visualize the lncRNA/circRNA-miRNA-mRNA network.

2.7. Statistical Analysis

All data were expressed as the mean ± standard error. Statistical analyses were performed using GraphPad Prism 6 (GraphPad Software, Inc., La Jolla, CA, USA). was considered to indicate a statistically significant difference.

3. Results

3.1. Differential Expression Analysis

A total of 312 DEMs (55 upregulated and 257 downregulated), 43 DELs (24 upregulated and 19 downregulated), and 85 DECs (4 upregulated and 81 downregulated) were identified between the CHD patients and control individuals. The complete upregulated and downregulated DEMs, DELs, and DECs are listed in supplementary materials: Tables S1–S3. Finally, based on the chosen criteria of value <0.05 and | log₂(FC) | > 0.5, we plotted the heat maps for the DEMs, DELs, and DECs, respectively, as shown in Figures 2(a)–2(c).

(a)

(b)

(c)

Figure 2: Hierarchical clustering and heat map analysis of differentially expressed profiles of exosomal RNAs. Note: (a) mRNAs, (b) lncRNAs, and (c) circRNAs. The color scale indicates the expression of differentially expressed exosomal RNAs. Red and blue indicate upand downregulation, respectively. Abbreviations are the same as in Figure 1.

3.2. PPI Network and Module Analyses

To further study the specific DEMs from an interactive perspective, a PPI network was constructed using the STRING database. A statistically significant PPI network consisted of 50 nodes and 189 edges. Nodes with a degree of ≥5 were regarded as hub mRNAs in the network. Two modules were formed in the PPI network with an MCODE score ≥7: module 1 with an MCODE score of 16.875 (nodes = 17) and module 2 with an MCODE score of 7.429 (nodes = 8). Hub mRNAs, namely, ubiquitin C (UBC) and 16 histone cluster family genes, were present in module 1, and cathepsin G, myeloperoxidase, cathelicidin antimicrobial peptide, defensin alpha 1, defensin alpha 3, matrix metallopeptidase 8, azurocidin 1, and defensin alpha 1B were present in module 2 (Figure 3).

Figure 3: Protein-protein interaction network and the significant modules. Note: the red nodes represent the upregulated mRNAs and the green nodes represent the downregulated mRNAs. The most significant module identified by MCODE had a score = 16.875; the second most significant module score = 7.429. Abbreviations: MCODE: molecular complex detection; other abbreviations are the same as in Figure 1.

3.3. Functional Enrichment Analyses

Functional enrichment analyses indicated that the DEMs were mainly enriched in “chromatin silencing at rDNA,” “telomere organization,” and “negative regulation of gene expression, epigenetic” for the BPs. CC analysis showed that the DEMs were significantly enriched for the “nucleosome,” “nuclear chromosome,” and “nuclear chromosome” components. For the MF category, the DEMs were enriched in “histone binding,” “protein heterodimerization activity,” and “poly(A) RNA binding” (Figure 4 and supplementary materials: Table S4). The enrichment analyses for hub genes showed that the biological function of UBC was enriched for a “cellular protein metabolic process,” its location was greatly enriched within the “extracellular exosome” compartment, and its molecular function was “poly(A) RNA binding” (Figure 5 and supplementary materials: Table S5).

Figure 4: The functional terms and pathways enriched for the differentially expressed mRNAs. Note: the red lines represent biological processes, the orange lines represent cellular components, the green lines represent molecular functions, the black lines represent KEGG pathways, and the gray lines represent −log10 ( value). Abbreviation: KEGG, Kyoto Encyclopedia of Genes and Genomes.

Figure 5: Functional and pathway enrichment analyses of hub genes selected from the protein-protein interaction network. Note: the red lines represent biological processes, the orange lines represent cellular components, the green lines represent molecular functions, the black lines represent KEGG pathways, and the gray lines represent −log10 ( value). Abbreviations are the same as in Figure 4.

3.4. LncRNA/circRNA-miRNA-mRNA ceRNA Network

Using the several construction tools described above, 175, 166, and 171 miRNAs were predicted to target DEMs, DELs, and DECs, respectively. The Venn diagram of predicted miRNAs is shown in supplementary materials: Figure S1. Subsequently, ceRNA network analyses were performed to unravel the functions of the identified differentially expressed exosomal-associated ceRNA network in CHD patients. This ceRNA network consisted of 48 DEMs, 73 predicted miRNAs, 10 DELs, and 15 DECs. From the ceRNA network, we identified lncRNA RPL7AP11 as competing for binding to hsa-miR-20a-5p and hsa-miR-17-5p, thereby affecting UBC expression (Figure 6 and supplementary materials: Table S6). These results suggest that the ceRNA networks we predict in this paper might be a key factor underlying the pathogenesis of CHD.

Figure 6: Competing endogenous RNA interaction network. Note: the red circle represents mRNAs, blue triangle represents miRNAs, brown hexagon represents lncRNAs, and purple diamond represents circRNAs. Abbreviations are the same as in Figure 1.

4. Discussion

Atherosclerotic disease and its thrombotic complications may lead to the development of CHD and, if untreated, can progress to myocardial infarction. Exosomes are a type of extracellular vesicle that contain molecular constituents (protein, DNA, and RNA) of the cells that secrete them [22].

Recently, dysregulated expression of RNAs (lncRNAs, circRNAs, miRNAs, and mRNAs) and regulation of their networks have been confirmed to affect the pathogenesis and progression of various tumors, such as endometrial carcinoma [23], ovarian cancer [24], and ependymoma [25]. In the cardiovascular system, several novel serum RNAs have been found to be associated with CHD [8, 26]. Although serum RNAs are often degraded by RNA enzymes and might not accurately reflect pathological differences, exosomes can protect them from degradation [9]. We, therefore, identified serum exosomal-associated RNAs and constructed a ceRNA network associated with CHD, revealing a new targeting axis in its pathogenesis. To our knowledge, this is the first study to explore exosomal-associated ceRNA networks in CHD.

In this study, we first identified 312 DEMs, 85 DECs, and 43 DELs involved in the pathogenesis of coronary heart disease. Enrichment analysis and PPI network construction were subsequently performed and revealed UBC to be one of the most important hub genes. After predicting miRNAs targeting mRNA, an exosomal-associated circRNA/lncRNA-miRNA-mRNA ceRNA network was constructed. Our results suggest specific ceRNA axes in the pathogenesis of CHD that may be promising targets for diagnosis of CHD.

UBC belongs to the ubiquitin family and is associated with protein degradation, DNA repair, kinase modification, autophagy, regulation of inflammation, and cell signaling pathways [27, 28]. Our enrichment analysis of both DEMs and hub genes showed that UBC participates in cellular protein metabolic processes. Ji et al. [29] showed that the expression of ubiquitin was significantly higher in CHD patients than in healthy individuals and the levels of ubiquitin varied with the severity of different classes of CHD. Our study further confirmed the function of UBC in the pathogenesis of CHD, suggesting its potential value as a noninvasive biomarker.

MiR-17-5p has been reported to regulate the cell cycle, proliferation, and apoptosis, and broad support has been provided for its role in regulation of cardiovascular diseases. Deletion of miR17 in neonatal mice is lethal, and overexpression of miR-17-5p can extend the life span of mice [30]. Liu et al. [31] showed that upregulation of miR-17-5p could contribute to hypoxia-induced proliferation of human pulmonary artery smooth muscle cells, leading to pulmonary hypertension. Yang et al. [32] found that miR-17-5p silencing protects heart function after acute myocardial infarction by decreasing the rate of apoptosis and repairing vascular injury. Moreover, recent studies have shown that circulating miR-17-5p could be a novel biomarker for diagnosis of acute myocardial infarction [33].

MiR-20b-5p has been shown to attenuate hypoxia-induced apoptosis in cardiomyocytes [34]. Also, Zhen et al. [35] found that overexpression of miR-20b-5p could increase cell viability and repress autophagy and apoptosis in human umbilical vein endothelial cells that had experienced hypoxia-reoxygenation injury. Since both hypoxia and hypoxia-reoxygenation models are similar to patients with myocardial infarction and subsequent revascularization, miR-20b-5p may play a role in regulating CHD. However, there has been little research on the function of miR-20b-5p in CHD patients, and confirmation of its role requires further research.

The lncRNA RPL7AP11 (ribosomal protein L7a pseudogene 11) is a pseudogene of ribosomal protein L7a (RPL7A). Zhang et al. [36] found that RPL7A was downregulated in a vascular endothelial cell line (ECV 304) induced by high-density lipoprotein. Pseudogenes, abundant in the human genome, were traditionally considered to be nonfunctional “junk” genes [37], although recent studies have demonstrated their role in various diseases. However, there has been limited research on RPL7AP11, and more evidence is needed.

Our study has shown that RPL7AP11 can sponge hsa-miR-17-5p and hsa-miR-20b-5p to upregulate UBC, thus regulating the pathogenesis of CHD through cellular protein metabolism.

There are several limitations to the present study. Firstly, the sample was not large. An additional validation cohort should be included in future studies to analyze the expression of these identified lncRNAs, circRNAs, miRNAs, and mRNAs. Secondly, how these novel exosomal-associated ceRNA axes participate in the development of CHD is still unclear. Further cell and animal experiments are needed to verify these findings. Moreover, further studies should focus on exploring the ceRNA networks between lncRNA/circRNA in serum exosomes and miRNA-mRNA in target recipient cells.

5. Conclusions

In conclusion, our comprehensive study identified several exosomal-associated lncRNA/circRNA-miRNA-mRNA interaction axes (RPL7AP11/hsa-miR-17-5p/UBC and RPL7AP11/hsa-miR-20b-5p/UBC) in the progression of CHD, which may be crucial targets for disease treatment.

Data Availability

The authors confirm that all data underlying the findings are fully available without restriction. All relevant data are accessible from the exoRBase database (http://www.exorbase.org/). Processed data are available from the corresponding author Jing Wang (e-mail: roubaobao09@163.com).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

JW, JM, and YFZ conceptualized and designed the study; JW and LCL acquired data; JM and PZ analysed and interpreted the data; RHR, YQW, and JW drafted the manuscript; and all authors approved the final version of the manuscript. Jia Mao, Yufei Zhou, and Licheng Lu contributed equally to this work.

Acknowledgments

The authors would like to show their gratitude to Research Square in which their manuscript was posted as the preprint version (https://www.researchsquare.com/article/rs-131078/v1). They would like to thank Dr. Yulu Hong (Department of Computer Science and Technology, Central South University, Changsha, Hunan, China), an expert in computer science, for the check of the database they used.

Supplementary Materials

Table S1: all upregulated and downregulated DEMs between the normal controls and the CHD samples. Table S2: all upregulated and downregulated DELs between the normal controls and the CHD samples. Table S3: all upregulated and downregulated DECs between the normal controls and the CHD samples. Table S4: enrichment analyses of differentially expressed mRNAs. Table S5: enrichment analyses of module hub genes. Table S6: complete list of competing endogenous RNA network pairs. Figure S1: Venn diagram of predicted microRNAs targeting DEMs, DELs, and DECs. Abbreviations: DEMs, differentially expressed mRNAs; DELs, differentially expressed lncRNAs; and DECs, differentially expressed circRNAs. (Supplementary Materials)

References

J. E. Dalen, J. S. Alpert, R. J. Goldberg, and R. S. Weinstein, “The epidemic of the 20th century: coronary heart disease,” The American Journal of Medicine, vol. 127, no. 9, pp. 807–812, 2014.

K. Saliminejad, H. R. Khorram Khorshid, S. Soleymani Fard, and S. H. Ghaffari, “An overview of microRNAs: biology, functions, therapeutics, and analysis methods,” Journal of Cellular Physiology, vol. 234, pp. 5451–5465, 2019.

J. Zhu, H. Fu, Y. Wu, and X. Zheng, “Function of lncRNAs and approaches to lncRNA-protein interactions,” Science China Life Sciences, vol. 56, pp. 876–885, 2013.

E. Lasda and R. Parker, “Circular RNAs: diversity of form and function,” RNA, vol. 20, pp. 1829–1842, 2014.

C. Jaquenod De Giusti, M. Santalla, and S. Das, “Exosomal non-coding RNAs (Exo-ncRNAs) in cardiovascular health,” Journal of Molecular and Cellular Cardiology, vol. 137, pp. 143–151, 2019.

R. Ahmadi, E. Heidarian, R. Fadaei, N. Moradi, M. Malek, and S. Fallah, “miR-342-5p expression levels in coronary artery disease patients and its association with inflammatory cytokines,” Clinical Laboratory, vol. 64, pp. 603–609, 2018.

H. Wang, N. Zhang, G. Li, and B. Xu, “Proinflammatory cytokine IFN-γ, lncRNA BANCR and the occurrence of coronary artery disease,” Life Sciences, vol. 231, Article ID 116510, 2019.

L. Li, L. Wang, H. Li et al., “Characterization of LncRNA expression profile and identification of novel LncRNA biomarkers to diagnose coronary artery disease,” Atherosclerosis, vol. 275, pp. 359–367, 2018.

Q. Ge, Y. Zhou, J. Lu, Y. Bai, X. Xie, and Z. Lu, “miRNA in plasma exosome is stable under different storage conditions,” Molecules, vol. 19, pp. 1568–1575, 2014.

A. T. Jan, S. Rahman, S. Khan, S. A. Tasduq, and I. Choi, “Biology, pathophysiological role, and clinical implications of exosomes: a critical appraisal,” Cells, vol. 8, 2019.

Y. Zhang, Y. W. Hu, L. Zheng, and Q. Wang, “Characteristics and roles of exosomes in cardiovascular disease,” DNA and Cell Biology, vol. 36, pp. 202–211, 2017.

S. Li, Y. Li, B. Chen et al., “exoRBase: a database of circRNA, lncRNA and mRNA in human blood exosomes,” Nucleic Acids Research, vol. 46, pp. D106–D112, 2018.

P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, pp. 2498–2504, 2003.

M. Ashburner, C. A. Ball, J. A. Blake et al., “Gene Ontology: tool for the unification of biology,” Nature Genetics, vol. 25, no. 1, pp. 25–29, 2000.

D. W. Huang, B. T. Sherman, Q. Tan et al., “The DAVID Gene Functional Classification Tool: a novel biological module-centric algorithm to functionally analyze large gene lists,” Genome Biology, vol. 8, p. R183, 2007.

C. Xie, X. Mao, J. Huang et al., “KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases,” Nucleic Acids Research, vol. 39, pp. W316–W322, 2011.

V. Agarwal, G. W. Bell, J. W. Nam, and D. P. Bartel, “Predicting effective microRNA target sites in mammalian mRNAs,” Elife, vol. 4, 2015.

Y. Chen and X. Wang, “miRDB: an online database for prediction of functional microRNA targets,” Nucleic Acids Research, vol. 48, pp. D127–d31, 2020.

A. Jeggari, D. S. Marks, and E. Larsson, “miRcode: a map of putative microRNA target sites in the long non-coding transcriptome,” Bioinformatics, vol. 28, pp. 2062-2063, 2012.

J.-H. Yang, J.-H. Li, P. Shao, H. Zhou, Y.-Q. Chen, and L.-H. Qu, “starBase: a database for exploring microRNA-mRNA interaction maps from Argonaute CLIP-Seq and Degradome-Seq data,” Nucleic Acids Research, vol. 39, pp. D202–D209, 2011.

L. Salmena, L. Poliseno, Y. Tay, L. Kats, and P. P. Pandolfi, “A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language?” Cell, vol. 146, pp. 353–358, 2011.

R. Kalluri and V. S. Lebleu, “The biology, function, and biomedical applications of exosomes,” Science, vol. 367, 2020.

C. Liu, Y. H. Zhang, Q. Deng et al., “Cancer-related triplets of mRNA-lncRNA-miRNA revealed by integrative network in uterine corpus endometrial carcinoma,” BioMed Research International, vol. 2017, Article ID 3859582, 7 pages, 2017.

Y. Zhou, X. Zheng, B. Xu, W. Hu, T. Huang, and J. Jiang, “The identification and analysis of mRNA-lncRNA-miRNA cliques from the integrative network of ovarian cancer,” Frontiers in Genetics, vol. 10, p. 751, 2019.

F. Liu, H. Dong, Z. Mei, and T. Huang, “Investigation of miRNA and mRNA Co-expression network in ependymoma,” Frontiers in Bioengineering and Biotechnology, vol. 8, p. 177, 2020.

L. Wang, C. Shen, Y. Wang et al., “Identification of circular RNA Hsa_circ_0001879 and Hsa_circ_0004104 as novel biomarkers for coronary artery disease,” Atherosclerosis, vol. 286, pp. 88–96, 2019.

G. Lopez-Castejon, “Control of the inflammasome by the ubiquitin system,” The FEBS Journal, vol. 287, pp. 11–26, 2020.

P. Grumati and I. Dikic, “Ubiquitin signaling and autophagy,” Journal of Biological Chemistry, vol. 293, pp. 5404–5413, 2018.

Y. Ji, J. Yao, Y. Zhao, J. Zhai, Z. Weng, and Y. He, “Extracellular ubiquitin levels are increased in coronary heart disease and associated with the severity of the disease,” Scandinavian Journal of Clinical and Laboratory Investigation, vol. 80, pp. 256–264, 2020.

H. Dellago, M. R. Bobbili, and J. Grillari, “MicroRNA-17-5p: at the crossroads of cancer and aging - a mini-review,” Gerontology, vol. 63, pp. 20–28, 2017.

G. Liu, P. Hao, J. Xu et al., “Upregulation of microRNA-17-5p contributes to hypoxia-induced proliferation in human pulmonary artery smooth muscle cells through modulation of p21 and PTEN,” Respiratory Research, vol. 19, p. 200, 2018.

S. Yang, T. Fan, Q. Hu et al., “Downregulation of microRNA-17-5p improves cardiac function after myocardial infarction via attenuation of apoptosis in endothelial cells,” Molecular Genetics and Genomics, vol. 293, pp. 883–894, 2018.

S. Xue, D. Liu, W. Zhu et al., “Circulating MiR-17-5p, MiR-126-5p and MiR-145-3p are novel biomarkers for diagnosis of acute myocardial infarction,” Frontiers in Physiology, vol. 10, p. 123, 2019.

Z. Zhou, S. Chen, Z. Tian et al., “miR-20b-5p attenuates hypoxia-induced apoptosis in cardiomyocytes via the HIF-1α/NF-κB pathway,” Acta Biochimica et Biophysica Sinica, vol. 52, pp. 927–934, 2020.

W. Zhen, D. Hui, S. Wenying, and S. Yulong, “microRNA-20b-5p regulates propofol-preconditioning-induced inhibition of autophagy in hypoxia-and-reoxygenation-stimulated endothelial cells,” Journal of Biosciences, vol. 45, 2020.

W. Zhang, B. Chen, W. Zeng, L. Bai, H. Xue, and G. Wu, “[Differentially expressed genes in vascular endothelial cell line ECV 304 induced by high density lipoprotein],” Zhonghua Yi Xue Za Zhi, vol. 82, pp. 921–923, 2002.

W. Lou, B. Ding, and P. Fu, “Pseudogene-derived lncRNAs and their miRNA sponging mechanism in human cancer,” Frontiers in Cell and Developmental Biology, vol. 8, p. 85, 2020.