Skip to main content

5-AZA-dC induces epigenetic changes associated with modified glycosylation of secreted glycoproteins and increased EMT and migration in chemo-sensitive cancer cells

Abstract

Background

Glycosylation, one of the most fundamental post-translational modifications, is altered in cancer and is subject in part, to epigenetic regulation. As there are many epigenetic-targeted therapies currently in clinical trials for the treatment of a variety of cancers, it is important to understand the impact epi-therapeutics have on glycosylation.

Results

Ovarian and triple negative breast cancer cells were treated with the DNA methyltransferase inhibitor, 5-AZA-2-deoxycytidine (5-AZA-dC). Branching and sialylation were increased on secreted N-glycans from chemo-sensitive/non-metastatic cell lines following treatment with 5-AZA-dC. These changes correlated with increased mRNA expression levels in MGAT5 and ST3GAL4 transcripts in ovarian cancer cell lines. Using siRNA transient knock down of GATA2 and GATA3 transcription factors, we show that these regulate the glycosyltransferases ST3GAL4 and MGAT5, respectively. Moreover, 5-AZA-dC-treated cells displayed an increase in migration, with a greater effect seen in chemo-sensitive cell lines. Western blots showed an increase in apoptotic and senescence (p21) markers in all 5-AZA-dC-treated cells. The alterations seen in N-glycans from secreted glycoproteins in 5-AZA-dC-treated breast and ovarian cancer cells were similar to the N-glycans previously known to potentiate tumour cell survival.

Conclusions

While the FDA has approved epi-therapeutics for some cancer treatments, their global effect is still not fully understood. This study gives insight into the effects that epigenetic alterations have on cancer cell glycosylation, and how this potentially impacts on the overall fate of those cells.

Graphic abstract

Background

Breast cancer is the leading cause of cancer-related deaths in women; triple negative breast cancer (TNBC) subtype, accounting for 10–20% of the cases [1], is more aggressive and harder to treat [2]. Although ovarian cancer is far less common, more than 70% of women present at a late stage, giving them a less than 40% overall survival rate [2]. DNA methylation is altered in cancer, and DNA methyltransferases are attractive targets for epi-therapies such as Vidaza (5-AZAcytidine) and Decitabine (5-AZA-2′-deoxycytidine) (https://clinicaltrials.gov/ct2/home) [2], which were approved by the FDA for the treatment of myelodysplastic syndrome [3]. Global hypomethylation and site-specific hypermethylation have been observed on the DNA of both breast and ovarian tumours [4]. Promoter hypermethylation results in gene silencing as a consequence of transcription factors (TFs) being unable to bind to a promoter/enhancer region of a gene [5]. The hypermethylation of genes involved in cell cycle control or tumour suppression causes these genes to be silenced and promotes cancer initiation and progression.

Glycosylation is the most common posttranslational modification of proteins and is altered in cancer and other pathologies such as chronic inflammatory diseases exemplified by rheumatoid arthritis [6]. In cancer, the predominant glycan changes reflect increases in highly sialylated and branched glycans, which are associated with inflammation, metastasis and disease progression [6]. In a previous publication, Saldova et al. [7] demonstrated increases in highly sialylated and branched glycans on the N-glycans from secreted glycoproteins from the ovarian cancer cell line OVCAR3 post-5-AZA-dC treatment [7]. Similar results were published by Chachadi et al. [8]. Specifically, they reported enhanced glycosylation changes post-5-AZA-dC treatment that were associated with enhanced metastatic potential [8]. Both papers not only highlight the role that DNA methylation plays in glycosylation, but also the need for a more in-depth look at the effects of 5-AZA-dC on glycosylation patterns in various tumour tissues, and on tumours of varying grade and stage [9].

The effects on cellular processes and end points when epigenetic alterations occur are well documented [10, 11]. However, with respect to altered glycosylation, processes such as cellular apoptosis and cellular senescence are less well understood. Research has been conducted, where cellular fate has been altered using various compounds to specifically target glycosylation. Gwak et al. [12] have shown that inhibiting N-linked protein glycosylation, using resveratrol, caused an accumulation of unfolded proteins in the endoplasmic reticulum (ER) of ovarian cancer cells and induced ER stress-mediated apoptosis [12]. In another study, NGI-1, an oligosaccharyltransferase (OST) inhibitor, was used to induce senescence in non-small cell lung carcinoma cells, by blocking the cell surface localisation of EGFR through altering glycosylation [13]. It is therefore warranted to consider that epigenetic tags and glycosylation should not be looked at as separate entities, but in combination, as they are concomitantly altered in cancer.

The current study follows up on the published work by Saldova et al. [7]. Our objective was to investigate whether changes in DNA methylation result in altered glycosylation, triggering differential cellular fate, in two chemo-resistant/chemo-sensitive pairs of ovarian and two triple negative breast cancer (TNBC) cell lines in-vitro. The impact of this treatment on apoptosis, senescence and epithelial to mesenchymal transition (EMT) is described. Survival analyses on online data repositories (kmplot.com) were undertaken to access the effect of glycan alterations in progression of ovarian and triple negative breast cancer (TNBC).

Results

5-AZA-dC treatment alters global DNA methylation, the cell cycle and cisplatin sensitivity

Two pairs of chemo-sensitive/chemo-resistant ovarian cancer (A2780/A2780cis, PEO1/PEO4), and two TNBC (MDA-MB-231 and MDA-MB-436) cell lines were chosen for this study, to determine what effect 5-AZA-dC treatment has on the secreted glycome of chemo-sensitive/chemo-resistant cancer cell lines. The inclusion of two TNBC cell lines (MDA-MB-231 and MDA-MB-436) allowed for comparison with another tumour type. Briefly, A2780 and A2780cis cell lines were treated with 1 μM 5-AZA-dC and the remaining cell lines were treated with 0.1 μM 5-AZA-dC. Following 5-AZA-dC treatment, cells were harvested for flow cytometry dot plot analyses (Fig. 1a). Global DNA methylation was assessed using an anti-5′-methylcytidine antibody. Analyses showed a reduction in DNA methylation in all 6 cell lines, ranging from 22.16% in PEO1 cells, to 66.26% in A2780cis (Fig. 1b). The demethylation was higher in the chemo-resistant (A2780cis, PEO4), compared to the chemo-sensitive (A2780, PEO1) cells, significantly in PEO1 compared to PEO4 (p = 0.033) (Fig. 1b).

Fig. 1.
figure 1

5-AZA-dC successfully demethylates global DNA. a Representative example of a flow cytometry dot plot for the PEO1 cell line, both untreated (UnT) and treated with 0.1 µM of 5-AZA-dC, measuring the fluorescent intensity (FI) of both PI and FITC. b Total percentage of DNA methylation post-5-AZA-dC treatment. c Histogram profiles, representing the percentage phases of the cell cycle (G1, S, and G2), before and after 5-AZA-dC treatment. All experiments are n = 3 biological replicates. Error bars represent ± SD. Significant changes are starred: *p value ≤ 0.05 or **p value ≤ 0.005 (T test)

Post-5-AZA-dC treatment, the phases of the cell cycle were significantly altered in the A2780cis, PEO1 and MDA-MB-231 treated cells, compared to untreated controls. The G1 phase was decreased for all 6 cell lines, reaching statistical significance for PEO1 and MDA-MB-231 treated cells. The S phase was significantly increased in PEO1 and MDA-MB-231 treated cells, and G2 phase was significantly increased in A2780cis, PEO1 and MDA-MB-231 treated cells (p ≤ 0.05) (Fig. 1c).

Chemo-sensitive A2780 cells treated with 5-AZA-dC and cisplatin showed an increase in cell numbers (p = 0.044), whereas the chemo-resistant A2780cis cells showed a decrease in cell numbers (p = 0.004). This suggests that 5-AZA-dC treatment increased resistance to cisplatin in the A2780 cell line while enhancing sensitivity to cisplatin in the A2780cis cell line (Additional file 5: Figure S1).

The effect of 5-AZA-dC treatment on glycosylation is cell line specific

The HILIC-UPLC chromatograms from secreted glycans were separated into 35 peaks (Fig. 2a) and from the cell glycans into 39 peaks (Additional file 5: Figure S2A). The major glycans present in each peak from the secreted glycoproteins are listed and matched to the main peaks in our previously published paper [14] Table 1. The major glycans were pooled to the common features in each cell line. Branching, galactosylation and sialylation were calculated Additional file 2: Table S1. The preliminary assignments of the cell glycans and pooling into the common features using HILIC-UPLC and exoglycosidase digestions are shown in Additional file 3: Table S2. Cell and secreted glycan had differences in profiles. Specifically, while secreted glycans mostly contained highly branched and sialylated glycans, the cell glycans contained mostly oligomannosylated glycans.

Fig. 2
figure 2

Glycosylation changes are cell line specific. a Representative UPLC chromatograms produced from secreted N-glycans of ovarian (A2780, A2780cis, PEO1 and PEO4) and triple negative breast cancer (TNBC) cell lines (MDA-MB-231 and MDA-MB-436) and their separation into 35 peaks. b Plotted peak areas from the secreted N-glycans of ovarian and breast cancer cell lines. The glycans in each peak (GP1–GP35) and features are listed in Additional file 2: Table S1. Significant changes (p  < 0.05) are starred: *P value ≤ 0.05 or **p value ≤ 0.005. (MANOVA). Heatmap histograms indicating fold changes in 5-AZA-dC treated compared to untreated cells, were created using Hierarchial Clustering Explorer HCE 3.5 software. Blue indicates decreases, and red indicates increases. The shade of colour corresponds to amounts of the decreases/increases

Table 1 Main glycans on secreted glycoproteins in all cell lines

When comparing chemo-sensitive with chemo-resistant cell lines, some differences in secreted glycans were observed (Additional file 5: Figure S3), namely, an increase in GP13 (M5A2BF1G2Lac1S(6)1) in both chemo-resistant cell lines, A2780cis compared to A2780 (p ≤ 0.005) and PEO4 compared to PEO1 (p ≤ 0.05), and an increase in GP14 (A2G2S2 isoforms) in PEO4 compared to PEO1 (p ≤ 0.005). There were no significant differences in cell glycans when comparing A2780 and A2780cis cell lines.

In the secreted glycans, the variation in glycans between 5-AZA-dC treated and untreated samples appeared to be cell line specific (Fig. 2b). Specifically, when looking at the individual glycan peaks (GPs), significant changes in 5-AZA-dC treated samples (*, **) are observed in the MDA-MB-436 cell line, in which there is a significant decrease in GP15, 18 (A2G2S2 isoforms) (p ≤ 0.05) and GP24, 26, 29 (A3G3S3 isoforms) (p ≤ 0.005) (up to 0.51fold) (Table 1). Although not statistically significant, only the chemo-sensitive cell lines (A2780 and PEO1) show an increased trend (fields in red, Fig. 2b) in peaks containing triantennary trigalactosylated glycans, mostly trisialylated, but also disialylated (GP22 = A3G3S2; GP24, 26, 29 = A3G3S3). Chemo-resistant cells (A2780cis and PEO4) and TNBC cells (MDA-MB-231 and MDA-MB-436), show a decrease in these highly branched and sialylated glycans. Significant changes (*, **) were observed in the TNBC cell line MDA-MB-436, with a significant increase in the core fucosylated glycan (CoreF), a decrease in tetragalactosylated glycans (G4), and a decrease in tri and tetra sialylated glycans (S3 and S4). Although not statistically significant, only the chemo-sensitive cell lines (A2780 and PEO1) showed decreases (in the blue colour) in core fucosylated glycans and glycans with polylactosamine extentions with outer arm fucose as well as increases (in the red colour) in trisialylated glycans (changes varied from 0.69fold decrease up to 1.75fold increase).

When branching, galactosylation and sialylation were summarised, and MDA-MB-436 cells showed a significant (*) decrease in galactosylation (p ≤ 0.05). Although not significant, chemo-sensitive cell lines showed increases in branching, galactosylation and sialylation, whereas the chemo-resistant ovarian (A2780cis and PEO4) and TNBC cells showed the opposite trend, specifically, a decrease in branching, galactosylation and sialylation) (Fig. 2b).

