- Open Access
Investigation of DNA damage response and apoptotic gene methylation pattern in sporadic breast tumors using high throughput quantitative DNA methylation analysis technology
Molecular Cancer volume 9, Article number: 303 (2010)
Sporadic breast cancer like many other cancers is proposed to be a manifestation of abnormal genetic and epigenetic changes. For the past decade our laboratory has identified genes involved in DNA damage response (DDR), apoptosis and immunesurvelliance pathways to influence sporadic breast cancer risk in north Indian population. Further to enhance our knowledge at the epigenetic level, we performed DNA methylation study involving 17 gene promoter regions belonging to DNA damage response (DDR) and death receptor apoptotic pathway in 162 paired normal and cancerous breast tissues from 81 sporadic breast cancer patients, using a high throughput quantitative DNA methylation analysis technology.
The study identified five genes with statistically significant difference between normal and tumor tissues. Hypermethylation of DR5 (P = 0.001), DCR1 (P = 0.00001), DCR2 (P = 0.0000000005) and BRCA2 (P = 0.007) and hypomethylation of DR4 (P = 0.011) in sporadic breast tumor tissues suggested a weak/aberrant activation of the DDR/apoptotic pathway in breast tumorigenesis. Negative correlation was observed between methylation status and transcript expression levels for TRAIL, DR4, CASP8, ATM, CHEK2, BRCA1 and BRCA2 CpG sites. Categorization of the gene methylation with respect to the clinicopathological parameters showed an increase in aberrant methylation pattern in advanced tumors. These uncharacteristic methylation patterns corresponded with decreased death receptor apoptosis (P = 0.047) and DNA damage repair potential (P = 0.004) in advanced tumors. The observation of BRCA2 -26 G/A 5'UTR polymorphism concomitant with the presence of methylation in the promoter region was novel and emerged as a strong candidate for susceptibility to sporadic breast tumors.
Our study indicates that methylation of DDR-apoptotic gene promoters in sporadic breast cancer is not a random phenomenon. Progressive epigenetic alterations in advancing tumors result in aberrant DDR-apoptotic pathway thereby promoting tumor development. We propose, since pathological epigenetic changes of the DDR-apoptotic genes are reversible modifications, these could further be targeted for therapeutic interventions.
Breast cancer is the leading cause of cancer mortality among women aged between 20 - 59 years; second leading cause of cancer mortality among all women . The sporadic form represents almost 90% of the total number of breast cancer cases, genetic etiology of which is least understood and the molecular mechanism underlying the onset and progression not clear. It is, however, believed to be a manifestation of abnormal genetic as well as epigenetic changes [2–4] along with the influence of dietary, environmental and physical factors . Previous studies from our laboratory have identified genes involved in DNA damage response (DDR), apoptosis and immunesurvelliance pathways such as polymorphisms in BRCA2, p53 , IFNG , TGFB1 , TRAIL  and mDNA , somatic mutations in IL6  and mDNA , copy number variation in H2AX  and aberrant expression of DDR pathway  to influence sporadic breast cancer risk in north Indian population. Since, studies have identified de-novo methyltransferases, DNMT3b, over-expressed in breast tumors , indicating an involvement of epigenetic modifications in oncogenesis, breast cancer susceptibility genes, identified under the categories of DNA damage response (DDR) and apoptosis related genes could be a major target of epigenetic inactivation in sporadic breast cancers. Also, numerous studies have shown that large scale methylation profiling of multiple CpG sites could prove essential to provide comprehensive information on DNA methylation changes occurring during neoplastic transformation [16–24].
In cancer cells, apoptosis induced by extrinsic pathway complements the intrinsic pathway . The extrinsic signal transduction pathway is activated by TRAIL that binds two types of receptors: DR4/5 and DCR1/2 . Binding of TRAIL to DR4 and/or DR5 results in receptor oligomerization and subsequent activation of CASP8 resulting in apoptosis . Similarly, CASP8 can also activate the intrinsic apoptotic pathway causing the release of CYCS from the mitochondria which serves to amplify the death receptor apoptotic signal . On the contrary, decoy receptors, FLIP and BCL2 block the apoptotic signal transduction and promote survival [25, 27]. Recent studies have also identified the involvement of death receptors in activating the DNA damage response (DDR) pathway. The activation of ATM, CHEK2 and H2AX in response to TRAIL acts as a positive feedback loop involving the activation of caspases [28, 29]. Additionally, RNF8 ubiquitinates H2AX thus, enhancing its activation . H2AX, ATM and/or CHEK2 result in phosphorylation and activation of P53 and regulating DNA repair as well as apoptosis . Although P53 is not essential for TRAIL mediated apoptosis [32, 33], it can cause transcriptional induction of pro-apoptotic genes such as TRAIL, DR4, DR5[34–36] and antagonize pro-survival genes, FLIP and BCL2[37, 38]. Further, acetylation of P53 by TIP60, also known as KAT5 promotes apoptosis by preferentially transactivating pro-apoptotic genes . Moreover, DNA damage response also results in activation of BRCA1 and BRCA2 that are involved in the repair of damaged DNA  (Figure 1).
Although, several methylation studies have been carried out for DR4, DR5, DCR1, DCR2, CASP8, TP53, BRCA1 and BRCA2 in various types of cancer, none of them provide a holistic view of the DDR - apoptotic pathway [21, 41–45]. Therefore, after investigating the breast tumor at the genetic and expression level, we for the first time provide an integrated information on methylation patterns of the DDR and death receptor apoptotic pathway genes (TRAIL, DR4, DR5, DCR1, DCR2, CASP8, FLIP, BCL2, CYCS, ATM, TP53, BRCA1, BRCA2, H2AX, RNF8, TIP60 and CHEK2) in sporadic breast tumors, using high throughput automated MALDI TOF [17–19, 22]. The study further dissects the role of methylation of candidate genes in sporadic breast tumorigenesis, its interaction with functional BRCA2 -26 SNP and the status of pro- and anti- apoptotic as well as DNA damage repair gene expression in breast tumor tissues.
Materials and methods
A total of 162 tissue samples (81 ductal carcinoma tissues and 81 adjacent normal tissues, both of parenchymal origin) were collected from 81 patients with sporadic breast cancer from Dharamshila Cancer Hospital and Research Centre, New Delhi; Rajiv Gandhi Cancer Institute and Research Center, New Delhi and Sanjay Gandhi Postgraduate Institute of Medical Sciences, Lucknow. The study samples adhered to the REMARK guidelines . The female patients ranged in the age group of 25 to 77 years, with a median of 48 years. None of the studied cases had a hereditary form of breast cancer. Prior approval was obtained from Jawaharlal Nehru University ethical committee and the informed consent taken of the concerned subjects for sample collection and study.
Differentiation between cancerous and normal tissue was based on magnetic resonance imaging, intra-operative gross surgical pathology and tissue histology. Breast tumor patients underwent imaging (mammography, ultra-sonography and MRI) to know the location, size and extent of tumor. Tumor was removed wide of its margin i.e. 1-2 cm envelope of normal breast tissue was left all around the tumor to ensure removal of any microscopic extensions of the tumor into normal breast parenchyma. What we considered as normal tissue was taken 2 cm away from the palpable tumor margin. This fact was further confirmed by frozen section histology and paraffin section histology before labeling the tumor and normal tissues for molecular biology experiments, thereby ensuring the absence of cross contamination.
Patients were classified on the basis of tumor size, nodal status, tumor stage, estrogen and progesterone receptor (ER and PR) status. At the time of diagnosis, 7 patients had stage I disease, 44 patients belonged to stage II, 28 patients to stage III and 2 patients had stage IV disease. ER/PR status and tumor grade could be obtained for 57 out of 81 pairs of tumor tissues studied. Immunohistochemical staining identified 27/24 ER/PR positive patients and 30/33 ER/PR negative patients. Out of 57 patients, 2 patients had grade I tumor, 30 had grade II tumor and 25 had grade III tumor (Additional file 1, Table S1). The samples collected were frozen immediately and stored at -80°C until use. DNA extraction was performed from 0.01 - 0.02 g of tissue sample, using the Genelute Mammalian genomic DNA isolation kit (Sigma Aldrich, St. Louis, Missouri, USA).
Bisulfite treatment and PCR
The EZ-96 DNA Methylation Kit (Zymo Research, Orange, CA, USA) was used for bisulfite conversion of the target sequences. The C/T conversion reaction was performed using the PCR program as follows: 95°C for 30 sec and 50°C for 15 min, which was repeated for 46 cycles. Primers were designed using epidesigner software (Sequenom, San Diego, CA, USA) to cover the regions with the most CpG sites (Additional file 2, Table S2). Our selected amplicons were mostly located in the promoter region of genes or started from the promoter and ended in the first exon. In PCR amplification, a T7-promoter tag was added to the reverse primer and a 10 mer-tag sequence was added to the forward primer to balance the PCR primer length. The bisulfite treated genomic DNA was amplified using Taq DNA polymerase (Roche Diagnostics, Mannheim, Germany) (4 min at 95°C followed by 45 cycles of 20 sec at 95°C, 30 sec at 62°C, and 1 min at 72°C with a 3 minute final extension). PCR products were analyzed further in Sequenom MassARRAY. Methylated and unmethylated positive control human DNA was procured from Sequenom. Fully methylated DNA was mixed with pure unmethylated DNA in a ratio of 100:0, 60:40, 40:60, and 0:100.
In vitro transcription, T cleavage assay and Mass Spectrometry
Unincorporated dNTPs were dephosphorylated by adding 1.7 μl H2O and 0.3U shrimp alkaline phosphatase (Sequenom) followed by incubation at 37°C for 40 min and heat inactivation of shrimp alkaline phosphatase at 85°C for 5 min. In general, 2 μl of the PCR was directly used as a template in a 5 μl transcription reaction. T7 RNA and DNA Polymerase (20U) (Sequenom) was used to incorporate dTTP in the transcripts. Ribonucleotides were used at 1 mmol/l and the dNTP substrate at 2.5 mmol/l. RNaseA enzyme (Sequenom) was added in the same step to cleave the in vitro transcripts (T-cleavage assay). The T cleavage reaction was carried out at 37°C for 3 hours and further diluted with H2O to a final volume of 27 μl. Conditioning of the phosphate backbone was achieved by adding 6 mg of Clean Resin (Sequenom) before performing MALDI-TOF MS.
22 nl of cleavage reaction was robotically dispensed onto a silicon matrix preloaded chips (SpectroCHIP; Sequenom), and the mass spectra obtained using a MassARRAY. Methylation ratios were generated by the MALDI-TOF and EpiTYPER software v1.0 (Sequenom). The assay was able to discriminate between the methylated and unmethylated components of the positive control according to the ratios.
Genotype and Expression Analysis
Sequence based analysis of the amplified 5'UTR region of the BRCA2 gene as described earlier  for genotype status was performed for 81 samples in which methylation study was carried out. Commercially available Taqman Gene expression Assay system for quantitating transcript level of TRAIL, DR4, DR5, DCR1, DCR2, CASP8, CASP8L, FLIPL, FLIPS, BCL2, CYCS, ATM, TP53, BRCA1, BRCA2, CHEK2[9, 14] and H2AX (Hs01573336_s1) genes belonging to the DDR and apoptotic pathway were used for studying mRNA expression in 40 representative tumor samples out of the 81 paired normal/tumor samples used in methylation study (Applied Biosystems, Foster City, CA, USA). GAPDH, B-Actin, PUM1, and MRPL19 (Applied Biosystems) were used as endogenous controls . Quantitative real-time PCR was carried out using an ABI Prism 7900 Sequence Detection System (Applied Biosystems). Threshold cycle (Ct) numbers were established by using SDS 1.1 RQ software (Applied Biosystems). All the reactions were carried out in duplicates. The normalization factor obtained from GeNorm software was used to compute normalized expression for the target genes .
Categorization of the results was done on the basis of pro- versus anti- apoptotic genes involved in death receptor apoptotic pathway [(TRAIL + DR4 + DR5 + CASP8 + ATM + H2AX + CHEK2 + CYCS)/(DCR1 + DCR2 + CASP8L + FLIPL + FLIPS + BCL2)] and DNA damage response pathway [ATM + CHEK2 + BRCA1 + BRCA2 + H2AX]. Earlier studies from our laboratory have shown that higher level of P53 expression in advanced breast tumors did not reflect in P53 activity in later stages of tumor development  and thus was not considered while calculating the DNA damage response signal.
Relative methylation was compared using the Wilcoxon signed-rank test, a nonparametric counterpart of the paired t-test to identify sites with statistically significant difference, between cancerous and paired normal tissues. Using the two-way hierarchical cluster analysis, the CpG fragments for each gene were clustered based on pair-wise Euclidean distances and linkage algorithm for all of the 81 pairs of tissue samples. Clustering of CpG units were plotted along the y axis and samples along the x axis. The procedure was carried out using the heatmap.2 function of the 'gplots' package using the R statistical software. Comparison of methylation between tumor and normal group was performed using Mann-Whitney U test for two groups http://udel.edu/~mcdonald/statkruskalwallis.html. Fisher's exact test and Spearman's correlation coefficient was calculated using the SPSS statistical package, version 13 (SPSS Inc., Chicago, IL, USA). P value was considered significant at and below ≤ 0.05.
Semiquantitative methylation profiles of 17 genes in cancerous and normal tissue
In this study, we analyzed the methylation patterns of 17 breast cancer genes in 162 (cancerous and normal) breast tissues from 81 breast tumor patients. For the 17 genes, one amplicon per gene was analyzed and the median amplicon length was 427 bp (range = 265-496 bp). In total we assessed 492 CpGs per sample with a total of 79,704 CpGs in 81 pairs of samples. We removed those CpG sites that could not be detected by MALDI-TOF MS and did not yield successful measurements because of generation of higher and lower mass fragments. CpGs which gave good results for > 90% of the samples were considered; these consisted of 227 CpGs per sample (139 CpG sites or units) equivalent to 36,774 CpGs in 81 pairs of samples equating to 46.14% of the total CpGs (Additional file 3, Table S3).
In normal tissues, 59% of all CpGs showed a mean value of < 5% methylation; 25% showed 5%-10% methylation and the remaining 16% displayed 10%-30% methylation. However, none of the normal tissues studied showed > 30% methylation, suggesting a lower level of methylation in our study group for the candidate genes (Figure 2a). This was in agreement with previous finding for normal tissue samples . Mann-Whitney U test showed average methylation to be significantly higher in tumors whose paired counterpart belonged to the mean methylation group of 5%-10% (P = 0.009). Moreover, average methylation turned out to be highly significant on comparison of whole sporadic breast tumor group with the normal group (P = 6.272 × 10-20) (Figure 2b).
Wilcoxon signed rank test followed by Benjamini Hochberg correction for individual CpG sites identified the sites that showed statistically significant difference between normal and tumor samples (P ≤ 0.05) (Table 1). Tumor tissue showed hypomethylation for DR4, FLIP, and RNF8 and hypermethylation for DR5, DCR1, DCR2, CASP8, CYCS, BRCA1, BRCA2, whereas H2AX, TRAIL, BCL2, ATM, CHEK2, TP53 and TIP60 did not show any significant difference between the breast tumor and the corresponding normal tissues. Categorization of the significant CpG units with respect to the transcription start site (TSS) revealed their presence in the promoter, 5'UTR, coding exon 1 and/or intron1-2 regions.
Further, unsupervised two-way hierarchical clustering with 139 CpG sites/units in 81 pair of samples showed separate grouping of normal and breast tumor samples with negligible intermixing (Figure 3) as confirmed by Fisher's exact test (P = 0.00000001) (Table 2). We again performed unsupervised two-way hierarchical clustering for the 33 CpG units that were found to be significant after supervised clustering in Wilcoxon test and Benjamini Hochberg correction. An overall low methylation was observed in the normal samples as compared to tumor samples (Additional file 4, Figure S1).
Mean methylation values (sample-wise average of methylation values of all CpG sites in one gene) reflected that DR5, DCR1, DCR2 and BRCA2 were hypermethylated whereas DR4 was hypomethylated in tumors as compared to adjacent normal tissues (Figure 4). Genes that did not show any difference between mean methylation pattern of normal and tumor tissues were: TRAIL, BCL2, CASP8, CYCS, FLIP, ATM, TP53, CHEK2, RNF8, TIP60, H2AX and BRCA1.
In order to understand the role of methylation in tumor progression, we compared the differential methylation pattern observed in normal versus tumor samples with respect to the various clinicopathological parameters (Table 3). ER negative tumors showed hypermethylation of DCR2 (87%; P = 0.00007) and hypomethylation of DR4 (77%; P = 0.001). Hypermethylation of DCR2 (75%; p = 0.001) in PR positive cases, hypomethylated DR4 (79%; p = 0.0002), hypermethylated DCR1 (64%; P = 0.004) and DCR2 (85%; P = 0.00003) in PR negative tumors was observed even after Benjamini Hochberg correction. Node negative tumors showed hypermethylation of DCR1 (63%; P = 0.003) and DCR2 (83%; P = 0.00001) while node positive displayed hypermethylation of DCR1 (66%, P = 0.001), DCR2 (73%; P = 0.00001) and BRCA2 (63%; P = 0.001). Collectively, tumor size 1 and 2 showed hypermethylation of DR5 (68%; P = 0.002), DCR1 (63%; P = 0.001) and DCR2 (80%; P = 0.0000001) as well as hypomethylation of DR4 (58%; P = 0.011), FLIP (61%; P = 0.007) and ATM (59%; P = 0.013) after Benjamini Hochberg correction. Similarly, tumor size 3 and 4 in combination showed hypermethylation of DCR1 (68%; P = 0.002) and DCR2 (73%; P = 0.001). Together, stage I and II tumors displayed hypermethylation of DCR1 (61%; P = 0.003) and DCR2 (82%; P = 0.000001), although hypomethylation was observed in some of the candidates but lost after Benjamini Hochberg correction. Likewise, stage III and IV tumors jointly showed hypermethylation in DCR1 (70%; P = 0.001) and DCR2 (70%; P = 0.0002) and BRCA2 (57%; P = 0.005). Grade I and II tumors exhibited hypermethylation in DCR2 (78%; P = 0.0002) while grade III tumors displayed DR5 (64%; P = 0.006), DCR1 (68%; P = 0.007) and DCR2 (84%; P = 0.0001) hypermethylation, and hypomethylation of DR4 (84%; P = 0.0003) and RNF8 (72%; P = 0.010).
Real time expression analysis of DDR - apoptotic pathway genes and their correlation with the methylation pattern
In order to understand the level of DNA damage repair and apoptosis with respect to tumor progression, a representative set of 40 tumor samples were studied for the real time expression of 15 genes (TRAIL, DR4, DR5, DCR1, DCR2, CASP8, CASP8L, FLIPL, FLIPS, BCL2, CYCS, ATM, TP53, BRCA1, BRCA2, CHEK2 and H2AX) belonging to the DDR and apoptotic pathway. The results of expression for these genes were categorized in pro- apoptotic versus anti- apoptotic signals and DNA damage response signals with respect to the tumor stages. The results showed that DNA damage response as well as apoptosis decreased with advancing sporadic breast tumors (P = 0.004; 0.047) (Figure 5) (Additional file 5, Figure S2).
Also, the relationship between promoter methylation and gene expression was evaluated. Spearman's correlation coefficient was calculated for each CpG unit. Negative correlation was observed for TRAIL -116 (P = 0.013; ρ = -0.350), DR4 +108 (P = 0.025; ρ = -0.331), CASP8 +571 (P = 0.046; ρ = -0.271), ATM -369 (P = 0.046; ρ = -0.270), ATM -314, -310 (P = 0.034; ρ = -0.291), CHEK2 -396, -393 (P = 0.040; ρ = -0.281), BRCA1 -314 (P = 0.030; ρ = -0.300), BRCA1 -251 (P = 0.004; ρ = -0.410), BRCA2 -101 (P = 0.036; ρ = -0.288) and BRCA2 +133, +136, +138 (P = 0.023; ρ = -0.318). However, positive correlation was observed for DR4 +23, +25, +28, +32 (P = 0.017; ρ = 0.335) and CYCS +241, +243 (P = 0.048; ρ = 0.266) (Table 4). On other hand, no correlation was observed for DR5, FLIP, BCL2, TP53 and H2AX. Insilico analysis using Alibaba2.1 revealed the presence of Sp1, ER and GBF2 transcription factor (TF) binding sites or absence of TF consensus sequence at these positions.
Involvement of BRCA2 in sporadic breast tumor
Earlier studies from our laboratory have identified BRCA2 -26 wild type GG and mutant AA genotype to provide risk whereas heterozygote GA to provide protection against sporadic breast tumors . Categorization of breast tumor samples on the basis of BRCA2 -26 G/A 5'UTR polymorphism, followed by methylation status of the candidate genes resulted in the identification of GG and AA genotypes to be associated with hypermethylated BRCA2 (63%;P = 0.008) , DR5 (69%;P = 0.001), DCR1 (67%; P = 0.0001), DCR2 (76%;P = 0.0000002), CASP8 (59%;P = 0.006) and hypomethylated FLIP (59%;P = 0.009). Whereas, the GA protector genotype was found to be associated with hypermethylated DCR2 (81%; P = 0.0008) (Table 5). Expression data for BRCA2 transcript compared with BRCA2 mean methylation level, revealed a negative correlation (P = 0.033; ρ = -293) (Figure 6).
Methylation alterations in cancer have been recognized for decades [47–50]. We analyzed quantitative changes in methylation of 17 promoter regions in 162 paired normal and cancerous breast tissues from 81 sporadic breast cancer patients using high throughput MALDI-TOF MS and evaluated their distribution, correlation and relationships to clinicopathological variables using common statistical methods. The genes selected belonged to the DNA damage response (DDR) and death receptor apoptotic pathways [30, 51–54]. Among the studied 17 genes, this is the first report of promoter methylation for CpG sites in H2AX, RNF8 and CYCS in human cancer. The limitation of the study was our ability to score only 46.14% of the total CpG of the studied region and the absence of screening of CpG site methylation for the entire promoter and the 5'UTR.
Although methylation levels observed in breast tissues were not very high, yet a significantly differential methylation pattern existed between normal and tumor samples, suggesting that the epigenetic control of the gene promoters was required for the maintenance of normal cellular homeostasis, deregulation of which could result in tumor development. Based on the DNA damage response and apoptotic pathway gene methylation profile, we were able to segregate the normal and tumor tissues. Our observation followed by statistical analysis of methylation status in breast tumors along with transcript expression revealed a negative correlation for TRAIL, DR4, CASP8, ATM, CHEK2, BRCA1 and BRCA2 CpG sites. Insilico analysis of these methylated CpG sites/units revealed the presence of consensus sequence for stimulatory protein 1 (Sp1) and estrogen receptor (ER) transcription factors. Studies have revealed that methylated CpGs in these recognition sites preclude the binding of the Sp1 transcription factor and thereby inhibit gene expression directly [55, 56]. DNA methylation has also been shown to be followed by binding of methyl-CpG binding domain proteins (MBDs) which contribute to gene repression by the recruitment of histone deacetylases (HDACs) to nucleosomes . The CpG sites within CHEK2 and BRCA2 were observed to be the consensus sequence for Sp1; whereas such sites in TRAIL, DR4, CASP8, ATM and BRCA1 did not have transcription factor binding sites. We believe methylation of these sites might inhibit/promote the binding of Sp1/MBDs, respectively to the methylated DNA, resulting in decreased gene expression.
For the first time, our study showed a collaborative involvement of DR5, DCR1, DCR2, CASP8, CYCS, BRCA1, BRCA2, H2AX hypermethylation and DR4, FLIP, RNF8 hypomethylation in sporadic breast tumor pathogenesis. Since apoptotic signaling through DR5 has been reported to be more potent than through DR4 , our data suggests that hypermethylation status of DR5 receptor facilitates tumor cells to evade apoptosis. Despite the hypomethylation state of DR4 receptor and hypermethylation state of anti-apoptotic decoy receptors (DCR1 and DCR2) [25, 42, 43], we propose that the final outcome in tumor pathogenesis depends on downstream signal transduction molecules, such as CASP8, FLIP and CYCS in established tumors. Hypermethylation of CASP8 and hypomethylation of FLIP inhibits the extrinsic apoptotic pathway, whereas hypermethylation of CYCS suggests a weak activation of intrinsic apoptotic pathway. Similarly, hypermethylation of BRCA1, BRCA2 and H2AX implies a decreased DNA damage repair. This is in accordance with the knowledge that inhibition of DNA damage response and apoptosis may contribute to tumor initiation, growth and metastasis in the pathogenesis of breast cancer [59, 60].
We believe decoy receptor down-regulation through hypermethylation in initial (stages I and II) as well as advanced stages (stage III and IV) represents a 'physiological' response to check the conversion of precancerous cell to a cancerous one by enhancing apoptosis, which can be evaded by established tumors via mechanisms such as aberrant methylation of DDR and apoptotic genes. The methylation profiles of the candidate genes in this pathway revealed increased apoptosis in initial tumors (T1 and T2) as a result of hypomethylated DR4 which initiates the death receptor signal transduction pathway in response to TRAIL. Although in initial tumors (T1 and T2) hypermethylation of DR5 and hypomethylation of FLIP is observed, we believe, the positive feedback loop due to hypomethylated ATM is strong enough to override the antiapoptotic effect of FLIP thereby enhancing DR4 mediated apoptosis. However, in advanced tumors (T3 and T4) most of the methylation profiles observed in initial tumor stages are lost consequently, resulting in decreased apoptosis. Hypermethylation of DR5 in grade III tumors also suggests a reduced apoptotic potential in more aggressive form of breast tumors. In the context of TRAIL mediating the activation of DNA damage response pathway [28, 29], hypermethylation of BRCA2 (stage III and IV) indicated an aberrant DNA damage repair mechanism in response to TRAIL mediated apoptosis. Thus our study showed a deregulation of death receptor apoptosis and DDR pathway genes through methylation, influencing breast tumorigenesis; also indicated by reduced apoptosis and DNA repair in advanced tumor. This is in accordance with the fact that suppression of apoptosis and DNA damage repair comprise the minimum common platform upon which all neoplastic evolution occurs . We suggest that besides genetic aberrations of DDR-apoptotic genes, their deregulated methylation also plays an intricate role in tumor development. On the whole, the study revealed cooperative involvement of the DDR and apoptotic genes at the level of methylation in assisting tumor cell survival and progression by apoptotic evasion and loss of DNA repair during sporadic breast tumor progression (Figure 7).
In addition, BRCA2 germline variation concomitant with the presence of methylation in the promoter region was novel and interesting and emerged as a strong candidate for susceptibility to sporadic breast tumors. Flanagan et.al.  showed in familial breast cancers that methylation profiles are defined by mutation status. On the basis of BRCA2 -26 GA heterozygotes observed earlier by us with a protector genotype, showing decreased LOH; and the risk providing GG/AA homozygotes showing increased LOH , we suggest that the presence of hypermethylated DR5, CASP8, BRCA2 and hypomethylated FLIP with BRCA2 -26 GG or AA genotype, might cause increased LOH, decreased DNA repair and apoptosis, providing a risk for sporadic breast tumor development and progression. These observations allow us to assign an important role to a concomitant presence of variant germline genetic background at 5'UTR of BRCA2 and the epigenetic modification in the process of oncogenesis.
No genome-wide evaluation have been carried out to identify altered DNA methylation patterns in the context of tumor initiation and/or progression  and their effects on DNA damage response - apoptotic pathway or gene networks. Our study indicates that promoter methylation of DDR-apoptotic genes in sporadic breast cancer is not a random phenomenon. It has two features: hypermethylation of DDR-apoptotic genes, and hypomethylation of anti-apoptotic/pro-survival genes. Progressive modification of the aberrant epigenetic alterations with advancing tumors results in deregulation of the DDR-apoptotic pathway thereby promoting tumor development. We also observe BRCA2 to play a major role at both genetic and epigenetic level in sporadic breast tumor pathogenesis. We propose, since pathological epigenetic changes of the DDR-apoptotic genes are reversible modifications , these could further be targeted for therapeutic interventions.
Ataxia Telangiectasia Mutated
B cell lymphoma 2
Breast Cancer 1 and 2, early onset proteins
Checkpoint Kinase 2
Decoy Receptors 1/2
Death Receptors 4/5
FLICE like inhibitory protein large/small
Histone H2A-member X
Ring finger protein 8
Lysine acetyltransferase 5
Tumor necrosis factor related apoptosis inducing ligand
Xu X, Gammon MD, Zhang Y, Bestor TH, Zeisel SH, Wetmur JG, Wallenstein S, Bradshaw PT, Garbowski G, Teitelbaum SL: BRCA1 promoter methylation is associated with increased mortality among women with breast cancer. Breast Cancer Res Treat. 2009, 115: 397-404. 10.1007/s10549-008-0075-5
Russo J, Yang X, Hu YF, Bove BA, Huang Y, Silva ID, Tahin Q, Wu Y, Higgy N, Zekri A, Russo IH: Biological and molecular basis of human breast cancer. Front Biosci. 1998, 3: D944-960.
Feinberg AP, Tycko B: The history of cancer epigenetics. Nat Rev Cancer. 2004, 4: 143-153. 10.1038/nrc1279
Baylin SB, Ohm JE: Epigenetic gene silencing in cancer - a mechanism for early oncogenic pathway addiction?. Nat Rev Cancer. 2006, 6: 107-116. 10.1038/nrc1799
Oldenburg RA, Meijers-Heijboer H, Cornelisse CJ, Devilee P: Genetic susceptibility for breast cancer: how many more genes to be found?. Crit Rev Oncol Hematol. 2007, 63: 125-149. 10.1016/j.critrevonc.2006.12.004
Gochhait S, Bukhari SI, Bairwa N, Vadhera S, Darvishi K, Raish M, Gupta P, Husain SA, Bamezai RN: Implication of BRCA2 -26G > A 5' untranslated region polymorphism in susceptibility to sporadic breast cancer and its modulation by p53 codon 72 Arg > Pro polymorphism. Breast Cancer Res. 2007, 9: R71- 10.1186/bcr1780
Saha A, Dhir A, Ranjan A, Gupta V, Bairwa N, Bamezai R: Functional IFNG polymorphism in intron 1 in association with an increased risk to promote sporadic breast cancer. Immunogenetics. 2005, 57: 165-171. 10.1007/s00251-005-0783-5
Saha A, Gupta V, Bairwa NK, Malhotra D, Bamezai R: Transforming growth factor-beta1 genotype in sporadic breast cancer patients from India: status of enhancer, promoter, 5'-untranslated-region and exon-1 polymorphisms. Eur J Immunogenet. 2004, 31: 37-42. 10.1111/j.1365-2370.2004.00442.x
Pal R, Gochhait S, Chattopadhyay S, Gupta P, Prakash N, Agarwal G, Chaturvedi A, Husain N, Husain SA, Bamezai RN: Functional implication of TRAIL -716 C/T promoter polymorphism on its in vitro and in vivo expression and the susceptibility to sporadic breast tumor. Breast Cancer Res Treat. 2010
Darvishi K, Sharma S, Bhat AK, Rai E, Bamezai RN: Mitochondrial DNA G10398A polymorphism imparts maternal Haplogroup N a risk for breast and esophageal cancer. Cancer Lett. 2007, 249: 249-255. 10.1016/j.canlet.2006.09.005
Saha A, Bairwa NK, Ranjan A, Gupta V, Bamezai R: Two novel somatic mutations in the human interleukin 6 promoter region in a patient with sporadic breast cancer. Eur J Immunogenet. 2003, 30: 397-400. 10.1111/j.1365-2370.2003.00423.x
Gochhait S, Bhatt A, Sharma S, Singh YP, Gupta P, Bamezai RN: Concomitant presence of mutations in mitochondrial genome and p53 in cancer development - a study in north Indian sporadic breast and esophageal cancer patients. Int J Cancer. 2008, 123: 2580-2586. 10.1002/ijc.23817
Srivastava N, Gochhait S, Gupta P, Bamezai RN: Copy number alterations of the H2AFX gene in sporadic breast cancer patients. Cancer Genet Cytogenet. 2008, 180: 121-128. 10.1016/j.cancergencyto.2007.09.024
Gochhait S, Dar S, Pal R, Gupta P, Bamezai RN: Expression of DNA damage response genes indicate progressive breast tumors. Cancer Lett. 2009, 273: 305-311. 10.1016/j.canlet.2008.08.009
Veeck J, Esteller M: Breast cancer epigenetics: from DNA methylation to microRNAs. J Mammary Gland Biol Neoplasia. 2010, 15: 5-17. 10.1007/s10911-010-9165-1
Stanssens P, Zabeau M, Meersseman G, Remes G, Gansemans Y, Storm N, Hartmer R, Honisch C, Rodi CP, Bocker S, van den Boom D: High-throughput MALDI-TOF discovery of genomic sequence polymorphisms. Genome Res. 2004, 14: 126-133. 10.1101/gr.1692304
Ehrich M, Nelson MR, Stanssens P, Zabeau M, Liloglou T, Xinarianos G, Cantor CR, Field JK, van den Boom D: Quantitative high-throughput analysis of DNA methylation patterns by base-specific cleavage and mass spectrometry. Proc Natl Acad Sci USA. 2005, 102: 15785-15790. 10.1073/pnas.0507816102
Ehrich M, Turner J, Gibbs P, Lipton L, Giovanneti M, Cantor C, van den Boom D: Cytosine methylation profiling of cancer cell lines. Proc Natl Acad Sci USA. 2008, 105: 4844-4849. 10.1073/pnas.0712251105
Radpour R, Kohler C, Haghighi MM, Fan AX, Holzgreve W, Zhong XY: Methylation profiles of 22 candidate genes in breast cancer using high-throughput MALDI-TOF mass array. Oncogene. 2009, 28: 2969-2978. 10.1038/onc.2009.149
Radpour R, Haghighi MM, Fan AX, Torbati PM, Hahn S, Holzgreve W, Zhong XY: High-throughput hacking of the methylation patterns in breast cancer by in vitro transcription and thymidine-specific cleavage mass array on MALDI-TOF silico-chip. Mol Cancer Res. 2008, 6: 1702-1709. 10.1158/1541-7786.MCR-08-0262
Flanagan JM, Cocciardi S, Waddell N, Johnstone CN, Marsh A, Henderson S, Simpson P, da Silva L, Khanna K, Lakhani S: DNA methylome of familial breast cancer identifies distinct profiles defined by mutation status. Am J Hum Genet. 2010, 86: 420-433. 10.1016/j.ajhg.2010.02.008
Ehrich M, Field JK, Liloglou T, Xinarianos G, Oeth P, Nelson MR, Cantor CR, van den Boom D: Cytosine methylation profiles as a molecular marker in non-small cell lung cancer. Cancer Res. 2006, 66: 10911-10918. 10.1158/0008-5472.CAN-06-0400
Ordway JM, Budiman MA, Korshunova Y, Maloney RK, Bedell JA, Citek RW, Bacher B, Peterson S, Rohlfing T, Hall J: Identification of novel high-frequency DNA methylation changes in breast cancer. PLoS One. 2007, 2: e1314- 10.1371/journal.pone.0001314
Holm K, Hegardt C, Staaf J, Vallon-Christersson J, Jonsson G, Olsson H, Borg A, Ringner M: Molecular subtypes of breast cancer are associated with characteristic DNA methylation patterns. Breast Cancer Res. 2010, 12: R36- 10.1186/bcr2590
Johnstone RW, Frew AJ, Smyth MJ: The TRAIL apoptotic pathway in cancer onset, progression and therapy. Nat Rev Cancer. 2008, 8: 782-798. 10.1038/nrc2465
Falschlehner C, Emmerich CH, Gerlach B, Walczak H: TRAIL signalling: decisions between life and death. Int J Biochem Cell Biol. 2007, 39: 1462-1475. 10.1016/j.biocel.2007.02.007
Ashkenazi A: Directing cancer cells to self-destruct with pro-apoptotic receptor agonists. Nat Rev Drug Discov. 2008, 7: 1001-1012. 10.1038/nrd2637
Solier S, Pommier Y: The apoptotic ring: a novel entity with phosphorylated histones H2AX and H2B and activated DNA damage response kinases. Cell Cycle. 2009, 8: 1853-1859.
Solier S, Sordet O, Kohn KW, Pommier Y: Death receptor-induced activation of the Chk2- and histone H2AX-associated DNA damage response pathways. Mol Cell Biol. 2009, 29: 68-82. 10.1128/MCB.00581-08
Harper JW, Elledge SJ: The DNA damage response: ten years after. Mol Cell. 2007, 28: 739-745. 10.1016/j.molcel.2007.11.015
Meek DW: Tumour suppression by p53: a role for the DNA damage response?. Nat Rev Cancer. 2009, 9: 714-723.
Ravi R, Bedi A: Sensitization of tumor cells to Apo2 ligand/TRAIL-induced apoptosis by inhibition of casein kinase II. Cancer Res. 2002, 62: 4180-4185.
Ray S, Shyam S, Fraizer GC, Almasan A: S-phase checkpoints regulate Apo2 ligand/TRAIL and CPT-11-induced apoptosis of prostate cancer cells. Mol Cancer Ther. 2007, 6: 1368-1378. 10.1158/1535-7163.MCT-05-0414
Liu X, Yue P, Khuri FR, Sun SY: p53 upregulates death receptor 4 expression through an intronic p53 binding site. Cancer Res. 2004, 64: 5078-5083. 10.1158/0008-5472.CAN-04-1195
Takimoto R, El-Deiry WS: Wild-type p53 transactivates the KILLER/DR5 gene through an intronic sequence-specific DNA-binding site. Oncogene. 2000, 19: 1735-1743. 10.1038/sj.onc.1203489
Kuribayashi K, Krigsfeld G, Wang W, Xu J, Mayes PA, Dicker DT, Wu GS, El-Deiry WS: TNFSF10 (TRAIL), a p53 target gene that mediates p53-dependent cell death. Cancer Biol Ther. 2008, 7: 2034-2038.
Fukazawa T, Fujiwara T, Uno F, Teraishi F, Kadowaki Y, Itoshima T, Takata Y, Kagawa S, Roth JA, Tschopp J, Tanaka N: Accelerated degradation of cellular FLIP protein through the ubiquitin-proteasome pathway in p53-mediated apoptosis of human cancer cells. Oncogene. 2001, 20: 5225-5231. 10.1038/sj.onc.1204673
Zhivotovsky B, Kroemer G: Apoptosis and genomic instability. Nat Rev Mol Cell Biol. 2004, 5: 752-762. 10.1038/nrm1443
Sykes SM, Mellert HS, Holbert MA, Li K, Marmorstein R, Lane WS, McMahon SB: Acetylation of the p53 DNA-binding domain regulates apoptosis induction. Mol Cell. 2006, 24: 841-851. 10.1016/j.molcel.2006.11.026
Sancar A, Lindsey-Boltz LA, Unsal-Kacmaz K, Linn S: Molecular mechanisms of mammalian DNA repair and the DNA damage checkpoints. Annu Rev Biochem. 2004, 73: 39-85. 10.1146/annurev.biochem.73.011303.073723
Yang Q, Kiernan CM, Tian Y, Salwen HR, Chlenski A, Brumback BA, London WB, Cohn SL: Methylation of CASP8, DCR2, and HIN-1 in neuroblastoma is associated with poor outcome. Clin Cancer Res. 2007, 13: 3191-3197. 10.1158/1078-0432.CCR-06-2846
van Noesel MM, van Bezouw S, Salomons GS, Voute PA, Pieters R, Baylin SB, Herman JG, Versteeg R: Tumor-specific down-regulation of the tumor necrosis factor-related apoptosis-inducing ligand decoy receptors DcR1 and DcR2 is associated with dense promoter hypermethylation. Cancer Res. 2002, 62: 2157-2161.
Shivapurkar N, Toyooka S, Toyooka KO, Reddy J, Miyajima K, Suzuki M, Shigematsu H, Takahashi T, Parikh G, Pass HI: Aberrant methylation of trail decoy receptor genes is frequent in multiple tumor types. Int J Cancer. 2004, 109: 786-792. 10.1002/ijc.20041
Kang JH, Kim SJ, Noh DY, Park IA, Choe KJ, Yoo OJ, Kang HS: Methylation in the p53 promoter is a supplementary route to breast carcinogenesis: correlation between CpG methylation in the p53 promoter and the mutation of the p53 gene in the progression from ductal carcinoma in situ to invasive ductal carcinoma. Lab Invest. 2001, 81: 573-579.
Cucer N, Taheri S, Ok E, Ozkul Y: Methylation status of CpG islands at sites -59 to +96 in exon 1 of the BRCA2 gene varies in mammary tissue among women with sporadic breast cancer. J Genet. 2008, 87: 155-158. 10.1007/s12041-008-0023-5
McShane LM, Altman DG, Sauerbrei W, Taube SE, Gion M, Clark GM: Reporting recommendations for tumor marker prognostic studies. J Clin Oncol. 2005, 23: 9067-9072. 10.1200/JCO.2004.01.0454
Suzuki MM, Bird A: DNA methylation landscapes: provocative insights from epigenomics. Nat Rev Genet. 2008, 9: 465-476. 10.1038/nrg2341
Ushijima T, Watanabe N, Shimizu K, Miyamoto K, Sugimura T, Kaneda A: Decreased fidelity in replicating CpG methylation patterns in cancer cells. Cancer Res. 2005, 65: 11-17.
Feinberg AP, Ohlsson R, Henikoff S: The epigenetic progenitor origin of human cancer. Nat Rev Genet. 2006, 7: 21-33. 10.1038/nrg1748
Laird PW: Principles and challenges of genome-wide DNA methylation analysis. Nat Rev Genet. 2010, 11: 191-203. 10.1038/nrg2732
Cotter TG: Apoptosis and cancer: the genesis of a research field. Nat Rev Cancer. 2009, 9: 501-507. 10.1038/nrc2663
Taylor RC, Cullen SP, Martin SJ: Apoptosis: controlled demolition at the cellular level. Nat Rev Mol Cell Biol. 2008, 9: 231-241. 10.1038/nrm2312
O'Driscoll M, Jeggo PA: The role of double-strand break repair - insights from human genetics. Nat Rev Genet. 2006, 7: 45-54. 10.1038/nrg1746
Okada H, Mak TW: Pathways of apoptotic and non-apoptotic death in tumour cells. Nat Rev Cancer. 2004, 4: 592-603. 10.1038/nrc1412
Li D, Da L, Tang H, Li T, Zhao M: CpG methylation plays a vital role in determining tissue- and cell-specific expression of the human cell-death-inducing DFF45-like effector A gene through the regulation of Sp1/Sp3 binding. Nucleic Acids Res. 2008, 36: 330-341. 10.1093/nar/gkm1028
Zelko IN, Mueller MR, Folz RJ: CpG methylation attenuates Sp1 and Sp3 binding to the human extracellular superoxide dismutase promoter and regulates its cell-specific expression. Free Radic Biol Med. 2010, 48: 895-904. 10.1016/j.freeradbiomed.2010.01.007
Veeck J, Esteller M: Breast cancer epigenetics: from DNA methylation to microRNAs. J Mammary Gland Biol Neoplasia. 2010, 15: 5-17. 10.1007/s10911-010-9165-1
Kelley RF, Totpal K, Lindstrom SH, Mathieu M, Billeci K, Deforge L, Pai R, Hymowitz SG, Ashkenazi A: Receptor-selective mutants of apoptosis-inducing ligand 2/tumor necrosis factor-related apoptosis-inducing ligand reveal a greater contribution of death receptor (DR) 5 than DR4 to apoptosis signaling. J Biol Chem. 2005, 280: 2205-2212. 10.1074/jbc.M410660200
Wu J: Apoptosis and angiogenesis: two promising tumor markers in breast cancer (review). Anticancer Res. 1996, 16: 2233-2239.
Fernandez Y, Gu B, Martinez A, Torregrosa A, Sierra A: Inhibition of apoptosis in human breast cancer cells: role in tumor progression to the metastatic state. Int J Cancer. 2002, 101: 317-326. 10.1002/ijc.10628
Evan GI, Vousden KH: Proliferation, cell cycle and apoptosis in cancer. Nature. 2001, 411: 342-348. 10.1038/35077213
Kalari S, Pfeifer GP: Identification of driver and passenger DNA methylation in cancer by epigenomic analysis. Adv Genet. 2010, 70: 277-308. full_text
Ramchandani S, Bhattacharya SK, Cervoni N, Szyf M: DNA methylation is a reversible biological signal. Proc Natl Acad Sci USA. 1999, 96: 6107-6112. 10.1073/pnas.96.11.6107
The authors thank the patient and control subjects for their participation. Ranjana Pal is grateful to the Council of Scientific and Industrial Research, New Delhi, India, for the pre-doctoral fellowship. Rameshwar N. K. Bamezai was supported by a research grant from the University Grants Commission to the National Centre of Applied Human Genetics, the University with Potential of Excellence project and the High-Throughput project from Department of Biotechnology (DBT), New Delhi.
The authors declare that they have no competing interests.
RP participated in the design of the experiments, experimental data acquisition, statistical analysis, interpretation of data and the writing of the draft. NS participated in the design of the experiments and interpretation of data. RC participated in experimental data acquisition. SG participated in the interpretation of data and gave critical suggestions. PG, NP and GA, as clinician and surgeon, participated in the acquisition of data. RB conceived the study, participated in the design of the experiments and the interpretation of data, revised the draft critically for intellectual content, and gave final approval of the version to be published. All authors have read and approved the final manuscript.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.