Comprehensive Analysis of mRNA Expression Profiling and Identification of Potential Diagnostic Biomarkers in Coronary Artery Disease

[Image: see text] The aim of this study is to investigate mRNA expression profiling by RNA sequencing (RNA-seq) in patients with coronary artery disease (CAD) and validate differentially expressed genes (DEGs) as novel biomarkers for CAD. Transcriptome-wide mRNA expression analysis of peripheral blood mononuclear cells was performed in five CAD patients and five controls. Functional enrichment analyses, protein–protein interaction network construction, and hub gene selection were further conducted. Relative expression levels of hub genes were validated by quantitative reverse transcription PCR in larger cohorts. Spearman correlation test and multiple linear regression analysis were applied to examine the relationship between confounding factors with severity of coronary artery atherosclerosis. Receiver operating characteristic (ROC) curve analysis was adopted to identify potentially diagnostic biomarkers for CAD. A total of 527 upregulated and 653 downregulated mRNAs were identified as DEGs in CAD patients. The relative expression levels of beta-transducin repeat containing E3 ubiquitin protein ligase (BTRC), F-box and leucine-rich repeat protein 4 (FBXL4), ubiquitin conjugating enzyme E2 D2 (UBE2D2), and ankyrin repeat and SOCS box containing 1 (ASB1) were significantly different between two groups (all p ≤ 0.05). The severity of coronary artery atherosclerosis was negatively associated with the BTRC gene relative expression level (r = −0.323, p < 0.001) and positively with UBE2D2 (r = 0.285, p < 0.001). ROC analysis of BTRC and UBE2D2 genes showed that the areas under the curve were 0.782 (95% CI: 0.720–0.845, p < 0.001) and 0.753 (95% CI: 0.681–0.824, p < 0.001), respectively. We described the characteristics of mRNA expression in the peripheral blood of CAD patients and controls by RNA-seq. Combined with Spearman correlation analysis and ROC analyses, BTRC and UBE2D2 genes had significantly diagnostic values, which may have potential to act as novel diagnostic biomarkers and therapeutic targets for CAD.