Disease Markers
Volume 2020 (2020), Article ID 3178642, 7 pages
https://doi.org/10.1155/2020/3178642

Research Article

Plasma-Derived Exosomal Circular RNA hsa_circ_0005540 as a Novel Diagnostic Biomarker for Coronary Artery Disease

Wei-peng Wu,1,2 Yan-hong Pan,1,2 Meng-yun Cai,¹ Jin-ming Cen,³ Can Chen,⁴ Lei Zheng,⁵ Xinguang Liu,1,2 and Xing-dong Xiong1,2

¹Guangdong Provincial Key Laboratory of Medical Molecular Diagnostics, Institute of Aging Research, Guangdong Medical University, Dongguan 523808, China
²Institute of Biochemistry & Molecular Biology, Guangdong Medical University, Zhanjiang 524023, China
³Department of Cardiovascular Disease, The First People’s Hospital of Foshan, Foshan 528000, China
⁴Department of Cardiovascular Disease, The Affiliated Hospital of Guangdong Medical University, Zhanjiang 524023, China
⁵Department of Laboratory Medicine, Nanfang Hospital, Southern Medical University, Guangzhou 510515, China

Correspondence should be addressed to Xing-dong Xiong; xiongxingdong@126.com

Received 18 November 2019; Revised 16 June 2020; Accepted 19 June 2020; Published 30 June 2020

Academic Editor: Wen-Jun Tu

Copyright © 2020 Wei-peng Wu 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. Exosomes exist in almost all body fluid and contain diverse biological contents which may be reflective of disease state. Circular RNAs (circRNAs) are stable in structure and have a long half-life in exosomes without degradation, thus making them reliable biomarkers. However, the potential of exosomal circRNAs as biomarkers of coronary artery disease (CAD) remains to be established. Here, we aimed to investigate the expression levels and the potential use of exosomal circRNAs as diagnostic biomarkers for CAD. Methods. CircRNA expression levels in exosomes obtained from three plasma samples of CAD patients and three paired controls were analyzed using RNA sequencing. Exosomal circRNAs obtained in the profiling phase were then verified in two-center validation cohorts. Finally, the ability of exosomal circRNAs, adjusting for Framingham Heart Study (FHS) risk factors, was determined to discriminate between CAD patients and non-CAD controls. Results. 355 circRNAs were differentially expressed between these two groups: 164 were upregulated, and 191 were downregulated. Here, we selected the potential circRNAs (, ) as candidate biomarkers for further validation. Our data showed that only hsa_circ_0005540 was significantly associated with CAD (). After adjustment for risk factors, hsa_circ_0005540 showed a high discriminatory power for CAD in ROC analyses (; , ). Conclusion. Our results suggest that plasma exosomal hsa_circ_0005540 can be used as a promising diagnostic biomarker of CAD.

1. Introduction

CAD is one of the most common forms of cardiovascular diseases and remains a leading cause of mortality worldwide. It poses a high economic burden to the national healthcare system mainly due to its late diagnosis and severe complications [1, 2]. Among the current clinical diagnosis methods, invasive coronary angiography (CAG) is the diagnostic “gold standard,” but its specialized technology and high cost limit it to a selective population [3]. Although the current standardized treatment strategies for CAD, which include coronary bypass surgery, balloon dilatation of coronary vessels, and percutaneous transluminal coronary angioplasty, have been developed, the survival of the patients has remained dismal and the prognosis ineligible [4, 5]. The main reason accounting for this discrepancy is mainly due to the current diagnostic methods unable to simultaneously achieve high sensitivity and convenience. Thus, highly sensitive and convenient diagnostic biomarkers are urgently needed for diagnosis of CAD.

