Overexpression of HDAC9 is associated with poor prognosis and tumor progression of breast cancer in Chinese females

Abstract

Background

Breast cancer represents a serious health issue among females. HDAC9 has been identified as an oncogene in human cancers. This study sought to assess the prognostic value and the biologic function of HDAC9 in breast cancer patients.

Methods

Expression of HDAC9 in breast cancer tissues and cells was evaluated by quantitative real-time polymerase chain reaction. Kaplan–Meier survival analysis and Cox regression assay were conducted to explore the prognostic significance of HDAC9. Cell experiments were performed to investigate the effects of HDAC9 on the biologic behaviors of breast cancer cells.

Results

Expression of HDAC9 was significantly upregulated in both cancerous tissues and cells compared with the normal controls (all P<0.05). Overexpression of HDAC9 was correlated with lymph node metastasis (P=0.021) and TNM stage (P=0.004). Patients with high HDAC9 had poor overall survival compared to those with low levels of HDAC9 (log-rank P<0.05). Elevated HDAC9 was found to be an independent prognostic factor for the patients (hazard ratio=2.996, 95% CI=1.611–5.572, P=0.001). According to the cell experiments, tumor cell proliferation, migration and invasion were suppressed by knockdown of HDAC9.

Conclusion

All data demonstrated that overexpression of HDAC9 serves as a prognostic biomarker and may be involved in the tumor progression of breast cancer.

Keywords:

• HDAC9
• prognosis
• progression
• breast cancer

Introduction

Breast cancer is the most common malignancy among females around the world.¹ Patients suffering from breast cancer usually have the symptoms of change in breast shape, fluid coming from the nipple, dimpling of skin or a red scaly patch of skin.^2,3 Every year, ~1,300,000 cases are diagnosed with breast cancer and about 465,000 deaths are estimated to occur worldwide.⁴ Some risk factors have been identified to be correlated with the occurrence of breast cancer, including obesity, drinking alcohol, lack of physical exercise, ionizing radiation, early age at first menstruation, family history and older age.⁵ Despite advances in surgery, chemotherapy and radiotherapy, the prognosis of breast cancer remains dismal.⁶ Data from previous studies revealed that the prognosis could be achieved by using some molecular biomarkers.⁷ It is, therefore, necessary to uncover more precise prognostic biomarkers for determining prognosis in patients with breast cancer.

It is generally considered that histone acetyltransferase and histone deacetylase (HDAC) play important roles during the regulation of gene transcription.^8,9 Currently, the relationship between HDACs and progression of cancer has attracted attention in different malignancies. Data in several previous studies reveal that an aberrant expression of HDACs has been detected in many tumor samples.^10,11 In addition, HDACs have been found to inhibit tumor suppressor expression by binding to the promoter region.^12,13 HDAC9, a subtype of HDAC, has been investigated in some types of human cancers. For example, Moreno et al demonstrated that the expression of HDAC9 was deregulated and it correlated with overall survival of patients with lymphoblastic leukemia.¹⁴ However, reports about the role of HDAC9 in breast cancer are currently limited.

To better understand the relationship of HDAC9 and breast cancer, this study examined the expression patterns and prognostic significance of HDAC9 in breast cancer patients. The effects of HDAC9 on biologic behaviors of cancer cells were also assessed.

Materials and methods

Patients and tissue sample collection

Tissue specimens used for the subsequent experiments were collected from 118 breast cancer patients who underwent surgery between 2007 and 2011 at the hospital and were verified by experienced pathologists. None of these patients had received any antitumor therapy before the sampling. The breast cancer tissues and the adjacent normal tissues were snap-frozen in liquid nitrogen after collection. Signed informed consent was obtained from each patient, and this study was approved by the Ethics Committee of Tenth People’s Hospital of Tongji University. Moreover, the clinicopathologic information on age, tumor size, estrogen receptor status, progesterone receptor status, human epidermal growth factor receptor 2 status, lymph node metastasis and TNM stage was recorded from the electronic medical records of the patients and are summarized in Table 1. After surgery, all the patients were enrolled in a 5-year follow-up survey. The survival information was obtained for the subsequent survival analysis.

