DNA methylation profiling distinguishes three clusters of breast cancer cell lines

Methylation change plays an important role in many cellular systems, including cancer development. During recent years, genome-wide or large-scale methylation data has become available thanks to rapid advances in high-throughput biotechnologies. So far, researchers have always used gene expression profiling to study disease subtypes and related therapies. In this study, we investigated methylation profiles in 30 breast cancer cell lines using methylation data generated by microarray technologies. Strong variation of the number of methylation peaks was found among these 30 cell lines; however, more peaks were found in the upstream regions than in downstream regions of genes. We further grouped the methylation profiles of these cell lines into three consensus clusters. Finally, we performed an integrative analysis of breast cancer cell lines using both methylation and gene-expression profiling data. There was no significant correlation between methylation-profiling subtypes and gene-expression profiling subtypes, suggesting the complex nature of methylation in the regulation of gene expression. However, we found basal B cell lines appeared exclusively in two methylation clusters. Although these results are preliminary, this study suggests that methylation profiling might be promising in disease subtype classification and the development of therapeutic strategies.     

DNA methylation profiling to assess pathogenicity of BRCA1 unclassified variants in breast cancer

Germline pathogenic mutations in BRCA1 increase risk of developing breast cancer. Screening for mutations in BRCA1 frequently identifies sequence variants of unknown pathogenicity and recent work has aimed to develop methods for determining pathogenicity. We previously observed that tumor DNA methylation can differentiate BRCA1-mutated from BRCA1-wild type tumors. We hypothesized that we could predict pathogenicity of variants based on DNA methylation profiles of tumors that had arisen in carriers of unclassified variants. We selected 150 FFPE breast tumor DNA samples [47 BRCA1 pathogenic mutation carriers, 65 BRCAx (BRCA1-wild type), 38 BRCA1 test variants] and analyzed a subset (n=54) using the Illumina 450K methylation platform, using the remaining samples for bisulphite pyrosequencing validation. Three validated markers (BACH2, C8orf31, and LOC654342) were combined with sequence bioinformatics in a model to predict pathogenicity of 27 variants (independent test set).  Predictions were compared with standard multifactorial likelihood analysis. Prediction was consistent for c.5194-12G>A (IVS 19-12 G>A) (P>0.99); 13 variants were considered not pathogenic or likely not pathogenic using both approaches. We conclude that tumor DNA methylation data alone has potential to be used in prediction of BRCA1 variant pathogenicity but is not independent of estrogen receptor status and grade, which are used in current multifactorial models to predict pathogenicity.     

CpG island hypermethylation of BRCA1 and loss of pRb as co-occurring events in basal/triple-negative breast cancer

Triple-negative breast cancer (TNBC) occurs in approximately 15% of all breast cancer patients, and the incidence of TNBC is greatly increased in BRCA1 mutation carriers. This study aimed to assess the impact of BRCA1 promoter methylation with respect to breast cancer subtypes in sporadic disease. Tissue microarrays (TMAs) were constructed representing tumors from 303 patients previously screened for BRCA1 germline mutations, of which a subset of 111 sporadic tumors had previously been analyzed with respect to BRCA1 methylation. Additionally, a set of eight tumors from BRCA1 mutation carriers were included on the TMAs. Expression analysis was performed on TMAs by immunohistochemistry (IHC) for BRCA1, pRb, p16, p53, PTEN, ER, PR, HER2, CK5/6, CK8, CK18, EGFR, MUC1, and Ki-67. Data on BRCA1 aberrations and IHC expression was examined with respect to breast cancer-specific survival. The results demonstrate that CpG island hypermethylation of BRCA1 significantly associates with the basal/triple-negative subtype. Low expression of pRb, and high/intense p16, were associated with BRCA1 promoter hypermethylation, and the same effects were seen in BRCA1 mutated tumors. The expression patterns of BRCA1, pRb, p16 and PTEN were highly correlated, and define a subgroup of TNBCs characterized by BRCA1 aberrations, high Ki-67 (≥ 40%) and favorable disease outcome. In conclusion, our findings demonstrate that epigenetic inactivation of the BRCA1 gene associates with RB/p16 dysfunction in promoting TNBCs. The clinical implications relate to the potential use of targeted treatment based on PARP inhibitors in sporadic TNBCs, wherein CpG island hypermethylation of BRCA1 represents a potential marker of therapeutic response.Key words: BRCA1, methylation, epigenetics, triple negative, breast cancer, retionblastoma tumor suppressor gene, pRb, p16     

Molecular Subtype-Specific Expression of MicroRNA-29c in Breast Cancer Is Associated with CpG Dinucleotide Methylation of the Promoter

