Histone demethylase GASC1 - a potential prognostic and predictive marker in invasive breast cancer

Background The histone demethylase GASC1 (JMJD2C) is an epigenetic factor suspected of involvement in development of different cancers, including breast cancer. It is thought to be overexpressed in the more aggressive breast cancer types based on mRNA expression studies on cell lines and meta analysis of human breast cancer sets. This study aimed to evaluate the prognostic and predictive value of GASC1 for women with invasive breast cancer. Methods All the 355 cases were selected from a cohort enrolled in the Kuopio Breast Cancer Project between April 1990 and December 1995. The expression of GASC1 was studied by immunohistochemistry (IHC) on tissue microarrays. Additionally relative GASC1 mRNA expression was measured from available 57 cases. Results In our material, 56% of the cases were GASC1 negative and 44% positive in IHC staining. Women with GASC1 negative tumors had two years shorter breast cancer specific survival and time to relapse than the women with GASC1 positive tumors (p=0.017 and p=0.034 respectively). The majority of GASC1 negative tumors were ductal cases (72%) of higher histological grade (84% of grade II and III altogether). When we evaluated estrogen receptor negative and progesterone receptor negative cases separately, there was 2 times more GASC1 negative than GASC1 positive tumors in each group (chi2, p= 0.033 and 0.001 respectively). In the HER2 positive cases, there was 3 times more GASC1 negative cases than GASC1 positives (chi2, p= 0.029). Patients treated with radiotherapy (n=206) and hormonal treatment (n=62) had better breast cancer specific survival, when they were GASC1 positive (Cox regression: HR=0.49, p=0.007 and HR=0.33, p=0.015, respectively). The expression of GASC1 mRNA was in agreement with the protein analysis. Conclusions This study indicates that the GASC1 is both a prognostic and a predictive factor for women with invasive breast cancer. GASC1 negativity is associated with tumors of more aggressive histopathological types (ductal type, grade II and III, ER negative, PR negative). Patients with GASC1 positive tumors have better breast cancer specific survival and respond better to radiotherapy and hormonal treatment.


Background
Breast cancer is a heterogeneous disease with different histopathological, molecular and clinical characteristics. Moreover, there is a wide variation in the progress of breast cancer in patients of the same age and with tumors of comparable clinical extent. The discovery of molecular markers, such as estrogen, progesterone and HER2 receptors has facilitated the classification of tumors and led to discovery of new cancer therapies. However, the selection of patients for appropriate adjuvant therapies still encounters difficulties [1]. Therefore there is an urgent need for novel diagnostic, prognostic and predictive markers which would make easier selection of patients for adjuvant therapies and possibly open novel perspectives for more efficient therapeutic strategies.
Nowadays epigenetics is making an increasingly important impact in cancer research [2][3][4]. Epigenetic research has not only provided novel insights into the molecular mechanisms of cancer, but it has also revealed useful diagnostic, prognostic and predictive biomarkers [5][6][7]. In addition to epigenetic modifications like DNA methylation and nucleosome positioning, histone modification patterns are altered in human tumors. However, methodological difficulties prevent use of altered histone modification profiles found in cancer as biomarkers [8]. Therefore, there is a considerable interest in understanding histone modifier genes and their products. One of these genes is GASC1 (gene amplified in squamous cell carcinoma 1; aliases: JMJD2C, JHDM3C, KDM4C) which codes a histone demethylase for di-and trimethylated lysine 9 and 36 on histone H3 (H3K9me3/2 and H3K36me3/2) [9,10]. H3K9me3/2 mark is generally associated with transcriptional repression and the formation of heterochromatin, while H3K36me3/2 is associated with transcriptionaly active genes and it is believed to play an important role in the suppression of incorrect transcription [10,11]. Because of its dual role in modifying H3 either by removing the repressive H3K9me3/2 or the active H3K36me3/2 factor, GASC1 has been considered to be a fine-tuning regulator of gene expression in normal development and differentiation as well as in cancer development and progression [12][13][14]. The involvement of GASC1 in development of breast cancer has been well documented in cell lines by Liu et al. [13] and Wu at al. [15]. They demonstrated that GASC1 is amplified and overexpressed in multiple breast cancer cell lines, it causes transformation of immortalized, non-transformed mammary epithelial cells, regulates expression of genes responsible for stem cells self-renewal and may be linked to the stem cell phenotypes in breast cancer. Additionally, they have found that the GASC1 is overexpressed in aggressive, basal-like breast cancers compared with non basallike breast cancers.
Moreover, an oncogenic role of GASC1 has been documented in prostate cancer where it enhances the transcription of androgen receptor-dependent genes and cell proliferation by interaction with ligand-bound androgen receptor [16]. GASC1 also plays an important role in normal development and differentiation by regulating expression of pluripotency genes, including NOTCH1, NANOG, Sox2 and Pou5 [17,18].
As far as we are aware, this marker has not been evaluated by immunohistochemistry in human breast tumors. Consequently, this study aimed at determining the relevance of GASC1 demethylase in the prognosis of invasive breast cancer progression and in the prediction of responses to particular adjuvant treatments in material from a large cohort of patients with a detailed clinical and histopathological classification of tumors and up to 20 years of follow-up [19]. Additionally, we investigated whether this marker could be utilized for more detailed classification of invasive breast cancer.

