Low methylation levels of the SFRP1 gene are associated with the basal-like subtype of breast cancer

  • Authors:
    • Young Ju Jeong
    • Hye Yeon Jeong
    • Jin Gu Bong
    • Sung Hwan Park
    • Hoon Kyu Oh
  • View Affiliations

  • Published online on: March 6, 2013     https://doi.org/10.3892/or.2013.2335
  • Pages: 1946-1954
Metrics: Total Views: 0 (Spandidos Publications: | PMC Statistics: )
Total PDF Downloads: 0 (Spandidos Publications: | PMC Statistics: )


Abstract

Epigenetic analyses have shown that aberrant DNA methylation signatures are associated with breast cancer molecular subtypes. In this study, we analyzed the methylation status of breast cancer-related genes in relation to the molecular subtypes and investigated whether the basal-like subtype displays distinct methylation profiles. By using pyrosequencing, we analyzed the DNA methylation status of 5 candidate genes in 60 breast cancer samples. We compared the methylation frequency across the molecular subtypes and analyzed the correlation between methylation levels and clinicopathological characteristics. A total of 59 cases displayed aberrant methylation. Amplification during polymerase chain reaction analysis failed in 1 case. The median methylation levels of the secreted frizzled-related protein 1 (SFRP1) gene were significantly lower in the basal-like subtype compared to the luminal A, luminal B and human epidermal growth factor receptor 2 (HER2) subtypes. Cadherin 13 (H-cadherin; CDH13) methylation levels were significantly higher in the HER2 tumors compared to the luminal A and basal-like subtypes. A comparison of the methylation status with clinicopathological characteristics revealed that the expression of bcl-2, progesterone receptor and epidermal growth factor receptor were associated with SFRP1 gene methylation status. Our results indicate that the basal-like subtype is associated with low methylation levels of the SFRP1 gene, suggesting that the methylation levels of specific breast cancer genes may potentially serve as epigenetic biomarkers and prognostic factors.

Introduction

Breast cancer is a molecularly, biologically and clinically heterogeneous disease. Previous microarray profiling studies on breast cancer have identified subtypes that are associated with different clinical outcomes (1,2). These subtypes are classified as luminal A, luminal B, human epidermal growth factor receptor 2 (HER2), basal-like and normal-like subtypes. The identification of gene expression-based breast cancer subtypes is considered a critical means of prognostication (3) and the basal-like subtype is an aggressive tumor that has a poor prognosis and no specific targeted therapy. Since it is not always feasible to obtain gene expression array information, a simple classification using a combination of immunohistochemical markers has been proposed and adopted in clinical practice (4,5). However, there is still a need for well-defined biomarker panels that allow breast cancer subtyping for clinical diagnostics and the management of the disease.

Breast carcinogenesis is a multistep process resulting from the accumulation of genetic alterations, as well as epigenetic changes, such as promoter methylation and histone modification (6,7). Promoter hypermethylation at specific gene loci leading to gene silencing is a major mechanism of epigenetic inactivation in cancer cells. A previous genome-wide study on breast cancer has led to the identification of a number of tumor suppressor genes that are inactivated by promoter hypermethylation (8), and epigenetic analyses have shown aberrant DNA methylation signatures associated with the molecular subtypes of breast cancer (9,10). However, limited information is available on the methylation status of candidate genes associated with each molecular subtype.

The purpose of this study was to analyze the methylation status of breast cancer-related genes according to the molecular subtypes found in Korean women, and to investigate whether the basal-like subtype displays distinct methylation profiles compared with the other subtypes. We included 5 genes that are involved in breast carcinogenesis and are commonly methylated in breast cancer [cadherin 13 (H-cadherin; CDH13), secreted frizzled-related protein 1 (SFRP1), fragile histidine triad (FHIT), Syk and retinoblastoma protein-interacting zinc-finger gene 1 (RIZ1)] and analyzed the methylation status of these genes using a sensitive and quantitative pyrosequencing assay.

Materials and methods

Patients and tumor characteristics

A total of 60 sporadic invasive ductal carcinoma (IDC) tissue samples were obtained from the Daegu Catholic University Hospital (Daegu, Korea). Each sample represented the 4 major molecular subtypes, encompassing the basal-like, HER2, luminal A and luminal B subtypes. All specimens were reviewed by an experienced pathologist. We subclassified the breast cancer samples according to immunohistochemical findings for the estrogen receptor (ER), progesterone receptor (PR), HER2 oncogene and Ki-67 labeling index (11). The basal-like subtype was defined as HER2-negative, ER- and PR-negative (triple-negative) breast cancer. The HER2 subtype was defined as HER2-positive, ER- and PR-negative. The luminal B subtype was defined as HER2-positive and ER- and/or PR-positive breast cancer (HER2-positive). The luminal A subtype was defined as ER- and/or PR-positive, HER2-negative breast cancer with a low Ki-67 index. ER- and/or PR-positive, HER2-negative breast cancer with a high Ki-67 status was classified as the luminal B (HER2-negative) subtype. Ethics approval for the study was obtained from the institutional review board at Daegu Catholic University Hospital.

