Aberrant methylation profiles and silencing of a small subset of tumor suppressors and cancer related genes involved in both the primary tumor and lymph node metastasis has been investigated for breast cancer [
11,
22-
24]. Several studies on breast cancer revealed some similarities and differences between promoter methylation pattern of the studied genes in primary and lymph node metastatic [
16,
19,
25]. These investigations highlight the important role of aberrant promoters’ methylation in the metastatic process. In the present study, we implemented the methylation signature of the 12 breast cancer candidate genes (
APC, BIN1, BMP6, BRCA1, CST6, ESR-b, GSTP1, P14, P16, P21, PTEN and
TIMP3) by comparing lymph nodes metastasis to their matched primary tumor tissues and normal tissues from the same breast cancer patients.
The quantitative methylation analysis of the 12 studied genes in the present cohort showed higher methylation proportion for the primary tumor tissue versus matched normal tissue and the differences were significant for
APC, BIN1, BMP6, BRCA1, CST6, ESR-b, P16, PTEN and
TIMP3 promoter regions (
P
<

0.05,
P
<

0.05,
P
<

0.01,
P
<

0.0001,
P
<

0.01,
P
<

0.01,
P
<

0.05,
P
<

0.05,
P
<

0.01; respectively). Among the significant methylated genes,
APC, BMP6, BRCA1 and
P16 represented higher methylation proportions in matched lymph node metastasis compared to those of the normal tissue (
P
<

0.05 and
P
<

0.01,
P
<

0.0001,
P
<

0.05; respectively). Present findings provided evidence of differences in methylation status between primary tumors and their corresponding matched lymph nodes and demonstrated methylation heterogeneity between primary tumors and metastatic lesion which are in line with previous reports described about primary tumor and metastasis in breast, gastric and colorectal cancers [
13,
26,
27]. The results also indicated that some of the cancer specific changes become altered over the metastasis procedure. The mechanism for losing or gaining methylation in lymph nodes metastatic is still not clear while the alteration of the methylation signature from primary to metastasis might be due to the adaptation response of the disseminated cells to the microenvironment at the site of colonization [
28].
The molecular function, biological processes and contribution of the studied genes to developing metastasis were analyzed by ResNet® 7 (Mammal). The pathway analysis revealed that
BMP6BRCA1 and
P16 have a role in prevention of neoplasm metastasis (Figure

). The relation of DNA methylation for
BRCA1 and
P16 with tumor recurrence has been reported with a high value in breast cancer patients [
29]. Moreover, association of
P16 hypermethylation with cancer progression and lymph node invasion has been shown in different studies [
30,
31]. Then, the aberrant methylation signatures of these genes found in the metastatic lymph node can provide a further clue to establish useful biomarkers for screening metastasis in breast cancer. Such biomarkers will need to be additionally explored for the potential opposing functional effects of specific hypermethylated CpG sites on transcriptional activity as well as for absolute percent methylation cutoffs for each breast tumor type that would enable reliable utility of them in the clinical laboratory setting. The possibility of using these hypermethylated genes as biomarkers in our previous study was investigated with the aim of developing a blood based panel for plasma and serum samples of breast cancer patients [
18]. Thus the detection of these three hypermethylated genes either as tissue specific biomarkers or as circulating biomarkers may give insight into the prognosis and therapeutic management of the breast cancer patients. Moreover, we identified significantly greater hypermethylation proportion of
BMP6 in the lymph node metastasis than their primary tumor tissue (
P
<

0.05).
BMP6 is a member of TGF-β super family and critically involved in many developmental processes [
32]. Recently, the close association of
BMP6 with progression of tumorigenesis and regulation of invasion for tumor cells has been reported [
33]. It has been assumed that demethylation of
BMP6 and re-expression of this gene might modulate metastasis and invasion in breast cancer [
34,
35].