Advertisement
Research Articles

Racial Variation in Breast Tumor Promoter Methylation in the Carolina Breast Cancer Study

Kathleen Conway, Sharon N. Edmiston, Chiu-Kit Tse, Christopher Bryant, Pei Fen Kuan, Brionna Y. Hair, Eloise A. Parrish, Ryan May and Theresa Swift-Scanlan
DOI: 10.1158/1055-9965.EPI-14-1228 Published June 2015

Abstract

Background: African American (AA) women are diagnosed with more advanced breast cancers and have worse survival than white women, but a comprehensive understanding of the basis for this disparity remains unclear. Analysis of DNA methylation, an epigenetic mechanism that can regulate gene expression, could help to explain racial differences in breast tumor clinical biology and outcomes.

Methods: DNA methylation was evaluated at 1,287 CpGs in the promoters of cancer-related genes in 517 breast tumors of AA (n = 216) or non-AA (n = 301) cases in the Carolina Breast Cancer Study (CBCS).

Results: Multivariable linear regression analysis of all tumors, controlling for age, menopausal status, stage, intrinsic subtype, and multiple comparisons [false discovery rate (FDR)], identified seven CpG probes that showed significant (adjusted P < 0.05) differential methylation between AAs and non-AAs. Stratified analyses detected an additional four CpG probes differing by race within hormone receptor–negative (HR−) tumors. Genes differentially methylated by race included DSC2, KCNK4, GSTM1, AXL, DNAJC15, HBII-52, TUSC3, and TES; the methylation state of several of these genes may be associated with worse survival in AAs. TCGA breast tumor data confirmed the differential methylation by race and negative correlations with expression for most of these genes. Several loci also showed racial differences in methylation in peripheral blood leukocytes (PBL) from CBCS cases, indicating that these variations were not necessarily tumor-specific.

Conclusions: Racial differences in the methylation of cancer-related genes are detectable in both tumors and PBLs from breast cancer cases.

Impact: Epigenetic variation could contribute to differences in breast tumor development and outcomes between AAs and non-AAs. Cancer Epidemiol Biomarkers Prev; 24(6); 921–30. ©2015 AACR.

Introduction

Breast cancer, the most common cancer among women in the United States (1), is a heterogeneous disease with multiple clinical, histopathologic, and molecular subtypes, exhibiting different therapeutic responses and prognoses (2–4). Racial differences exist in its presentation and outcome, with African-American (AA) women, especially younger women, being diagnosed with more advanced cancers (5) and having worse survival than white women even after controlling for known prognostic factors or treatment (5–7). Prior work using protein markers to classify breast tumor intrinsic subtypes in the Carolina Breast Cancer Study (CBCS) found a higher prevalence of basal-like breast tumors among AA women compared with white women (3); however, AA women also have worse outcomes across all intrinsic subtypes (8).

In an effort to better understand the molecular factors contributing to breast cancer development, outcomes, and racial disparities, we evaluated DNA promoter methylation profiles in invasive breast tumors from the CBCS. An epigenetic modification to DNA that does not alter the nucleotide sequence (9), DNA methylation usually occurs as the addition of a methyl group to a cytosine within a CpG dinucleotide. Aberrant hypermethylation of CpG islands in tumor-suppressor genes can result in their silencing in cancer, while hypomethylation can lead to increased oncogene expression (9, 10). Such methylation changes and their effects on gene expression have the potential to influence breast tumor phenotypes and clinical outcomes (10). Several prior studies using targeted methylation approaches for a few genes have suggested that racial variation exists in methylation between AAs and non-AAs (or Caucasians; refs. 11–14), and a recent whole-genome study using the Illumina 450K platform supports these earlier reports (15).

In this study, use of a cancer-focused promoter methylation array to assess methylation in a large population-based series of breast tumors from AA and non-AA breast cancer cases in the CBCS revealed differential tumor methylation of several genes, and most of these were also differentially methylated in white blood cells, thereby enabling a more informed interpretation of which epigenetic differences might potentially contribute to breast cancer development, influence breast tumor chemosensitivity or the disparity in outcomes.

Materials and Methods

CBCS cases and specimens

The CBCS is a population-based, case–control study of incident invasive breast cancer in North Carolina. Details of the study design have been described previously (16). Randomized recruitment was used to oversample younger and AA cases to ensure that they comprised roughly half the study sample. Race/ethnicity was self-reported. Among non-AA cases, 97% self-reported as Caucasian (n = 291) and the other 3% included 4 cases with Hispanic ethnicity, 3 American Indians, 6 Asian/Native Islanders, and 1 other; therefore, we refer to this group as non-AA.

Blood samples were obtained for DNA extraction from peripheral blood leukocytes (PBL). Ancestry informative markers (AIM) to estimate the proportion of African versus European ancestry were previously evaluated from PBL DNA in the CBCS (17), AIMs genotypes were derived from 144 single-nucleotide polymorphisms (SNP), and an AIMs score ranging from 0 to 1.0 was generated for each case representing lower to higher African ancestry.

Formalin-fixed paraffin-embedded (FFPE) breast tumors were sectioned and mounted on slides, histopathologically reviewed, and tumor areas macrodissected for DNA extraction (18). Clinical data were obtained from medical records or histopathologic review of tumor tissue. Breast tumor intrinsic subtypes were identified using estrogen receptor (ER), progesterone receptor (PR), HER2, CK5, CK6, and EGFR protein markers (3). Hormone receptor (HR) positivity was considered to be ER+ and/or PR+, while tumors designated HR-negative (HR−) were both ER− and PR−.

Methylation analysis