In the cell glycans, peaks GP37 (FA4G4S(6,6,6)3), GP39 (FA4G4S(6,6,6,6)4) and tetrasialylated glycans (S4) were significantly decreased in treated A2780 cells (Additional file 5: Figure S2B).

5-AZA-dC treatment increases EMT and migration in chemo-sensitive ovarian cancer cell lines A2780 and PEO1

Having identified changes in the secreted N-glycan structures post-5-AZA-dC treatment, we next investigated what impact 5-AZA-dC treatment had on tumourigenic phenotypes. Firstly, markers related to the EMT specifically, E-cadherin, which is downregulated in EMT and N-cadherin and Vimentin, upregulated in EMT were investigated (Fig. 3a) [15]. A2780, MDA-MB-231 and MDA-MB-436 cell lines expressed high levels of N-cadherin with no E-cadherin detectable before the treatment suggesting a marked mesenchymal phenotype (Fig. 3a). A2780 and PEO1 both showed a significant increase in vimentin (p ≤ 0.05, p < 0.005, respectively) and PEO1 had a trend decrease in E-cadherin (p = 0.051) post-5-AZA-dC treatment compared to untreated controls. This indicated an increase in EMT in these chemo-sensitive cell lines (Fig. 3a and Additional file 5: Figure S4). The chemo-resistant cell line PEO4 showed a significant increase in E-cadherin compared to un-treated controls (p ≤ 0.05), suggesting a decrease in EMT (Fig. 3a and Additional file 5: Figure S4). The other chemo-resistant cell lines (A2780cis and TNBC cells lines MDA-MB-231 and MDA-MB-436) demonstrated no significant changes with 5-AZA-dC treatment (Fig. 3A, Additional file 5: Figure S4).

Fig. 3
figure 3

EMT and cell migration increases in chemo-sensitive cell lines. a Fold change of the protein expression of E-cadherin, vimentin and N-cadherin normalised to α-tubulin as determined by densitometry (Image J) (Western blots are in Additional file 5: Figure S4). Each condition represents three biological replicates. b Quantification of the absorbance readings from the Oris™ migration assay. A comparison of the absorbance from untreated vs 5-AZA-dC treated cells. Error bars were calculated from 3 independent experiments with 4 technical replicates per experiment. Significant changes are starred: *p value ≤ 0.05 or **p value ≤ 0.005 (T test)

