#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Breast Cancer DNA Methylation Profiles Are Associated with Tumor Size and Alcohol and Folate Intake


Although tumor size and lymph node involvement are the current cornerstones of breast cancer prognosis, they have not been extensively explored in relation to tumor methylation attributes in conjunction with other tumor and patient dietary and hormonal characteristics. Using primary breast tumors from 162 (AJCC stage I–IV) women from the Kaiser Division of Research Pathways Study and the Illumina GoldenGate methylation bead-array platform, we measured 1,413 autosomal CpG loci associated with 773 cancer-related genes and validated select CpG loci with Sequenom EpiTYPER. Tumor grade, size, estrogen and progesterone receptor status, and triple negative status were significantly (Q-values <0.05) associated with altered methylation of 209, 74, 183, 69, and 130 loci, respectively. Unsupervised clustering, using a recursively partitioned mixture model (RPMM), of all autosomal CpG loci revealed eight distinct methylation classes. Methylation class membership was significantly associated with patient race (P<0.02) and tumor size (P<0.001) in univariate tests. Using multinomial logistic regression to adjust for potential confounders, patient age and tumor size, as well as known disease risk factors of alcohol intake and total dietary folate, were all significantly (P<0.0001) associated with methylation class membership. Breast cancer prognostic characteristics and risk-related exposures appear to be associated with gene-specific tumor methylation, as well as overall methylation patterns.


Published in the journal: . PLoS Genet 6(7): e32767. doi:10.1371/journal.pgen.1001043
Category: Research Article
doi: https://doi.org/10.1371/journal.pgen.1001043

Summary

Although tumor size and lymph node involvement are the current cornerstones of breast cancer prognosis, they have not been extensively explored in relation to tumor methylation attributes in conjunction with other tumor and patient dietary and hormonal characteristics. Using primary breast tumors from 162 (AJCC stage I–IV) women from the Kaiser Division of Research Pathways Study and the Illumina GoldenGate methylation bead-array platform, we measured 1,413 autosomal CpG loci associated with 773 cancer-related genes and validated select CpG loci with Sequenom EpiTYPER. Tumor grade, size, estrogen and progesterone receptor status, and triple negative status were significantly (Q-values <0.05) associated with altered methylation of 209, 74, 183, 69, and 130 loci, respectively. Unsupervised clustering, using a recursively partitioned mixture model (RPMM), of all autosomal CpG loci revealed eight distinct methylation classes. Methylation class membership was significantly associated with patient race (P<0.02) and tumor size (P<0.001) in univariate tests. Using multinomial logistic regression to adjust for potential confounders, patient age and tumor size, as well as known disease risk factors of alcohol intake and total dietary folate, were all significantly (P<0.0001) associated with methylation class membership. Breast cancer prognostic characteristics and risk-related exposures appear to be associated with gene-specific tumor methylation, as well as overall methylation patterns.

Introduction

Breast cancer is the most common non-skin cancer among American women. The American Cancer Society's estimates indicate approximately 1.3 million new cases of invasive breast cancer were diagnosed globally in 2007; and nearly 500,000 women died from the disease [1]. Currently, there are over 2.5 million breast cancer survivors in the US, and an estimated $8.1 billion dollars is spent each year on treatment of breast cancer [2].

The principal prognostic indicator currently in clinical use for breast cancer is the tumor-node-metastasis (TNM) stage [3], [4]. Morphological attributes of malignant tumors that influence disease prognosis are the size of the primary tumor (T), presence and extent of regional lymph node involvement (N) and presence of distant metastases (M). Molecular attributes of tumors are also considered in clinical decision-making; loss of hormone receptor expression [5] and increased expression of ERBB2 [6] have each been associated with poor prognosis. Although numerous recent studies have demonstrated that alterations of DNA methylation in breast cancers are common and may be important etiologic and prognostic markers [7]-[14], large gaps in our knowledge remain. There is a notable lack of studies examining tumor DNA methylation in relation to breast cancer risk factors such as diet or reproductive factors in conjunction with other important tumor markers. Patient exposures such as alcohol and folate intake have potentially strong mechanistic links to epigenetic dysregulation [15]. In addition, recent work in-vitro and in animal models suggest that long term exposure to estrogen may lead to epigenetic effects and altered profiles of DNA methylation [16], [17]. To explore associations of tumor methylation with important tumor and patient characteristics, we analyzed tumors from breast cancer patients in the Kaiser Permanente Division of Research Pathways Study using a large scale methylation array.

Results

Unsupervised clustering and locus-by-locus analysis

Table 1 shows the patient demographic, hormonal, dietary and tumor characteristics for the 162 women overall (and stratified by menopausal status in Table S1). Results of unsupervised hierarchical clustering of the 750 most variable CpG loci indicate the epigenetic heterogeneity of these tumors (Figure 1).

Unsupervised clustering heatmap of CpG methylation in breast carcinomas.
Fig. 1. Unsupervised clustering heatmap of CpG methylation in breast carcinomas.
Unsupervised hierarchical clustering heat map based on Manhattan distance and average linkage of the 750 autosomal CpG loci with the highest variance. Samples are in rows (n = 162), and CpG loci are in columns. Blue indicates methylated and yellow indicates unmethylated.