Exosomes are a class of nanoscale (30-150 nm) extracellular vesicles that can be generated and released by most cell types [6]. They are generated within the endosomal compartment and released into the extracellular environment after the fusion of multivesicular bodies (MVBs) with the plasma membrane [7–9]. These vesicles have been proved to be accessible in nearly all body fluids such as the blood [10], urine [11, 12], and saliva [13]. Exosomes can transfer functional RNAs, proteins, and lipids to adjacent cells and serve as mediators of intercellular communication [14–16]. In addition, exosomes have a relatively stable structure which can protect their cargos from destruction, thus making them potential candidates for noninvasive biomarkers [17].

CircRNAs are a large novel type of endogenous transcripts and represent a recent research hotspot in the field of RNA. They form covalently closed loop structures with neither 5-3 polarities nor a polyadenylated tail, which makes them much more stable than linear RNAs [18–20]. Recent studies showed that circRNAs could function as miRNA sponges, bind partners of proteins, regulators of alternative splicing, or even be translated into proteins driven by N6-methyladenosine modification [21–25]. Additionally, high abundance, relative stability, tissue/developmental-stage-specific expression, and evolutionary conservation among species make circRNAs very attractive for clinical application [18, 26, 27]. Due to those characteristics, circRNAs have a great potential to serve as significant biomarkers to diagnose diseases.

Recent studies confirmed that circRNAs are enriched and stable in exosomes and can continually play their roles after exosomes are secreted from the parent cells into the bodily fluid [28, 29]. There are growing evidences of exosomal circRNA biomarkers for various cancers [29–31]. For instance, it was reported that serum exosomal circKLDHC10 could distinguish patients with colorectal cancer from healthy controls, suggesting its significant potential as a circulating biomarker for cancer diagnosis [29]. However, the functional roles of exosomal circRNAs in CAD remain elusive. Thus, the current study is aimed at assessing the differential expressions of exosomal circRNAs in the plasma of CAD patients and non-CAD controls and to determine their potential diagnostic value for CAD.

2. Materials and Methods

2.1. Study Subjects

Sixty-one CAD patients and thirty-eight non-CAD controls enrolled from center 1 (the First People’s Hospital of Foshan) between March 2018 and August 2018 formed the profiling and internal validation phases. Forty-seven CAD patients and 51 non-CAD controls recruited at center 2 (the Affiliated Hospital of Guangdong Medical University) between January 2019 and May 2019 constituted the external validation phase. The diagnosis of CAD was based on the percentage narrowing of each coronary artery segment, with 50% stenosis of at least one of the major coronaries (the left main coronary trunk, anterior descending branch, circumflex artery, and right coronary artery). Patients with impaired ejection fraction, heart failure, or acute myocardial injury were excluded. Age-, sex-, and ethnicity-matched subjects were used as non-CAD controls. All the controls have normal electrocardiogram (ECG) records and no evidence of CAD. General exclusion criteria were a known history of leukopenia, thrombocytopenia, or severe hepatic or renal dysfunction, as well as evidence for inflammatory or malignant disease. Written informed consent was obtained from each of the participants. This study was approved by the Ethics Committee of the First People’s Hospital of Foshan and the Affiliated Hospital of Guangdong Medical University under the guidance of the Helsinki Declaration. The characteristics of the study subjects are shown in Table 1.

Table 1: Characteristics of the study populations.

Characteristics
Profiling phase
Internal validation phase
External validation phase

Non-CAD ()
CAD ()
Non-CAD ()
CAD ()
Non-CAD ()
CAD ()

Age (years)

Sex (male)
2 (66.67%)
1 (33.33%)
13 (37.14%)
36 (62.06%)
19 (37.25%)
35 (74.47%)

Smoking
0 (0%)
1 (33.33%)
5 (14.28%)
26 (44.83%)
0 (0%)
8 (17.02%)

Hypertension
1 (33.33%)
3 (100.00%)
12 (34.29%)
37 (63.79%)
9 (17.65%)
16 (34.04%)

Diabetes
0 (0%)
2 (66.67%)
1 (2.86%)
18 (31.03%)
5 (9.80%)
16 (34.04%)

Systolic BP (mmHg)

FPG (mM)

TC (mM)

HDL (mM)

Coronary artery disease

    1 vessel

