Construction of a fatty acid metabolism-related gene signature for predicting prognosis and immune response in breast cancer

Abstract

Background: Breast cancer has the highest incidence among malignant tumors in women, and its prevalence ranks first in global cancer morbidity.Aim: This study aimed to explore the feasibility of a prognostic model for patients with breast cancer based on the differential expression of genes related to fatty acid metabolism.Methods: The mRNA expression matrix of breast cancer and paracancer tissues was downloaded from The Cancer Genome Atlas database. The differentially expressed genes related to fatty acid metabolism were screened in R language. The TRRUST database was used to predict transcriptional regulators related to hub genes and construct an mRNA–transcription factor interaction network. A consensus clustering approach was used to identify different fatty acid regulatory patterns. In combination with patient survival data, Lasso and multivariate Cox proportional risk regression models were used to establish polygenic prognostic models based on fatty acid metabolism. The median risk score was used to categorize patients into high- and low-risk groups. Kaplan–Meier survival curves were used to analyze the survival differences between both groups. The Cox regression analysis included risk score and clinicopathological factors to determine whether risk score was an independent risk factor. Models based on genes associated with fatty acid metabolism were evaluated using receiver operating characteristic curves. A comparison was made between risk score levels and the fatty acid metabolism-associated genes in different subtypes of breast cancer. The differential gene sets of the Kyoto Encyclopedia of Genes and Genomes for screening high- and low-risk populations were compared using a gene set enrichment analysis. Furthermore, we utilized CIBERSORT to examine the abundance of immune cells in breast cancer in different clustering models.Results: High expression levels of ALDH1A1 and UBE2L6 prevented breast cancer, whereas high RDH16 expression levels increased its risk. Our comprehensive assessment of the association between prognostic risk scoring models and tumor microenvironment characteristics showed significant differences in the abundance of various immune cells between high- and low-risk breast cancer patients.Conclusions: By assessing fatty acid metabolism patterns, we gained a better understanding of the infiltration characteristics of the tumor microenvironment. Our findings are valuable for prognosis prediction and treatment of patients with breast cancer based on their clinicopathological characteristics.