Tab. 1. Patient demographic, hormonal, dietary, and tumor characteristics.
Patient demographic, hormonal, dietary, and tumor characteristics.

In array-wide locus-by-locus analysis the strongest associations of methylation of individual loci (Q-values <0.05) were observed for tumor grade (loci n = 209), tumor size (loci n = 74), estrogen receptor status (loci n = 183), progesterone receptor status (loci n = 69), and triple negative status (tumors negative for both estrogen and progesterone receptors as well as ERBB2; loci n = 130; Table S2). Together with tumor size, patient lymph node status is used in tumor staging. Among five CpG loci whose methylation was significantly associated (Q<0.05) with lymph node status, four (two in COL1A2, and one each in LOX and P2RX7) were also associated with tumor size (Q<0.05). Additionally, there was a trend of increased methylation associated with increased tumor size: for all 74 CpG loci that were significantly associated with tumor size (Q<0.05) methylation increased with larger tumor size. Similarly, all five CpGs associated with disease-positive lymph nodes had increased methylation in tumors in women with disease-positive lymph nodes. Details of locus-by-locus analyses for tumor grade, size, hormone receptor, and triple negative status (loci with Q<0.05) are given in Table S3.

Array validation

Methylation array validation was performed at CpGs with highly ranked associations from locus-by-locus analysis. The array CpG whose methylation was most significantly increased with increasing tumor stage was in the FES gene (Table S3) and array methylation was significantly correlated with Sequenom methylation (rho = 0.68, P = 1.1E-12, n = 85; Figure 2A). Promoter CpGs in P2RX7 and HSD17B12 had significantly increased methylation (Q<0.0001, and Q = 0.01 respectively) with increasing tumor size (Table S3) and array methylation at these CpGs were significantly correlated with Sequenom methylation (P2RX7; rho = 0.65, P = 8.6E-12, n = 88; HSD17B12; rho = 0.34, P = 5.4E-05, n = 137; Figure 2B and 2C). A promoter CpG in GSTM2 had significantly increased methylation with increasing tumor grade (Table S3) and array methylation was significantly correlated with Sequenom methylation (rho = 0.83, P<2.2E-16, n = 140; Figure 2D). Additionally, in all cases, Sequenom methylation values were significantly associated with respective covariates; tumor stage with FES methylation (P = 0.05), tumor size with P2RX7 (P<0.005) and HSD17B12 methylation (P<0.02), and tumor grade with GSTM2 methylation (P<0.001). Furthermore, relative mRNA expression of GSTM2 was significantly decreased among tumors with high array methylation at both CpGs associated with tumor grade (P<0.001 and P<0.03, Figure S1).

Array methylation is validated by Sequenom EpiTYPER.
Fig. 2. Array methylation is validated by Sequenom EpiTYPER.
Results from GoldenGate array methylation values are plotted versus Sequenom EpiTYPER quantitative methylation values. (A) Sequenom FES methylation is significantly correlated with GoldenGate methylation average β at the coordinate array CpG (Spearman correlation rho = 0.68, P = 1.1E-12, n = 85). (B) Sequenom P2RX7 methylation is significantly correlated with GoldenGate methylation average β at the coordinate array CpG (rho = 0.65, P = 8.6E-12, n = 88). (C) Sequenom HSD17B12 methylation is significantly correlated with GoldenGate methylation average β at the coordinate array CpG (rho = 0.34, P = 5.4E-05, n = 137). (D) Sequenom GSTM2 methylation is significantly correlated with GoldenGate methylation average β at the coordinate array CpG (rho = 0.83, P<2.2E-16, n = 140).

Clustering of DNA methylation patterns with RPMM

In order to explore overall methylation profiles of these tumors and their potential relationships with patient demographic, tumor and exposure characteristics we applied a modified model-based form of unsupervised clustering known as recursively partitioned mixture modeling (RPMM) [18]. The RPMM resulted in the eight methylation classes (average methylation profiles shown in Figure 3). Patient race was significantly associated with methylation class membership (P = 0.015, Table 2), with the majority of African Americans (54%) residing in class 2, and 40% of Hispanic cases residing in class 4. An association between methylation class membership and alcohol consumption approached statistical significance (P = 0.07, ever vs. never drinker, Table 2). Both supplemental folic acid intake (µg/day) and total dietary folate (µg/day) had associations with methylation class membership that approached statistical significance (P = 0.06 and P = 0.08 respectively; Table 2). For both folate variables, cases in methylation class 4 had the lowest intake and cases in methylation class 6 had the highest intake. Of the tumor characteristic variables, only tumor size was significantly associated with overall methylation profile (P = 0.0006, Table 2).