1 (33.33%)

25 (43.10%)

17 (36.17%)

    2 vessels

1 (33.33%)

24 (41.38%)

19 (40.43%)

    3 vessels

1 (33.33%)

9 (15.52%)

11 (23.40%)

Data are summarized by either . Abbreviations: BP: blood pressure; FPG: fasting plasma glucose; TC: total cholesterol; HDL: high-density lipoprotein.

2.2. Exosomal RNA Isolation

For exosomal RNA extraction, an exoRNeasy Serum/Plasma Midi kit (Qiagen, CA, USA) was used to isolate total exosome RNA from 1 ml prefiltered plasma for validation cohorts according to the manufacturer’s protocol. Briefly, prefiltered plasma was mixed 1 : 1 with 2x binding buffer (XBP) and added to the exoEasy membrane affinity column to bind the exosomes to the membrane. After centrifugation, the flow-through was discarded and the wash buffer (XWP) was added to the column to wash off nonspecifically retained material. After another centrifugation and discarding of the flow-through, the plasma exosomes were lysed by adding QIAzol to the spin column, and the lysate was collected by centrifugation. Following the addition of chloroform, thorough mixing and centrifugation to separate organic and aqueous phases, the aqueous phase was recovered and mixed with ethanol. The sample-ethanol mixture was added to the RNeasy MinElute spin column and centrifuged. The column was washed once with buffer RWT, and then twice with buffer RPE followed by elution of RNA in water. Collected RNA was then used in the subsequent experiments.

2.3. Exosomal CircRNA Sequence Analysis

Total RNA for RNA-seq analysis was extracted using an exoRNeasy Serum/Plasma Maxi kit (Qiagen, CA, USA) from 4 ml prefiltered plasma according to the manufacturer’s protocol. The purity and concentration of RNA were determined from the OD260/280 readings using a spectrophotometer (NanoDrop ND-1000). The RNA integrity was assessed using Agilent Bioanalyzer 2100 (Agilent Technologies). Total RNA from each sample was subjected to the RiboMinus Eukaryote Kit (Qiagen, CA, USA) to remove ribosomal RNA before the construction of RNA-seq libraries. Strand-specific RNA-seq libraries were prepared using NEBNext Ultra Directional RNA Library Prep Kit for Illumina (NEB) following the manufacturer’s instructions. Briefly, approximately 50 ng of ribosome-depleted RNA samples were fragmented and then used for first- and second-strand cDNA syntheses with random hexamer primers. For second-strand cDNA synthesis, dUTP mix (without dTTP) was used which allows the removal of the second strand. An End-It DNA End Repair Kit was used to repair the ends of the double-stranded cDNA fragments, which were then modified by the Klenow fragment so that an A was added to the 3 end of the DNA fragments; the fragments were finally ligated to adapters. The ligated cDNA products were purified and treated with uracil DNA glycosylase to remove the second-strand cDNA. Purified first-strand cDNA was subjected to 14 cycles of PCR amplification, followed by library analysis with a Bioanalyser 2100 (Agilent, CA, USA). The strand-specific RNA-seq libraries were sequenced using Illumina HiSeq™ 2500 platform.

2.4. Quantitative Real-Time RT-PCR

Total RNAs were converted to cDNA using PrimeScript™ RT Reagent Kit (Takara, Otsu, Japan). The differentially expressed exosomal circRNAs selected from the profiling data were measured by quantitative real-time RT-PCR using the LightCycler® 96 detection system at 95°C for 3 min and amplified by 40 cycles of denaturing at 95°C for 10 s and 60°C for 30 s. The CT value was the fractional cycle number at which the fluorescence exceeded the given threshold. ACTIN was used to normalize the RNA preparation. The relative fold-change was calculated using the 2-ΔΔCT method. Primer sequences are listed in TABLE S1 (see TABLE S1 in the Supplementary Material).

2.5. Statistical Analysis

