DNA methylation is required to maintain both DNA replication timing precision and 3D genome organization integrity.

DNA replication timing and three-dimensional (3D) genome organization are associated with distinct epigenome patterns across large domains. However, whether alterations in the epigenome, in particular cancer-related DNA hypomethylation, affects higher-order levels of genome architecture is still unclear. Here, using Repli-Seq, single-cell Repli-Seq, and Hi-C, we show that genome-wide methylation loss is associated with both concordant loss of replication timing precision and deregulation of 3D genome organization. Notably, we find distinct disruption in 3D genome compartmentalization, striking gains in cell-to-cell replication timing heterogeneity and loss of allelic replication timing in cancer hypomethylation models, potentially through the gene deregulation of DNA replication and genome organization pathways. Finally, we identify ectopic H3K4me3-H3K9me3 domains from across large hypomethylated domains, where late replication is maintained, which we purport serves to protect against catastrophic genome reorganization and aberrant gene transcription. Our results highlight a potential role for the methylome in the maintenance of 3D genome regulation.

DNA Hypermethylation Encroachment at CpG Island Borders in Cancer Is Predisposed by H3K4 Monomethylation Patterns.

Promoter CpG islands are typically unmethylated in normal cells, but in cancer a proportion are subject to hypermethylation. Using methylome sequencing we identified CpG islands that display partial methylation encroachment across the 5' or 3' CpG island borders. CpG island methylation encroachment is widespread in prostate and breast cancer and commonly associates with gene suppression. We show that the pattern of H3K4me1 at CpG island borders in normal cells predicts the different modes of cancer CpG island hypermethylation. Notably, genetic manipulation of Kmt2d results in concordant alterations in H3K4me1 levels and CpG island border DNA methylation encroachment. Our findings suggest a role for H3K4me1 in the demarcation of CpG island methylation borders in normal cells, which become eroded in cancer.

Replication timing and epigenome remodelling are associated with the nature of chromosomal rearrangements in cancer.

DNA replication timing is known to facilitate the establishment of the epigenome, however, the intimate connection between replication timing and changes to the genome and epigenome in cancer remain largely uncharacterised. Here, we perform Repli-Seq and integrated epigenome analyses and demonstrate that genomic regions that undergo long-range epigenetic deregulation in prostate cancer also show concordant differences in replication timing. A subset of altered replication timing domains are conserved across cancers from different tissue origins. Notably, late-replicating regions in cancer cells display a loss of DNA methylation, and a switch in heterochromatin features from H3K9me3-marked constitutive to H3K27me3-marked facultative heterochromatin. Finally, analysis of 214 prostate and 35 breast cancer genomes reveal that late-replicating regions are prone to cis and early-replication to trans chromosomal rearrangements. Together, our data suggests that the nature of chromosomal rearrangement in cancer is related to the spatial and temporal positioning and altered epigenetic states of early-replicating compared to late-replicating loci.

Enduring epigenetic landmarks define the cancer microenvironment.

The growth and progression of solid tumors involves dynamic cross-talk between cancer epithelium and the surrounding microenvironment. To date, molecular profiling has largely been restricted to the epithelial component of tumors; therefore, features underpinning the persistent protumorigenic phenotype of the tumor microenvironment are unknown. Using whole-genome bisulfite sequencing, we show for the first time that cancer-associated fibroblasts (CAFs) from localized prostate cancer display remarkably distinct and enduring genome-wide changes in DNA methylation, significantly at enhancers and promoters, compared to nonmalignant prostate fibroblasts (NPFs). Differentially methylated regions associated with changes in gene expression have cancer-related functions and accurately distinguish CAFs from NPFs. Remarkably, a subset of changes is shared with prostate cancer epithelial cells, revealing the new concept of tumor-specific epigenome modifications in the tumor and its microenvironment. The distinct methylome of CAFs provides a novel epigenetic hallmark of the cancer microenvironment and promises new biomarkers to improve interpretation of diagnostic samples.

Comprehensive evaluation of genome-wide 5-hydroxymethylcytosine profiling approaches in human DNA.

