Identification of potential M2 macrophage-associated diagnostic biomarkers in coronary artery disease

M2 macrophages have been reported to be important in the progression of coronary artery disease (CAD). Thus, the present study aims at exploring the diagnostic value of M2 macrophage-associated genes in CAD.

Transcriptome profile of CAD and control samples were downloaded from Gene Expression Omnibus database. The proportion of immune cells was analyzed using cell type identification by estimating relative subsets of RNA transcripts. Weighted Gene Co-expression Network Analysis (WGCNA) was carried out to screen the relevant module associated with M2 macrophages. Differential CAD and control samples of expressed genes (DEGs) were identified by the limma R package. Functional enrichment analysis by means of the clusterProfiler R package. Least absolute shrinkage and selection operator (LASSO) and random forest (RF) algorithms were carried out to select signature genes. Receiver operating curves (ROC) were plotted to evaluate the diagnostic value of selected signature genes. The expressions of potential diagnostic markers were validated by RT-qPCR. The ceRNA network of diagnostic biomarkers was constructed via miRwalk and Starbase database. CMap database was used to screen candidate drugs in the treatment of CAD by targeting diagnostic biomarkers.

A total of 166 M2 macrophage-associated genes were identified by WGCNA. By intersecting those genes with 879 DEGs, 53 M2 macrophage-associated DEGs were obtained in the present study. By LASSO, RF, and ROC analyses, C1orf105, CCL22, CRYGB, FRK, GAP43, REG1P, CALB1, and PTPN21 were identified as potential diagnostic biomarkers. RT-qPCR showed the consistent expression patterns of diagnostic biomarkers between GEO dataset and clinical samples. Perhexiline, alimemazine and mecamylamine were found to be potential drugs in the treatment of CAD.

We identified eight M2 macrophage-associated diagnostic biomarkers and candidate drugs for the CAD treatment.