DNA lysates prepared from breast tumors were sodium bisulfate–treated using the EZ DNA Methylation Gold Kit (Zymo Research). Methylation profiling was accomplished using the Illumina GoldenGate Cancer Panel I array (Illumina) on 517 breast tumors and 69 PBLs from cases, 61 of which were matched to individual cases. Methylation data were preprocessed using GenomeStudio Methylation software (Illumina). Methylation is represented by β, an estimate of the fraction of methylated DNA, and ranges from 0 (unmethylated) to 1.0 (fully methylated). The array interrogated 1,505 CpG sites located in the upstream regulatory region (promoter or exon 1) of 807 genes. Of the 1,505 CpG probes, a total of 1,287 CpGs were included in the final dataset after filtering (18), including probes previously identified as overlapping an underlying sequence variant (19). Compared with 163 cases not evaluated due to inadequate quantity or quality of DNA (of 680 total with tumors in CBCS phase I), the 517 breast cancer cases evaluated were younger (P = 0.03) but did not differ in other characteristics. Array data were deposited to Gene Expression Omnibus under accession number GSE51557.

Statistical analyses

Statistical analyses were carried out using the R package (http://www.r-project.org/) or SAS 9.3. CpGs differentially methylated between breast tumors of AAs and non-AAs were identified via generalized linear models (GLM), modeling the methylation value with logit link function, and including race (AA vs. non-AA) as a predictor, while also adjusting for age, menopausal status, stage, and intrinsic subtype (20). P values were adjusted for multiple comparisons via the Benjamini–Hochberg method for controlling the false discovery rate (FDR). A GLM model to compare methylation by race was performed first for all breast tumors, and then was repeated with stratification of breast tumors on HR status. The box and whisker plots were constructed to display the distribution of methylation β values for each group, and included the mean, median, and interquartile range of methylation values, and the group minimum and maximum values. Univariate analysis of methylation by race in PBLs was determined by the Student t test. Correlation of methylation in paired tumors and PBLs from cases was determined by the Spearman rank-correlation coefficient. The Kaplan–Meier plots and the log-rank P values were used to illustrate disease-specific survival according to methylation level (high, intermediate, and low) for select probes within the AA or non-AA case groups. For each probe, cutoff points were based on the tertiles in the combined AA/non-AA dataset.

Validation studies

Independent validation of racial differences in breast tumor methylation and correlations of methylation with gene expression were conducted using breast tumor data from The Cancer Genome Atlas (TCGA; ref. 21). Methylation analysis in TCGA was performed using the Illumina Infinium 450K methylation array, while gene expression data were generated using RNA sequencing. Only 371 of the 1,505 CpG probes interrogated on the Illumina GoldenGate methylation platform exactly match those on the 450K methylation array; however, prior analysis of 450K probes showed that methylation levels were concordant for matched probes between distinct methylation platforms (22). Therefore, we compared methylation between AA and white patients at probes from the 450K array in TCGA that interrogated direct match CpG sites when available, or those closest within 200 bp upstream or downstream of the GoldenGate probes of interest. We also tested associations between methylation and gene expression for these probes. Two separate subsets of TCGA data were used for these analyses. For comparisons of breast tumor methylation, β values for each probe of interest from the 450K array were compared between 42 AA women and 291 white women in the TCGA dataset. T tests were used to identify probes that differed by race at P < 0.05. TCGA data from 581 breast cancer patients were used to examine relationships between methylation and gene expression. Pearson correlation coefficients were calculated in all 581 TCGA breast tumors, or within HR+ or HR− subsets, based on RNAseq (Illumina) log2 RSEM gene normalized expression values with methylation β values for 450K CpG probes, with significance set at P < 0.05. Independent validation of racial differences in methylation was also performed using an Illumina 450K methylation dataset in Gene Expression Omnibus on lymphoblastoid cell lines from healthy female AA (n = 80) and Caucasian American subjects (n = 49; accession #GSE36369; ref. 23). Moreover, technical validation was conducted to compare gene methylation from the GoldenGate array with that from quantitative methylation–specific polymerase chain reaction (Q-MSP; ref. 24).

Results

Characteristics of breast cancer cases

Characteristics of AA and non-AA cases and their tumors evaluated for DNA methylation are detailed in Table 1. Breast cancer cases were mostly early stage (>86% stages I or II). AA cases were older, more frequently postmenopausal, and more likely to have HR−, high-grade, or basal-like breast tumors compared with non-AAs, as previously reported (3). PBL samples derived from a subset of AA and non-AA cases did not differ in age (P = 0.44) or menopausal status (P = 0.12).

View this table:
Table 1.

Demographic, tumor, and clinical characteristics of the AA and non-AA breast cancer cases evaluated for DNA methylation

Racial differences in breast tumor methylation

Comparison of β values for 1,287 CpG probes in 517 breast tumors using GLM while controlling for age, menopausal status, stage, intrinsic subtype, and multiple comparisons (FDR) identified a total of 24 CpG probes that differed significantly at q < 0.05 between AA (n = 216) and non-AA (n = 301) cases (Table 2). The majority of these probes (n = 17), however, were considered technically ambiguous or unreliable because they were previously reported to overlap a site of a potential SNP or repeat (19), or based on review of updated probe target sequence information in NCBI databases (e.g., Ensembl, dbSNP, and Blast) were subsequently found to overlap underlying sequence variants or show ancestral population differences of >10% (as detailed in Supplementary Table S1). In total, seven probes showing differential methylation by race were retained for further analysis, including four that had no known underlying variant (DSC2_E90_F, KCNK4_E3_F, GSTM1_P266_F, and HBII-52_E142_F) and three that overlapped a single SNP but that was unlikely to impact probe performance due to its location toward the end of the target sequence and the lack of evidence for ancestral differences in allele frequencies (AXL_P223_R, DNAJC15_E26_R, and TES_P182_F; Table 2). Although the GSTM1 gene is polymorphic for a deletion variant, this deletion occurs within the coding sequence and does not involve the probe target region in the promoter (25). Importantly, racial variation in methylation for these seven loci was unlikely to be appreciably related to the known racial differences in breast tumor subtype distribution or differences in methylation patterns between the major intrinsic subtypes because we adjusted for IHC-based subtype, stage, age, and menopausal status in the GLM. Racial variation in breast tumor methylation was evident for these top CpG probes whether race was self-reported or defined by AIMs (Supplementary Fig. S1).

View this table:
Table 2.

CpG probes differentially methylated by race in multivariable GLM analysis of all breast tumors

Stratified GLM was subsequently performed to determine whether the racial differences initially observed in CpG methylation were evident within HR-defined subtypes. As shown in Supplementary Table S2, four additional probes that did not previously meet the threshold of q < 0.05 showed significant racial differences among cases with HR− tumors, including TUSC3_E29_R, RAF1_P330_F, SMARCA3_P17_R, and IMPACT_P234_R. Notably, more CpG loci showed racial variation within HR− tumors than HR+ tumors. Only HBII-52_E142_F showed significant differential methylation within HR+ tumors of AAs versus non-AAs. The box and whisker plots summarizing the distribution of methylation β values for probes significantly varying by race (at q < 0.05), overall or within HR-defined tumor subsets, are shown in Fig. 1A–C. Technical validation comparing methylation obtained from the GoldenGate array with Q-MSP confirmed the comparability of the methylation measurements obtained by these two methods (Supplementary Fig. S2).

Figure 1.
Figure 1.

The box and whisker plots showing CpG probes exhibiting racial differences in breast tumor methylation. CpG probes shown differed significantly (at q < 0.05) in tumor methylation between AA and non-AA cases in GLM after controlling for age, menopausal status, stage, intrinsic subtype, and multiple comparisons (FDR). Each box plot shows the median β-value (dark bar within box), the interquartile range (outer boundaries of box), and the whiskers extend to the group minimum and maximum β values. Probes shown either do not overlap a sequence variant or overlap a single SNP that is likely to have negligible impact on probe performance. Probes show significant racial differences in methylation in all CBCS breast tumors (n = 517; A), HR+ breast tumors (n = 337; B), and HR− breast tumors (n = 163; C). Shown are adjusted P values and q values from GLM. D, validation of racial variation in methylation for candidate genes in TCGA breast tumors. Data are shown for 450K CpG probes that exactly matched or were closest to (within 200 bp) probes on the GoldenGate array. P values were determined by the t test.

Although there were too few genes differentially methylated by racial group to perform formal gene ontology or pathway analyses, several of the genes showing epigenetic racial differences have roles in DNA repair, transcription, or mediate other DNA interactions (GSTM1, DNAJC15, and SMARCA3/HLTF), are involved in cell adhesion (DSC2 and TES), or are kinases important in signal transduction (AXL and RAF1).

Because AAs experience more adverse breast cancer outcomes than non-AAs, we determined whether the top candidate CpG markers might be differentially associated with disease-specific survival among AAs versus non-AAs. The Kaplan–Meier plots in Supplementary Fig. S3 indicate that a low or intermediate (vs. high) level of methylation at several probes (HBII-52_E142_F, TES_P182_F, DNAJC15_P65_F, DNAJC15_E26_R, and DSC2_E90_F) was marginally associated with worse survival in AAs but not in non-AAs.

Validation in TCGA breast tumors

Breast tumor methylation and gene expression data from TCGA (21) were used to validate both our methylation findings and to test for correlations with expression for genes differentially methylated by race in the CBCS. Comparison of 450K methylation β values for probes that either directly matched or were located within 200 bp of GoldenGate probes in 42 AA and 291 white TCGA patients confirmed our findings of racial differences in methylation for AXL, GSTM1, KCNK4, and DNAJC15 (Table 3 and Fig. 1D). Methylation was also significantly inversely correlated with gene expression for the majority of our candidate genes (Supplementary Table S3).

View this table:
Table 3.

TCGA validation of differential breast tumor methylation in genes showing racial variation in the CBCS

Racial differences in methylation in PBLs or matched tumor/PBL pairs

Methylation patterns are considered to be tissue-specific, yet, if methylation at certain loci differs by race rather than tissue type, we expect that these differences might also be detectable in other normal cells, even if derived from a different tissue. To determine whether the racial differences in CpG methylation were restricted to tumors or occurred more broadly, we compared methylation in PBLs from 29 AAs and 40 non-AA cases matched on age and menopausal status. As shown in Table 4, the majority of CpG markers that were differentially methylated by race in tumors also differed significantly in PBLs, and the directionality of the difference in PBLs reflected that of the tumors. For example, similar to the pattern in tumors, DSC2_E90_F, GSTM1_P266_F, KCNK4_E3_F, and AXL_P223_R all showed higher methylation in PBLs of AAs than in non-AAs even though the absolute levels of methylation varied somewhat between tumor and PBL. Differential methylation by race was also observed for CpG sites in AXL, GSTM1, DSC2, and DNAJC15 in an independent series of lymphoblastoid cell lines from 80 female AAs compared with 49 female Caucasian Americans described by Heyn and colleagues (Supplementary Table S4; ref. 23).

View this table:
Table 4.

Differential methylation in PBLs of AA versus non-AA cases in the CBCS

We also compared methylation between matched pairs of breast tumors and PBLs within individual AA and non-AA cases (n = 61) for several probes exhibiting the largest methylation differences by race. Methylation in tumor/PBL pairs was significantly correlated (with P < 0.05) for the majority of probes (not shown). The box and whisker plots in Fig. 2A show that among the 61 cases for whom tumor and PBL samples are matched, the methylation differences observed in tumors are also evident within PBLs. The bar chart in Fig. 2B provides a comparison of individual β values in matched tumors and PBLs from breast cancer cases, further highlighting the differential methylation detected between AAs and non-AAs.

Figure 2.
Figure 2.

Methylation in matched pairs of breast tumors and PBLs from AA and non-AA cases in the CBCS. A, the box and whisker plots showing the distribution of methylation β values in 61 matched pairs of breast tumors and PBLs from AA and non-AA cases. B, bar charts showing individual methylation β values of matched tumor (black) and PBL (gray) pairs by race (26 AA, 35 non-AA) and age group (40, <50 years; 21, 50+ years).

Discussion

In this study, using a methylation array targeted to the promoters of genes important in cancer, we identified several genes differing in tumor DNA methylation between AAs and non-AAs even after controlling for potential confounders, including intrinsic subtype distributions that are known to vary by race (3). Previous studies, including the CBCS, have shown that breast tumor methylation varies between intrinsic subtypes (10, 18, 21, 26) or HR status (27), and by clinicopathologic characteristics such as tumor size and grade (18). Importantly, most of the CpG loci differing by race in CBCS breast tumors were also similarly differentially methylated in normal white blood cells of cases, although the absolute levels of methylation were not necessarily the same. A recent study that used the whole-genome Illumina 450K methylation array to assess racial differences in breast tumor methylation excluded probes that overlapped SNPs or that were otherwise ambiguous (15); however, we opted to retain such probes but to note their status. The impact of a probe that overlaps a target sequence potentially containing a SNP is difficult to predict, as it would depend on the SNP frequency in a given population as well as the probe sequence and precise location of the underlying SNP. Because many probes on the Illumina methylation arrays often overlap multiple CpG sites, the effect of an additional SNP-based single-nucleotide mismatch on probe binding and methylation detection may be minimal. Consistent with this, we have previously observed in the CBCS that some probes overlapping a SNP performed similarly in their detection of methylation (i.e., generated similar β values) to nearby probes in the same gene that do not (22). Genotyping for SNPs located within gene promoters was not performed in the CBCS and thus we could not rule out that racial variation in methylation might be related to the effects of underlying SNPs for certain probes; however, for most probes that overlapped a SNP, the variant alleles were relatively rare (minor allele frequency <10%) but in some cases differed by ancestry. Thus, to focus on CpG probes that would most reliably detect methylation differences (rather than underlying SNPs), we reviewed all probe target sequences against NCBI databases, assessed reported racial differences in allelic frequencies if available, and then removed those for which the underlying feature was likely to impede a reliable measurement of methylation. In the end, a total of 11 CpG probes were retained that were likely to provide a reasonably accurate assessment of racial variation in methylation.

Several prior studies using targeted methylation approaches to examine a small number of candidate genes known to be involved in breast cancer reported differences in breast tumor methylation between AAs and those of European ancestry. Consistent with our findings of more methylation differences by race within HR− tumors, Mehrotra and colleagues (14) found significantly higher methylation in AAs than in Caucasian woman for four genes within a five-gene marker panel (HIN-1, TWIST1, CCND2, RASSF1A, RARB); these differences were only evident within ER− tumors and among women diagnosed before age 50 years. Similarly, Wang and colleagues (13) examined an eight-marker panel (p16, RASSF1A, RARβ2, ESR1, CDH13, HIN1, SFRP1, and LINE1) and reported differential methylation for CDH13 between ER− tumors from younger AAs and European-American (EA) women. Using the Illumina 450K whole-genome methylation array, the same array platform used in TCGA, Ambrosone and colleagues (15) also reported that many more loci were differentially methylated in AAs compared with EAs in ER− than in ER+ breast cancers. The genes exhibiting differential promoter methylation by race in the CBCS are unique, and not among the candidate markers examined previously or the top loci reported from the 450K analysis in the Ambrosone study (15). Importantly, it should be noted that our analyses controlled for differences in intrinsic subtype, stage, and age, whereas these prior studies did not. Therefore, it is unclear whether some loci reported previously as varying by race or ancestry may actually reflect differences in tumor intrinsic subtype distributions between the racial groups, such as within ER− tumors that are comprised the basal-like, claudin-low, and HER2+/ER− subtypes. Intriguingly, several genes showing racial variation in tumor methylation (HBII-52, TES, DNAJC15, and DSC2) may be associated with differences in survival among AAs, while no such differences were observed among non-AAs. Further work is needed to confirm these results and to assess the possibility that underlying copy number or other changes contributed to these findings.

In the CBCS, we also observed racial differences in methylation in normal PBL samples from cases, with the overall patterns of tumor and PBL methylation for these probes being correlated, although the absolute levels of methylation in tumor versus PBL differed for some loci. Differences in the quality of DNA derived from FFPE tumors versus PBLs could potentially contribute to such methylation differences; however, PBL samples did not consistently show higher or lower methylation levels than in tumors. Although variation in PBL methylation profiles primarily reflect the admixture of white blood cell populations (28, 29), we still observed epigenetic differences by race. The finding of racial variation in methylation in both tumor and PBLs suggests underlying ancestral differences in the epigenetic state of some genes, with about half of differentially methylated loci being potentially explained by known sequence variations that could disrupt probe binding and methylation measurements. This is consistent with several prior studies noting that such genetic variants contributed to a substantial portion of racial differences in DNA methylation (23, 30–32). Racial differences in DNA methylation patterns in normal tissues have also been detected in PBLs from women in a multiethnic New York City Birth Cohort (33), in umbilical cord blood from newborns (34), and in breast tissue from reduction mammoplasty (11), raising the question of whether such normal epigenetic variation contributes to cancer risk. In a study of normal human prostate tissue and prostate cancer, several genes exhibited differential methylation between AAs and EAs (12).

The genes displaying racial variation in tumor methylation in the CBCS included several involved in transcription or interactions with DNA (DNAJC15, SMARCA3, or HLTF), signal transduction (AXL and RAF1), carcinogen detoxification (GSTM1), cell adhesion (DSC2), or are known tumor suppressors (TUSC3 and TES). Desmocollin 2 (DSC2), a membrane glycoprotein involved in cell–cell adhesion and maintenance of normal epithelial architecture is an independent prognostic marker for esophageal (35) and pancreatic cancers (36), and high DSC2 expression may be a marker of basal-like breast tumors (37). The carcinogen detoxifying enzyme, GSTM1, exhibits an intragenic deletion variant resulting in loss of enzyme activity (25), elevated cancer risk, and shows racial variation between blacks and whites (38); however, whether the variant is linked with promoter methylation state is unclear. Testin (TES) and TUSC3 are tumor-suppressor genes that are silenced, in part, via promoter methylation (39); loss of expression of TES is an independent poor prognostic marker in breast cancer (40), while silencing of TUSC3 may be a poor prognostic factor in ovarian cancer (41). AXL is a receptor tyrosine kinase and an epithelial-to-mesenchymal transition (EMT)–induced regulator of breast cancer metastasis and patient survival (42–45). AXL expression is regulated, in part, by Sp1 transcription factor binding and methylation (42, 44, 46). In contrast to most other genes showing epigenetic racial variation in the CBCS, AXL methylation was not inversely correlated with mRNA expression in the TCGA tumor set, most likely because the 450K AXL CpG probes analyzed in relation to gene expression appear to be outside of these critical Sp1 binding sites. Intriguingly, AXL was reported to be aberrantly hypermethylated in association with prenatal tobacco smoke exposure (47), suggesting that epigenetic racial variation may be associated with lifestyle or environmental exposures and could potentially contribute to cancer risk.

The strengths of this study include the large and well-characterized population-based series of mostly early-stage breast cancer cases, providing the power to detect even modest differences in methylation. Epigenetic racial differences for several of the genes in the CBCS were independently confirmed in TCGA breast tumors although the number of AAs in this dataset was limited. Moreover, the availability of demographic, clinical and subtype information allowed us to control for factors that are known to vary by patient or tumor subsets. However, our study was limited by use of the targeted cancer-focused array rather than a whole-genome methylation array.

The implications of the racial differences in methylation discovered in this study for breast cancer development, progression, or outcomes in AA or non-AA women are unclear, and further work is needed to confirm these findings. It is of interest that the epigenetic differences were most evident within HR− tumors despite controlling for intrinsic subtype; however, it is possible that these racial variations signify different proportions of HR− subtypes in AA and non-AA women. Interestingly, most of the genes exhibiting differential tumor methylation by race in the CBCS (DSC2, AXL, DNAJC15, TES, and TUSC3) are among the genes defining subclasses of triple-negative breast cancer as reported by Lehmann and colleagues (48). Whether the prevalence of triple-negative subtypes, which were not defined within the CBCS, vary between AAs and non-AAs is presently unknown. The results of this study highlight the possibility that methylation patterns of breast tumors, particularly the more aggressive HR− tumors, may differ by ancestry, and that the racial differences may not be a tumor-specific phenomenon, suggesting that such variations could also contribute to cancer risk. Clearly, more research is needed to validate these findings and determine whether they contribute to the known racial disparity in breast cancer survival.

Disclosure of Potential Conflicts of Interest

No potential conflicts of interest were disclosed.

In Memoriam

This work is dedicated to the memory of our colleague, the late Dr. Robert Millikan, principal investigator of the CBCS.

Authors' Contributions

Conception and design: K. Conway, S.N. Edmiston

Development of methodology: K. Conway, S.N. Edmiston

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): K. Conway, S.N. Edmiston, T. Swift-Scanlan

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): K. Conway, S.N. Edmiston, C.-K. Tse, C. Bryant, P.F. Kuan, B.Y. Hair, R. May, T. Swift-Scanlan