Construction of tissue microarrays (TMAs)

Representative paraffin tumor blocks were selected according to the primary evaluation of hematoxylin and eosin (H&E)-stained slides before they were prepared for TMA analysis. Two tumor tissue cores (1 mm in diameter) were taken from each of the donor breast cancer tissue blocks using a manual punch arrayer (Quick-Ray™; Uni-Tech Science, Seoul, Korea). The cores were placed in a new recipient paraffin block that ultimately contained 72–96 tissue cores. Each array block contained both tumor and control tissue samples. Multiple sections (5-μm-thick) were cut from the TMA blocks and then mounted onto microscope slides. The TMA H&E-stained sections were reviewed under a light microscope to confirm the presence of representative tumor areas.

Immuohistochemical staining and interpretation

Immunohistochemical analysis was performed on 5-μm-thick TMA tissue sections using the Bond Polymer Intense Detection system (Leica Microsystems, Mount Waverley, Victoria, Australia) according to the manufacturer's instructions with minor modifications. Briefly, the 5-μm-thick sections of formalin-fixed and paraffin-embedded TMA tissues were deparaffinized with Bond Dewax Solution (Leica Microsystems), and an antigen retrieval procedure was performed using Bond ER Solution (Leica Microsystems) for 30 min at 100°C. The endogenous peroxidase was quenched by a 5-min incubation with hydrogen peroxide. Sections were incubated for 15 min at an ambient temperature with commercially available primary monoclonal antibodies for ER (1:100, clone 6F11; Novocastra), PR (1:100, clone 16; Novocastra), HER2 (1:250, A0485; Dako), Ki-67 (1:200, MM1-L; Novocastra), Bcl-2 (1:4, clone 124; Dako), p53 (1:200, BP53.12; Zymed) and epidermal growth factor receptor (EGFR) (1:100, clone EGFR.25; Novocastra) using a biotin-free polymeric horseradish peroxidase-linker antibody conjugate system in a Bond-Max automatic slide stainer (Leica Microsystems).

A cut-off value of 10% for the stained nuclei was used to define ER and PR positivity. Cytoplasmic staining of any intensity in >10% of the tumor cells was scored as positive for Bcl-2. Membranous staining for HER-2 with strong complete staining in 10% of the tumor cells was regarded as HER-2 overexpression. p53 staining was scored positive if >10% of the cells were stained with a strong intensity. The Ki-67 labeling index was expressed as a percentage and was graded as ‘high’ if the number of positive cells was ≥14%.

DNA extraction and sodium bisulfate treatment

For DNA extraction, 8 tissue sections (5–10-μm-thick) were obtained from the paraffin-embedded primary breast cancer tissues. Genomic DNA was isolated using the QIAamp DNA FFPE Tissue kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. The purified DNA was quantified using a ND-1000 spectrophotometer (NanoDrop Technologies, Inc., Wilmington, DE, USA). The quality of the DNA was verified by gel electrophoresis. Sodium bisulfate modification of 200–500 ng genomic DNA was performed using the EZ DNA Methylation-Gold kit (Zymo Research, Orange, CA, USA) according to the manufacturer's instructions.

Candidate selection

Over 100 individual candidate genes have been reported to be commonly hypo- and hypermethylated in breast cancer. We carried out literature searches on PubMed http://www.ncbi.nlm.nih.gov/pubmed for the keywords (breast cancer, cancer and methylation) and searched serial analysis of gene expression (SAGE) data (GeneCards, http://www.genecards.org/cgi-bin/cardsearch). We selected 5 genes, cadherin 13 (H-cadherin; CDH13), secreted frizzled-related protein 1 (SFRP1), FHIT, Syk and RIZ1, and functional annotation of the candidate genes was carried out using the functional annotation table function in the DAVID database http://david.abcc.ncifcrf.gov (12).

Pyrosequencing

Methylation analysis was carried out using pyrosequencing. Primers were designed using the PyroMark Assay Design program version 2.0.1.15 (Qiagen) and the sequences are presented in Table I. Polymerase chain reaction (PCR) was carried out using bisulfate-treated DNA under the following conditions: 95°C for 5 min; 45 cycles of 95°C for 30 sec, 55°C for 30 sec and 72°C for 30 sec; and a final extension of 5 min at 72°C. PCR was conducted using a PCR premix (Enzynomics, Daejeon, Korea), and the quality and quantity of the PCR product was confirmed by performing agarose gel (2%) electrophoresis with loading 4 μl of 20 PCR products. Pyrosequencing was performed using the Pyro Gold kit and PSQ 96MA instrument (Qiagen) as instructed by the manufacturer. The methylation index (MtI) of each gene in each sample was calculated as the average value of mC/(mC + C) for all examined CpG sites in target regions. All experiments included a negative control without a template.

Table I