Cell migration (Fig. 3b) of the 5-AZA-dC treated cell lines compared to untreated controls was then investigated. The Oris migration assay results showed a significant increase in migration in A2780 (p ≤ 0.005) post-5-AZA-dC treatment (Fig. 3b). None of the other cell lines showed any difference in migration post-5-AZA-dC treatment. However, although not significant, the PEO1 cell line had an increasing trend in migration (Fig. 3b). To see how the migration relates to the proliferation, faster growing cells (A2780, A2780cis and MB-MDA-231) were seeded at half the cell concentrations (2.5 × 104 cells) compared to the slow growing cells (5 × 104 cells) (PEO1, PEO4 and MB-MDA-436). Additionally, A2780cis cells, having a relatively high proportion of cells in S-phase (Fig. 1c), were treated with the solubilised dye after the migration results were measured, to see the amounts of the cells in each 5-AZA-dC treated and untreated wells. These results showed that there was no contribution of proliferation to the migration results (Additional file 5: Figure S5).

Evaluation of senescence and apoptosis post-5-AZA-dC treatment

To determine what potential impact 5-AZA-dC had on inducing senescence or apoptosis, PARP was analysed as a marker of apoptosis and p16, retinoblastoma (Rb) and p21 for cellular senescence. 5-AZA-dC treatment induced a significant increase in p21 expression in all cell lines (p ≤ 0.05) compared to untreated controls. For Rb, there was a significant decrease in expression in PEO1 and PEO4 cells compared to untreated controls. In relation to apoptosis, there was a significant increase in PARP cleavage in the MDA-MB-231 (p ≤ 0.05) post-5-AZA-dC treatment. There were no significant differences found for p16 across all cell lines and treatments (Fig. 4 and Additional file 5: Figure S6).

Fig. 4.
figure 4

5-AZA-dC treatment results in increased expression of the cellular senescence marker p21 and increased PARP-cleavage. Protein expression of p16, p21, Rb and PARP normalised to α-tubulin as determined by densitometry (Image J) (Western blots are represented in Additional file 5: Figure S6). Each condition was undertaken on two biological replicates. Significant changes are starred: *p value  ≤ 0.05 or **p value  ≤ 0.005 (T test)

GATA2 and 3 transcription factors regulate expression of MGAT5 and ST3GAL4 glycosyltransferases

To further investigate the glycosylation changes seen in the UPLC profiles, the expression levels of the various glycosyltransferases and sugar donor enzymes were analysed (Fig. 5a). Preliminary results on one biological replicate (2–3 technical replicates) showed alterations in expression in all the enzymes analysed in all cell lines. Consistent with increases in highly branched and sialylated glycans, the expression of MGAT5 and ST4GAL3, glycosyltransferases responsible for branching and sialylation, was also increased in chemo-sensitive cell lines A2780 and PEO1 post-5-AZA-dC treatment. Therefore, this suggested that these enzymes could be potentially regulated by promoter methylation. However, through preliminary analyses, only one of the two MGAT5 RefSeq curated transcript variants has a CpG island (variant 1) and ST3GAL4 is not subject to promotor DNA methylation (Supporting data1). Briefly, using data from the ENCODE project, we identified transcription factors (TFs) with binding sites in the proximal promoter regions of both MGAT5 and ST3GAL4. While GATA1-3 were not the only TFs with binding sites in the promoters of both MGAT5 and ST3GAL4, the finding fact that all three were present in both promoters suggested to us that these TFs may play an important role MGAT5 and ST3GAL4 transcription and therefore warranted further investigation [14]. GATA1 showed very low expression in all cell lines (data not shown). However, increased GATA2 levels correlated with increased ST3GAL4 levels (coef = 0.943; p = 0.005), and increased GATA3 levels correlated with increased MGAT5 levels (coef = 0.886) p ≤ 0.05) (Fig. 5b).

Fig. 5
figure 5

UPLC N-glycan changes are associated with the mRNA expression levels of glycosyltransferases. a Histogram representation of the relative expression of both the MGAT5 and ST3GAL4 glycosyltransferases and enzymes involved in the sugar nucleotide donor pathway (MPI, TSTA3 and GMDS) in 5-AZA-dC treated compared to the untreated controls. b A comparison of the expression of the two transcription factors, GATA2 and 3 with the glycosyltransferases, ST3GAL4 and MGAT5. The relative expression level of each gene was calculated according to the ddCt method normalised to TBP. Each condition was undertaken as one biological replicate with 2–3 technical replicates. Spearman’s correlation was used

Subsequent siRNA knockdown of the GATA2 and GATA3 was performed to determine the impact this might have on the expression levels of the glycosyltransferases MGAT5 and ST3GAL4 (Fig. 6, Additional file 5: Figure S7). For these experiments, the chemo-sensitive/chemo-resistant pairs A2780/A2780cis were chosen for their high expression levels of GATA2 and PEO1/PEO4 for their high expression levels of GATA3.

Fig. 6
figure 6

Expression of GATA2 and 3 affects the expression of MGAT5 and ST3GAL4. a Quantitative RT-PCR analysis of MGAT5 and ST3GAL4 and their transcription factors (GATA2 and GATA3). The relative mRNA expression level of each gene was calculated according to the ddCt method, normalised to TBP. Each condition represents two biological replicates. b Densitometric (Image J) fold change of the protein expression of GATA2, GATA3, MGAT5 and ST3GAL4 following siRNA knockdown of GATA2 or 3 (Western blots are in Additional file 5: Figure S7). Each condition represents two or three biological replicates. *p value  ≤ 0.05 or **p value ≤ 0.005 (fold expression in particular conditions was compared as follows: untreated vs 5-AZA-dC treated and both conditions with vs without siGATA2 treatments). (Two-way ANOVA and Spearman’s correlation) All p values for this figure are in Additional file 4: Table S3 including the correlation between mRNA and protein expression)

Following successful knockdown of GATA2 in the A2780/A2780cis pair, there was a significant decrease in the mRNA levels of ST3GAL4 and MGAT5 compared to non-transfected cells (Fig. 6a). For 5-AZA-dC-treated A2780 cells, successful GATA2 knockdown demonstrated a significant mRNA decrease in MGAT5 (p ≤ 0.005). 5-AZA-dC-treated A2780cis cells showed no successful GATA2 knockdown, and a significant increase in MGAT5 and ST3GAL4 was observed (p ≤ 0.05) (Fig. 6a). Following successful knockdown of GATA3 in the PEO1/PEO4 isogenic pair, there was a significant decrease in the mRNA levels of MGAT5 in PEO1 compared to non-transfected cells (p ≤ 0.05) (Fig. 6a). For 5-AZA-dC treated PEO1/PEO4 cells, there was no decrease in GATA3, and ST3GAL4 was significantly increased in PEO1 (p ≤ 0.05) (Fig. 6a).

At the protein level, MGAT5 and ST3GAL4 were significantly decreased following siRNA knock down of GATA2 and 3, but only in untreated cells (Fig. 6b).

Finally, increased mRNA levels of GATA2 were significantly associated with increases in the mRNA levels of MGAT5 and ST3GAL4. In addition, increased MGAT5 protein levels correlated with increased protein levels of ST3GAL4 in all cell lines. Moreover, increased MGAT5 gene levels correlated with increased gene levels of ST3GAL4 in A2780 and A2780cis cell lines (ρ = 0.483–0.911) (p ≤ 0.05) (Additional file 4: Table S3).

In silico survival analysis identifies increased ST3GAL4 mRNA expression as a marker of poor prognosis.

