Decreased brain-expressed X-linked 4 (BEX4) expression promotes growth of oral squamous cell carcinoma

Background Brain-expressed X-linked (BEX) 4 is a member of BEX family. The functional role of BEX4 in oral squamous cell carcinoma (OSCC) remains unknown. Methods Expression level of BEX family members (BEX1-5) in OSCC tissues and the paired normal epithelial were examined. Functions of epigenetic changes (DNA methylation and histone modifications) on BEX4 suppression in OSCC were examined by zebularine and trichostatin A (TSA) treatment on OSCC cell lines. Lentivector containing full-length BEX4 was used to generate OSCC cell lines with stable BEX4 expression. Effects of BEX4 expression on OSCC proliferation were monitored with xCELLigence RTCA real-time cell analyzer. BEX4-overexpressing CAL27 was implanted into nude mice to evaluate the effects on tumor growth in vivo. The signaling pathways regulated by BEX4 in OSCC was explored using human whole-transcript expression microarray. Results Among the 5 BEX family members, BEX1 and BEX4 showed significant down-regulation in OSCC (P < 0.001). BEX3, in comparison, was overexpressed in the primary tumor. BEX4 expression in OSCC cell lines was re-activated after zebularine and TSA treatment. High BEX4 expression could suppress proliferation of OSCC in vitro. Subcutaneous tumor volume of BEX4-overexpressing CAL27 was remarkably reduced in nude mice. Microarray experiment showed that S100A family members (S100A7, S100A7A, S100A8, S100A9 & S100A12) might be the downstream targets of BEX4 in OSCC. Conclusions BEX4 functions as tumor suppressor by inhibiting proliferation and growth of oral cancer. Decreased BEX4 contributes to the increased proliferative propensity of OSCC.


Background
Epithelial cancer derived from the oral mucosa is a common head and neck cancer [1]. Histologically, oral squamous cell carcinoma (OSCC) is the major subtype contributing to over 90 % of the cases [2]. Despite the advances in cancer detection techniques and treatment methods in the past decade, the incidence of OSCC remains high and the 5-year survival rate of OSCC patients remains unsatisfactory [3]. OSCC is an aggressive head and neck cancer due to the rapid proliferative capability and high invasive nature [4]. The survival rate of OSCC patients showed close association with the tumor size [5]. It has been shown that tumor size is an independent predictive factor for survival of OSCC patients [6]. Therefore, a better understanding on the underlying mechanisms involve in OSCC growth is important to improve the prognosis of OSCC patients.
Silencing key tumor suppressor genes by epigenetic mechanisms are essential for the carcinogenic development of OSCC. Identifying key methylated genes with functional implications in the development of OSCC would be useful in diagnosis and as therapeutic targets [7]. High-throughput parallel screening using expression and methylation microarrays had shown that brainexpressed X-linked (BEX) members were susceptible to epigenetic inactivation in OSCC [8]. Human BEX members (BEX1, BEX2, BEX3, BEX4, BEX5) are all located on chromosome Xq22. They are highly conserved genes among different species (including chimp, mouse, rat, dog and human) implying their importance in cellular functions [9]. At present, little is known about the functional implications of BEX family members in the progression of human malignancies. It has been reported that BEX family members could regulate cell cycle and apoptotic signaling [10,11].
BEX family members displayed distinct tissue distribution and have varying roles dependent on different cellular context. Thus, in the present study, we first explored the expression patterns of the human BEX family members in OSCC. Given that BEX4 expression was remarkably reduced in our cohort, we hypothesized that BEX4 could possibly play a tumor-suppressing role in OSCC. Further, whether DNA methylation or chromatin remodeling would be a possible cause of BEX4 suppression in OSCC will also be investigated. The functional role of BEX4 on tumor growth was examined using OSCC cell lines and xenograft models. As there was no information about the functional role of BEX4 in human malignancies, global gene expression change of BEX4-overexpressing OSCC cell line will be explored using microarray.