The differences of the demographic, clinical pathological characteristics and circRNA expression levels between the CAD patients and non-CAD controls were estimated using Student’s -test and test for continuous variables and categorical variables, respectively. Data were expressed as percentages for the qualitative variables and (SD) for the quantitative variables. In the scatter plot depicting the circRNA expression, the horizontal lines represent the medians. The area under the receiver operator characteristic (ROC) curve (AUC) was used to evaluate the diagnostic value of exosomal circRNAs for CAD. The corresponding sensitivity and specificity were identified through the ROC curve analysis. Multivariate logistic regression models were constructed with stepwise addition of FHS risk factors (age, gender, smoking history, diabetes, hypertension, and TC and HDL levels) for each circRNA. All statistical analyses were performed using SPSS software. was considered statistically significant.

3. Results

3.1. Profiling of Exosomal CircRNAs in CAD Patients

To identify alterations in circRNA expression in circulating exosomes between CAD patients and non-CAD controls, we first performed circRNA screening. Specifically, the circRNA profiles in isolated circulating exosomes from CAD patients () and non-CAD controls () were assessed using an Illumina HiSeq™ 2500. As illustrated in the hierarchical clustering analysis (Figure 1(a)), circRNA expression profiles in CAD patients were distinctly different from those in non-CAD controls. The volcano plot showed the statistical significance of differentially expressed circRNAs between CAD patients and non-CAD controls and identified 355 circRNAs whose levels changed substantially (), including 164 upregulated and 191 downregulated circRNAs (Figure 1(b)). We then verified the expression of hsa_circ_0005540, hsa_circ_0000676, and hsa_circ_0007385, which were differentially expressed in the profile phase (, ).

(a)

(b)

Figure 1: Differentially expressed exosomal circRNAs between CAD patients and non-CAD controls. (a) Hierarchical clustering analysis of exosomal circRNA profile in CAD patients () and non-CAD controls (). The expression of exosomal circRNAs is hierarchically clustered on the -axis; CAD patients and non-CAD controls are hierarchically clustered on the -axis. Expression values are presented in red and blue to indicate upregulation and downregulation, respectively. (b) Volcano plot: -axis: log₂(fold change); -axis: -log₁₀( value). The red points indicate exosomal circRNAs 1.5-fold upregulated significantly and the green points represent exosomal circRNAs 1.5-fold downregulated with statistical significance.

3.2. Internal Validation of Differentially Expressed Exosomal CircRNAs

To technically and biologically validate these exosomal circRNAs from the screening stage, we performed quantitative real-time RT-PCR measurements in 58 CAD patients and 35 non-CAD controls, referred to as internal validation cohort. The results showed that hsa_circ_0005540 and hsa_circ_0000676 were significantly increased in CAD patients as compared with non-CAD controls (, Figures 2(a) and 2(b)). However, there was no significant difference for hsa_circ_0007385 between CAD patients and non-CAD controls (, Figure 2(c)).

(a)

(b)

(c)

Figure 2: Expression levels of exosomal circRNAs were quantified by qPCR in the internal validation cohort. (a–c) Expression levels of exosomal circRNAs: (a) hsa_circ_0005540, (b) hsa_circ_0000676, and (c) hsa_circ_0007385.

3.3. External Validation of Differentially Expressed Exosomal CircRNAs

To further explore the applicability of these exosomal circRNAs as potential diagnostic biomarkers for CAD, we tested these circRNAs in an external validation cohort consisting of 47 CAD patients and 51 non-CAD controls. As shown in Figure 3(a), only hsa_circ_0005540 was significantly increased in patients with CAD as compared with non-CAD controls (), whereas there was no significant difference for hsa_circ_0000676 or hsa_circ_0007385 between the two groups (Figures 3(b) and 3(c)).

(a)

(b)

(c)

Figure 3: Expression levels of exosomal circRNAs were quantified by qPCR in the external validation cohort. (a–c) Expression levels of exosomal circRNAs: (a) hsa_circ_0005540, (b) hsa_circ_0000676, and (c) hsa_circ_0007385.