Recursively partitioned mixture model of CpG methylation in breast carcinomas.
Fig. 3. Recursively partitioned mixture model of CpG methylation in breast carcinomas.
The figure depicts the results of RPMM. Columns represent CpG sites and rows represent methylation classes. The height of each row is proportional to the number of observations residing in the class, total n = 162. Blue indicates methylated and yellow indicates unmethylated. Methylation classes are numbered one through eight on the left. The color of the columns within each class represents the average methylation of the CpG for that class.

Tab. 2. RPMM methylation calss membership by patient demographic and tumor characteristic covariates.
RPMM methylation calss membership by patient demographic and tumor characteristic covariates.
*Running 10,000 permutations.

Trends of DNA methylation related to alcohol and folate intake

Associations between alcohol intake and dietary folate and methylation class membership approached statistical significance. While methylation of only one CpG locus (in IL17RB) was significantly associated with folate intake in locus-by-locus tests (Q<0.05), regression coefficients from univariate locus-by-locus analysis plotted against their respective P-values revealed trends in the pattern of methylation for both alcohol and folate intake. Figure 4A illustrates the strong trend for patients with increasing alcohol intake to have negative regression coefficients, indicative of decreased methylation. In contrast, the trend for patients with increasing total dietary folate shows a strong shift to positive regression coefficients, indicative of increased methylation (Figure 4B).

There is an opposite trend for direction of association between breast carcinoma CpG methylation and alcohol intake compared to folate intake.
Fig. 4. There is an opposite trend for direction of association between breast carcinoma CpG methylation and alcohol intake compared to folate intake.
P-values for alcohol intake (g/day) and total dietary folate (µg/day) are plotted versus regression coefficients from locus-by-locus analysis of CpG methylation. Horizontal blue dotted line intercepts the y-axis at 0.05 to illustrate significance (before correction for multiple comparisons). The vertical solid black like intercepts the x-axis at zero to illustrate the contrasting trends. (A) There is a trend toward decreased methylation with increasing alcohol intake. (B) There is a trend toward increased methylation with increasing dietary folate.

Multivariate modeling of RPMM classes

The relationships between methylation classes and several covariates of interest were then modeled together using multinomial logistic regression in order to adjust for other factors in the model. Patient age, alcohol consumption, total dietary folate, and tumor size were each strongly associated with methylation class membership when controlling for all modeled variables (all Wald P-values <0.0001) and complete model details are given in Table S4. Figure 5 displays an illustration of the model results for covariates significantly associated with methylation classes. As alcohol consumption increased, there was an increased probability of cases residing in methylation classes 3 and 8, and a concomitant decrease in the probability of cases residing in classes 2 and 4 (Figure 5B). Increasing total dietary folate intake imparted a striking increase in the probability of membership in class 6, and a decreased probability of class membership in classes 1, 3, 4, and 7 (Figure 5C). The strong association between tumor size and methylation class membership remained after controlling for potential confounders, with the probability of patients being in class 2 increasing from about 20% to about 60% across the span of tumor size from 0 mm to 80+mm (Figure 5D). Accompanying this trend for tumor size were simultaneous decreases in the probability of cases with increasingly large tumors residing in classes 1 and 5–8, while tumor size had less influence on the probability for residing in classes 3 or 4 (Figure 5D).

Probability of methylation class membership is significantly associated with tumor size, patient age, alcohol intake, and dietary folate when controlling for potential confounders in a multinomial logistic regression model.
Fig. 5. Probability of methylation class membership is significantly associated with tumor size, patient age, alcohol intake, and dietary folate when controlling for potential confounders in a multinomial logistic regression model.
Results from a multinomial logistic regression plot the probability of methylation class membership versus covariates controlled for age, race, alcohol consumption, total dietary folate, tumor stage (low vs. high), tumor grade, tumor size, estrogen receptor status, and histology. The referent class (methylation class 3) is on the bottom of the plot in blue-green, remaining classes are plotted in numeric order from bottom to top as shown in the legend. (A) Patient age is significantly associated with methylation class membership (Overall Wald P<0.0001), and all methylation classes except class 4 are individually significantly associated with patient age. (B) Alcohol consumption is significantly associated with methylation class membership (Wald P<0.0001), and methylation classes 2, 4, 5, and 8 are individually significantly associated with alcohol intake. (C) √Total dietary folate intake is significantly associated with methylation class membership (Wald P<0.0001), and all methylation classes are individually significantly associated with total dietary folate. (D) Tumor size is significantly associated with methylation class membership (Wald P<0.0001), and all methylation classes except class 4 are individually significantly associated with tumor size.

Hormone receptor status among postmenopausal cases

Although neither estrogen nor progesterone receptor status were significantly associated with RPMM methylation profiles, large numbers of specific CpG loci had significant methylation associations with these tumor characteristics in locus-by-locus analysis (Table S2 and Table S3). Compared to the overall population of women diagnosed with breast cancer in the Kiaser Permanente Northern California cancer registry from 200–2009, this surgical cohort has a higher prevalence of hormone receptor positivity (78% overall vs. 88% here), particularly among pre-menopausal women's tumors (74% overall vs. 95% here). We therefore stratified on menopausal status, running RPMM on methylation data from post menopausal patients' tumors only (n = 117). This model resulted in eleven methylation classes (Figure S2) and methylation class membership was significantly associated with estrogen receptor status (P<0.03), and the association for triple negative tumors approached significance (P = 0.07) detailed results available in Table S5.