Cell cultures
Human tongue squamous cell carcinoma cell line CAL27 was obtained from American Type Culture Collection. Buccal squamous cell carcinoma cell line YD-38 was obtained from Korean cell line bank. Both cell lines were maintained in RPMI 1640 medium supplemented with 10 % fetal bovine serum, 100 U/ml penicillin and 100 μg/ ml streptomycin. Drug treatment started after cell seeding for 24 h. For zebularine treatment, the cells were harvested after 72 h. For trichostatin A (TSA) treatment, the cells were harvested after 24 h.

Clinical samples
Sixty-five OSCC and the paired normal tissues were examined. Tissue were obtained from patients underwent surgery in Division of Head and Neck, Queen Mary Hospital, Hong Kong. All the cancer tissues were squamous cell carcinoma. Consent of tissue donation for research purpose was obtained before sample collection. The study was approved by the ethnic review board of our institute (registered number: UW12-123).
RNA extraction and real-time quantitative polymerase chain reaction (qPCR) Total RNA was extracted using Trizol reagent (Invitrogen). First-strand cDNA synthesis was performed using High Capacity cDNA Reverse Transciption Kit (ABI). Transcript level was evaluated by qPCR using FastStart Universal Probe Master (Roche Applied Science) on LightCycler® 480 (Roche Applied Science). Table 1 listed the primers and probes used in this study. Reactions were performed at 95°C for 10 min followed by 45 cycles of 95°C for 15 s and 60°C for 1 min. Expression levels were evaluated using the comparative threshold cycle methods.

Immunocytochemistry
OSCC cells were seeded on chamber slides and treated with zebularine and TSA. Cells were fixed with 4 % paraformaldehyde and permeabilized using 0.2 % Triton X-100 in PBS. After blocking, cells were incubated with anti-BEX4 antibodies (Abcam) and CF488-conjugated secondary antibodies (Biotium). Cells were visualized by fluorescent microscope (Nikon).

Immunohistochemistry
Tissue samples embedded into paraffin blocks were cut into 4-μm sections. After deparaffinization and rehydration, sections were microwaved in 10 mM sodium citrate buffer for antigen retrieval and stained with anti-BEX4 antibodies (Abcam) at room temperature for 1 h. DAKO EnVision + System, HRP (DAKO) was used for visualization of immunoreaction. Then, sections were stained with Mayer's hematoxylin, dehydrated and photographed under light microscope. BEX4 protein expression level was graded as low, moderate or high depending on the staining intensity. Association between BEX4 expression and the clinicopathological parameters of OSCC patients was evaluated by Chi-square test or Fisher's exact test.

Lentivirus vector construction, lentivirus production and infection
Full-length BEX4 coding sequence was cloned into lentiviral expression vector pCDH (System Biosciences) to generate the BEX4-expressing vector (pCDH-BEX4 vector). Vector sequences were verified by direct sequencing. Virus packaging was performed by transient transfection into 293 T cells using Lipofectamine® 2000 Transfection Reagent (Invitrogen). OSCC cell lines CAL27 and YD-38 were transduced with medium containing lentivirus.

SiRNA transfection
BEX4 siRNA-1, BEX4 siRNA-2 and negative control siRNA (QIAGEN) were transfected into OSCC cells using HiPerFect transfection reagent (QIAGEN) according to the manufacturer's protocol. After 72 h, cells were collected and the efficiency of BEX4 silencing was determined by qPCR and immunocytochemistry.

Western blot
His-tag BEX4 protein was obtained from Novus Biologicals. Protein was extracted from cells using RIPA lysis buffer. Protein concentration was measured by BCA protein assay system (Pierce). Proteins separated by sodium dodecyl sulphate-poly-acrylamide gel electrophoresis were transferred to PVDF membranes and incubated with anti-BEX4 antibodies (Abcam). After incubation with horseradish peroxidase-conjugated goat anti-rabbit secondary antibody, hybridization signals were developed using ECL Plus Western Blotting Detection Reagents (Amersham Biosciences).