Basal-like breast cancer is a molecularly distinct subtype of breast cancer that is highly aggressive and has a poor prognosis. MicroRNA-29c (miR-29c) has been shown to be significantly down-regulated in basal-like breast tumors and to be involved in cell invasion and sensitivity to chemotherapy. However, little is known about the genetic and regulatory factors contributing to the altered expression of miR-29c in basal-like breast cancer. We here report that epigenetic modifications at the miR-29c promoter, rather than copy number variation of the gene, may drive the lower expression of miR-29c in basal-like breast cancer. Bisulfite sequencing of CpG sites in the miR-29c promoter region showed higher methylation in basal-like breast cancer cell lines compared to luminal subtype cells with a significant inverse correlation between expression and methylation of miR-29c. Analysis of primary breast tumors using The Cancer Genome Atlas (TCGA) dataset confirmed significantly higher levels of methylation of the promoter in basal-like breast tumors compared to all other subtypes. Furthermore, inhibition of CpG methylation with 5-aza-CdR increases miR-29c expression in basal-like breast cancer cells. Flourescent In Situ Hybridization (FISH) revealed chromosomal abnormalities at miR-29c loci in breast cancer cell lines, but with no correlation between copy number variation and expression of miR-29c. Our data demonstrated that dysregulation of miR-29c in basal-like breast cancer cells may be in part driven by methylation at CpG sites. Epigenetic control of the miR-29c promoter by epigenetic modifiers may provide a potential therapeutic target to overcome the aggressive behavior of these cancers.     

Genetic and Epigenetic Regulation of TOX3 Expression in Breast Cancer

Genome wide association studies (GWAS) have identified low penetrance and high frequency single nucleotide polymorphisms (SNPs) that contribute to genetic susceptibility of breast cancer. The SNPs at 16q12, close to the TOX3 and CASC16 genes, represent one of the susceptibility loci identified by GWAS, showing strong evidence for breast cancer association across various populations. To examine molecular mechanisms of TOX3 regulation in breast cancer, we investigated both genetic and epigenetic factors using cell lines and datasets derived from primary breast tumors available through The Cancer Genome Atlas (TCGA). TOX3 expression is highly up-regulated in luminal subtype tumors compared to normal breast tissues or basal-like tumors. Expression quantitative trait loci (eQTL) analyses revealed significant associations of rs3803662 and rs4784227 genotypes with TOX3 expression in breast tumors. Bisulfite sequencing of four CpG islands in the TOX3 promoter showed a clear difference between luminal and basal-like cancer cell lines. 5-Aza-2’-deoxycytidine treatment of a basal-like cancer cell line increased expression of TOX3. TCGA dataset verified significantly lower levels of methylation of the promoter in luminal breast tumors with an inverse correlation between methylation and expression of TOX3. Methylation QTL (mQTL) analyses showed a weak or no correlation of rs3803662 or rs4784227 with TOX3 promoter methylation in breast tumors, indicating an independent relationship between the genetic and epigenetic events. These data suggest a complex system of TOX3 regulation in breast tumors, driven by germline variants and somatic epigenetic modifications in a subtype specific manner.     

CpG island hypermethylation of BRCA1 and loss of pRb as co-occurring events in basal/triple-negative breast cancer

Triple-negative breast cancer (TNBC) occurs in approximately 15% of all breast cancer patients, and the incidence of TNBC is greatly increased in BRCA1 mutation carriers. This study aimed to assess the impact of BRCA1 promoter methylation with respect to breast cancer subtypes in sporadic disease. Tissue microarrays (TMAs) were constructed representing tumors from 303 patients previously screened for BRCA1 germline mutations, of which a subset of 111 sporadic tumors had previously been analyzed with respect to BRCA1 methylation. Additionally, a set of eight tumors from BRCA1 mutation carriers were included on the TMAs. Expression analysis was performed on TMAs by immunohistochemistry (IHC) for BRCA1, pRb, p16, p53, PTEN, ER, PR, HER2, CK5/6, EGFR, MUC1 and Ki-67. Data on BRCA1 aberrations and IHC expression was examined with respect to breast cancer-specific survival. The results demonstrate that CpG island hypermethylation of BRCA1 significantly associates with the basal/triple-negative subtype. Low expression of pRb, and high/intense p16, were associated with BRCA1 promoter hypermethylation, and the same effects were seen in BRCA1 mutated tumors. The expression patterns of BRCA1, pRb, p16 and PTEN were highly correlated, and define a subgroup of TNBCs characterized by BRCA1 aberrations, high Ki-67 (≥40%) and favorable disease outcome. In conclusion, our findings demonstrate that epigenetic inactivation of the BRCA1 gene associates with RB/p16 dysfunction in promoting TNBCs. The clinical implications relate to the potential use of targeted treatment based on PARP inhibitors in sporadic TNBCs, wherein CpG island hypermethylation of BRCA1 represents a potential marker of therapeutic response.     

Epigenetic and genetic burden measures are associated with tumor characteristics in invasive breast carcinoma