BACKGROUND: The discovery that 5-methylcytosine (5mC) can be oxidized to 5-hydroxymethylcytosine (5hmC) by the ten-eleven translocation (TET) proteins has prompted wide interest in the potential role of 5hmC in reshaping the mammalian DNA methylation landscape. The gold-standard bisulphite conversion technologies to study DNA methylation do not distinguish between 5mC and 5hmC. However, new approaches to mapping 5hmC genome-wide have advanced rapidly, although it is unclear how the different methods compare in accurately calling 5hmC. In this study, we provide a comparative analysis on brain DNA using three 5hmC genome-wide approaches, namely whole-genome bisulphite/oxidative bisulphite sequencing (WG Bis/OxBis-seq), Infinium HumanMethylation450 BeadChip arrays coupled with oxidative bisulphite (HM450K Bis/OxBis) and antibody-based immunoprecipitation and sequencing of hydroxymethylated DNA (hMeDIP-seq). We also perform loci-specific TET-assisted bisulphite sequencing (TAB-seq) for validation of candidate regions. RESULTS: We show that whole-genome single-base resolution approaches are advantaged in providing precise 5hmC values but require high sequencing depth to accurately measure 5hmC, as this modification is commonly in low abundance in mammalian cells. HM450K arrays coupled with oxidative bisulphite provide a cost-effective representation of 5hmC distribution, at CpG sites with 5hmC levels >~10%. However, 5hmC analysis is restricted to the genomic location of the probes, which is an important consideration as 5hmC modification is commonly enriched at enhancer elements. Finally, we show that the widely used hMeDIP-seq method provides an efficient genome-wide profile of 5hmC and shows high correlation with WG Bis/OxBis-seq 5hmC distribution in brain DNA. However, in cell line DNA with low levels of 5hmC, hMeDIP-seq-enriched regions are not detected by WG Bis/OxBis or HM450K, either suggesting misinterpretation of 5hmC calls by hMeDIP or lack of sensitivity of the latter methods. CONCLUSIONS: We highlight both the advantages and caveats of three commonly used genome-wide 5hmC profiling technologies and show that interpretation of 5hmC data can be significantly influenced by the sensitivity of methods used, especially as the levels of 5hmC are low and vary in different cell types and different genomic locations.

Genome-scale methylation assessment did not identify prognostic biomarkers in oral tongue carcinomas.

BACKGROUND: DNA methylation profiling of heterogeneous head and neck squamous cell carcinoma (HNSCC) cohorts has been reported to predict patient outcome. We investigated if a prognostic DNA methylation profile could be found in tumour tissue from a single uniform subsite, the oral tongue. The methylation status of 109 comprehensively annotated oral tongue squamous cell carcinoma (OTSCC) formalin-fixed paraffin-embedded (FFPE) samples from a single institution were examined with the Illumina HumanMethylation450K (HM450K) array. Data pre-processing, quality control and analysis were performed using R packages. Probes mapping to SNPs, sex chromosomes and unreliable probes were accounted for prior to downstream analyses. The relationship between methylation and patient survival was examined using both agnostic approaches and feature selection. The cohort was enlarged by incorporation of 331 The Cancer Genome Atlas (TCGA) HNSCC samples, which included 91 TCGA OTSCC samples with HM450K and survival data available. RESULTS: Given the use of FFPE-derived DNA, we defined different cohorts for separate analyses. Overall, similar results were found between cohorts. With an unsupervised approach, no distinct hypermethylated group of samples was identified and nor was a prognostic methylation profile identified. The use of multiple downstream feature selection approaches, including a linear models for microarray data (LIMMA), centroid feature selection (CFS), and recursive feature elimination (RFE) support vector machines, similarly failed to identify a significant methylation signature informative for patient prognosis or any clinicopathological data available. Furthermore, we were unable to confirm the prognostic methylation profiles or specific prognostic loci reported within the literature for HNSCC. CONCLUSIONS: With genome-scale assessment of DNA methylation using HM450K in one of the largest OTSCC cohorts to date, we were unable to identify a hypermethylated group of tumours or a prognostic methylation signature. This suggests that either DNA methylation in isolation is not likely to be of prognostic value or larger cohorts are required to identify such a biomarker for OTSCC.

Genome-wide DNA methylation profiling in triple-negative breast cancer reveals epigenetic signatures with important clinical value.

Analysis of cancer methylomes has dramatically changed our concept of the potential of diagnostic and prognostic methylation biomarkers in disease stratification. Through whole-genome methylation capture sequencing of triple-negative breast cancers (TNBCs) we recently identified differentially methylated regions with diagnostic and prognostic value that promise to stratify TNBCs for more personalized management.

Three-dimensional disorganization of the cancer genome occurs coincident with long-range genetic and epigenetic alterations.

