VDAC1 Conversely Correlates with Cytc Expression and Predicts Poor Prognosis in Human Breast Cancer Patients

Chen, Fangfang; Yin, Shuai; Luo, Bin; Wu, Xiaoyan; Yan, Honglin; Yan, Dandan; Chen, Chuang; Guan, Feng; Yuan, Jingping

doi:https://doi.org/10.1155/2021/7647139

Oxidative Medicine and Cellular Longevity

On this page

Abstract Introduction Methods Results Discussion Conclusion Data Availability Ethical Approval Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Special Issue

Role of Redox Homeostasis in Cancer Biology and Anticancer Therapy

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 7647139 | https://doi.org/10.1155/2021/7647139

VDAC1 Conversely Correlates with Cytc Expression and Predicts Poor Prognosis in Human Breast Cancer Patients

Fangfang Chen,¹Shuai Yin,²Bin Luo,¹Xiaoyan Wu,¹Honglin Yan,¹Dandan Yan,¹Chuang Chen,³Feng Guan,¹and Jingping Yuan¹

Academic Editor: Bin Duan

Received28 Jul 2020

Revised02 Jan 2021

Accepted04 Feb 2021

Published17 Feb 2021

Abstract

Aim. The main objective of this article was to evaluate the association of voltage-dependent anion channel 1 (VDAC1) with Cytochrome C (Cytc) expression, various clinicopathological features, and prognosis in breast cancer (BC) patients. Meanwhile, the correlation of Cytc expression with various clinical features and 5-year disease-free survival (5-DFS) of BC was also investigated. Methods. In vivo, expression of VDAC1 and Cytc was examined in 219 BC tissues and 100 benign breast lesions by immunohistochemical (IHC) analysis. In vitro, MTT and wound healing migration assay were performed to detect the effect of VDAC1 on BC cells. Results. Expression of VDAC1 is conversely associated with Cytc in BC (), especially in triple-negative breast cancer (TNBC) (). Knockdown of VDAC1 inhibited proliferation () and migration () of MCF-7 cells. High expression of VDAC1 and low expression of Cytc had a significant association with multiple clinicopathological parameters () and poor 5-DFS () in BC. Conclusion. VDAC1 was elevated in BC tissues and conversely associated with Cytc. Detection of VDAC1 may provide guidance for the poor prognosis of BC, especially TNBC.

1. Introduction

Breast cancer (BC) is by far the most common female malignant tumor in the world, with a high mortality rate, and is the leading cause of cancer-related death in women [1]. In China, 367900 new cases of BC were diagnosed in 2018, which accounts for 19.2% of all newly diagnosed female cancers [2]. In recent years, although breast conserving surgery, targeted therapy, endocrine therapy, and immunotherapy have made progress in breast cancer treatment, the incidence rate of the disease has increased [3]. And only limited success has been achieved in cases of advanced cancer [4]. The burden of this disease is heavy, and it remains a serious health threat for women relative to the large Chinese population [5]. Therefore, it is extremely important to find new prognostic markers and early diagnosis of BC, so that patients can receive better treatments.

Cancer cells have a series of common characteristics, including high proliferation and antiapoptosis [6, 7]. The formation of cancer is related to cell metabolic reprogramming, such as enhanced aerobic glycolysis (Warburg effect) [8, 9]. And the “Warburg effect” is postulated to be closely associated with the closure of voltage-dependent anion channel (VDAC) [10, 11]. VDAC has been identified as a 31 kDa pore-forming protein, which is present in the mitochondrial outermembrane (OMM) of all eukaryotes [12]. Three homologous genes encode three VDAC isoforms: VDAC1, VDAC2, and VDAC3 [13]. Among them, VDAC1 is the most abundantly expressed and best characterized one, which forms a channel for the entry and exit of mitochondrial metabolites including ATP and NADH across the OMM, thereby regulating the activity of mitochondria [14]. VDAC1 plays a key role in the regulation of apoptosis by both interaction with proapoptotic (Bcl2 and Bcl-xL) proteins and antiapoptotic proteins (Bax, Bak, and Bim) [15, 16]. It has been proved that VDAC1 is involved in several diseases, such as neurodegenerative, cardiac injury, and neoplastic [17]. Several studies have demonstrated that VDAC1 is remarkable expressed in malignant tumors such as uterine cervical cancer [18], hepatocellular carcinoma [19], and cholangiocellular carcinoma [20], indicating that it plays a significant role in high energy-demanding cancer cells. Furthermore, VDAC1 has also been reported to promote tumor growth and play a controversial role in the prognosis of different cancers [21]. A number of anticancer agents have taken VDAC1 as a novel pharmacological target [12]. However, the correlation of VDAC1 expression with tumorigenesis and progression of solid tumor of BC and the role of VDAC1 acts in the prognosis of BC patients has not been well-documented.