Table 1 Association of HDAC9 expression with the clinical features of breast cancer patients

Features	Total (N=118)	HDAC9 expression		P-value
		Low (n=54)	High (n=64)
Age (years)				0.763
≤50	42	20	22
>50	76	34	42
Tumor size (cm)				0.360
≤3	45	23	22
>3	73	31	42
ER status				0.264
Negative	46	24	22
Positive	72	30	42
PR status				0.129
Negative	48	26	22
Positive	70	28	42
HER2 status				0.360
Negative	45	23	22
Positive	73	31	42
LN metastasis				0.021
Negative	52	30	22
Positive	66	24	42
TNM stage				0.004
I–II	55	33	22
III–IV	63	21	42

Abbreviations: ER, estrogen receptor; HER2, human epidermal growth factor receptor 2; LN, lymph node; PR, progesterone receptor.

Cell lines and transfection

Human breast cancer cell lines MCF-7, BT474 and normal human breast epithelial cell line MCF-10A were purchased from the American Type Culture Collection (ATCC, Manassas, VA, USA). All these cells were cultured in Roswell Park Memorial Institute (RPMI)-1640 medium supplemented with 10% fetal bovine serum and kept in a humidified incubator with 5% CO₂ at cell transfection was conducted by using Lipofectamine 2000 Reagent (Thermo Fisher Scientific, Waltham, MA, USA) following the instructions of the manufacturer. Cancer cells transfected with HDAC9 siRNA were defined as the experimental group and the other cells treated with control siRNA or Lipofectamine 2000 were used as the control group.

HDAC9 inhibition

In addition to siRNA, vorinostat (suberoylanilide hydroxamic acid [SAHA]) was also used to suppress the expression of HDAC9 in breast cancer cells. MCF-7 and BT474 were seeded in a 96-well culture plate and treated with 10 μM SAHA (vorinostat; MedChem Express Co., Monmouth Junction, NJ, USA) or vehicle (dimethyl sulfoxide [DMSO], 1:1,000) for 48 h. The drug and the medium were replenished every 24 h.

RNA extraction and quantitative real-time polymerase chain reaction (qRT-PCR)

Total RNA was extracted from the tissues and cells using TRIzol reagent (Thermo Fisher Scientific) as per the manufacturer’s instructions. Pure RNA was obtained by calculating the ratio of OD A260/280, which was used to indicate the purity of RNA. Reverse transcription was performed to synthesize cDNA from RNA by Transcriptor First Strand cDNA Synthesis Kit (Roche, Vilvoorde, Belgium). Expression of HDAC9 was investigated by using qRT-PCR, which was carried out with the SYBR green I Master Mix kit (Thermo Fisher Scientific) and was run on the 7300 Real-Time PCR System (Thermo Fisher Scientific). The primer sequences of HDAC9 were as follows: forward: 5′-AACTGGAGCAGCAGAGGCAAG-3′, reverse: 5′-TACTTCTGTACTTGCCACTGCC-3′. Besides, GAPDH was used as the internal control with the following primers: forward: 5′-GGCCTCCAAGGAGTAAGACC-3′, reverse: 5′-AGGGGTCTACATGGCAACTG-3′. The final relative expression of HDAC9 was calculated with 2^−ΔΔCt method and normalized to GAPDH.

Cell proliferation assay

In order to examine the effects of HDAC9 on cell proliferation of breast cancer cells, the colorimetric MTT analysis was carried out in this study. MCF-7 and BT474 cells were seeded in two 96-well culture plates. Cells in one plate were transfected with HDAC9 siRNA or control vectors (control siRNA or Lipofectamine 2000). Meanwhile, the cells in another plate were treated with SAHA or DMSO for 48 h. After the treatment, each well was added with 10 μL MTT (5 mg/mL; Sigma-Aldrich) and incubated at 37°C for 4 h. Then, 100 μL DMSO (Sigma-Aldrich) was added to the wells to dissolve the formazan crystals. The absorbance value was measured at 490 nm with a spectrophotometer (Multiskan MK3; Thermo Fisher Scientific). Experiments were repeated in triplicate.