Increased ST3GAL4 mRNA levels were associated with worse progression free survival in patients with a diagnosis of ovarian cancer when separated by the median (logrank p = 0.022, hazard ratio = 0.86) (Fig. 7a) (kmplot.com). High mRNA levels of ST3GAL4 were also associated with shorter recurrence free survival in patients with a diagnosis of lymph node positive TNBC when separated by the median (logrank p = 0.022, hazard ratio = 2.01) (Fig. 7b). ST33GAL4 expression did not influence recurrence free survival in this cohort. Survival analysis of available datasets showed that MGAT5 expression did not significantly correlate with progression free survival and overall survival of patients with ovarian cancer or recurrence free survival of patients with TNBC [16].

Fig. 7
figure 7

Survival analysis indicates the poor prognostic significance of increased ST3GAL4 expression. ST3GAL4 mRNA in a ovarian cancer patients in progression free survival separated by the median and in b lymph node positive triple negative breast cancer patients in recurrence free survival when separated by the median. HR hazard ratio

Discussion

The epigenetic link between glycosylation and cancer is a relatively new concept [2]. This study looked at the impact 5-AZA-dC had on the N-glycosylation of TNBC and ovarian cancer cell lines. DNA methylation was decreased in all cells following pharmacological demethylation, more in the chemo-resistant cells from the chemo-resistant/chemo-sensitive pair. 5-AZA-dC increased resistance to cisplatin in chemo-sensitive cells while enhancing sensitivity in chemo-resistant cell lines.

There were two glycans, namely hybrid pentamanosylated bisected outer arm fucosylated biantennary digalactosylated monosialylated glycan with polylactosamine extension and biantennary digalactosylated disialylated glycan, increased on secreted glycoproteins in chemo-resistant compared to chemo-sensitive cell lines.

Although not reaching significance, there was an observed increase in highly branched, galactosylated and sialylated glycans on the secreted glycoproteins, in the chemo-sensitive ovarian cancer cell lines A2780 and PEO1 consistent with the findings of Saldova et al. [7], where the ovarian cancer cell line OVCAR3 was studied [7]. The chemo-resistant ovarian cancer and TNBC breast cancer cell lines showed decreases in these types of glycans. Chakraborty et al. [17] and Saldova et al. [7] published conflicting results in relation to the changes seen in branched glycans post-5-AZA-dC treatment, Saldova et al. [7] accredited this to the variation in cell line models used [7, 17] and confirmed by this study. The impact these glycan alterations potentially have on tumourigenesis is significant, as increased altered glycosylation such as branching and sialylation is widely implicated in immune evasion, drug resistance, metastasis and EMT [18,19,20]. The markers analysed for EMT were epithelial (E)-cadherin, N-cadherin and Vimentin. E-cadherin is a protein necessary for apical-basal polarity [21], and N-cadherin and Vimentin are required for enabling cell motility [21]. This combination of EMT markers would then identify a cadherin switch typical of cells undergoing EMT. A2780, MDA-MB-231 and MDA-MB-436 cell lines expressed high levels of N-cadherin in the absence of E-cadherin before 5-AZA treatment suggesting a marked mesenchymal phenotype. Increases in EMT were shown only in the chemo-sensitive ovarian cancer cell lines A2780 and PEO1; the same cell lines with increased branching, galactosylation and sialylation on cell secreted glycans. These results correlate with the migratory capacity of these cells, as a significant increase in the migration of A2780 and a trend in an increase in migration in PEO1 cells was seen after 5-AZA-dC treatment. E-cadherin is one of the most important molecules in cell–cell adhesion in epithelial tissues and has been found to be shed from the cell surface into the secreted fraction in ovarian cancer cells [22]. MGAT5, responsible for increased branching of N-linked glycans, when overexpressed in a gastric cancer cell line, induced an error in the trafficking of the E-cadherin protein, resulting in its mis/re-location of E-cadherin from the cell membrane to the cytoplasm, resulting in reduction in cell to cell adhesion [2]. This is consistent with the finding that increased branching of the N-linked glycans on E-cadherin results in tumours with increased mobility and metastatic potential [23].

Differential markers which are well established with other cellular events such as apoptosis and senescence were examined. Firstly, PARP cleavage, a marker of apoptosis, showed an increase in all 6 cell lines post-5-AZA-dC treatment. 5-AZA-dC causes DNA damage inducing strand breaks and this recruits PARP, which functions in the DNA damage response and regulates DNA repair, which forms covalent adducts that induce apoptosis [24]. DNA methylation downregulates PARP expression [25], therefore, our results are consistent with the increased expression of PARP following 5-AZA-dC treatment in our cell lines. A marker of cellular senescence, p21, was significantly increased in all cell lines following 5-AZA-dC treatment. P21 inhibits G1 phase [26], consistent with our cells decreasing numbers in the G1 phase, following 5-AZA-dC treatment. P21 is responsible for growth initiation, local tumour cell invasion and aggressiveness [27].

Another marker of senescence, p16, showed no significant variation in expression. P16 is inactivated in tumours and is highly expressed in senescent cells [28], Rb is activated upon senescence by p21 or p16 [28], preventing premature senescence though DNMT1. DNA methylation inactivates p21 [29] and p16 [30]. Therefore, it is consistent with our significant increase in p21 in all cell lines. 5-AZA-dC is a 5-methylcytosine nucleoside analogue which can replace cytosine after conversion to its triphosphate form. Its main mode of action though is in the trapping of DNMTs and subsequently reducing their ability to methylate DNA [31]. This process takes numerous cell replications before the effect is pronounced. Therefore, p16, which is activated later for the stabilisation of the senescence phenotype, may need a longer time period before differential levels of p16 expression would be observed [32]. DNA methylation down regulates Rb [33]; however, we saw a significant decrease in Rb in PEO1 and PEO4 cell lines post-5-AZA-dC treatment. These cell lines have numerous mutations in BRCA 1 and 2 proteins [34], and the Rb gene is frequently inactivated by gene disruption in BRCA [35]. Interestingly, A2780 has BRCA1 and 2 mutations also [34].

Expression of glycosyltransferases responsible for altered glycosylation was investigated. We have looked at the glycosyltransferases which were previously found altered and responsible for changes on secreted N-glycans in ovarian cancer [7] and found that MGAT5 and ST3GAL4 correlated with increases in highly branched and sialylated structures on secreted glycans in the chemo-sensitive ovarian cancer cell lines. Our in silico survival analysis indicated the poor prognostic significance of increased ST3GAL4 expression in ovarian and lymph node positive TNBC patients, in agreement with the fact, that increase in ST3GAL4 and associated sialylation is associated with cancer progression and chemo-resistance [2, 23, 36]. Expression of other glycosyltransferases, such as MGAT3, which increased expression was found to be responsible for a decrease in core fucosylated tetra-antennary secreted N-glycans post-5-AZA-dC treatment in HepG2 cell lines [37], could be potentially also further investigated. Our in-silico analysis identified potential TFs with binding sites in the proximal promoter regions of both MGAT5 and ST3GAL4, GATA1-3, which may regulate the expression of these glycosyltransferases. While GATA3 expression is required for normal development of the mammary gland, it has recently been reported to be overexpressed in many metastatic breast cancers [38]. Increased GATA2/3 expression is associated with a more aggressive phenotype and associated with metastatic disease [39, 40]. Both TFs are highly dependent on their methylation status for regulation [41, 42]. Knockdown of either GATA2 or GATA3, resulted in a significant decrease in MGAT5, ST3GAL4, or both, on gene expression and protein levels, as confirmed by RT-qPCR and western blots. Interestingly, following 5-AZA-dC treatment, there was only successful knockdown of GATA2 in the A2780 cell line with an associated decrease in MGAT5. In the A2780cis and PEO1 cell lines, the expression of the ST3GAL4 and MGAT5 was significantly increased after 5-AZA-dC treatment. Demethylation of GATA TFs may cause their significant increase following 5-AZA-dC, and therefore, this effect could mask the effect of the siRNA silencing. Decreases in the mRNA and protein levels of GATA2/GATA3 correlated with decreases in mRNA and protein levels of ST3GAL4/MGAT5. The GATA family of TFs contain six proteins, GATA1-6. It has been reported that aberrant DNA methylation on GATA2 affects GATA6 expression levels in gastric cancer progression [43]. 5-AZA-dC treatment may induce DNA methylation changes on GATA TFs and further contribute to disease progression. Considering only GATA1-3 was assessed, it may be possible that GATA 4–6 also have binding motifs on the glycosyltransferases and warrants further investigation.