Cytochrome C (Cytc) has been proposed to be a proapoptotic factor which is located in the inner membrane of mitochondria and plays an irreplaceable role in mitochondrial electron transport and intrinsic type II apoptosis [22, 23]. Cytc release forms an essential step in the apoptotic cascade, upon binding with apoptosis protease-activating factor 1 (Apaf-1), dATP, and procaspase-9 [24, 25]. A number of studies have shown that Cytc can induce apoptosis of cancer cells [26, 27]. Jamsheed et al. [28] demonstrated that in non-small-cell lung cancer patients, Cytc level was lower than healthy individuals, and the lower expression of Cytc were associated with advanced stages, high grade histological differentiation, and shorter overall survival. Rahul et al. [29] explored that in prostate cancer, the deficiency of Cytc contributed to tumor invasiveness and therapeutic resistance and led to faster recurrence. It has been widely reported that VDAC1 is involved in the release of Cytc, which is able to signal Cytc to initiate the mitochondrial-mediated cell death cascade [30]. In contrast, VDAC1 has also been reported to interact with antiapoptotic proteins such as Bcl-2 and hexokinase (HK) to control the release of Cytc [15]. In several melanoma and prostate cancer cell lines, the expression level of VDAC1 is related to the induction of Cytc release, which provides more possibilities for the pharmacological treatment of tumors [31]. However, the association of VDAC1 expression with Cytc in BC is still elusive and rarely addressed, thereby indicating further investigation is needed.

The main objectives of this article were as follows: (i) evaluate the association of VDAC1 expression with Cytc in BC; (ii) investigate the correlation of various clinical features and 5-year disease-free survival (5-DFS) of BC with VDAC1 and Cytc, respectively, and (iii) explore the effects of VDAC1 on cell proliferation and migration in BC cell line.

2. Material and Methods

2.1. Patient Tissue Samples

A total of 219 formalin-fixed, paraffin-embedded primary invasive breast cancer tissue samples were collected from the patients diagnosed with BC through histopathologic evaluation on surgical tissue specimens. 100 cases of benign breast lesions were collected as controls. All the patients underwent surgical treatment at Renmin Hospital of Wuhan University between August 2009 and December 2010. There were no any previous chemotherapies, radiotherapies, or other treatments before surgery in these patients. The patients’ written informed consent was obtained before the operation, and the study was approved by the Ethics Committee of Renmin Hosptial of Wuhan University. Patients were all followed up for 5 years. The follow-up data was calculated as the period from the date of surgery to the end of follow-up or death. We followed up all the patients by telephone interviews or outpatient clinic visits.

2.2. Cell Culture

The breast cancer cell line MCF-7 was purchased from the Cell Bank of Chinese Academy of Sciences (Shanghai, China). These cells were incubated in Dulbecco’s modified Eagle’s medium (DMEM, Invitrogen, Carlsbad, CA, USA) supplemented with 10% fetal bovine serum (FBS, Invitrogen), routinely maintained at 37°C and in an atmosphere of 5% CO2.

2.3. Immunohistochemistry (IHC)

A tissue array was used which included two tumor samples from each patient. The paraffin-embedded tissues were cut into 4 μm thick sections, then deparaffinized and dehydrated following standard procedures. Subsequently, paraffin sections were rinsed with PBS () and then blocked with 3% hydrogen peroxide at 37°C for endogenous peroxidase ablation for 10 min. Antigen retrieval was conducted by microwave heating with citrate buffer (pH 6.0) for 20 min. Then, the samples were exposed to normal goat serum at 37°C for 20 min to decrease nonspecific antibody binding. The tissue sections were incubated overnight at 4°C with the primary antibody (anti-VDAC1, 1 : 1000, ab15895, Abcam, UK; anti-Cytc, 1 : 50, sc-13561, Santa Cruz Biotechnology, Inc., Santa Cruz, CA, USA). After rinsing in PBS, the tissue sections were incubated with horseradish peroxidase-labeled anti-rabbit antibodies at 37^∘C for 20 min. Then, the tissue sections were rinsed with PBS for 4 times and then dripped with freshly prepared 3,3-diaminobenzidine (DAB). Microscopically, the staining was terminated when the tissue sections were brown-yellow or brown. Subsequently, all the tissue sections were restained with hematoxylin for about 1 min. Finally, the slices were dehydrated with ethanol and toluene and then sealed with neutral gum. PBS was used to replace the primary antibody as negative control.

2.4. Evaluation of Immunohistochemical Staining

The slides were viewed via Olympus BX53 (Tokyo, Japan) microscope. IHC staining was evaluated independently by two pathologists under the double-blind condition. VDAC1 was mainly expressed in membrane of tumor cells. VDAC1 immunohistochemical staining in tumor cells was evaluated semiquantitatively as follows: (1) staining intensity: 0 (no staining), 1 (weak staining), 2 (moderate staining), and 3 (strong staining); (2) the extent of staining: 0 (0%), 1 (1-20%), 2 (21-50%), and 3 (>50%). Five most representative fields of high magnification (400x) were selected to calculate the final score. The final immunohistochemical score was multiplied by staining intensity and extent, theoretically from 0 to 9. Scores less than or equal to 3 were defined as low expression, and scores greater than or equal to 4 were described as high expression. Cytc protein was predominantly located in membrane and cytoplasm of tumor cells. The staining intensity was classified as four grades as follows: 0 (no staining), 1 (light yellow), 2 (brown-yellow), and 3 (dark brown). The percentage of positive cells was classified as five grades as follows: 0 (0%), 1 (≤30%), 2 (31-50%), 3 (51-80%), and 4(≥80%). Five most representative fields of high magnification (400x) were selected to calculate the final score. The final immunohistochemical score was the product of staining intensity and extent, theoretically from 0 to 12. Scores less than 4 were defined as low expression, and scores greater than or equal to 4 were described as high expression.

