Summary We analyzed primary breast cancers by genomic DNA copy number arrays, DNA methylation, exome sequencing, mRNA arrays, microRNA sequencing and reverse phase protein arrays. Our ability to integrate information across platforms provided key insights into previously-defined gene expression subtypes and demonstrated the existence of four main breast cancer classes when combining data from five platforms, each of which shows significant molecular heterogeneity. Somatic mutations in only three genes (TP53, PIK3CA and GATA3) occurred at > 10% incidence across all breast cancers; however, there were numerous subtype-associated and novel gene mutations including the enrichment of specific mutations in GATA3, PIK3CA and MAP3K1 with the Luminal A subtype. We identified two novel protein expression-defined subgroups, possibly contributed by stromal/microenvironmental elements, and integrated analyses identified specific signaling pathways dominant in each molecular subtype including a HER2/p-HER2/HER1/p-HER1 signature within the HER2-Enriched expression subtype. Comparison of Basal-like breast tumors with high-grade Serous Ovarian tumors showed many molecular commonalities, suggesting a related etiology and similar therapeutic opportunities. The biologic finding of the four main breast cancer subtypes caused by different subsets of genetic and epigenetic abnormalities raises the hypothesis that much of the clinically observable plasticity and heterogeneity occurs within, and not across, these major biologic subtypes of breast cancer.

Colorectal cancer (CRC) is a heterogeneous disease in which unique subtypes are characterized by distinct genetic and epigenetic alterations. Here we performed comprehensive genome-scale DNA methylation profiling of 125 colorectal tumors and 29 adjacent normal tissues. We identified four DNA methylation-based subgroups of CRC using model-based cluster analyses. Each subtype shows characteristic genetic and clinical features, indicating that they represent biologically distinct subgroups. A CIMP-high (CIMP-H) subgroup, which exhibits an exceptionally high frequency of cancer-specific DNA hypermethylation, is strongly associated with MLH1 DNA hypermethylation and the BRAF V600E mutation. A CIMPlow (CIMP-L) subgroup is enriched for KRAS mutations and characterized by DNA hypermethylation of a subset of CIMP-H-associated markers rather than a unique group of CpG islands. Non-CIMP tumors are separated into two distinct clusters. One non-CIMP subgroup is distinguished by a significantly higher frequency of TP53 mutations and frequent occurrence in the distal colon, while the tumors that belong to the fourth group exhibit a low frequency of both cancer-specific DNA hypermethylation and gene mutations and are significantly enriched for rectal tumors. Furthermore, we identified 112 genes that were down-regulated more than twofold in CIMP-H tumors together with promoter DNA hypermethylation. These represent~7% of genes that acquired promoter DNA methylation in CIMP-H tumors. Intriguingly, 48/112 genes were also transcriptionally down-regulated in non-CIMP subgroups, but this was not attributable to promoter DNA hypermethylation. Together, we identified four distinct DNA methylation subgroups of CRC and provided novel insight regarding the role of CIMP-specific DNA hypermethylation in gene silencing.

Chromatin alterations are fundamental hallmarks of cancer. To study chromatin alterations in primary gastric adenocarcinomas, we perform nanoscale chromatin immunoprecipitation sequencing of multiple histone modifications in five gastric cancers and matched normal tissues. We identify hundreds of somatically altered promoters and predicted enhancers. Many cancer-associated promoters localize to genomic sites lacking previously annotated transcription start sites (cryptic promoters), driving expression of nearby genes involved in gastrointestinal cancer, embryonic development and tissue specification. Cancer-associated promoters overlap with embryonic stem cell regions targeted by polycomb repressive complex 2, exhibiting promoter bivalency and DNA methylation loss. We identify somatically acquired elements exhibiting germline allelic biases and non-coding somatic mutations creating new promoters. Our findings demonstrate the feasibility of profiling chromatin from solid tumours with limited tissue to identify regulatory elements, transcriptional patterns and regulatory genetic variants associated with cancer.

Targeted therapies against the BCR-ABL1 kinase have revolutionized treatment of chronic phase (CP) chronic myeloid leukemia (CML). In contrast, management of blast crisis (BC) CML remains challenging because BC cells acquire complex molecular alterations that confer stemness features to progenitor populations and resistance to BCR-ABL1 tyrosine kinase inhibitors. Comprehensive models of BC transformation have proved elusive because of the rarity and genetic heterogeneity of BC, but are important for developing biomarkers predicting BC progression and effective therapies. To better understand BC, we performed an integrated multiomics analysis of 74 CP and BC samples using whole-genome and exome sequencing, transcriptome and methylome profiling, and chromatin immunoprecipitation followed by high-throughput sequencing. Employing pathway-based analysis, we found the BC genome was significantly enriched for mutations affecting components of the polycomb repressive complex (PRC) pathway. While transcriptomically, BC progenitors were enriched and depleted for PRC1- and PRC2-related gene sets respectively. By integrating our data sets, we determined that BC progenitors undergo PRC-driven epigenetic reprogramming toward a convergent transcriptomic state. Specifically, PRC2 directs BC DNA hypermethylation, which in turn silences key genes involved in myeloid differentiation and tumor suppressor function via so-called epigenetic switching, whereas PRC1 represses an overlapping and distinct set of genes, including novel BC tumor suppressors. On the basis of these observations, we developed an integrated model of BC that facilitated the identification of combinatorial therapies capable of reversing BC reprogramming (decitabine+PRC1 inhibitors), novel PRC-silenced tumor suppressor genes (NR4A2), and gene expression signatures predictive of disease progression and drug resistance in CP.

Activation of estrogen receptor ␣ (ER␣) results in both induction and repression of gene transcription; while mechanistic details of estrogen induction are well described, details of repression remain largely unknown. We characterized several ER␣-repressed targets and examined in detail the mechanism for estrogen repression of Reprimo (RPRM), a cell cycle inhibitor. Estrogen repression of RPRM is rapid and robust and requires a tripartite interaction between ER␣, histone deacetylase 7 (HDAC7), and FoxA1. HDAC7 is the critical HDAC needed for repression of RPRM; it can bind to ER␣ and represses ER␣'s transcriptional activity-this repression does not require HDAC7's deacetylase activity. We further show that the chromatin pioneer factor FoxA1, well known for its role in estrogen induction of genes, is recruited to the RPRM promoter, is necessary for repression of RPRM, and interacts with HDAC7. Like other FoxA1 recruitment sites, the RPRM promoter is characterized by H3K4me1/me2. Estrogen treatment causes decreases in H3K4me1/me2 and release of RNA polymerase II (Pol II) from the RPRM proximal promoter. Overall, these data implicate a novel role for HDAC7 and FoxA1 in estrogen repression of RPRM, a mechanism which could potentially be generalized to many more estrogen-repressed genes and hence be important in both normal physiology and pathological processes.