- Open Access
Overexpressed VDAC1 in breast cancer as a novel prognostic biomarker and correlates with immune infiltrates
World Journal of Surgical Oncology volume 20, Article number: 211 (2022)
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.
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.
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 of VDAC1 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.
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.
Breast cancer (BC) is a public health problem that is plaguing women around the world over the past few decades. In 2020, the global incidence of BC has surpassed lung cancer as the most common malignant tumor in women with approximately 2.26 million new cases . With the advancement of diagnostic and therapeutic strategies, the 5-year survival rate in female patients diagnosed with BC has improved, but the mortality rate of BC is still highest in female cancers [2, 3]. Therefore, early diagnosis and prognosis prediction are crucial for physicians to monitor the disease progression and develop an individualized treatment plan. Currently, early detection of BC mainly relies on some serum biomarkers such as carcinoembryonic antigen (CEA) and carbohydrate antigen 153 (CA153), but their sensitivity and specificity are still low [4, 5]. In addition, there are no effective biomarkers for prognosis prediction of BC currently. Therefore, it is an urgent need for us to identify some novel biomarkers with high sensitivity and specificity for early diagnosis and effective biomarkers for prognosis prediction that is important for optimal treatment planning.
It is well known that mitochondria exist in almost all eukaryotic cells as energy-converting organelles that provide energy to cells through oxidative phosphorylation. Mitochondria are essential to the regulation of cellular energy metabolism, biosynthesis, and cell death, and the dysfunction of which is closely related to the incidence and development of many diseases, including cancers . The voltage-dependent anion channel (VDAC) protein is the most abundant pore-forming protein located in the outer mitochondrial membrane (OMM) of eukaryotic cells working as a gatekeeper of mitochondria that regulate the entry and exit of metabolites, Ca2+, fatty acid ions, and reactive oxygen species across the OMM, as well as a hub protein that interacts with other proteins from cytosol and endoplasmic reticulum to regulate cellular metabolism and apoptosis [7,8,9]. The VDAC protein has been identified as three isoforms encoded by three homologous genes: VDAC1, VDAC2, and VDAC3. VDAC1 has the most abundant expression level among them, VDAC2 is known as an anti-apoptotic protein, and VDAC3 is involved in the ciliary disassembly . It has been reported that the overexpression of VDAC1 is related to many diseases including neurodegeneration, cardiovascular diseases, type 2 diabetes, and different types of cancers. VDAC1 might be served as a novel pharmacological target for anti-cancer therapeutics [8, 9, 11].
However, the evidence about the correlation between the expression level of VDAC1 and BC has been rarely reported and remains unclear, which deserves further exploration. Therefore, the purpose of this study was to analyze the relationships between VDAC1 expression and clinicopathological features, diagnosis value, and prognosis value in BC patients utilizing various online databases. In addition, the co-expression genes of VDAC1 and the relationships between VDAC1 expression and tumor-infiltrating immune cells, as well as their corresponding gene markers, would be analyzed, which would explore the potential mechanism of VDAC1’s role involved in the incidence and development of BC.
Materials and methods
Data collection and gene expression analysis
The expression levels of VDAC1 in different types of cancers were analyzed by the Tumor Immune Estimation Resource (TIMER) database (https://cistrome.shinyapps.io/timer/)  with the data from The Cancer Genome Atlas (TCGA). The data of the VDAC1 expression in BC patients were downloaded from TCGA (https://portal.gdc.cancer.gov/), including 1065 tumor samples and 111 normal samples, and then the data was converted to as log2 Transcripts Per Million (TPM). Three datasets (GSE21422 , GSE33447 , and GSE31192 ) were selected as validation sets from the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/). The expression levels of VDAC1 between normal and tumor groups were analyzed by T-test. All of the p values in our study less than 0.05 was considered statistically significant. In addition, the figures regarding the immunohistochemistry of VDAC1 expression in human BC tissues and normal tissues were obtained from the Human Protein Atlas (https://www.proteinatlas.org/).
Analysis of the clinicopathological features and diagnostic value
Data of the clinicopathological features and expression levels of VDAC1 in BC patients were downloaded from TCGA and divided into the low-expression group and high-expression group according to the median expression value of VDAC1. The clinicopathological features are depicted in Table 1. T-test or Kruskal-Wallis test was performed to analyze the different expression levels of VDAC1 among different groups. Moreover, the receiver operating characteristics (ROC) curves were then performed to identify the diagnosis value of VDAC1 to distinguish the normal group and the tumor group, as well as the groups with different types of BC.
Overall survival analysis
As noted above, the samples were divided into two groups according to the expression level of VDAC1. The overall survival (OS) curves were visualized by Kaplan-Meier analysis and log-rank test. Then, we further downloaded two datasets (GSE1456  and GSE159956) from the GEO database to validate the OS curves. In addition, univariate and multivariate Cox regression analyses were used to calculate death hazard ratios of clinicopathological features and VDAC1 expression and identify whether VDAC1 could be served as an independent prognostic factor for BC.
Identification and enrichment analysis of co-expression genes
The co-expression genes of VDAC1 in the data of BC from TCGA were identified via the database LinkedOmics (http://www.linkedomics.org/login.php)  and Pearson correlation test. Then, we selected the top 200 co-expression genes to perform Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis by the online database Metascape (http://metascape.org)  to explore the VDAC1-related molecular mechanisms. We set the minimum counts larger than 3, p-value less than 0.05, and minimum enrichment factors larger than 1.5 as thresholds. Moreover, the top 200 co-expression genes were uploaded to the STRING database (https://cn.string-db.org/)  for the construction of the protein-protein interaction (PPI) network with a minimum required interaction score of 0.9 and then visualized by the gene-networking tool Cytoscape (version 3.8.2).
Immune cell analysis
The single sample GSEA method from the R package “GSVA”  was applied to present infiltration enrichment of 24 common immune cells and the relation between VDAC1 expression with immune cell infiltration was analyzed by the Spearman test. The Microenvironment Cell Populations-counter (MCP-counter) algorithm  and the TIMER algorithm were also employed to analyze the correction between VDAC1 expression and immune cell infiltration. The Estimation of Stromal and Immune cells in Malignant Tumor tissues using Expression (ESTIMATE) algorithm  was used to calculate and compare stromal scores, immune scores, and ESTIMATE scores between low and high VDAC1 expression groups of BC samples from TCGA dataset by Wilcoxon rank-sum test. Moreover, we further analyzed the association of the expression of VDAC1 in BC with multiple marker genes of immune cells by the TIMER database.
Immune checkpoint gene analysis and immunity therapy response prediction
Patients with high expression of immune checkpoint inhibitors (ICIs) will receive greater benefits from ICI therapy . We analyzed the correction between VDAC1 expression and eight immune checkpoint genes (ICGs) by the Spearman test, which was SIGLEC15, TIGIT, CD274, HAVCR2, PDCD1, CTLA4, LAG3, and PDCD1LG2. 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 by the Wilcoxon rank-sum test to predict ICI therapy response.
VDAC1 exhibited higher expression levels in BC than that in normal tissues
Initially, we evaluated the expression levels of VDAC1 in different types of human malignant tumors from the TCGA RNA-seq data via the TIMER (Fig. 1A). The expression level of VDAC1 in the BC tissue group was remarkably higher than that in the normal tissue group. Moreover, the expression of VDAC1 was especially higher in human epidermal growth factor receptor 2 (HER2)-enriched BC than that in luminal and basal-like BC tissues. In addition, the result revealed that VDAC1 exhibited significantly higher expression levels in many types of malignant tumors such as bladder cancer (BLCA), cholangiocarcinoma (CHOL), chromophobe renal cell carcinoma (KICH), liver hepatocellular carcinoma (LIHC), lung adenocarcinoma (LUAD), and stomach adenocarcinomas (STAD). Analysis of the expression of VDAC1 expression between BC tissue group and normal tissue group, as well as matched normal tissues, suggested that VDAC1 exhibit remarkably higher expression level in BC (Fig. 1B, C) (p < 0.001, respectively). The result was further validated by three datasets from the GEO database (Fig. 1D–F) (p < 0.001; p = 0.042; p = 0.001, respectively). Correspondingly, the immunohistochemistry from the Human Protein Atlas revealed that VDAC1 protein was significantly expressed in BC tissues, especially in cytoplasmic and membranous of BC cells, and rarely expressed in normal tissues (Fig. 1G, H).
Expression of VDAC1 was associated with ER, PR, and HER2
The clinicopathological features of the BC patients from the TCGA dataset are shown in Table 1. We analyzed the expression of VDAC1 in BC patients with different clinicopathological features and the samples with indeterminate information were excluded. The result of our analysis shows that there was no significant difference in age (p = 0.681), T stages (p = 0.328), N stages (p = 0.374), M stages (p = 0.980), and pathologic stages (p = 0.282). The high expression of VDAC1 was correlated with the positive expression of estrogen receptor (ER) (p = 0.004), progesterone receptor (PR) (p = 0.033), and HER2 (p = 0.001) (Fig. 2).
VDAC1 could be served as a biomarker for BC diagnosis
The ROC curves were used to identify the effectiveness of VDAC1 mRNA expression level to distinguish BC tissues from the normal tissues. The area under the curve (AUC) was 0.854 [95% confidence interval (CI): 0.825–0.882], with 76.0% sensitivity and 82.9% specificity (Fig. 3A). The expression of VDAC1 was validated to be associated with the status of ER, PR, and HER2 in the above findings; thus, we further evaluate the diagnostic value of VDAC1 in negative and positive expression groups of ER, PR, and HER2, respectively. In the comparison of ER-negative expression and normal tissue groups, the AUC was 0.812 (95% CI: 0.768–0.857), with 70.5% sensitivity and 82.9% specificity (Fig. 3B). In the comparison of ER-positive expression and normal tissue groups, the AUC was 0.866 (95% CI: 0.837–0.894), with 76.3% sensitivity and 83.8% specificity (Fig. 3C). The AUC of VDAC1 was 0.837 (95% CI: 0.799–0.874 ) with 72.2% sensitivity and 83.8% specificity (Fig. 3D) in the comparison between PR-negative expression and normal tissue groups and the AUC was 0.863 in the comparison between PR-positive and normal tissue groups (95% CI: 0.833–0.892) with 77.2% sensitivity and 82.9% specificity (Fig. 3E). The AUC of VDAC1 was 0.849 (95% CI: 0.816–0.882) with 74.3% sensitivity and 82.9% specificity in the comparison between HER-negative expression and normal tissue groups (Fig. 3F). Intriguingly, we found that VDAC1 exhibits the best effectiveness to distinguish HER2-positive tumor tissues from the normal tissues, the AUC of which was 0.900 (95% CI: 0.863–0.938) with 74.5% sensitivity and 95.5% specificity (Fig. 3G).
VDAC1 was an independent factor for prognosis prediction in BC
To identify the correlation between VDAC1 expression and the OS of BC patients from TCGA, the survival curves were visualized by Kaplan-Meier analysis and log-rank test. As indicated in Fig. 4A, higher VDAC1 mRNA expression was remarkably related to shorter OS (HR = 1.76, p = 0.001). The result was validated by two datasets from GEO (Fig. 4B, C).
Cox regression analysis was further performed to evaluate the prognostic value of VDAC1 in BC. The p values of the clinicopathological features less than 0.1 in the univariate Cox model were considered statistically significant and the features were then further put into multivariate Cox analysis. As indicated in Fig. 4D, OS of BC patients was significantly correlated with age, T stages, N stages, pathologic stages, ER status, HER2 status, and VDAC1 expression. In the multivariate Cox model, the clinicopathological features including age, T stage, and ER expression, as well as VDAC1 expression, were independent predictors for poor prognosis prediction in BC (Fig. 4E). Overall, our results revealed that VDAC1 was an independent factor of unfavorable prognosis in BC patients.
Enrichment analysis and network establishment for co-expression genes of VDAC1 in BC
To explore the possible mechanism and biological function of VDAC1’s role involved in BC, we identified the co-expression genes of VDAC1 in BC via the online database LinkedOmics. As shown in Fig. 5A, a total of 10,651 genes were significantly corrected with the expression of VDAC1, including 3651 positively corrected genes (red dots) and 7000 negatively corrected genes (green dots), and the false discovery rate (FDR) was less than 0.01. The top 50 positively and negatively corrected genes were shown in two heatmaps (Fig. 5B, C), respectively.
The top 200 co-expression genes and VDAC1 were then selected for further GO and KEGG pathway enrichment analysis by the Matascape. The results of GO analysis suggested that VDAC1 and its related genes are remarkably enriched in protein folding, ATP metabolic process, mitochondrion organization, etc. in the biological process (BP) ontology; mitochondrial envelope, intracellular protein-containing complex, respirasome, etc. in the cellular component (CC) ontology; and unfolded protein binding, etc. in the molecular function (MF) ontology (Fig. 5D–G). KEGG pathway enrichment analysis showed that Parkinson disease, proteasome, and cell cycle were the most enriched pathway.
In addition, the top 200 co-expression genes and VDAC1 were uploaded to STRING for construction of the PPI network and then were visualized via the gene-networking tool Cytoscape (Fig. 5H). The PPI network contained 84 nodes and 207 edges, in which VDAC1 was the hub gene that related to another 4 genes.
VDAC1 was correlated with tumor-infiltrating immune cells in BC
Firstly, we used the ssGSEA method to present infiltration enrichment of 24 common types of immune cells in BC. Subsequently, the relation between VDAC1 expression with immune cell infiltration was analyzed by Spearman’ analysis. A total of 24 immune cells and their relationships with VDAC1 expression are shown in Fig. 6A and 17 types of immune cells were significantly correlated with VDAC1 expression. There were 15 types of immune cells negatively corrected with VDAC1 expression, which were plasmacytoid DCs (pDCs) (r = −0.377, p = 2.60e−37), CD8 + T cells (r = −0.232, p = 2.34e−14), cytotoxic cells (r = −0.215, p = 1.71e−12), natural killer (NK) cells (r = −0.206, p = 1.24e−11), T effector memory (Tem) cells (r = −0.203, p = 2.30e−11), B cells (r = −0.191, p = 3.53e−10), T cells (r = −0.182, p = 2.45e−09), dendritic cells (DCs) (r = −0.180, p = 3.33e−09), immature DCs (iDCs) (r = −0.180, p = 3.45e−09), NK 56-cells (r = −0.123, p = 5.84e−05), T follicular helper (Tfh) cells (r = −0.117, p = 1.23e−04), type 17 Th (Th17) cells (r = −0.117, p = 1.26e−04), neutrophils (r = −0.114, p = 2.08e−04), type 1 Th (Th1) cells (r = −0.086, p = 0.005), and macrophages (r = −0.071, p = 0.020). Two types of immune cells positively corrected with VDAC1 expression were type 2 Th (Th2) cells (r = 0.329, p = 0.000), and T helper (Th) cells (r = 0.077, p = 0.013) (Fig. 6B). The results of the MCP-counter method suggested that VDAC1 expression was negatively corrected with seven kinds of immune cells (Fig. 6C) and the TIMER algorithm suggested that VDAC1 expression was negatively corrected with infiltrating levels of six common immune cells (Fig. 6D). The ESTIMATE algorithm suggested that the high-VDAC1 group had lower immune scores (p < 0.001), stromal scores (p < 0.001), and ESTIMATE scores (p < 0.001) than the low-VDAC1 group (Fig. 6E).
Correction between the expression of VDAC1 and marker genes of immune cells in BC
We further analyzed the association between the expression of VDAC1 and various immune signatures of various immune cells via TIMER databases. A total of 60 immune marker genes were analyzed, which were well accepted as corresponding markers of different types of immune cells, such as B cells, macrophages, neutrophils, monocytes, NK cells, dendritic cells, and different functional T cells. Our results revealed that most immune markers of immune cells were remarkably associated with the expression of VDAC1 in BC (Table 2).
Correction between the expression of VDAC1 and ICGs in BC
At last, we analyzed the association between the expression of VDAC1 and eight ICGs. We found that the expression of VDAC1 was negatively corrected with the expression of PDCD1 (r = −0.160, p < 0.001), CTLA4 (r = −0.084, p = 0.006), LAG3 (r = −0.111, p < 0.001), SIGLEC15 (r = −0.062, p = 0.044), and TIGIT (r = −0.089, p = 0.003) (Fig. 7A). Further analysis showed that there was no difference in the TIDE score between the high-VDAC1 group and the low-VDAC1 group (p = 0.798) (Fig. 7B) despite the proportion of patients who responded to ICI therapy in the high-expression group (37.1%) was higher than in the low-expression group (33.2%).
Breast cancer is the most common malignant tumor in women. The incidence of young patients diagnosed with BC has been rising over the past decade . Many biomarkers have been proven to be involved in the BC progression and associated with the prognosis of BC patients [26,27,28,29]. The prognosis role of cancer driver genes and methylated genes in BC have also been reported [30, 31]. Interestingly, a high systemic immune-inflammation index has been reported to predict the poor prognosis of BC patients as a promising indicator . In addition, absolute lymphocyte count and insulin resistance have been reported to be served as predictors for chemotherapy response [33, 34]. Tang et al. have successfully developed a nomogram to predict pathological complete response after neoadjuvant chemotherapy of ER-positive and HER2-negative BC patients . All in all, in the future, more and more tools would be used as predictors in BC treatment.
In the present study, we explored the clinical significance of a mitochondrial porin, VDAC1, in BC via the datasets obtained from TCGA. We also found that increased expression of VDAC1 was related to the positive status of ER, PR, and especially HER2. We observed that the expression level of VDAC1 was remarkably increased in BC tissues compared to normal controls, which was consistent with a recent study that detected the expression of VDAC1 protein in BC tissues and benign breast lesions via immunohistochemistry  and validated by 3 datasets obtained from GEO database. In addition, we found that the expression of VDAC1 was higher in many types of tumor tissues than that in corresponding normal tissues, which was consistent with previous studies [37,38,39] and suggested that VDAC1 might act as an oncogene that is related to tumorigenesis and tumor progression. Currently, early detection is critical to the treatment of cancers. In our study, we performed ROC analysis and found that VDAC1 was an effective biomarker for the diagnosis of BC patients, especially HER2-positive BC patients. Moreover, in our study, we revealed that high expression of VDAC1 was closely related to a shorter OS rate of BC patients. Cox regression analysis indicated that VDAC1 could be severed as an independent risk factor for poor prognosis of BC patients, which is consistent with previous research . VDAC1 has also been identified with poor outcomes in other types of malignant tumors, such as pancreatic cancer, hepatocellular carcinoma, and cervical cancer [40,41,42]. These results might support that VDAC1 plays an important role in the development of cancers by different potential mechanisms.
We further attempted to explore the potential mechanisms and functions of VDAC1 involved in BC; thus, we performed GO and KEGG pathway enrichment analysis of co-expression of VDAC1. GO analysis suggested that VDAC1 and its co-expression showed significant enrichment in mitochondrial energy metabolism and protein modification. VDAC1 is a multi-functional channel protein and served as a hub that is involved in the control of cell metabolism, oxidative stress, apoptosis, mtDNA release, and more . Overexpression of VDAC1 has been validated to be associated with cancers [8, 9, 37,38,39,40], but its role involved in BC remains unclear. In recent years, more and more evidence suggests that cancer is a mitochondrial metabolic disease, and dysfunction of mitochondria is critical to tumorigenesis . As a gatekeeper of mitochondria, VDAC1 can decide the fate of cancer cells by regulating metabolic and energetic functions. The overexpressed VDAC1 in the cancer cells can interact with hexokinase (HK), a rate-limiting enzyme of glycolysis, and then promote mitochondrial ATP to coupling to glucose, which contributes to cancer cell metabolism [8, 9, 11]. VDAC1 has also been found to be involved in anti-apoptosis by interacting with anti-apoptotic proteins overexpressed in cancer such as Bcl-2, Bcl-xL, HK, and mediating the release of Cytc to prevent the cancer cells from apoptosis [8, 9, 11]. A recent study suggested that the high expression of VDAC1 in BC could be inhibited by the bromodomain inhibitor (JQ1) and associated with the resistance to JQ1 . It has been reported that mitochondrial fission factor (MFF) can bind to VDAC1 and the MFF-VDAC1 complex can be severed as an actionable therapeutic target in BC . Moreover, microRNA-7 was validated to decrease the expression of VDAC1 to inhibit hepatocellular carcinoma proliferation and metastasis . In addition, VDAC1 was found to be upregulated by microRNA-320a to promote proliferation and invasion of non-small cell lung cancer . These findings reveal that VDAC1 could be severed as a novel drug target for the treatment of BC. In our pathway enrichment analysis, the hypoxia-inducible factor 1 (HIF-1) signaling pathway was enriched by VDAC1 co-expression genes. HIF-1 signaling is known to play an important role in the tumor microenvironment; thus, HIF-1 signaling is a promising target for the treatment of cancers .
In our study, we also explored the underlying relationship between VDAC1 expression and immune cell infiltration. Our study suggested that VDAC1 expression was negatively corrected with the infiltration levels of most types of immune cells, such as DCs and CD8 + T cells, as well as the expression of their corresponding markers. It is well accepted now that CD8 + T cells are central in mediating anti-cancer immunity and activated by DCs to initiate anti-cancer immunity . VDAC1 expression was negatively corrected with neutrophils and macrophages, which play important roles in anti-cancer immunity . Further analysis of the gene markers of immune cells showed that markers of M2 macrophages such as CD163 and MS4A4A had weak correlations with BICC1 expression. M2 macrophage marker NOS2 and TAM marker CD80 showed moderate relationships with VDAC1 expression, which suggested that VDAC1 might involve in the regulation of the polarization of macrophages. NK cells are a specialized type of immune cells that can kill adjacent cells with surface markers related to oncogenic transformation. In the past decades, the fields of NK cell-based cancer immunotherapy have grown exponentially . The negative correction between VDAC1 expression and NK cell infiltration, as well as their markers CD7, indicated that VDAC1 might inhibit the activation of NK cells. These results suggested that overexpressed VDAC1 seemed to dampen tumor immunity thus contributing to the tumorigenesis and development of BC.
With the development of immunotherapy, ICI therapy has exhibited great efficacy in breast cancer treatment but only a small number of patients benefit from it . In our study, we analyzed the relationship between eight ICGs and VDAC1 expression, among which PDCD1, CTLA4, LAG3, SIGLEC15, and TIGIT were negatively corrected with VDAC1 expression in BC. The results suggested that BC patients with low VDAC1 expression might receive greater benefits from ICI therapy with a better prognosis. However, TIDE scores between low and high expression groups showed no difference. VDAC1 was not a reliable biomarker to predict the response to ICI therapy.
All in all, our study demonstrated the overexpression of VDAC1 in BC, which might be served as a novel biomarker for the diagnosis of BC patients. High expression of VDAC1 was associated with poor prognosis and VDAC1 was an independent factor for poor outcome prediction of BC patients. High expression of VDAC1 was closely associated with low infiltration levels of most types of immune cells. Our study revealed that VDAC1 might inhibit tumor immunity and might be served as a novel therapeutic target in BC.
Carbohydrate antigen 153
Voltage-dependent anion channel
Outer mitochondrial membrane
Tumor Immune Estimation Resource
The Cancer Genome Atlas
Transcripts Per Million
Gene Expression Omnibus
Receiver operating characteristics
Kyoto Encyclopedia of Genes and Genomes
Estimation of Stromal and Immune cells in Malignant Tumor tissues using Expression
Single sample GSEA
Microenvironment Cell Populations-counter
Immune checkpoint inhibitor
Immune checkpoint genes
Tumor immune dysfunction and exclusion
Human epidermal growth factor receptor 2
Area under the curve
Ferlay J, Colombet M, Soerjomataram I, et al. Cancer statistics for the year 2020: an overview. Int J Cancer. 2021. https://doi.org/10.1002/ijc.33588.
DeSantis CE, Ma J, Goding Sauer A, et al. Breast cancer statistics, 2017, racial disparity in mortality by state. CA Cancer J Clin. 2017;67(6):439–48. https://doi.org/10.3322/caac.21412.
Bray F, Ferlay J, Soerjomataram I, et al. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424. https://doi.org/10.3322/caac.21492.
Tang S, Zhou F, Sun Y, et al. CEA in breast ductal secretions as a promising biomarker for the diagnosis of breast cancer: a systematic review and meta-analysis. Breast Cancer. 2016;23(6):813–9. https://doi.org/10.1007/s12282-016-0680-9.
Li X, Xu Y, Zhang L. Serum CA153 as biomarker for cancer and noncancer diseases. Prog Mol Biol Transl Sci. 2019;162:265–76. https://doi.org/10.1016/bs.pmbts.2019.01.005.
Annesley SJ, Fisher PR. Mitochondria in health and disease. Cells. 2019;8(7):680. https://doi.org/10.3390/cells8070680.
Shoshan-Barmatz V, De S, Meir A. The mitochondrial voltage-dependent anion channel 1, Ca2+ transport, apoptosis, and their regulation. Front Oncol. 2017;7:60. https://doi.org/10.3389/fonc.2017.00060.
Shoshan-Barmatz V, Shteinfer-Kuzmine A, Verma A. VDAC1 at the intersection of cell metabolism, apoptosis, and diseases. Biomolecules. 2020;10(11):1485. https://doi.org/10.3390/biom10111485.
Shoshan-Barmatz V, Krelin Y, Shteinfer-Kuzmine A, et al. Voltage-dependent anion channel 1 as an emerging drug target for novel anti-cancer therapeutics. Front Oncol. 2017;7:154. https://doi.org/10.3389/fonc.2017.00154.
De Pinto V, Reina S, Gupta A, et al. Role of cysteines in mammalian VDAC isoforms’ function. Biochim Biophys Acta. 2016;1857(8):1219–27. https://doi.org/10.1016/j.bbabio.2016.02.020.
Magrì A, Reina S, De Pinto V. VDAC1 as pharmacological target in cancer and neurodegeneration: focus on its role in apoptosis. Front Chem. 2018;6:108. https://doi.org/10.3389/fchem.2018.00108.
Li T, Fan J, Wang B, et al. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77(21):e108–10. https://doi.org/10.1158/0008-5472.CAN-17-0307.
Kretschmer C, Conradi A, Kemmner W, et al. Latent transforming growth factor binding protein 4 (LTBP4) is downregulated in mouse and human DCIS and mammary carcinomas. Cell Oncol (Dordr). 2011;34(5):419–34. https://doi.org/10.1007/s13402-011-0023-y.
Lian ZQ, Wang Q, Li WP, et al. Screening of significantly hypermethylated genes in breast cancer using microarray-based methylated-CpG island recovery assay and identification of their expression levels. Int J Oncol. 2012;41(2):629–38. https://doi.org/10.3892/ijo.2012.1464.
Harvell DM, Kim J, O'Brien J, et al. Genomic signatures of pregnancy-associated breast cancer epithelia and stroma and their regulation by estrogens and progesterone. Horm Cancer. 2013;3:140–53. https://doi.org/10.1007/s12672-013-0136-z.
Hall P, Ploner A, Bjöhle J, et al. Hormone-replacement therapy influences gene expression profiles and is associated with breast-cancer prognosis: a cohort study. BMC Med. 2006;4:16. https://doi.org/10.1186/1741-7015-4-16.
Vasaikar SV, Straub P, Wang J, et al. LinkedOmics: analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res. 2018;46(D1):D956–63. https://doi.org/10.1093/nar/gkx1090.
Zhou Y, Zhou B, Pache L, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10(1):1523. https://doi.org/10.1038/s41467-019-09234-6.
Szklarczyk D, Gable AL, Lyon D, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13. https://doi.org/10.1093/nar/gky1131.
Hänzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics. 2013;14:7. https://doi.org/10.1186/1471-2105-14-7.
Becht E, Giraldo NA, Lacroix L, et al. Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression. Genome Biol. 2016;17:218. https://doi.org/10.1186/s13059-016-1070-5.
Yoshihara K, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612. https://doi.org/10.1038/ncomms3612.
Yuasa T, Masuda H, Yamamoto S, Numao N, Yonese J. Biomarkers to predict prognosis and response to checkpoint inhibitors. Int. J. Clin. Oncol. 2017;22:629–34. https://doi.org/10.1007/s10147-017-1122-1.
Jiang P, et al. Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response. Nat Med. 2018;24(10):1550–8. https://doi.org/10.1038/s41591-018-0136-1.
Li Y, Lu S, Zhang Y, et al. Loco-regional recurrence trend and prognosis in young women with breast cancer according to molecular subtypes: analysis of 1099 cases. World J Surg Oncol. 2021;19(1):113. https://doi.org/10.1186/s12957-021-02214-5.
Kim J. In silico analysis of differentially expressed genesets in metastatic breast cancer identifies potential prognostic biomarkers. World J Surg Oncol. 2021;19(1):188. https://doi.org/10.1186/s12957-021-02301-7.
Zhang YP, Na WT, Dai XQ, et al. Over-expression of SRD5A3 and its prognostic significance in breast cancer. World J Surg Oncol. 2021;19(1):260. https://doi.org/10.1186/s12957-021-02377-1.
Du Y, Miao Z, Wang K, et al. Expression levels and clinical values of miR-92b-3p in breast cancer. World J Surg Oncol. 2021;19(1):239. https://doi.org/10.1186/s12957-021-02347-7.
Zhou X, Xiao C, Han T, et al. Prognostic biomarkers related to breast cancer recurrence identified based on Logit model analysis. World J Surg Oncol. 2020;18(1):254. https://doi.org/10.1186/s12957-020-02026-z.
Mao XH, Ye Q, Zhang GB, et al. Identification of differentially methylated genes as diagnostic and prognostic biomarkers of breast cancer. World J Surg Oncol. 2021;19(1):29. https://doi.org/10.1186/s12957-021-02124-6.
Du XW, Li G, Liu J, et al. Comprehensive analysis of the cancer driver genes in breast cancer demonstrates their roles in cancer prognosis and tumor microenvironment. World J Surg Oncol. 2021;19(1):273. https://doi.org/10.1186/s12957-021-02387-z.
Ji Y, Wang H. Prognostic prediction of systemic immune-inflammation index for patients with gynecological and breast cancers: a meta-analysis. World J Surg Oncol. 2020;18(1):197. https://doi.org/10.1186/s12957-020-01974-w.
Morisaki T, Kashiwagi S, Asano Y, et al. Prediction of survival after eribulin chemotherapy for breast cancer by absolute lymphocyte counts and progression types. World J Surg Oncol. 2021;19(1):324. https://doi.org/10.1186/s12957-021-02441-w.
Alan O, Akin Telli T, Aktas B, et al. Is insulin resistance a predictor for complete response in breast cancer patients who underwent neoadjuvant treatment? World J Surg Oncol. 2020;18(1):242. https://doi.org/10.1186/s12957-020-02019-y.
Tang L, Shu X, Tu G. Exploring the influencing factors of the pathologic complete response in estrogen receptor-positive, HER2-negative breast cancer after neoadjuvant chemotherapy: a retrospective study. World J Surg Oncol. 2022;20(1):27. https://doi.org/10.1186/s12957-022-02492-7.
Chen F, Yin S, Luo B, et al. VDAC1 conversely correlates with Cytc expression and predicts poor prognosis in human breast cancer patients. Oxid Med Cell Longev. 2021;2021:7647139. https://doi.org/10.1155/2021/7647139.
Pittala S, Krelin Y, Shoshan-Barmatz V. Targeting liver cancer and associated pathologies in mice with a mitochondrial VDAC1-based peptide. Neoplasia. 2018;20(6):594–609. https://doi.org/10.1016/j.neo.2018.02.012.
Zhang G, Jiang G, Wang C, et al. Decreased expression of microRNA-320a promotes proliferation and invasion of non-small cell lung cancer cells by increasing VDAC1 expression. Oncotarget. 2016;7(31):49470–80. https://doi.org/10.18632/oncotarget.9943.
Gao W, Xua J, Wang F, et al. Mitochondrial proteomics approach reveals voltage-dependent anion channel 1 (VDAC1) as a potential biomarker of gastric cancer. Cell Physiol Biochem. 2015;37(6):2339–54. https://doi.org/10.1159/000438588.
Wang W, Zhang T, Zhao W, et al. A single talent immunogenic membrane antigen and novel prognostic predictor: voltage-dependent anion channel 1 (VDAC1) in pancreatic cancer. Sci Rep. 2016;6:33648. https://doi.org/10.1038/srep33648.
Wu CH, Lin YW, Wu TF, et al. Clinical implication of voltage-dependent anion channel 1 in uterine cervical cancer and its action on cervical cancer cells. Oncotarget. 2016;7(4):4210–25. https://doi.org/10.18632/oncotarget.6704.
Wang F, Qiang Y, Zhu L, et al. MicroRNA-7 downregulates the oncogene VDAC1 to influence hepatocellular carcinoma proliferation and metastasis. Tumour Biol. 2016;37(8):10235–46. https://doi.org/10.1007/s13277-016-4836-1.
Seyfried TN. Cancer as a mitochondrial metabolic disease. Front Cell Dev Biol. 2015;3:43. https://doi.org/10.3389/fcell.2015.00043.
Yang G, Zhou D, Li J, et al. VDAC1 is regulated by BRD4 and contributes to JQ1 resistance in breast cancer. Oncol Lett. 2019;18(3):2340–7. https://doi.org/10.3892/ol.2019.10534.
Seo JH, Chae YC, Kossenkov AV, et al. MFF regulation of mitochondrial cell death is a therapeutic target in cancer. Cancer Res. 2019;79(24):6215–26. https://doi.org/10.1158/0008-5472.CAN-19-1982.
Ikeda H, Kakeya H. Targeting hypoxia-inducible factor 1 (HIF-1) signaling with natural products toward cancer chemotherapy. J Antibiot (Tokyo). 2021;74(10):687–95. https://doi.org/10.1038/s41429-021-00451-0.
Fu C, Jiang A. Dendritic cells and CD8 T cell immunity in tumor microenvironment. Front Immunol. 2018;9:3059. https://doi.org/10.3389/fimmu.2018.03059.
Qu X, Tang Y, Hua S. Immunological approaches towards cancer and inflammation: a cross talk. Front Immunol. 2018;9:563. https://doi.org/10.3389/fimmu.2018.00563.
Myers JA, Miller JS. Exploring the NK cell platform for cancer immunotherapy. Nat Rev Clin Oncol. 2021;18(2):85–100. https://doi.org/10.1038/s41571-020-0426-7.
Setordzi P, Chang X, Liu Z, et al. The recent advances of PD-1 and PD-L1 checkpoint signaling inhibition for breast cancer immunotherapy. Eur J Pharmacol. 2021;895:173867. https://doi.org/10.1016/j.ejphar.2021.173867.
This work was supported by the Science and Technology Project of Shantou, China (Grant Nos. 170828211930352;120), funds from the Clinical Research Fund Project of Wu Jieping Medical Foundation, and funds from the Guangdong Provincial Key Laboratory of Breast Cancer Diagnosis and Treatment Research.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Fang, Y., Liu, J., Zhang, Q. et al. Overexpressed VDAC1 in breast cancer as a novel prognostic biomarker and correlates with immune infiltrates. World J Surg Onc 20, 211 (2022). https://doi.org/10.1186/s12957-022-02667-2
- Breast cancer
- Immune cells