Cancer Research AACR Conference on Molecular Diagnostics - 2008  Translational Medicine Conference in Israel
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Cancer Research Clinical Cancer Research
Cancer Epidemiology Biomarkers & Prevention Molecular Cancer Therapeutics
Molecular Cancer Research Cancer Prevention Research
Cancer Prevention Journals Portal Cancer Reviews Online
Annual Meeting Education Book Meeting Abstracts Online

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Data
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Widschwendter, M.
Right arrow Articles by Laird, P. W.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Widschwendter, M.
Right arrow Articles by Laird, P. W.
[Cancer Research 64, 3807-3813, June 1, 2004]
© 2004 American Association for Cancer Research


Regular Articles

Association of Breast Cancer DNA Methylation Profiles with Hormone Receptor Status and Response to Tamoxifen

Martin Widschwendter5, Kimberly D. Siegmund2, Hannes M. Müller5, Heidi Fiegl5, Christian Marth5, Elisabeth Müller-Holzner5, Peter A. Jones1,3 and Peter W. Laird3,4

Departments of 1 Urology, 2 Preventive Medicine, 3 Biochemistry and Molecular Biology, and 4 Surgery, University of Southern California, Keck School of Medicine, Norris Comprehensive Cancer Center, Los Angeles, California, and 5 Department of Obstetrics and Gynecology, University Hospital, University of Innsbruck, Innsbruck, Austria


    ABSTRACT
 Top
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
We have generated DNA methylation profiles of 148 human breast tumors and found significant differences in hormone receptor (HR) status between clusters of DNA methylation profiles. Of 35 DNA methylation markers analyzed, the ESR1 gene, encoding estrogen receptor {alpha}, proved to be the best predictor of progesterone receptor status, whereas methylation of the PGR gene, encoding progesterone receptor, was the best predictor of estrogen receptor status. ESR1 methylation outperformed HR status as a predictor of clinical response in patients treated with the antiestrogen tamoxifen, whereas promoter methylation of the CYP1B1 gene, encoding a tamoxifen- and estradiol-metabolizing cytochrome P450, predicted response differentially in tamoxifen-treated and nontamoxifen-treated patients. High levels of promoter methylation of the ARHI gene, encoding a RAS-related small G-protein, were strongly predictive of good survival in patients who had not received tamoxifen therapy. Our results reveal an as yet unrecognized degree of interaction between DNA methylation and HR biology in breast cancer cells and suggest potentially clinically useful novel DNA methylation predictors of response to hormonal and non-hormonal breast cancer therapy.


    INTRODUCTION
 Top
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Breast cancer is the most common malignancy among females in most Western countries, who have an overall lifetime risk of >10% of developing invasive breast cancer (1) . The presence of estrogen receptor (ER) and/or progesterone receptor (PR) is an important diagnostic feature of breast cancer, reflective of disease etiology (2) and predictive of response to treatment with the antiestrogen tamoxifen (3 , 4) . Recent advances in molecular profiling by gene expression (cDNA) microarrays have led to a further refinement of the subclassification of breast cancer into five major subtypes and to the identification of gene expression signatures associated with prognosis (5, 6, 7, 8, 9) . Molecular profiling in breast cancer has thus far focused primarily on the use of cDNA microarrays, which are limited by the innate instability of RNA and are poorly compatible with the formalin fixation and paraffin embedding of tumor tissues used in routine histopathology. Therefore, we have explored the use of DNA methylation markers as an alternative approach to molecular profiling (10) . Hypermethylation of promoter CpG islands, which is frequently observed in breast cancer (11, 12, 13) , is often associated with transcriptional silencing of the associated gene, thus providing a DNA-based surrogate marker for expression status (14) . Microarray-based methods of DNA methylation analysis are hampered by modest quantitative accuracy, poor sensitivity to low levels of CpG island hypermethylation, and technical challenges in target DNA preparation, which requires either bisulfite PCR amplification of each individual locus (15) or the use of restriction enzyme digestion (11) , which is not consistently reliable with formalin-fixed tissues. As an alternative, we have used a moderate-throughput, fluorescence-based, semiautomated quantitative technique called MethyLight (16) to screen a panel of 35 methylation markers in 148 cases of breast cancer. Interestingly, we found that among these 35 markers, the best predictor of ER status was methylation of the PR gene (PGR). Conversely, the best predictor of PR status was methylation of the ER gene (ESR1).

Hormone receptor (HR) status, defined as ER and/or PR positivity, has been shown to predict response to tamoxifen treatment (3 , 4) . Interestingly, although tamoxifen is thought to act through the ER, PR status is an independent factor predictive of adjuvant endocrine treatment benefit (4) . Tamoxifen, which is a selective ER modulator, has been shown to dramatically reduce the risk of breast cancer (17) and of breast cancer recurrence (18) . Since its introduction more than 25 years ago, tamoxifen has been the mainstay of the endocrine adjuvant treatment of breast cancer, has become the most widely used anticancer drug, and may be considered one of the first targeted therapies (18) . In this study we found that ESR1 methylation predicts survival only in tamoxifen-treated patients and that ARHI methylation predicts survival only in non-tamoxifen-treated patients, whereas CYP1B1 methylation predicts survival differentially in tamoxifen-treated and nontreated patients. We propose that these differences in DNA methylation profiles reflect alternative pathways of tumorigenesis associated with differences in HR status, possibly due to different originating cell types (9) and/or disease etiology (2) .


    MATERIALS AND METHODS
 Top
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Tissues.
Tumor samples were retrieved from the tissue bank of the Department of Obstetrics and Gynecology, Innsbruck University Hospital (Innsbruck, Austria). Clinical, pathological, and follow-up data are stored in a database in accordance with hospital privacy rules. Specimens were brought to the pathologist (E. M-H.) immediately after resection, and part of the tissue was placed in liquid nitrogen and stored at –80°C until lyophilization. A total of 148 patients with breast cancer treated at the Department of Obstetrics and Gynecology, University of Innsbruck between 1989 and 2000 were included in this study. Patient characteristics are provided separately in Data Supplement 1.