Methods
Our initial material consisted of 392 breast cancer cases selected from a cohort enrolled into the Kuopio Breast Cancer Project in the Kuopio University Hospital, Kuopio, Finland between April 1990 and December 1995 [20,21].
The tumor samples from these patients were fixed in 10% buffered formalin and embedded in paraffin. The histological diagnosis was confirmed by reviewing one to four original sections of the primary tumor. From the total material (392 tumors), we excluded 37 benign and in situ cases. In the remaining 355 cases of invasive breast tumors without distant metastases, we evaluated GASC1 expression by immunohistochemical (IHC) staining in the nuclei of the tumor epithelial cells (Figure 1). Further we analyzed how the GASC1 status would influence the breast cancer specific survival and the time to relapse.
In the analyzed by IHC material, 198 (55.8%) cases were GASC1 negative and 157 (44.2%) GASC1 positive. In this cohort, 105 patients had died of breast cancer, 106 patients had died from other causes than breast cancer and 144 patients were still alive at the time of analysis. The mean follow up time at the cut-off point in February 2011 was 10.6 years, ranging from 0.1 to 20.4 years. Sixty eight patients had undergone resection, 285 were treated with mastectomy, and two patients did not undergo surgery. Postoperative radiotherapy was given to 206 patients. Altogether 62 patients had received adjuvant tamoxifen, and 69 patients were treated with adjuvant chemotherapy, mainly the intravenous CMF regimen (500 mg/m 2 , methotrexate 40 mg/m 2 , 5-fluorouracil 500 mg/m).
The tissue microarray (TMA) was constructed as described previously [22]. The diagnosis of the cases was based on the World Health Organization (WHO), classification of breast and female genital organs [23]. The presence of metastases was determined at the time of the operation. The collection of the material and the clinical features of the patients have been described in a previous study [21]. The research was approved by the ethical committee of University of Kuopio/University of Eastern Finland and Kuopio University Hospital.