The development and progression of invasive breast cancer is characterized by alterations to the genome and epigenome. However, the relationship between breast tumor characteristics, disease subtypes, and patient outcomes with the cumulative burden of these molecular alterations are not well characterized. We determined the average departure of tumor DNA methylation from adjacent normal breast DNA methylation using Illumina 450K methylation data from 700 invasive breast tumors and 90 adjacent normal breast tissues in The Cancer Genome Atlas. From this we generated a novel summary measure of altered DNA methylation, the DNA methylation dysregulation index (MDI), and examined the relation of MDI with tumor characteristics and summary measures that quantify cumulative burden of genetic mutation and copy number alterations. Our analysis revealed that MDI was significantly associated with tumor stage (P = 0.017). Across invasive breast tumor subtypes we observed significant differences in genome-wide DNA MDIs (P = 4.9E–09) and in a fraction of the genome with copy number alterations (FGA) (P = 4.6E–03). Results from a linear regression adjusted for subject age, tumor stage, and estimated tumor purity indicated a positive significant association of MDI with both MCB and FGA (P = 0.036 and P < 2.2E–16). A recursively partitioned mixture model of all 3 somatic alteration burden measures resulted in classes of tumors whose epigenetic and genetic burden profile were associated with the PAM50 subtype and mutations in TP53, PIK3CA, and CDH1. Together, our work presents a novel framework for characterizing the epigenetic burden and adds to the understanding of the aggregate impact of epigenetic and genetic alterations in breast cancer.     

Deletion of Chromosomes 13q and 14q Is a Common Feature of Tumors with BRCA2 Mutations

Introduction

Germline BRCA1 or BRCA2 mutations account for 20–30% of familial clustering of breast cancer. The main indication for BRCA2 screening is currently the family history but the yield of mutations identified in patients selected this way is low.

Methods

To develop more efficient approaches to screening we have compared the gene expression and genomic profiles of BRCA2-mutant breast tumors with those of breast tumors lacking BRCA1 or BRCA2 mutations.

Results

We identified a group of 66 genes showing differential expression in our training set of 7 BRCA2-mutant tumors and in an independent validation set of 19 BRCA2-mutant tumors. The differentially expressed genes include a prominent cluster of genes from chromosomes 13 and 14 whose expression is reduced. Gene set enrichment analysis confirmed that genes in specific bands on 13q and 14q showed significantly reduced expression, suggesting that the affected bands may be preferentially deleted in BRCA2-mutant tumors. Genomic profiling showed that the BRCA2-mutant tumors indeed harbor deletions on chromosomes 13q and 14q. To exploit this information we have created a simple fluorescence in situ hybridization (FISH) test and shown that it detects tumors with deletions on chromosomes 13q and 14q.

Conclusion

Together with previous reports, this establishes that deletions on chromosomes 13q and 14q are a hallmark of BRCA2-mutant tumors. We propose that FISH to detect these deletions would be an efficient and cost-effective first screening step to identify potential BRCA2-mutation carriers among breast cancer patients without a family history of breast cancer.     

BRCA1 and BRCA2 missense variants of high and low clinical significance influence lymphoblastoid cell line post-irradiation gene expression

The functional consequences of missense variants in disease genes are difficult to predict. We assessed if gene expression profiles could distinguish between BRCA1 or BRCA2 pathogenic truncating and missense mutation carriers and familial breast cancer cases whose disease was not attributable to BRCA1 or BRCA2 mutations (BRCAX cases). 72 cell lines from affected women in high-risk breast ovarian families were assayed after exposure to ionising irradiation, including 23 BRCA1 carriers, 22 BRCA2 carriers, and 27 BRCAX individuals. A subset of 10 BRCAX individuals carried rare BRCA1/2 sequence variants considered to be of low clinical significance (LCS). BRCA1 and BRCA2 mutation carriers had similar expression profiles, with some subclustering of missense mutation carriers. The majority of BRCAX individuals formed a distinct cluster, but BRCAX individuals with LCS variants had expression profiles similar to BRCA1/2 mutation carriers. Gaussian Process Classifier predicted BRCA1, BRCA2 and BRCAX status, with a maximum of 62% accuracy, and prediction accuracy decreased with inclusion of BRCAX samples carrying an LCS variant, and inclusion of pathogenic missense carriers. Similarly, prediction of mutation status with gene lists derived using Support Vector Machines was good for BRCAX samples without an LCS variant (82–94%), poor for BRCAX with an LCS (40–50%), and improved for pathogenic BRCA1/2 mutation carriers when the gene list used for prediction was appropriate to mutation effect being tested (71–100%). This study indicates that mutation effect, and presence of rare variants possibly associated with a low risk of cancer, must be considered in the development of array-based assays of variant pathogenicity.     

Genome-wide CpG island methylation and intergenic demethylation propensities vary among different tumor sites

The epigenetic landscape of cancer includes both focal hypermethylation and broader hypomethylation in a genome-wide manner. By means of a comprehensive genomic analysis on 6637 tissues of 21 tumor types, we here show that the degrees of overall methylation in CpG island (CGI) and demethylation in intergenic regions, defined as ‘backbone’, largely vary among different tumors. Depending on tumor type, both CGI methylation and backbone demethylation are often associated with clinical, epidemiological and biological features such as age, sex, smoking history, anatomic location, histological type and grade, stage, molecular subtype and biological pathways. We found connections between CGI methylation and hypermutability, microsatellite instability, IDH1 mutation, 19p gain and polycomb features, and backbone demethylation with chromosomal instability, NSD1 and TP53 mutations, 5q and 19p loss and long repressive domains. These broad epigenetic patterns add a new dimension to our understanding of tumor biology and its clinical implications.