A three-dimensional chromatin state underpins the structural and functional basis of the genome by bringing regulatory elements and genes into close spatial proximity to ensure proper, cell-type-specific gene expression profiles. Here, we performed Hi-C chromosome conformation capture sequencing to investigate how three-dimensional chromatin organization is disrupted in the context of copy-number variation, long-range epigenetic remodeling, and atypical gene expression programs in prostate cancer. We find that cancer cells retain the ability to segment their genomes into megabase-sized topologically associated domains (TADs); however, these domains are generally smaller due to establishment of additional domain boundaries. Interestingly, a large proportion of the new cancer-specific domain boundaries occur at regions that display copy-number variation. Notably, a common deletion on 17p13.1 in prostate cancer spanning the TP53 tumor suppressor locus results in bifurcation of a single TAD into two distinct smaller TADs. Change in domain structure is also accompanied by novel cancer-specific chromatin interactions within the TADs that are enriched at regulatory elements such as enhancers, promoters, and insulators, and associated with alterations in gene expression. We also show that differential chromatin interactions across regulatory regions occur within long-range epigenetically activated or silenced regions of concordant gene activation or repression in prostate cancer. Finally, we present a novel visualization tool that enables integrated exploration of Hi-C interaction data, the transcriptome, and epigenome. This study provides new insights into the relationship between long-range epigenetic and genomic dysregulation and changes in higher-order chromatin interactions in cancer.

DNA methylation of oestrogen-regulated enhancers defines endocrine sensitivity in breast cancer.

Expression of oestrogen receptor (ESR1) determines whether a breast cancer patient receives endocrine therapy, but does not guarantee patient response. The molecular factors that define endocrine response in ESR1-positive breast cancer patients remain poorly understood. Here we characterize the DNA methylome of endocrine sensitivity and demonstrate the potential impact of differential DNA methylation on endocrine response in breast cancer. We show that DNA hypermethylation occurs predominantly at oestrogen-responsive enhancers and is associated with reduced ESR1 binding and decreased gene expression of key regulators of ESR1 activity, thus providing a novel mechanism by which endocrine response is abated in ESR1-positive breast cancers. Conversely, we delineate that ESR1-responsive enhancer hypomethylation is critical in transition from normal mammary epithelial cells to endocrine-responsive ESR1-positive cancer. Cumulatively, these novel insights highlight the potential of ESR1-responsive enhancer methylation to both predict ESR1-positive disease and stratify ESR1-positive breast cancer patients as responders to endocrine therapy.

MicroRNA profiling of the pubertal mouse mammary gland identifies miR-184 as a candidate breast tumour suppressor gene.

INTRODUCTION: The study of mammalian development has offered many insights into the molecular aetiology of cancer. We previously used analysis of mammary morphogenesis to discover a critical role for GATA-3 in mammary developmental and carcinogenesis. In recent years an important role for microRNAs (miRNAs) in a myriad of cellular processes in development and in oncogenesis has emerged. METHODS: microRNA profiling was conducted on stromal and epithelial cellular subsets microdissected from the pubertal mouse mammary gland. miR-184 was reactivated by transient or stable overexpression in breast cancer cell lines and examined using a series of in vitro (proliferation, tumour-sphere and protein synthesis) assays. Orthotopic xenografts of breast cancer cells were used to assess the effect of miR-184 on tumourigenesis as well as distant metastasis. Interactions between miR-184 and its putative targets were assessed by quantitative PCR, microarray, bioinformatics and 3' untranslated region Luciferase reporter assay. The methylation status of primary patient samples was determined by MBD-Cap sequencing. Lastly, the clinical prognostic significance of miR-184 putative targets was assessed using publicly available datasets. RESULTS: A large number of microRNA were restricted in their expression to specific tissue subsets. MicroRNA-184 (miR-184) was exclusively expressed in epithelial cells and markedly upregulated during differentiation of the proliferative, invasive cells of the pubertal terminal end bud (TEB) into ductal epithelial cells in vivo. miR-184 expression was silenced in mouse tumour models compared to non-transformed epithelium and in a majority of breast cancer cell line models. Ectopic reactivation of miR-184 inhibited the proliferation and self-renewal of triple negative breast cancer (TNBC) cell lines in vitro and delayed primary tumour formation and reduced metastatic burden in vivo. Gene expression studies uncovered multi-factorial regulation of genes in the AKT/mTORC1 pathway by miR-184. In clinical breast cancer tissues, expression of miR-184 is lost in primary TNBCs while the miR-184 promoter is methylated in a subset of lymph node metastases from TNBC patients. CONCLUSIONS: These studies elucidate a new layer of regulation in the PI3K/AKT/mTOR pathway with relevance to mammary development and tumour progression and identify miR-184 as a putative breast tumour suppressor.