2.5. Transfection

Small interfering RNA (siRNA) duplexes targeting human VDAC1 (si-VDAC1-1: GTCTAGGACTGGAATTTCA, si-VDAC1-2: GGAGACCGCTGTCAATCTT, and si-VDAC1-3: GCTGCGACATGGATTTCGA) were designed and synthesized by Guangzhou Ruibo Biotechnology Co., Ltd (Guangzhou, China). SiRNA duplexes with nonspecific sequences were used as siRNA negative control. The transfections were carried out using Lipofectamine-RNAi MAX (Invitrogen) following the manufacturer’s protocol. The siRNA-targeted human VDAC1 is designed from the messenger RNA (mRNA) sequences of human VDAC1 gene (RefSeq: NM_003374). The effect of silencing was verified at the protein level.

2.6. Western Blot Analysis

Cultured cells were collected and lysed to harvest protein. Denatured protein samples were electrophoretically separated on 10% SDS-polyacrylamide gel (PAGE), and then, the proteins were transferred to polyvinylidenefluoride (PVDF) membranes. After transfer, the membranes were incubated with a blocking buffer consisting of 50 mm Tris-HCl (pH 7.5), 100 mm NaCl, 5% nonfat dry milk, and 0.05% Tween 20 for 2 h. Next, the membranes were incubated with primary antibodies against VDAC1 (Abcam Ab15895, 1 : 1000) and Cytc (Santa Cruz sc-13561, 1 : 100) at 4°C overnight. After being washed for three times, the membranes were further incubated with horseradish peroxidase-conjugated sheep antihuman (Amersham Biosciences, Piscataway, NJ) IgG antibodies at a dilution of 1 : 5000 for 2 h at room temperature. Protein was visualized using an enhanced chemiluminescence system (ECL) reagent (KeyGEN BioTECH, China). The volumes of target bands were normalized to GAPDH. Then, the membranes were detected by automatic chemiluminescence analyzer, and the gray value of related bands was analyzed by Image J software (Broken Symmetry Software). Each experiment was repeated at least three times.

2.7. MTT Assay

Cell proliferation viability of MCF-7 cells was measured using the MTT assay. In brief, cells were seeded after transfection with siRNA for the indicated time in triplicate in 96-well plates at densities of . The cells were incubated under standard conditions for 48 hours. Then, cells were treated by MTT dye (5 mg/ml 10 μl per well) for 4 h. The absorbance values were measured using a microplate reader at a wavelength of 490 nm. The growth assays were repeated three times and reported as percentage changes compared to the controls.

2.8. Wound Healing Migration Assay

A total of were plated into a 6-well plate, and complete convergence was allowed. Wounds were created by scraping confluent cell monolayers with 200 μl pipette tips. The cells were washed with PBS twice, after which cells were cultured in 10% FBS DMEM. Three digital images were taken at 0 and 48 h after wounding under a microscope (Olympus, Japan). The assays were accomplished in triplicate and repeated three times. The change of wound width was measured, and wound area filled was calculated by Image J software (Broken Symmetry Software).

2.9. Statistical Analysis

SPSS software version 17.0 was used to carry out all the statistical analyses. Chi-square test and Student’s test were used to analyze the association between categorical variables. Spearman’s rank correlation analysis performed to evaluate correlations between variables. The survival curves were disposed by using the Kaplan-Meier method and log-rank test. We performed univariate and multivariate survival analysis through Cox proportional hazard regression model to assess the independent prognostic factors in BC patients. Hazard ratios (HRs) and their 95% confidence intervals (CIs) were calculated for both univariate and multivariate analyses. Two-tailed values of <0.05 were considered statistically significant. Graph-Pad Prism 8.0.1 (GraphPad, San Diego, CA) software was used to present graphs.

3. Results

3.1. The Characteristics of Patients

The patients were composed of 219 females with a median age of 50 (age range, 29–78) years. Among the 219 cases, 165 (75.34%) were classified into G1 and G2 stages, 54 (24.66%) were defined with G3 on the basis of histological differentiation. 153 (69.86%) patients were classified as stages I and II, and 66 (30.14%) as III on the basis of tumor node metastasis (TNM) stage. Other basic clinicopathological characteristics, including age, menopausal status, lymph node metastasis (LNM), estrogen receptor (ER), progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2), and recurrence, were presented in Table 1.

3.2. Expression of VDAC1 Is Significantly Higher in BC Tissues than in Benign Breast Lesions