Primer sequences used for PCR and pyrosequencing.

Table I

Primer sequences used for PCR and pyrosequencing.

Primer namePrimer sequence (5′→3′)
CDH13_F TAAGGAAAATATGTTTAGTGTAGT
CDH13_R AAATTCTCCACTACATTTTATCC
CDH13_S GTGTAGTAGAGTGTATGAATGAAAA
SFRP1_F TTTTAGGAGGTTTTTGGAAGT
SFRP1_R ACTCTACCCCCTATTCTCC
SFRP1_S AGGTTTTTGGAAGTTTG
FHIT_F GGGAGGTAAGTTTAAGTGGAATATTG
FHIT_R CCACTAAACTCCCAAATAATAACCTAAC
FHIT_S GTAAGTTTAAGTGGAATATTGT
Syk_F TTAGTAGGGAGGGTTAGGG
Syk_R CTCATTTTAAACAACTTCCTTAAC
Syk_S ATATTGGGAGGAAGTG
RIZ1_F AGTAAGTTTTTTAAGGGTAGGATTAT
RIZ1_R CCCTAATACCCAAAAACAATAACCAA
RIZ1_S GTTTTTTAAGGGTAGGATTATTAT

[i] CDH13, cadherin 13 (H-cadherin); SFRP1, secreted frizzled-related protein 1; RIZ1, retinoblastoma protein-interacting zinc-finger gene 1. F, forward primer for PCR, R, reverse primer for PCR; S, primer used for pyrosequencing.

Statistical analysis

Statistical analyses were carried out using SPSS version 15.0 software (SPSS Inc., Chicago, IL, USA). A one-sample Kolmogorov-Smirnov test was used to evaluate the normal distribution fit of continuous parameters. The clinicopathological characteristics were compared across the 4 different breast cancer subtypes using the χ2 test or Fisher's exact test for categorical data, and ANOVA or the non-parametric Kruskal-Wallis test for continuous data. A comparison of the mean methylation frequencies across the subtypes was performed using ANOVA or the Kruskal-Wallis test, and distributions of methylation levels across the different subtypes were depicted for each gene using box plots. Associations between methylation status and clinicopathological characteristics were assessed using the Student's t-test or the non-parametric Mann-Whitney U test for categorical variables, and correlation between 2 continuous variables was assessed using correlation analysis. All tests were two-sided and a P-value <0.05 was considered to indicate a statistically significant difference.

Clustering analyses were performed using GeneSpring GX version 7.3 (Agilent Technologies Inc., Santa Clara, CA, USA) on the basis of the mean methylation levels of genes and the CpG sites. Hierarchical clustering was performed using Pearson's correlation distance and average linkage.

Results

Clinicopathological characteristics

The patient characteristics are presented in Table II. The average age of the 60 patients with invasive breast cancer was 51.77±13.22 years (range, 26–90 years). A total of 15 cases were included for each molecular subtype in the 60 breast cancer samples. TNM staging was as follows: stage I, 29 patients (48.3%); stage II, 21 patients (35.0%); stage III, 6 patients (10.0%); and stage IV, 4 patients (6.7%).

Table II

General patient characteristics.

Table II

General patient characteristics.

CharacteristicsValue
Age (years), mean (range)51.77±13.22 (26–90)
Menopausal status, n (%)
 Pre-menopausal27 (45.8)
 Post-menopausal32 (54.2)
Tumor size (cm), mean (range)1.80±0.93 (0.10–4.50)
Histological grade, n (%)
 I13 (21.7)
 II11 (18.3)
 III36 (60.0)
Nodal involvement, n (%)
 Negative40 (69.0)
 Positive18 (31.0)
Distant metastasis, n (%)
 Negative58 (96.7)
 Positive2 (3.3)
Molecular subtype, n (%)
 Luminal A15 (25.0)
 Luminal B15 (25.0)
 HER215 (25.0)
 Basal-like15 (25.0)
Methylation level of candidate gene, mean %
CDH1315.66±13.84
SFRP115.67±11.21
FHIT3.43±0.97
Syk8.73±5.32
RIZ148.30±11.55
Lymphovascular invasion, n (%)
 Negative39 (66.1)
 Positive20 (33.9)
ER, n (%)
 Negative31 (51.7)
 Positive29 (48.3)
PR, n (%)
 Negative33 (55.0)
 Positive27 (45.0)
HER2 overexpression, n (%)
 Negative30 (50.0)
 Positive30 (50.0)
Ki-67, n (%)
 <14%25 (41.7)
 ≥14%35 (58.3)

[i] HER2, human epidermal growth factor receptor 2; CDH13, cadherin 13 (H-cadherin); SFRP1, secreted frizzled-related protein 1; RIZ1, retinoblastoma protein-interacting zinc-finger gene 1; ER, estrogen receptor PR, progesterone receptor.