Discussion

It is becoming increasingly common to include data on molecular alterations from patient tumor samples into routine clinical practice as a means of improving prognosis and evaluating the predictive power of alterations of interest. As technology improves and population-based studies and clinical trials are conducted, medicine is being ushered into a new era of molecular characterization of disease. Tumor-node-metastasis (TNM) stage is the current prognostic indicator for breast cancer, though several clinical trials are currently under way to investigate the utility of molecular markers [19], and as more patients elect neoadjuvant therapy (specifically pre-operative chemotherapy), improved clinical staging and additional staging tools are poised to have great impact. Most current studies and one commercially available tool (Oncotype DX) are focused on gene expression markers, though the inherent instability of mRNA may make implementation of these strategies challenging outside of major surgical centers or centralized commercial laboratories. In contrast, DNA methylation is a stable mechanism of control of transcription, and the stability of DNA makes it an attractive target for accurate and reproducible assessment. Here we reported that tumor size, a cornerstone of breast cancer prognosis, is associated with tumor DNA methylation profile. In addition, we found that alcohol and folate intake, exposures related to disease risk, are independently associated with tumor DNA methylation profiles. This work sheds light on the relationship between important etiologic exposures and molecular subclasses of disease, extends the evidence for the utility of molecular characterization in tumor staging, and can be accomplished with minimal tissue in a pre-operative context.

The recently updated American Joint Committee on Cancer (AJCC) staging manual for breast cancer does not include additional molecular markers, though the committee acknowledged their consideration of markers such as hormone receptor status and stated that TNM staging “may play increasingly less important roles than understanding the biology of the cancer” [4]. Examining TNM variables we found that overall DNA methylation profile and methylation alterations in dozens of individual CpG loci were significantly associated with tumor size (all increased methylation). In contrast, methylation alterations of only five CpG loci (two in COL1A2, and one each in FAS, LOX, and P2RX7) were significantly associated with disease-positive lymph nodes. However, methylation of four of five lymph-node-positive associated CpGs (excepting FAS) were also significantly associated with tumor size, suggesting that these phenotypes are mechanistically related, and at least in part manifest via epigenetic alterations. As FAS encodes a TNF-receptor involved in regulating apoptosis it is not surprising that methylation-induced silencing of this receptor is associated with disease-positive lymph node status. In addition, hypermethylation of COL1A2 (collagen type I, alpha 2) has been associated with both proliferation and migration activity in bladder cancer [20], LOX is involved in the control of normal collagen deposition [21], and P2RX7 loss has been linked to morphologic changes in stroma related to altered collagen fibril alignment [22]. Collectively these data suggest that perturbations in collagen and collagen-related genes promote tumor growth and invasion, perhaps by altering the architecture of connective tissues in the tumor microenvironment. In support of this hypothesis, recent work in a mouse model has shown that altered mammary stromal tissue collagen expression significantly increases tumor formation and invasiveness potential [23]. Additionally, Chernov et al. showned that epigenetic alterations in collagen and collagen-related genes allows the deposition of an invasion-promoting collagen matrix in both breast and brain tumor cell lines [24].

The primary objective of TNM staging is to provide a standard prognosis nomenclature for patient care [4], and our results suggest that methylation markers may be a robust proxy for tumor size. Importantly, broader application of neoadjuvant therapy complicates breast cancer staging since chemotherapy can considerably decrease tumor size prior to surgical treatment, and it is still unclear whether clinical or pathologic stage best informs prognosis and treatment decisions [19]. The AJCC has added methodology (yc or ypTNM) for differentiating clinical and pathologic staging; in part, this is from recognition of the increasing use of neoadjuvant therapy for patients with operable, early stage disease [4], [25], [26]. Our data illustrate the promise of tumor DNA methylation for augmenting tumor staging. However, additional study of the relationship between tumor methylation and size in both pretreatment and postoperative samples is necessary. Specifically, the value of methylation to act as an additional marker of size in the neoadjuvant setting should be evaluated in future studies that compare both imaging and pathologically based size determination.

In order to evaluate the predictive power of DNA methylation profiles and individual loci for disease prognosis and recurrence, these patients continue to be followed for these events. Associations between DNA methylation and patient survival have been reported for individual genes such as GSTP1 and PITX [7], [8], [10], though overall DNA methylation profiles, or patterns of methylation at selected CpG loci or genes, may improve predictive power. Well recognized molecular subtypes of breast cancer such as hormone receptor negative and ERBB2 over-expressing tumors are known to be associated with reduced survival [27], and it will be necessary to extensively examine methylation markers stratified by commonly used molecular tumor markers. However, we did not find significant associations between ERBB2 status and CpG methylation in our analysis. Nonetheless, other well recognized molecular subtype markers; estrogen receptor, progesterone receptor, and triple negative status were among the covariates with the highest number of significant CpGs from array-wide locus-by-locus analysis. However, hormone receptor status and triple negativity were not associated with methylation profile when modeling all cases. Premenopausal patients' tumors in our surgical cohort had a higher prevalence of hormone receptor positivity compared to the overall population of premenopausal patients diagnosed with breast cancer. In order to address the potential bias this introduced we modeled the methylation profiles of postmenopausal patients' tumors separately and found a significant association between estrogen receptor status and methylation class. Additional study will be needed to better understand the role of hormone receptor and growth factor receptor expression in these tumors as they relate to methylation profile in the context of a patient's menopausal status.

