Overexpressed VDAC1 in breast cancer as a novel prognostic biomarker and correlates with immune infiltrates

Background: More and more evidence suggests that cancer is a mitochondrial metabolic disease recently and mitochondria dysfunction is critical to tumorigenesis. As a gatekeeper of mitochondria, the voltage-dependent anion channel 1 (VDAC1) is associated with the development of breast cancer (BC). However, its potential mechanism and clinical significance remain unclear; thus, in this research, we aimed to explore it. Methods: VDAC1 expression in BC tissues and normal tissues was obtained from The Cancer Genome Atlas (TCGA) and validated by datasets from the gene expression omnibus (GEO) database. Then, the relationships between VDAC1 expression and clinicopathological features were analyzed. Receiver operating characteristics (ROC) curves were used to identify the diagnostic value of VDAC1. The prognostic value was evaluated by Kaplan-Meier survival curves and Cox regression analysis. VDAC1 with its co-expression genes were subjected to enrichment analysis to explore potential mechanisms in BC and the protein-protein interaction (PPI) network was constructed. At last, the association between VDAC1 expression and infiltration levels of immune cell infiltration by various methods, as well as their corresponding markers, was analyzed. We also analyzed the correction between VDAC1 expression and eight immune checkpoint genes and the tumor immune dysfunction and exclusion (TIDE) scores of each BC sample in TCGA were calculated and the differences between high and low VDAC1 expression groups were analyzed. Results: VDAC1 expression was remarkably elevated in BC (p < 0.001), and high expression of VDAC1 was associated with the positive expression of ER (p = 0.004), PR (p = 0.033), and HER2 (p = 0.001). ROC analysis suggested that VDAC1 had diagnosed value in BC. The Kaplan-Meier analysis suggested that higher expression of VDAC1 was associated with shorter overall survival (OS), and further Cox regression analysis revealed that VDAC1 was an independent factor of unfavorable prognosis in BC patients. Enrichment analysis ofVDAC1 and its co-expression suggested that VDAC1 was related to the regulation of mitochondrial energy metabolism and protein modification, and the HIF-1 singing pathway might be the potential mechanism in BC. Notably, we found that VDAC1 expression was infiltration levels of most types of immune cells, as well as the expression of marker genes of immune cells. The ICGs PDCD1, CTLA4, LAG3,SIGLEC15, and TIGIT were negatively corrected with VDAC1 expression in BC. TIDE scores between the low and high expression groups showed no difference. Conclusion: Overexpressed VDAC1 in BC could be severed as a novel biomarker for diagnosis and VDAC1 was an independent factor for adverse prognosis prediction. Our study revealed that VDAC1 might inhibit tumor immunity and might be a novel therapeutic target in BC.