Histopathological Analyses.
All breast cancer specimens were reviewed by a single pathologist (E. M-H.). HR positivity was defined as presence of ER and/or PR in >10% of tumor cells (immunohistochemistry was done for the 106 breast cancers) or >=15 fmol/mg protein (biochemical assays were performed for 42 breast cancer specimens, which were obtained before immunohistochemistry had been established in our laboratory).

DNA Methylation Analyses.
Genomic DNA was isolated using a QIAmp tissue kit (Qiagen, Hilden, Germany). Sodium bisulfite conversion of genomic DNA was performed as described previously (16) . DNA methylation analysis was performed by MethyLight (16 , 19) . Three sets of primers and probes designed specifically for bisulfite-converted DNA were used [a methylated set for the gene of interest and two reference sets, ß-actin (ACTB) and collagen 2A1 (COL2A1)] to normalize for input DNA. The specificity of the reactions for methylated DNA was confirmed separately using SssI (New England Biolabs)-treated human peripheral blood lymphocyte DNA (Promega), which results in near complete methylation of this reference DNA (16) . The percentage of fully methylated molecules at a specific locus was calculated by dividing the GENE:ACTB ratio of a sample by the GENE:ACTB ratio of SssI-treated sperm DNA and multiplying by 100 and calculated separately by dividing the GENE:COL2A1 ratio of a sample by the GENE:COL2A1 ratio of SssI-treated sperm DNA and multiplying by 100. The mean of these two resulting values was used in subsequent statistical analyses. We use the abbreviation PMR (percentage of fully methylated reference) to indicate this measurement (20) . The initial 64 methylation markers and the final panel of 35 markers were selected based on published reports demonstrating a role for DNA methylation in breast cancer or due to the fact that they are involved in HR action. Primer and probe sequences are shown in Data Supplement 2.

Expression Analyses.
RNA isolation and expression analyses were performed as described previously (19) . Before cDNA synthesis, RNA samples were treated with DNase to ensure removal of contaminating genomic DNA. TATA box-binding protein served as the reference gene. Primer and probe sequences are shown in Data Supplement 3.

Statistics.
To cluster samples and DNA methylation markers, we used agglomerative hierarchical cluster analysis in SPLUS 2000 (Insightful Corp.) Because many of the CpG regions had undetectable methylation, we categorized the PMR values into quartiles (coded 1–4). If >25% but <50% of the samples had undetectable methylation, this resulted in scores of 1 (undetectable methylation), 2 (detectable methylation and <=50th percentile), 3 (51st–75th percentile), and 4 (>75th percentile). If >50% but <75% of the samples had undetectable methylation, the scores were 1, 3, and 4. If >75% of the samples had undetectable methylation, the score was either 1 for undetectable methylation or 4 for detectable methylation. Manhattan distance, the sum of the absolute deviations across methylation markers, was used to measure dissimilarity. The dissimilarity between clusters was measured by the group average method (21) . We tested the association between (categorized) PMR values and HR status using logistic regression. Separate analyses were conducted for each gene. Multiple linear regression was used to study the relationship between DNA methylation and ESR1 gene expression. A total of 75 samples had ESR1 gene expression measured, but one was omitted from the analyses as an outlier (ESR1 upstream A expression > 3 SDs above the mean). Comparisons between groups of samples and between different genes were simplified by expressing all expression data relative to the mean for the entire set of 74 samples. We used Cox regression to study the association between PMR values and overall (and disease-free) survival, treating PMR quartiles as ordered categorical variables. Using an interaction model, we tested whether the association of PMR values and survival varied by treatment with tamoxifen therapy (received tamoxifen treatment versus did not receive tamoxifen treatment; tamoxifen treatment was defined as 20 mg of tamoxifen daily for 5 years or until recurrence of disease). The analyses were adjusted for nodal status (0, 1–3, and >3) and tumor stage (I, II, and III/IV). Nodal status was coded using indicator variables for two of the three levels. All analyses were age-adjusted.


    RESULTS
 Top
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
We prescreened 65 DNA methylation markers on a limited set of pilot samples (8 breast cancer cell lines and 8 breast carcinomas) to identify markers with sufficiently high methylation frequencies and/or methylation levels. From this initial set, 35 informative markers were selected for MethyLight analysis on 148 primary breast carcinomas obtained from the University Hospital, University of Innsbruck (a summary of clinical characteristics is shown in Data Supplement 1). A semiautomated MethyLight platform was used to execute these reactions, and data were obtained for 4978 of the 5180 analyses performed (96% success rate; a summary of the data is shown in Data Supplement 4).

Two-dimensional unsupervised hierarchical clustering analysis of cases versus methylation markers revealed that the tumors segregated naturally into groups of cases with distinct methylation profiles (Fig. 1)Citation . These two major clusters differed significantly in their HR status [P = 0.0011 for cluster 1 (indicated in green; n = 87) versus cluster 2 (indicated in red; n = 56) for ER+ versus ER–; P = 0.0013 for PR+ versus PR–; P = 0.0011 for HR+ versus HR–] and in age (P = 0.0080 for cluster 1 versus cluster 2; mean age, 57 versus 63 years, respectively). We adjusted for age in all subsequent analyses. We did not detect significant clustering of cases by HER2 status (22) , menopausal status, relapse, death, grade, nodes, stage, or tumor diameter (data not shown).



View larger version (86K):
[in this window]
[in a new window]
 