Table III presents the clinicopathological characteristics according to the 4 breast cancer subtypes. The basal subtype was characterized by a high histologicical grade (P<0.001), low extensive intraductal component (EIC) (P<0.001), the presence of necrosis (P<0.001), a high Ki-67 level (P<0.001) and a positive expression of EGFR (P<0.001).

Table III

Clinicopathological characteristics according to breast cancer subtype.

Table III

Clinicopathological characteristics according to breast cancer subtype.

Subtype

CharacteristicsLuminal ALuminal BHER2Basal-likeP-value
Age (years), mean ± SD57.3±15.246.9±7.150.3±11.852.6±15.90.273
Menopausal status, n (%)
 Pre-menopausal7 (46.7)8 (57.1)7 (46.7)5 (33.3)0.643
 Post-menopausal8 (53.3)6 (42.9)8 (53.3)10 (66.7)
Tumor size (cm), mean ± SD1.6±0.71.6±1.01.7±0.82.3±1.00.082
Histological grade, n (%)
 I9 (60.0)2 (13.3)2 (15.4)0 (0.0)<0.001
 II5 (33.3)4 (26.7)2 (13.3)0 (0.0)
 III1 (6.7)9 (60.0)11 (73.3)15 (100.0)
Nodal involvement, n (%)
 Negative10 (66.7)10 (66.7)12 (85.7)8 (57.1)0.432
 Positive5 (33.3)5 (33.3)2 (14.3)6 (42.9)
Distant metastasis, n (%)
 Negative15 (100.0)15 (100.0)14 (93.3)14 (93.3)1.000
 Positive0 (0.0)0 (0.0)1 (6.7)1 (6.7)
Lymphovascular invasion, n (%)
 Negative10 (66.7)10 (66.7)11 (73.3)8 (57.1)0.844
 Positive5 (33.3)5 (33.3)4 (26.7)6 (42.9)
EIC (%), mean ± SD13.9±19.731.7±38.812.9±37.73.6±9.5<0.001
Necrosis, n (%)
 Negative14 (100.0)11 (73.3)6 (40.0)1 (12.5)<0.001
 Positive0 (0.0)4 (26.7)9 (60.0)7 (87.5)
Ki-67, n (%)
 <14%15 (100.0)6 (40.0)4 (26.7)0 (0.0)<0.001
 ≥14%0 (0.0)9 (60.0)11 (73.3)15 (100.0)
ER, n (%)
 Negative0 (0.0)1 (6.7)15 (100.0)15 (100.0)<0.001
 Positive15 (100.0)14 (93.3)0 (0.0)0 (0.0)
PR, n (%)
 Negative1 (6.7)2 (13.3)15 (100.0)15 (100.0)<0.001
 Positive14 (93.3)13 (86.7)0 (0.0)0 (0.0)
HER2, n (%)
 Negative15 (100.0)0 (0.0)0 (0.0)15 (100.0)<0.001
 Positive0 (0.0)15 (100.0)15 (100.0)0 (0.0)
Bcl-2, n (%)
 Negative13 (86.7)15 (27.3)13 (86.7)14 (93.3)0.740
 Positive2 (13.3)0 (0.0)2 (13.3)1 (6.7)
p53, n (%)
 Negative2 (13.3)1 (6.7)3 (20.0)5 (33.3)0.353
 Positive13 (86.7)14 (93.3)12 (80.0)10 (66.7)
EGFR, n (%)
 Negative14 (93.3)14 (93.3)7 (50.0)0 (90.0)<0.001
 Positive1 (6.7)1 (6.7)7 (50.0)15 (100.0)

[i] EIC, extensive intraductal component; ER, estrogen receptor; PR, progesterone receptor; HER2, HER2, human epidermal growth factor receptor 2; EGFR, epidermal growth factor receptor.

Methylation levels in different molecular subtypes of breast cancer

A total of 59 cases showed aberrantly methylated genes. Amplification of 1 sample failed during PCR analysis. The mean methylation levels of CDH13, SFRP1, FHIT, Syk and RIZ1 were 15.66±13.84, 15.67±11.21, 3.43±0.97, 8.73±5.32 and 48.30%±11.55%, respectively. The methylation status of CDH13 and SFRP1 was significantly different according to breast cancer molecular subtypes (Table IV, Fig. 1). The CDH13 methylation level was significantly higher in HER2 tumors compared to the luminal and basal-like subtypes (P=0.006 and P=0.012, respectively) (Table V). The median methylation level and average methylation ratio of the SFRP1 gene were significantly lower in the basal-like subtype compared to the luminal A, luminal B and HER2 subtypes (P=0.002, P<0.001 and P=0.003, respectively). FHIT, Syk and RIZ1 methylation levels showed no significant differences among the molecular subtypes.

Table IV

Methylation levels in breast cancer subtypes.

Table IV

Methylation levels in breast cancer subtypes.

Methylation level (mean ± SD)