Real-time proliferation assay
Real-Time cell kinetic analyzer xCELLigence RTCA (ACEA Biosciences). was used to monitor the dynamic changes of cell proliferation. Data analysis was performed using RTCA Control Unit and the preinstalled RTCA software. For real-time proliferation assay, E-plate 16 was used. Cells were seeded directly onto E-plate. Changes in baseline impedance resulting from the increase of cell numbers were monitored by gold micro-electrodes located at the bottom of E-plate. The proportional changes in impedance were recorded continuously and expressed as cell index (CI). Change of CI with time was monitored continuously for 72 h.

Cell cycle analysis
Apoptosis was determined by propidium iodide (PI) staining. Cells were harvested and fixed with 70 % ethanol for 24 h at 4°C. The cells were stained with PI and RNase A in PBS for 30 min in the dark. Cell-cycle distribution was analyzed with a FACS flow cytometer (BD Biosciences).

Animal models
Athymic nu/nu mice (5 weeks old, weight range: 18-22 g) were used. CAL27 with pCDH-BEX4 vector or empty vector pCDH were injected subcutaneously into the right flank of mice. Four animals were used in each group. Tumor growth was measured daily with calipers in two dimensions. After 30 days, mice were euthanized. The animals were maintained under pathogen-free conditions with controlled temperature and humidity.

Microarray
Affymetrix whole-transcript expression microarray Human Exon 1.0 ST was used. All the RNA samples were tested with Agilent 2100 Bioanalyser (Agilent Technologies). All the samples RNA had RIN value over 8. All the tests were performed in the Centre for Genomic Sciences, the University of Hong Kong. Data analysis was evaluated using GeneSpring version 13 (Agilent Technologies). Genes with more than 1.5-fold difference and p value below 0.01 were selected for qPCR validation.

Statistical analysis
Data were expressed as mean ± standard deviation. Ttest or Chi square test were used to compare difference between groups as appropriate. Data analysis was performed using SPSS® for Windows version 14.0 (SPSS, Inc., Chicago, IL, USA). All the tests were two-sided. Pvalues <0.05 were considered as statistically significant. SurvExpress (Interface v2.0, Database Update: Sep 30, 2015) was used to perform survival analysis on oral cancer dataset GSE26549 as described [12]. The dataset contains global gene expression data from 86 oral leukoplakia tissue samples examined with Affymetrix Human Gene 1.0 ST Array [13].

Aberrant expression of BEX4 in oral cancers
To identify the key BEX family members involved in OSCC, we used qPCR to quantify BEX1-5 transcript levels in 65 primary OSCC tissues and compared with the paired normal epithelia. In comparison with the normal counterpart, expression of BEX1 (p < 0.001), BEX4 (p < 0.001) and BEX5 (p = 0.006) were significantly reduced in the OSCC tissues (Fig. 1a). In contrast, BEX3 expression showed significant elevation in OSCC (p = 0.026). For BEX2, no significant difference was observed between the OSCC and paired normal tissues (p = 0.457). The characteristics of the 65 patients examined were shown in Table 2.
Low BEX4 expression was associated with poor outcome in patients with oral pre-neoplastic lesions We analyzed the association between BEX family member expression and the outcome of 86 oral leukoplakia patients from GSE26549 using SurvExpress. Low BEX4 expression was associated with poor cancer-free survival of oral leukoplakia patients (Hazard ratio = 2.76, p = 0.0072) (Fig. 1b, c). The correlation of low BEX4 expression with poor cancer-free survival in leukoplakia patients was more significant than other BEX family members including BEX1, BEX2, NGFRAP1 (BEX3) and BEX5 (Fig. 1c).

