Therefore, in this study, weighted gene co-expression network analysis (WGCNA), single-sample gene set enrichment analysis (ssGSEA), multivariate COX analysis, least absolute shrinkage, and selection operator (LASSO), and support vector machine-recursive feature elimination (SVM-RFE) algorithm, along with a series of analyses were performed, and four immune-related genes (APOD, CXCL14, IL33, and LIFR) were identified as biomarkers correlated with breast cancer prognosis. The gene discussed is IL33; the disease is breast carcinoma.