We found significant, independent associations between both alcohol and folate intake and overall tumor DNA methylation profiles when controlling for potential confounders. Folate is a B vitamin that donates its methyl group for homocysteine remethylation to methionine as part of one-carbon metabolism. In turn, methionine is the methyl donor for DNA methylation via S-adenosyl methionine. However, alcohol is known to interfere with folate absorption in the intestine and hepatic release of folate, and hence, supply to tissues [28]. In fact, strong evidence of an etiologic role for alcohol in breast cancer has been reported in multiple meta-analyses of prospective and case-control studies with an excess risk for each alcoholic drink per day of about 10% [29], [30]. In contrast, meta-analysis of prospective studies has not provided clear support for an overall protective association between folate intake and breast cancer risk [31]. Yet, meta-analysis of case control studies of dietary folate, including results from the Shanghai Breast Cancer Study (whose participants are not regular alcohol drinkers) generally support a protective role for folate [31], [32].

While there have been numerous studies of alcohol and folate in relation to risk of breast cancer, investigations of the relationship between these exposures and epigenetic alterations in tumors themselves are scarce. Tao et al. reported that the prevalence of breast tumor methylation at CDKN2A, CDH1, and RARB did not differ by folate intake or lifetime alcohol consumption in genotype strata of one-carbon metabolism enzymes methylenetetrahydrofolate reductase (MTHFR) and methionine synthase (MTR) [33]. Consistent with these findings (and perhaps the lack of similar null results in the literature), we too did not find associations between alcohol or folate and methylation of CpG loci in CDKN2A, CDH1, and RARB. Further, after correcting for multiple comparisons, no CpG loci had significant alcohol-related methylation, and only one CpG locus (in the IL17RB promoter) was associated with folate intake. Alone, these results suggested that folate and alcohol intake do not influence tumor DNA methylation. However, plots of regression coefficients indicated strong independent trends for increased folate and reduced alcohol intake associations with increased CpG methylation. Since global, low-level effects of alcohol and folate intake on CpG methylation may not be detectable at individual CpGs in a genome-wide context, we examined the global relationships between alcohol or folate intake and DNA methylation using RPMM methylation classes. Modeling both exposures together revealed highly significant, independent associations between alcohol and folate and DNA methylation profile. Another human cancer for which alcohol is an important etiologic factor is head and neck squamous cell carcinoma, and previous work from our group demonstrated a similar relationship between DNA methylation profiles of these tumors and alcohol consumption [34]. Taken together with the weak mutagenic potential of alcohol [35], these results suggest that a major carcinogenic mechanism of action of alcohol is interference with epigenetic regulation through disruption of one-carbon metabolism.

In summary, we found tumor DNA methylation associated with tumor characteristics predictive of prognosis, and DNA methylation and patient exposures known to be related to disease risk. Additional study is needed to determine the prognostic value of DNA methylation markers. However, the potential clinical utility of tumor-size-related DNA methylation is apparent.

Materials and Methods

Study population

The Pathways Study is a prospective cohort study of breast cancer survival actively recruiting women diagnosed with invasive breast cancer from the Kaiser Permanente Northern California (KPNC) patient population since January 2006. Further study details are provided elsewhere [36]. Written informed consent is obtained from all participants before they are enrolled in the study. The study was approved by the IRB of KPNC and all collaborating sites.

Demographic, hormonal, and dietary factors

During the in-person baseline interview, participants were asked detailed information on family history of cancer and reproductive history, including: age at first full-term pregnancy, number of biological children, breastfeeding, and menopausal status. Additional information was collected on smoking, alcohol use, hormone use (oral contraceptives, hormone replacement therapy), and demographics (age at breast cancer diagnosis, race/ethnicity, household income, education). Self-reported height and weight around diagnosis was obtained to calculate body mass index (BMI, kg/m2). Any missing values were supplemented by concurrent information from KPNC electronic medical records.

Tumor characteristics

Data on estrogen and progesterone receptor status and ERBB2 expression were obtained from the KPNC Cancer Registry [37]. Tumor size was measured in a uniform manner by participating study pathologists. Data are collected, coded, and added to the KPNC registry approximately four months post-diagnosis to allow for the completion of treatment. For all breast surgical specimens, hormone receptor status and ERBB2 expression is routinely determined by IHC at the KPNC regional IHC lab, and if the IHC staining for ERBB2 expression is equivocal (less than 30% strong staining, but more than 10% weak staining), by fluorescence in situ hybridization at the KPNC regional cytogenetics lab.