Low BEX4 expression was associated with advanced tumor stage
To confirm the specificity of BEX4 antibody, Western blot was performed on the his-tagged recombinant human BEX4 protein. Signal of the immunoreaction increased in parallel with the concentration of recombinant BEX4 (Fig. 2a). In tissues obtained from OSCC patients, majority of normal tissues (96.9 %) showed moderate or high expression of BEX4. In the OSCC tissues, low expression of BEX4 was observed in 53.3 % tumor samples from patients with OSCC, while it was only found in 3.3 % paired normal tissues (Fig. 2b). Moreover, 75 % tumor tissues with T3-4 stage exhibited low BEX4 expression. In contrast, only 10 % tumor tissues with T1-2 stage displayed low BEX4 expression (Table 3). Low BEX4 level was associated with advanced stage tumor (p = 0.001, Table 3).
Representative images of BEX4 staining on oral tissues were shown in Fig. 2c.

BEX4 was epigenetically silenced in OSCC
Next, we sought to examine whether DNA methylation and histone acetylation are the major epigenetic alterations leading to BEX4 suppression in OSCC. DNA methylation was controlled by DNA methyltransferases that catalyze the addition of a methyl group to DNA. Zebularine inhibited DNA methyltransferases through the formation of a covalent complex with them [14]. Histone deacetylases (HDACs) catalyze the removal of acetyl groups of histone, leading to chromatin condensation and transcriptional suppression. TSA was a specific inhibitor of HDAC activity [15]. Following treatment with zebularine and TSA, dose-depending increase in BEX4 expression was observed in CAL27 cells (Fig. 3).
The protein expression level of BEX4 was also enhanced following zebularine and TSA treatment (Fig. 3).

BEX4 reduced the proliferation and tumor growth of OSCC
Having established that BEX4 expression was suppressed in OSCC, we hypothesized that BEX4 has a tumor suppressing function in OSCC. We generated BEX4overexpressing OSCC cell lines (Fig. 4a-c) and measured the change of proliferation rate. BEX4-overexpressing cells showed significant reduction in OSCC proliferation (Fig. 4d, e). Apoptotic cells were determined by PI staining using flow cytometry. Flow cytometry analysis showed that BEX4 overexpression in OSCC did not induce apoptosis of OSCC (Fig. 4f). To confirm our results, two different BEX4-specific siRNAs were used to silence the expression of BEX4 in OSCC cells (Fig. 5a, b). Silence of BEX4 resulted in a significant increase in the proliferation of CAL27 and YD-38 cells (Fig. 5c, d). In the animal model, BEX4-overexpressing xenograft developed from CAL27 cells had a significantly smaller tumor volume and slower growth rate (Fig. 6). Tumor volume of control xenograft increased by 4.7-fold from day 17 to day 53. In contrast, tumor volume of BEX4overexpressing xenograft increased by only 1.7-fold (Fig. 6). The data is consistent with the cell line models suggesting that BEX4 was involved in controlling OSCC growth.