GeneLuminal ALuminal BHER2Basal-likeP-value
CDH139.2±5.918.2±17.425.1±16.211.1±7.7.80.015
SFRP115.5±7.824.3±10.718.0±10.74.0±2.0<0.001
FHIT3.4±1.13.2±1.13.8±1.03.3±0.60.367
Syk7.7±5.38.2±5.79.0±5.210.4±5.20.626
RIZ153.0±11.548.7±10.946.2±11.444.5±12.10.257

[i] HER2, human epidermal growth factor receptor 2; CDH13, cadherin 13 (H-cadherin); SFRP1, secreted frizzled-related protein 1; RIZ1, retinoblastoma protein-interacting zinc-finger gene 1.

Table V

Comparison of methylation levels between breast cancer subtypes.

Table V

Comparison of methylation levels between breast cancer subtypes.

Basal-like vs. LuminalBasal-like vs. HER2Luminal vs. HER2



GeneMean differenceP-valueMean differenceP-valueMean differenceP-value
CDH13−2.6110.533−14.0220.006−11.4110.012
SFRP1−16.089<0.001−14.0760.0012.0140.927
FHIT−0.0470.879−0.5410.137−0.4940.121
Syk2.3760.1981.3920.510−0.9840.571
RIZ1−6.3130.113−1.7440.7044.5700.236

[i] HER2, human epidermal growth factor receptor 2; CDH13, cadherin 13 (H-cadherin); SFRP1, secreted frizzled-related protein 1; RIZ1, retinoblastoma protein-interacting zinc-finger gene 1.

Gene-specific patterns of methylation for classification of breast cancer

The varying methylation frequencies of the CDH13 and SFRP1 genes in the different breast cancer molecular subtypes provides evidence for subtype-specific methylation profiling in breast cancer classification. In particular, a low frequency of SFRP1 methylation was significantly associated with the basal-like subtype (P<0.001). We determined the area under the curve (AUC) of the receiver operating characteristic (ROC) curve for the mean methylation level of the SFRP1 gene. For the cut-off value of 7.50%, the AUC was 0.941, with a sensitivity of 91.7% and a specificity of 84.2%.

Unsupervised hierarchical clustering based on the 5 genes and the variable methylated CpG loci of each gene resulted in the formation of multiple small clusters that had no similar biological patterns. As shown in Fig. 2, the results showed classes that were not well separated and that did not appear to be associated with a molecular subtype.

Methylation status and clinicopathological parameters

A comparison of the methylation status with the clinicopathologic characteristics revealed that the negative expression of ER, PR and HER-2, and the positive expression of Bcl-2 and EGFR were associated with a low level of SFRP1 gene methylation (P=0.002, P=0.015, P<0.001, P=0.001 and P<0.001, respectively) (Table VI). The correlations between the methylation status of RIZ1, Syk and FHIT and the clinicopathological characteristics were dissimilar, RIZ1 was associated with menopausal status and age (P=0.015 and P=0.042, respectively). Syk was associated with histological grade (P=0.016).

Table VI

Association between clinicopathological characteristics and DNA methylation levels (P-value).

Table VI

Association between clinicopathological characteristics and DNA methylation levels (P-value).

CDH13SFRP1FHITSykRIZ1
Age0.9500.8230.5120.3450.042
Menopausal status0.8590.6330.3160.8890.015
Stage0.5190.2990.2030.8480.797
Tumor size0.3740.2270.7810.8420.256
LN(+)0.5450.1730.4760.2940.855
Metastasis0.3070.1800.1870.2840.237
Histological grade0.1580.1680.7580.0160.807
LVI0.6080.3550.7880.8960.983
ER0.2480.0020.7680.1700.174
PR0.1380.0150.4280.1760.380
HER20.006<0.0010.4980.8190.637
Ki-670.1940.1620.7590.0270.125
Bcl-20.3080.0010.5450.5500.933
p530.2710.3470.8020.7610.945
EGFR0.973<0.0010.8580.9740.209
Necrosis0.5340.0640.7150.6760.751

[i] CDH13, cadherin 13 (H-cadherin); SFRP1, secreted frizzled-related protein 1; RIZ1, retinoblastoma protein-interacting zinc-finger gene 1; LN(+), lymph node status; LVI, lymphovascular invasion; ER, estrogen receptor; PR, progesterone receptor; HER2, HER2, human epidermal growth factor receptor 2; EGFR, epidermal growth factor receptor.

Discussion

Gene expression microarray analysis has made it possible to identify distinct molecular subtypes that are associated with different clinical outcomes (13). Previous studies have documented the aberrant methylation of CpG islands in gene promoters in breast carcinogenesis (7,8), and have indicated that specific DNA methylation patterns may be associated with some of the known breast cancer subtypes (9,10).