Study samples

162 tumor specimens from the initial diagnostic biopsy were obtained from the KPNC tumor biorepository for methylation analysis. All tumor specimens were from patients who did not receive neoadjuvant chemotherapy.

Methylation analysis

FFPE tissue DNA was extracted using the QIAamp DNA mini kit according to the manufacturer's protocol (Qiagen, Valencia, CA). DNA was treated with sodium bisulfite to convert unmethylated cytosines to uracil using the EZ DNA Methylation Kit (Zymo Research, Orange, CA) according to the manufacturer's protocol. Illumina GoldenGate methylation bead arrays were used to simultaneously interrogate 1505 CpG loci associated with 803 cancer-related genes. Bead arrays have a similar sensitivity as quantitative methylation-specific PCR and were run at the UCSF Institute for Human Genetics, Genomics Core Facility according to the manufacturer's protocol and as described by Bibikova et al [38]. GoldenGate array methylation data are publicly available on the Gene Expression Omnibus archive, accession GSE22290.

Array methylation validation by Sequenom EpiTYPER mass spectroscopy

Array methylation was validated with Sequenom EpiTYPER base-specific cleavage and MALDI-TOF MS of bisulfite treated DNA [39]. EpiTYPER assays were designed for CpG loci both with significant associations between methylation and tumor characteristic variables as well as a high standard deviation of methylation values across samples. One assay (for COL1A2) failed the design process. Samples were processed at the UCSF Institute for Human Genetics, Genomics Core Facility. Briefly, PCR with primers located on either side of the CpG sites of interest are transcribed into an RNA transcript and cleaved base specifically. The cleavage products are analyzed by MALDI-TOF MS, and a characteristic mass signal pattern that distinguishes methyl-cytosine from thymine is obtained.

Gene expression by RT–PCR

Messenger RNA expression was measured using RT-PCR with preamplification using a validated approach [40]. RNA extraction was performed using the RecoverAll (Ambion), with a 16 hour tissue digestion and yields were determined using a Nanodrop spectrophotometer. Samples were concentration-normalized and reverse-transcribed with iScript cDNA synthesis kit (BioRad). Following cDNA synthesis, we performed linear, gene specific preamplification of samples and controls using the TaqMan preamp protocol (Applied Biosystems). Relative expression was measured using a HT7900 real time PCR instrument (Applied Biosystems).

Statistical analysis

Data assembly

Data were assembled with BeadStudio methylation software from Illumina (SanDiego, CA). All array data points are represented by fluorescent signals from both methylated (Cy5) and unmethylated (Cy3) alleles, and methylation level is given by β = (max(Cy5, 0))/(|Cy3|+|Cy5|+100), the average methylation (β) value is derived from the ∼30 replicate methylation measurements. Raw average β values were analyzed without normalization as recommended by Illumina. At each locus for each sample the detection P-value was used to determine sample performance; all samples, had detection P-values <1e-5 at more than 75% of CpG loci and passed performance criteria. CpG loci with a median detection P-value >0.05 (n = 8, 0.5%), were eliminated from analysis. All CpG loci on the X chromosome were excluded from analysis. The final dataset contained 1413 CpG loci associated with 773 genes.

Unsupervised clustering

Subsequent analyses were carried out using the R software [41]. For exploratory and visualization purposes, hierarchical clustering was performed using R function hclust with Manhattan metric and average linkage. To discern and describe the relationships between CpG methylation and patient and tumor covariates a modified model-based form of unsupervised clustering known as recursively partitioned mixture modeling (RPMM) was used as described in [18] and as used in [42]. Permutation tests (running 10,000 permutations) were used to test for association with methylation class by generating a distribution of the test statistic for the null distribution for comparison to the observed distribution. For continuous variables, the permutation test was run with the Kruskal-Wallis test statistic. For categorical variables we used the standard chi-square statistic for testing association between two categorical variables.

Locus-by-locus analysis

Associations between covariates and methylation at individual CpG loci were tested with a generalized linear model. The β-distribution of average β values was accounted for with a quasi-binomial logit link with an estimated scale parameter constraining the mean between 0 and 1, in a manner similar to that described by Hsuing et al. [43]. Array-wide scanning for CpG loci associations with sample type or covariate used false discovery rate estimation and Q-values computed by the qvalue package in R [44].

Multinomial logistic regression

Multinomial logistic regression was used to model methylation class while controlling for potential confounders. Referent class selection does not affect the underlying interpretation of the model and as class three was neither the largest, nor the smallest class, and was relatively hypomethylated it was chosen as the referent class. Because of the potentially large number of methylation classes, logistic regression coefficients were regularized using a ridge (L2) penalty, with coefficients for a common (non-intercept) covariate across outcome levels shrunk toward zero [34] the tuning parameter was selected by minimizing Bayesian information criterion.

Sequenom EpiTYPER methylation and RT–PCR