Discussion
Our results indicated that BEX4 functions as tumor suppressor in controlling proliferation and growth of OSCC. At present, little is known about the expression status and functions of brain expressed X-linked family in OSCC. Thus, we first examined the mRNA levels of BEX1, BEX2, BEX3, BEX4 and BEX5 in our OSCC tissues and compared with the paired normal epithelia. We identified that BEX4 suppression is a common event in our OSCC and is potentially associated with the poor outcome of patients. The association between low BEX4 with poor outcome of patients with oral pre-neoplastic lesions suggested that reduced BEX4 is an early event in the development of OSCC. OSCC is a male dominant disease. Recent epidemiology study revealed that there is an increasing number of female OSCC cases [16]. In our cohort, reduced   b Immunostaining analysis showed that BEX4 protein increased in CAL27 and YD-38 after zebularine and TSA treatment Members of BEX family have been implicated in regulating apoptosis in various human cancer cells and normal cells. In breast cancer, BEX2 could modulate ceramide-induced apoptosis via protein phosphatase 2A [17,18]. BEX1 interacted with B-cell lymphoma 2 (BCL-2) and inhibited the formation of BCL-2/BCL-2-associated X protein (BAX) complex, leading to enhancement of imatinib-induced apoptosis in human leukemic cell line K562 [19]. BEX3 could interact with the death domain of p75 neurotrophin receptor (p75NTR) and induce apoptosis upon nerve growth factor (NGF) treatment in HEK293 cells [20]. In light of the regulating role of BEX family in apoptosis, we investigated whether BEX4 At present, the functional roles of BEX family members in cancer remain unclear. In OSCC, it has been reported that BEX1 is involved in modulating nuclear factor-κB (NF-κB) signaling pathway [8]. BEX2 expression could possibly be modulated by the aberrantly activated mechanistic target of rapamycin (mTOR) in cancer cells [21]. BEX protein family is also known as intrinsically disordered proteins (IDPs), which exist as dynamic ensemble of conformation with no discrete 3D structure. IDPs are suggested to play an important role in tumorigenesis and are involved in modulating multiple gene regulating signaling pathways by their ability to recognize multiple interaction partners [22]. Our results revealed that BEX family member, BEX4, is one of the IDP involved in the pathogenesis of OSCC.
Epigenetic silencing of tumor suppressor gene is a characteristic feature of OSCC. Existing data suggests that DNA methylation and chromatin remodeling play an important role in modulating BEX family expression in human malignancies [23,24]. In glioma, BEX1 and BEX2 are silenced and the expression could be activated by treatment with DNA methyltransferase inhibitor and histone deacetylase inhibitor [25]. Our results are consistent with the current knowledge suggesting that BEX4 is epigenetically silenced in OSCC.
OSCC overexpressing BEX4 showed significant upregulation of S100A family members (S100A7, S100A8, S100A9, S100A12, S100A15 or S100A7A). S100A is a family of calcium binding protein involved in multiple processes of cancer development [26]. Of which, S100A7 dysregulation is believed to be involved in a wide range of pathological processes associated with OSCC. S100A7 (also known as psoriasin) was expressed in early stage oral cancer. S100A7 expression was reduced in tumor in late stage [27]. Additionally, highly differentiated OSCC usually had high expression of S100A7 as compared to the moderate or poorly differentiated counterparts [28]. Forced expression of S100A7 is accompanied by repressed cell proliferation and tumor progression in orthotopic tongue tumor model [27]. Apart from S100A7, aberrant expression of other S100A family members is emerging as an important event in head and neck cancers. In head and neck squamous cell carcinoma, S100A8/A9 (calprotectin) could suppress cell growth by inducing G2/M cell cycle checkpoint arrest [29]. The mRNA and protein expression levels of S100A8/ A9 are decreased in head and neck cancers. S100A8/A9 can modulate G2/M Cdc2/cyclin B1 complex activation resulting in G2/M checkpoint arrest [29]. For S100A12, high level in tumor tissue is a good prognostic indicator for oropharyngeal squamous cell carcinoma patient [30]. S100A7A (koebnerisin) are highly homologous to S100A7. Both S100A7A and S100A7 could serve as oncogenic inactivator in human malignancies [31].
The molecular mechanisms underlying S100A proteins dysregulation remains unclear. In OSCC, it has been shown that overexpression of β-catenin decreased the expression of S100A7, while silence of β-catenin enhanced its expression, indicating that β-catenin signaling negatively modulated S100A7 [27]. In addition, c-Myc suppressed S100A7 expression by directly binding to its promoter [32].

Conclusions
Based on the current results, we suggested that BEX4 controls OSCC proliferation. Although the mechanisms underlying BEX4 suppression in OSCC remains unresolved, induced expression of BEX4 in the presence of epigenetic drugs suggested that BEX4 is epigenetically silenced in OSCC. Taken together, we suggested that BEX4 functions as a novel tumor suppressor gene involved in OSCC.