We observed a significantly lower methylation level at a CpG site in the SFRP1 gene in the basal-like breast cancer subtype when compared with other molecular subtypes, consistent with the results of a previous study (13). SFRP1 is a member of the SFRP family (14) and a putative inhibitor of Wnt signaling (15,16) that plays a key role in embryonic development, cell differentiation and proliferation, as well as tumor development and progression (17). It has been described that SFRP1 expression is frequently lost or downregulated in breast cancer (18,19), and the loss of SFRP1 expression is associated with poor prognosis, indicating a putative tumor suppressor gene function of SFRP1(19). Furthermore, a recent study indicated that SFRP1 has potential as a novel biomarker for the triple-negative breast cancer phenotype (20). Thus, SFRP1 hypermethylation may contribute to breast carcinogenesis, and may be useful as a novel prognostic marker of breast cancer, particularly basal-like tumors. In our study, 83.3% of the breast cancer cases had an aberrant methylation of the SFRP1 gene. This result suggests that promoter hypermethylation is the predominant mechanism of SFRP1 gene silencing in human breast cancer, as shown in a previous study (21), although SFRP1 protein expression analysis was not included in this study. Of the 60 cases that we analyzed, 15 were of the basal-like subtype and tended to have decreased methylation of the SFRP1 gene, suggesting that the DNA methylation of this gene may contribute to the phenotype of breast cancer subtypes, particularly the basal-like subtype.

The basal-like subtype of breast cancer has a triple-negative phenotype with a poor prognosis and no specific targeted therapy, despite an increased response to chemotherapy compared to other breast cancer subtypes. Of note, a recent study reported that, in triple-negative breast cancer, SFRP1 expression significantly correlates with an increased sensitivity to neoadjuvant chemotherapy with taxane/anthracycline, and preliminary experiments with siRNA-mediated knockdown of SFRP1 showed a significantly decreased sensitivity to paclitaxel, doxorubicin and cisplatin (20). In our study, the low frequency of SFRP1 methylation in the basal-like subtype may account for the higher expression of SFRP1 in basal-like tumors compared to the other breast cancer subtypes and the associated increase in response to chemotherapy in triple-negative breast cancer compared to the other subtypes. These results suggest that SFRP1 signaling and methylation have potential as biomarkers tailored for the basal-like subtype, and may be used to predict chemotherapy response, as well as for treatment selection.

In addition, analysis of genes related with subtype-specific methylation revealed that CDH13 was specifically hypermethylated in HER2 tumors when compared with the other molecular subtypes. This result is in accordance with a previous study stating that the CDH13 gene is highly methylated in HER2/neu-positive breast cancers (22). The CDH13 gene, coding for H-cadherin, is a member of the cadherin family and a putative mediator of cell-cell interaction and cancer cell invasion and metastasis. It is considered as a tumor suppressor gene and may contribute to the enhancement of tumor progression and invasion (23). However, we did not observe a correlation between CDH13 methylation and stage, tumor grade, lymph node status and metastasis, which may reflect the limited sample size used in this study.

Several genes have previously been shown to be aberrantly methylated in breast cancer (8,24,25). Furthermore, breast cancer subtype-specific epigenotypes have been investigated through candidate gene approaches, as well as genome-wide DNA methylation analysis. However, these observations require further confirmation as the methylation frequency of a candidate gene and subtype-specific methylation patterns vary between independent studies. For example, Bediaga et al(9) and Holm et al(10) who used the array-based DNA methylation-profiling approach, found that each molecular subtype of breast cancer displayed specific methylation profiles (9,10). However, these studies did not identify and validate specific genes, such as SFRP1 and CDH13, whose methylation status was used to discriminate between the basal-like and HER2 subtypes in our study. Nevertheless, our findings are consistent with those of previous studies, stating that the HER2 subtype is associated with the preferential hypermethylation of several genes, and that basal-like tumors have several genes with low methylation levels (9,10,22,26). Taken together, our results support the findings that methylation may play a significant role in different breast tumor phenotypes.

In the present study, we used quantitative DNA methylation analysis to measure the methylation frequency in 5 genes known to be involved in breast carcinogenesis and are commonly methylated in breast cancer. Pyrosequencing is uniquely capable of quantifying methylation in explicit sequence context, thereby enabling several consecutive CpG sites to be quantified individually in a single assay. Although approaches for genome-wide DNA methylation analysis hold promise for identifying novel epigenetic targets, pyrosequencing is effective for identifying and quantifying the aberrant methylation of breast cancer genes and determining the criteria that may characterize and discriminate specific molecular subtypes.

To detect possible common patterns of methylation associated with the breast cancer molecular subtypes, hierarchical clustering of the 5 genes was performed. However, the discriminating ability of clustering analysis was poor compared with using the methylation status of the SFRP1 gene alone. This may be due to the relatively small sample size and the small number of candidate genes, which were insufficient to produce meaningful results (27). Furthermore, 3 of the 5 candidate genes, FHIT, Syk and RIZ1, displayed dissimilar methylation frequencies, and the methylation status of these genes may not be useful for the molecular classification of breast cancer.

