Abstract
BACKGROUND: Immune cell infiltration (ICI) has a close relationship with the progression of atherosclerosis (AS). Therefore, the current study was aimed to explore the role of genes related to ICI and to investigate potential mechanisms in AS. METHODS: Single-sample gene set enrichment analysis (ssGSEA) was applied to explore immune infiltration in AS and controls. Genes related to immune infitration were mined by weighted gene co-expression network analysis (WGCNA). The function of those genes were analyzed by enrichment analyses of the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO). The interactions among those genes were visualized in the protein-protein interaction (PPI) network, followed by identification of hub genes through Cytoscape software. A receiver operating characteristic (ROC) plot was generated to assess the performance of hub genes in AS diagnosis. The expressions of hub genes were measured by reverse transcription quantitative real-time PCR (RT-qPCR) in human leukemia monocyticcell line (THP-1) derived foam cells and macrophages, which mimic AS and control, respectively. RESULTS: We observed that the proportions of 27 immune cells were significantly elevated in AS. Subsequent integrative analyses of differential expression and WGCNA identified 99 immune cell-related differentially expressed genes (DEGs) between AS and control. Those DEGs were associated with tryptophan metabolism and extracellular matrix (ECM)-related functions. Moreover, by constructing the PPI network, we found 11 hub immune cell-related genes in AS. The expression pattern and receiver ROC analyses in two independent datasets showed that calsequestrin 2 (CASQ2), nexilin F-Actin binding protein (NEXN), matrix metallopeptidase 12 (MMP12), C-X-C motif chemokine ligand 10 (CXCL10), phospholamban (PLN), heme oxygenase 1 (HMOX1), ryanodine receptor 2 (RYR2), chitinase 3 like 1 (CHI3L1), matrix metallopeptidase 9 (MMP9), actin alpha cardiac muscle 1 (ACTC1) had good performance in distinguishing AS from control samples. Furthermore, those biomarkers were shown to be correlated with angiogenesis and immune checkpoints. In addition, we found 239 miRNAs and 47 transcription factor s (TFs), which may target those biomarkers and regulate their expressions. Finally, we found that RT-qPCR results were consistent with sequencing results.