Background {#Sec1}
==========

Gliomas are the most common central nervous system tumor in adults. Glioblastoma, also known as glioblastoma multiform (GBM), is the most aggressive, malignant, and lethal glioma. Despite radiation and temozolomide therapy, the prognosis of patients with GBM remains extremely poor, with a median survival time of only 14.5 months from diagnosis to death \[[@CR1]-[@CR3]\]. Recently, DNA methylation alterations have been widely reported in human gliomas, including global hypomethylation and promoter-associated CpG island hypermethylation \[[@CR4]\]. To better understand the molecular mechanisms and pathological effects of genomic methylation abnormality on GBM, we performed genome-wide DNA methylation analysis on GBM (n =  42) and normal brain tissues (n  = 8). We found that the *testis derived transcript* (*TES*) gene promoter was broadly hypermethylated in GBM patients.

The *TES* gene encodes a protein (Tes) of 421 amino acids which plays a role in protein--protein interactions \[[@CR5],[@CR6]\]. TES has an important role in cell adhesions and largely affects cell motility in a vast number of cancers \[[@CR7]\]. According to recent studies, TES is widely considered as a putative tumor suppressor gene in malignances such as head and neck squamous cell carcinomas \[[@CR8]\], ovarian cancer \[[@CR9]\], primary gastric cancer \[[@CR10]\], and prostate cancer \[[@CR11]\]. However, TES methylation status and pathological function in GBM still remain largely unclear \[[@CR12],[@CR13]\].

Herein, we examined the expression of TES in human GBMs and investigated the biological role of TES in the pathogenesis of this tumor. We found that TES was down-regulated in both glioblastoma cell lines and GBM tissues. TES expression level correlated inversely with the clinical outcome. Our results also indicate that down-regulation of TES protein expression is anti-apoptotic to GBM cell line.

Methods {#Sec2}
=======

Tissue samples {#Sec3}
--------------

All samples were collected from the Department of Neurosurgery, Beijing Tiantan Hospital, between June 2009 and June 2012. Tissues were frozen in liquid nitrogen immediately after surgery; some were stored at −80°C following the extraction of RNA, whereas others were formalin-fixed and paraffin-embedded for immunohistochemistry. Informed consent from patients and ethics approval from the ethics committee was achieved. This study have been approved by the Review board of Beijing Tiantan Hospital, Capital Medical University.

DNA methylation profiling {#Sec4}
-------------------------

Methylation profiling was performed using the Illumina Infinium Methylation27 Bead Array. The BeadChip contains 27,578 highly informative CpG sites covering more than 14,000 human RefSeq genes. This allows searches to interrogate all these sites per sample at a single nucleotide resolution. Bisulfite modification of DNA, chip processing, and data analysis were carried out by following the manufacturer's manual. DNA was extracted from peripheral-blood buffy coats using the QIAmp DNA mini kit according to the manufacturer's protocol (Qiagen, Valencia, CA, USA); the quality of DNA samples was assessed by electrophoresis in a 1% agarose gel. DNA were then bisulfite-modified using the EZ DNA methylation kit (Zymo Research) and hybridized according to the manufacturer's instructions.

Immunohistochemistry {#Sec5}
--------------------

For immunohistochemistry, we used a standard avidin-biotin-peroxidase complex method. Briefly, 4-μm thick tissue sections were dewaxed and the endogenous peroxidase activity was blocked with 3% hydrogen peroxide in deionized water. For antigen retrieval, sections were immersed in an antigen repair liquid kit and then rinsed in phosphate-buffered saline (PBS; pH 7.2). The sections were incubated at 37°C for 1 h with anti-TES monoclonal antibody. After rinsing in PBS, sections were incubated with a goat anti-mouse secondary antibody for 30 min, and washed in PBS. A DAB coloration kit (ZSGB-BIO ORIGENE, ZLI-9018) was then used; PBS was used as a negative control. The percentage of positive tumor cells was counted under a bright-field Olympus microscope. Sections with fewer than 25% labeled cells indicated low expression of TES. Sections with labeling of ≥25% indicated high expression of TES.

Cell lines and treatment with 5-Aza-2′-deoxycytidine (5-aza-dC) treatment {#Sec6}
-------------------------------------------------------------------------

Glioblastoma cell lines (U251) were routinely cultured in Dulbecco's modified eagle medium (Hyclone, Thermo Scientific, USA). The medium was supplemented with 1% L-glutamine (Hyclone), 10% fetal bovine serum (Gibco, Invitrogen Corporation, Carlsbad, CA, USA), 1% penicillin streptomycin (Hyclone), and grown at 37°C in a humidified atmosphere of 5% CO~2~/95% air. 5-aza-dC (Sigma Aldrich Corporation, St Louis, MO, USA) was dissolved in PBS (Hyclone) and filtered with a 0.22 μm filter membrane.

Western blotting {#Sec7}
----------------

U251 cells were lysed in 1% NP40 buffer (150 mM, 1% NP40, 50 Mm Tris, pH 8.0) following treatment with 5-aza-dC; 30 μg of protein samples were loaded onto 10% SDS-PAGE gels and then transferred to polyvinylidene fluoride membranes (Millipore). The membranes were then blocked with skimmed milk in Tris-buffered saline Tween-20 for 2 h. The following primary antibodies were used: anti-TES (Abcam, Cambridge, MA, USA) and anti-b-tubulin (Sigma-Aldrich), followed by reaction with a goat anti-mouse second antibody.

Flow cytometry analysis of cell apoptosis {#Sec8}
-----------------------------------------

Cell apoptosis was analyzed by an Annexin V-FITC⁄propidium iodide (PI) kit according to the manufacturer's instructions. Briefly, harvested cells were resuspended in 100 μL Annexin V-FITC binding buffer and adjusted to about 1 × 10^6^ ⁄mL, then 5 μL Annexin V-FITC and 5 μL PI (20 mg ⁄mL) were added and incubated for 15 min at room temperature in the dark. Flow cytometry was conducted on a fluorescence-activated cell sorter (FACSCalibur™ BD Biosciences, Hercules, CA, USA). Each experiment was performed in triplicate.

Statistical analysis {#Sec9}
--------------------

Cox analysis was performed using Matlab2009. Survival curves were calculated according to the Kaplan-Meier method and differences between curves were assessed using the log-rank test. The difference between high and low TES-expressing tissues was assessed using the χ^2^ test. Significant differences among groups were determined using the Student's *t*-test. *P* \<0.05 was considered to be significant. All statistical analyses were carried out using the software GraphPad Prism (GraphPad Software, La Jolla, CA, USA) and SPSS version 16.0 (SPSS, Chicago, IL, USA).

Results and discussion {#Sec10}
======================

Identification of promoter hypermethylation of *TES* gene in GBM {#Sec11}
----------------------------------------------------------------

We performed genome-wide DNA methylation analysis on 42 GBM and 8 normal brain tissues. To identify genes that had different methylation status between pGBM and normal tissues, all 50 methylation data were subjected to Significant Analysis of Microarray and 42 GBM were analyzed by Survival Analysis (Cox analysis) Metlab2009. Overall, 21 differential genes displayed correlation between promoter methylation and the overall survival in glioblastomas: *TUBA2*, *CKMT2*, *AMT*, *HOXC11*, *TRIM58*, *LCE1B*, *KLF14*, *SERPINB12*, *CYP2A7*, *TREML2*, *UNQ467*, *SPRR3*, *TES*, *RGN*, *DEFA1*, *AJAP1*, *PDE4C*, *SPRR2D*, *WNT7B*, *UBQLN3*, and *LCE1F*. Independent *t*-test and Kaplan-Meier method were performed on these genes to identify the differentiation and prognostic value in a univariate manner. Finally, eight genes were identified: *WNT7B*, *CKMT2*, *AMT*, *UNQ467*, *DEFA1*, *TUBA2*, *TES*, and *UBQLN3*. We selected the *TES* gene as the best candidate gene for further study (Figure [1](#Fig1){ref-type="fig"}A). To verify the hypermethylation of *TES* in GBM, treatment with demethylation reagent 5-aza-dC on the GBM cell line U251 (low TES expression) can greatly increase TES expression (Figure [1](#Fig1){ref-type="fig"}B). This result confirms the microarray data that *TES* is hypermethylated in GBM.Figure 1***TES*** **promoter methylation status in glioblastoma (GBM, n = 42) and normal (n = 8) samples analyzed by DNA promoter methylation microarray. (A)** Comparison of *TES* promoter methylation status in GBM and normal samples by microarray analysis (up) and independent *t*-test analysis (down). **(B)** Western blot of TES in U251 cell line before and after 5-aza-2-deoxycytidine treatment. Anti-tublin was used as a protein loading control.

Down-regulation of TES protein expression is in GBM tissues {#Sec12}
-----------------------------------------------------------

TES protein expression in GBM was surveyed by immunohistochemistry using anti-TES monoclonal antibody on 37 independent GBM and 10 non-tumor samples. A reduced TES protein level was found in 33 (89.2%) tumor tissues and in 2 (20%) non-tumor tissues (*P* \<0.001, Figure [2](#Fig2){ref-type="fig"}A). For evaluation of TES protein level, classification standards were as follows: high, \>25% tumors cells were positively stained, low, no staining or \<25% tumor cells showed negatively staining. Collectively, these results demonstrated a significant down regulation of TES protein expression in GBM (*P* \<0.001, Figure [2](#Fig2){ref-type="fig"}B).Figure 2**TES expression in GBM. (A)** Immunohistochemical staining of TES protein expression in normal tissue (left) and a GBM specimen (right) (magnification × 400). **(B)** Collective summary of TES staining intensity from all samples.

TES enhances apoptosis in glioblastoma cells {#Sec13}
--------------------------------------------

To evaluate pathogenic effect of low expression of TES in GBM, apoptotic effect was assessed on 5 μm 5-aza-dC treated and untreated U251 cells. The percentage of apoptotic cells was calculated by flow cytometry using the Annexin-V/PI double-staining assay (Figure [3](#Fig3){ref-type="fig"}A). In the untreated control group (no 5-aza-dC was given), the total rate of apoptosis was 8.75 (Figure [3](#Fig3){ref-type="fig"}B). However, in U251 cells treated with 5-aza-dC, many more apoptotic cells were detected: the rate raised to 12.25, suggesting that TES promotes apoptosis of glioblastoma cells.Figure 3**Cell apoptosis induced by treatment of 5-aza-2-deoxycytidine.** The percentage of apoptotic cells was calculated by flow cytometry using the Annexin-V/PI double-staining assay. Control vs. 5AZA (\**P* \<0.05). **(A)** (Control) Apoptosis of U251 cells in the group without 5-aza-dC treatment (8.75); **(B)** (5AZA) Apoptosis of U251 cells in the 5 mL 5-aza-2-dC group (15.25).

TES protein expression is highly correlated with worse outcome of GBM patients {#Sec14}
------------------------------------------------------------------------------

After TES promoter hypermethylation was identified in some GBM patients, we further estimated the association of this abnormality with GBM patients' survival using Kaplan-Meier plot, as shown in Figure [4](#Fig4){ref-type="fig"}A. Results showed that the overall survival (OS) of patients with TES hypermethylation is significantly shorter than those with hypomethylation (median survival, 236.5 vs. 565 days *P* \<0.05, log-rank test). Patient outcome was also evaluated independently from TES protein expression level from results of immunohistochemistry analysis. Patients with high TES expression had a median OS of 221 days, while those with low expression had an OS of 345 days (median survival, 221 vs. 345 days *P* \<0.05; Figure [4](#Fig4){ref-type="fig"}B). These results suggest that the down-regulation of TES predicts worse outcome of GBM patients and that promoter hypermethylation should be the major reason for TES down-regulation, although mutation, loss of heterozygosity, etc. cannot be excluded.Figure 4**Survial analysis according to TES expression and methylation status.** Kaplan-Meier survival estimates of overall survival according to the TES methylation status in 42 GBM patients **(A)** and according to the TES expression level in GBM patients **(B)**.

The elucidation biomarkers of glioblastoma, the most common and lethal of brain tumors, has been reported since the 1990's. Somatic mutations have been demonstrated for *TP53* on 17p \[[@CR14]-[@CR16]\], *PTEN* on 10q \[[@CR17]-[@CR19]\], *RB1* on 13q \[[@CR20],[@CR21]\], and homozygous deletions of the *CDKN2A/p16/p14* locus on *9p* have been documented \[[@CR16],[@CR22]\]. Other putative tumor-suppressor genes alternated at all other loci remain elusive. Methylation of promoter regions inducing silencing of gene expression, has been implicated as an alternative mechanism for tumor-suppressor gene inactivation. In glioblastomas, frequent promoter hypermethylation has been reported for *PRDX1* \[[@CR23],[@CR24]\], *WNK2* \[[@CR24]\], *SGNE1* \[[@CR25]\], and *ALDH1A3* \[[@CR26]\]. Blough et al. demonstrated that 1p Loss silence NHE-1 in oligodendroglioma \[[@CR27]\]. Some studies showed the methylation of promoter regions with other alterations. Chou et al. identified hypermethylation of *RB1* with IDH1 and IDH2 mutation in gliomas \[[@CR28]\].

The *TES* gene was previously identified as located in the fragile chromosomal region FRA7G at 7q31.2 \[[@CR29]\]. Loss of TES expression has been frequently found in various cancers \[[@CR5],[@CR29]\]. Non-synonymous mutation homozygous deletions have not been observed, consistent with CpG promoter hypermethylation being a mechanism of *TES* gene inactivation \[[@CR5]\]. Previous studies have demonstrated that TES inhibited the growth of breast and uterine cells and promoted ovarian cancer cell apoptosis through caspase-dependent and caspase-independent processes \[[@CR9],[@CR30]\]. Although *TES* could be a candidate tumor suppressor gene in glioma cells, the mechanism of this role is still unknown.

TES, a highly conserved protein, consists of three C-terminal LIM domains and a PET (prickle, espinas, and testin) domain of unknown function. Group 3 LIM domain proteins, including the TES domain, are a component of the focal adhesion complex and localize to cell-matrix adhesions and cell-cell contacts \[[@CR31]\]. Over-expression of TES increases cell spreading and decreased cell motility \[[@CR32]\]. Adenoviral transduction of the *TES* gene into breast and uterine cancer cell lines promotes apoptosis and tumor reduction \[[@CR30]\]. Further functional analysis demonstrated that *TES* acted as a tumor suppressor gene *in vivo*, as *TES* knockout mice showed an increased susceptibility to induced gastric cancer \[[@CR6]\]. Our results in this study showed that *TES* is widely down-regulated in GBM. Elevation of TES expression can increase apoptosis of GBM cells.

TES can be silenced either by loss of heterozygosity or by promoter hypermethylation in several tumor types, including head and neck squamous cell carcinoma, ovarian cancer, primary gastric cancer, prostate cancer, breast cancer, among others \[[@CR8]-[@CR12]\].

Our data suggests promoter hypermethylation is the major cause of down-regulation or silencing of *TES* in GBM, even if other mechanisms cannot be excluded at this point.

Conclusions {#Sec15}
===========

Using whole-genome DNA methylation microarray we found that the *TES* gene is a novel tumor suppressor gene and might represent a valuable prognostic marker for glioblastoma, indicating a potential future target for GBM therapy.

**Competing interests**

The authors declare that they have no competing interests.

**Authors' contributions**

YB participated in the design of the study, sample collection, performed Western blotting, flow cytometry analysis of cell apoptosis, and immunohistochemistry analysis, and drafted the manuscript. QGZ performed the statistical analysis, carried out the immunohistochemistry analysis, and helped to revise the manuscript. XHW conceived of the study and participated in its design and coordination, and helped to draft the manuscript. All authors read and approved the final manuscript.

This work was kindly supported in part by the following grants: National Natural Science Foundation of China (81071626, 30872933).