Immunohistochemistry for GASC1
Immunohistochemical staining was performed on 4 μmthick sections cut, from TMA block. After deparaffinization and rehydration, the sections were heated in a microwave oven for 3 × 5 min in citrate buffer (pH 6.0). Then they were treated for 5 min with 5% hydrogen peroxide to block endogenous peroxidase. Next, the sections were incubated for 35 min at room temperature in 1.5% normal serum diluted in PBS to block non-specific binding. After that, the sections were incubated overnight at 4°C with the mouse monoclonal anti -GASC1 antibody (Origene, TA 500587) at dilution 1:100. The slides were then incubated with a biotinylated secondary antibody (35 min) and avidin-biotin-peroxidase complex (45 min) (ABC Vectastain Mouse Elite Kit, Vector Laboratories, Burlingame, CA, USA). After each step of the immunostaining procedure the slides were rinsed with PBS. The color was developed with diaminobenzidine tetrahydrochloride (DAB) (Sigma, St. Louis, MO, USA). The slides were counterstained with Mayer's haematoxylin, washed, dehydrated, cleared and mounted with Depex (BDH, Poole, UK). In the negative controls, the primary antibody was omitted.
The immunoreactivity for GASC1 was analyzed in the nuclei of epithelial tumor cells taking into account the number of positively stained nuclei and intensity of staining. The number of positively stained nuclei was semiquantified as follows: The intensity of nuclear staining was evaluated on a four grade scale: The evaluation was performed on two separate array cores (A and B) by two pathologists (YS, BB) blinded to outcome ( Figure 1). The results in series A and B were similar (substantial inter-series agreement was achieved; kappa for both nuclear number and intensity was 0.7; p=0.000). The final score was obtained by combining these results into four groups. Tumors were designated as negative if their scores from both series A and B were 0. Tumors were designated as positive if their score from at least one series was positive.

mRNA expression analysis
For this study from this cohort there were 57 samples available with invasive breast tumors. RNA was isolated from fresh frozen tissue stored at −70°C using Protein and RNA Isolation System for Small RNAs (mirVana TM Paris TM ). We used High Capacity cDNA Reverse Transcription Kit to synthesize cDNA (Applied Biosystems, Foster City, USA) following the manufacturer's instructions. The endogenous control gene was chosen by the investigation conducted by McNeill et al. [24] where they determined PPIA to be the best choice for breast cancer mRNA analysis. Next we performed quantitative real-time RT-PCR using gene specific TaqMan W Gene Expression Assays (Applied Biosystems). Relative gene expression values were calculated as the ratio between the target gene and the endogenous control PPIA, obtained for each sample from the standard curves. Finally, the values greater than the mean were designated as positives and lower than the mean as negatives.

Statistical analysis
The expression of GASC1 in different groups was compared using chi-squared test. In breast cancer specific survival analysis, the end point was death from breast cancer with deaths from other causes being censored, whereas in time to relapse analysis, the end point was breast cancer recurrence, either local or distant. Overall survival was calculated from the date of diagnosis to the date of death or the last follow-up date. Kaplan-Meier analysis was applied to estimate breast cancer specific survival and overall survival according to the adjuvant treatment; different groups were compared with the logrank test. Multivariate analyses were conducted with the Cox regression model. The statistical analyses were performed using SPSS version 17.0 (SPSS Inc., Chicago, IL, USA). p < 0.05 was considered statistically significant.

GASC1 IHC negativity is an independent prognostic factor of poorer breast cancer specific survival
Overall, women with GASC1 negative tumors (n=198) had two years shorter breast cancer specific survival than the women with GASC1 positive tumors (n=157; Table 1, Figure 2). Stratification according to clinical parameters revealed that the GASC1 negative women, over 55 years of age, with tumor size T2, T3 or T4, with positive nodal status and with clinical stage II, III or IV had a significantly poorer survival than the GASC1 positive ones. In patients aged 55 years or younger and patients with clinical stage I, the GASC1 status did not influence the survival time (Table 1).
After including the above clinical parameters as covariates in the Cox regression analysis, we confirmed that GASC1 negativity was an independent factor predicting poorer breast cancer specific survival, equal to the positive nodal status (p=0.001, Table 2, Figure 3). In this analysis tumor size (T2, T3 and T4) and clinical stage (II, III and IV) had no effect on breast cancer specific survival and time to relapse possibly, because there were only 34 patients (9.6%) with T3 and T4 tumors and 42 patients (12%) at Stage III and Stage IV in our material. Probably, this number of cases was not sufficient to show a significant influence on survival in multivariate analysis. However, in univariate analysis the patients with T1 and at Stage I survived significantly better and had significantly longer time to relapse than the patients with more advanced disease. Similarly, in bivariate analysis, where as the second variable in addition to tumor size or Stage was entered GASC1 status, tumor size or stage and GASC1 status had a significant effect on the breast cancer specific survival and time to relapse.