Study limitations and future perspectives To further investigate cell-line specific patterns of glycan changes, future investigations would include whole genome methylation and transcriptome analysis post-5-AZA-dC treatment. To better characterise the roles of GATA2 and GATA3 in MGAT5 and ST3GAL4 expression, chromatin immunoprecipitation could confirm that both GATA2 and GATA3 bind to MGAT5 and ST3GAL4 genes thereby regulating their expression. The possibility that 5-AZA-dC treatment results in the hypomethylation these GATA TFs could be confirmed by bisulfite pyrosequencing. Additionally, binding of GATA may be methylation-sensitive. GATA1 TF binds to CGATA elements only if cytosine is unmethylated [44].

Based on our results, the 5-AZA-dC treatment, currently in the clinical trials, could be potentially beneficial for patients with chemo-resistant cancer but potentially harmful for chemo-sensitive cancer presentations; therefore, patients should be carefully selected for this treatment. Glycomics changes in glycosylation on serum proteins have been altered with cancer treatment and prognosis [6, 36, 45] and could be further investigated.

Conclusion

Our results show that the alterations produced by 5-AZA-dC treatment are cell line specific. This study is the first to connect increases in branching and sialylation on N-glycome from secreted glycoproteins in chemo-sensitive cell lines to increase migration post-5-AZA-dC treatment. This was directly related to the altered transcription of the glycosyltransferases ST3GAL4 and MGAT5, regulated in part by GATA2 and GATA3 TFs. Increased expression of ST3GAL4 was associated poor recurrence free survival in ovarian and lymph node positive TNBC patients. While more investigation is required, there appears to be a direct and novel link between the GATA TFs and these glycosyltransferases. 5-AZA-dC also triggers a therapeutic-induced senescence (TIS) and an EMT phenotype in the chemo-sensitive ovarian cancer cell lines A2780 and PEO1, with an associated increase in cellular migration. Based on our results, the 5-AZA-dC treatment could be potentially beneficial for patients with chemo-resistant cancer but harmful for chemo-sensitive patients. Therefore, our study is important for patient selection for this treatment.

Methods

Tissue culture and 5-AZA-2′-deoxycytidine treatment

The ovarian cancer cell lines A2780, PEO1, PEO4, and the TNBC cell lines MDA-MB-231, MDA-MB-436 were obtained and cultured as described in Greville et al. [14]. The cells were treated with 5-AZA-dC, specifically, A2780 and A2780cis cells were treated with 1 μM 5-AZA-dC, while PEO1, PEO4, MDA-MB-231 and MDA-MB-436 were treated with 0.1 μM 5-AZA-dC (Sigma)) every 24 h for 3 days (72 h). A dose response curve of 0.1, 1 and 2 μM 5-AZA-dC treatment was carried out for each of the cell lines with cell viability and DNA demethylation assessed. For A2780 and A2780cis, the optimum concentration was found to be 1 μM and for other cell lines 0.1 μM 5-AZA-dC was optimal. During the last 5-AZA-dC treatment, complete growth medium was replaced with serum free medium. Following treatment, the serum-free medium was collected for analysis of secreted glycoproteins. To investigate the potential impact of 5-AZA-2′-deoxycytidine treatment on cisplatin sensitivity, cells treated with cisplatin and combination of cisplatin and 5-AZA-dC were cultured at 37 °C with 5% CO2 for 3 days and cells were counted and compared.

Flow cytometry to assess DNA methylation status was performed as described in Greville et al. [14].

Harvesting of secreted and cell glycoproteins

N-glycans from the secreted glycoproteins

Supernatants from the cultured cells were collected and concentrated using an Amicon Ultra-15 10 K ultrafiltration (Millipore) to a final volume of ~ 200 μL. Proteins were precipitated with a half volume of 50:50 TCA: acetone (w/v) on ice. The mixture was then incubated for 45 min on ice and centrifuged at 13,000 rpm for 5 min. The resultant pellet was washed with cold acetone and centrifuged again at 13,000 rpm for 5 min.

N-glycans from the cell glycoproteins

Cultured cells were harvested and resuspended in 250μL sample buffer (2% SDS, 62.5 mM TRIS pH 6.6) and left on ice for 20 min, pipetted again and left on ice further 20 min. The whole content was then drawn up and down 10times through a 21 gauge needle and centrifuged at 13,400 rpm at 4 °C for 20 min. The resultant supernatant was removed, and the pellets were dried.

These final pellets of secreted and cell glycoproteins were dried and resuspended in sample buffer for subsequent glycan analysis.

Glycan analysis

N-glycan analyses including glycan release, labelling, HILIC-UPLC and exoglycosidase digestions were performed as described in Greville et al. [14]. Briefly, N-glycans were released from glycoproteins in samples by in situ digestion with Peptide N-glycosidaseF (PNGaseF; Prozyme) in-gel blocks, and fluorescently labelled with 2-aminobenzamide (2AB) by reductive amination. HILIC-UPLC was carried out on a BEH Glycan 1.7 μM 2.1 × 150 mm column (Waters) on an Acquity UPLC H-Class (Waters) coupled with an Acquity fluorescence detector using 30 min method and calibrated using a dextran ladder. The 2AB-labelled oligosaccharides were digested using arrays of the following exoglycosidase enzymes: Arthrobacter ureafaciens sialidase (ABS), Streptococcus pneumoniae sialidase (NAN1), bovine testes β-galactosidase (BTG), bovine kidney α-fucosidase (BKF), β-N-acetylglucosaminidase cloned from Streptococcus pneumonia, expressed in E. coli (GUH), and jack bean α-mannosidase (JBM), almond meal α-fucosidase (AMF).

Feature analysis

Glycan peaks were pooled based on similar structural or compositional features of the peak glycan members. Features pertaining to a peak were determined based on the major glycan members of that peak (Additional file 2: Table S1 for secreted glycans and Additional file 3: Table S2 for cell glycans).

Electrophoresis and western blot analysis

Proteins extracted from trypsinised cells using RIPA lysis buffer (BioRad) were separated by SDS-PAGE using a 4–15% precast TGX gels (BioRad) and transferred onto PVDF membranes using the Trans-Blot®Turbo™ system (Biorad). Blots were blocked and incubated with rabbit monoclonal antibodies targeting PARP, p21, p16, N-Cadherin, Vimentin, GATA2 and GATA3; rabbit polyclonal antibody targeting LC3, mouse monoclonal antibodies targeting E-Cadherin, retinoblastoma (Rb), ST3GAL4 and MGAT5 as described in Greville et al. [14]. The membranes were then incubated with a secondary goat anti-mouse or anti-rabbit antibody. The blots were developed using TMB for enzymatic colourimetric detection. To analyse protein loading, mouse monoclonal, α-Tubulin antibody (1:10,000; Santa Cruz, CA, USA) was used. Western blots were quantified using ImageJ™ software (FIJI).