Fig. 1. Unsupervised hierarchical clustering of DNA methylation markers and breast carcinomas. Agglomerative hierarchical analysis was used to cluster samples (columns) and DNA methylation gene markers (rows) as listed on the right. Quartiles of DNA methylation measures (percentage of fully methylated reference values) are indicated by color, ranging from dark green indicating the lowest methylation (first quartile), and light green, light red, and dark red indicating the second, third, and fourth (highest methylation) quartiles, respectively. Failed analyses are indicated in white. The two major sample clusters are shown in green (cluster 1) and red (cluster 2). Hormone receptor (HR) status is indicated by black ovals (HR–) or white ovals (HR+).

 
The hierarchical clustering data suggest that HR status is associated with profiles of multiple methylation markers, but it does not reveal which markers contribute most to the clustering. We therefore investigated the association of each of the 35 markers with HR status individually, and we ranked them according to the strength of their association (Table 1)Citation . Fifteen of the 35 genes yielded Ps < 0.05 using methylation values (PMR measurements as described in "Materials and Methods") as predictors of HR (ER+ and/or PR+) status (Table 1Citation , HR Status Predictors). Since we tested 35 different markers simultaneously, we adjusted for the multiple comparisons by controlling the false discovery rate, which is the expected proportion of false positive tests among all positive tests (see "Materials and Methods"). After this adjustment, three markers (SOCS1, RASSF1A, and BCL2) were significantly associated with HR status (Table 1Citation , HR Status Predictors). Interestingly, SOCS1 deficiency results in accelerated mammary gland development in mice (23) and is known to be methylated in human tumors (24) . Promoter methylation of RASSF1A has previously been reported to be methylated in a large percentage of human breast cancers (25) , and this methylation can even be detected in epithelial hyperplasia and intraductal papillomas (26) . Here, we show that RASSF1A promoter methylation is associated with HR status in advanced breast tumors. Our observation of a significant negative association between BCL2 methylation and HR status is consistent with earlier reports that BCL2 expression is positively associated with HR positivity (27) and that its down-regulation is negatively associated with HR positivity (28) .


View this table:
[in this window]
[in a new window]
 
Table 1 DNA methylation markers as predictors of HRa status

The 35 DNA methylation markers are ranked according to the strength of their association with HR as determined by logistic regression. HR status (HR, ER, and PR, in which HR refers to ER and/or PR positive) was used as the outcome, and quartiled DNA methylation (PMR) values were used as predictors, with adjustment for age. Significant associations after adjustment for multiple comparisons are indicated in bold. In the ER analysis, the data were adjusted for PR status, whereas in the PR analysis, the data were adjusted for ER status. "Association" refers to the direction of the association (+ indicates a positive relationship between DNA methylation and HR status, and – indicates an inverse relationship).

 
Remarkably, of all 35 markers, DNA methylation of the ESR1 gene encoding the ER{alpha} was the least associated with HR status. PGR methylation was also not significantly associated with HR status after adjustment for multiple comparisons. Breast tumors are often concordant for ER and PR status. In our study, 127 of 148 breast cancer specimens were either double receptor (ER and PR) positive (n = 86) or double receptor negative (n = 41), whereas only 21 tumors were positive for either just ER (n = 12) or PR (n = 9). This highly significant association (P = 2.6 x 10–16 by {chi}2) between ER and PR status is attributed to induction of PGR gene expression by activated ER (29 , 30) . This makes it difficult to separate the effects of the two receptors. We addressed this problem by investigating which methylation markers best predict the status of ER and PR individually, while adjusting for the status of the other receptor and for age (Table 1Citation , ER Status Predictors and PR Status Predictors). Interestingly, methylation of the 5' CpG island of the PGR gene turned out to be the best predictor of ER status (Table 1Citation , ER Status Predictors). Positive ER status is inversely associated with PGR CpG island methylation, consistent with the well-established induction of PGR gene expression by activated ER (31) . On the other hand, ESR1 methylation turned out to be the best predictor of PR status (Table 1Citation , PR Status Predictors), even though it is the least significant predictor of overall HR status (Table 1Citation , HR Status Predictors). Thus, whereas neither ESR1 nor PGR methylation marker was a good predictor of overall HR status, each is the best predictor of the status of the other receptor, but not of their own cognate receptor. It is not a priori clear why this would be the case, or why there would be a positive association between PR status and ESR1 methylation (Table 1Citation , PR Status Predictors) because activated ER is known to induce PGR gene expression. This raises the question of whether ESR1 methylation is truly reflective of reduced ESR1 expression. We analyzed ESR1 expression by quantitative real-time reverse transcription-PCR in 74 samples for which we had frozen tissue available for RNA analysis. We did not find a clear inverse relationship between ESR1 expression levels and quartiled ESR1 methylation levels when the tumors were analyzed collectively (Fig. 2)Citation . However, PR– tumors did show a statistically significant inverse trend between ESR1 gene expression levels and ESR1 methylation levels (Fig. 2)Citation . This suggests that PR+ status may confer resistance of the ESR1 gene expression to ESR1 DNA methylation. We have further explored this interesting observation, as described in "Discussion," but we have not yet resolved the mechanism. It should be noted that the levels of methylation that we observed at the ESR1 locus are quite low, with a median PMR of 0.8 (see Data Supplement 4), and may therefore be more of a reflection of chromatin status at the ESR1 locus rather than a major driving force in silencing ESR1 gene expression. Regardless, it is yet another intriguing finding of an interaction between DNA methylation and HR biology.



View larger version (13K):
[in this window]
[in a new window]
 