GASC1 IHC negative cases have a shorter time to relapse than the GASC1 positive cases
Overall, there were 132 cases with relapse and breast cancer was the cause of death of 105 women. Eleven women with a relapse died from other causes and 16 were still alive at the time of analysis. Women with GASC1 negative tumors had a shorter time to relapse than the women with GASC1 positive tumors (Table 3, Figure 4). After stratification according to the clinical parameters, the results were similar to those obtained in breast cancer specific survival analysis, which suggests that among subjects with more advanced stages of the disease, GASC1 negative cases had poorer prognosis (Table 3). Cox regression analysis verified that the GASC1 negativity is an independent factor predicting a shorter time to relapse in women with invasive breast cancer (p=0.002, Table 2, Figure 5).
GASC1 IHC negative cases are more likely to have a relapse of breast cancer and to suffer from more aggressive tumors than the GASC1 positive cases Among patients with a relapse (n=132) there were 26% more GASC1 negative cases (n=83) than GASC1 positive (n=49), while among patients without relapse (n=223) the difference between GASC1 negative (n=115) and positive (n=108) cases was only 4% (chi2, p=0.038). The majority of GASC1 negative tumors were ductal (72 %) of higher histological grade (84% of grade II and III altogether). The GASC1 negative tumors significantly differ in terms of histological type and grade, estrogen receptor (ER), progesterone receptor (PR) and HER2 status from GASC1 positive tumors (p=0.005, p=0.000, p=0.033, p=0.001 and p=0.029 respectively; Table 4). The aforementioned results pointed to the possibility that GASC1 negativity might only be a marker for aggressive tumor subtype but not an independent marker of breast cancer specific survival. To check this possibility, we performed the Kaplan-Meier survival nnumber of cases, HRhazard ratio.

Cum Survival
GASC1 IHC positivity is an independent marker for better prognosis in patients treated with radiotherapy or tamoxifen We also analyzed the predictive value of GASC1 staining according to the adjuvant treatments that the patients were given. There were 206 patients who were treated with adjuvant radiotherapy. Kaplan-Meier analysis detected a better relapse-free survival (Log Rank 0.017) and breast   cancer specific survival (Log Rank 0.021) in the patients with GASC1 positive tumors. In the Cox regression analysis after adjusting for age, stage, chemotherapy, hormonal therapy and hormone receptor status, GASC1 positivity was statistically significantly associated with improved relapse-free and better breast cancer specific survival ( Table 6). Overall survival was not affected by the GASC1 staining intensity.
Altogether sixty two patients with ER positive tumors were treated with tamoxifen as their only adjuvant    (Table 6). GASC1 staining status did not have any effect on the survival of patients given chemotherapy as their only medical adjuvant treatment (n=69). Furthermore, there were no significant differences in overall survival in any of the treatment groups (Table 6).

GASC1 mRNA expression is in line with the immunohistochemical data
The expression of GASC1 mRNA was evaluated in 57 available cases from the material used in tissue microarrays (TMA). The Kaplan-Meier survival curves were similar to those obtained from protein analysis: cases with low expression of GASC1 mRNA had 2.3 years shorter breast cancer specific survival than cases with high GASC1 mRNA expression (p=0.132, Log Rank; Figure 6). GASC1 mRNA expression in GASC1 negative cases was significantly lower than in GASC1 positive ones (Mann-Whitney: p=0.004). The expression of GASC1 mRNA was significantly lower in grade II and III tumors compared with its expression in grade I tumors (Figure 7). This supports the previous finding from TMA analysis showing that the great majority (84%) of tumors of grade II and III were GASC1 negative.
The cases with negative or weak PR expression had a significantly lower relative level of GASC1 mRNA than the cases with high PR expression (Mann-Whitney: p=0.016). In contrast, HER2 negative cases showed significantly higher GASC1 mRNA expression than the HER2 positive counterparts (Mann-Whitney: p=0.004), which was in line with the protein staining results (Figure 8).