Spearman correlation coefficients and test P-values are reported for correlation between array and Sequenom methylation values. Tests for association between methylation and mRNA expression used relative mRNA expression versus array methylation average β stratified into two groups around 0.5 with the Kruskal-Wallis test statistic.

Supporting Information

Attachment 1

Attachment 2

Attachment 3

Attachment 4

Attachment 5

Attachment 6

Attachment 7


Zdroje

1. 2007 Global Cancer Facts & Figures 2007. Atlanta American Cancer Society

2. 2009 Breast Cancer Facts & Figures 2009–2010. Atlanta American Cancer Society

3. SingletarySE

ConnollyJL

2006 Breast cancer staging: working with the sixth edition of the AJCC Cancer Staging Manual. CA Cancer J Clin 56 37 47; quiz 50–31

4. EdgeSB

ByrdDR

ComptonCC

FritzAG

GreenFL

2009 AJCC Cancer Staging Atlas: Springer

5. McGuireWL

1978 Steroid receptors in human breast cancer. Cancer Res 38 4289 4291

6. SlamonDJ

ClarkGM

WongSG

LevinWJ

UllrichA

1987 Human breast cancer: correlation of relapse and survival with amplification of the HER-2/neu oncogene. Science 235 177 182

7. HarbeckN

NimmrichI

HartmannA

RossJS

CuferT

2008 Multicenter study using paraffin-embedded tumor tissue testing PITX2 DNA methylation as a marker for outcome prediction in tamoxifen-treated, node-negative breast cancer patients. J Clin Oncol 26 5036 5042

8. HartmannO

SpyratosF

HarbeckN

DietrichD

FassbenderA

2009 DNA methylation markers predict outcome in node-positive, estrogen receptor-positive breast cancer with adjuvant anthracycline-based chemotherapy. Clin Cancer Res 15 315 323

9. NovakP

JensenT

OshiroMM

WattsGS

KimCJ

2008 Agglomerative epigenetic aberrations are a common event in human breast cancer. Cancer Res 68 8616 8625

10. RonnebergJA

TostJ

SolvangHK

AlnaesGI

JohansenFE

2008 GSTP1 promoter haplotypes affect DNA methylation levels and promoter activity in breast carcinomas. Cancer Res 68 5562 5571

11. SinhaS

SinghRK

AlamN

RoyA

RoychoudhuryS

2008 Frequent alterations of hMLH1 and RBSP3/HYA22 at chromosomal 3p22.3 region in early and late-onset breast carcinoma: clinical and prognostic significance. Cancer Sci 99 1984 1991

12. SoaresJ

PintoAE

CunhaCV

AndreS

BaraoI

1999 Global DNA hypomethylation in breast carcinoma: correlation with prognostic factors and tumor progression. Cancer 85 112 118

13. VeeckJ

BektasN

HartmannA

KristiansenG

HeindrichsU

2008 Wnt signalling in human breast cancer: expression of the putative Wnt inhibitor Dickkopf-3 (DKK3) is frequently suppressed by promoter hypermethylation in mammary tumours. Breast Cancer Res 10 R82

14. VeeckJ

NoetzelE

BektasN

JostE

HartmannA

2008 Promoter hypermethylation of the SFRP2 gene is a high-frequent alteration and tumor-specific epigenetic marker in human breast cancer. Mol Cancer 7 83

15. MahoneyMC

BeversT

LinosE

WillettWC

2008 Opportunities and strategies for breast cancer prevention through risk reduction. CA Cancer J Clin 58 347 371

16. BredfeldtTG

GreathouseKL

SafeSH

HungMC

BedfordMT

2010 Xenoestrogen-Induced Regulation of EZH2 and Histone Methylation via Estrogen Receptor Signaling to PI3K/AKT. Mol Endocrinol 24 993 1006

17. Starlard-DavenportA

TryndyakVP

JamesSR

KarpfAR

LatendresseJR

2010 Mechanisms of epigenetic silencing of the Rassf1a gene during estrogen-induced breast carcinogenesis in ACI rats. Carcinogenesis 31 376 381

18. HousemanEA

ChristensenBC

MarsitCJ

KaragasMR

WrenschMR

2008 Model-based clustering of DNA methylation array data: a recursive-partitioning algorithm for high-dimensional data arising as a mixture of beta distributions. BMC Bioinformatics 9

19. JerussJS

MittendorfEA

TuckerSL

Gonzalez-AnguloAM

BuchholzTA

2008 Staging of breast cancer in the neoadjuvant setting. Cancer Res 68 6477 6481

20. MoriK

EnokidaH

KagaraI

KawakamiK

ChiyomaruT

2009 CpG hypermethylation of collagen type I alpha 2 contributes to proliferation and migration activity of human bladder cancer. Int J Oncol 34 1593 1602

21. HongHH

PischonN

SantanaRB

PalamakumburaAH

ChaseHB

2004 A role for lysyl oxidase regulation in the control of normal collagen deposition in differentiating osteoblast cultures. J Cell Physiol 200 53 62

22. MayoC

RenR

RichC

SteppMA

