Investigating the role of FADS family members in breast cancer based on bioinformatic analysis and experimental validation

Abstract

Breast cancer (BC) is the most common malignant tumor in women worldwide. Emerging evidence indicates the significance of fatty acid metabolism in BC. Fatty acid desaturase (FADS) is closely associated with cancer occurrence and development. Here, bioinformatic analysis and experimental validation were applied to investigate the potential functions of FADS in BC. Several public databases, including TCGA, GEO, HPA, Kaplan–Meier plotter, STRING, DAVID, cBioPortal, TIMER, TRRUST, and LinkedOmics were used to determine mRNA/protein expression levels, prognostic significance, functional enrichment, genetic alterations, association with tumor-infiltrating immune cells, and related transcription factors and kinases. BC tissues showed higher and lower mRNA expression of FADS2/6/8 and FADS3/4/5, respectively. FADS1/2/6 and FADS3/4/5 showed higher and lower protein expression levels, respectively, in BC tissues. Moreover, FADS1/7 up- and FADS3/8 down-regulation predicted poor overall and recurrence-free survival, while FADS2/5 up- and FADS4 down-regulation were associated with poor recurrence-free survival. Receiver operating characteristic curves revealed that FADS2/3/4/8 were indicative diagnostic markers. FADS family members showing differential expression levels were associated with various clinical subtypes, clinical stages, lymph node metastasis status, copy number variants, DNA methylation, and miRNA regulation in BC. The mRNA expression level of FADS1/2/3/4/5/7/8 was observed to be significantly negatively correlated with DNA methylation. FADS1/2 upregulation was significantly correlated with clinical stages. FADS1/4 expression was obviously lower in BC patients with higher lymph node metastasis than lower lymph node metastasis, while FADS7/8 expression was obviously higher in BC patients with higher lymph node metastasis than lower lymph node metastasis. FADS family members showed varying degrees of genetic alterations, and Gene Ontology and KEGG pathway enrichment analyses suggested their involvement in lipid metabolism. Their expression level was correlated with immune cell infiltration levels. FADS2 was chosen for further validation analyses. We found FADS2 to be significantly over-expressed in clinical BC tissue samples. The proliferation, migration, and invasion abilities of MDA-MB-231 and BT474 cells were significantly reduced after FADS2 knockdown. Furthermore, FADS2 may promote the occurrence and development of BC cells via regulating the epithelial–mesenchymal transition (EMT) pathway. Altogether, our results suggest that FADS1/2/3/4 can serve as potential therapeutic targets, prognostic indicators, and diagnostic markers in patients with BC.