RT-qPCR (reverse transcription quantitative PCR)

Eight gene transcripts were analysed comprising of glycosyltransferases mannosyl-(α1,6-)-glycoprotein β1,6-N-acetyl-glucosaminyltransferase (MGAT5), β-galactoside alpha-2,3-sialyltransferase 4 (ST3GAL4), GDP-mannose-4,6-dehydratase (GMDS), mannose phosphate isomerase (MPI), tissue specific transplantation antigen (TSTA3) GATA1, GATA2 and GATA3 as described in Greville et al. [14]. TBP was selected as a reference gene [46,47,48]. At least two independent technical replicates were performed for each sample. Samples were analysed by triplicate in each experiment. Results were expressed as the mean ± SD values.

Migration assay

The Oris™ Cell Migration Assay (Platypus Technologies) was used to assess migration in all untreated and treated cell lines as described in Greville et al. [14]. 2.5 × 104 (for quicker growing cells A2780, A2780cis and MB-MDA-231) or 5 × 104 cells/well/100 μL (for slower growing cells PEO1, PEO4 and MB-MDA-436) were added into stopper-loaded wells to normalise for the effect of proliferation. Proliferation controls were additionally treated with 1% SDS in water to solubilise the colour and incubated 1 hr at RT on a rocking plate. A reading of each well was taken at 595 nm.

Transient GATA knockdown

Cells were transiently transfected with 100 nM of either siRNA targeting GATA2 (A2780, A2780cis), or GATA3 (PEO1, PEO4) (Dharmacon) transcripts as described in Greville et al. [14].

Statistical analysis

All data are expressed as the means ± the standard deviation (SD). Statistical analyses were performed using SPSS statistical software for Windows (version 24.0; SPSS Inc.). For significances in glycan data, the HILIC-UPLC data were logit transformed and then used MANOVA and Tukey test. For Western blot densitometry, proliferation and migration assays, parametric T test was used. For RT-qPCR, a 2-way ANOVA was used with a post-hoc T test. p values were adjusted for multiple testing using the Benjamini–Hochberg method. The criterion for significance was set at *p value ≤ 0.05 or **p value ≤ 0.005.

Survival analysis

The prognostic significance of ST3GAL4 and MGAT5 was evaluated using the Kaplan–Meier method (kmplotter.com). Progression free survival was chosen as the outcome of interest for ovarian cancer, while recurrence free survival was chosen as the outcome of patients with lymph node positive triple negative breast cancer. Expression levels were divided above and below the median into two cohorts; high vs low expression. Differences in survival was determined by the log rank test and expressed as a logrank p value. Univariate analysis was carried out using Cox proportional hazard models, and hazard ratios (HR) and 95% confidence intervals (CI) were calculated.

Availability of data and materials

All data generated or analysed during this study are included in this published article (and its supplementary information files).

Abbreviations

5-AZA-dC:

5-AZA-2-deoxycytidine

5′MeC:

Anti-5′methylCytidine

CI:

Confidence intervals

DNMTi:

DNA methylatransferase inhibitors

EMT:

Epithelial to mesenchymal transition

ER:

Endoplasmic reticulum

GMDS:

GDP-mannose-4,6-dehydratase

HILIC:

Hydrophilic interaction liquid chromatography

HR:

Hazard ratio

GU:

Glucose unit

MGAT5:

Glycosyltransferases mannosyl-(α1,6-)-glycoprotein β1,6-N-acetyl-glucosaminyltransferase

MPI:

Mannose phosphate isomerase

OST:

Oligosaccharyltransferase

PI:

Propidium iodide

PNGase F:

Peptide N-glycosidase F

RT-qPCR:

Reverse transcription quantitative PCR

ST3GAL4:

β-Galactoside alpha-2,3-sialyltransferase 4

STR:

Short tandem repeat

TF:

Transcription factor

TNBC:

Triple negative breast cancer

TSTA3:

Tissue specific transplantation antigen

UPLC:

Ultra performance liquid chromatography