Fig. 2. A, ESR1 gene expression and ESR1 DNA methylation. ESR1 gene expression from exon 1A was measured by real-time quantitative reverse transcription-PCR. The resulting values were divided by the expression values for the control gene TATA box-binding protein and normalized to a mean ratio of 1.0 across all 74 samples. Methylation levels were measured immediately downstream of the transcription start site of exon 1A (amplicon located from +14 to +114) by MethyLight, and the resulting percentage of fully methylated reference values was divided into quartiles. In the two panels on the right, the four quartiles were further subdivided into progesterone receptor-positive and progesterone receptor-negative tumors.

 
Hormone therapy with the antiestrogen tamoxifen is frequently used as an adjuvant treatment in breast cancer. HR status has been shown to predict response to tamoxifen treatment (3 , 4) . Our patients showed a similar trend using a Cox proportional hazards model to test whether the association between HR status and survival differed by treatment with tamoxifen therapy (mean follow-up, 5.2 years). However, the interaction with treatment response was not statistically significant (Table 2)Citation . Our results described above suggest a link between HR status and cellular DNA methylation profiles in breast cancer cells. Therefore, we tested whether any of the 35 DNA methylation markers would be better predictors of response to tamoxifen treatment than HR status. We ranked our 35 DNA methylation markers according to their ability to predict response to tamoxifen therapy as measured by a test for interaction in a Cox model with adjustments for treatment-specific effects of HR status. The model was also adjusted for age, stage, and number of positive nodes. Three markers were statistically significant predictors of response to tamoxifen therapy (Table 3)Citation . Interestingly, we had shown above that ESR1 methylation is a good predictor of PR status. Here, we show that ESR1 methylation outperforms PR status as a predictor of response to tamoxifen (Table 3)Citation . We further analyzed the relationship between ESR1 methylation and tamoxifen response by comparing the survival curves of patients with above or below median ESR1 methylation levels who either received tamoxifen therapy or did not receive tamoxifen therapy (Fig. 3)Citation . High ESR1 methylation was a significant predictor of better survival in the tamoxifen-treated group but showed no significant predictive value for the non-tamoxifen-treated group (Fig. 3)Citation .


View this table:
[in this window]
[in a new window]
 
Table 2 Cox regression analysis of the association between HRa status and overall survival or disease-free survival in patients treated with TAM or not treated with TAM

 

View this table:
[in this window]
[in a new window]
 
Table 3 Cox regression analysis of the association between DNA methylation PMRa values (grouped by quartiles) and overall survival or disease-free survival

Significant associations are indicated in bold.

 


View larger version (45K):
[in this window]
[in a new window]
 
Fig. 3. Survival curves. Kaplan-Meier survival curves are shown for the three genes that were statistically significant predictors of response to tamoxifen therapy in a Cox proportional hazards interaction model (Table 2)Citation . Above median methylation levels are indicated by HIGH mC for ESR1 and ARHI, whereas below median methylation levels are indicated by LOW mC for these two genes. For CYP1B1 methylation, POS mC refers to any detectable methylation, whereas NEG mC refers to a lack of detectable methylation. The high or positive methylation survival curves are shown in red, with squares indicating censored events. The low or negative methylation survival curves are shown in green, with circles indicating censored events.

 
ARHI promoter methylation was a highly significant predictor of survival in patients who had not received tamoxifen therapy but showed no predictive value for patients treated with tamoxifen therapy (Table 3Citation ; Fig. 3Citation ). Finally, CYP1B1 methylation was a highly significant predictor of tamoxifen response in the interaction model (Table 3)Citation . The survival curves reveal that this is due to a differential predictive behavior of this marker in tamoxifen-treated versus nontreated patients (Fig. 3)Citation . These results show that DNA methylation markers can outperform HR status as predictors of response to tamoxifen therapy and suggest that DNA methylation markers may be of clinical use in directing hormone therapy in breast cancer patients.


    DISCUSSION
 Top
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Significant progress has been made in recent years toward the implementation of DNA methylation markers as clinical tools in cancer detection and diagnosis (10) . Here, we explore their utility in classifying breast cancer tumors and in predicting response to hormonal therapy. Although our initial goal was open-ended and was not directed toward HR status or response to hormone therapy, unsupervised clustering of the data indicated a relationship between DNA methylation and HR status. To our surprise, methylation levels of the genes encoding the HRs were not the best predictors of HR status. This suggests that DNA methylation markers are not perfect inverse predictors of gene expression status and also that they may contain relevant information that is independent of gene expression status. Transcriptional repression by promoter DNA methylation is thought to be mediated through changes in chromatin structure. The association between DNA methylation and gene expression may therefore show threshold effects, rather than a simple linear relationship. Indeed, robust ESR1 expression in PR– tumors is seen only in the lowest quartile of ESR1 methylation (Fig. 2)Citation . The lack of observed effect of ESR1 methylation on ESR1 expression in PR+ tumors (Fig. 2)Citation is interesting. As mentioned earlier, the PMR values obtained for ESR1 methylation are very low, with a median PMR of 0.8, which may not be sufficiently high to cause gene silencing. Moreover, our ESR1 methylation assay is located immediately downstream of the transcription start site for exon 1A of the ESR1 gene (+14 to +114), rather than in the promoter region itself, due to limitations imposed by MethyLight primer design criteria. However, PR– tumors did show the expected inverse relationship between ESR1 methylation and expression with these same assays (Fig. 2)Citation . Others have shown that ESR1 gene expression can be activated in ER– cell lines by DNA methyltransferase and histone deacetylase inhibitors (32) . One hypothesis that could reconcile these disparate observations is that PR+ tumors rely more extensively on expression initiated at upstream exons of the ESR1 gene (33) . Expression driven from these upstream promoters would not be expected to be affected by the DNA methylation that we measure at exon 1A (34) . Indeed, we identified several putative progesterone response elements located near upstream exon 1C (data not shown). However, we found that the relative utilization of ESR1 exons 1A, 1B, 1C, 1D, and 1E was similar in PR+ and PR– tumors. A summary of this analysis is shown in Data Supplement 5. This interesting difference between PR+ and PR– tumors in the relationship between ESR1 methylation and expression will require further investigation.

