Abstract
BACKGROUND: Coronary artery disease (CAD) is the leading cause of death worldwide, and aberrant phagocytosis may be involved in its development. Understanding this aspect may provide new avenues for prompt CAD diagnosis. METHODS: CAD-related information was obtained from Gene Expression Omnibus datasets GSE66360, GSE113079, and GSE59421. We identified 995 upregulated and 1086 downregulated differentially expressed genes (DEGs) in GSE66360. Weighted gene co-expression network analysis revealed a module of 503 genes relevant to CAD. Using clusterProfiler, we revealed 32 CAD-related PRFs. Eight candidate genes were identified in a protein-protein interaction network. Machine learning algorithms identified CAD biomarkers that underwent gene set enrichment analysis, immune cell analysis with CIBERSORT, microRNA (miRNA) prediction using the miRWalk database, transcription factor (TF) level predication through ChEA3, and drug prediction with DGIdb. Cytoscape visualized the miRNA -mRNA- TF, miRNA-single nucleotide polymorphism-mRNA, and biomarker-drug networks. RESULTS: IL1B, TLR2, FCGR2A, SYK, FCER1G, and HCK were identified as CAD biomarkers. The area under the curve of a diagnostic model based on the six biomarkers was > 0.7 for the GSE66360 and GSE113079 datasets. Gene set enrichment analysis revealed differences in their biological pathways. CIBERSORT revealed that 10 immune cell types were differentially expressed between the CAD and control groups. The TF-mRNA-miRNA network showed that has-miR-1207-5p regulates HCK and FCER1G expression and that RUNX1 and SPI may be important TFs. Ninety-five drugs were predicted, including aspirin, which influenced ILIB and FCERIG. CONCLUSION: In this study, six biomarkers (IL1B, TLR2, FCGR2A, SYK, FCER1G, and HCK) related to CAD phagocytic regulatory factors were identified, and their expression regulatory relationships in CAD were further studied, providing a deeper understanding of the pathogenesis, diagnosis, and potential treatment strategies of CAD.