Our study had several limitations. First, our study included a relatively small number of samples and candidate genes. These factors limit the application of our results to clinical settings. Further studies with a larger number of breast cancer cases are required to provide additional evidence. In addition, an analysis of additional candidate genes may help define subtype-specific methylation profiling for breast cancer classification. Second, we did not include normal breast tissue in the present study. Studies have shown that candidate genes are hypermethylated in breast tumors, whereas matched normal breast tissues have very low or no methylation levels (21,22,2830). However, direct comparisons of methylation levels in tumor tissue with those of normal tissue would be required to prove that epigenetic mechanisms may play a key role in the development of breast cancer. Furthermore, comparisons of gene expression with methylation levels in tumor and normal tissue would be helpful to address this issue further.

In conclusion, our study revealed that the basal-like subtype of breast cancer displays distinct methylation profiles compared to other subtypes. We found that the basal-like subtype is associated with low methylation levels of the SFRP1 gene, and that the hypermethylation of CDH13 differs among the molecular subtypes of breast cancer. These results suggest that the analysis of molecular markers, such as gene hypermethylation are useful for the improved characterization of molecular subtypes, and that the methylation levels of specific genes in breast cancer may potentially serve as epigenetic biomarkers and prognostic factors. Further studies with larger sample sizes and more comprehensive DNA methylation profiling with validation are required to clarify the predictive and prognostic value of gene methylation patterns in breast cancer.

Acknowledgements

The present research has been supported by Korea Breast Cancer Foundation.

References

1 

Sørlie T, Perou CM, Tibshirani R, et al: Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci USA. 98:10869–10874. 2001.PubMed/NCBI

2 

Sørlie T, Tibshirani R, Parker J, et al: Repeated observation of breast tumor subtypes in independent gene expression data sets. Proc Natl Acad Sci USA. 100:8418–8423. 2003.PubMed/NCBI

3 

Parker JS, Mullins M, Cheang MC, et al: Supervised risk predictor of breast cancer based on intrinsic subtypes. J Clin Oncol. 27:1160–1167. 2009. View Article : Google Scholar : PubMed/NCBI

4 

Nielsen TO, Hsu FD, Jensen K, et al: Immunohistochemical and clinical characterization of the basal-like subtype of invasive breast carcinoma. Clin Cancer Res. 10:5367–5374. 2004. View Article : Google Scholar : PubMed/NCBI

5 

Cheang MCU, Chia SK, Voduc D, et al: Ki67 index, HER2 status, and prognosis of patients with luminal B breast cancer. J Natl Cancer Inst. 101:736–750. 2009. View Article : Google Scholar

6 

Jaenisch R and Bird A: Epigenetic regulation of gene expression: how the genome integrates intrinsic and environmental signals. Nat Genet. 33:245–254. 2003. View Article : Google Scholar : PubMed/NCBI

7 

Dworkin AM, Huang TH and Toland AE: Epigenetic alterations in breast: Implications for breast cancer detection, prognosis and treatment. Semin Cancer Biol. 19:165–171. 2009. View Article : Google Scholar : PubMed/NCBI

8 

Widschwendter M and Jones PA: DNA methylation and breast carcinogenesis. Oncogene. 21:5462–5482. 2002. View Article : Google Scholar : PubMed/NCBI

9 

Bediaga N, Acha-Sagredo A, Guerra I, et al: DNA methylation epigenotypes in breast cancer molecular subtypes. Breast Cancer Res. 12:R772010. View Article : Google Scholar : PubMed/NCBI

10 

Holm K, Hegardt C, Staaf J, et al: Molecular sybtypes of breast cancer are associated with characteristic DNA methylation patterns. Breast Cancer Res. 12:R362010. View Article : Google Scholar : PubMed/NCBI

11 

Goldhirsch A, Wood WC, Coates AS, et al: Strategies for subtypes - dealing with the diversity of breast cancer: highlights of the St. Gallen International Expert Consensus on the Primary Therapy of Early Breast Cancer 2011. Ann Oncol. 22:1736–1747. 2011. View Article : Google Scholar : PubMed/NCBI

12 

Huang DW, Sherman BT and Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources. Nat Protoc. 4:44–57. 2009.PubMed/NCBI

13 

Wang S, Dorsey TH, Terunuma A, Kittles RA, Ambs S and Kwabi-Addo B: Relationship between tumor DNA methylation status and patient characteristics in African-American and European-American women with breast cancer. PLoS One. 7:e379272012.PubMed/NCBI

14 

Rattner A, Hsieh JC, Smallwood PM, et al: A family of secreted proteins contains homology to the cysteine-rich ligand-binding domain of frizzled receptors. Proc Natl Acad Sci USA. 94:2859–2863. 1997. View Article : Google Scholar : PubMed/NCBI

15 

Finch PW, He X, Kelley MJ, et al: Purification and molecular cloning of a secreted, Frizzled-related antagonist of Wnt action. Proc Natl Acad Sci USA. 94:6770–6775. 1997. View Article : Google Scholar : PubMed/NCBI

16 

Cadigan KM and Nusse R: Wnt signaling: a common theme in animal development. Genes Dev. 11:3286–3305. 1997. View Article : Google Scholar : PubMed/NCBI