3.4. Exosomal hsa_circ_0005540 Combined with FHS Risk Factors

As the above results were obtained in two different cohorts, we further assessed hsa_circ_0005540 as a biomarker for CAD in all individuals, including 105 CAD patients and 86 non-CAD controls. Levels of hsa_circ_0005540 were significantly increased in CAD patients as compared with the non-CAD controls (, Figure 4(a)). Then, we examined the correlation between hsa_circ_0005540 levels and the FHS risk factors by Pearson’s correlation test. Our data revealed that hsa_circ_0005540 were weakly correlated with the FHS risk factors (, , Supplementary Figure S1). After adjustment for risk factors, hsa_circ_0005540 showed a high discriminatory power for CAD in ROC analyses (; -0.906, , Figure 4(b)), with the sensitivity of 0.810 and the specificity of 0.765. These results imply that the exosomal hsa_circ_0005540 could serve as a potential diagnostic biomarker for CAD.

(a)

(b)

Figure 4: Expression level and ROC curve analysis of exosomal hsa_circ_0005540 in two-center validation cohorts. (a) Expression level of hsa_circ_0005540. (b) ROC curve analysis of exosomal hsa_circ_0005540 combined with FHS RFs for diagnosis of CAD.

4. Discussion

Currently, CAD can be diagnosed by either noninvasive or invasive methods, but each has its limitations. Therefore, a highly sensitive and specific biomarker is urgently needed to facilitate the diagnosis of CAD. In this study, we first explored the differentially expressed exosomal circRNAs in CAD patients and non-CAD controls, and potential circRNAs were selected for verification. Our results from clinical samples demonstrated that the levels of exosomal hsa_circ_0005540 in plasma from patients with CAD were significantly elevated compared with those non-CAD controls. The present work has led us to conclude that exosomal circRNAs can be clinically practicable biomarkers for CAD diagnosis.

To date, clinical studies of circRNAs expression in CAD are limited and have produced varying results. Several cross-sectional studies have reported differences in some circRNAs among patients with CAD compared to controls [32–34]. A study from Zhao et al. was the first to investigate the circRNA profile in the peripheral blood of CAD patients and determine its correlation with the severity of CAD, which suggests that hsa_circ_0124644 might be a sensitive and specific biomarker for diagnosing CAD [34]. Another study offered a transcriptome-wide overview of aberrantly expressed circRNAs in CAD patients and identified hsa_circ_0001879 and hsa_circ_0004104 as novel circRNA biomarkers to diagnose CAD [33]. However, no exosomal circRNAs have been reported to be associated with CAD. The aim of this study was to investigate the potential of exosomal circRNAs in plasma as diagnostic biomarkers for CAD. Our results suggest that plasma exosomal hsa_circ_0005540 can be used as a promising diagnostic biomarker for CAD.

hsa_circ_0005540 is located at chr5:94204037-94248681, and the name of its source gene is MCTP1, which encodes multiple C2 domain transmembrane protein 1. There is currently no definitive evidence demonstrating the biological function of hsa_circ_0005540. CircRNAs can regulate gene expression by acting as miRNA sponges. There are some binding sites for miR-221 and miR-145 in hsa_circ_0005540, and these miRNAs were associated with CAD [35–38]. Moreover, miR-221 was closely related to endothelial cell survival, migration, and capillary tube formation [35]. miR-145 was explored for its potential as a powerful biomarker for diagnosis of CAD and had been shown to be associated with the severity of CAD [36, 38]. Therefore, we speculate that hsa_circ_0005540 is involved in the progression of CAD.

In addition to CAG, many methods are currently used to diagnose CAD, including routine ECG, Holter monitoring, TET, and CTA. The sensitivities of these methods have been shown to be 0.29, 0.65, 0.79, and 0.92, respectively [39–42]. The specificities are 0.67, 0.90, 0.80, and 0.75, respectively [39–42]. In our study, we found that the sensitivity and specificity of exosomal hsa_circ_0005540 combined with FHS risk factors were 0.810 and 0.765, respectively. Based on these comparisons, the diagnostic value of exosomal hsa_circ_0005540 combined with FHS risk factors was greater than that of routine ECG and TET, which was approximately equal to Holter monitoring. When considering the cost and convenience of diagnostic methods, exosomal hsa_circ_0005540 might ameliorate the diagnosis of CAD.

