Abstract
Background: Although obesity is associated with breast cancer incidence and prognosis, the underlying mechanisms are poorly understood. Identification of obesity-associated epigenetic changes in breast tissue may advance mechanistic understanding of breast cancer initiation and progression. The goal of this study, therefore, was to investigate associations between obesity and gene methylation in breast tumors.
Methods: Using the Illumina GoldenGate Cancer I Panel, we estimated the association between body mass index (BMI) and gene methylation in 345 breast tumor samples from phase I of the Carolina Breast Cancer Study, a population-based case–control study. Multivariable linear regression was used to identify sites that were differentially methylated by BMI. Stratification by tumor estrogen receptor (ER) status was also conducted.
Results: In the majority of the 935 probes analyzed (87%), the average beta value increased with obesity (BMI ≥ 30). Obesity was significantly associated with differential methylation (FDR q < 0.05) in just two gene loci in breast tumor tissue overall and in 21 loci among ER-positive tumors. Obesity was associated with methylation of genes that function in immune response, cell growth, and DNA repair.
Conclusions: Obesity is associated with altered methylation overall, and with hypermethylation among ER-positive tumors in particular, suggesting that obesity may influence the methylation of genes with known relevance to cancer. Some of these differences in methylation by obese status may influence levels of gene expression within breast cells.
Impact: If our results are validated, obesity-associated methylation sites could serve as targets for prevention and treatment research. Cancer Epidemiol Biomarkers Prev; 24(3); 580–6. ©2015 AACR.
Introduction
Larger body size, including the body mass index (BMI)-defined categories of overweight and obese, are associated with both the incidence and prognosis of breast cancer, although the risk varies by menopausal status and tumor subtype. In premenopausal women, an inverse association has generally been found between BMI and breast cancer incidence (1, 2), whereas increased BMI has been shown to increase the risk among postmenopausal women (3), with an estimated 3% increase in risk with each one unit gain in BMI (4). There is also evidence that obesity may increase incidence of estrogen receptor (ER)-positive but not ER-negative tumors, although these relationships differ by age and menopausal status (5, 6).
Beyond associations between BMI and breast cancer incidence, obesity is also associated with worse prognosis and outcome in both pre- and postmenopausal women (4). Obese patients with breast cancer are more likely to be diagnosed with advanced stage at diagnosis (3), larger tumor size (1), nodal involvement (1), higher grade (4), and higher mitotic cell count (4). They are also at higher risk of experiencing recurrence (7), metastasis (1), and mortality (7).
The biologic mechanisms underlying obesity's effects on carcinogenesis are poorly understood. Gene methylation, a mechanism for controlling the expression of genes (8), has the potential to influence the carcinogenic process during both tumor formation and progression. In normal tissues, methylation plays a role in whether cells are growing or senescent and how cells differentiate; if these patterns become abnormal, cancer may result (9). Abnormal methylation can also influence tumor invasiveness and metastasis (10).
Although obesity has been associated with increased methylation in several cancer-related genes (11, 12), previous research has been limited by the use of white blood cells (WBC) instead of breast tissue (11, 12) or by focusing on only candidate genes (13). We sought to test whether obesity is associated with methylation in a panel of cancer-related genes in the tumor tissue of women diagnosed with breast cancer. Identifying genes that are differentially methylated by BMI can provide evidence for the underlying biologic processes through which obesity influences breast carcinogenesis and may facilitate the development of chemoprevention and treatment targets.
Materials and Methods
Study participants
Data from phase I of the Carolina Breast Cancer Study (CBCS) was used to assess gene methylation in tumor tissue collected from women with breast cancer. The CBCS study design data collection methods have been previously described (14). Briefly, women diagnosed with breast cancer were recruited with randomized sampling between 1993 and 1996 from the North Carolina Central Cancer Registry. African-American women and women younger than 50 years of age were oversampled to comprise approximately 50% of the study population. Participation in the study was limited to English-speaking women ages 20 to 74 years. Participants were interviewed in-person by trained nurses to collect demographic and breast cancer risk factor data. The University of North Carolina at Chapel Hill Institutional Review Board approved this study.
Tissue specimens
Only participants with sufficient tumor block material for methylation analyses were included in the current study. CBCS pathologists verified the cancer diagnosis of all participants. The breast tumors were grossly dissected after review by a pathologist, who encircled areas of malignant cellularity. While some intratumoral stroma may be present among tumor cells, this sampling method ensured minimal contamination of our methylation signatures by profiles of adipose, stroma, or normal epithelium. As a result, our findings are not likely to be driven by between-sample differences in tumor composition.
Gene methylation assessment
The Illumina Cancer Panel I platform (Illumina Inc.), which consists of 1,505 gene sites with known relevance to cancer (15), measures DNA methylation at cytosine–guanine dinucleotide (CpG) sites. The methylation level of each CpG site on the panel was measured as a beta value, which was calculated by comparing the fluorescent signal ratio of the methylated allele to the sum of the signals from the methylated and unmethylated alleles. The beta values range from 0 to 1, representing the fraction of methylated DNA, with 0 indicating no methylation and 1 indicating complete methylation. Of the 1,505 CpG sites on the platform, 570 were removed from the analysis due to poor performance or because the sites overlapped regions with SNPs or copy number variation, which may render them unreliable (16). The remaining 935 CpG sites (17) were assessed for association with BMI. The input sequence of all probes that were found to be significantly associated with BMI were checked against target sequences in the NCBI Blast database to determine whether any of the probes ambiguously mapped to the genome (18). The array data have been deposited in Gene Expression Omnibus under accession number GSE51557.
In a previous study, the demographic and tumor characteristics were compared for those who were eligible for participation in the study and for those whose methylation values were successfully assessed (17). The only characteristic that reached statistical significance (P < 0.05) was age, where the tumors of younger women were more likely to be successfully assayed. Age was included as a confounder in our analyses. There were no significant differences in tumor size or clinical stage.
Statistical analyses
Linear regression, using limma with an empirical Bayes approach (19), was conducted to assess which CpG sites were associated with prediagnostic BMI. The R statistical package (www.r-project.org/) was used to conduct all regression analyses. Prediagnosis BMI was self-reported at the time of interview and was categorized as normal weight (18.5 ≤ BMI < 25), overweight (25 ≤ BMI < 30), or obese (BMI ≥ 30). Comparisons were 2-category contrasts (obese vs. normal weight, with overweight women excluded) for regression modeling based on previous studies and our own exploratory analyses that showed that the overweight class was highly heterogeneous with regard to obesity-dependent gene expression (20). Therefore, to enhance our ability to detect obesity-associated methylation, the intermediate overweight phenotype group was excluded. Use of a categorical rather than continuous coding for BMI also minimized potential bias due to underestimation of weight (21).
Statistical analyses were restricted to those of African-American or white race. Race and age (treated as a continuous variable) were considered confounders a priori and were included as covariates in all analyses. Other variables assessed for confounding included menopausal status, age at first full-term birth, parity, alcohol use, smoking status, physical activity, fruit and vegetable intake, income, and education status. Among these variables, only menopausal status was associated with both BMI (P < 0.05) and gene methylation (FDR q < 0.05), resulting in a final adjustment set that included age, race, and menopausal status. In addition to testing the association between BMI and gene methylation in breast tumor tissue overall, the relationship was also assessed among ER-positive tumors because prior research has shown that this subgroup of tumors may be more stable genomically than other tumor subtypes (22, 23). In a sensitivity analysis restricted to ER-negative tumors, there were no CpG sites associated with obesity in a regression model adjusting for age, race, and menopausal status.
An FDR q value of <0.05 was used as the statistical significance cutoff value for identifying differentially methylated genes.
Although we did not have data on gene expression for the tissue analyzed in this study, we assessed the correlation between methylation and expression in breast tumors that were collected through The Cancer Genome Atlas (TCGA; http://cancergenome.nih.gov/) for the methylation probes that were significantly associated with BMI in our study. TCGA determined methylation in breast tumors using the Illumina HumanMethylation 450 Beadchip panel (24), whereas the present study used the Illumina GoldenGate Cancer I Panel. Only 7 of the BMI-associated probes in the present study matched directly with probes in the 450 Beadchip panel. For the remaining BMI-associated probes, we identified probes on the 450 Beadchip array that were within 200 base pairs, either up- or downstream, of the GoldenGate probe. The correlation between methylation and expression was examined for all of the breast tumor tissue available in the TCGA and also for hormone receptor–positive tumors (the closest approximation to ER-positive tumors for which we had data access).
Results
Methylation sites associated with BMI
Of the 517 patients eligible for study participation, 492 had both prediagnosis BMI data and methylation data; of these, 345 were either obese or of normal weight (Table 1). The median age for the 492 study participants was 47 years. In unadjusted analyses, 30 CpG sites were differentially methylated by BMI status (obese vs. normal) at q < 0.05; after adjustment for age, race, and menopausal status, 2 CpG sites remained significantly differentially methylated: SH3BP2_P771_R on the promoter region of the SH3BP2 gene, which regulates transcriptional activity in immune cells, and Xist_seq_80_S95_R on the non–protein-coding XIST gene. For both sites, the average methylation beta values were higher in the tissue of obese patients, when compared with those of normal weight.
Selected characteristics of eligible breast cancer cases in the CBCS, phase I
. | Totala . | Normalb . | Overweightb . | Obeseb . |
---|---|---|---|---|
. | n (%) . | n (%) . | n (%) . | n (%) . |
Total | 492 (100) | 195 (100) | 147 (100) | 150 (100) |
Race | ||||
African-American | 200 (40.7) | 37 (19.0) | 61 (41.5) | 102 (68.0) |
White | 282 (57.3) | 152 (77.9) | 83 (56.5) | 47 (31.3) |
Other | 10 (2.0) | 6 (3.1) | 3 (2.0) | 1 (0.7) |
Menopausal status | ||||
Premenopausal | 262 (53.3) | 108 (55.4) | 78 (53.1) | 76 (50.7) |
Postmenopausal | 230 (46.7) | 87 (44.6) | 69 (46.9) | 74 (49.3) |
Alcohol use | ||||
Never | 141 (28.7) | 41 (21.0) | 45 (30.6) | 55 (36.7) |
Ever | 351 (71.3) | 154 (79.0) | 102 (69.4) | 95 (63.3) |
Smoking status | ||||
Never | 251 (51.0) | 93 (47.7) | 79 (53.7) | 79 (52.7) |
Former | 146 (29.7) | 57 (29.2) | 42 (28.6) | 47 (31.3) |
Current | 95 (19.3) | 45 (23.1) | 26 (17.7) | 24 (16.0) |
Income, $ | ||||
<15,000 | 119 (24.2) | 31 (17.2) | 30 (22.6) | 58 (40.3) |
15,000–30,000 | 99 (20.1) | 38 (21.1) | 27 (20.3) | 34 (23.6) |
30,000–50,000 | 113 (23.0) | 50 (27.8) | 33 (24.8) | 30 (20.8) |
≥50,000 | 126 (25.6) | 61 (33.9) | 43 (32.3) | 22 (15.3) |
Educationc | ||||
<High school | 90 (18.3) | 18 (9.2) | 26 (17.7) | 46 (30.7) |
High school | 269 (54.7) | 107 (54.9) | 87 (59.2) | 75 (50.0) |
College | 133 (27.0) | 70 (35.9) | 34 (23.1) | 29 (19.3) |
ER-positive | ||||
Yes | 285 (59.6) | 118 (62.1) | 77 (54.2) | 90 (61.6) |
No | 193 (40.4) | 72 (37.9) | 65 (45.8) | 56 (38.4) |
Staged | ||||
I | 168 (36.7) | 73 (38.8) | 49 (36.8) | 46 (33.6) |
II | 238 (52.0) | 99 (52.7) | 70 (52.6) | 69 (50.4) |
III and IV | 52 (11.4) | 16 (8.5) | 14 (10.5) | 22 (16.1) |
. | Totala . | Normalb . | Overweightb . | Obeseb . |
---|---|---|---|---|
. | n (%) . | n (%) . | n (%) . | n (%) . |
Total | 492 (100) | 195 (100) | 147 (100) | 150 (100) |
Race | ||||
African-American | 200 (40.7) | 37 (19.0) | 61 (41.5) | 102 (68.0) |
White | 282 (57.3) | 152 (77.9) | 83 (56.5) | 47 (31.3) |
Other | 10 (2.0) | 6 (3.1) | 3 (2.0) | 1 (0.7) |
Menopausal status | ||||
Premenopausal | 262 (53.3) | 108 (55.4) | 78 (53.1) | 76 (50.7) |
Postmenopausal | 230 (46.7) | 87 (44.6) | 69 (46.9) | 74 (49.3) |
Alcohol use | ||||
Never | 141 (28.7) | 41 (21.0) | 45 (30.6) | 55 (36.7) |
Ever | 351 (71.3) | 154 (79.0) | 102 (69.4) | 95 (63.3) |
Smoking status | ||||
Never | 251 (51.0) | 93 (47.7) | 79 (53.7) | 79 (52.7) |
Former | 146 (29.7) | 57 (29.2) | 42 (28.6) | 47 (31.3) |
Current | 95 (19.3) | 45 (23.1) | 26 (17.7) | 24 (16.0) |
Income, $ | ||||
<15,000 | 119 (24.2) | 31 (17.2) | 30 (22.6) | 58 (40.3) |
15,000–30,000 | 99 (20.1) | 38 (21.1) | 27 (20.3) | 34 (23.6) |
30,000–50,000 | 113 (23.0) | 50 (27.8) | 33 (24.8) | 30 (20.8) |
≥50,000 | 126 (25.6) | 61 (33.9) | 43 (32.3) | 22 (15.3) |
Educationc | ||||
<High school | 90 (18.3) | 18 (9.2) | 26 (17.7) | 46 (30.7) |
High school | 269 (54.7) | 107 (54.9) | 87 (59.2) | 75 (50.0) |
College | 133 (27.0) | 70 (35.9) | 34 (23.1) | 29 (19.3) |
ER-positive | ||||
Yes | 285 (59.6) | 118 (62.1) | 77 (54.2) | 90 (61.6) |
No | 193 (40.4) | 72 (37.9) | 65 (45.8) | 56 (38.4) |
Staged | ||||
I | 168 (36.7) | 73 (38.8) | 49 (36.8) | 46 (33.6) |
II | 238 (52.0) | 99 (52.7) | 70 (52.6) | 69 (50.4) |
III and IV | 52 (11.4) | 16 (8.5) | 14 (10.5) | 22 (16.1) |
aColumns may not sum to total due to missing data.
bNormal weight: 18.5 ≤ BMI < 25; overweight: 25 ≤ BMI < 30; obese: BMI ≥ 30.
cHigh school: high school completion (or equivalent), technical or business school, and some college; college: college completion and post-graduate school or professional degree.
dAccording to AJCC breast tumor staging guidelines.
Methylation in ER-positive tumors
In analyses restricted to ER-positive tumors (n = 208), 21 probes were significantly associated with BMI (obese vs. normal) after adjustment for age, race, and menopausal status (Table 2). Twelve of the 21 significant sites were located in CpG islands. Three CpG sites mapped to chromosome 4 (SH3BP2_P771_R, SH3BP2_E18_F, and PKD2_P336_R), 3 mapped to the X chromosome (Xist_seq_80_S95_R, MCF2_P1024_R, and BTK_P105_F), and 2 each mapped to chromosome 10 (ERCC6_P698_R, BMPR1A_E88_F), chromosome 11 (TSG101_P257_R, SPI1_E205_F), and chromosome 12 (IGFBP6_E47_F, ARHGDIB_P148_R). Two of the probes (BMPR1A_E88_F and TDGF1_P428_R) mapped to more than one chromosome; as such, these probes should be interpreted with caution. Figure 1 displays a volcano plot, a visual representation of the regression analysis wherein each circle represents a methylation site; circles with positive coefficients signify higher average beta values in the obese compared with those of normal weight, whereas circles with negative coefficients signify lower beta values in the obese. The majority of the sites were hypermethylated in the obese group, as were all of the probes that were significantly associated with obesity (circles above the horizontal line).
Volcano plot of 935 methylation sites in ER-positive tumors. The x-axis represents coefficient from model regressing obesity on CpG site, adjusting for age, race, and menopausal status and y-axis represents −log10 of P value from model regressing obesity on CpG site, adjusting for age, race, and menopausal status. Horizontal line indicates FDR q-value cutoff of 0.05. Each circle represents a CpG site.
Volcano plot of 935 methylation sites in ER-positive tumors. The x-axis represents coefficient from model regressing obesity on CpG site, adjusting for age, race, and menopausal status and y-axis represents −log10 of P value from model regressing obesity on CpG site, adjusting for age, race, and menopausal status. Horizontal line indicates FDR q-value cutoff of 0.05. Each circle represents a CpG site.
Methylation sites associated with BMI (FDR < 0.05) in the CBCS among patients with ER-positive tumors
Probe ID . | Gene . | Gene description . | CHRa . | Islandb . | Coefficient . | P . | q . | Normalc . | Overweightc,d . | Obesec . |
---|---|---|---|---|---|---|---|---|---|---|
DNMT3B_P352_R | DNMT3B | DNA cytosine-5 methyltransferase 3 beta isoform 2 | 20q11 | N | 0.0230 | 0.00017 | 0.03818 | 15.5% | 17.3% | 17.6% |
IGFBP6_E47_F | IGFBP6 | Insulin-like growth factor binding protein 6 | 12q13 | N | 0.0447 | 0.00032 | 0.03818 | 37.5% | 40.5% | 43.3% |
ERCC6_P698_R | ERCC6 | Excision repair cross-complementing rodent repair deficiency, complementation group 6 | 10q11 | Y | 0.0480 | 0.00037 | 0.03818 | 43.1% | 46.8% | 47.3% |
SH3BP2_P771_Re | SH3BP2 | SH3 domain–binding protein 2 | 4p16 | Y | 0.0324 | 0.00048 | 0.03818 | 17.4% | 20.8% | 22.7% |
PLAGL1_E68_R | PLAGL1 | Pleiomorphic adenoma gene-like 1 isoform 1 | 6q24–25 | Y | 0.0422 | 0.00049 | 0.03818 | 69.8% | 72.5% | 75.4% |
FZD9_E458_F | FZD9 | Frizzled 9 | 7q11 | Y | 0.0535 | 0.00096 | 0.04640 | 49.2% | 50.0% | 55.4% |
BMPR1A_E88_F | BMPR1A | Bone morphogenetic protein receptor, type IA precursor | 10q22 | Y | 0.0320 | 0.00101 | 0.04640 | 65.0% | 66.7% | 69.2% |
TSG101_P257_R | TSG101 | Tumor susceptibility gene 101 | 11p15 | Y | 0.0400 | 0.00132 | 0.04640 | 25.4% | 29.0% | 29.8% |
RIPK4_E166_F | RIPK4 | Receptor-interacting serine-threonine kinase 4 | 21q22 | Y | 0.0171 | 0.00133 | 0.04640 | 7.2% | 9.6% | 10.1% |
Xist_seq_80_S95_Re | XIST | X inactive–specific transcript (non–protein-coding) | Xq13 | N | 0.0299 | 0.00143 | 0.04640 | 74.3% | 76.3% | 78.9% |
TDGF1_P428_R | TDGF1 | Teratocarcinoma-derived growth factor 1 | 3p21 | N | 0.0412 | 0.00147 | 0.04640 | 40.7% | 41.8% | 44.6% |
MCF2_P1024_R | MCF2 | MCF.2 cell line–derived transforming sequence | Xq27 | N | 0.0393 | 0.00151 | 0.04640 | 77.6% | 80.6% | 82.3% |
SPI1_E205_F | SPI1 | Spleen focus-forming virus (SFFV) proviral integration oncogene spi1 | 11p11 | Y | 0.0372 | 0.00154 | 0.04640 | 54.1% | 57.8% | 58.8% |
SH3BP2_E18_F | SH3BP2 | SH3 domain–binding protein 2 | 4p16 | N | 0.0284 | 0.00208 | 0.04934 | 25.3% | 28.6% | 29.9% |
ABO_E110_F | ABO | ABO blood group (transferase A, alpha 1-3-N-acetylgalactosaminyltransferase; transferase B, alpha 1–3-galactosyltransferase) | 9q34 | Y | 0.0373 | 0.00225 | 0.04934 | 15.8% | 18.0% | 20.6% |
CHI3L2_P226_F | CHI3L2 | Chitinase 3-like 2 isoform c | 1p13 | N | 0.0352 | 0.00227 | 0.04934 | 74.8% | 75.2% | 74.6% |
BTK_P105_F | BTK | Bruton agammaglobulinemia tyrosine kinase | Xq21–22 | N | 0.0363 | 0.00240 | 0.04934 | 39.6% | 42.4% | 44.8% |
SNURF_P78_F | SNURF | SNRPN upstream reading frame protein | 15q11–12 | Y | 0.0381 | 0.00254 | 0.04934 | 69.5% | 71.0% | 74.3% |
ARHGDIB_P148_R | ARHGDIB | Rho GDP dissociation inhibitor (GDI) beta | 12p12 | N | 0.0340 | 0.00255 | 0.04934 | 44.5% | 45.5% | 47.0% |
SEZ6L_P299_F | SEZ6L | Seizure-related 6 homolog (mouse)-like precursor | 22q12 | Y | 0.0482 | 0.00255 | 0.04934 | 17.5% | 18.0% | 21.3% |
PKD2_P336_R | PKD2 | Polycystin 2 | 4q22 | Y | 0.0260 | 0.00265 | 0.04934 | 13.1% | 14.6% | 16.5% |
Probe ID . | Gene . | Gene description . | CHRa . | Islandb . | Coefficient . | P . | q . | Normalc . | Overweightc,d . | Obesec . |
---|---|---|---|---|---|---|---|---|---|---|
DNMT3B_P352_R | DNMT3B | DNA cytosine-5 methyltransferase 3 beta isoform 2 | 20q11 | N | 0.0230 | 0.00017 | 0.03818 | 15.5% | 17.3% | 17.6% |
IGFBP6_E47_F | IGFBP6 | Insulin-like growth factor binding protein 6 | 12q13 | N | 0.0447 | 0.00032 | 0.03818 | 37.5% | 40.5% | 43.3% |
ERCC6_P698_R | ERCC6 | Excision repair cross-complementing rodent repair deficiency, complementation group 6 | 10q11 | Y | 0.0480 | 0.00037 | 0.03818 | 43.1% | 46.8% | 47.3% |
SH3BP2_P771_Re | SH3BP2 | SH3 domain–binding protein 2 | 4p16 | Y | 0.0324 | 0.00048 | 0.03818 | 17.4% | 20.8% | 22.7% |
PLAGL1_E68_R | PLAGL1 | Pleiomorphic adenoma gene-like 1 isoform 1 | 6q24–25 | Y | 0.0422 | 0.00049 | 0.03818 | 69.8% | 72.5% | 75.4% |
FZD9_E458_F | FZD9 | Frizzled 9 | 7q11 | Y | 0.0535 | 0.00096 | 0.04640 | 49.2% | 50.0% | 55.4% |
BMPR1A_E88_F | BMPR1A | Bone morphogenetic protein receptor, type IA precursor | 10q22 | Y | 0.0320 | 0.00101 | 0.04640 | 65.0% | 66.7% | 69.2% |
TSG101_P257_R | TSG101 | Tumor susceptibility gene 101 | 11p15 | Y | 0.0400 | 0.00132 | 0.04640 | 25.4% | 29.0% | 29.8% |
RIPK4_E166_F | RIPK4 | Receptor-interacting serine-threonine kinase 4 | 21q22 | Y | 0.0171 | 0.00133 | 0.04640 | 7.2% | 9.6% | 10.1% |
Xist_seq_80_S95_Re | XIST | X inactive–specific transcript (non–protein-coding) | Xq13 | N | 0.0299 | 0.00143 | 0.04640 | 74.3% | 76.3% | 78.9% |
TDGF1_P428_R | TDGF1 | Teratocarcinoma-derived growth factor 1 | 3p21 | N | 0.0412 | 0.00147 | 0.04640 | 40.7% | 41.8% | 44.6% |
MCF2_P1024_R | MCF2 | MCF.2 cell line–derived transforming sequence | Xq27 | N | 0.0393 | 0.00151 | 0.04640 | 77.6% | 80.6% | 82.3% |
SPI1_E205_F | SPI1 | Spleen focus-forming virus (SFFV) proviral integration oncogene spi1 | 11p11 | Y | 0.0372 | 0.00154 | 0.04640 | 54.1% | 57.8% | 58.8% |
SH3BP2_E18_F | SH3BP2 | SH3 domain–binding protein 2 | 4p16 | N | 0.0284 | 0.00208 | 0.04934 | 25.3% | 28.6% | 29.9% |
ABO_E110_F | ABO | ABO blood group (transferase A, alpha 1-3-N-acetylgalactosaminyltransferase; transferase B, alpha 1–3-galactosyltransferase) | 9q34 | Y | 0.0373 | 0.00225 | 0.04934 | 15.8% | 18.0% | 20.6% |
CHI3L2_P226_F | CHI3L2 | Chitinase 3-like 2 isoform c | 1p13 | N | 0.0352 | 0.00227 | 0.04934 | 74.8% | 75.2% | 74.6% |
BTK_P105_F | BTK | Bruton agammaglobulinemia tyrosine kinase | Xq21–22 | N | 0.0363 | 0.00240 | 0.04934 | 39.6% | 42.4% | 44.8% |
SNURF_P78_F | SNURF | SNRPN upstream reading frame protein | 15q11–12 | Y | 0.0381 | 0.00254 | 0.04934 | 69.5% | 71.0% | 74.3% |
ARHGDIB_P148_R | ARHGDIB | Rho GDP dissociation inhibitor (GDI) beta | 12p12 | N | 0.0340 | 0.00255 | 0.04934 | 44.5% | 45.5% | 47.0% |
SEZ6L_P299_F | SEZ6L | Seizure-related 6 homolog (mouse)-like precursor | 22q12 | Y | 0.0482 | 0.00255 | 0.04934 | 17.5% | 18.0% | 21.3% |
PKD2_P336_R | PKD2 | Polycystin 2 | 4q22 | Y | 0.0260 | 0.00265 | 0.04934 | 13.1% | 14.6% | 16.5% |
aCHR = chromsome.
bIsland indicates whether the probe is located in a CpG Island; N = no; Y = yes.
cUnadjusted % methylation; normal weight: 18.5 ≤ BMI < 25; overweight: 25 ≤ BMI < 30; obese BMI ≥ 30.
dOverweight individuals were not included in the regression analysis, but are included here for descriptive purposes.
eAssociated with BMI in overall breast tumor tissue regression analysis.
The methylation sites that showed the most significant changes in the beta values between obese versus normal BMI patients were on the DNMT3B gene, which codes for a DNA methyltransferase; the IGFBP6 gene, which codes for an IGF-binding protein; the ERCC6 gene, which codes for a protein that functions in excision repair; and the SH3BP2 gene.
Correlation with gene expression
The 21 methylation sites that were associated with BMI in the overall or ER-positive analyses were mapped to breast tumors collected by TCGA to assess the correlation between gene methylation and expression. Of the 21 sites, 20 matched to probes that were at least within 200 base pairs of the probe used in TCGA data and 3 sites matched to more than one probe (Supplementary Table S1). Approximately 60% of the 20 BMI-associated methylation sites were significantly and inversely associated with gene expression (correlation coefficient ≤ −0.20) in the full TCGA breast tumor tissue data set, and 55% were significantly and inversely associated with expression in the hormone receptor–positive tumors. Although probes that were correlated with gene expression in both the full TCGA dataset and among hormone-positive tumors were more likely to be on CpG islands than probes that were not associated with gene expression, the associations were not significant at an α of 0.05 level (P = 0.36 and 0.17, respectively).
Discussion
Because methylation is a mechanism for controlling gene expression (8), exploring the relationship between obesity and methylation can provide insight on BMI-associated expression of cancer-related genes. Our regression analyses suggest that BMI is associated with gene methylation in cancer-related genes overall and especially among ER-positive breast tumors. When assessing the correlation between methylation and expression for our BMI-associated methylation sites (from both our overall and ER-positive analyses) in breast tumors collected through TCGA, about 60% showed a significant inverse correlation between methylation and expression.
Previous research in breast tissue has demonstrated that BMI is associated with the differential expression of genes involved in pathways related to immune response and IGF availability (20, 25). In our study, BMI was associated with methylation of genes involved not only in immune response and IGF pathways but also in DNA methylation and DNA repair. Notably, SH3BP2 was the most significant BMI-associated site in our overall analysis, whereas SH3BP2, IGFBP6, DNMT3, and ERCC6 were the most statistically significant BMI-associated sites among ER-positive tumors. All 4 of these genes were hypermethylated with obesity. The SH3BP2 gene codes for a protein that acts in the signal transduction of various immune response pathways (26). It has been implicated in responses involving neutrophils (27), B cells (28–30), and T cells (26), and it may function as a tumor suppressor gene (26). Our finding that SH3BP2 is hypermethylated in the tumor tissue of obese patients with breast cancer suggests that obesity may confer an abnormal immune response in breast tissue. This is interesting in light of previous research also showing increased infiltration of macrophages in the breast tissue of obese women (20, 31). The IGFBP6 gene codes for a protein that binds with high affinity to IGFII (32), which plays a role in cellular proliferation, differentiation, and migration, and is regulated by the IGFBP6 protein (32). Expression of IGFBP6 may be associated with cancer incidence (32), with studies finding lower expression of the protein in the serum of patients with breast cancer than in the serum of those with benign breast disease (33, 34). Our finding of hypermethylation of IGFBP6 in obese patients suggests that obesity may lead to reduced expression in breast tissue via gene methylation. BMI may also impact broader mechanisms controlling methylation. For example, DNMT3 is a member of a family of DNA methyltransferases that code for enzymes involved in the methylation of promoter regions of genes (35), and DNMT3 mRNA has been shown to be overexpressed in breast tumor tissue (35, 36) and cell lines (37). Finally, the ERCC6 protein, also known as the Cockayne syndrome complementation group B (CSB) protein, is involved in transcription-coupled nucleotide excision DNA repair (38).
Our study was strengthened by examining the association between BMI and methylation in breast tumor tissue rather than WBCs, as in most previous studies. Nonetheless, prior research has shown that differences in methylation between those with and without cancer may be driven by changes in the distribution of WBC populations (39). Obesity may also be associated with changes in the distribution of WBC populations, resulting in methylation differences by BMI status. In one study of WBCs, methylation was increased in those with higher BMI for 3 of 4 differentially methylated CpG sites (MMP9, PM20D1, and PRKG1; ref. 12). Of those sites, only MMP9 was included in the methylation panel used in this study, and it was not associated with obesity either overall or among ER-positive tumors. In another study of methylation in WBCs, BMI was associated with increased methylation of several sites on HIF3A, a gene that encodes for a protein involved in responses to hypoxia (11). However, in the latter study, the magnitude of change was not replicated when the HIF3A loci were examined in adipose tissue, underscoring the importance of tissue specific methylation. There were no CpG sites on the HIF3A gene in the methylation panel used for this study.
Given previous reports of differentially methylated genes associated with BMI in both genome-wide and candidate loci studies (11, 12), it is interesting that studies examining global methylation in WBCs are inconclusive regarding the direction of methylation in response to increasing BMI. In a study of global methylation in 85 women, overweight and obese patients had nonstatistically significant higher levels of methylation than those with a normal BMI after adjustment for variables including childhood smoking exposure, nulliparity, and age at first birth (40). When methylation of the repetitive element, LINE-1, was measured in the peripheral blood of men and women, those with a normal BMI had the lowest mean methylation percentage and those who were obese had the highest methylation percentages, although the results were not statistically significant in crude analyses (41). Other studies reviewed in Terry and colleagues (42) did not find significant correlations between methylation values in the LINE-1 and ALU repetitive elements and BMI (43, 44) and one study examining BMI in women of childbearing age reported less LINE-1 methylation in those with higher BMI (45).
These aforementioned studies suggest that obesity may be weakly associated with an increase in methylation in WBCs but do not address whether BMI may affect methylation in breast tissue. Given our emphasis on the tumor tissue, the most relevant previous research for our purposes are studies of breast tumor methylation and BMI (13, 46). In one study, Tao and colleagues examined the association between obesity and methylation in 3 candidate genes, E-cadherin, p16 (CDKN2A), and RAR-β(2), in breast tumor tissue and showed no statistically significant differences in percentage methylation by BMI status in their study population; regression analyses stratified by ER status were not reported (13). Although the 3 genes examined by Tao and colleagues (13) were represented on our methylation panel, we were unable to ascertain the exact CpG loci on which methylation values were measured. Given that DNA methylation levels can vary greatly within short base pair distances (47), we were unable to meaningfully compare our gene methylation results with those of Tao and colleagues. In another study, Naudshad and colleagues examined the association between BMI and methylation in Ec-SOD (SOD3), RASSF1, BRCA1, and BNIP3; all of the gene loci except for BNIP3 were significantly, increasingly methylated with higher BMI, whereas BNIP3 was significantly, inversely associated with BMI (46). All but the BNIP3 gene examined by Naushad and colleagues (46) were represented on our methylation panel. However, there were no exact loci matches between the probes available on our methylation panel and the Ec-SOD and BRCA1 sites examined by Naushad and colleagues; we were unable to determine the location of the RASSF1 probe. Because gene methylation is loci-specific (47), we were unable to compare their reported methylation results to our findings.
Our study is novel in that it assesses the relationship between BMI and methylation in breast tumor tissue using a panel of cancer-associated genes. This enabled us to discern relationships between obesity and methylation that might not be apparent in WBCs or that might be overlooked by focusing on a small number of candidate genes. Nonetheless, because the methylation array was limited to cancer-relevant genes, only a subset of genes that may be differentially methylated by obese status were analyzed. Thus, our results are biased toward methylation sites with established associations with cancer. Furthermore, there is evidence that genes with known relevance to cancer can be either hypermethylated or hypomethylated relative to normal breast tissue (10, 48). However, the goal of our study was to determine the association between obesity and methylation status, regardless of the genes' hyper- or hypomethylation status relative to normal tissue. Our analysis demonstrates that obesity is associated with the degree of methylation in some cancer-relevant genes, providing insight into which genes should be targeted by research for further exploration.
Although our study was strengthened by assessing the correlation between methylation and expression for our BMI-associated sites, we were unable to assess whether these differences in expression correlated with functional differences. In addition, we were unable to validate our results in an independent dataset because we were not able to identify public data sources where both methylation and BMI were measured and reported. TCGA, for example, does not collect data on BMI. Our study, however, contributes data on both BMI and methylation status to the public domain, and future studies can seek to replicate our findings using these data.
Our analysis was also strengthened by analyzing ER-positive tumors, which have been shown to feature fewer chromosomal aberrations than other tumor subtypes (22, 23). However, there were some limitations associated with our analysis. Although tissue from ER-positive tumors may be more stable than that from other tumor subtypes, in the analysis of any tumor tissue, it is difficult to determine which methylation changes influence the carcinogenic process and which occur as a result of it (49). Furthermore, there may be heterogeneity within ER-positive tumors (i.e., some tumors are luminal A and others luminal B) that may have influenced our power to detect associations between methylation and obesity. A further limitation of the study design was the lack of power to adequately assess the association between BMI and gene methylation by menopausal status. Given the variability in the obesity–breast cancer association by menopausal status (1, 2), it is possible that different or additional genes would be associated with methylation by menopausal status. Finally, the motivating goal of this study was to explore whether there might be epigenetic biomarkers of the effect of obesity on epithelial tumor cells that can be used to direct research of targeted breast cancer therapies. However, a focus on normal breast tissue is also important for future research, as it can give insight into how obesity may influence cancer initiation and promotion.
In summary, our study provides evidence that obesity may be associated with increased gene methylation in breast tumor tissue and may play a role in the expression of genes associated with immune response, IGF availability, and other pathways involved in carcinogenesis. If obesity does act through methylation to control expression of carcinogenesis-related genes, differentially methylated sites could serve as targets for the chemoprevention and treatment of breast cancer in women with obesity as a risk factor.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Authors' Contributions
Conception and design: B.Y. Hair, M.A. Troester, K. Conway
Development of methodology: B.Y. Hair, S.N. Edmiston, K. Conway
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): B.Y. Hair, S.N. Edmiston, K. Conway
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): B.Y. Hair, M.A. Troester, M.C. Wu, A.F. Olshan, T. Swift-Scanlan, K. Conway
Writing, review, and/or revision of the manuscript: B.Y. Hair, M.A. Troester, S.N. Edmiston, W.R. Robinson, A.F. Olshan, T. Swift-Scanlan, K. Conway
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): B.Y. Hair, S.N. Edmiston, E.A. Parrish, K. Conway
Study supervision: M.A. Troester, K. Conway
Grant Support
This research was funded, in part, by the University Cancer Research Fund of North Carolina, the National Cancer Institute Specialized Program of Research Excellence (SPORE) in Breast Cancer (NIH/NCI P50-CA58223), the Lineberger Comprehensive Cancer Center Core Grant (NIH/NCI P30-CA16086). Support for K. Conway was provided by a grant from Susan G. Komen for the Cure (#KG081397). Support for W.R. Robinson was provided by the National Cancer Institute (1K01CA172717-01). Support for M.A. Troester was provided by the National Cancer Institute/National Institute of Environmental Health Sciences Breast Cancer and the Environment Research Program (U01-ES019472).
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.