- Research article
- Open Access
Decreased methylation in the SNAI2 and ADAM23 genes associated with de-differentiation and haematogenous dissemination in breast cancers
BMC Cancer volume 18, Article number: 875 (2018)
In breast cancer (BC), deregulation of DNA methylation leads to aberrant expressions and functions of key regulatory genes. In our study, we investigated the relationship between the methylation profiles of genes associated with cancer invasivity and clinico-pathological parameters. In detail, we studied differences in the methylation levels between BC patients with haematogenous and lymphogenous cancer dissemination.
We analysed samples of primary tumours (PTs), lymph node metastases (LNMs) and peripheral blood cells (PBCs) from 59 patients with sporadic disseminated BC. Evaluation of the DNA methylation levels of six genes related to invasivity, ADAM23, uPA, CXCL12, TWIST1, SNAI1 and SNAI2, was performed by pyrosequencing.
Among the cancer-specific methylated genes, we found lower methylation levels of the SNAI2 gene in histologic grade 3 tumours (OR = 0.61; 95% CI, 0.39–0.97; P = 0.038) than in fully or moderately differentiated cancers. We also evaluated the methylation profiles in patients with different cancer cell dissemination statuses (positivity for circulating tumour cells (CTCs) and/or LNMs). We detected the significant association between reduced DNA methylation of ADAM23 in PTs and presence of CTCs in the peripheral blood of patients (OR = 0.45; 95% CI, 0.23–0.90; P = 0.023).
The relationships between the decreased methylation levels of the SNAI2 and ADAM23 genes and cancer de-differentiation and haematogenous dissemination, respectively, indicate novel functions of those genes in the invasive processes. After experimental validation of the association between the lower values of SNAI2 and ADAM23 methylation and clinical features of aggressive BCs, these methylation profiles could improve the management of metastatic disease.
Breast cancer (BC) is one of the most common tumours occurring in women worldwide, with high mortality rates . In recent population study, it was observed that of 13,785 BC patients undergoing adjuvant chemotherapy, 24.4% developed metastatic disease, with median survival of 20 months. In contrast to previous studies, no general improvement of survival in metastatic recurrent BC patients after chemotherapy has been evidenced over the last three decades . Therefore, to assist in the development of effective diagnostic and therapeutic strategies for BC, the identification of more precise biomarkers is necessary.
It is generally accepted that during the initiation and progression of cancer, deregulation of epigenetic processes, including DNA methylation, occurs and results in aberrant expression and function of a number of key regulatory genes [3, 4]. Changes in DNA methylation in primary tumours (PTs) compared with normal breast tissues have been partially identified by numerous researchers; however, whole-genome bisulfite sequencing technology is now available, which allows the comprehensive analysis of normal and BC methylomes .
The metastatic cascade involves multiple steps enabling the detachment of cancer cells from PTs, their intravasation to the circulatory system, their avoidance of the immune reaction and their migration through the circulation. Some tumour cells can extravasate from vessels and invade distant organs. In a new microenvironment, tumour cells can persist in a dormant stage or form micrometastases, which can later develop into macrometastases [6,7,8,9].
The critical processes underlying the spread of cancer cells include the degradation of the basement membrane (BM) and the extracellular matrix (ECM). The BM is formed by dense interconnected sheets of ECM proteins, and the loss of tissue architecture and the disruption of structural boundaries between the epithelium and the stroma are the typical features of advanced stages of malignancy . ADAMs (a disintegrin and metalloproteases) have been defined as transmembrane proteins with protease and adhesion functions . ADAM23 belongs to the ADAM protein family although it does not have metalloprotease activity; however, its disintegrin-like domain facilitates its interaction with integrins, resulting in integrin-dependent cell adhesion [12,13,14]. Cancer cells not expressing ADAM23 have higher migration capacities in cell line models . Moreover, promoter hypermethylation of the ADAM23 gene in breast PTs is significantly associated with more advanced grade, higher proliferation of cancer cells, poor prognosis and lower survival rates [15,16,17].
Cancer cells contribute to the degradation of the ECM by secreting various proteolytic enzymes, such as urokinase plasminogen activator (uPA). uPA is a serine protease that contributes to the conversion of inactive plasminogen to plasmin and that, in cooperation with the uPA receptor (uPAR), degrades the fundamental ECM components, such as laminin and fibronectin . Several in vitro and in vivo studies have documented the essential role of uPA in tumour growth, the detachment of cancer cells from PT and haematogenous dissemination [19, 20]. Moreover, the prognostic impact of uPA and plasminogen activator inhibitor-1 (PAI-1) expression levels has been identified predominantly in patients with HER2+ tumours . Experimental studies have documented that promoter methylation is the dominant mechanism of uPA silencing in BC progression [22, 23]. Previously, the complete demethylation of the uPA gene and expression of uPA mRNA was found only in highly invasive oestrogen receptor-negative BC cell lines; in normal and low invasivity lines, uPA hypermethylation with reduced mRNA expression was observed . In BC patients, uPA hypomethylation has been associated with metastatic disease, but the uPA hypermethylation in non-selected BC was evidenced using the more complex methylation analyses of multiple CpG islands [24, 25].
Invasion by cancer cells occurs at the interface between tumour and host tissue, and the intensive exchange of many enzymes, chemokines and cytokines contributes to the destruction of the BM, the local modulation of the ECM and the stimulation of cancer cell migration . Chemokine CXC motif ligand 12 (CXCL12), which is produced at high levels by cells of the tumour microenvironment and target tissues including lymph nodes, interacts with its signal receptor, CXCR4, expressed on the surface of cancer cells. Using the chemotactic gradient of CXCL12-CXCR4, CXCL12 participates in the regulation of the metastatic process [27, 28]. Both encoding genes, CXCL12 and CXCR4, are epigenetically regulated during the modulation of metastatic potential in BC. In patients with histologically advanced disease and worse survival, CXCL12 methylation and CXCR4 overexpression without methylation have been found [29,30,31]. In our previous study we documented the relationship between CXCL12 hypermethylation and increased risk for the development of LNMs in BC patients .
For the infiltration of epithelial cancer cells into circulation, it is necessary for them to shift into a more mesenchymal phenotype through the epithelial-mesenchymal transition (EMT). In normal physiological development, EMT is active in embryonal and organ development. In epithelial tumours, EMT contributes to the loss of cancer cell polarity, cell-cell contacts and cell-BM interactions. Owing to these mesenchymal characteristics, cancer cells obtain the capacity to migrate and invade, and they became be more resistant to apoptosis . The EMT process is regulated by three main families of transcription factors (TFs), Zeb, Snail and Twist. Of these TFs, ZEB1 and ZEB2 (zinc finger E-box binding homeobox 1/2), SNAI1 and SNAI2 (snail family zinc finger 1/2) and TWIST1 and TWIST2 (basic helix-loop-helix transcription factor 1/2) are known as direct regulators of the adhesion molecule E-cadherin . Using immunohistochemical analyses, changes in Twist, Snail and Slug expressions have been observed in BC. In patients with metastatic disease, increasing levels of TWIST protein are associated with cancer-caused death, and higher levels of SLUG (SNAI2) protein are present in cases with higher tumour grades . More prevalent TWIST1 hypermethylation has been observed in BC compared to normal breast tissue; however, no direct evidence of the epigenetic regulation of TWIST1 mRNA or the relevant protein expression has been found [35,36,37]. Higher TWIST1 methylation levels have been evidenced in non-triple negative (TN) BC than in TN BC and in ductal compared to lobular invasive carcinomas [38, 39]. The dynamic changes in DNA methylation in genes encoding other two EMT-TFs, SNAI1 and SNAI2, have been observed in a EMT/MET cell-line model, where increased methylation was associated with decreased transcription in both genes . In a previous in vitro study, regardless of the significantly higher expression levels of both SNAI2 and ZEB2 genes in metastatic BC cell lines, only moderate DNA methylation differences were found in SNAI2 between highly and poorly metastatic lines, which was in contrast to the ZEB2 gene, which was markedly hypomethylated in highly metastatic cell lines . These results indicate more complex interactions in the epigenetic regulation of cancer invasivity including EMT.
Cancer cells disseminate from PTs either via blood vessels after intensive neo-vascularisation (haematogenous spread) or via the lymphatic system, named lymphogenous spread, which occurs after neo-lymphangiogenesis [42, 43]. In cancers, the active growth of new blood vessels from pre-existing ones is triggered by tumour cells expressing vascular endothelial growth factor A (VEGFA) in the hypoxic microenvironment of tumour tissue. New angiogenic blood vessels are abnormal, with a lack of pericytes and BM connection and with some fenestrations, allowing intravasation and the haematogenous spread of cancer cells . In BC, it has been found that endothelial cells produce soluble factors, such as hepatocyte growth factor (HGF), which contribute to EMT . Several studies have shown that in numerous cancer types including BC, the neo-vascularization associated with the presence of circulating tumour cells (CTCs) in blood vessels and the higher number of CTCs are correlated with poor prognosis and decreased overall survival [46,47,48,49]. On the other hand, tumour cells can express and secrete several lymphangiogenic factors, including VEGFC/D, to promote the formation of new lymphatic vessels in the tumour stroma . From the lymphatic vasculature, cancer cells can enter the blood circulation via the thoracic duct or form metastases in lymph nodes. Moreover, cancer cell motility and spread are affected by the density of blood and lymphatic vessels, the interstitial blood pressure, tumour hypoxia, regional lymph node metastases (LMNs) and the recently discovered primo vascular system . The preferred dissemination route of cancer cells depends on cancer types and microenvironmental conditions, but it has already been reported that carcinomas, such as BC, and melanomas more often spread through the lymphatic vessels and develop LNMs, in comparison to soft tissue sarcomas, which predominantly result in systematic metastasis [51, 52].
In our previous study of methylation profiles in 206 BC patients we showed that the risk for LNMs development and higher proliferation of cancer cells was increased by hypermethylation of CXCL12 and ADAM23 genes, respectively . For present study, we selected 59 BC patients with haematogenous and/or lymphogenous cancer dissemination to investigate more deeply the relationship between the DNA methylation patterns and the clinico-pathological parameters of more aggressive BC. In addition to CXCL12, ADAM23 and uPA genes contributing to the inhibition of cell adhesion, BM degradation and ECM remodelling, TWIST1, SNAI1 and SNAI2 genes, the protein products of which are known as key inducers of EMT, were involved. The changes in methylation levels of patients with haematogenous and/or lymphogenous dissemination of BC could help to identify potential biomarkers for metastatic potential of BC alone and also for different routes of metastatic progression.
Patients and samples
In the present study we analysed PT tissues and LNM samples, both isolated from formalin fixed paraffin-embedded tissues (FFPE) and peripheral blood cells (PBCs) from 59 selected patients with disseminated BC. Their tumours were considered as sporadic, because these patients did not fulfil criteria for hereditary BC testing. PBC samples from 53 healthy women and 10 normal mammary gland samples from mammoplasties were used as controls. PT and LNM specimens were collected at the National Cancer Institute in Bratislava, Slovakia between 2012 and 2014, and the control samples were collected in previous studies. This study was approved by the Institutional Review Board of the National Cancer Institute of Slovakia, and written informed consent was obtained from all patients and controls. The clinical and histopathological characteristics were obtained from patient records, and the tumours were defined according to TNM classification. The clinico-histopathological data of evaluated patients are summarized in Table 1 . No patient underwent preoperative radiotherapy or chemotherapy before sample collection. The age of the women with BC ranged from 37 to 79 years (mean 60.19 ± 10.38 years). The PBC samples were obtained from control women aged 40 to 85 years (mean 58.17 ± 10.50 years). The age of the normal breast tissue donors ranged from 31 to 56 years (mean 46.90 ± 9.06 years). No control persons had signs or symptoms of cancer or other serious diseases.
DNA extraction and sodium bisulfite DNA modification
The representative samples of PTs (minimally 80% of cancer cells), LNMs and healthy tissues from mammary glands were selected from FFPE blocks. Tissue samples identified under light microscope were marked and 3 mm diameter cores of tissue were removed from donor blocks into the recipient master block using a tissue microarray method, as described previously . DNA was extracted by the MagneSil Genomic Fixed Tissue System (Promega, Madison, Wisconsin, USA). DNA from PBCs was isolated using the FlexiGene DNA Kit (Qiagen, Hilden, Germany). All extraction methods were performed according to the manufacturer’s instructions. A NanoDrop 1000 spectrophotometer (Thermo Fisher Scientific, Bremen, Germany), was used to measure DNA concentration. DNA bisulfite modification was performed by the EpiTect Bisulfite Kit (Qiagen, Hilden, Germany) and the CpGenome DNA Modification Kit (Chemicon, Temecula, California, USA) for FFPEs (2 μg of isolated DNA) and PBC samples (1 μg of isolated DNA), respectively. The samples of modified DNA were divided into aliquots and stored at − 18 °C. The principle of DNA bisulfite modification is the conversion of unmethylated cytosines to uracils, while methylated cytosines remain unchanged.
PCR and methylation analyses
For the evaluation of the DNA methylation profiles in the six selected cancer-associated genes (ADAM23, uPA, CXCL12, TWIST1, SNAI1 and SNAI2), the quantitative pyrosequencing method was used. Between 3 and 8 CpGs were analysed in each gene in the CpG islands of the promoter regions flanking the transcription start site. The primers for PCR and the pyrosequencing reactions were designed by PyroMark Assay Design software 2.0 (Qiagen, Hilden, Germany). The primers are shown in Table 2. All the designed assays were validated according to the manufacturer’s instructions. The PCR reactions were performed by the PyroMark PCR Kit (Qiagen, Hilden, Germany) following manufacturer’s instructions. The concentrations of the PCR primers were 0.32 μM for uPA, CXCL12, TWIST1, SNAI1, and SNAI2 and 0.40 μM for ADAM23, with the annealing conditions of 56 °C for 30 s except for the ADAM23 gene, for which 52 °C for 30 s was used. For pyrosequencing, the PyroGold Reagent Kit (Qiagen, Hilden, Germany) and PyroMark Q24 System were used, and the results were evaluated by PyroMark Q24 2.0.6 software (Qiagen, Hilden, Germany). Methylation levels were defined as the percentage of average methylation in all CpG sites in each individual gene.
CTCs detection in peripheral blood
CTCs were detected in 5 mL of peripheral blood depleted of CD45-positive (CD45+) cells for CTC enrichment using a quantitative real-time polymerase chain reaction (qRT-PCR) assay, as described previously [55, 56]. Patient samples with higher levels of KRT19 gene transcripts than those of healthy donors were scored as epithelial CTC-positive, while patient samples with higher mRNA levels of EMT-transcription factors (TWIST1, SNAI1, SLUG (SNAI2) and ZEB1) than those of healthy donors were scored as mesenchymal CTC-positive. The expression of at least one of the markers (either epithelial or mesenchymal) at levels above the defined cut off was sufficient to define a sample as CTC-positive.
The highest expression levels of the KRT19 and EMT-inducing TF gene transcripts relative to that of GAPDH were 3.4 × 10− 3 (median 2.8 × 10− 6, range: 0–3.4 × 10− 3) for KRT19, 7.5 × 10− 4 (median 0, range: 0–7.5 × 10− 4) for TWIST1, 3.8 × 10− 2 (median 3.1 × 10− 3, range: 5.0 × 10− 4 - 3.8 × 10− 2) for SNAIL1 and 1.7 × 10− 1 (median 1.4 × 10− 2, range: 2.2 × 10− 3 - 1.7 × 10− 1) for ZEB1, while SLUG transcripts were not detected in any of the samples from healthy donor. These highest expression values in healthy donors were used as “cutoff” to determine CTCs positivity .
SPSS statistics 15.0 software was used for the statistical analyses of the data. The normality of distribution was tested by the Shapiro-Wilk test. Normally distributed samples were tested by Student’s t-test or analysis of variance (ANOVA) with Bonferroni’s corrections. Nonparametric Mann-Whitney U or Kruskal-Wallis H tests were used for non-normally distributed data. According to the normality of data, Pearson’s or Spearman’s correlations were used. Univariate analyses were performed for categorical variables using χ2 or Fisher’s exact test. A statistical two-tailed significance was regarded as P value < 0.05. A logistic regression adjusted for age was used to determine the effect of independent categorical variables and DNA methylation on the clinical status and tumour dissemination.
DNA methylation in various sample types in BC patients and controls
In this study, we evaluated specific DNA methylation levels by the quantitative pyrosequencing method in three types of samples isolated from 59 BC patients, namely, in PT and LNM tissues and PBCs. From the six analysed genes, higher DNA methylation levels were found in PTs for ADAM23, uPA, CXCL12 and TWIST1 genes, with values of 9.98%, 14.08%, 13.41% and 21.88%, respectively. Similar methylation levels of these four genes were found in LMN samples, and a positive correlation between PT and LNM tissues was observed (r = 0.598, P = < 0.001; r = 0.498, P = 0.011; r = 0.380, P = 0.013; r = 0.428, P = 0.026; respectively). In BC patients, increased methylation levels were detected in PT and LNM samples compared with PBCs for all studied genes. Methylation results in different types of samples are summarized in Table 3.
Furthermore, in univariate analysis, we found significantly different methylation levels in PBC samples between BC patients and controls for ADAM23, uPA, SNAI1 and SNAI2, with P values ≤0.001 for all. Similarly, significant differences in DNA methylation between PT and normal breast tissues were observed for CXCL12, TWIST1 and SNAI2 genes (P = 0.001, P = 0.006 and P = 0.009, respectively). However, these results cannot be correctly interpreted for ageing processes interference with DNA methylation in patient’s tumour tissue samples.
DNA methylation and clinico-histopathological characteristics
In the group of 59 BC patients, the relationships between the DNA methylation levels of all six evaluated genes and the clinico-histopathological features such as tumour histology, tumour size, LNM status, TNM, grade (G), hormonal receptor (HR) and HER2 status, Ki-67 proliferative index, tumour multifocality and CTC presence were analysed in PT, LNM and PBC samples. We identified significant differences in methylation levels between patients with G1 and G2 BC compared to those with G3 BC. Compared to patients with G3 BC, patients with lower grades presented higher methylation levels in three genes: CXCL12 in PBCs, SNAI1 in LNMs and SNAI2 in PT samples. The mean methylation levels in G1 and G2 compared to G3 for those three genes were 3.06% vs. 2.35%, 5.39% vs. 3.42%, and 7.38% vs. 4.00% with P values of 0.021, 0.007, and 0.003, respectively (Fig. 1). Methylation results for all six genes in patients with different histological grade are summarized in Additional file 1.
The relationship between grade and SNAI2 methylation in tumour tissue was confirmed by multivariate logistic regression analysis. The influence of the methylation of EMT-associated genes (TWIST1, SNAI1 and SNAI2) and the clinico-histopathological features of invasive BC, namely, age, HR and HER2 status, tumour histology and Ki-67 proliferative index on tumour grade were assessed in a multivariate model. The methylation of the SNAI2 gene negatively associated with higher tumour grade (grade 3) (OR = 0.61; 95% CI, 0.39–0.97; P = 0.038) and higher grade positively associated with the high Ki-67 proliferative index (OR = 17.34; 95% CI, 2.89–103.88; P = 0.002) (Table 4).
DNA methylation and dissemination of tumour cells
In the present study, we investigated the differences in methylation profiles between patients with haematogenous dissemination of cancer cells, indicated by the presence of CTCs in peripheral blood, and/or lymphogenous spread, represented by the presence of LNMs.
To evaluate the association between DNA methylation and haematogenous spread, we analysed DNA methylation levels in patients with epithelial CTCs (epi CTC+) and in those undergoing EMT process, as indicated by the presence of mesenchymal CTCs (mes CTC+) in their blood. We observed significantly lower methylation levels for ADAM23, TWIST1 and SNAI2 in the PTs of patients with epi CTC+ than in the PTs of patients with epi CTC- (4.45% vs. 11.46%, P = 0.003; 13.09% vs. 24.36%, P = 0.045 and 4.31% vs. 7.02%, P = 0.026, respectively) (Fig. 2). In BC patients with mes CTC+ in circulation, significant lower methylation levels of the ADAM23 gene were identified in PT and LNM tissues than in the tissues of patients with mes CTC- (4.79% vs. 11.89%, P = 0.023 for PTs and 4.29% vs. 11.31%, P = 0.024 for LNMs) (Fig. 3). Additionally, we analysed methylation profiles in three groups of patients with epi CTC+, mes CTC+ and without any CTC (CTC-) and we found the significant differences in methylation levels in SNAI1 gene between epi CTC+ and mes CTC+ patients, in SNAI2 gene between epi CTC+ and CTC- patients. For ADAM23 gene the similar values were observed between patients with epi CTC+ and mes CTC+, where DNA methylation was significantly lower for both these groups when compared to CTC- patients (Table 5).
We also evaluated DNA methylation in three groups of patients with different statuses of cancer cell dissemination, namely, CTC-negative patients with LNM (CTC-LNM+, n = 29), CTC-positive patients without LNM (CTC + LNM-, n = 16) and those positive for both parameters (CTC + LNM+, n = 13). For those analyses, one BC patient with micrometastasis in lymph node (LN) was excluded. We found significant differences in methylation levels among these three groups for ADAM23 in PT and LNM tissues (Fig. 4) and TWIST1 in PBCs.
The analysis of methylation levels in PTs and PBC samples between patients with exclusively haematogenous dissemination (CTC + LNM-) and those with only lymphogenous spread (CTC-LNM+) identified significant differences between methylation levels of the ADAM23 gene in PTs and the CXCL12 and TWIST1 genes in PBC samples (4.69% vs. 13.39%, P = 0.002; 2.31% vs. 3.03%, P = 0.036; 3.31% vs. 4.48%, P = 0.008). Using multivariate logistic regression, we confirmed the inverse association between the ADAM23 methylation level in PTs and the presence of CTCs in the peripheral blood of patients (OR = 0.45; 95% CI, 0.23–0.90; P = 0.023) (Table 6).
The high mortality rate associated with metastatic disease is a serious medical problem. To improve screening strategies and preventative treatment of metastatic progression in BC patients, detailed knowledge is needed about the molecular mechanisms that result in aggressive tumour phenotypes, including invasiveness and metastasis. In addition to many other molecular markers, DNA methylation biomarkers, particularly the specific alterations in methylation, offer tools for early detection, diagnosis and more effective management of advanced BC.
In the present study, we investigated promoter methylation in six genes associated with tumour invasivity (ADAM23, uPA, CXCL12, TWIST1, SNAI1 and SNAI2). We found significantly different methylation levels for CXCL12, TWIST1 and SNAI2 genes in tumours of BC patients compared with healthy breast tissues, with mean values of 13.4%, 21.9% and 6.37%, respectively; however, these relations could be influenced by the fact that BC patients were older than healthy donors. With the process of aging, especially after the menopause, mammary glands undergo involution. Both glandular parenchyma and connective tissue atrophy with reduction of adipose tissue. These morphological changes are accompanied with molecular alterations including DNA methylation. In normal breast tissues, the age-related methylation in 204 CpG loci was observed, which progressively alterated in breast tumours . Therefore, for evaluation of cancer-specific methylation profiles, the younger persons cannot be accepted as accurate controls to BC patients, who undergone the aging processes.
In 206 patients, the methylation level of the CXCL12 gene was 12.04% ; however, in another laboratory, quantitative methylation analyses of CXCL12 gene were performed in only three BC samples, with variable methylation levels detected, ranging from 6 to 89% . Other authors observed slightly lower mean levels of TWIST1 methylation than we evidenced in our study, namely, 17.5% and 15.3% [35, 36], but Gort and colleagues determined a higher methylation level of TWIST1 in 76 samples of invasive BC, with a mean of 34% . The median of our results for TWIST1 methylation (18.5%) is approximately two-fold higher than the median of 8.85% measured in BC patients by other researchers using the same pyrosequencing method . In BC cell lines, an inverse association between DNA methylation and the transcription levels of the SNAI1 and SNAI2 genes has been documented . These in vitro results together with ours in BC patients indicate that in addition to CXCL12 and TWIST1 methylation, the aberrant methylation profile of the SNAI2 gene could also contribute to the invasivity of BC.
For the appropriate staging of BC, the tumour size and lymph node status are considered to be the most important clinical categories. However, the histologic grade, which characterizes the degree of cancer cell differentiation, has not yet been integrated into the TNM classification. In a recent meta-analysis, the prognostic importance of histologic grade for overall survival was documented, without any association with tumour size or nodal status . Moreover, integrative bioinformatics analyses have identified the 22-gene tumour aggressiveness grading classifier that indicate the variable features of cancer cells in tumour tissues with different levels of de-differentiation because of individual sets of aberrant genetic changes . To the best of our knowledge, no study has been published that investigated the association between changes in whole-genome methylation profiles and cancer de-differentiation processes.
It is generally accepted that the SLUG (SNAI2) transcription factor, encoded by the SNAI2 gene, induces EMT, initiating cancer dissemination processes. This fact was supported by a previous study, in which the higher expression levels of the SLUG protein and SLUG mRNA were observed in the tumours of patients with metastatic BC or disease recurrence . Moreover, the results of an in vitro study documented epigenetic regulation through promoter methylation in both the SNAI1 and SNAI2 genes during EMT and reverse mesenchymal-epithelial transition . However, our results indicate the other possible role of the SNAI2 gene in BC, which is the de-differentiation of cancer cells. SNAI2 involvement in cell differentiation has been shown in human epidermal progenitor cells. Overexpressed SNAI2 inhibits the differentiation of these cells by binding to the differentiation and adhesion genes across the genome. On the other hand, low levels of SNAI2, and therefore weak binding during epidermal differentiation, lead to the full expression of the differentiation programme. Similar processes could be expected in cancer cells in which high levels of SNAI2 result in more intensive binding of SNAI2 to the genomic targets, determining the differentiation status of epithelial cells. Cancer cells then trigger EMT and the inhibition of cell differentiation . Based on the assumption that decreased methylation could allow the up-regulation of SNAI2 gene expression, reduced differentiation in tumours could be expected, as was observed in our study by the inverse relationship between SNAI2 methylation and histologic grade.
Two important roles in cancer cell spread are played by ADAM proteins performing proteinase activities through the metalloproteinase domain and regulating cell adhesion by their interaction with integrins, which participate in tumour growth and metastasis as well as in tumour angiogenesis . Among the ADAM proteins, proteolytically inactive ADAM23 negatively regulates cancer cell migration by binding αvβ3 integrin to its disintegrin domain . In addition to cell migration, the active integrins promote the interaction between CTCs and thrombocytes, thereby mediating the protection of cancer cells against degradation in the vasculature . In BC cell lines, Costa and colleagues found that ADAM23 promoter hypermethylation was strongly associated with the reduction in mRNA and protein expression, and they also observed higher methylation levels in PTs with more advanced grade . In our study, we investigated whether there are differences in methylation levels in BC patients with different routes of cancer dissemination. We found significantly lower ADAM23 methylation levels in tumours of CTC-positive patients, regardless of their epithelial or mesenchymal phenotype, than in tumours of patients with LNM. Our results indicate that the higher expression level of the ADAM23 gene could contribute to its newly identified function in the haematogenous dissemination of BC.
The dissemination of cancer cells is performed via blood and lymphatic vessels after neo-vascularization and neo-lymphangiogenesis in tumour tissue . In a very interesting study, the molecular features of endothelial cells (EC) located in tumour-associated vasculature were studied to evaluate a new strategy to inhibit tumour growth and cancer dissemination. In this model of EC purified and cultured from ovarian cancer and normal adrenal glands as control tissues, the authors identified 158 highly expressed transcripts. Of these mRNAs, the ADAM23, FAP, GPMNB and PRSS3 were found in the tumour-derived endothelium, but no expression was observed in tumour cells. Moreover, ADAM23, GPMNB and PRSS3 expression only occurred in the blood vessels of human cancer samples . These findings seem to not be organ-specific but rather to be exclusive to tumour-derived EC because expression of several genes that are typical for tumour vasculature was also observed in other types of malignancies in the colon, brain and breast [67,68,69]. In the context of our results, we consider that decreased methylation of the ADAM23 gene could partially represent the molecular profile of EC from the tissue surrounding the tumour body rather than from the tumour cells. Therefore, patients with higher ADAM23 expression levels, which may be the result of decreased promoter methylation, could have contributed to the dense development of tumour blood vessels, resulting in the haematogenous rather than the lymphogenous spread of cancer cells.
In our study, among the selected genes related to the partial processes of cancer cell invasivity, we identified higher cancer-specific methylation levels in the CXCL12, TWIST1 and SNAI2 genes. In contrast to CXCL12 and TWIST1 methylation, the epigenetic regulation of the SNAI2 gene was previously determined only in BC cell line models, not in tumour samples. Furthermore, the decreased methylation of the SNAI2 gene in tumours with histologic grade 3 indicates a new function of a typical EMT gene in cancer cell de-differentiation. The differences in ADAM23 methylation profiles in BC patients with various route of cancer cells dissemination suggest that ADAM23 could participate in haematogenous spread. Our research uncovers new relationships between aberrant methylation profiles and clinical characteristics of advanced BC. However, experimental investigation of the associated molecular mechanisms is needed before the SNAI2 and ADAM23 methylation profiles can contribute to the management of metastatic progression in BC.
A disintegrin and metalloproteases
Analysis of variance
Circulating tumour cells
- epi CTC:
Epithelial circulating tumour cell
Formalin fixed paraffin-embedded tissue
Hepatocyte growth factor
Lymph node metastases
- mes CTC:
Mesenchymal circulating tumour cell
Peripheral blood cells
Quantitative real-time polymerase chain reaction
Vascular endothelial growth factor
Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics. 2012. CA Cancer J Clin. 2015;65:87–108.
Tevaarwerk AJ, Gray RJ, Schneider BP, Smith ML, Wagner LI, Fetting JH, et al. Survival in patients with metastatic recurrent breast cancer after adjuvant chemotherapy: little evidence of improvement over the past 30 years. Cancer. 2013;119:1140–8.
Jones PA, Baylin SB. The epigenomics of cancer. Cell. 2007;128:683–92.
Esteller M. Cancer epigenomics: DNA methylomes and histone-modification maps. Nat Rev Genet. 2007;8:286–98.
Lin IH, Chen DT, Chang YF, Lee YL, Su CH, Cheng C, et al. Hierarchical clustering of breast cancer methylomes revealed differentially methylated and expressed breast cancer genes. PLoS One. 2015;10:e0118453.
Gupta GP, Massagué J. Cancer metastasis: building a framework. Cell. 2006;127:679–95.
Seyfried TN, Huysentruyt LC. On the origin of cancer metastasis. Crit Rev Oncog. 2013;18:43–73.
Páez D, Labonte MJ, Bohanes P, Zhang W, Benhanim L, Ning Y, et al. Cancer dormancy: a model of early dissemination and late cancer recurrence. Clin Cancer Res. 2012;18:645–53.
Chaffer CL, Weinberg RA. A perspective on cancer cell metastasis. Science. 2011;331:1559–64.
Tanner K. Regulation of the basement membrane by epithelia generated forces. Phys Biol. 2012;9:065003.
Duffy MJ, Mullooly M, O'Donovan N, Sukor S, Crown J, Pierce A, et al. The ADAMs family of proteases: new biomarkers and therapeutic targets for cancer? Clin Proteomics. 2011;8:9.
Giebeler N, Zigrino P. A Disintegrin and Metalloprotease (ADAM): Historical Overview of Their Functions. Toxins (Basel). 2016;8:122.
Cal S, Freije JM, López JM, Takada Y, López-Otín C. ADAM 23/MDC3, a human disintegrin that promotes cell adhesion via interaction with the alphavbeta3 integrin through an RGD-independent mechanism. Mol Biol Cell. 2000;11:1457–69.
Wang L, Hoggard JA, Korleski ED, Long GV, Ree BC, Hensley K, et al. Multiple non-catalytic ADAMs are novel integrin α4 ligands. Mol Cell Biochem. 2018;442:29-38.
Verbisck NV, Costa ET, Costa FF, Cavalher FP, Costa MD, Muras A, et al. ADAM23 negatively modulates alpha(v)beta(3) integrin activation during metastasis. Cancer Res. 2009;69:5546–52.
Costa FF, Verbisck NV, Salim AC, Ierardi DF, Pires LC, Sasahara RM, et al. Epigenetic silencing of the adhesion molecule ADAM23 is highly frequent in breast tumors. Oncogene. 2004;23:1481–8.
Fridrichova I, Smolkova B, Kajabova V, Zmetakova I, Krivulcik T, Mego M, et al. CXCL12 and ADAM23 hypermethylation are associated with advanced breast cancers. Transl Res. 2015;165:717–30.
Santibanez JF, Krstic J. Transforming Growth Factor-Beta and Urokinase Type Plasminogen Interplay in Cancer. Curr Protein Pept Sci. 2017; doi.org/10.2174/1389203718666171030103801
Frandsen TL, Holst-Hansen C, Nielsen BS, Christensen IJ, Nyengaard JR, Carmeliet P, et al. Direct evidence of the importance of stromal urokinase plasminogen activator (uPA) in the growth of an experimental human breast cancer using a combined uPA gene-disrupted and immunodeficient xenograft model. Cancer Res. 2001;61:532–7.
Bekes EM, Deryugina EI, Kupriyanova TA, Zajac E, Botkjaer KA, Andreasen PA, et al. Activation of pro-uPA is critical for initial escape from the primary tumor and hematogenous dissemination of human carcinoma cells. Neoplasia. 2011;13:806–21.
Witzel I, Milde-Langosch K, Schmidt M, Karn T, Becker S, Wirtz R, et al. Role of urokinase plasminogen activator and plasminogen activator inhibitor mRNA expression as prognostic factors in molecular subtypes of breast cancer. Onco Targets Ther. 2014;7:2205–13.
Guo Y, Pakneshan P, Gladu J, Slack A, Szyf M, Rabbani SA. Regulation of DNA methylation in human breast cancer. Effect on the urokinase-type plasminogen activator gene production and tumor invasion. J Biol Chem. 2002;277:41571–9.
Pakneshan P, Szyf M, Farias-Eisner R, Rabbani SA. Reversal of the hypomethylation status of urokinase (uPA) promoter blocks breast cancer growth and metastasis. J Biol Chem. 2004a;279:31735–44.
Pakneshan P, Têtu B, Rabbani SA. Demethylation of urokinase promoter as a prognostic marker in patients with breast carcinoma. Clin Cancer Res. 2004b;10:3035–41.
Chen CM, Chen HL, Hsiau TH, Hsiau AH, Shi H, Brock GJ, et al. Methylation target array for rapid analysis of CpG island hypermethylation in multiple tissue genomes. Am J Pathol. 2003;163:37–45.
Chang C, Werb Z. The many faces of metalloproteases: cell growth, invasion, angiogenesis and metastasis. Trends Cell Biol. 2001;11:S37–43.
Sun Y, Mao X, Fan C, Liu C, Guo A, Guan S, et al. CXCL12-CXCR4 axis promotes the natural selection of breast cancer cell metastasis. Tumour Biol. 2014;35:7765–73.
Ray P, Stacer AC, Fenner J, Cavnar SP, Meguiar K, Brown M, et al. CXCL12-γ in primary tumors drives breast cancer metastasis. Oncogene. 2015;34:2043–51.
Zhou W, Jiang Z, Liu N, Xu F, Wen P, Liu Y, et al. Down-regulation of CXCL12 mRNA expression by promoter hypermethylation and its association with metastatic progression in human breast carcinomas. J Cancer Res Clin Oncol. 2009;135:91–102.
Ramos EA, Camargo AA, Braun K, Slowik R, Cavalli IJ, Ribeiro EM, et al. Simultaneous CXCL12 and ESR1 CpG island hypermethylation correlates with poor prognosis in sporadic breast cancer. BMC Cancer. 2010;10:23.
Ramos EA, Grochoski M, Braun-Prado K, Seniski GG, Cavalli IJ, Ribeiro EM, et al. Epigenetic changes of CXCR4 and its ligand CXCL12 as prognostic factors for sporadic breast cancer. PLoS One. 2011;6:e29461.
Kalluri R, Weinberg RA. The basics of epithelial-mesenchymal transition. J Clin Invest. 2009;119:1420–8.
Sánchez-Tilló E, Liu Y, de Barrios O, Siles L, Fanlo L, Cuatrecasas M, et al. EMT-activating transcription factors in cancer: beyond EMT and tumor invasiveness. Cell Mol Life Sci. 2012;69:3429–56.
Martin TA, Goyal A, Watkins G, Jiang WG. Expression of the transcription factors snail, slug, and twist and their clinical significance in human breast cancer. Ann Surg Oncol. 2005;12:488–96.
Cho YH, Yazici H, Wu HC, Terry MB, Gonzalez K, Qu M, et al. Aberrant promoter hypermethylation and genomic hypomethylation in tumor, adjacent normal tissues and blood from breast cancer patients. Anticancer Res. 2010;30:2489–96.
Cho YH, Shen J, Gammon MD, Zhang YJ, Wang Q, Gonzalez K, et al. Prognostic significance of gene-specific promoter hypermethylation in breast cancer patients. Breast Cancer Res Treat. 2012;131:197–205.
Gort EH, Suijkerbuijk KP, Roothaan SM, Raman V, Vooijs M, van der Wall E, et al. Methylation of the TWIST1 promoter, TWIST1 mRNA levels, and immunohistochemical expression of TWIST1 in breast cancer. Cancer Epidemiol Biomark Prev. 2008;17:3325–30.
Branham MT, Marzese DM, Laurito SR, Gago FE, Orozco JI, Tello OM, et al. Methylation profile of triple-negative breast carcinomas. Oncogenesis. 2012;1:e17.
Fackler MJ, McVeigh M, Evron E, Garrett E, Mehrotra J, Polyak K, et al. DNA methylation of RASSF1A, HIN-1, RAR-beta, cyclin D2 and twist in in situ and invasive lobular breast carcinoma. Int J Cancer. 2003;107:970–5.
Chen Y, Wang K, Qian CN, Leach R. DNA methylation is associated with transcription of snail and slug genes. Biochem Biophys Res Commun. 2013;430:1083–90.
Rodenhiser DI, Andrews J, Kennette W, Sadikovic B, Mendlowitz A, Tuck AB, et al. Epigenetic mapping and functional analysis in a breast cancer metastasis model using whole-genome promoter tiling microarrays. Breast Cancer Res. 2008;10:R62.
Paduch R. The role of lymphangiogenesis and angiogenesis in tumor metastasis. Cell Oncol (Dordr). 2016;39:397–410.
Stacker SA, Baldwin ME, Achen MG. The role of tumor lymphangiogenesis in metastatic spread. FASEB J. 2002;16:922–34.
Lee E, Pandey NB, Popel AS. Crosstalk between cancer cells and blood endothelial and lymphatic endothelial cells in tumour and organ microenvironment. Expert Rev Mol Med. 2015;17:e3.
Sigurdsson V, Hilmarsdottir B, Sigmundsdottir H, Fridriksdottir AJ, Ringnér M, Villadsen R, et al. Endothelial induced EMT in breast epithelial cells with stem cell properties. PLoS One. 2011;6:e23833.
Scher HI, Heller G, Molina A, Attard G, Danila DC, Jia X, et al. Circulating tumor cell biomarker panel as an individual-level surrogate for survival in metastatic castration-resistant prostate cancer. J Clin Oncol. 2015;33:1348–55.
Romiti A, Raffa S, Di Rocco R, Roberto M, Milano A, Zullo A, et al. Circulating tumor cells count predicts survival in colorectal cancer patients. J Gastrointestin Liver Dis. 2014;23:279–84.
Hou JM, Krebs MG, Lancashire L, Sloane R, Backen A, Swain RK, et al. Clinical significance and molecular characteristics of circulating tumor cells and circulating tumor microemboli in patients with small-cell lung cancer. J Clin Oncol. 2012;30:525–32.
Zhang L, Riethdorf S, Wu G, Wang T, Yang K, Peng G, et al. Meta-analysis of the prognostic value of circulating tumor cells in breast cancer. Clin Cancer Res. 2012;18:5701–10.
Skobe M, Hawighorst T, Jackson DG, Prevo R, Janes L, Velasco P, et al. Induction of tumor lymphangiogenesis by VEGF-C promotes breast cancer metastasis. Nat Med. 2001;7:192–8.
Wong SY, Hynes RO. Lymphatic or hematogenous dissemination: how does a metastatic tumor cell decide? Cell Cycle. 2006;5:812–7.
Nathanson SD, Kwon D, Kapke A, Alford SH, Chitale D. The role of lymph node metastasis in the systemic dissemination of breast cancer. Ann Surg Oncol. 2009;16:3396–405.
Wolff AC, Hammond ME, Hicks DG, Dowsett M, McShane LM, Allison KH, et al. Recommendations for human epidermal growth factor receptor 2 testing in breast cancer: American Society of Clinical Oncology/College of American Pathologists clinical practice guideline update. Arch Pathol Lab Med. 2014;138:241–56.
Zmetakova I, Danihel L, Smolkova B, Mego M, Kajabova V, Krivulcik T, et al. Evaluation of protein expression and DNA methylation profiles detected by pyrosequencing in invasive breast cancer. Neoplasma. 2013;60:635–46.
Cierna Z, Mego M, Janega P, Karaba M, Minarik G, Benca J, et al. Matrix metalloproteinase 1 and circulating tumor cells in early breast cancer. BMC Cancer. 2014;14:472.
Mego M, Karaba M, Minarik G, Benca J, Sedlackova T, Tothova L, et al. Relationship between circulating tumor cells, blood coagulation, and urokinase-plasminogen- activator system in early breast cancer patients. Breast J. 2015;21:155–60.
Mego M, Cholujova D, Minarik G, Sedlackova T, Gronesova P, Karaba M, et al. CXCR4-SDF-1 interaction potentially mediates trafficking of circulating tumor cells in primary breast cancer. BMC Cancer. 2016;16:127.
Johnson KC, Koestler DC, Cheng C, Christensen BC. Age-related DNA methylation in normal breast tissue and its relationship with invasive breast tumor methylation. Epigenetics. 2014;9:268–75.
Cheol Kim D, Thorat MA, Lee MR, Cho SH, Vasiljević N, Scibior-Bentkowska D, et al. Quantitative DNA methylation and recurrence of breast cancer: a study of 30 candidate genes. Cancer Biomark. 2012;11:75–88.
Schwartz AM, Henson DE, Chen D, Rajamarthandan S. Histologic grade remains a prognostic factor for breast cancer regardless of the number of positive lymph nodes and tumor size: a study of 161 708 cases of breast cancer from the SEER program. Arch Pathol Lab Med. 2014;138:1048–52.
Aswad L, Yenamandra SP, Ow GS, Grinchuk O, Ivshina AV, Kuznetsov VA. Genome and transcriptome delineation of two major oncogenic pathways governing invasive ductal breast cancer development. Oncotarget. 2015;6:36652–74.
Mistry DS, Chen Y, Wang Y, Zhang K, Sen GL. SNAI2 controls the undifferentiated state of human epidermal progenitor cells. Stem Cells. 2014;32:3209–18.
Reiss K, Ludwig A, Saftig P. Breaking up the tie: disintegrin-like metalloproteinases as regulators of cell migration in inflammation and invasion. Pharmacol Ther. 2006;111:985–1006.
Felding-Habermann B, O'Toole TE, Smith JW, Fransvea E, Ruggeri ZM, Ginsberg MH, et al. Integrin activation controls metastasis in human breast cancer. Proc Natl Acad Sci U S A. 2001;98:1853–8.
Witte MH, Dellinger MT, McDonald DM, Nathanson SD, Boccardo FM, Campisi CC, et al. Lymphangiogenesis and hemangiogenesis: potential targets for therapy. J Surg Oncol. 2011;103:489–500.
Ghilardi C, Chiorino G, Dossi R, Nagy Z, Giavazzi R, Bani M. Identification of novel vascular markers through gene expression profiling of tumor-derived endothelium. BMC Genomics. 2008;9:201.
St Croix B, Rago C, Velculescu V, Traverso G, Romans KE, Montgomery E, et al. Genes expressed in human tumor endothelium. Science. 2000;289:1197–202.
Madden SL, Cook BP, Nacht M, Weber WD, Callahan MR, Jiang Y, et al. Vascular gene expression in nonneoplastic and malignant brain. Am J Pathol. 2004;165:601–8.
Parker BS, Argani P, Cook BP, Liangfeng H, Chartrand SD, Zhang M, et al. Alterations in vascular gene expression in invasive breast carcinoma. Cancer Res. 2004;64:7857–66.
The authors would like to thank the surgeons Marian Karaba, Juraj Benca and Martin Bohac for providing the collection of breast tumour and normal breast tissue samples.
This publication is the result of the implementation of several research projects funded by the Scientific Grant Agency of the Ministry of Education, science, research and sport of the Slovak Republic and Slovak Academy of Sciences (grants VEGA 2/0092/15 and VEGA 2/0102/17), the European Fund for Regional Development through Operation Program of Research and Development (project code ITMS 26240220058), and the Slovak Research and Development Agency (grants APVV-0076-10, APVV-14-0327 and APVV-16-0010). These agencies had no role in the design of the study and collection, analysis, and interpretation of data and in writing and in the decision to submit the manuscript for publication.
Availability of data and materials
All data generated or analysed during this study are included in this published article.
Ethics approval and consent to participate
This study was approved by the Institutional Review Board of the National Cancer Institute of Slovakia. All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. Informed consent was obtained from all individual participants included in the study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
DNA methylation levels in various types of samples from BC patients with different histological grade. A table containing methylation results for all six genes and statistical evaluation between two groups of patients with G1 + G2 and G3 tumours in samples of peripheral blood cells, primary tumours and lymph node metastases. (DOCX 15 kb)
About this article
Cite this article
Kalinkova, L., Zmetakova, I., Smolkova, B. et al. Decreased methylation in the SNAI2 and ADAM23 genes associated with de-differentiation and haematogenous dissemination in breast cancers. BMC Cancer 18, 875 (2018). https://doi.org/10.1186/s12885-018-4783-x
- SNAI2 methylation
- ADAM23 methylation
- Breast cancer
- Cancer cell de-differentiation
- Haematogenous dissemination