To examine the expression level of VDAC1 protein, we performed IHC on 219 cases of BC tissues and 100 cases of the benign breast lesions. As shown in Figure 1, VDAC1 protein was mainly expressed in the membrane of breast cancer cells. Table 2 showed the result of IHC staining of VDAC1 protein. Of 219 BC samples, 162 (73.97%) showed high expression of VDAC1 protein, and 57 (26.03%) showed low expression. In benign breast lesions, 55 (55%) showed high expression of VDAC1 protein, and 45 (45%) showed low expression. The expression of VDAC1 was significantly higher in BC tissues (, ) as determined by chi-square test.

(a)

(b)

(c)

(d)

3.3. Expression of Cytc Is Lower in BC Tissues than in Benign Breast Lesions

We examined the expression of Cytc by IHC in 219 cases of BC tissue and 100 cases of the benign breast lesions. As shown in Figure 2, Cytc protein was mainly expressed in the membrane and cytoplasm of BC tumor cells. The result of IHC staining of Cytc protein was summarized in Table 2, which showed that the rate of high Cytc expression was 29.68% (65/219) and the rate of low Cytc expression was 70.32% (154/219) in all BC samples. While in benign breast lesions, the rate of high Cytc expression was 47% (47/100), and the rate of low Cytc expression was 53% (53/100). Thus, the expression of Cytc was significantly lower in BC tissues than benign breast lesions (, ).

(a)

(b)

(c)

(d)

3.4. VDAC1 Protein Expression Correlates with Cytc Protein Expression and Clinicopathological Parameters in BC Tissues, but Not in Benign Breast Lesions

In the benign breast lesions, there was no significant association between VDAC1 expression and Cytc protein expression (, ; Table 3). Conversely, as shown in Table 4, high expression of VDAC1 protein was inversely associated with Cytc protein expression in BC tissues (, , ), which was also shown in Figure 3. Furthermore, as shown in Table 5, high expression of VDAC1 protein was detected in 73.97% (162/219) of BC tissues, which was significantly associated with advanced TNM stage (, ), higher histological grade (, ), recurrence (, ), HER2 gene amplification (, ), and lymph node metastasis (, ), but not with other examined clinicopathological parameters, including age (, ), ER status (, ), PR status (, ), and menopause (, ).

(a)

(b)

(c)

(d)

3.5. Cytc Expression Is Also Correlated with Various Clinical Features in BC Tissues

As is shown in Table 5, low expression of Cytc was associated with higher histological grade (, ), ER status (, ), PR status (, ), and recurrence (, ), but not age (, ), TNM stage (, ), lymph node metastasis (, ), HER2 gene amplification (, ), or menopause (, ) in tumor samples.

3.6. Expression Level of Cytc Protein Is Increased after VDAC1 Knockdown In Vitro

Firstly, we conducted transient silencing of VDAC1 in the MCF-7 cell line with three different siRNAs and verified the knockdown efficiency with Western blot. Si-VDAC1-1 achieved more effective knockdown efficiency, so it was used in all subsequent experiments (, Figure 4(a)).

(a)

(b)

(c)

As VDAC1 was conversely associated with Cytc expression in BC tissues, we also explored the expression level of Cytc in vitro after knockdown of VDAC1. Consistently, silencing of VDAC1 led to Cytc expression upregulation at the protein level in MCF-7 cells (, Figure 4(b)).

3.7. VDAC1 Promoted the Proliferation and Migration of BC Cells In Vitro

Additionally, we explored whether VDAC1 could affect the cell proliferation and migration. After transfection with si-VDAC1-1, a significant reduction in the proliferation rate of MCF-7 cells was observed with the MTT assay compared with the control (, Figure 4(c)). Afterwards, wound healing migration assay was performed. As a result, we validated that VDAC1 silencing obviously inhibited MCF-7 cell migration compared with the controls (, Figure 5). Therefore, knockdown of VDAC1 can suppress the proliferative and migrative ability of MCF-7 cells in vitro.

(a)

(b)

3.8. Correlation Analysis of the 5-DFS with the Expression of VDAC1, Cytc Protein, and Other Parameters in BC

Kaplan-Meier analysis showed that low VDAC1 protein expression in BC predicted a better survival and lower mortality rate (Figure 6(a), ). Similarly, high expression of Cytc also predicted a better outcome of BC patients (Figure 6(b), ). Univariate analysis of predictive factors for the 5-DFS in BC patients was performed by Cox proportional hazards regression model (Table 6). In univariate analysis, histological grade (), TNM stage (), ER status (), PR status (), and lymph node metastasis () were also significantly correlated with the 5-DFS of BC patients. However, age, HER2 gene amplification, and menopause had no significant association with 5-DFS in BC patients () (Table 6). In order to analyze whether the above univariate was an independent prognostic factor, multivariate Cox proportional hazard model for 5-DFS was performed. The results indicated that both the expressions of VDAC1 and Cytc were independent prognostic parameters for 5-DFS of BC patients [HR: 3.982 (1.723-9.207), ; HR: 0.542 (0.307-0.959), , respectively, Table 6]. Concurrently, TNM stage, PR status, lymph node metastasis, and histological grade were also independent predictors regarding the 5-DFS of breast cancer patients.

(a)

(b)

(c)

(d)

3.9. High Expression of VDAC1 Is Associated with the Poor Prognosis of HER2-Negative Breast Cancer