Clinical and epidemiological studies in the past have suggested that breast cancer is composed of at least two distinct groups (2 , 35) . More recently, molecular profiling of breast cancer using gene expression profiles has revealed five distinct clusters composed of one basal-like subgroup, one ERBB2-overexpressing subgroup, two luminal-like subgroups, and one normal breast tissue-like subgroup (9) . Because we have not performed gene expression microarray experiments on our group of breast tumors, we cannot directly compare our DNA methylation clustering results with the five major groups identified by gene expression profiles. However, it seems likely that the DNA methylation cluster 2, which contains mostly HR+ tumors (Fig. 1)Citation , overlaps with the two luminal-like subgroups, which contain ER+ tumors (5) . The DNA methylation cluster 1 contains the majority of HR– tumors and likely overlaps with the other three gene expression subtypes, which tend to be ER– (5) . It seems likely that the gene expression profile subgroupings represent a much more stable subgrouping because these analyses are based on a much larger number of samples and genes (9) . Nevertheless, the undirected clustering of our methylation data led us to the identification of an interesting link between DNA methylation patterns and HR biology.

We chose to use MethyLight technology for this study, rather than methylation-specific PCR or the methylation microarray technologies currently under development. One of the unique features of MethyLight technology is that the resulting data are composed of a mixture of discrete and variable measures. The discrete measures arise from the large number of data points with undetectable methylation (PMR values of 0) versus the data points with positive detection of methylation. This type of data structure is similar to that obtained with methylation-specific PCR analysis. On the other hand, the quantitative nature of MethyLight also generates continuous measures for the samples with detectable levels of DNA methylation. We show here that useful information can be extracted from both types of measures. For example, among the methylation markers predictive of response to tamoxifen therapy, CYP1B1 was used as a discrete measure of positive versus negative DNA methylation, similar to methylation-specific PCR analysis. However, a methylation-specific PCR-based approach for the other two markers predictive of treatment response would have been noninformative because ESR1 and ARHI are positive in 100% and 99.3% of the samples, respectively (see Data Supplement 4). The quantitative aspect of MethyLight analysis was required to reveal the association of these methylation markers with response to tamoxifen therapy.

Of 35 DNA methylation markers tested, three genes showed the potential to serve as independent predictors of clinical response to systemic hormonal therapy with tamoxifen. Two of the three genes (ESR1 and CYP1B1) are known to be intimately involved in the function and metabolism of estradiol. This lends credence to the biological relevance of DNA methylation changes in breast tumors. The third gene (ARHI) encodes a RAS-related small G-protein, which may play a role in the regulation of breast cancer cell growth (36) . We found that patients with high levels of ARHI methylation had better survival than patients with low levels of ARHI methylation. However, this effect was completely obliterated in the tamoxifen-treated group (Fig. 3)Citation . This may be due to the ability of antiestrogens such as tamoxifen to block growth factor-induced mitogenesis, possibly involving pathways regulated by ARHI (36 , 37) .

ESR1 encodes the ER{alpha}. Patients treated with tamoxifen who had high levels of tumor ESR1 methylation showed better survival than tamoxifen-treated patients with low levels of ESR1 methylation. The survival benefit in patients with high levels of ESR1 methylation may be due, in part, to the positive association between ESR1 methylation and PR status (Table 1Citation , PR Status Predictors). PR status appears to be a better predictor of response to tamoxifen than ER status (Table 2)Citation .

CYP1B1 encodes cytochrome P450 1B1, which catalyzes the conversion of 17-ß-estradiol (E2) to the catechol estrogen metabolites 2-OH-E2 and 4-OH-E2. The 2-hydroxylated form of E2 has been shown to have weak ER agonist or antagonist properties (38) . CYP1B1 is also the principal catalyst of 4-hydroxytamoxifen trans-cis-isomerization, which converts the primary potent antiestrogen trans-4-hydroxytamoxifen to the weak estrogen agonist cis-4-hydroxytamoxifen (39) . We have not investigated gene expression levels of CYP1B1 as a function of CYP1B1 methylation, and our measured levels of methylation are quite low (see Data Supplement 4). Nevertheless, if patients with positive CYP1B1 methylation do indeed have reduced CYP1B1 expression, then these patients would be expected to have lower rates of 4-hydroxytamoxifen trans-cis-isomerization and would thus retain higher levels of active antiestrogen. This would be consistent with the better survival of these patients in the tamoxifen-treated group (Fig. 3)Citation . Conversely, in the patients who did not receive tamoxifen therapy, patients with tumor CYP1B1 methylation would have a reduced capacity for conversion of E2 to its weaker catechol derivatives. This would be consistent with the observation that these patients show a worse survival among the group not receiving tamoxifen therapy (Fig. 3)Citation .

Our results show a level of interaction between DNA methylation changes in breast cancer and HR status or response to hormonal therapy that was not previously appreciated. Because DNA methylation markers rely on DNA as an analyte, as opposed to the more chemically labile RNA molecule, these results suggest exciting opportunities for the development of robust assays for clinical diagnosis and for predicting response to antiestrogen therapy in the adjuvant setting.


    ACKNOWLEDGMENTS
 
We thank Drs. Mihaela Velicescu and Daniel J. Weisenberger for critical evaluation of the manuscript. We are grateful to Tiffany I. Long for technical advice and assistance.


    FOOTNOTES
 
