Integrating Single-Cell RNA-Seq and Bulk RNA-Seq Data to Explore the Key Role of Fatty Acid Metabolism in Breast Cancer

Cancer immune escape is associated with the metabolic reprogramming of the various infiltrating cells in the tumor microenvironment (TME), and combining metabolic targets with immunotherapy shows great promise for improving clinical outcomes. Among all metabolic processes, lipid metabolism, especially fatty acid metabolism (FAM), plays a major role in cancer cell survival, migration, and proliferation. However, the mechanisms and functions of FAM in the tumor immune microenvironment remain poorly understood. We screened 309 fatty acid metabolism-related genes (FMGs) for differential expression, identifying 121 differentially expressed genes. Univariate Cox regression models in The Cancer Genome Atlas (TCGA) database were then utilized to identify the 15 FMGs associated with overall survival. We systematically evaluated the correlation between FMGs’ modification patterns and the TME, prognosis, and immunotherapy. The FMGsScore was constructed to quantify the FMG modification patterns using principal component analysis. Three clusters based on FMGs were demonstrated in breast cancer, with three patterns of distinct immune cell infiltration and biological behavior. An FMGsScore signature was constructed to reveal that patients with a low FMGsScore had higher immune checkpoint expression, higher immune checkpoint inhibitor (ICI) scores, increased immune microenvironment infiltration, better survival advantage, and were more sensitive to immunotherapy than those with a high FMGsScore. Finally, the expression and function of the signature key gene NDUFAB1 were examined by in vitro experiments. This study significantly demonstrates the substantial impact of FMGs on the immune microenvironment of breast cancer, and that FMGsScores can be used to guide the prediction of immunotherapy efficacy in breast cancer patients. In vitro experiments, knockdown of the NDUFAB1 gene resulted in reduced proliferation and migration of MCF-7 and MDA-MB-231 cell lines.