As is shown in Table 5, the expression of VDAC1 was significantly associated with HER2 gene, but not with ER and PR status. As targeted therapy can be implemented for HER2-positive breast cancer, we analyzed the correlation of VDAC1 protein expression with the prognosis of HER2-positive breast cancer patients and HER2-negative breast cancer patients, respectively. Kaplan-Meier analysis showed that VDAC1 expression had no correlation with 5-DFS in HER2-positive breast cancer patients (). However, a significant relevance was found between high VDAC1 expression and shorter 5-DFS () in HER2-negative breast cancer patients (Figure 6(c)). As triple-negative breast cancer (TNBC) is a special molecular subtype of HER2-negative breast cancer, which is defined by the absence of the ER, PR, and HER2 genes, with no standard treatment at present, we then explored the correlation of VDAC1 expression with the prognosis of TNBC. Interestingly, VDAC1 was also associated with reduced 5-DFS in TNBC () (Figure 6(d)), predicting a poorer survival and higher mortality rate. To further validate the prognostic significance of VDAC1, Cox proportional hazards model analysis was also performed in TNBC. As presented in Table 7, VDAC1 was an independent predictor of poor 5-DFS [HR: 4.018 (1.415-11.407), ]. Meanwhile, VDAC1 protein expression was inversely associated with Cytc in TNBC (, , , Table 8, Figure 7). In summary, high expression of VDAC1 is associated with the prognosis of TNBC and inversely correlates with Cytc.

(a)

(b)

(c)

(d)

4. Discussion

In this study, we observed that VDAC1 was elevated while Cytc was decreased in breast carcinoma patients compared with benign breast lesions. VDAC1 protein expression was conversely associated with Cytc in BC, especially in TNBC. Meanwhile, knockdown of VDAC1 inhibits BC cell proliferation and migration in vitro. Furthermore, both high VDAC1 and low Cytc protein expression had significant positive correlation with poor prognosis. VDAC1 expression can also be an independent prognostic factor of BC, especially TNBC.

VDAC1 participates in cancer metabolism via its modulatory roles in the transport of various metabolites [32]. In many types of cancer, the interaction of VDAC1 with HK, especially HK II, directly accesses to mitochondrial ATP for phosphorylation of glucose to glucose-6-phosphate and contributes to the cancer cells unrestricted growth and the inhibiting of apoptosis [21, 33]. VDAC1 has been found to be involved in tumor proliferation, migration, metastasis, and invasion [13]. It has been demonstrated that VDAC1 acts a controvertial role in the prognosis of different malignant tumors. For example, in uterine cervical cancer, high expression of VDAC1 was associated with exhibited deeper stromal invasion, larger tumor size, higher recurrence, and poorer overall survival [18]. Conversely, in cholangiocellular carcinoma, the low expression of VDAC1 correlated with higher cancer stage classification, lymph node involvement, and reduced survival [20]. These conflicting consequences demonstrate the differential effects of VDAC1 expression in different kinds of cancer and may need further exploration. In the present study, 219 cases of primary invasive breast cancer tissues were collected, and it was found that VDAC1 protein was primarily located in the membrane of breast cancer cells. The expression of VDAC1 protein in breast cancer solid tumors was significantly higher than that in benign breast lesions, and high expression of VDAC1 protein correlated with advanced TNM stage, higher histological grade, recurrence, lymph node metastasis, and HER2 gene amplification, thereby suggesting that high VDAC1 expression may play a role in promoting tumorigenesis and progression of BC. Interestingly, our present results in breast cancer were consistent with the reports in pancreatic cancer [34] and colorectal cancer [35], which indicated VDAC1 expression was upregulated in tumor and promoted the growth and invasion of cancer cells. In addition, knockdown of VDAC1 in multiple types of cancer cell lines, including colon and lung cancer cells, has been demonstrated to block proliferation and migration of the cancer cells in vitro [36–38]. As is shown above, our studies achieved the similar results.

Furthermore, by the multivariate analysis, we found that overexpression of VDAC1 protein in BC tissues could be as an independent poor prognostic factor. Consistent with this result, Chih-Hsien and Eiran’s studies also showed that high expression of VDAC1 was correlated with poorer prognosis in uterine cervical cancer and hepatocellular carcinoma [18, 19]. Therefore, these data suggest that VDAC1 has the potential to be a poor prognostic marker in BC. As TNBC showed more progressively malignant manifestation with worse clinical outcomes and is the most insensitive subtype of breast cancer to drug treatment, we next explored the correlation of VDAC1 expression with the prognosis of TNBC. Amazingly, our result demonstrated that high level of VDAC1 protein was also associated with reduced 5-DFS and acted as an independent predictor of poor prognosis in TNBC, suggesting more potential use of VDAC1 should be exploited in prognostic marker and therapeutic target.