Writing, review, and/or revision of the manuscript: K. Conway, S.N. Edmiston, C. Bryant, B.Y. Hair, T. Swift-Scanlan

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): K. Conway, S.N. Edmiston, E.A. Parrish

Study supervision: K. Conway

Grant Support

This research was supported by grants to K. Conway from Susan G. Komen Foundation (grant #KG081397) and from the University Cancer Research Fund of North Carolina. The CBCS was also funded by the University Cancer Research Fund of North Carolina and the National Cancer Institute Specialized Program of Research Excellence (SPORE) in Breast Cancer (NIH/NCI P50-CA58223).

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.

Acknowledgments

The authors thank Dr. Andy Olshan for critical review of this article and the participants and study staff of the CBCS for their continuing support of the study.

Footnotes

  • Note: Supplementary data for this article are available at Cancer Epidemiology, Biomarkers & Prevention Online (http://cebp.aacrjournals.org/).

  • Received October 31, 2014.
  • Revision received February 11, 2015.
  • Accepted March 3, 2015.

References

  1. 1.↵
    1. DeSantis C,
    2. Siegel R,
    3. Bandi P,
    4. Jemal A
    . Breast cancer statistics, 2011. CA Cancer J Clin 2011;61:409–18.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Perou CM,
    2. Sorlie T,
    3. Eisen MB,
    4. van de Rijn M,
    5. Jeffrey SS,
    6. Rees CA,
    7. et al.
    Molecular portraits of human breast tumors. Nature 2000;406:747–52.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Carey LA,
    2. Perou CM,
    3. Livasy CA,
    4. Dressler LG,
    5. Cowan D,
    6. Conway K,
    7. et al.
    Race, breast cancer subtypes, and survival in the Carolina Breast Cancer Study. JAMA 2006;295:2492–502.
    OpenUrlCrossRefPubMed
  4. 4.↵
    1. Carey LA,
    2. Dees EC,
    3. Sawyer L,
    4. Gatti L,
    5. Moore DT,
    6. Collichio F,
    7. et al.
    The triple negative paradox: primary tumor chemosensitivity of breast cancer subtypes. Clin Cancer Res 2007;13:2329–34.
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    1. Smigal C,
    2. Jemal A,
    3. Ward E,
    4. Cokkinides V,
    5. Smith R,
    6. Howe H,
    7. et al.
    Trends in breast cancer by race and ethnicity: update 2006. CA Cancer J Physicians 2006;56:168–83.
    OpenUrl
  6. 6.↵
    1. Maskarinec G,
    2. Sen C,
    3. Koga K,
    4. Conroy SM
    . Ethnic differences in breast cancer survival: status and determinants. Womens Health 2011;7:677–87.
    OpenUrlCrossRef
  7. 7.↵
    1. Shavers VL,
    2. Brown ML
    . Racial and ethnic disparities in the receipt of cancer treatment. J Natl Cancer Inst 2002;94:334–57.
    OpenUrlAbstract/FREE Full Text
  8. 8.↵
    1. O'Brien KM,
    2. Cole SR,
    3. Tse CK,
    4. Perou CM,
    5. Carey LA,
    6. Foulkes WD,
    7. et al.
    Intrinsic breast tumor subtypes, race, and long-term survival in the Carolina Breast Cancer Study. Clin Cancer Res 2010;16:6100–10.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Herman JG
    . Epigenetic changes in cancer and preneoplasia. Cold Spring Harb Symp Quant Biol 2005;70:329–33.
    OpenUrlAbstract/FREE Full Text
  10. 10.↵
    1. Kamalakaran S,
    2. Varadan V,
    3. Giercksky Russnes HE,
    4. Levy D,
    5. Kendall J,
    6. Janevski A,
    7. et al.
    DNA methylation patterns in luminal breast cancers differ from non-luminal subtypes and can identify relapse risk independent of other clinical variables. Mol Oncol 2011;5:77–92.
    OpenUrlCrossRefPubMed
  11. 11.↵
    1. Dumitrescu RG,
    2. Marian C,
    3. Krishnan SS,
    4. Spear SL,
    5. Kallakury BVS,
    6. Perry DJ,
    7. et al.
    Familial and racial determinants of tumour suppressor genes promoter hypermethylation in breast tissues from healthy women. J Cell Mol Med 2010;14:1468–75.
    OpenUrlPubMed
  12. 12.↵
    1. Kwabi-Addo B,
    2. Wang S,
    3. Chung W,
    4. Jelinek J,
    5. Patierno SR,
    6. Wang B,
    7. et al.
    Identification of differentially methylated genes in normal prostate tissues from African American and Caucasian men. Clin Cancer Res 2010;16:3539–47.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    1. Wang S,
    2. Dorsey TH,
    3. Terunuma A,
    4. Kittles RA,
    5. Ambs S,
    6. Kwabi-Addo B
    . Relationship between tumor DNA methylation status and patient characteristics in African-American and European-American women with breast cancer. PLoS ONE 2012;7:e37928.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Mehrotra J,
    2. Ganpat MM,
    3. Kanaan Y,
    4. Fackler MJ,
    5. McVeigh M,
    6. Lahti-Domenici J,
    7. et al.
    Estrogen receptor/progesterone receptor-negative breast cancers of young African-American women have a higher frequency of methylation of multiple genes than those of Caucasian women. Clin Cancer Res 2004;10:2052–7.
    OpenUrlAbstract/FREE Full Text
  15. 15.↵
    1. Ambrosone CB,
    2. Young AC,
    3. Sucheston LE,
    4. Wang D,
    5. Yan L,
    6. Liu S,
    7. et al.
    Genome-wide methylation patterns provide insight into differences in breast tumor biology between American women of African and European ancestry. Oncotarget 2014;5:237–48.
    OpenUrlPubMed
  16. 16.↵
    1. Newman B,
    2. Moorman PG,
    3. Millikan R,
    4. Qaqish BF,
    5. Geradts J,
    6. Aldrich TE,
    7. et al.
    The Carolina Breast Cancer Study: integrating population-based epidemiology and molecular biology. Breast Cancer Res Treat 1995;35:51–60.
    OpenUrlCrossRefPubMed
  17. 17.↵
    1. Barnholtz-Sloan JS,
    2. Shetty PB,
    3. Guan X,
    4. Nyante SJ,
    5. Luo J,
    6. Brennan DJ,
    7. et al.
    FGFR2 and other loci identified in genome-wide association studies are associated with breast cancer in African-American and younger women. Carcinogenesis 2010;31:1417–23.
    OpenUrlAbstract/FREE Full Text
  18. 18.↵
    1. Conway K,
    2. Edmiston SN,
    3. May R,
    4. Kuan PF,
    5. Chu H,
    6. Bryant C,
    7. et al.
    DNA methylation profiling in the Carolina Breast Cancer Study defines cancer subclasses differing in clinicopathologic characteristics and survival. Breast Cancer Res 2014;16:450.
    OpenUrlPubMed
  19. 19.↵
    1. Byun HM,
    2. Siegmund KD,
    3. Pan F,
    4. Weisenberger DJ,
    5. Kanel G,
    6. Laird PW,
    7. et al.
    Epigenetic profiling of somatic tissues from human autopsy specimens identifies tissue- and individual-specific DNA methylation patterns. Hum Mol Genet 2009;18:4808–17.
    OpenUrlAbstract/FREE Full Text
  20. 20.↵
    1. Wedderburn,
    2. Robert WM
    . Quasi-likelihood functions, generalized linear models, and the Gauss-Newton method. Biometrika 1974;61:439–47.
    OpenUrlAbstract/FREE Full Text
  21. 21.↵
    Cancer Genome Atlas Network. Comprehensive molecular portraits of human breast tumours. Nature 2012;490:61–70.
    OpenUrlCrossRefPubMed
  22. 22.↵
    1. Bardowell SA,
    2. Parker J,
    3. Fan C,
    4. Crandell J,
    5. Perou CM,
    6. Swift-Scanlan T
    . Differential methylation relative to breast cancer subtype and matched normal tissue reveals distinct patterns. Breast Cancer Res Treat 2013;142:365–80.
    OpenUrlPubMed
  23. 23.↵
    1. Heyn H,
    2. Moran S,
    3. Hernando-Herraez I,
    4. Sayols S,
    5. Gomez A,
    6. Sandoval J,
    7. et al.
    DNA methylation contributes to natural human variation. Genome Res 2013;23:1363–72.
    OpenUrlAbstract/FREE Full Text
  24. 24.↵
    1. Thomas N E,
    2. Slater NA,
    3. Edmiston SN,
    4. Zhou X,
    5. Kuan PF,
    6. Groben PA,
    7. et al.
    DNA methylation profiles in primary cutaneous melanomas are associated with clinically significant pathologic features. Pigment Cell Melanoma Res 2014;27:1097–1105.
    OpenUrlPubMed
  25. 25.↵
    1. Xu S,
    2. Wang Y,
    3. Roe B,
    4. Pearson WR
    . Characterization of the human class Mu glutathione S-transferase gene cluster and the GSTM1 deletion. J Biol Chem 1998;273:3517–27.
    OpenUrlAbstract/FREE Full Text
  26. 26.↵
    1. Rønneberg JA,
    2. Fleischer T,
    3. Solvang HK,
    4. Nordgard SH,
    5. Edvardsen H,
    6. Potapenko I,
    7. et al.
    Methylation profiling with a panel of cancer related genes: association with estrogen receptor, TP53 mutation status and expression subtypes in sporadic breast cancer. Mol Oncol 2011;5:61–76.
    OpenUrlCrossRefPubMed
  27. 27.↵
    1. Fackler MJ,
    2. Umbricht CB,
    3. Williams D,
    4. Argani P,
    5. Cruz LA,
    6. Merino VF,
    7. et al.
    Genome-wide methylation analysis identifies genes specific to breast cancer hormone receptor status and risk of recurrence. Cancer Res 2011;71:6195–207.
    OpenUrlAbstract/FREE Full Text
  28. 28.↵
    1. Houseman EA,
    2. Accomando WP,
    3. Koestler DC,
    4. Christensen BC,
    5. Marsit CJ,
    6. Nelson HH,
    7. et al.
    DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics 2012;13:86.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Koestler DC,
    2. Christensen B,
    3. Karagas MR,
    4. Marsit CJ,
    5. Langevin SM,
    6. Kelsey KT,
    7. et al.
    Blood-based profiles of DNA methylation predict the underlying distribution of cell types: a validation analysis. Epigenetics 2013;8:816–26.
    OpenUrlCrossRefPubMed
  30. 30.↵
    1. Zhang X,
    2. Moen EL,
    3. Liu C,
    4. Mu W,
    5. Gamazon ER,
    6. Delaney SM,
    7. et al.
    Linking the genetic architecture of cytosine modifications with human complex traits. Hum Mol Genet 2014;23:5893–905.
    OpenUrlAbstract/FREE Full Text
  31. 31.↵
    1. Shoemaker R,
    2. Deng J,
    3. Wang W,
    4. Zhang K
    . Allele-specific methylation is prevalent and is contributed by CpG-SNPs in the human genome. Genome Res 2010;20:883–9.
    OpenUrlAbstract/FREE Full Text
  32. 32.↵
    1. Moen EL,
    2. Zhang X,
    3. Mu W,
    4. Delaney SM,
    5. Wing C,
    6. McQuade J,
    7. et al.
    Genome-wide variation of cytosine modifications between European and African populations and the implications for complex traits. Genetics 2013;194:987–96.
    OpenUrlAbstract/FREE Full Text
  33. 33.↵
    1. Terry MB,
    2. Ferris JS,
    3. Pilsner R,
    4. Flom JD,
    5. Tehranifar P,
    6. Santella RM,
    7. et al.
    Genomic DNA methylation among women in a multiethnic New York City birth cohort. Cancer Epidemiol. Biomarkers Prev 2008;17:2306–10.
    OpenUrlAbstract/FREE Full Text
  34. 34.↵
    1. Adkins RM,
    2. Krushkal J,
    3. Tylavsky FA,
    4. Thomas F
    . Racial differences in gene-specific DNA methylation levels are present at birth. Birth Defects Research Part A. Clin Mol Teratol 2011;91:728–36.
  35. 35.↵
    1. Fang WK,
    2. Gu W,
    3. Li EM,
    4. Wu ZY,
    5. Shen ZY,
    6. Shen JH,
    7. et al.
    Reduced membranous and ectopic cytoplasmic expression of DSC2 in esophageal squamous cell carcinoma: an independent prognostic factor. Hum Pathol 2010;41:1456–65.
    OpenUrlCrossRefPubMed
  36. 36.↵
    1. Hamidov Z,
    2. Altendorf-Hofmann A,
    3. Chen Y,
    4. Settmacher U,
    5. Petersen I,
    6. Knösel T
    . Reduced expression of desmocollin 2 is an independent prognostic biomarker for shorter patients survival in pancreatic ductal adenocarcinoma. J Clin Pathol 2011;64:990–4.
    OpenUrlAbstract/FREE Full Text
  37. 37.↵
    1. Culhane AC,
    2. Quackenbush J
    . Confounding effects in “A six-gene signature predicting breast cancer lung metastasis”. Cancer Res 2009;69:7480–5.
    OpenUrlAbstract/FREE Full Text
  38. 38.↵
    1. Bell DA,
    2. Taylor JA,
    3. Paulson DF,
    4. Robertson CN,
    5. Mohler JL,
    6. Lucier GW
    . Genetic risk and carcinogen exposure: a common inherited defect of the carcinogen-metabolism gene glutathione S-transferase M1 (GSTM1) that increases susceptibility to bladder cancer. J Natl Cancer Inst 1993;85:1159–64.
    OpenUrlAbstract/FREE Full Text
  39. 39.↵
    1. Tobias ES,
    2. Hurlstone AF,
    3. MacKenzie E,
    4. McFarlane R,
    5. Black DM
    . The TES gene at 7q31.1 is methylated in tumours and encodes a novel growth-suppressing LIM domain protein. Oncogene 2001;20:2844–53.
    OpenUrlCrossRefPubMed
  40. 40.↵
    1. Zhu J,
    2. Li X,
    3. Kong X,
    4. Moran MS,
    5. Su P,
    6. Haffty BG,
    7. et al.
    Testin is a tumor suppressor and prognostic marker in breast cancer. Cancer Sci 2012;103:2092–101.
    OpenUrlPubMed
  41. 41.↵
    1. Pils D,
    2. Horak P,
    3. Vanhara P,
    4. Anees M,
    5. Petz M,
    6. Alfanz A,
    7. et al.
    Methylation status of TUSC3 is a prognostic factor in ovarian cancer. Cancer 2013;119:946–54.
    OpenUrlPubMed
  42. 42.↵
    1. Choufani S,
    2. Shapiro JS,
    3. Susiarjo M,
    4. Butcher DT,
    5. Grafodatskaya D,
    6. Lou Y,
    7. et al.
    A novel approach identifies new differentially methylated regions (DMRs) associated with imprinted genes. Genome Res 2011;21:465–76.
    OpenUrlAbstract/FREE Full Text
  43. 43.↵
    1. Breton CV,
    2. Salam MT,
    3. Gilliland FD
    . Heritability and role for the environment in DNA methylation in AXL receptor tyrosine kinase. Epigenetics 2011;6:895–8.
    OpenUrlCrossRefPubMed
  44. 44.↵
    1. Gjerdrum C,
    2. Tiron C,
    3. Høiby T,
    4. Stefansson I,
    5. Haugen H,
    6. Sandal T,
    7. et al.
    Axl is an essential epithelial-to-mesenchymal transition-induced regulator of breast cancer metastasis and patient survival. Proc Natl Acad Sci U S A 2010;107:1124–9.
    OpenUrlAbstract/FREE Full Text
  45. 45.↵
    1. Vajkoczy P,
    2. Knyazev P,
    3. Kunkel A,
    4. Capelle HH,
    5. Behrndt S,
    6. von Tengg-Kobligk H,
    7. et al.
    Dominant-negative inhibition of the Axl receptor tyrosine kinase suppresses brain tumor cell growth and invasion and prolongs survival. Proc Natl Acad Sci U S A 2006;103:5799–804.
    OpenUrlAbstract/FREE Full Text
  46. 46.↵
    1. Mudduluru G,
    2. Allgayer H
    . The human receptor tyrosine kinase Axl gene—promoter characterization and regulation of constitutive expression by Sp1, Sp3 and CpG methylation. Biosci Rep 2008;28:161–76.
    OpenUrlAbstract/FREE Full Text
  47. 47.↵
    1. Breton CV,
    2. Byun HM,
    3. Wenten M,
    4. Pan F,
    5. Yang A,
    6. Gilliland FD
    . Prenatal tobacco smoke exposure affects global and gene-specific DNA methylation. Am J Respir Crit Care Med 2009;180:462–7.
    OpenUrlCrossRefPubMed
  48. 48.↵
    1. Lehmann BD,
    2. Bauer JA,
    3. Chen X,
    4. Sanders ME,
    5. Chakravarthy AB,
    6. Shyr Y,
    7. et al.
    Identification of human triple-negative breast cancer subtypes and preclinical models for selection of targeted therapies. J Clin Invest 2011;121:2750–67.
    OpenUrlCrossRefPubMed
Back to top
Cancer Epidemiology Biomarkers & Prevention: 24 (6)
June 2015
Volume 24, Issue 6

Sign up for alerts

View this article with LENS

Open full page PDF
Article Alerts
Email Article
Citation Tools
Racial Variation in Breast Tumor Promoter Methylation in the Carolina Breast Cancer Study
Kathleen Conway, Sharon N. Edmiston, Chiu-Kit Tse, Christopher Bryant, Pei Fen Kuan, Brionna Y. Hair, Eloise A. Parrish, Ryan May and Theresa Swift-Scanlan
Cancer Epidemiol Biomarkers Prev June 1 2015 (24) (6) 921-930; DOI: 10.1158/1055-9965.EPI-14-1228
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo

Jump to section

Advertisement