Grant support: Austrian Science Foundation Grants J2024 and P15995-B05 (M. Widschwendter).

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Note: Supplementary data for this article are available at Cancer Research Online (http://cancerres.aacrjournals.org).

Requests for reprints: Peter W. Laird, University of Southern California/Norris Cancer Center, Room 6418, 1441 Eastlake Avenue, Los Angeles, CA 90089-9176. Phone: (323) 865-0650; Fax: (323) 865-0158; E-mail: plaird{at}usc.edu

Received 12/ 9/03. Revised 2/26/04. Accepted 3/23/04.


    REFERENCES
 Top
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 

  1. Feuer EJ, Wun LM, Boring CC, et al The lifetime risk of developing breast cancer. J Natl Cancer Inst (Bethesda), 85: 892-7, 1993.[Abstract/Free Full Text]
  2. Potter JD, Cerhan JR, Sellers TA, et al Progesterone and estrogen receptors and mammary neoplasia in the Iowa Women’s Health Study: how many kinds of breast cancer are there?. Cancer Epidemiol Biomark Prev, 4: 319-26, 1995.[Abstract]
  3. Anonymous. Tamoxifen for early breast cancer: an overview of the randomised trials. Early Breast Cancer Trialists’ Collaborative Group. Lancet, 351: 1451-67, 1998.[CrossRef][Medline]
  4. Bardou VJ, Arpino G, Elledge RM, Osborne CK, Clark GM. Progesterone receptor status significantly improves outcome prediction over estrogen receptor status alone for adjuvant endocrine therapy in two large breast cancer databases. J Clin Oncol, 21: 1973-9, 2003.[Abstract/Free Full Text]
  5. Perou CM, Sorlie T, Eisen MB, et al Molecular portraits of human breast tumours. Nature (Lond), 406: 747-52, 2000.[CrossRef][Medline]
  6. Sorlie 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-74, 2001.[Abstract/Free Full Text]
  7. van’t Veer LJ, Dai H, van de Vijver MJ, et al Gene expression profiling predicts clinical outcome of breast cancer. Nature (Lond), 415: 530-6, 2002.[CrossRef][Medline]
  8. van de Vijver MJ, He YD, van’t Veer LJ, et al A gene-expression signature as a predictor of survival in breast cancer. N Engl J Med, 347: 1999-2009, 2002.[Abstract/Free Full Text]
  9. Sorlie 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-23, 2003.[Abstract/Free Full Text]
  10. Laird PW. Early detection: the power and the promise of DNA methylation markers. Nat Rev Cancer, 3: 253-66, 2003.[CrossRef][Medline]
  11. Yan PS, Chen CM, Shi H, et al Dissecting complex epigenetic alterations in breast cancer using CpG island microarrays. Cancer Res, 61: 8375-80, 2001.[Abstract/Free Full Text]
  12. Yang X, Yan L, Davidson NE. DNA methylation in breast cancer. Endocr Relat Cancer, 8: 115-27, 2001.[Abstract]
  13. Widschwendter M, Jones PA. DNA methylation and breast carcinogenesis. Oncogene, 21: 5462-82, 2002.[CrossRef][Medline]
  14. Jones PA, Baylin SB. The fundamental role of epigenetic events in cancer. Nat Rev Genet, 3: 415-28, 2002.[Medline]
  15. Adorjan P, Distler J, Lipscher E, et al Tumour class prediction and discovery by microarray-based DNA methylation analysis. Nucleic Acids Res, 30: E21 2002.
  16. Eads CA, Danenberg KD, Kawakami K, et al MethyLight: a high-throughput assay to measure DNA methylation. Nucleic Acids Res, 28: E32 2000.
  17. Powles TJ. Anti-oestrogenic prevention of breast cancer: the make or break point. Nat Rev Cancer, 2: 787-94, 2002.[CrossRef][Medline]
  18. Jordan VC. Tamoxifen: a most unlikely pioneering medicine. Nat Rev Drug Discov, 2: 205-13, 2003.[CrossRef][Medline]
  19. Eads CA, Danenberg KD, Kawakami K, et al CpG island hypermethylation in human colorectal tumors is not associated with DNA methyltransferase overexpression. Cancer Res, 59: 2302-6, 1999.[Abstract/Free Full Text]
  20. Eads CA, Lord RV, Wickramasinghe K, et al Epigenetic patterns in the progression of esophageal adenocarcinoma. Cancer Res, 61: 3410-8, 2001.[Abstract/Free Full Text]
  21. Kaufman L, Rousseeuw PJ. . Finding groups in data: an introduction to cluster analysis, Wiley Interscience New York 1990.
  22. Slamon DJ, Clark GM, Wong SG, et al Human breast cancer: correlation of relapse and survival with amplification of the HER-2/neu oncogene. Science (Wash DC), 235: 177-82, 1987.[Abstract/Free Full Text]
  23. Lindeman GJ, Wittlin S, Lada H, et al SOCS1 deficiency results in accelerated mammary gland development and rescues lactation in prolactin receptor-deficient mice. Genes Dev, 15: 1631-6, 2001.[Abstract/Free Full Text]
  24. Yoshikawa H, Matsubara K, Qian GS, et al SOCS-1, a negative regulator of the JAK/STAT pathway, is silenced by methylation in human hepatocellular carcinoma and shows growth- suppression activity. Nat Genet, 28: 29-35, 2001.[CrossRef][Medline]
  25. Dammann R, Yang G, Pfeifer GP. Hypermethylation of the cpG island of Ras association domain family 1A (RASSF1A), a putative tumor suppressor gene from the 3p21.3 locus, occurs in a large percentage of human breast cancers. Cancer Res, 61: 3105-9, 2001.[Abstract/Free Full Text]
  26. Lehmann U, Langer F, Feist H, et al Quantitative assessment of promoter hypermethylation during breast cancer development. Am J Pathol, 160: 605-12, 2002.[Abstract/Free Full Text]
  27. Castiglione F, Sarotto I, Fontana V, et al Bcl2, p53 and clinical outcome in a series of 138 operable breast cancer patients. Anticancer Res, 19: 4555-63, 1999.[Medline]
  28. Park SH, Kim H, Song BJ. Down regulation of bcl2 expression in invasive ductal carcinomas is both estrogen- and progesterone-receptor dependent and associated with poor prognostic factors. Pathol Oncol Res, 8: 26-30, 2002.[Medline]
  29. Savouret JF, Bailly A, Misrahi M, et al Characterization of the hormone responsive element involved in the regulation of the progesterone receptor gene. EMBO J, 10: 1875-83, 1991.[Medline]
  30. Clark GM, Osborne CK, McGuire WL. Correlations between estrogen receptor, progesterone receptor, and patient characteristics in human breast cancer. J Clin Oncol, 2: 1102-9, 1984.[Abstract]
  31. Nardulli AM, Greene GL, O’Malley BW, Katzenellenbogen BS. Regulation of progesterone receptor messenger ribonucleic acid and protein levels in MCF-7 cells by estradiol: analysis of estrogen’s effect on progesterone receptor synthesis and degradation. Endocrinology, 122: 935-44, 1988.[Abstract]
  32. Yang X, Phillips DL, Ferguson AT, et al Synergistic activation of functional estrogen receptor (ER)-alpha by DNA methyltransferase and histone deacetylase inhibition in human ER-alpha-negative breast cancer cells. Cancer Res, 61: 7025-9, 2001.[Abstract/Free Full Text]
  33. Kos M, Reid G, Denger S, Gannon F. Minireview: genomic organization of the human ERalpha gene promoter region. Mol Endocrinol, 15: 2057-63, 2001.[Abstract/Free Full Text]
  34. Jones PA. The DNA methylation paradox. Trends Genet, 15: 34-7, 1999.[CrossRef][Medline]
  35. Fox MS. On the diagnosis and treatment of breast cancer. JAMA, 241: 489-94, 1979.[Abstract]
  36. Luo RZ, Fang X, Marquez R, et al ARHI is a Ras-related small G-protein with a novel N-terminal extension that inhibits growth of ovarian and breast cancers. Oncogene, 22: 2897-909, 2003.[CrossRef][Medline]
  37. Clarke R, Leonessa F, Welch JN, Skaar TC. Cellular and molecular pharmacology of antiestrogen action and resistance. Pharmacol Rev, 53: 25-71, 2001.[Abstract/Free Full Text]
  38. Gupta M, McDougal A, Safe S. Estrogenic and antiestrogenic activities of 16alpha- and 2-hydroxy metabolites of 17beta-estradiol in MCF-7 and T47D human breast cancer cells. J Steroid Biochem Mol Biol, 67: 413-9, 1998.[CrossRef][Medline]
  39. Crewe HK, Notley LM, Wunsch RM, Lennard MS, Gillam EM. Metabolism of tamoxifen by recombinant human cytochrome P450 enzymes: formation of the 4-hydroxy, 4'-hydroxy and N-desmethyl metabolites and isomerization of trans-4-hydroxytamoxifen. Drug Metab Dispos, 30: 869-74, 2002.[Abstract/Free Full Text]



This article has been cited by other articles:


Home page
Cancer Epidemiol. Biomarkers Prev.Home page
J. C. Baker Jr., J. H. Ostrander, S. Lem, G. Broadwater, G. R. Bean, N. C. D'Amato, V. K. Goldenberg, C. Rowell, C. Ibarra-Drendall, T. Grant, et al.
ESR1 Promoter Hypermethylation Does Not Predict Atypia in RPFNA nor Persistent Atypia after 12 Months Tamoxifen Chemoprevention
Cancer Epidemiol. Biomarkers Prev., August 1, 2008; 17(8): 1884 - 1890.
[Abstract] [Full Text] [PDF]


Home page
Cancer Res.Home page
E. M. Wolff, G. Liang, C. C. Cortez, Y. C. Tsai, J. E. Castelao, V. K. Cortessis, D. D. Tsao-Wei, S. Groshen, and P. A. Jones
RUNX3 Methylation Reveals that Bladder Tumors Are Older in Patients with a History of Smoking
Cancer Res., August 1, 2008; 68(15): 6208 - 6214.
[Abstract] [Full Text] [PDF]


Home page
J. Mol. Diagn.Home page
A. J. Iafrate and D. N. Louis
"MGMT for pt mgmt": Is Methylguanine-DNA Methyltransferase Testing Ready for Patient Management?
J. Mol. Diagn., July 1, 2008; 10(4): 308 - 310.
[Abstract] [Full Text] [PDF]


Home page
Toxicol SciHome page
H.-C. Yin, H.-P. Tseng, H.-Y. Chung, C.-Y. Ko, W.-S. Tzou, D. R. Buhler, and C.-H. Hu
Influence of TCDD on Zebrafish CYP1B1 Transcription during Development
Toxicol. Sci., May 1, 2008; 103(1): 158 - 168.
[Abstract] [Full Text] [PDF]


Home page
GutHome page
S. Ogino, T. Kawasaki, G. J Kirkner, Y. Suemoto, J. A Meyerhardt, and C. S Fuchs
Molecular correlates with MGMT promoter methylation and silencing support CpG island methylator phenotype-low (CIMP-low) in colorectal cancer
Gut, November 1, 2007; 56(11): 1564 - 1571.
[Abstract] [Full Text] [PDF]


Home page
Cancer Res.Home page
K. H. Taylor, R. S. Kramer, J. W. Davis, J. Guo, D. J. Duff, D. Xu, C. W. Caldwell, and H. Shi
Ultradeep Bisulfite Sequencing Analysis of DNA Methylation Patterns in Multiple Gene Promoters by 454 Sequencing
Cancer Res., September 15, 2007; 67(18): 8511 - 8518.
[Abstract] [Full Text] [PDF]


Home page
J. Mol. Diagn.Home page
S. Ogino, T. Kawasaki, G. J. Kirkner, P. Kraft, M. Loda, and C. S. Fuchs
Evaluation of Markers for CpG Island Methylator Phenotype (CIMP) in Colorectal Cancer by a Large Population-Based Sample
J. Mol. Diagn., July 1, 2007; 9(3): 305 - 314.
[Abstract] [Full Text] [PDF]


Home page
Mol Cancer ResHome page
M. F.G. de Maat, N. Umetani, E. Sunami, R. R. Turner, and D. S.B. Hoon
Assessment of Methylation Events during Colorectal Tumor Progression by Absolute Quantitative Analysis of Methylated Alleles
Mol. Cancer Res., May 1, 2007; 5(5): 461 - 471.
[Abstract] [Full Text] [PDF]


Home page
Cancer Res.Home page
V. Muthusamy, S. Duraisamy, C. M. Bradbury, C. Hobbs, D. P. Curley, B. Nelson, and M. Bosenberg
Epigenetic Silencing of Novel Tumor Suppressors in Malignant Melanoma
Cancer Res., December 1, 2006; 66(23): 11187 - 11193.
[Abstract] [Full Text] [PDF]


Home page
Clin. Cancer Res.Home page
K. Visvanathan, S. Sukumar, and N. E. Davidson
Epigenetic Biomarkers and Breast Cancer: Cause for Optimism.
Clin. Cancer Res., November 15, 2006; 12(22): 6591 - 6593.
[Full Text] [PDF]


Home page
J. Mol. Diagn.Home page
S. Ogino, T. Kawasaki, G. J. Kirkner, M. Loda, and C. S. Fuchs
CpG Island Methylator Phenotype-Low (CIMP-Low) in Colorectal Cancer: Possible Associations with Male Sex and KRAS Mutations
J. Mol. Diagn., November 1, 2006; 8(5): 582 - 588.
[Abstract] [Full Text] [PDF]


Home page
GutHome page
S Ogino, M Cantor, T Kawasaki, M Brahmandam, G J Kirkner, D J Weisenberger, M Campan, P W Laird, M Loda, and C S Fuchs
CpG island methylator phenotype (CIMP) of colorectal cancer is best characterised by quantitative DNA methylation analysis and prospective cohort studies
Gut, July 1, 2006; 55(7): 1000 - 1006.
[Abstract] [Full Text] [PDF]


Home page
J. Mol. Diagn.Home page
S. Ogino, T. Kawasaki, M. Brahmandam, M. Cantor, G. J. Kirkner, D. Spiegelman, G. M. Makrigiorgos, D. J. Weisenberger, P. W. Laird, M. Loda, et al.
Precision and Performance Characteristics of Bisulfite Conversion and Real-Time PCR (MethyLight) for Quantitative DNA Methylation Analysis
J. Mol. Diagn., May 1, 2006; 8(2): 209 - 217.
[Abstract] [Full Text] [PDF]


Home page
Mol Cancer ResHome page
T. M. Sissung, D. K. Price, A. Sparreboom, and W. D. Figg
Pharmacogenetics and Regulation of Human Cytochrome P450 1B1: Implications in Hormone-Mediated Tumor Metabolism and a Novel Target for Therapeutic Intervention
Mol. Cancer Res., March 1, 2006; 4(3): 135 - 150.
[Abstract] [Full Text] [PDF]


Home page
Genome ResHome page
M. Bibikova, Z. Lin, L. Zhou, E. Chudin, E. W. Garcia, B. Wu, D. Doucet, N. J. Thomas, Y. Wang, E. Vollmer, et al.
High-throughput DNA methylation profiling using universal bead arrays
Genome Res., March 1, 2006; 16(3): 383 - 393.
[Abstract] [Full Text] [PDF]


Home page
Cancer Res.Home page
H. Fiegl, S. Millinger, G. Goebel, E. Muller-Holzner, C. Marth, P. W. Laird, and M. Widschwendter
Breast Cancer DNA Methylation Profiles in Cancer Cells and Tumor Stroma: Association with HER-2/neu Status in Primary Breast Cancer
Cancer Res., January 1, 2006; 66(1): 29 - 33.
[Abstract] [Full Text] [PDF]


Home page
Nucleic Acids ResHome page
D. J. Weisenberger, M. Campan, T. I. Long, M. Kim, C. Woods, E. Fiala, M. Ehrlich, and P. W. Laird
Analysis of repetitive element DNA methylation by MethyLight
Nucleic Acids Res., December 2, 2005; 33(21): 6823 - 6836.
[Abstract] [Full Text] [PDF]


Home page
Mol. Endocrinol.Home page
D. Sharma, J. Blum, X. Yang, N. Beaulieu, A. R. Macleod, and N. E. Davidson
Release of Methyl CpG Binding Proteins and Histone Deacetylase 1 from the Estrogen Receptor {alpha} (ER) Promoter upon Reactivation in ER-Negative Human Breast Cancer Cells
Mol. Endocrinol., July 1, 2005; 19(7): 1740 - 1751.
[Abstract] [Full Text] [PDF]


Home page
Cancer Res.Home page
J. W.M. Martens, I. Nimmrich, T. Koenig, M. P. Look, N. Harbeck, F. Model, A. Kluth, J. Bolt-de Vries, A. M. Sieuwerts, H. Portengen, et al.
Association of DNA Methylation of Phosphoserine Aminotransferase with Response to Endocrine Therapy in Patients with Recurrent Breast Cancer
Cancer Res., May 15, 2005; 65(10): 4101 - 4117.
[Abstract] [Full Text] [PDF]


Home page
Cancer Epidemiol. Biomarkers Prev.Home page
K. Woodson, D. J. Weisenberger, M. Campan, P. W. Laird, J. Tangrea, L. L. Johnson, A. Schatzkin, and E. Lanza
Gene-Specific Methylation and Subsequent Risk of Colorectal Adenomas among Participants of the Polyp Prevention Trial
Cancer Epidemiol. Biomarkers Prev., May 1, 2005; 14(5): 1219 - 1223.
[Abstract] [Full Text] [PDF]