Cytc release from mitochondria is the driving force for apoptosome leading to apoptotic cell death in several malignant tumor [39], and VDAC1 has been widely reported to be interacted with pro- or antiapoptotic proteins such as Bcl-2 and HK, whereby mediating the release of Cytc [15, 40]. In many types of cancer cells, when HK2 binds to VDAC1, the interaction between VDAC1 and Bcl-2 protein family will be blocked, resulting in a decrease in Cytc release, thus protecting tumor cells from apoptosis [21]. In our present study, VDAC1 protein expression was inversely associated with Cytc in BC, especially in TNBC. Simultaneously, silencing of VDAC1 upregulated Cytc expression at the protein level in MCF-7 cells. Furthermore, Cytc was lower expressed in BC compared with benign breast lesions and low expression of Cytc was correlated with higher histological grade, ER status, PR status, and recurrence. A similar correlation was found in prostate cancer which Cytc deficiency contributed to tumor invasiveness and faster recurrence [29]. Importantly, we also found that decreased expression of Cytc was an independent prognostic factor and played a pivotal role in the poorer 5-DFS of that cohort of BC patients. As a result, contrary to VDAC1, Cytc has the potential to be an improved prognostic marker in BC.

It should be acknowledged that there are also some limitations in this study. Firstly, due to the limited number of patients in this study, a larger cohort is required. Secondly, we took the benign breast tissues as control, which was not paired comparison, so associations should be interpreted with cognizance of these possible differences. Besides, more investigations need to be conducted to explore the mechanisms by which VDAC1 reduces Cytc expression in BC. Afterwards, the prognostic significance of VDAC1 in BC was discussed only at histological level. The exact role of VDAC1 in BC, especially in TNBC, still needs to be evaluated in follow-up mechanistic investigations. Finally, our study evaluated prognosis by 5-DFS rather than OS. Since DFS was not correlated with OS sometimes, the influence of VDAC1 expression on OS is still a topic for future research.

5. Conclusion

Our study showed for the first time that VDAC1 was elevated in BC tissues. Meanwhile, our findings firstly revealed that VDAC1 expression was conversely associated with Cytc, and knockdown of VDAC1 inhibited malignant behavior of BC. Besides, high VDAC1 level was associated with reduced 5-DFS of BC patients. By the multivariate analysis, we found that overexpression of VDAC1 could be employed as an independent poor prognostic factor in breast cancer, including TNBC, which was intractable in clinical. All in all, VDAC1 can be exploited as a potential prognostic marker and therapeutic target in BC, especially TNBC.

Data Availability

The datasets used and analyzed during the current study are available from the corresponding authors on reasonable request.

Ethical Approval

This study was approved by the Ethical Committee of Renmin Hospital of Wuhan University (WDRY2019-K010). The written informed consents were obtained from all the patients.

Conflicts of Interest

The authors declare that they have no competing interests.

Authors’ Contributions

FC and JY made contributions to the conception and design of the study; CC and XW collected the samples and performed immunohistochemistry experiments; DY and HY made contributions to statistical analysis; BL and FG evaluated the immunohistochemical staining of all sections; SY made contributions to critical revision of the manuscript. All authors read and approved the final manuscript.

Acknowledgments

This work was supported by grants from the Science and Technology Planning Project of Wuhan (Grant numbers 2017060201010172) and Guidance Foundation of Renmin Hospital of Wuhan University (Grant numbers RMYD2018M27).

.

References