17 

Polakis P: Wnt signaling and cancer. Genes Dev. 14:1837–1851. 2000.

18 

Ugolini F, Adélaïde J, Charafe-Jauffret E, et al: Differential expression assay of chromosome arm 8p genes identifies Frizzled-related (FRP1/FRZB) and Fibroblast Growth Factor Receptor 1 (FGFR1) as candidate breast cancer genes. Oncogene. 18:1903–1910. 1999. View Article : Google Scholar

19 

Klopocki E, Kristiansen G, Wild PJ, et al: Loss of SFRP1 is associated with breast cancer progression and poor prognosis in early stage tumors. Int J Oncol. 25:641–649. 2004.PubMed/NCBI

20 

Liedtke C, Ruckert C, Goette M, et al: Secreted frizzled receptor protein 1 (sFRP-1) as both a potential novel biomarker of triple negative breast cancer (TNBC), and its sensitivity against taxane/anthracycline containing neoadjuvant chemotherapy. Cancer Res. 69(24 Suppl): 40472009. View Article : Google Scholar

21 

Veeck J, Niederacher D, An H, Klopocki E, et al: Aberrant methylation of the Wnt antagonist SFRP1 in breast cancer is associated with unfavourable prognosis. Oncogene. 25:3479–3488. 2006. View Article : Google Scholar : PubMed/NCBI

22 

Fiegl H, Millinger S, Goebel G, et al: Breast cancer DNA methylation profiles in cancer cells and tumor stroma: association with HER-2/neu status in primary breast cancer. Cancer Res. 66:29–33. 2006. View Article : Google Scholar : PubMed/NCBI

23 

Celebiler CA, Kilic Y, Saydam S, et al: Predicting invasive phenotype with CDH1, CDH13, CD44, and TIMP3 gene expression in primary breast cancer. Cancer Sci. 100:2341–2345. 2009. View Article : Google Scholar : PubMed/NCBI

24 

Hinshelwood RA and Clark SJ: Breast cancer epigenetics: normal human mammary epithelial cells as a model system. J Mol Med. 86:1315–1328. 2008. View Article : Google Scholar : PubMed/NCBI

25 

Lo PK and Sukumar S: Epigenomics and breast cancer. Pharmacogenomics. 9:1879–1902. 2008. View Article : Google Scholar : PubMed/NCBI

26 

Feng W, Shen L, Wen S, et al: Correlation between CpG methylation profiles and hormone receptor status in breast cancers. Breast Cancer Res. 9:R572007. View Article : Google Scholar : PubMed/NCBI

27 

Eisen MB, Spellman PT, Brown PO and Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 95:14863–14868. 1998. View Article : Google Scholar : PubMed/NCBI

28 

Yuan Y, Mendez R, Sahin A and Dai JL: Hypermethylation leads to silencing of the SYK gene in human breast cancer. Cancer Res. 61:5558–5561. 2001.PubMed/NCBI

29 

Raish M, Dhillon VS, Ahmad A, et al: Promoter hypermethylation in tumor suppressing genes p16 and FHIT and their relationship with estrogen receptor and progesterone receptor status in breast cancer patients from Northern India. Transl Oncol. 2:264–270. 2009. View Article : Google Scholar

30 

Du Y, Carling T, Fang W, Piao Z, Sheu JC and Huang S: Hypermethylation in human cancers of the RIZ1 tumor suppressor gene, a member of a histone/protein methyltransferase superfamily. Cancer Res. 61:8094–8099. 2001.PubMed/NCBI

Related Articles

Journal Cover

May 2013
Volume 29 Issue 5

Print ISSN: 1021-335X
Online ISSN:1791-2431

Sign up for eToc alerts

Recommend to Library

Copy and paste a formatted citation
x
Spandidos Publications style
Jeong YJ, Jeong HY, Bong JG, Park SH and Oh HK: Low methylation levels of the SFRP1 gene are associated with the basal-like subtype of breast cancer. Oncol Rep 29: 1946-1954, 2013.
APA
Jeong, Y.J., Jeong, H.Y., Bong, J.G., Park, S.H., & Oh, H.K. (2013). Low methylation levels of the SFRP1 gene are associated with the basal-like subtype of breast cancer. Oncology Reports, 29, 1946-1954. https://doi.org/10.3892/or.2013.2335
MLA
Jeong, Y. J., Jeong, H. Y., Bong, J. G., Park, S. H., Oh, H. K."Low methylation levels of the SFRP1 gene are associated with the basal-like subtype of breast cancer". Oncology Reports 29.5 (2013): 1946-1954.
Chicago
Jeong, Y. J., Jeong, H. Y., Bong, J. G., Park, S. H., Oh, H. K."Low methylation levels of the SFRP1 gene are associated with the basal-like subtype of breast cancer". Oncology Reports 29, no. 5 (2013): 1946-1954. https://doi.org/10.3892/or.2013.2335