Abstract
Circular RNAs (circRNAs) are involved in many physiological functions. Whether circulating circRNAs serve as markers for coronary artery disease (CAD) is unknown. Seven CAD-related microarray datasets were downloaded from the Gene Expression Omnibus (GEO) database and were analyzed using clustering and functional enrichment to identify hub mRNAs and miRNAs. StarBase V3.0 and circinteractome databases were used to predict interactions between circRNAs and miRNAs whereas miRwalk and DIANA TOOLS were used to predict interactions between miRNAs and mRNAs. Altogether, this helped establish a circRNA-miRNA-mRNA triple network for diagnosis of CAD. Five non-coding RNAs (ncRNAs) were identified in our study population with the use of quantitative real-time PCR (RT-PCR). The prognostic values of circYOD1, hsa-miR-21-3p and hsa-miR-296-3p were evaluated using a receiver operating characteristic (ROC) curve. A CAD circRNA-miRNA-mRNA network was established from our analyses containing one circRNA, four miRNAs and thirteen mRNAs. After performing RT-PCR validation between CAD and non-CAD samples, only three ncRNAs of five ncRNAs showed significance for further analysis. The area under ROC curve (AUC) of circ-YOD1 was 0.824, the AUC of hsa-miR-21-3p was 0.731 and hsa-miR-296-3p was 0.776. The pairwise comparison results showed that circ-YOD1 had statistical significance (PYOD1-21 < 0.01 and PYOD1-296 < 0.05). The results of functional enrichment analysis of interacting genes and microRNAs showed that the shared circ-YOD1 may act as a new biomarker for CAD. Our investigation of the triple regulatory networks of circRNA-miRNA-mRNA in CAD revealed circ-YOD1 as a potential biomarker for CAD.