Global Burden of Disease Cancer Collaboration, C. Fitzmaurice, D. Abate et al., “Global, regional, and national cancer incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 29 cancer groups, 1990 to 2017: a systematic analysis for the global burden of disease study,” JAMA Oncology, vol. 5, no. 12, pp. 1749–1768, 2019.
View at: Publisher Site | Google Scholar
F. Bray, J. Ferlay, I. Soerjomataram, R. L. Siegel, L. A. Torre, and A. Jemal, “Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries,” CA: a Cancer Journal for Clinicians, vol. 68, no. 6, pp. 394–424, 2018.
View at: Publisher Site | Google Scholar
C. E. DeSantis, J. Ma, M. M. Gaudet et al., “Breast cancer statistics, 2019,” CA: a Cancer Journal for Clinicians, vol. 69, no. 6, pp. 438–451, 2019.
View at: Publisher Site | Google Scholar
A. Costa, Y. Kieffer, A. Scholer-Dahirel et al., “Fibroblast heterogeneity and immunosuppressive environment in human breast cancer,” Cancer Cell, vol. 33, no. 3, pp. 463–479.e10, 2018, e10.
View at: Publisher Site | Google Scholar
P. Ji, Y. Gong, C. C. Jiang, X. Hu, G. H. di, and Z. M. Shao, “Association between socioeconomic factors at diagnosis and survival in breast cancer: a population-based study,” Cancer Medicine, vol. 9, no. 5, pp. 1922–1936, 2020.
View at: Publisher Site | Google Scholar
D. Hanahan and R. A. Weinberg, “Hallmarks of cancer: the next generation,” Cell, vol. 144, no. 5, pp. 646–674, 2011.
View at: Publisher Site | Google Scholar
K. Meehan and L. J. Vella, “The contribution of tumour-derived exosomes to the hallmarks of cancer,” Critical Reviews in Clinical Laboratory Sciences, vol. 53, no. 2, pp. 121–131, 2015.
View at: Publisher Site | Google Scholar
W. H. Koppenol, P. L. Bounds, and C. V. Dang, “Otto Warburg's contributions to current concepts of cancer metabolism,” Nature Reviews. Cancer, vol. 11, no. 5, pp. 325–337, 2011.
View at: Publisher Site | Google Scholar
J. Lu, M. Tan, and Q. Cai, “The Warburg effect in tumor progression: mitochondrial oxidative metabolism as an anti-metastasis mechanism,” Cancer Letters, vol. 356, no. 2, pp. 156–164, 2015.
View at: Publisher Site | Google Scholar
V. V. Lemeshko, “VDAC electronics: 1. VDAC-hexo(gluco)kinase generator of the mitochondrial outer membrane potential,” Biochimica et Biophysica Acta, vol. 1838, no. 5, pp. 1362–1371, 2014.
View at: Publisher Site | Google Scholar
V. V. Lemeshko, “VDAC electronics: 2. A new, anaerobic mechanism of generation of the membrane potentials in mitochondria,” Biochimica et Biophysica Acta, vol. 1838, no. 7, pp. 1801–1808, 2014.
View at: Publisher Site | Google Scholar
V. Shoshan-Barmatz, Y. Krelin, A. Shteinfer-Kuzmine, and T. Arif, “Voltage-dependent anion channel 1 as an emerging drug target for novel anti-cancer therapeutics,” Frontiers in Oncology, vol. 7, p. 154, 2017.
View at: Publisher Site | Google Scholar
V. Shoshan-Barmatz, D. Ben-Hail, L. Admoni, Y. Krelin, and S. S. Tripathi, “The mitochondrial voltage-dependent anion channel 1 in tumor cells,” Biochimica et Biophysica Acta (BBA) - Biomembranes, vol. 1848, no. 10, pp. 2547–2575, 2015.
View at: Publisher Site | Google Scholar
V. Shoshan-Barmatz, S. Pittala, and D. Mizrachi, “VDAC1 and the TSPO: expression, interactions, and associated functions in health and disease states,” International Journal of Molecular Sciences, vol. 20, 13, no. 13, p. 3348, 2019.
View at: Publisher Site | Google Scholar
A. K. S. Camara, Y. F. Zhou, P. C. Wen, E. Tajkhorshid, and W. M. Kwok, “Mitochondrial VDAC1: a key gatekeeper as potential therapeutic target,” Frontiers in Physiology, vol. 8, p. 460, 2017.
View at: Publisher Site | Google Scholar
V. Shoshan-Barmatz, Y. Krelin, and A. Shteinfer-Kuzmine, “VDAC1 functions in Ca²⁺ homeostasis and cell life and death in health and disease,” Cell Calcium, vol. 69, pp. 81–100, 2018.
View at: Publisher Site | Google Scholar
V. Shoshan-Barmatz, E. Nahon-Crystal, A. Shteinfer-Kuzmine, and R. Gupta, “VDAC1, mitochondrial dysfunction, and Alzheimer's disease,” Pharmacological Research, vol. 131, pp. 87–101, 2018.
View at: Publisher Site | Google Scholar
C. H. Wu, Y. W. Lin, T. F. Wu, J. L. Ko, and P. H. Wang, “Clinical implication of voltage-dependent anion channel 1 in uterine cervical cancer and its action on cervical cancer cells,” Oncotarget, vol. 7, no. 4, pp. 4210–4225, 2016.
View at: Publisher Site | Google Scholar
F. Wang, Y. Qiang, L. Zhu et al., “MicroRNA-7 downregulates the oncogene VDAC1 to influence hepatocellular carcinoma proliferation and metastasis,” Tumour Biology, vol. 37, no. 8, pp. 10235–10246, 2016.
View at: Publisher Site | Google Scholar
R. G. Feichtinger, D. Neureiter, R. Kemmerling, J. A. Mayr, T. Kiesslich, and B. Kofler, “Low VDAC1 expression is associated with an aggressive phenotype and reduced overall patient survival in cholangiocellular carcinoma,” Cell, vol. 8, no. 6, p. 539, 2019.
View at: Publisher Site | Google Scholar
A. Magri, S. Reina, and V. De Pinto, “VDAC1 as pharmacological target in cancer and neurodegeneration: focus on its role in apoptosis,” Frontiers in Chemistry, vol. 6, p. 108, 2018.
View at: Publisher Site | Google Scholar
M. Zhang, J. Zheng, R. Nussinov, and B. Ma, “Release of cytochrome C from Bax pores at the mitochondrial membrane,” Scientific Reports, vol. 7, no. 1, p. 2635, 2017.
View at: Publisher Site | Google Scholar
F. J. Bock and S. W. G. Tait, “Mitochondria as multifaceted regulators of cell death,” Nature Reviews. Molecular Cell Biology, vol. 21, no. 2, pp. 85–100, 2020.
View at: Publisher Site | Google Scholar
R. Santucci, F. Sinibaldi, P. Cozza, F. Polticelli, and L. Fiorucci, “Cytochrome c: an extreme multifunctional protein with a key role in cell fate,” International Journal of Biological Macromolecules, vol. 136, pp. 1237–1246, 2019.
View at: Publisher Site | Google Scholar
M. Zhou, Y. Li, Q. Hu et al., “Atomic structure of the apoptosome: mechanism of cytochrome c- and dATP-mediated activation of Apaf-1,” Genes & Development, vol. 29, no. 22, pp. 2349–2361, 2015.
View at: Publisher Site | Google Scholar
N. Yadav, R. Gogada, J. O'Malley et al., “Molecular insights on cytochrome c and nucleotide regulation of apoptosome function and its implication in cancer,” Biochimica et Biophysica Acta (BBA) - Molecular Cell Research, vol. 1867, no. 1, article 118573, 2020.
View at: Publisher Site | Google Scholar
C. M. Figueroa, B. N. Suarez, A. M. Molina, J. C. Fernandez, Z. Torres, and K. Griebenow, “Smart release nano-formulation of cytochrome C and hyaluronic acid induces apoptosis in cancer cells,” Journal of Nanomedicine & Nanotechnology, vol. 8, no. 1, 2017.
View at: Publisher Site | Google Scholar
J. Javid, R. Mir, P. K. Julka, P. C. Ray, and A. Saxena, “Extracellular cytochrome c as a biomarker for monitoring therapeutic efficacy and prognosis of non-small cell lung cancer patients,” Tumour Biology, vol. 36, no. 6, pp. 4253–4260, 2015.
View at: Publisher Site | Google Scholar
R. Kumar, T. A. Bhat, E. M. Walsh et al., “Cytochromec deficiency confers apoptosome and mitochondrial dysfunction in African-American men with prostate cancer,” Cancer Research, vol. 79, no. 7, pp. 1353–1368, 2019.
View at: Publisher Site | Google Scholar
V. Shoshan-Barmatz, S. De, and A. Meir, “The mitochondrial voltage-dependent anion channel 1, Ca²⁺ transport, apoptosis, and their regulation,” Frontiers in Oncology, vol. 7, p. 60, 2017.
View at: Publisher Site | Google Scholar
J. C. Lai, W. Tan, L. Benimetskaya, P. Miller, M. Colombini, and C. A. Stein, “A pharmacologic target of G3139 in melanoma cells may be the mitochondrial VDAC,” Proceedings of the National Academy of Sciences of the United States of America, vol. 103, no. 19, pp. 7494–7499, 2006.
View at: Publisher Site | Google Scholar
V. Shoshan-Barmatz, Y. Krelin, and Q. Chen, “VDAC1 as a player in mitochondria-mediated apoptosis and target for modulating apoptosis,” Current Medicinal Chemistry, vol. 24, no. 40, pp. 4435–4446, 2017.
View at: Publisher Site | Google Scholar
V. Shoshan-Barmatz and M. Golan, “Mitochondrial VDAC1: function in cell life and death and a target for cancer therapy,” Current Medicinal Chemistry, vol. 19, no. 5, pp. 714–735, 2012.
View at: Publisher Site | Google Scholar
W. Wang, T. Zhang, W. Zhao et al., “A single talent immunogenic membrane antigen and novel prognostic predictor: voltage-dependent anion channel 1 (VDAC1) in pancreatic cancer,” Scientific Reports, vol. 6, no. 1, article 33648, 2016.
View at: Publisher Site | Google Scholar
X. Liu, B. He, T. Xu et al., “MiR-490-3p functions as a tumor suppressor by inhibiting oncogene VDAC1 expression in colorectal cancer,” Journal of Cancer, vol. 9, no. 7, pp. 1218–1230, 2018.
View at: Publisher Site | Google Scholar
T. Arif, L. Vasilkovsky, Y. Refaely, A. Konson, and V. Shoshan-Barmatz, “Silencing VDAC1 expression by siRNA inhibits cancer cell proliferation and tumor growth in vivo,” Molecular therapy - Nucleic acids, vol. 8, p. 493, 2017.
View at: Publisher Site | Google Scholar
S. Freitas, R. Martins, M. Costa et al., “Hierridin B isolated from a marine cyanobacterium alters VDAC1, mitochondrial activity, and cell cycle genes on HT-29 colon adenocarcinoma cells,” Marine Drugs, vol. 14, no. 9, p. 158, 2016.
View at: Publisher Site | Google Scholar
A. Maimaiti, A. Aili, H. Kuerban, and X. Li, “VDAC1 mediated anticancer activity of gallic acid in human lung adenocarcinoma A549 cells,” Anti-Cancer Agents in Medicinal Chemistry, vol. 18, no. 2, pp. 255–262, 2018.
View at: Publisher Site | Google Scholar
G. Ichim, J. Lopez, S. U. Ahmed et al., “Limited mitochondrial permeabilization causes DNA damage and genomic instability in the absence of cell death,” Molecular Cell, vol. 57, no. 5, pp. 860–872, 2015.
View at: Publisher Site | Google Scholar
S. Xu and H. R. Herschman, “A tumor agnostic therapeutic strategy for hexokinase 1-null/hexokinase 2-positive cancers,” Cancer Research, vol. 79, no. 23, pp. 5907–5914, 2019.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Fangfang Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

544

Downloads

883

Citations