Discussion
This is the first study investigating GASC1 status in a relatively large group of clinical samples (altogether 355 cases of invasive breast cancer). We show that GASC1 negativity is an independent prognostic factor of poorer    breast cancer specific survival, not only overall but also in groups of patients treated with different adjuvant therapies, evidence of its prognostic as well as predictive value for women with invasive breast cancer.
GASC1 negative women had poorer breast cancer specific survival, when they were in more advanced stages of the disease (T2 -T4, positive lymph nodes, clinical stages II -IV). However, the survival of women in initial stages of the disease did not differ between GASC1 positive and negative cases, even though the percentage of GASC1 positive tumors in these women was similar to the percentage in women in more advanced stages of the disease (results not shown). This observation points to a modulating role for GASC1 in the regulation of gene expression rather than a clear-cut activating or inhibitory role. An interesting finding from the survival analysis was that HER2 negative patients survived significantly better when they had GASC1 positive tumors regardless of their clinical stage. This observation might point the search for new therapeutic possibilities in the direction towards the GASC1 signaling pathway. Identification of molecules involved in this pathway and elucidating their role in the breast cancer pathophysiology could well be beneficial especially for those patients with triple negative tumors.
Moreover, we observed that GASC1 negative cases overall and in advanced clinical stages had a shorter time to relapse than the GASC1 positive patients, which was in agreement with the breast cancer specific survival analysis. Further, we demonstrated that the GASC1 positivity was statistically significantly associated with improved relapse-free and breast cancer specific survival in patients who were treated with adjuvant radiotherapy or with adjuvant tamoxifen which supported the above results. Those patients with ER-positive, HER2-negative disease are the group in whom decisions about adjuvant chemotherapy are most difficult. The relative indications for chemoendocrine therapy and endocrine therapy alone are given by the St Gallen International Expert Consensus [1]. However, there are no guidelines for those ER-positive, HER2-negative patients who are node positive with grade II tumors of size 2.1 -5 cm (T2-T3). Our results indicate that the patients with HER2 negative, grade II tumors of ductal type have better prognosis when they are GASC1 positive (Table 5). Moreover GASC1 positive patients respond better to hormonal treatment. The results suggest that GASC1 positivity in these patients might be an indication for endocrine therapy alone, however this will need to be confirmed in a larger group of patients.
The majority of GASC1 negative tumors were ductal cases of higher histological grade, more often ER and PR negative and HER2 positive. We confirmed these results by analysis of GASC1 mRNA expression and we also excluded the possibility that GASC1 is only a marker for aggressive tumor subtype. As far as we are aware, there are no studies of IHC GASC1 expression in human breast cancer and therefore, we analyzed the list of genes which were affected by GASC1 depletion in study of Loh at al. [17] (Supplementary Table 1 by Loh at al. [17]). Several genes which were up-regulated by GASC1 depletion in this study have also been found to be upregulated in certain human tumors. Gene Sp5 (fold change 4.8) was reported to be overexpressed in several human tumors including hepatocellular carcinoma, gastric cancer, and colon cancer [25,26]. Annexin A1 A B (ANXA1, fold change 7.1) has been shown to have oncogenic potential in breast cancer progression and metastasis [27] and was significantly correlated with unfavorable prognostic features of breast cancer [28].
Insulin-like growth factor-2 (IGF-2, fold change 4.5) and H19 (fold change 4.2) are located together within an imprinted domain of chromosome 11p15.5. Loss of imprinting of IGF2 resulting in increased expression of IGF-2 is a common genetic alteration in human malignancies and aberrant methylation of IGF2/H19 locus has been detected in multiple human cancers including breast cancer [29][30][31]. Overexpression of IGF-2 was also positively correlated with increased risk of breast cancer, nodal positivity and higher tumor grade [32]. Forkhead box Q1 (Foxq1, fold change 6.7) gene expression was correlated with high-grade basal-like breast cancers where it was associated with metastases and poor clinical outcomes [33,34].
The above evidence suggests that the absence of GASC1 might cause over-expression of several genes responsible for carcinogenesis and a more aggressive tumor phenotype in different cancers including breast cancer, which supports our finding that GASC1 negativity is correlated with more aggressive breast tumors and poorer prognosis. Based on these results, we speculate that GASC1 and possibly its up-and/or downstream factors might help to define breast tumors with better or poorer prognosis, especially in those which are considered to be more aggressive. However this hypothesis needs confirmation in a larger group of patients. GASC1 possesses enzymatic activity i.e. it specifically demethylates tri-and dimethylated (me3 and me2) lysine (K) in residues 9 and 36 on histone H3 (H3K9me3/2 and H3K36me3/2) [9,16,35] and together with histone methyltransferases, it dynamically modulates the methylation status of H3. The biological significance of H3 methylation depends on the modified residue (K9 or K36), the degree of methylation (me1, me2, me3) and the genomic position of H3 (promoter region or coding region) [36,37]. In general, H3K9me3/2 is found in promoter regions of inactive genes and demethylation in this site triggers promoter activation, whereas H3K36me3/2 is enriched within the body of active genes and demethylation in this site is related to the termination of transcription [10,38]. However, the increase of H3K9me3 inactive mark in coding region of the gene has also been associated with active gene expression [39].
The transcriptional repressive effect of the H3K9 methylation can be attributed to the association with the repressive protein HP1 (heterochromatin protein 1) [11]. Cloos at al. [9] reported that GASC1 could effectively demethylate H3K9me3 and H3K9me2, release the repressive protein HP1 and reduce heterochromatin in vivo. One could speculate that in the absence of GASC1 HP1 would be recruited to H3K9me3/2 and stabilize heterochromatin. Recent data have shown that HP1α is over-expressed in numerous cancers and that HP1α over-expression is associated with increased cell proliferation, most likely through silencing of genes inhibiting cell proliferation. Moreover, the same authors demonstrated that HP1α overexpression in breast cancer patient samples correlated with an increased risk of death [40]. These findings could explain, at least to some extent, why GASC1 negativity is associated with more aggressive tumors and poorer breast cancer specific survival.
H3K9me3 status is also mediated by different histone methyltransferases. Suv39H1 (suppressor of variegation 3-9 homolog 1) is a major methyltransferase responsible for H3K9me3 that is intimately linked to DNA methylation. Dong at al. (2012) have found that H3K9me3 and DNA methylation on the E-cadherin promoter were higher in basal-like breast cancer cell lines. Furthermore, they showed that knockdown of Suv39H1 restored Ecadherin expression by blocking H3K9me3 and DNA methylation and resulted in an inhibition of cell migration, invasion and metastasis of basal-like breast cancer [41]. GASC1 activity eliminates H3K9me3, which in view of the above data, might partly explain the protective role of GASC1 evident in our data.
The role of GASC1 in prostate cancer has been studied by Wissmann at al. [16]. GASC1 colocalizes with androgen receptor (AR) in both normal prostate and prostate carcinomas. It enhances the transcription of AR responsive genes, and this property may contribute to tumor cell proliferation. In primary breast carcinomas, including triple-negative breast cancer, AR expression was interestingly associated with a significantly better disease free interval and overall survival [42][43][44]. Moreover, AR positive breast tumors have been shown to be smaller, more often node negative, lower in histological grade and clinical stage than their AR negative counterparts [42][43][44][45]. The clinico-pathological characteristics of AR positive tumors in the above studies resemble the characteristics of the GASC1 positive tumors in our study. This suggests that the GASC1 and the AR might be functionally connected also in breast cancers. However, the effect of AR stimulation in breast tissue is reported to be opposite to AR stimulation in prostate [46].
Regardless of repressing or activating AR dependent genes, GASC1 as an epigenetic factor can have a more general influence on gene expression and it also may modulate the expression of genes which are not dependent on the AR. GASC1 demethylase activity has been shown to regulate the expression of genes critical for stem cell self-renewal, including NOTCH1 and NANOG, and possibly to be linked to the stem cell phenotypes in breast cancer [13,17,18]. The mammary gland is a subject to many tissue remodeling events occurring during puberty, pregnancy, lactation and menopause. Therefore, there needs to be a compartment of normal stem cells or early progenitor cells with high proliferative potential and differentiation abilities in order to maintain the mammary gland function. It has been proposed that breast tumors arise from normal stem cells or early progenitor cells through deregulation of normal self-renewal [47]. It is possible that GASC1 maintains the normal population of stem cells through mechanisms connected with HP1 regulation as described above. This function of GASC1 might be responsible for a lower recurrence rate in GASC1 positive cases compared with the GASC1 negative and the better outcome of GASC1 positive patients treated with radiotherapy.
The results obtained from clinical material are not consistent with in vitro results from breast cancer cell lines. Liu at al. [13] reported that GASC1 expression was significantly higher in those cell lines representing aggressive basal-like breast cancers compared with the cell lines representing nonbasal-like breast cancers. Moreover, they demonstrated that GASC1 was able to induce transformed phenotypes when overexpressed in immortalized, non-transformed ER negative mammary epithelial MCF10A cells. In our clinical material, GASC1 positivity was associated with less aggressive tumors. The reason for this discrepancy might be that cell lines only generally reflect the pathophysiological environment of the tumor without including all the up-and downstream factors present in the organism that influence the tumor behavior. Liu and coworkers [13] linked the above oncogenic properties of GASC1 with induction of NOTCH1 by GASC1. However, there is evidence that NOTCH proteins can act either as oncogenes or as tumor suppressors depending on the cellular context [48,49]. One reason for the difference between our results and the report of Liu at al. [13] is that the latter did not take into account the cellular context of GASC1 action i.e. stromal or physiological interactions were absent in the cell cultures.
In addition to cell culture experiments, Liu at al. [13] meta-analyzed data sets concerning gene expression in human breast cancer from ONCOMINE. The results concerning GASC1 mRNA expression from this metaanalysis were consistent with their results obtained from cell culture, but did not support our results from IHC analysis of GASC1 expression. The discrepancy might be attributed to differences in the studied material. The data set from Fiank at al. [50] concerned gene expression in tumor stroma, while we assessed GASC1 expression in the epithelial compartment of the tumor. The analysis which was based on material used by Kreike at al. [51] consisted of 58% triple negative/basal like carcinomas, while in our material, triple negative tumors constituted only 11% of the total [51]. The discrepancy might also come from differences in methodology. We performed immunohistochemistry and compared numbers of GASC1 positive and negative tumors in different histopathological and clinical groups, while Liu at al. based their work on GASC1 mRNA and protein expression analysis in different experimental setups. Our sample size available for GASC1 mRNA expression study was too small to permit a detailed statistical analysis, which prohibits a relevant comparison between our GASC1 expression results with the GASC1 expression evaluated by others.

Conclusions
In summary, this study shows that the histone demethylase, GASC1, can have clinical significance as a prognostic as well as a predictive factor, indicating better prognosis and better response to radiotherapy and adjuvant treatment. GASC1 can possess either tumor suppressing or oncogenic properties, which may be context and organ specific.