The underlying molecular mechanisms and biomarkers of plaque vulnerability based on bioinformatics analysis

The study aimed to identify the underlying differentially expressed genes (DEGs) and mechanism of unstable atherosclerotic plaque using bioinformatics methods.

GSE120521, which includes four unstable samples and four stable atherosclerotic samples, was downloaded from the GEO database. DEGs were identified using LIMMA. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of DEGs were performed using the Database for metascape Visualization online tool. Based on the STRING database, protein-protein interactions (PPIs) network among DEGs were constructed. Regulatory networks were visualized using Cytoscape. We use the xCell to analyze the different immune cell subtypes.

A total of 1626 DEGs (1034 up-regulated and 592 down-regulated DEGs) were identified between unstable and stable samples. I pulled 62 transcription factors (34 up-regulated TFs and 28 down-regulated TFs) from the Trust database. The up-regulated TFs were mainly enrichment in positive regulation of myeloid leukocyte differentiation, and the down-regulated TFs were mainly enrichment in connective tissue development. In the PPI network, RB1, CEBPA, PPARG, BATF was the most significantly up-regulated gene in ruptured atherosclerotic samples. The immune cell composition enriched in CD cells and macrophages in the unstable carotid plaque.

Upregulated RB1, CEBPA, PPARG, BATF and down-regulated SRF, MYOCD, HEY2, GATA6 might perform critical promotional roles in atherosclerotic plaque rupture, furthermore, number and polarization of macrophages may play an important role in vulnerable plaques.