References

  1. Reddy SM, Barcenas CH, Sinha AK, Hsu L, Moulder SL, Tripathy D, et al. Long-term survival outcomes of triple-receptor negative breast cancer survivors who are disease free at 5 years and relationship with low hormone receptor positivity. Br J Cancer. 2018;118(1):17–23.

    Article  CAS  PubMed  Google Scholar 

  2. Greville G, McCann A, Rudd PM, Saldova R. Epigenetic regulation of glycosylation and the impact on chemo-resistance in breast and ovarian cancer. Epigenetics. 2016;11(12):845–57.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Steensma DP. Can hypomethylating agents provide a platform for curative therapy in myelodysplastic syndromes? Best Pract Res Clin Haematol. 2012;25(4):443–51.

    Article  CAS  PubMed  Google Scholar 

  4. Sandoval J, Esteller M. Cancer epigenomics: beyond genomics. Curr Opin Genet Dev. 2012;22(1):50–5.

    Article  CAS  PubMed  Google Scholar 

  5. Medvedeva YA, Khamis AM, Kulakovskiy IV, Ba-Alawi W, Bhuyan MSI, Kawaji H, et al. Effects of cytosine methylation on transcription factor binding sites. BMC Genomics. 2014;15:119.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  6. Marino K, Saldova R, Adamczyk B, Rudd PM. Changes in serum N-glycosylation profiles: functional significance and potential for diagnostics. In: Rauter AP, editor. Carbohydrate chemistry: chemical and biological approaches, vol. 37. Cambridge: RSC Publishing; 2012.

    Google Scholar 

  7. Saldova R, Dempsey E, Perez-Garay M, Marino K, Watson JA, Blanco-Fernandez A, et al. 5-AZA-2’-deoxycytidine induced demethylation influences N-glycosylation of secreted glycoproteins in ovarian cancer. Epigenetics. 2011;6(11):1362–72.

    Article  CAS  PubMed  Google Scholar 

  8. Chachadi VB, Cheng H, Klinkebiel D, Christman JK, Cheng P-W. 5-Aza-2′-deoxycytidine increases sialyl Lewis X on MUC1 by stimulating β-galactoside:α2,3-sialyltransferase 6 gene. Int J Biochem Cell Biol. 2011;43(4):586–93.

    Article  CAS  PubMed  Google Scholar 

  9. Saldova R, McCann A, Rudd PM. Commentary on paper: 5-Aza-2’-deoxycytidine increases sialyl Lewis X on MUC1 by stimulating beta-galactoside:alpha2,3-sialyltransferase 6 gene (Chachadi et al.). Int J Biochem Cell Biol. 2012;44(5):737.

    Article  CAS  PubMed  Google Scholar 

  10. Sidler C, Kovalchuk O, Kovalchuk I. Epigenetic regulation of cellular senescence and aging. Front Genetics. 2017;8:138.

    Article  CAS  Google Scholar 

  11. Paschall AV, Liu K. Epigenetic regulation of apoptosis and cell cycle regulatory genes in human colon carcinoma cells. Genomics Data. 2015;5:189–91.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Gwak H, Kim S, Dhanasekaran DN, Song YS. Resveratrol triggers ER stress-mediated apoptosis by disrupting N-linked glycosylation of proteins in ovarian cancer cells. Cancer Lett. 2016;371(2):347–53.

    Article  CAS  PubMed  Google Scholar 

  13. Lopez-Sambrooks C, Shrimal S, Khodier C, Flaherty DP, Rinis N, Charest JC, et al. Oligosaccharyltransferase inhibition induces senescence in RTK-driven tumor cells. Nat Chem Biol. 2016;12:1023.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Greville G, Llop E, Huang C, Creagh-Flynn J, Pfister S, O'Flaherty R, et al. Hypoxia alters epigenetic and N-glycosylation profiles of ovarian and breast cancer cell lines in-vitro. Front Oncol. 2020 (in press).

  15. Steinestel K, Eder S, Schrader AJ, Steinestel J. Clinical significance of epithelial-mesenchymal transition. Clin Transl Med. 2014;3:17.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Gyorffy B, Lanczky A, Eklund AC, Denkert C, Budczies J, Li Q, et al. An online survival analysis tool to rapidly assess the effect of 22,277 genes on breast cancer prognosis using microarray data of 1,809 patients. Breast Cancer Res Treat. 2010;123(3):725–31.

    Article  PubMed  CAS  Google Scholar 

  17. Chakraborty AK, de Frietas SJ, Chakraborty D, Funasaka Y, Bhattacharya M, Chatterjee A, et al. GnT-V expression and metastatic phenotypes in macrophage–melanoma fusion hybrids is down-regulated by 5-Aza-dC: evidence for methylation sensitive, extragenic regulation of GnT-V transcription. Gene. 2006;374:166–73.

    Article  CAS  PubMed  Google Scholar 

  18. Jandus C, Boligan KF, Chijioke O, Liu H, Dahlhaus M, Demoulins T, et al. Interactions between Siglec-7/9 receptors and ligands influence NK cell-dependent tumor immunosurveillance. J Clin Invest. 2014;124(4):1810–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Lesniak D, Sabri S, Xu Y, Graham K, Bhatnagar P, Suresh M, et al. Spontaneous epithelial-mesenchymal transition and resistance to HER-2-targeted therapies in HER-2-positive luminal breast cancer. PLoS ONE. 2013;8:e71987.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Nangia-Makker P, Balan V, Raz A. Galectin-3-binding and metastasis. Methods Mol Biol (Clifton, NJ). 2012;878:251–66.

    Article  CAS  Google Scholar 

  21. Kalluri R, Weinberg RA. The basics of epithelial-mesenchymal transition. J Clin Investig. 2009;119(6):1420–8.

    Article  CAS  PubMed  Google Scholar 

  22. Faca VM, Ventura AP, Fitzgibbon MP, Pereira-Faca SR, Pitteri SJ, Green AE, et al. Proteomic analysis of ovarian cancer cells reveals dynamic processes of protein secretion and shedding of extra-cellular domains. PLoS ONE. 2008;3(6):e2425.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. da Fonseca LM, da Silva VA, Freire-de-Lima L, Previato JO, Mendonça-Previato L, Capella MAM. Glycosylation in cancer: interplay between multidrug resistance and epithelial-to-mesenchymal transition? Front Oncol. 2016;6:158.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Orta ML, Calderon-Montano JM, Dominguez I, Pastor N, Burgos-Moron E, Lopez-Lazaro M, et al. 5-Aza-2’-deoxycytidine causes replication lesions that require Fanconi anemia-dependent homologous recombination for repair. Nucleic Acids Res. 2013;41(11):5827–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Gong C, Tao G, Yang L, Liu J, Liu Q, Li W, et al. Methylation of PARP-1 promoter involved in the regulation of nano-SiO2-induced decrease of PARP-1 mRNA expression. Toxicol Lett. 2012;209(3):264–9.

    Article  CAS  PubMed  Google Scholar 

  26. Ferrandina G, Stoler A, Fagotti A, Fanfani F, Sacco R, De Pasqua A, et al. p21WAF1/CIP1 protein expression in primary ovarian cancer. Int J Oncol. 2000;17(6):1231–5.

    CAS  PubMed  Google Scholar 

  27. Dai M, Al-Odaini AA, Fils-Aime N, Villatoro MA, Guo J, Arakelian A, et al. Cyclin D1 cooperates with p21 to regulate TGFbeta-mediated breast cancer cell migration and tumor local invasion. Breast Cancer Res BCR. 2013;15(3):R49.

    Article  CAS  PubMed  Google Scholar 

  28. Ben-Porath I, Weinberg RA. The signals and pathways activating cellular senescence. Int J Biochem Cell Biol. 2005;37(5):961–76.

    Article  CAS  PubMed  Google Scholar 

  29. Fang JY, Lu YY. Effects of histone acetylation and DNA methylation on p21( WAF1) regulation. World J Gastroenterol. 2002;8(3):400–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Shibata KR, Aoyama T, Shima Y, Fukiage K, Otsuka S, Furu M, et al. Expression of the p16INK4A gene is associated closely with senescence of human mesenchymal stem cells and is potentially silenced by DNA methylation during in vitro expansion. Stem cells. 2007;25(9):2371–82.

    Article  CAS  PubMed  Google Scholar 

  31. Gowher H, Jeltsch A. Mechanism of inhibition of DNA methyltransferases by cytidine analogs in cancer therapy. Cancer Biol Ther. 2004;3(11):1062–8.

    Article  CAS  PubMed  Google Scholar 

  32. Campisi J, d’AddadiFagagna F. Cellular senescence: when bad things happen to good cells. Nat Rev Mol Cell Biol. 2007;8:729.

    Article  CAS  PubMed  Google Scholar 

  33. Mileo AM, Mattarocci S, Matarrese P, Anticoli S, Abbruzzese C, Catone S, et al. Hepatitis C virus core protein modulates pRb2/p130 expression in human hepatocellular carcinoma cell lines through promoter methylation. J Exp Clin Cancer Res CR. 2015;34:140.

    Article  PubMed  CAS  Google Scholar 

  34. Stordal B, Timms K, Farrelly A, Gallagher D, Busschots S, Renaud M, et al. BRCA1/2 mutation analysis in 41 ovarian cell lines reveals only one functionally deleterious BRCA1 mutation. Mol Oncol. 2013;7(3):567–79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Jonsson G, Staaf J, Vallon-Christersson J, Ringner M, Gruvberger-Saal SK, Saal LH, et al. The retinoblastoma gene undergoes rearrangements in BRCA1-deficient basal-like breast cancer. Cancer Res. 2012;72(16):4028–36.

    Article  PubMed  CAS  Google Scholar 

  36. Saldova R, Piccard H, Perez-Garay M, Harvey DJ, Struwe WB, Galligan MC, et al. Increase in sialylation and branching in the mouse serum N-glycome correlates with inflammation and ovarian tumour progression. PLoS ONE. 2013;8(8):e71159.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Klasic M, Kristic J, Korac P, Horvat T, Markulin D, Vojta A, et al. DNA hypomethylation upregulates expression of the MGAT3 gene in HepG2 cells and leads to changes in N-glycosylation of secreted glycoproteins. Sci Rep. 2016;6:24363.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. De Lara S, Parris TZ, Werner Ronnerman E, Helou K, Kovacs A. GATA3 as a putative marker of breast cancer metastasis-A retrospective immunohistochemical study. Breast J. 2018;24(2):184–8.

    Article  PubMed  CAS  Google Scholar 

  39. Wang Y, He X, Ngeow J, Eng C. GATA2 negatively regulates PTEN by preventing nuclear translocation of androgen receptor and by androgen-independent suppression of PTEN transcription in breast cancer. Hum Mol Genet. 2012;21(3):569–76.

    Article  CAS  PubMed  Google Scholar 

  40. Chen HJ, Huang RL, Liew PL, Su PH, Chen LY, Weng YC, et al. GATA3 as a master regulator and therapeutic target in ovarian high-grade serous carcinoma stem cells. Int J Cancer. 2018;143:3106–19.

    Article  CAS  PubMed  Google Scholar 

  41. Rodriguez-Bravo V, Carceles-Cordon M, Hoshida Y, Cordon-Cardo C, Galsky MD, Domingo-Domenech J. The role of GATA2 in lethal prostate cancer aggressiveness. Nat Rev Urol. 2017;14(1):38–48.

    Article  CAS  PubMed  Google Scholar 

  42. Zhu Y, Qiu Y, Yu H, Yi S, Su W, Kijlstra A, et al. Aberrant DNA methylation of GATA binding protein 3 (GATA3), interleukin-4 (IL-4), and transforming growth factor-β (TGF-β) promoters in Behcet’s disease. Oncotarget. 2017;8(38):64263–72.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Song SH, Jeon MS, Nam JW, Kang JK, Lee YJ, Kang JY, et al. Aberrant GATA2 epigenetic dysregulation induces a GATA2/GATA6 switch in human gastric cancer. Oncogene. 2018;37(8):993–1004.

    Article  CAS  PubMed  Google Scholar 

  44. Yang L, Chen Z, Stout ES, Delerue F, Ittner LM, Wilkins MR, et al. Methylation of a CGATA element inhibits binding and regulation by GATA-1. Nat Commun. 2020;11(1):2560.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Saldova R, Haakensen VD, Rodland E, Walsh I, Stockmann H, Engebraaten O, et al. Serum N-glycome alterations in breast cancer during multimodal treatment and follow-up. Mol Oncol. 2017;11(10):1361–79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Shaw RJ, Omar MM, Rokadiya S, Kogera FA, Lowe D, Hall GL, et al. Cytoglobin is upregulated by tumour hypoxia and silenced by promoter hypermethylation in head and neck cancer. Br J Cancer. 2009;101(1):139–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Fink T, Lund P, Pilgaard L, Rasmussen JG, Duroux M, Zachar V. Instability of standard PCR reference genes in adipose-derived stem cells during propagation, differentiation and hypoxic exposure. BMC Mol Biol. 2008;9:98.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  48. Rasmussen LM, Frederiksen KS, Din N, Galsgaard E, Christensen L, Berchtold MW, et al. Prolactin and oestrogen synergistically regulate gene expression and proliferation of breast cancer cells. Endocr Relat Cancer. 2010;17(3):809–22.

    Article  CAS  PubMed  Google Scholar 