Cell migration and invasion analysis

To uncover the effects of HDAC9 on cell migration and invasion, the Transwell analysis was carried out with a 24-well Transwell chamber. The cells transfected with siRNA or treated with SAHA were added in the upper compartment with a concentration of 1×10⁵ per well and then incubated in serum-free RPMI-1640 medium at 37°C for 24 h. The lower compartment contained 300 μL RPMI-1640 medium supplemented with 20% fetal bovine serum, which was used as the chemotactic factor. After incubation for 24 h, the cells that migrated to the lower compartment were stained with 0.1% crystal violet and counted by a microscope. For the invasion assay, the upper chamber was coated with Matrigel (BD, Bedford, MA, USA).

Statistical analysis

All the statistical analyses were conducted by SPSS software (SPSS Inc., Chicago, IL, USA), and the data used were expressed as mean±SD. Differences between the two groups were examined by Student’s t-test. Association of HDAC9 with clinicopathologic features was assessed with chi-square test. Survival analysis was conducted using Kaplan–Meier method and log-rank test. Cox regression analysis was adopted to confirm the prognostic performance of HDAC9 for breast cancer patients. Statistical differences with P<0.05 were considered as statistically significant.

Results

HDAC9 expression in tissue specimens and cells

In this study, the expression of HDAC9 in the tissues and cell lines was estimated by qRT-PCR. The analysis results showed that HDAC9 expression was significantly higher in breast cancer tissues than that in the paired normal tissues (P<0.001, Figure 1A). To confirm this result, HDAC9 expression was also evaluated in breast cancer cells, and it was found that its expression was upregulated in breast cancer cells compared with the normal cells (all P<0.05, Figure 1B), which was in accordance with the results in tissue samples.

Figure 1 The mRNA expression of HDAC9 measured by qRT-PCR.

Notes: (A) HDAC9 expression in breast cancer tissues and paired normal tissues. The expression of HDAC9 was higher in breast cancer tissues than that in the matched normal tissues (***P<0.001). (B) HDAC9 expression in cells. The expression of HDAC9 was increased in breast cancer cells compared with the normal cells (**P<0.01, ***P<0.001).
Abbreviation: qRT-PCR, quantitative real-time polymerase chain reaction.

Association of HDAC9 with the clinical characteristics of breast cancer patients

Chi-square test was used to estimate the relationship between HDAC9 expression and cancer patients’ clinicopathologic data. To facilitate this assay, a cutoff value of mean HDAC9 expression was chosen to be used to classify the patients into low HDAC9 expression group (n=54) and high HDAC9 expression group (n=64). All the analysis results are detailed in Table 1 and reveal that the expression of HDAC9 was influenced by lymph node metastasis (P=0.021) and TNM stage (P=0.004). However, no relationship was found between HDAC9 expression and age, tumor size, estrogen receptor status, progesterone receptor status or human epidermal growth factor receptor 2 status (all P>0.05).

Prognostic significance of HDAC9 for breast cancer

In this study, the relationship between HDAC9 and overall survival of breast cancer patients was also investigated. Based on the survival information obtained from the 5-year follow-up survey, the Kaplan–Meier survival analysis was performed for breast cancer patients. Survival curves shown in Figure 2 reveal that patients with high HDAC9 expression levels had shorter survival time compared to those with low HDAC9 levels (log-rank P<0.05). Furthermore, Cox regression analysis was carried out to examine the factors that might influence the overall survival. The univariate and multivariate Cox analyses results detailed in Table 2 demonstrate that upregulated HDAC9 was closely correlated with poor overall survival and could be used as an independent prognostic factor for patients with breast cancer (hazard ratio=2.996, 95% CI=1.611–5.572, P=0.001).

Table 2 Univariate and multivariate Cox analysis for HDAC9 in breast cancer patients