It should be pointed out that there were some limitations in our study. Although 2 independent cohorts for validation had been performed in our study, a larger clinical population is warranted to further validate the performance of the association of exosomal hsa_circ_0005540 as a potential blood-based signature in CAD. Furthermore, since the levels of exosomal circRNAs may be affected by multiple parameters such as the change in expression in the tissue and the release of exosomes by the cells into the circulation, further studies are necessary to explore the mechanisms underlying the dysregulation and the putative impact of the changes of exosomal hsa_circ_0005540 level in the pathophysiology.

5. Conclusion

In conclusion, our study reported a novel signature of exosomal hsa_circ_0005540, allowing to distinguish patients with CAD from non-CAD controls, and provided first insights into the levels of exosomal circRNAs in patients with CAD that could be translated into noninvasive blood-based biomarker.

Data Availability

Raw and normalized data files for the RNA-seq analysis have been deposited in the NCBI Gene Expression Omnibus under accession number GSE152498. You may view the GSE152498 study at https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE152498.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

Wei-peng Wu and Yan-hong Pan contributed equally to this work.

Acknowledgments

The study was supported by grants from the National Natural Science Foundation of China (81871120), the Natural Science Foundation of Guangdong Province (2019A1515010334, 2017KZDXM039), the Yangfan Training Program of Guangdong Province (4YF16006G), and the Construction of Basic Medical Disciplines in Guangdong Medical University (4SG19047G).

Supplementary Materials

Table S1: sequences of all real-time PCR primers. Figure S1: the Pearson correlation coefficient between hsa_circ_0005540 and FHS risk factors. (Supplementary Materials)

References

E. J. Benjamin, P. Muntner, A. Alonso et al., “Heart disease and stroke statistics-2019 update: a report from the American Heart Association,” Circulation, vol. 139, no. 10, pp. e56–528, 2019.

Y. Hao, J. Liu, J. Liu et al., “Sex differences in in-hospital management and outcomes of patients with acute coronary syndrome,” Circulation, vol. 139, no. 15, pp. 1776–1785, 2019.

S. Achenbach and W. G. Daniel, “Noninvasive coronary angiography--an acceptable alternative?” The New England Journal of Medicine, vol. 345, no. 26, pp. 1909-1910, 2001.

GBD 2013 Mortality and Causes of Death Collaborators, “Global, regional, and national age-sex specific all-cause and cause-specific mortality for 240 causes of death, 1990-2013: a systematic analysis for the Global Burden of Disease Study 2013,” The Lancet, vol. 385, no. 9963, pp. 117–171, 2015.

S. W. Waldo, E. A. Secemsky, C. O’Brien et al., “Surgical ineligibility and mortality among patients with unprotected left main or multivessel coronary artery disease undergoing percutaneous coronary intervention,” Circulation, vol. 130, no. 25, pp. 2295–2301, 2014.

R. Kalluri, “The biology and function of exosomes in cancer,” The Journal of Clinical Investigation, vol. 126, no. 4, pp. 1208–1215, 2016.

M. Mittelbrunn and F. Sanchez-Madrid, “Intercellular communication: diverse structures for exchange of genetic information,” Nature Reviews Molecular Cell Biology, vol. 13, no. 5, pp. 328–335, 2012.

G. Raposo and W. Stoorvogel, “Extracellular vesicles: exosomes, microvesicles, and friends,” The Journal of Cell Biology, vol. 200, no. 4, pp. 373–383, 2013.

C. Thery, M. Ostrowski, and E. Segura, “Membrane vesicles as conveyors of immune responses,” Nature Reviews Immunology, vol. 9, no. 8, pp. 581–593, 2009.