Trinkaus-RandallV

2008 Regulation by P2X7: epithelial migration and stromal organization in the cornea. Invest Ophthalmol Vis Sci 49 4384 4391

23. ProvenzanoPP

InmanDR

EliceiriKW

KnittelJG

YanL

2008 Collagen density promotes mammary tumor initiation and progression. BMC Med 6 11

24. ChernovAV

BaranovskayaS

GolubkovVS

WakemanDR

SnyderEY

2010 Microarray-based transcriptional and epigenetic profiling of matrix metalloproteinases, collagens and related genes in cancer. J Biol Chem 285 19647 19659

25. FisherB

BrownA

MamounasE

WieandS

RobidouxA

1997 Effect of preoperative chemotherapy on local-regional disease in women with operable breast cancer: findings from National Surgical Adjuvant Breast and Bowel Project B-18. J Clin Oncol 15 2483 2493

26. FisherB

BryantJ

WolmarkN

MamounasE

BrownA

1998 Effect of preoperative chemotherapy on the outcome of women with operable breast cancer. J Clin Oncol 16 2672 2685

27. CleatorS

HellerW

CoombesRC

2007 Triple-negative breast cancer: therapeutic options. Lancet Oncol 8 235 244

28. HillmanRS

SteinbergSE

1982 The effects of alcohol on folate metabolism. Annu Rev Med 33 345 354

29. HamajimaN

HiroseK

TajimaK

RohanT

CalleEE

2002 Alcohol, tobacco and breast cancer--collaborative reanalysis of individual data from 53 epidemiological studies, including 58,515 women with breast cancer and 95,067 women without the disease. Br J Cancer 87 1234 1245

30. KeyJ

HodgsonS

OmarRZ

JensenTK

ThompsonSG

2006 Meta-analysis of studies of alcohol and breast cancer with consideration of the methodological issues. Cancer Causes Control 17 759 770

31. LarssonSC

GiovannucciE

WolkA

2007 Folate and risk of breast cancer: a meta-analysis. J Natl Cancer Inst 99 64 76

32. ShrubsoleMJ

JinF

DaiQ

ShuXO

PotterJD

2001 Dietary folate intake and breast cancer risk: results from the Shanghai Breast Cancer Study. Cancer Res 61 7136 7141

33. TaoMH

ShieldsPG

NieJ

MarianC

AmbrosoneCB

2009 DNA promoter methylation in breast tumors: no association with genetic polymorphisms in MTHFR and MTR. Cancer Epidemiol Biomarkers Prev 18 998 1002

34. MarsitCJ

ChristensenBC

HousemanEA

KaragasMR

WrenschMR

2009 Epigenetic profiling reveals etiologically distinct patterns of DNA methylation in head and neck squamous cell carcinoma. Carcinogenesis 30 416 422

35. DumitrescuRG

ShieldsPG

2005 The etiology of alcohol-induced breast cancer. Alcohol 35 213 225

36. KwanML

AmbrosoneCB

LeeMM

BarlowJ

KrathwohlSE

2008 The Pathways Study: a prospective study of breast cancer survivorship within Kaiser Permanente Northern California. Cancer Causes Control 19 1065 1076

37. OehrliMD

QuesenberryCP

LeydenW

2006 Annual Report on Trends, Incidence, and Outcomes. Kaiser Permanente Northern California Cancer Registry

38. BibikovaM

LinZ

ZhouL

ChudinE

GarciaEW

2006 High-throughput DNA methylation profiling using universal bead arrays. Genome Res 16 383 393

39. EhrichM

NelsonMR

StanssensP

ZabeauM

LiloglouT

2005 Quantitative high-throughput analysis of DNA methylation patterns by base-specific cleavage and mass spectrometry. Proc Natl Acad Sci U S A 102 15785 15790

40. LiJ

SmythP

CahillS

DenningK

FlavinR

2008 Improved RNA quality and TaqMan Pre-amplification method (PreAmp) to enhance expression analysis from formalin fixed paraffin embedded (FFPE) materials. BMC Biotechnol 8 10

41. R Development CT 2007 R: A Language and Environment for Statistical Computing. Vienna, Austria R Foundation for Statistical Computing

42. ChristensenBC

HousemanEA

GodleskiJJ

MarsitCJ

LongackerJL

2009 Epigenetic profiles distinguish pleural mesothelioma from normal pleura and predict lung asbestos burden and clinical outcome. Cancer Res 69 227 234

43. HsiungDT

MarsitCJ

HousemanEA

EddyK

FurnissCS

2007 Global DNA methylation level in whole blood as a biomarker in head and neck squamous cell carcinoma. Cancer Epidemiol Biomarkers Prev 16 108 114

44. StoreyJ

TaylorJ

SiegmundD

2004 Strong control, conservative point estimation, and simultaneous conservative consistency of false discovery rates: A unified approach. J Royal Stat Soc Series B 187 205

Štítky
Genetika Reprodukční medicína

Článek vyšel v časopise

PLOS Genetics


2010 Číslo 7
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy Podcasty Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#