Download references

Funding

Authors acknowledge funding from the Science foundation Ireland Starting Investigator Research grant (SFI SIRG) under grant number 13/SIRG/2164, Spanish Ministry of Economy and Competitivity (grant BIO 2015–66356-R), the University of Girona (grant MPCUdG2016/028) and the University College Dublin (UCD PhD Advance Core Scheme).

Author information

Authors and Affiliations

Authors

Contributions

PMR, AMC and RS conceived and initiated this study. GG, EL, JH, ASP and RS carried out the experiments. GG, EL, ASP and RS generated and interpreted the results. EL, RP, SM and RS conducted the statistical analyses. GG, EL, RP, PMR, AMC and RS drew the conclusions. GG and RS with significant input from AMC wrote the manuscript which was reviewed by all co-authors. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Radka Saldova.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1:

Supporting data: Glycosylation gene promoter CpG island analysis

Additional file 2:

 Table S1. N-glycans and features from secreted glycoproteins from all cell lines.

Additional file 3:

 Table S2. Detailed assignments of N-glycans from cell glycoproteins of A2780 cell line using UPLC and exoglycosidase digestions.

Additional file 4:

 Table S3. P-values for siRNA knock down experiments (Fig. 6)

Additional file 5:

 Figure S1: 5-AZA-dC treatment increases resistance to cisplatin in chemo-sensitive and decreases in chemo-resistant cell lines. % viability of A2780 and A2780cis cells treated with 1 μM 5-AZA-dC and 1 μM 5-AZA-dC in combination with 1 μM cisplatin. * = P-value ≤ 0.05 or ** = P-value ≤ 0.005 (T-test). Figure S2: Glycosylation changes of cell glycans on A2780 and A2780cis chemo-sensitive/chemo-resistant pair. (A) Representative UPLC chromatograms produced from secreted N-glycans of ovarian chemo-sensitive- chemo-resistant pair (A2780, A2780cis) and their separation into 39 peaks. (B) Plotted peak areas from the cell N-glycans of these cell lines. The glycans in each peak (GP1-GP39) and features are listed in Table S2. Significant changes (p < 0.05) are starred: * = P-value ≤ 0.05 or ** = P-value ≤ 0.005. (MANOVA). Heatmap histograms indicating fold changes in 5-AZA-dC treated compared to untreated cells were created using Hierarchial Clustering Explorer HCE 3.5 software. Blue indicates decreases, and red indicates increases. The shade of colour corresponds to amounts of the decreases/increases. Figure S3: Secreted glycans differ in chemo-resistant comparing to chemo-sensitive cell lines. Plotted peak areas of GP13 and 14 from the secreted N-glycans of chemo-sensitive and chemo-resistant ovarian cancer cell lines. Significant changes (p < 0.05) are starred: * = P-value ≤ 0.05 or ** = P-value ≤ 0.005. (MANOVA). Figure S4. Representative Western blots of EMT markers in the 4 ovarian cancer cell lines and 2 TNBC cell lines post-5-AZA-dC treatment (T) compared to non-treated (UT) controls. Figure S5: Migration results are not attributable to proliferation. Migration (A) and proliferation (B) or the representative 1 uM 5-AZA-dC treated relative to the untreated A2780cis cells. Figure S6. Western blot analyses of the senescence markers p16, p21 and Rb and of cellular apoptosis (PARP cleavage) markers, post-5-AZA-dC treatment. Representative Western blots of the senescence associated proteins p21, p16, Rb and the apoptosis marker PARP of 4 ovarian cancer cell lines and 2 TNBC cell lines in 5-AZA-dC treated (T) compared to untreated (UT) controls. Figure S7. SiRNA knockdown of GATA2/3, impacts on the protein expression levels of MGAT5 and ST3GAL4. Western blot analyses of siRNA knockdown of GATA2/3 in 4 ovarian cancer cell lines (A2780, A2780cis, PEO1 and PEO4). GATA2 was knocked down in the A2780 and A2780cis cell lines and GATA3 in the PEO1 and PEO4 cell lines. Results are shown comparing (i) untreated controls (UT) with siRNA GATA knockdown UT (siG2/3) (left) with (ii) 5-AZA-dC treated (T) with siRNA GATA knockdown treated T (siG + T). All results are in triplicates (n = 3, biological replicates).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Greville, G., Llop, E., Howard, J. et al. 5-AZA-dC induces epigenetic changes associated with modified glycosylation of secreted glycoproteins and increased EMT and migration in chemo-sensitive cancer cells. Clin Epigenet 13, 34 (2021). https://doi.org/10.1186/s13148-021-01015-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13148-021-01015-7

Keywords