Coordinated epigenetic remodelling of transcriptional networks occurs during early breast carcinogenesis.

BACKGROUND: Dysregulation of the epigenome is a common event in malignancy; however, deciphering the earliest cancer-associated epigenetic events remains a challenge. Cancer epigenome studies to date have primarily utilised cancer cell lines or clinical samples, where it is difficult to identify the initial epigenetic lesions from those that occur over time. Here, we analysed the epigenome of human mammary epithelial cells (HMEC) and a matched variant cell population (vHMEC) that have spontaneously escaped senescence and undergone partial carcinogenic transformation. Using this model of basal-like breast carcinogenesis, we provide striking new insights into the very first epigenetic changes that occur during the initial stages of malignancy. RESULTS: The first phase of malignancy is defined by coordinated changes in the epigenome. At the chromatin level, this is embodied in long-range epigenetic deregulation, which involves the concomitant but atypical acquisition or loss of active and repressive histone modifications across large regional blocks. Changes in DNA methylation also occurs in a highly coordinated manner. We identified differentially methylated regions (DMRs) in the very earliest passages of vHMECs. Notably, we find that differential methylation targets loci regulated by key transcription factors including p53, AHR and E2F family members suggesting that epigenetic deregulation of transcription factor binding is a key event in breast carcinogenesis. Interestingly, DMRs identified in vHMEC are extensively methylated in breast cancer, with hypermethylation frequently encroaching into neighbouring regions. A subset of vHMEC DMRs exhibited a strong basal-like cancer specific hypermethylation. CONCLUSIONS: Here, we generated epigenome-wide maps of the earliest phase of breast malignancy and show long-range epigenetic deregulation and coordinated DNA hypermethylation targets loci regulated by key transcription factors. These findings support a model where induction of breast cancer occurs through epigenetic disruption of transcription factor binding leading to deregulation of cancer-associated transcriptional networks. With their stability and very early occurrence, vHMECs hypermethylated loci could serve as excellent biomarkers for the initial detection of basal breast cancer.

Methylome sequencing in triple-negative breast cancer reveals distinct methylation clusters with prognostic value.

Epigenetic alterations in the cancer methylome are common in breast cancer and provide novel options for tumour stratification. Here, we perform whole-genome methylation capture sequencing on small amounts of DNA isolated from formalin-fixed, paraffin-embedded tissue from triple-negative breast cancer (TNBC) and matched normal samples. We identify differentially methylated regions (DMRs) enriched with promoters associated with transcription factor binding sites and DNA hypersensitive sites. Importantly, we stratify TNBCs into three distinct methylation clusters associated with better or worse prognosis and identify 17 DMRs that show a strong association with overall survival, including DMRs located in the Wilms tumour 1 (WT1) gene, bi-directional-promoter and antisense WT1-AS. Our data reveal that coordinated hypermethylation can occur in oestrogen receptor-negative disease, and that characterizing the epigenetic framework provides a potential signature to stratify TNBCs. Together, our findings demonstrate the feasibility of profiling the cancer methylome with limited archival tissue to identify regulatory regions associated with cancer.

A bioengineered microenvironment to quantitatively measure the tumorigenic properties of cancer-associated fibroblasts in human prostate cancer.

Stromal-epithelial cell interactions play an important role in cancer and the tumor stroma is regarded as a therapeutic target. In vivo xenografting is commonly used to study cellular interactions not mimicked or quantified in conventional 2D culture systems. To interrogate the effects of tumor stroma (cancer-associated fibroblasts or CAFs) on epithelia, we created a bioengineered microenvironment using human prostatic tissues. Patient-matched CAFs and non-malignant prostatic fibroblasts (NPFs) from men with moderate (Gleason 7) and aggressive (Gleason 8-9 or castrate-resistant) prostate cancer were cultured with non-tumorigenic BPH-1 epithelial cells. Changes in the morphology, motility and phenotype of BPH-1 cells in response to CAFs and NPFs were analyzed using immunofluorescence and quantitative cell morphometric analyses. The matrix protein gene expression of CAFs, with proven tumorigenicity in vivo, had a significantly different gene expression profile of matrix proteins compared to patient matched NPFs. In co-culture with CAFs (but not NPFs), BPH-1 cells had a more invasive, elongated phenotype with increased motility and a more directed pattern of cell migration. CAFs from more aggressive tumors (Gleason 8-9 or CRPC) were not quantitatively different to moderate grade CAFs. Overall, our bioengineered microenvironment provides a novel 3D in vitro platform to systematically investigate the effects of tumor stroma on prostate cancer progression.