Characteristics	Univariate analysis			Multivariate analysis
	HR	95% CI	P-value	HR	95% CI	P-value
HDAC9	2.996	1.611–5.572	0.001	2.996	1.611–5.572	0.001
Age	1.117	0.646–1.930	0.692	–	–	–
Tumor size	1.243	0.726–2.129	0.428	–	–	–
ER status	0.996	0.582–1.704	0.989	–	–	–
PR status	0.916	0.537–1.562	0.748	–	–	–
HER2 status	1.291	0.744–2.240	0.363	–	–	–
Lymph node metastasis	1.236	0.710–2.151	0.454	–	–	–
TNM stage	0.870	0.509–1.489	0.612	–	–	–

Abbreviations: ER, estrogen receptor; HER2, human epidermal growth factor receptor 2; HR, hazard ratio; PR, progesterone receptor.

Figure 2 Kaplan–Meier survival curves for breast cancer patients based on the expression of HDAC9.

Notes: P-value was calculated using the log-rank test. Patients with high HDAC9 expression had shorter survival time than those with low HDAC9 expression (log-rank P<0.05).

HDAC9 reduction inhibits proliferation, migration and invasion of breast cancer cells

To investigate the functional role of HDAC9 in PCa, its effects on tumor cell proliferation, migration and invasion were examined in this study. Two breast cancer cell lines MCF-7 and BT474 were transfected with HDAC9 siRNA or treated with SAHA to suppress the expression of HDAC9. The results of qRT-PCR showed that the expression of HDAC9 in breast cancer cells transfected with HDAC9 siRNA was significantly lower than that in the cells with control vectors (P<0.01, Figure 3). Similarly, markedly reduced HDAC9 was also detected in the cells treated with SAHA compared to that in those with DMSO (P<0.01, Figure 3). These data indicate that the HDAC9 expression was successfully reduced by siRNA and SAHA. MTT assay was adopted to analyze the cell proliferation, which showed that breast cancer cell proliferation was suppressed in HDAC9-knockdown cells compared with the control cells (P<0.05, Figure 4A and B). In addition to proliferation, cell migration and invasion were also assessed using Transwell analysis. According to the analysis results, the tumor cell migration and invasion were both found to be inhibited in the cells with reduction of HDAC9 compared with the controls (P<0.05, Figure 4C–F).

Figure 3 Expression of HDAC9 in cells treated with siRNA or SAHA.

Notes: (A) HDAC9 expression was significantly decreased in MCF-7 transfected with HDAC9 siRNA and MCF-7 treated with SAHA compared with the controls (**P<0.01, ***P<0.001). (B) HDAC9 expression was significantly decreased in BT474 transfected with HDAC9 siRNA and BT474 treated with SAHA compared with the controls (**P<0.01).
Abbreviations: DMSO, dimethyl sulfoxide; SAHA, suberoylanilide hydroxamic acid.

Figure 4 Effects of HDAC9 expression on cell proliferation, migration and invasion of breast cancer cells.

Notes: (A, B) Cell proliferation was suppressed by HDAC9 knockdown in MCF-7 and BT474 cells (*P<0.05, **P<0.01). (C, D) Number of migrated cells was decreased by a reduction of HDAC9 in breast cancer cells (**P<0.01, ***P<0.001). (E, F) Cell invasion was inhibited by the downregulation of HDAC9 in both MCF-7 and BT474 cells (**P<0.01).
Abbreviations: DMSO, dimethyl sulfoxide; SAHA, suberoylanilide hydroxamic acid.

Discussion