M. P. Caby, D. Lankar, C. Vincendeau-Scherrer, G. Raposo, and C. Bonnerot, “Exosomal-like vesicles are present in human blood plasma,” International Immunology, vol. 17, no. 7, pp. 879–887, 2005.

E. J. Hoorn, T. Pisitkun, R. Zietse et al., “Prospects for urinary proteomics: exosomes as a source of urinary biomarkers,” Nephrology, vol. 10, no. 3, pp. 283–290, 2005.

T. Pisitkun, R. F. Shen, and M. A. Knepper, “Identification and proteomic profiling of exosomes in human urine,” Proceedings of the National Academy of Sciences of the United States of America, vol. 101, no. 36, pp. 13368–13373, 2004.

K. Boriachek, M. N. Islam, A. Möller et al., “Biological functions and current advances in isolation and detection strategies for exosome nanovesicles,” Small, vol. 14, no. 6, p. 1702153, 2018.

A. Montecalvo, A. T. Larregina, W. J. Shufesky et al., “Mechanism of transfer of functional microRNAs between mouse dendritic cells via exosomes,” Blood, vol. 119, no. 3, pp. 756–766, 2012.

M. Mittelbrunn, C. Gutiérrez-Vázquez, C. Villarroya-Beltri et al., “Unidirectional transfer of microRNA-loaded exosomes from T cells to antigen- presenting cells,” Nature Communications, vol. 2, no. 1, 2011.

J. Skog, T. Würdinger, S. van Rijn et al., “Glioblastoma microvesicles transport RNA and proteins that promote tumour growth and provide diagnostic biomarkers,” Nature Cell Biology, vol. 10, no. 12, pp. 1470–1476, 2008.

J. A. Tickner, A. J. Urquhart, S.-A. Stephenson, D. J. Richard, and K. J. Oâ€™Byrne, “Functions and therapeutic roles of exosomes in cancer,” Frontiers in Oncology, vol. 4, 2014.

S. Memczak, M. Jens, A. Elefsinioti et al., “Circular RNAs are a large class of animal RNAs with regulatory potency,” Nature, vol. 495, no. 7441, pp. 333–338, 2013.

H. Suzuki and T. Tsukahara, “A view of pre-mRNA splicing from RNase R resistant RNAs,” International Journal of Molecular Sciences, vol. 15, no. 6, pp. 9331–9342, 2014.

L. L. Chen and L. Yang, “Regulation of circRNA biogenesis,” RNA Biology, vol. 12, no. 4, pp. 381–388, 2015.

Z. Li, C. Huang, C. Bao et al., “Exon-intron circular RNAs regulate transcription in the nucleus,” Nature Structural & Molecular Biology, vol. 22, no. 3, pp. 256–264, 2015.

Y. Zhang, X. O. Zhang, T. Chen et al., “Circular intronic long noncoding RNAs,” Molecular Cell, vol. 51, no. 6, pp. 792–806, 2013.

R. Ashwal-Fluss, M. Meyer, N. R. Pamudurti et al., “circRNA biogenesis competes with pre-mRNA splicing,” Molecular Cell, vol. 56, no. 1, pp. 55–66, 2014.

T. B. Hansen, T. I. Jensen, B. H. Clausen et al., “Natural RNA circles function as efficient microRNA sponges,” Nature, vol. 495, no. 7441, pp. 384–388, 2013.

D. Barbagallo, A. Caponnetto, D. Brex et al., “CircSMARCA5 regulates VEGFA mRNA splicing and angiogenesis in glioblastoma multiforme through the binding of SRSF1,” Cancers, vol. 11, no. 2, p. 194, 2019.

J. Salzman, R. E. Chen, M. N. Olsen, P. L. Wang, and P. O. Brown, “Cell-type specific features of circular RNA expression,” PLoS Genetics, vol. 9, no. 9, article e1003777, 2013.