As the most prevalent cancer occurring in women, breast cancer has received lots of attention on its progression and treatment.¹⁵ It can be diagnosed at all age groups and has been found to be a great threat to healthy life.¹⁶ So far, great progress has been made in therapeutic methods, and the mortality of breast cancer has been reduced.¹⁷ However, the prognosis of some cancer cases is not satisfactory mainly due to the advanced stage of the tumors.¹⁸ Therefore, prognosis is urgently needed to be improved for breast cancer patients. Data from recent studies show that cancer prognosis has improved due to using related biomarkers, which play crucial roles during tumor progression.^19,20 In breast cancer, some prognostic biomarkers have also been identified. For instance, Fu et al demonstrated that SOX17 expression was downregulated and correlated with the poor prognosis of breast cancer.²¹ MicroRNA-106b, another example, has been proved to be involved in the recurrence of breast cancer and to act as an independent prognostic factor in patients with breast cancer.²² Jerzak et al reported that the thyroid hormone receptor α (THRα) represented an efficient prognostic biomarker in breast cancer patients.²³ All these data indicated the pivotal role of cancer-related molecules for breast cancer prognosis. Consequently, more molecular biomarkers should be identified for prediction of prognosis in breast cancer.

During the development of breast cancer, numerous epigenetic changes take place, such as methylation of DNA and diverse histone modifications, including phosphorylation, methylation, sumoylation, ubiquitination and acetylation.²⁴ HDACs represent a series of proteins which play pivotal roles in acetylation.²⁵ Eighteen members of the HDAC family have been identified in human beings. Recent studies have demonstrated that HDACs can be used as therapeutic targets and are involved in tumor progression in different cancers, including breast cancer.²⁶ The altered expression of HDAC1, HDAC5, HDAC6 and HDAC8 has been identified as a reliable prognostic biomarker for breast cancer patients.^27–30 HDAC2 and HDAC3 were found to be associated with the aggressive behavior of breast cancer.³¹ HDAC4 has been reported to mediate the antitumor effects of microRNA-125a-5p in breast cancer.³² HDAC9, also a member of HDAC family, has been found to be deregulated in some tumor samples, such as medulloblastomas and lung cancer.^33,34 Moreover, Lapierre et al focused on the functional role of HDAC9 in proliferation of breast cancer cells and found that overexpression of HDAC9 could promote cell proliferation.³⁵ However, the clinical significance of HDAC9 has been rarely reported in patients with breast cancer. To better understand the role of HDAC9, its expression patterns and prognostic significance were assessed in patients with breast cancer.

In this study, the expression of HDAC9 was measured in the breast cancer tissues and cells using qRT-PCR. The analysis results showed that HDAC9 expression was remarkably upregulated in breast cancer tissues and cell lines compared with the normal controls. Moreover, the overexpression of HDAC9 was found to be correlated with lymph node metastasis and TNM stage. Thus, we considered that HDAC9 might act as an oncogene in breast cancer and is involved in the development of this malignancy. Additionally, the clinical significance of HDAC9 in breast cancer prognosis was also investigated. From the Kaplan–Meier survival curves, we found that patients with high HDAC9 expression levels had shorter survival time than those with low HDAC9 levels, suggesting the increased HDAC9 was correlated with poor overall survival. Data in the Cox regression analysis revealed that the upregulated HDAC9 was an independent prognostic factor for patients with breast cancer.

To further confirm the functional role of HDAC9 during breast cancer progression, its effects on biologic behaviors of breast cancer cells were also examined in the current study. In the breast cancer cells, siRNA and SAHA were separately used to silence the expression of HADC9. The analysis results demonstrated that cell proliferation, migration and invasion were suppressed by the knockdown of HDAC9. All the above data suggested that HDAC9 reduction could inhibit tumor progression in breast cancer. In a previous study performed by Lapierre et al,³⁵ overexpression of HDAC9 was reported to promote cell proliferation of breast cancer cells, and it might exert its antitumor effects by targeting sex-determining region Y-box 9 protein (SOX9), which has been described as a crucial molecule during tumor progression of breast cancer by regulating Wnt/β-catenin pathway.³⁶ However, the precise molecular mechanisms underlying the role of HDAC9 in breast cancer need to be confirmed in further studies. The results of this study might be limited by the small sample size, and therefore, further studies are necessary with larger research cohort.

Conclusion

In summary, all data in this study revealed that overexpression of HDAC9 was correlated with the progression of breast cancer and it could be used as a candidate prognostic biomarker.

Disclosure

The authors report no conflicts of interest in this work.