W. R. Jeck, J. A. Sorrentino, K. Wang et al., “Circular RNAs are abundant, conserved, and associated with ALU repeats,” RNA, vol. 19, no. 2, pp. 141–157, 2013.

C. Preußer, L.-H. Hung, T. Schneider et al., “Selective release of circRNAs in platelet-derived extracellular vesicles,” Journal of Extracellular Vesicles, vol. 7, no. 1, p. 1424473, 2018.

Y. Li, Q. Zheng, C. Bao et al., “Circular RNA is enriched and stable in exosomes: a promising biomarker for cancer diagnosis,” Cell Research, vol. 25, no. 8, pp. 981–984, 2015.

X. Dai, C. Chen, Q. Yang et al., “Exosomal circRNA_100284 from arsenite-transformed cells, via microRNA-217 regulation of EZH2, is involved in the malignant transformation of human hepatic cells by accelerating the cell cycle and promoting cell proliferation,” Cell Death & Disease, vol. 9, no. 5, p. 454, 2018.

H. Zhang, T. Deng, S. Ge et al., “Exosome circRNA secreted from adipocytes promotes the growth of hepatocellular carcinoma by targeting deubiquitination-related USP7,” Oncogene, vol. 38, no. 15, pp. 2844–2859, 2019.

X. Li, Z. Zhao, D. Jian, W. Li, H. Tang, and M. Li, “Hsa-circRNA11783-2 in peripheral blood is correlated with coronary artery disease and type 2 diabetes mellitus,” Diabetes & Vascular Disease Research, vol. 14, no. 6, pp. 510–515, 2017.

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.

Z. Zhao, X. Li, C. Gao et al., “Peripheral blood circular RNA hsa_circ_0124644 can be used as a diagnostic biomarker of coronary artery disease,” Scientific Reports, vol. 7, no. 1, 2017.

Q. Zhang, I. Kandic, and M. J. Kutryk, “Dysregulation of angiogenesis-related microRNAs in endothelial progenitor cells from patients with coronary artery disease,” Biochemical and Biophysical Research Communications, vol. 405, no. 1, pp. 42–46, 2011.

H. Gao, R. R. Guddeti, Y. Matsuzawa et al., “Plasma levels of microRNA-145 are associated with severity of coronary artery disease,” PLoS One, vol. 10, no. 5, article e0123477, 2015.

Y. Minami, M. Satoh, C. Maesawa et al., “Effect of atorvastatin on microRNA 221/222 expression in endothelial progenitor cells obtained from patients with coronary artery disease,” European Journal of Clinical Investigation, vol. 39, no. 5, pp. 359–367, 2009.

Z. Zhao, X. Li, C. Gao et al., “Circulating miR-155, miR-145 and let-7c as diagnostic biomarkers of the coronary artery disease,” Scientific Reports, vol. 7, no. 1, 2017.

J. L. Sun, R. Han, J. H. Guo, X. Y. Li, X. L. Ma, and C. Y. Wang, “The diagnostic value of treadmill exercise test parameters for coronary artery disease,” Cell Biochemistry and Biophysics, vol. 65, no. 1, pp. 69–76, 2013.

A. Arbab-Zadeh, M. F. D. Carli, R. Cerci et al., “Accuracy of computed tomographic angiography and single-photon emission computed tomography-acquired myocardial perfusion imaging for the diagnosis of coronary artery disease,” Circulation: Cardiovascular Imaging, vol. 8, no. 10, p. e003533, 2015.

Y. Jiang, J. P. Tian, H. Wang, B. X. Chen, and F. H. Du, “Diagnostic value of combined parameters derived from ambulatory electrocardiography for detecting coronary artery disease in non-active chest pain patients,” Pakistan Journal of Medical Sciences, vol. 30, no. 6, pp. 1331–1335, 2014.

A. Firouzi, B. Ghadrdoost, and M. Haghjoo, “Accuracy of cardiogoniometry compared with electrocardiography in the diagnosis of coronary artery disease,” Research in Cardiovascular Medicine, vol. 4, no. 1, p. 1, 2015.